Researcher profile

Stephen H. Lubow

Stephen H. Lubow contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

8 published item(s)

preprint2022arXiv

Accretion onto a binary from a polar circumbinary disc

We present hydrodynamical simulations to model the accretion flow from a polar circumbinary disc onto a high eccentricity ($e=0.78$) binary star system with near unity mass ratio ($q=0.83$), as a model for binary HD 98800 BaBb. We compare the polar circumbinary disc accretion flow with the previously studied coplanar case. In the coplanar case, the circumbinary disc becomes eccentric and the accretion alternates from being dominant onto one binary member to the other. For the polar disc case involving a highly eccentric binary, we find that the circumbinary disc retains its initially low eccentricity and that the primary star accretion rate is always about the same as the secondary star accretion rate. Recent observations of the binary HD 98800 BaBb, which has a polar circumbinary disc, have been used to determine the value of the $\rm Hα$ flux from the brighter component. From this value, we infer that the accretion rate is much lower than for typical T Tauri stars. The eccentric orbit of the outer companion HD 98800 A increases the accretion rate onto HD 98800 B by $\sim 20$ per cent after each periastron passage. Our hydrodynamical simulations are unable to explain such a low accretion rate unless the disc viscosity parameter is very small, $α< 10^{-5}$. Additional observations of this system would be useful to check on this low accretion rate.

preprint2022arXiv

Circumbinary disk evolution in the presence of an outer companion star

We consider a hierarchical triple system consisting of an inner eccentric binary with an outer companion. A highly misaligned circumbinary disk around the inner binary is subject to two competing effects: (i) nodal precession about the inner binary eccentricity vector that leads to an increase in misalignment (polar alignment) and (ii) Kozai-Lidov (KL) oscillations of eccentricity and inclination driven by the outer companion that leads to a reduction in the misalignment. The outcome depends upon the ratio of the timescales of these effects. If the inner binary torque dominates, then the disk aligns to a polar orientation. If the outer companion torque dominates, then the disk undergoes KL oscillations. In that case, the highly eccentric and misaligned disk is disrupted and accreted by the inner binary, while some mass is transferred to the outer companion. However, when the torques are similar, the outer parts of the circumbinary disk can undergo large eccentricity oscillations while the inclination remains close to the polar orientation. The range of initial disk inclinations that evolve to a polar orientation is smaller in the presence of the outer companion. Disk breaking is also more likely, at least temporarily, during the polar alignment process. The stellar orbits in HD 98800 have parameters such that polar alignment of the circumbinary disk is expected. In the absence of the gas, solid particles are unstable at much smaller radii than the gas disk inner tidal truncation radius because KL driven eccentricity leads to close encounters with the binary.

preprint2022arXiv

Eccentric dust ring formation in Kozai-Lidov gas disks

A highly misaligned gas disk around one component of a binary star system can undergo global Kozai-Lidov (KL) oscillations for which the disk inclination and eccentricity are exchanged. With hydrodynamical simulations of a gas and dust disk we explore the effects of these oscillations on the dust density distribution. For dust that is marginally coupled to the gas (${\rm St} \approx 1$), we find that the dust undergoes similar dynamical behaviour to the gas disk but the radial distribution of dust may be very different to the gas. The inward radial drift of the dust is faster in an eccentric disk leading to a smaller outer dust disk radius. The dust breaks into multiple narrow eccentric rings during the highly eccentric disk phase. Eccentric dust ring formation may have significant implications for the formation of planets in misaligned disks. We suggest that multiple dust rings may generally occur within gas disks that have sufficiently strong eccentricity peaks at intermediate radii.

preprint2022arXiv

Improvements to Pan-STARRS1 Astrometry: II. Corrections for Differential Chromatic Refraction

In a previous paper, we applied the Gaia DR2 catalog to improve the astrometric accuracy of about 1.7 billion objects in Pan-STARRS1 Data Release 2 (PS1 DR2). We report here on further improvements made by utilizing Gaia EDR3 and correcting for the effects of differential chromatic refraction (DCR) in declination. We extend the correction algorithm in Paper 1 by iteratively subtracting color- and declination-dependent PS1/Gaia EDR3 declination residuals. We determine the astrometric improvement for ~440 million reference objects that are point-like and cross-match to Gaia EDR3. For this set of objects, Gaia EDR3 provides a ~3% improvement in PS1 astrometry over Gaia DR2, and DCR corrections provide an additional ~5% improvement. DCR corrections increase substantially for objects observed away from the zenith. DCR corrections lead to an astrometric improvement of ~30% for blue objects (0<g-i<1) that are 50° away from the zenith. The amplitude of systematic astrometric errors from these effects is substantially reduced to less than 1 mas for objects with PS1 colors in the range 0 < g-i < 4.5, which makes this a useful astrometric reference catalog in fields where there are few Gaia stars. The improved astrometric data will be available through the Mikulski Archive for Space Telescopes PS1 catalog interfaces.

preprint2020arXiv

Formation of the polar debris disc around 99 Herculis

We investigate the formation mechanism for the observed nearly polar aligned (perpendicular to the binary orbital plane) debris ring around the eccentric orbit binary 99 Herculis. An initially inclined nonpolar debris ring or disc will not remain flat and will not evolve to a polar configuration, due to the effects of differential nodal precession that alter its flat structure. However, a gas disc with embedded well coupled solids around the eccentric binary may evolve to a polar configuration as a result of pressure forces that maintain the disc flatness and as a result of viscous dissipation that allows the disc to increase its tilt. Once the gas disc disperses, the debris disc is in a polar aligned state in which there is little precession. We use three-dimensional hydrodynamical simulations, linear theory, and particle dynamics to study the evolution of a misaligned circumbinary gas disc and explore the effects of the initial disc tilt, mass, and size. We find that for a wide range of parameter space, the polar alignment timescale is shorter than the lifetime of the gas disc. Using the observed level of alignment of 3 deg. from polar, we place an upper limit on the mass of the gas disc of about 0.014 M_sun at the time of gas dispersal. We conclude that the polar debris disc around 99 Her can be explained as the result of an initially moderately inclined gas disc with embedded solids. Such a disc may provide an environment for the formation of polar planets.

preprint2020arXiv

Polar planets around highly eccentric binaries are the most stable

We study the orbital stability of a non-zero mass, close-in circular orbit planet around an eccentric orbit binary for various initial values of the binary eccentricity, binary mass fraction, planet mass, planet semi--major axis, and planet inclination by means of numerical simulations that cover $5 \times 10^4$ binary orbits. For small binary eccentricity, the stable orbits that extend closest to the binary (most stable orbits) are nearly retrograde and circulating. For high binary eccentricity, the most stable orbits are highly inclined and librate near the so-called generalised polar orbit which is a stationary orbit that is fixed in the frame of the binary orbit. For more extreme mass ratio binaries, there is a greater variation in the size of the stability region (defined by initial orbital radius and inclination) with planet mass and initial inclination, especially for low binary eccentricity. For low binary eccentricity, inclined planet orbits may be unstable even at large orbital radii (separation $> 5 \,a_{\rm b}$). The escape time for an unstable planet is generally shorter around an equal mass binary compared with an unequal mass binary. Our results have implications for circumbinary planet formation and evolution and will be helpful for understanding future circumbinary planet observations.

preprint2020arXiv

Sustained Kozai-Lidov oscillations in misaligned circumstellar gas disks

A disk around one component of a binary star system with sufficiently high inclination can undergo Kozai-Lidov (KL) oscillations during which the disk inclination and disk eccentricity are exchanged. Previous studies show that without a source of accretion, KL unstable disks exhibit damped oscillations, due to viscous dissipation, that leave the disk stable near or below the critical inclination for KL oscillations. With three-dimensional hydrodynamical simulations we show that a highly misaligned circumbinary disk that flows onto the binary components forms highly inclined circumstellar disks around each component. We show that a continuous infall of highly inclined material allows the KL oscillations to continue. The KL disk oscillations produce shocks and eccentricity growth in the circumstellar disks that affect the conditions for planet formation.

preprint2019arXiv

Multi-planet disc interactions in binary systems

We investigate the evolution of a multi--planet--disc system orbiting one component of a binary star system. The planet--disc system is initially coplanar but misaligned to the binary orbital plane. The planets are assumed to be giants that open gaps in the disc. We first study the role of the disc in shaping the mutual evolution of the two planets using a secular model for low initial tilt. In general we find that the planets and the disc do not remain coplanar, in agreement with Lubow & Martin (2016) for the single planet case. Instead, the planets and the disc undergo tilt oscillations. A high mass disc between the two planets causes the planets and the disc to nodally precess at the same average rate but they are generally misaligned. The amplitude of the tilt oscillations between the planets is larger while the disc is present. We then consider higher initial tilts using hydrodynamical simulations and explore the possibility of the formation of eccentric Kozai-Lidov (KL) planets. We find that the inner planet&#39;s orbit undergoes eccentricity growth for a large range of disc masses and initial misalignments. For a low disc mass and large initial misalignment both planets and the disc can undergo KL oscillations. Furthermore, we find that sufficiently massive discs can cause the inner planet to increase its inclination beyond $90^{\circ}$ and therefore to orbit the binary in a retrograde fashion. The results have important implications for the explanation of very eccentric planets and retrograde planets observed in multi-planet systems.