-
Decretion disc size in Be/X-ray binaries depends upon the disc aspect ratio
Authors:
Rebecca G. Martin,
Stephen H. Lubow,
Philip J. Armitage,
Daniel J. Price
Abstract:
With three-dimensional hydrodynamical simulations we show that the size of the decretion disc and the structure of the accretion flow onto the neutron star in a Be/X-ray binary strongly depends upon the disc aspect ratio, $H/R$. We simulate a Be star disc that is coplanar to the orbit of a circularly or moderately eccentric neutron star companion, thereby maximising the effects of tidal truncation…
▽ More
With three-dimensional hydrodynamical simulations we show that the size of the decretion disc and the structure of the accretion flow onto the neutron star in a Be/X-ray binary strongly depends upon the disc aspect ratio, $H/R$. We simulate a Be star disc that is coplanar to the orbit of a circularly or moderately eccentric neutron star companion, thereby maximising the effects of tidal truncation. For low disc aspect ratio, $H/R\lesssim 0.1$, the disc is efficiently tidally truncated by the neutron star. Most material that escapes the Roche lobe of the Be star is accreted by the neutron star through tidal streams. For larger disc aspect ratio, the outflow rate through the Be star disc is higher, tidal truncation becomes inefficient, the disc fills the Roche lobe and extends to the orbit of the companion. Some material escapes the binary as a gas stream that begins near the L2 point. While the accretion rate onto the neutron star is higher, the fraction of the outflow that is accreted by the neutron star is smaller. Low density Be star discs are expected to be approximately isothermal, such that $H/R$ increases with radius. Tidal truncation is therefore weaker for larger separation binaries, and lower mass primaries.
△ Less
Submitted 27 April, 2024;
originally announced April 2024.
-
Inclination instability of circumbinary planets
Authors:
Stephen H. Lubow,
Anna C. Childs,
Rebecca G. Martin
Abstract:
We analyze a tilt instability of the orbit of an outer planet in a two planet circumbinary system that we recently reported. The binary is on an eccentric orbit and the inner circumbinary planet is on a circular polar orbit that causes the the binary to undergo apsidal precession. The outer circumbinary planet is initially on a circular or eccentric orbit that is coplanar with respect to the binar…
▽ More
We analyze a tilt instability of the orbit of an outer planet in a two planet circumbinary system that we recently reported. The binary is on an eccentric orbit and the inner circumbinary planet is on a circular polar orbit that causes the the binary to undergo apsidal precession. The outer circumbinary planet is initially on a circular or eccentric orbit that is coplanar with respect to the binary. We apply a Hamiltonian in quadrupole order of the binary potential to show that the tilt instability is the result of a secular resonance in which the apsidal precession rate of the binary matches the nodal precession rate of the outer planet. Resonance is possible because the polar inner planet causes the apsidal precession of the binary to be retrograde. The outer planet periodically undergoes large tilt oscillations for which we analytically determine the initial evolution and maximum inclination. Following a typically relatively short adjustment phase, the tilt grows exponentially in time at a characteristic rate that is of order the absolute value of the binary apsidal precession rate. The analytic results agree well with numerical simulations. This instability is analogous to the Kozai-Lidov instability, but applied to a circumbinary object. The instability fails to operate if the binary mass ratio is too extreme. The instability occurs even if the outer planet is instead an object of stellar mass and involves tilt oscillations of the inner binary.
△ Less
Submitted 15 April, 2024;
originally announced April 2024.
-
Relativistic Effects on Circumbinary Disk Evolution: Breaking the Polar Alignment around Eccentric Black Hole Binary Systems
Authors:
Anna C. Childs,
Rebecca G. Martin,
C. J. Nixon,
Aaron M. Geller,
Stephen H. Lubow,
Zhaohuan Zhu,
Stephen Lepp
Abstract:
We study the effects of general relativity (GR) on the evolution and alignment of circumbinary disks around binaries on all scales. We implement relativistic apsidal precession of the binary into the hydrodynamics code {\sc phantom}. We find that the effects of GR can suppress the stable polar alignment of a circumbinary disk, depending on how the relativistic binary apsidal precession timescale c…
▽ More
We study the effects of general relativity (GR) on the evolution and alignment of circumbinary disks around binaries on all scales. We implement relativistic apsidal precession of the binary into the hydrodynamics code {\sc phantom}. We find that the effects of GR can suppress the stable polar alignment of a circumbinary disk, depending on how the relativistic binary apsidal precession timescale compares to the disk nodal precession timescale. Studies of circumbinary disk evolution typically ignore the effects of GR which is an appropriate simplification for low mass or widely separated binary systems. In this case, polar alignment occurs providing that the disks initial misalignment is sufficiently large. However, systems with a very short relativistic precession timescale cannot polar align and instead move toward coplanar alignment. In the intermediate regime where the timescales are similar, the outcome depends upon the properties of the disk. Polar alignment is more likely in the wavelike disk regime (where the disk viscosity parameter is less than the aspect ratio, $α<H/r$) since the disk is in good radial communication. In the viscous disk regime disk breaking is more likely. Multiple rings can destructively interact with one another resulting in short disk lifetimes, and the disk moving towards coplanar alignment. Around main-sequence star or stellar mass black hole binaries, polar alignment may be suppressed far from the binary but in general the inner parts of the disk can align to polar. Polar alignment may be completely suppressed for disks around supermassive black holes for close binary separations.
△ Less
Submitted 14 December, 2023;
originally announced December 2023.
-
Retrograde discs around one component of a binary are unstable to tilting
Authors:
Madeline Overton,
Rebecca G. Martin,
Stephen H. Lubow,
Stephen Lepp
Abstract:
With hydrodynamic simulations we show that a coplanar disc around one component of a binary can be unstable to global tilting when the disc orbits in a retrograde direction relative to the binary. The disc experiences the largest inclination growth relative to the binary orbit in the outermost radii of the disc, closest to the companion. This tilt instability also occurs for test particles. A retr…
▽ More
With hydrodynamic simulations we show that a coplanar disc around one component of a binary can be unstable to global tilting when the disc orbits in a retrograde direction relative to the binary. The disc experiences the largest inclination growth relative to the binary orbit in the outermost radii of the disc, closest to the companion. This tilt instability also occurs for test particles. A retrograde disc is much larger than a prograde disc since it is not tidally truncated and instead spreads outwards to the orbit of the companion. The coplanar retrograde disc remains circular while a coplanar prograde disc can become eccentric. We suggest that the inclination instability is due to a disc resonance caused by the interaction of the tilt with the tidal field of the binary. This model is applicable to Be/X-ray binaries in which the Be star disc may be retrograde relative to the binary orbit if there was a sufficiently strong kick from the supernova that formed the neutron star companion. The accretion on to the neutron star and the resulting X-ray outbursts are weaker in the retrograde case compared to the prograde case.
△ Less
Submitted 17 November, 2023;
originally announced November 2023.
-
AC Her: Evidence of the first polar circumbinary planet
Authors:
Rebecca G. Martin,
Stephen H. Lubow,
David Vallet,
Narsireddy Anugu,
Douglas R. Gies
Abstract:
We examine the geometry of the post-asymptotic giant branch (AGB) star binary AC Her and its circumbinary disk. We show that the observations describe a binary orbit that is perpendicular to the disk with an angular momentum vector that is within $9^\circ$ of the binary eccentricity vector, meaning that the disk is close to a stable polar alignment. The most likely explanation for the very large i…
▽ More
We examine the geometry of the post-asymptotic giant branch (AGB) star binary AC Her and its circumbinary disk. We show that the observations describe a binary orbit that is perpendicular to the disk with an angular momentum vector that is within $9^\circ$ of the binary eccentricity vector, meaning that the disk is close to a stable polar alignment. The most likely explanation for the very large inner radius of the dust is a planet within the circumbinary disk. This is therefore both the first reported detection of a polar circumbinary disk around a post-AGB binary and the first evidence of a polar circumbinary planet. We consider the dynamical constraints on the circumbinary disk size and mass. The polar circumbinary disk feeds circumstellar disks with gas on orbits that are highly inclined with respect to the binary orbit plane. The resulting circumstellar disk inclination could be anywhere from coplanar to polar depending upon the competition between the mass accretion and binary torques.
△ Less
Submitted 30 October, 2023;
originally announced October 2023.
-
Tilted circumbinary planetary systems as efficient progenitors of free-floating planets
Authors:
Cheng Chen,
Rebecca G. Martin,
Stephen H. Lubow,
C. J. Nixon
Abstract:
The dominant mechanism for generating free-floating planets has so far remained elusive. One suggested mechanism is that planets are ejected from planetary systems due to planet-planet interactions. However, instability around a single star requires a very compactly spaced planetary system. We find that around binary star systems instability can occur even with widely separated planets that are on…
▽ More
The dominant mechanism for generating free-floating planets has so far remained elusive. One suggested mechanism is that planets are ejected from planetary systems due to planet-planet interactions. However, instability around a single star requires a very compactly spaced planetary system. We find that around binary star systems instability can occur even with widely separated planets that are on tilted orbits relative to the binary orbit due to combined effects of planet-binary and planet-planet interactions, especially if the binary is on an eccentric orbit. We investigate the orbital stability of planetary systems with various planet masses and architectures. We find that the stability of the system depends upon the mass of the highest mass planet. The order of the planets in the system does not significantly affect stability but, generally, the most massive planet remains stable and the lower mass planets are ejected. The minimum planet mass required to trigger the instability is about that of Neptune for a circular orbit binary and a super-Earth of about $10$ Earth masses for highly eccentric binaries. Hence, we suggest that planet formation around misaligned binaries can be an efficient formation mechanism for free-floating planets. While most observed free-floating planets are giant planets, we predict that there should be more low-mass free floating planets that are as yet unobserved than higher mass planets.
△ Less
Submitted 24 October, 2023;
originally announced October 2023.
-
Tidal truncation of circumplanetary disks fails above a critical disk aspect ratio
Authors:
Rebecca G. Martin,
Philip J. Armitage,
Stephen H. Lubow,
Daniel J. Price
Abstract:
We use numerical simulations of circumplanetary disks to determine the boundary between disks that are radially truncated by the tidal potential, and those where gas escapes the Hill sphere. We consider a model problem, in which a coplanar circumplanetary disk is resupplied with gas at an injection radius smaller than the Hill radius. We evolve the disk using the PHANTOM Smoothed Particle Hydrodyn…
▽ More
We use numerical simulations of circumplanetary disks to determine the boundary between disks that are radially truncated by the tidal potential, and those where gas escapes the Hill sphere. We consider a model problem, in which a coplanar circumplanetary disk is resupplied with gas at an injection radius smaller than the Hill radius. We evolve the disk using the PHANTOM Smoothed Particle Hydrodynamics code until a steady-state is reached. We find that the most significant dependence of the truncation boundary is on the disk aspect ratio $H/R$. Circumplanetary disks are efficiently truncated for $H/R \lesssim 0.2$. For $H/R \simeq 0.3$, up to about half of the injected mass, depending on the injection radius, flows outwards through the decretion disk and escapes. As expected from analytic arguments, the conditions ($H/R$ and Shakura-Sunyaev $α$) required for tidal truncation are independent of planet mass. A simulation with larger $α=0.1$ shows stronger outflow than one with $α=0.01$, but the dependence on transport efficiency is less important than variations of $H/R$. Our results suggest two distinct classes of circumplanetary disks: tidally truncated thin disks with dust-poor outer regions, and thicker actively decreting disks with enhanced dust-to-gas ratios. Applying our results to the PDS 70c system, we predict a largely truncated circumplanetary disk, but it is possible that enough mass escapes to support an outward flow of dust that could explain the observed disk size.
△ Less
Submitted 30 June, 2023;
originally announced June 2023.
-
On the frequencies of circumbinary discs in protostellar systems
Authors:
Daniel Elsender,
Matthew R. Bate,
Ben S. Lakeland,
Eric L. N. Jensen,
Stephen H. Lubow
Abstract:
We report the analysis of circumbinary discs formed in a radiation hydrodynamical simulation of star cluster formation. We consider both pure binary stars and pairs within triple and quadruple systems. The protostellar systems are all young (ages < $10^5$ yrs). We find that the systems that host a circumbinary disc have a median separation of $\approx 11$ au, and the median characteristic radius o…
▽ More
We report the analysis of circumbinary discs formed in a radiation hydrodynamical simulation of star cluster formation. We consider both pure binary stars and pairs within triple and quadruple systems. The protostellar systems are all young (ages < $10^5$ yrs). We find that the systems that host a circumbinary disc have a median separation of $\approx 11$ au, and the median characteristic radius of the discs is $\approx 64$ au. We find that $89$ per cent of pure binaries with semi-major axes $a<1$ au have a circumbinary disc, and the occurrence rate of circumbinary discs is bimodal with log-separation in pure binaries with a second peak at $a \approx 50$ au. Systems with $a>100$ au almost never have a circumbinary disc. The median size of a circumbinary disc is between $\approx 5-6\ a$ depending on the order of the system, with higher order systems having larger discs relative to binary separation. We find the underlying distribution of mutual inclinations between circumbinary discs and binary orbit of both observed and simulated discs to not differ statistically.
△ Less
Submitted 9 June, 2023;
originally announced June 2023.
-
Orbital stability of two circumbinary planets around misaligned eccentric binaries
Authors:
Cheng Chen,
Stephen H. Lubow,
Rebecca G. Martin,
C. J. Nixon
Abstract:
With $n$-body simulations we investigate the stability of tilted circumbinary planetary systems consisting of two nonzero mass planets. The planets are initially in circular orbits that are coplanar to each other, as would be expected if they form in a flat but tilted circumbinary gas disc and decouple from the disc within a time difference that is much less than the disc nodal precession period.…
▽ More
With $n$-body simulations we investigate the stability of tilted circumbinary planetary systems consisting of two nonzero mass planets. The planets are initially in circular orbits that are coplanar to each other, as would be expected if they form in a flat but tilted circumbinary gas disc and decouple from the disc within a time difference that is much less than the disc nodal precession period. We constrain the parameters of stable multiple planet circumbinary systems. Both planet-planet and planet-binary interactions can cause complex planet tilt oscillations which can destabilise the orbits of one or both planets. The system is considerably more unstable than the effects of these individual interactions would suggest, due to the interplay between these two interactions. The stability of the system is sensitive to the binary eccentricity, the orbital tilt and the semi-major axes of the two circumbinary planets. With an inner planet semi-major axis of $5\,a_{\rm b}$, where $a_{\rm b}$ is semi-major axis of the binary, the system is generally stable if the outer planet is located at $\gtrsim 8\,a_{\rm b}$, beyond the 2:1 mean motion resonance with the inner planet. For larger inner planet semi-major axis the system is less stable because the von-Zeipel--Kozai--Lidov mechanism plays a significant role, particularly for low binary-eccentricity cases. For the unstable cases, the most likely outcome is that one planet is ejected and the other remains bound on a highly eccentric orbit. Therefore we suggest that this instability is an efficient mechanism for producing free-floating planets.
△ Less
Submitted 13 March, 2023; v1 submitted 9 March, 2023;
originally announced March 2023.
-
Coplanar circumbinary planets can be unstable to large tilt oscillations in the presence of an inner polar planet
Authors:
Anna C. Childs,
Rebecca G. Martin,
Stephen Lepp,
Stephen H. Lubow,
Aaron M. Geller
Abstract:
Mutually misaligned circumbinary planets may form in a warped or broken gas disc or from later planet-planet interactions. With numerical simulations and analytic estimates we explore the dynamics of two circumbinary planets with a large mutual inclination. A coplanar inner planet causes prograde apsidal precession of the binary and the stationary inclination for the outer planet is higher for lar…
▽ More
Mutually misaligned circumbinary planets may form in a warped or broken gas disc or from later planet-planet interactions. With numerical simulations and analytic estimates we explore the dynamics of two circumbinary planets with a large mutual inclination. A coplanar inner planet causes prograde apsidal precession of the binary and the stationary inclination for the outer planet is higher for larger outer planet orbital radius. In this case a coplanar outer planet always remains coplanar. On the other hand, a polar inner planet causes retrograde apsidal precession of the binary orbit and the stationary inclination is smaller for larger outer planet orbital radius. For a range of outer planet semi-major axes, an initially coplanar orbit is librating meaning that the outer planet undergoes large tilt oscillations. Circumbinary planets that are highly inclined to the binary are difficult to detect -- it is unlikely for a planet to have an inclination below the transit detection limit in the presence of a polar inner planet. These results suggest that there could be a population of circumbinary planets that are undergoing large tilt oscillations.
△ Less
Submitted 21 February, 2023;
originally announced February 2023.
-
Precession and polar alignment of accretion discs in triple (or multiple) stellar systems
Authors:
Simone Ceppi,
Cristiano Longarini,
Giuseppe Lodato,
Nicolás Cuello,
Stephen H. Lubow
Abstract:
We investigate the mechanism of polar alignment for accretion discs in hierarchical systems (HSs) with more than two stars. In eccentric binary systems, low mass discs that are sufficiently tilted to the binary orbit align in a polar configuration with respect to the binary plane by aligning their angular momentum to the binary eccentricity vector. In HSs, secular evolution of the orbital paramete…
▽ More
We investigate the mechanism of polar alignment for accretion discs in hierarchical systems (HSs) with more than two stars. In eccentric binary systems, low mass discs that are sufficiently tilted to the binary orbit align in a polar configuration with respect to the binary plane by aligning their angular momentum to the binary eccentricity vector. In HSs, secular evolution of the orbital parameters makes the eccentricity vector of the system precess with time. This precession undermines the stability of the polar orbit for accretion discs hosted in HSs. We analytically show that the binary criteria for polar alignment derived in the literature are necessary but not sufficient conditions for polar alignment in HSs. Then, we derive an analytical criterion for polar alignment in HSs. In general, we find that discs orbiting the innermost level of a HS can go polar. Conversely, radially extended discs orbiting the outer levels of a HS cannot polarly align and evolve as orbiting around a circular binary. We confirm our findings through detailed numerical simulations. Also, our results are compatible with the observed distribution of disc-orbit mutual inclination. Finally, we compare the observed distribution of disc inclinations in the binary and in the HS populations. Binaries host mainly coplanar discs, while HSs show a wide range of disc inclinations. We suggest that the wider range of inclinations in HSs results from the secular oscillation of their orbital parameters (such as Kozai-Lidov oscillations), rather than from a different initial condition or evolution between HSs and binaries.
△ Less
Submitted 7 February, 2023;
originally announced February 2023.
-
Formation of polar circumstellar discs in binary star systems
Authors:
Jeremy L. Smallwood,
Rebecca G. Martin,
Stephen H. Lubow
Abstract:
We investigate the flow of material from highly misaligned and polar circumbinary discs that feed the formation of circumstellar discs around each binary component. With three-dimensional hydrodynamic simulations we consider equal mass binaries with low eccentricity. We also simulate inclined test particles and highly-misaligned circumstellar discs around one binary component for comparison. Durin…
▽ More
We investigate the flow of material from highly misaligned and polar circumbinary discs that feed the formation of circumstellar discs around each binary component. With three-dimensional hydrodynamic simulations we consider equal mass binaries with low eccentricity. We also simulate inclined test particles and highly-misaligned circumstellar discs around one binary component for comparison. During Kozai-Lidov (KL) cycles, the circumstellar disc structure is altered through exchanges of disc eccentricity with disc tilt. Highly inclined circumstellar discs and test particles around individual binary components can experience very strong KL oscillations. The continuous accretion of highly misaligned material from the circumbinary disc allows the KL oscillations of circumstellar discs to be long-lived. In this process, the circumbinary material is continuously delivered with a high inclination to the lower inclination circumstellar discs. We find that the simulation resolution is important for modeling the longevity of the KL oscillations. An initially polar circumbinary disc forms nearly polar, circumstellar discs that undergo KL cycles. The gas steams accreting onto the polar circumstellar discs vary in tilt during each binary orbital period, which determines how much material is accreted onto the discs. The long-lived KL cycles in polar circumstellar discs may lead to the formation of polar S-type planets in binary star systems.
△ Less
Submitted 27 January, 2023;
originally announced January 2023.
-
Grid-Based Simulations of Polar Circumbinary Disks: Polar Alignment and Vortex Formation
Authors:
Ian Rabago,
Zhaohuan Zhu,
Rebecca G. Martin,
Stephen H. Lubow
Abstract:
We describe the first grid-based simulations of the polar alignment of a circumbinary disk. We simulate the evolution of an inclined disk around an eccentric binary using the grid-based code ATHENA++. The use of a grid-based numerical code allows us to explore lower disk viscosities than have been examined in previous studies. We find that the disk aligns to a polar orientation when the $α$ viscos…
▽ More
We describe the first grid-based simulations of the polar alignment of a circumbinary disk. We simulate the evolution of an inclined disk around an eccentric binary using the grid-based code ATHENA++. The use of a grid-based numerical code allows us to explore lower disk viscosities than have been examined in previous studies. We find that the disk aligns to a polar orientation when the $α$ viscosity is high, while disks with lower viscosity nodally precess with little alignment over 1000 binary orbital periods. The timescales for polar alignment and disk precession are compared as a function of disk viscosity, and are found to be in agreement with previous studies. At very low disk viscosities (e.g. $α= 10^{-5}$), anticyclonic vortices are observed along the inner edge of the disk. These vortices can persist for thousands of binary orbits, creating azimuthally localized overdensities as well as multiple pairs of spiral arms. The vortex is formed at $\sim 3-4$ times the binary semi-major axis, close to the inner edge of the disk, and orbits at roughly the local Keplerian speed. The presence of a vortex in the disk may play an important role in the evolution of circumbinary systems, such as driving episodic accretion and accelerating the formation of polar circumbinary planets.
△ Less
Submitted 18 January, 2023;
originally announced January 2023.
-
Polar circumtriple planets and disks can only form close to a triple star
Authors:
Stephen Lepp,
Rebecca G. Martin,
Stephen H. Lubow
Abstract:
Observations of protoplanetary disks around binary and triple star systems suggest that misalignments between the orbital plane of the stars and the disks are common. Motivated by recent observations of polar circumbinary disks, we explore the possibility for polar circumtriple disks and therefore polar circumtriple planets that could form in such a disk. With n-body simulations and analytic metho…
▽ More
Observations of protoplanetary disks around binary and triple star systems suggest that misalignments between the orbital plane of the stars and the disks are common. Motivated by recent observations of polar circumbinary disks, we explore the possibility for polar circumtriple disks and therefore polar circumtriple planets that could form in such a disk. With n-body simulations and analytic methods we find that the inclusion of the third star, and the associated apsidal precession, significantly reduces the radial range of polar orbits so that circumtriple polar disks and planets can only be found close to the stellar system. Outside of a critical radius, that is typically in the range of 3-10 times the outer binary separation depending upon the binary parameters, the orbits behave the same as they do around a circular orbit binary. For some observed systems that have shorter period inner binaries, the critical radius is considerably larger. If polar circumtriple planets can form, we suggest that it is likely that they form in a disk that was subject to breaking.
△ Less
Submitted 3 January, 2023;
originally announced January 2023.
-
Polar alignment of a massive retrograde circumbinary disc around an eccentric binary
Authors:
Charles P. Abod,
Cheng Chen,
Jeremy Smallwood,
Ian Rabago,
Rebecca G. Martin,
Stephen H. Lubow
Abstract:
A test particle orbit around an eccentric binary has two stationary states in which there is no nodal precession: coplanar and polar. Nodal precession of a misaligned test particle orbit centres on one of these stationary states. A low mass circumbinary disc undergoes the same precession and moves towards one of these states through dissipation within the disc. For a massive particle orbit, the st…
▽ More
A test particle orbit around an eccentric binary has two stationary states in which there is no nodal precession: coplanar and polar. Nodal precession of a misaligned test particle orbit centres on one of these stationary states. A low mass circumbinary disc undergoes the same precession and moves towards one of these states through dissipation within the disc. For a massive particle orbit, the stationary polar alignment occurs at an inclination less than $90^{\circ}$, this is the prograde-polar stationary inclination. A sufficiently high angular momentum particle has an additional higher inclination stationary state, the retrograde-polar stationary inclination. Misaligned particle orbits close to the retrograde-polar stationary inclination are not nested like the orbits close to the other stationary points. We investigate the evolution of a gas disc that begins close to the retrograde-polar stationary inclination. With hydrodynamical disc simulations, we find that the disc moves through the unnested crescent shape precession orbits and eventually moves towards the prograde-polar stationary inclination thus increasing the parameter space over which circumbinary discs move towards polar alignment. If protoplanetary discs form with an isotropic orientation relative to the binary orbit, then polar discs may be more common than coplanar discs around eccentric binaries, even for massive discs. This has implications for the alignment of circumbinary planets.
△ Less
Submitted 28 September, 2022;
originally announced September 2022.
-
Forced Eccentricity in Circumbinary Discs
Authors:
Stephen H. Lubow
Abstract:
We analyze the eccentric response of a low mass coplanar circumbinary disc to secular tidal forcing by a Keplerian eccentric orbit central binary. The disc acquires a forced eccentricity whose magnitude depends on the properties of the binary and disc. The largest eccentricities occur when there is a global apsidal resonance in the disc. The driving frequency by the binary is its apsidal frequency…
▽ More
We analyze the eccentric response of a low mass coplanar circumbinary disc to secular tidal forcing by a Keplerian eccentric orbit central binary. The disc acquires a forced eccentricity whose magnitude depends on the properties of the binary and disc. The largest eccentricities occur when there is a global apsidal resonance in the disc. The driving frequency by the binary is its apsidal frequency that is equal to zero. A global resonance occurs when the disc properties permit the existence of a zero apsidal frequency free eccentric mode. Resonances occur for different free eccentric modes that differ in the number of radial nodes. For a disc not at resonance, the eccentricity distribution has somewhat similar form to the eccentricity distributions in discs at resonance that have the closest matching disc aspect ratios. For higher disc aspect ratios, the forced eccentricity distribution in a 2D disc is similar to that of the fundamental free mode. The forced eccentricity distribution in a 3D disc is similar to that of higher order free modes, not the fundamental mode, unless the disc is very cool. For parameters close to resonance, large phase shifts occur between the disc and binary eccentricities that are locked in phase. Forced eccentricity may play an important role in the evolution of circumbinary discs and their central binaries.
△ Less
Submitted 13 September, 2022;
originally announced September 2022.
-
Improvements to Pan-STARRS1 Astrometry: II. Corrections for Differential Chromatic Refraction
Authors:
Richard L. White,
Stephen H. Lubow,
Bernie Shiao
Abstract:
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…
▽ More
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.
△ Less
Submitted 20 June, 2022;
originally announced June 2022.
-
Accretion onto a binary from a polar circumbinary disc
Authors:
Jeremy L. Smallwood,
Stephen H. Lubow,
Rebecca G. Martin
Abstract:
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 accre…
▽ More
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.
△ Less
Submitted 18 May, 2022;
originally announced May 2022.
-
Circumbinary disk evolution in the presence of an outer companion star
Authors:
Rebecca G. Martin,
Stephen Lepp,
Stephen H. Lubow,
Matthew A. Kenworthy,
Grant M. Kennedy,
David Vallet
Abstract:
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…
▽ More
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.
△ Less
Submitted 14 February, 2022;
originally announced February 2022.
-
Eccentric dust ring formation in Kozai-Lidov gas disks
Authors:
Rebecca G. Martin,
Stephen H. Lubow
Abstract:
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 t…
▽ More
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.
△ Less
Submitted 9 January, 2022;
originally announced January 2022.
-
Orbital dynamics of two circumbinary planets around misaligned eccentric binaries
Authors:
Cheng Chen,
Stephen H. Lubow,
Rebecca G. Martin
Abstract:
We investigate the orbital dynamics of circumbinary planetary systems with two planets around a circular or eccentric orbit binary. The orbits of the two planet are initially circular and coplanar to each other, but misaligned with respect to the binary orbital plane. The binary-planet and planet-planet interactions result in complex planet tilt oscillations. We use analytic models and numerical s…
▽ More
We investigate the orbital dynamics of circumbinary planetary systems with two planets around a circular or eccentric orbit binary. The orbits of the two planet are initially circular and coplanar to each other, but misaligned with respect to the binary orbital plane. The binary-planet and planet-planet interactions result in complex planet tilt oscillations. We use analytic models and numerical simulations to explore the effects of various values of the planet semi-major axes, binary eccentricity, and initial inclination. Around a circular orbit binary, secular tilt oscillations are driven by planet-planet interactions and are periodic. In that case, planets undergo mutual libration if close together and circulation if far apart with an abrupt transition at a critical separation. Around an eccentric orbit binary, secular tilt oscillations are driven by both planet-planet interactions and binary-planet interactions. These oscillations generally display more than one frequency and are generally not periodic. The transition from mutual planet libration to circulation is not sharp and there is a range of separations for which the planets are on orbits that are sometimes mutually librating and sometimes circulating. In addition, at certain separations, there are resonances for which tilt oscillations are complicated but periodic. For planets that are highly misaligned with respect to an eccentric orbit binary, there are stationary (non-oscillating) tilt configurations that are generalisations of polar configurations for the single planet case. Tilt oscillations of highly inclined planets occur for initial tilts that depart from the stationary configuration.
△ Less
Submitted 29 November, 2021;
originally announced November 2021.
-
Interactions of Young Binaries with Disks
Authors:
Stephen H. Lubow,
Pawel Artymowicz
Abstract:
The environment of a binary star system may contain two circumstellar disks, one orbiting each of the stars, and a circumbinary disk orbiting about the entire binary. The disk structure and evolution are modified by the presence of the binary. Resonances emit waves and open disk gaps. The binary's total mass and mass ratio as well as orbital elements can be modified by the disks. Signatures of the…
▽ More
The environment of a binary star system may contain two circumstellar disks, one orbiting each of the stars, and a circumbinary disk orbiting about the entire binary. The disk structure and evolution are modified by the presence of the binary. Resonances emit waves and open disk gaps. The binary's total mass and mass ratio as well as orbital elements can be modified by the disks. Signatures of these interactions provide observational tests of the dynamical models. The interaction of young planets with protoplanetary disks circularizes the orbits of Jupiter-mass planets and may produce much more massive extrasolar planets on eccentric orbits.
△ Less
Submitted 14 November, 2021;
originally announced November 2021.
-
An Analytic Solution to the Kozai-Lidov Evolution Equations
Authors:
Stephen H. Lubow
Abstract:
A test particle in a noncoplanar orbit about a member of a binary system can undergo Kozai-Lidov oscillations in which tilt and eccentricity are exchanged. An initially circular highly inclined particle orbit can reach high eccentricity. We consider the nonlinear secular evolution equations previously obtained in the quadrupole approximation. For the important case that the initial eccentricity of…
▽ More
A test particle in a noncoplanar orbit about a member of a binary system can undergo Kozai-Lidov oscillations in which tilt and eccentricity are exchanged. An initially circular highly inclined particle orbit can reach high eccentricity. We consider the nonlinear secular evolution equations previously obtained in the quadrupole approximation. For the important case that the initial eccentricity of the particle orbit is zero, we derive an analytic solution for the particle orbital elements as a function of time that is exact within the quadrupole approximation. The solution involves only simple trigonometric and hyperbolic functions. It simplifies in the case that the initial particle orbit is close to being perpendicular to the binary orbital plane. The solution also provides an accurate description of particle orbits with nonzero but sufficiently small initial eccentricity. It is accurate over a range of initial eccentricity that broadens at higher initial inclinations. In the case of an initial inclination of pi/3, an error of 1% at maximum eccentricity occurs for initial eccentricities of about 0.1.
△ Less
Submitted 20 July, 2021;
originally announced July 2021.
-
Sustained Kozai-Lidov oscillations in misaligned circumstellar gas disks
Authors:
Jeremy L. Smallwood,
Rebecca G. Martin,
Stephen H. Lubow
Abstract:
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 f…
▽ More
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.
△ Less
Submitted 20 December, 2020;
originally announced December 2020.
-
Improvements to Pan-STARRS1 Astrometry Using Gaia
Authors:
Stephen H. Lubow,
Richard L. White,
Bernie Shiao
Abstract:
We use the Gaia DR2 catalog to improve the astrometric accuracy of about 1.7 billion objects in Pan-STARRS1 Data Release 2 (PS1 DR2). We also obtain proper motions for these PS1 objects. The cross match between Gaia and PS1 reveals residuals that are correlated on a scale of about 1 arcmin. We apply a spatially adaptive correction algorithm for all PS1 objects having more than two detections to re…
▽ More
We use the Gaia DR2 catalog to improve the astrometric accuracy of about 1.7 billion objects in Pan-STARRS1 Data Release 2 (PS1 DR2). We also obtain proper motions for these PS1 objects. The cross match between Gaia and PS1 reveals residuals that are correlated on a scale of about 1 arcmin. We apply a spatially adaptive correction algorithm for all PS1 objects having more than two detections to reduce these residuals and align the object positions to Gaia. For point-like PS1 objects that cross match to Gaia, the algorithm reduces PS1/Gaia residuals by 33% in position (median value of 13.5 mas reduced to 9.0 mas) and by 24% in proper motion (median value of 6.3 mas/yr reduced to 4.8 mas/yr). The residuals for the corrected positions are smallest for objects with the most point-like morphologies and with intermediate magnitudes of about 17 mag. The residual errors in declination are systematically larger than those in right ascension; the declination errors increase with zenith angle in proportion to the air mass of the observations. Declination positional residuals at a given declination generally vary with color and are consistent with the effects of differential atmospheric refraction. In principle, these residuals could be reduced further by taking into account object color.
△ Less
Submitted 19 October, 2020;
originally announced October 2020.
-
Polar planets around highly eccentric binaries are the most stable
Authors:
Cheng Chen,
Stephen H. Lubow,
Rebecca G. Martin
Abstract:
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…
▽ More
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.
△ Less
Submitted 15 April, 2020;
originally announced April 2020.
-
Formation of the polar debris disc around 99 Herculis
Authors:
Jeremy L. Smallwood,
Alessia Franchini,
Cheng Chen,
Eric Becerril,
Stephen H. Lubow,
Chao-Chin Yang,
Rebecca G. Martin
Abstract:
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…
▽ More
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.
△ Less
Submitted 3 March, 2020;
originally announced March 2020.
-
Multi-planet disc interactions in binary systems
Authors:
Alessia Franchini,
Rebecca G. Martin,
Stephen H. Lubow
Abstract:
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 ge…
▽ More
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'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.
△ Less
Submitted 1 November, 2019;
originally announced November 2019.
-
Robust Registration of Astronomy Catalogs with Applications to the Hubble Space Telescope
Authors:
Fan Tian,
Tam�s Budav�ri,
Amitabh Basu,
Stephen H. Lubow,
Richard L. White
Abstract:
Astrometric calibration of images with a small field of view is often inferior to the internal accuracy of the source detections due to the small number of accessible guide stars. One important experiment with such challenges is the Hubble Space Telescope (HST). A possible solution is to cross-calibrate overlapping fields instead of just relying on standard stars. Following the approach of \citet{…
▽ More
Astrometric calibration of images with a small field of view is often inferior to the internal accuracy of the source detections due to the small number of accessible guide stars. One important experiment with such challenges is the Hubble Space Telescope (HST). A possible solution is to cross-calibrate overlapping fields instead of just relying on standard stars. Following the approach of \citet{2012ApJ...761..188B}, we use infinitesimal 3D rotations for fine-tuning the calibration but devise a better objective that is robust to a large number of false candidates in the initial set of associations. Using Bayesian statistics, we accommodate bad data by explicitly modeling the quality, which yields a formalism essentially identical to an $M$-estimation in robust statistics. Our results on simulated and real catalogs show great potentials for improving the HST calibration, and those with similar challenges.
△ Less
Submitted 28 August, 2019;
originally announced August 2019.
-
Orbital dynamics of circumbinary planets
Authors:
Cheng Chen,
Alessia Franchini,
Stephen H. Lubow,
Rebecca G. Martin
Abstract:
We investigate the dynamics of a nonzero mass, circular orbit planet around an eccentric orbit binary for various values of the binary eccentricity, binary mass fraction, planet mass, and planet semi--major axis by means of numerical simulations. Previous studies investigated the secular dynamics mainly by approximate analytic methods. In the stationary inclination state, the planet and binary pre…
▽ More
We investigate the dynamics of a nonzero mass, circular orbit planet around an eccentric orbit binary for various values of the binary eccentricity, binary mass fraction, planet mass, and planet semi--major axis by means of numerical simulations. Previous studies investigated the secular dynamics mainly by approximate analytic methods. In the stationary inclination state, the planet and binary precess together with no change in relative tilt. For both prograde and retrograde planetary orbits, we explore the conditions for planetary orbital libration versus circulation and the conditions for stationary inclination. As was predicted by analytic models, for sufficiently high initial inclination, a prograde planet's orbit librates about the stationary tilted state. For a fixed binary eccentricity, the stationary angle is a monotonically decreasing function of the ratio of the planet--to--binary angular momentum $j$. The larger $j$, the stronger the evolutionary changes in the binary eccentricity and inclination. We also calculate the critical tilt angle that separates the circulating from the librating orbits for both prograde and retrograde planet orbits. The properties of the librating orbits and stationary angles are quite different for prograde versus retrograde orbits. The results of the numerical simulations are in very good quantitative agreement with the analytic models. Our results have implications for circumbinary planet formation and evolution.
△ Less
Submitted 17 August, 2019;
originally announced August 2019.
-
Circumbinary disk inner radius as a diagnostic for disk--binary misalignment
Authors:
Alessia Franchini,
Stephen H. Lubow,
Rebecca G. Martin
Abstract:
We investigate the misalignment of the circumbinary disk around the binary HD 98800 BaBb with eccentricity $e\simeq 0.8$. \cite{Kennedy2019} observed the disk to be either at an inclination of $48^{\circ}$ or polar aligned to the binary orbital plane. Their simulations showed that alignment from $48^\circ$ to a polar configuration can take place on a shorter timescale than the age of this system.…
▽ More
We investigate the misalignment of the circumbinary disk around the binary HD 98800 BaBb with eccentricity $e\simeq 0.8$. \cite{Kennedy2019} observed the disk to be either at an inclination of $48^{\circ}$ or polar aligned to the binary orbital plane. Their simulations showed that alignment from $48^\circ$ to a polar configuration can take place on a shorter timescale than the age of this system. We perform hydrodynamical numerical simulations in order to estimate the cavity size carved by the eccentric binary for different disk inclinations as an independent check of polar alignment. Resonance theory suggests that torques on the inner parts of a polar disk around such a highly eccentric binary are much weaker than in the coplanar case, indicating a significantly smaller central cavity than in the coplanar case. We show that the inferred inner radius (from carbon monoxide measurements) of the accretion disk around BaBb can exclude the possibility of it being mildly inclined with respect to the binary orbital plane and therefore confirm the polar configuration. This study constitutes an important diagnostic for misaligned circumbinary disks, since it potentially allows us to infer the disk inclination from observed gas disk inner radii.
△ Less
Submitted 11 July, 2019;
originally announced July 2019.
-
Polar alignment of a protoplanetary disc around an eccentric binary III: Effect of disc mass
Authors:
Rebecca G. Martin,
Stephen H. Lubow
Abstract:
Martin & Lubow (2017) found that an initially sufficiently misaligned low mass protoplanetary disc around an eccentric binary undergoes damped nodal oscillations of tilt angle and longitude of ascending node. Dissipation causes evolution towards a stationary state of polar alignment in which the disc lies perpendicular to the binary orbital plane with angular momentum aligned to the eccentricity v…
▽ More
Martin & Lubow (2017) found that an initially sufficiently misaligned low mass protoplanetary disc around an eccentric binary undergoes damped nodal oscillations of tilt angle and longitude of ascending node. Dissipation causes evolution towards a stationary state of polar alignment in which the disc lies perpendicular to the binary orbital plane with angular momentum aligned to the eccentricity vector of the binary. We use hydrodynamic simulations and analytic methods to investigate how the mass of the disc affects this process. The simulations suggest that a disc with nonzero mass settles into a stationary state in the frame of the binary, the generalised polar state, at somewhat lower levels of misalignment with respect to the binary orbital plane, in agreement with the analytic model. Provided that discs settle into this generalised polar state, the observational determination of the misalignment angle and binary properties can be used to determine the mass of a circumbinary disc. We apply this constraint to the circumbinary disc in HD 98800. We obtain analytic criteria for polar alignment of a circumbinary ring with mass that approximately agree with the simulation results. Very broad misaligned discs undergo breaking, but the inner regions at least may still evolve to a polar state. The long term evolution of the disc depends on the evolution of the binary eccentricity that we find tends to decrease. Although the range of parameters required for polar alignment decreases somewhat with increasing disc mass, such alignment appears possible for a broad set of initial conditions expected in protostellar circumbinary discs.
△ Less
Submitted 24 September, 2019; v1 submitted 25 April, 2019;
originally announced April 2019.
-
Alignment of a circumbinary disc around an eccentric binary with application to KH 15D
Authors:
Jeremy L. Smallwood,
Stephen H. Lubow,
Alessia Franchini,
Rebecca G. Martin
Abstract:
We analyse the evolution of a mildly inclined circumbinary disc that orbits an eccentric orbit binary by means of smoother particle hydrodynamic (SPH) simulations and linear theory. We show that the alignment process of an initially misaligned circumbinary disc around an eccentric orbit binary is significantly different than around a circular orbit binary and involves tilt oscillations. The more e…
▽ More
We analyse the evolution of a mildly inclined circumbinary disc that orbits an eccentric orbit binary by means of smoother particle hydrodynamic (SPH) simulations and linear theory. We show that the alignment process of an initially misaligned circumbinary disc around an eccentric orbit binary is significantly different than around a circular orbit binary and involves tilt oscillations. The more eccentric the binary, the larger the tilt oscillations and the longer it takes to damp these oscillations. A circumbinary disc that is only mildly inclined may increase its inclination by a factor of a few before it moves towards alignment. The results of the SPH simulations agree well with those of linear theory. We investigate the properties of the circumbinary disc/ring around KH 15D. We determine disc properties based on the observational constraints imposed by the changing binary brightness. We find that the inclination is currently at a local minimum and will increase substantially before setting to coplanarity. In addition, the nodal precession is currently near its most rapid rate. The recent observations that show a reappearance of Star B impose constraints on the thickness of the layer of obscuring material. Our results suggest that disc solids have undergone substantial inward drift and settling towards to disc midplane. For disc masses $\sim 0.001 M_\odot$, our model indicates that the level of disc turbulence is low $α\ll 0.001$. Another possibility is that the disc/ring contains little gas.
△ Less
Submitted 1 April, 2019;
originally announced April 2019.
-
Generalised warped disk equations
Authors:
Rebecca G. Martin,
Stephen H. Lubow,
J. E. Pringle,
Alessia Franchini,
Zhaohuan Zhu,
Stephen Lepp,
Rebecca Nealon,
C. J. Nixon,
David Vallet
Abstract:
The manner in which warps in accretion disks evolve depends on the magnitude of the viscosity. For small viscosity $(α< H/R)$ the warp evolves in a wave-like manner; for large viscosity $H/R<α\ll 1$ it evolves diffusively. Here $α$ is the viscosity parameter and $H/R$ the disk aspect ratio. Currently there is no simple set of equations which describes the evolution in both regimes. In this paper w…
▽ More
The manner in which warps in accretion disks evolve depends on the magnitude of the viscosity. For small viscosity $(α< H/R)$ the warp evolves in a wave-like manner; for large viscosity $H/R<α\ll 1$ it evolves diffusively. Here $α$ is the viscosity parameter and $H/R$ the disk aspect ratio. Currently there is no simple set of equations which describes the evolution in both regimes. In this paper we describe a possible solution to this problem and introduce a set of one-dimensional equations that describe the evolution of a warped disk that are applicable in both high and low viscosity regimes for arbitrary tilts, but small warps.
△ Less
Submitted 28 February, 2019;
originally announced February 2019.
-
Misaligned accretion disc formation via Kozai-Lidov oscillations
Authors:
Alessia Franchini,
Rebecca G. Martin,
Stephen H. Lubow
Abstract:
We investigate the formation and evolution of misaligned accretion discs around the secondary component of a binary through mass transfer driven by Kozai-Lidov oscillations of the circumprimary disc's eccentricity and inclination. We perform SPH simulations to study the amount of mass transferred to the secondary star as a function of both the disc and binary parameters. For the range of parameter…
▽ More
We investigate the formation and evolution of misaligned accretion discs around the secondary component of a binary through mass transfer driven by Kozai-Lidov oscillations of the circumprimary disc's eccentricity and inclination. We perform SPH simulations to study the amount of mass transferred to the secondary star as a function of both the disc and binary parameters. For the range of parameters we explore, we find that increasing the disc aspect ratio, viscosity parameter and initial inclination as well as decreasing the binary mass ratio leads to larger amount of mass transfer, up to a maximum of about ten per cent of the initial mass of the primary disc. The circumsecondary disc forms with a high eccentricity and a high inclination and is also able to undergo KL oscillations. The circumsecondary disc oscillations have a shorter period than those in the disc around the primary. We find that some of the material that escapes the Roche-lobe of the two components forms a narrow misaligned circumbinary accretion disc. This study has implications for disc evolution in young binary star systems.
△ Less
Submitted 11 February, 2019;
originally announced February 2019.
-
White Dwarf Pollution by Asteroids from Secular Resonances
Authors:
Jeremy L. Smallwood,
Rebecca G. Martin,
Mario Livio,
Stephen H. Lubow
Abstract:
In the past few decades, observations have revealed signatures of metals polluting the atmospheres of white dwarfs. The diffusion timescale for metals to sink from the atmosphere of a white dwarf is of the order of days for a hydrogen-dominated atmosphere. Thus, there must be a continuous supply of metal-rich material accreting onto these white dwarfs. We investigate the role of secular resonances…
▽ More
In the past few decades, observations have revealed signatures of metals polluting the atmospheres of white dwarfs. The diffusion timescale for metals to sink from the atmosphere of a white dwarf is of the order of days for a hydrogen-dominated atmosphere. Thus, there must be a continuous supply of metal-rich material accreting onto these white dwarfs. We investigate the role of secular resonances that excite the eccentricity of asteroids allowing them to reach star-grazing orbits leading them to tidal disruption and the formation of a debris disc. Changes in the planetary system during the evolution of the star lead to a change in the location of secular resonances. In our Solar System, the engulfment of the Earth will cause the $ν_6$ resonance to shift outwards which will force previously stable asteroids to undergo secular resonant perturbations. With analytic models and $N$--body simulations we show that secular resonances driven by two outer companions can provide a source of continuous pollution. Secular resonances are a viable mechanism for the pollution of white dwarfs in a variety of exoplanetary system architectures.
△ Less
Submitted 7 July, 2018;
originally announced July 2018.
-
Polar alignment of a protoplanetary disc around an eccentric binary II: Effect of binary and disc parameters
Authors:
Rebecca G. Martin,
Stephen H. Lubow
Abstract:
In a recent paper Martin & Lubow showed that a circumbinary disc around an eccentric binary can undergo damped nodal oscillations that lead to the polar (perpendicular) alignment of the disc relative to the binary orbit. The disc angular momentum vector aligns to the eccentricity vector of the binary. We explore the robustness of this mechanism for a low mass disc (0.001 of the binary mass) and it…
▽ More
In a recent paper Martin & Lubow showed that a circumbinary disc around an eccentric binary can undergo damped nodal oscillations that lead to the polar (perpendicular) alignment of the disc relative to the binary orbit. The disc angular momentum vector aligns to the eccentricity vector of the binary. We explore the robustness of this mechanism for a low mass disc (0.001 of the binary mass) and its dependence on system parameters by means of hydrodynamic disc simulations. We describe how the evolution depends upon the disc viscosity, temperature, size, binary mass ratio, orbital eccentricity and inclination. We compare results with predictions of linear theory. We show that polar alignment of a low mass disc may occur over a wide range of binary-disc parameters. We discuss the application of our results to the formation of planetary systems around eccentric binary stars.
△ Less
Submitted 21 June, 2018;
originally announced June 2018.
-
Linear Analysis of the Evolution of Nearly Polar Low Mass Circumbinary Discs
Authors:
Stephen H. Lubow,
Rebecca G. Martin
Abstract:
Martin & Lubow (2017) showed through simulations that an initially tilted disc around an eccentric binary can evolve to polar alignment in which the disc lies perpendicular to the binary orbital plane. We apply linear theory to show both analytically and numerically that a nearly polar aligned low mass circumbinary disc evolves to polar alignment and determine the alignment timescale. Significant…
▽ More
Martin & Lubow (2017) showed through simulations that an initially tilted disc around an eccentric binary can evolve to polar alignment in which the disc lies perpendicular to the binary orbital plane. We apply linear theory to show both analytically and numerically that a nearly polar aligned low mass circumbinary disc evolves to polar alignment and determine the alignment timescale. Significant disc evolution towards the polar state around moderately eccentric binaries can occur for typical protostellar disc parameters in less than a typical disc lifetime for binaries with orbital periods of order 100 years or less. Resonant torques are much less effective at truncating the inner parts of circumbinary polar discs than the inner parts of coplanar discs. For polar discs, they vanish for a binary eccentricity of unity. The results agree with the simulations in showing that discs can evolve to a polar state. Circumbinary planets may then form in such discs and reside on polar orbits.
△ Less
Submitted 5 October, 2017;
originally announced October 2017.
-
Kozai-Lidov Disc Instability
Authors:
Stephen H. Lubow,
Gordon I. Ogilvie
Abstract:
Recent results by Martin et al. (2014) showed in 3D SPH simulations that tilted discs in binary systems can be unstable to the development of global, damped Kozai-Lidov (KL) oscillations in which the discs exchange tilt for eccentricity. We investigate the linear stability of KL modes for tilted inviscid discs under the approximations that the disc eccentricity is small and the disc remains flat.…
▽ More
Recent results by Martin et al. (2014) showed in 3D SPH simulations that tilted discs in binary systems can be unstable to the development of global, damped Kozai-Lidov (KL) oscillations in which the discs exchange tilt for eccentricity. We investigate the linear stability of KL modes for tilted inviscid discs under the approximations that the disc eccentricity is small and the disc remains flat. By using 1D equations, we are able to probe regimes of large ratios of outer to inner disc edge radii that are realistic for binary systems of hundreds of AU separations and are not easily probed by multi-dimensional simulations. For order unity binary mass ratios, KL instability is possible for a window of disc aspect ratios H/r in the outer parts of a disc that roughly scale as (n_b/n)^2 < H/r < n_b/n, for binary orbital frequency n_b and orbital frequency n at the disc outer edge. We present a framework for understanding the zones of instability based on the determination of branches of marginally unstable modes. In general, multiple growing eccentric KL modes can be present in a disc. Coplanar apsidal-nodal precession resonances delineate instability branches. We determine the range of tilt angles for unstable modes as a function of disc aspect ratio. Unlike the KL instability for free particles that involves a critical (minimum) tilt angle, disc instability is possible for any nonzero tilt angle depending on the disc aspect ratio.
△ Less
Submitted 24 April, 2017;
originally announced April 2017.
-
Polar alignment of a protoplanetary disk around an eccentric binary
Authors:
Rebecca G. Martin,
Stephen H. Lubow
Abstract:
We use three-dimensional hydrodynamical simulations to show that an initially mildly misaligned circumbinary accretion disk around an eccentric binary can evolve to an orientation that is perpendicular to the orbital plane of the binary (polar alignment). As the disk evolves to the perpendicular state, it undergoes nodal libration oscillations of the tilt angle and the longitude of the ascending n…
▽ More
We use three-dimensional hydrodynamical simulations to show that an initially mildly misaligned circumbinary accretion disk around an eccentric binary can evolve to an orientation that is perpendicular to the orbital plane of the binary (polar alignment). As the disk evolves to the perpendicular state, it undergoes nodal libration oscillations of the tilt angle and the longitude of the ascending node. Dissipation within the disk causes the oscillations to damp. The process operates above a critical initial misalignment angle that depends upon the eccentricity of the binary and the mass of the disk. For binary eccentricity of 0.5, the process operates typically for disk masses smaller than a few percent of the binary mass and initial tilt angle of more than 40 degrees. This evolution has important implications for planet formation around eccentric binary star systems.
△ Less
Submitted 2 February, 2017;
originally announced February 2017.
-
Fragmentation of Kozai-Lidov Disks
Authors:
Wen Fu,
Stephen H. Lubow,
Rebecca G. Martin
Abstract:
We analyze the gravitational instability (GI) of a locally isothermal inclined disk around one component of a binary system. Such a disk can undergo global Kozai-Lidov (KL) cycles if the initial disk tilt is above the critical KL angle (of about 40 degrees). During these cycles, an initially circular disk exchanges its inclination for eccentricity, and vice versa. Self-gravity may suppress the cyc…
▽ More
We analyze the gravitational instability (GI) of a locally isothermal inclined disk around one component of a binary system. Such a disk can undergo global Kozai-Lidov (KL) cycles if the initial disk tilt is above the critical KL angle (of about 40 degrees). During these cycles, an initially circular disk exchanges its inclination for eccentricity, and vice versa. Self-gravity may suppress the cycles under some circumstances. However, with hydrodynamic simulations including self-gravity we show that for a sufficiently high initial disk tilts and for certain disk masses, disks can undergo KL oscillations and fragment due to GI, even when the Toomre Q value for an equivalent undisturbed disk is well within the stable regime (Q > 2). We suggest that KL triggered disk fragmentation provides a mechanism for the efficient formation of giant planets in binary systems and may enhance fragmentation of disks in massive black hole binaries.
△ Less
Submitted 22 December, 2016;
originally announced December 2016.
-
Turbulence, Transport and Waves in Ohmic Dead Zones
Authors:
Daniel Gole,
Jacob B. Simon,
Stephen H. Lubow,
Philip J. Armitage
Abstract:
We use local numerical simulations to study a vertically stratified accretion disk with a resistive mid-plane that damps magnetohydrodynamic (MHD) turbulence. This is an idealized model for the dead zones that may be present at some radii in protoplanetary and dwarf novae disks. We vary the relative thickness of the dead and active zones to quantify how forced fluid motions in the dead zone change…
▽ More
We use local numerical simulations to study a vertically stratified accretion disk with a resistive mid-plane that damps magnetohydrodynamic (MHD) turbulence. This is an idealized model for the dead zones that may be present at some radii in protoplanetary and dwarf novae disks. We vary the relative thickness of the dead and active zones to quantify how forced fluid motions in the dead zone change. We find that the residual Reynolds stress near the mid-plane decreases with increasing dead zone thickness, becoming negligible in cases where the active to dead mass ratio is less than a few percent. This implies that purely Ohmic dead zones would be vulnerable to episodic accretion outbursts via the mechanism of Martin & Lubow (2011). We show that even thick dead zones support a large amount of kinetic energy, but this energy is largely in fluid motions that are inefficient at angular momentum transport. Confirming results from Oishi & Mac Low (2009), the perturbed velocity field in the dead zone is dominated by an oscillatory, vertically extended circulation pattern with a low frequency compared to the orbital frequency. This disturbance has the properties predicted for the lowest order r mode in a hydrodynamic disk. We suggest that in a global disk similar excitations would lead to propagating waves, whose properties would vary with the thickness of the dead zone and the nature of the perturbations (isothermal or adiabatic). Flows with similar amplitudes would buckle settled particle layers and could reduce the efficiency of pebble accretion.
△ Less
Submitted 9 June, 2016;
originally announced June 2016.
-
Planet-disc evolution and the formation of Kozai-Lidov planets
Authors:
Rebecca G. Martin,
Stephen H. Lubow,
Chris Nixon,
Philip J. Armitage
Abstract:
With hydrodynamical simulations we determine the conditions under which an initially coplanar planet-disc system that orbits a member of a misaligned binary star evolves to form a planet that undergoes Kozai-Lidov (KL) oscillations once the disc disperses. These oscillations may explain the large orbital eccentricities, as well as the large misalignments with respect to the spin of the central sta…
▽ More
With hydrodynamical simulations we determine the conditions under which an initially coplanar planet-disc system that orbits a member of a misaligned binary star evolves to form a planet that undergoes Kozai-Lidov (KL) oscillations once the disc disperses. These oscillations may explain the large orbital eccentricities, as well as the large misalignments with respect to the spin of the central star, observed for some exoplanets. The planet is assumed to be massive enough to open a gap in the disc. The planet's tilt relative to the binary orbital plane is subject to two types of oscillations. The first type, present at even small inclination angles relative to the binary orbital plane, is due to the interaction of the planet with the disc and binary companion and is amplified by a secular resonance. The second type of oscillation is the KL oscillation that operates on both the planet and disc at larger binary inclination angles. We find that for a sufficiently massive disc, even a relatively low inclination planet-disc system can force a planet to an inclination above the critical KL angle, as a consequence of the first type of tilt oscillation, allowing it to undergo the second type of oscillation. We conclude that the hydrodynamical evolution of a sufficiently massive and inclined disc in a binary system broadens the range of systems that form eccentric and misaligned giant planets to include a wide range of initial misalignment angles (20 to 160 degrees).
△ Less
Submitted 9 March, 2016;
originally announced March 2016.
-
Version 1 of the Hubble Source Catalog
Authors:
Bradley C. Whitmore,
Sahar S. Allam,
Tamas Budavari,
Stefano Casertano,
Ronald A. Downes,
Thomas Donaldson,
S. Michael Fall,
Stephen H. Lubow,
Lee Quick,
Louis-Gregory Strolger,
Geoff Wallace,
Richard L. White
Abstract:
The Hubble Source Catalog is designed to help optimize science from the Hubble Space Telescope by combining the tens of thousands of visit-based source lists in the Hubble Legacy Archive into a single master catalog. Version 1 of the Hubble Source Catalog includes WFPC2, ACS/WFC, WFC3/UVIS, and WFC3/IR photometric data generated using SExtractor software to produce the individual source lists. The…
▽ More
The Hubble Source Catalog is designed to help optimize science from the Hubble Space Telescope by combining the tens of thousands of visit-based source lists in the Hubble Legacy Archive into a single master catalog. Version 1 of the Hubble Source Catalog includes WFPC2, ACS/WFC, WFC3/UVIS, and WFC3/IR photometric data generated using SExtractor software to produce the individual source lists. The catalog includes roughly 80 million detections of 30 million objects involving 112 different detector/filter combinations, and about 160 thousand HST exposures. Source lists from Data Release 8 of the Hubble Legacy Archive are matched using an algorithm developed by Budavari & Lubow (2012). The mean photometric accuracy for the catalog as a whole is better than 0.10 mag, with relative accuracy as good as 0.02 mag in certain circumstances (e.g., bright isolated stars). The relative astrometric residuals are typically within 10 mas, with a value for the mode (i.e., most common value) of 2.3 mas. The absolute astrometric accuracy is better than $\sim$0.1 arcsec for most sources, but can be much larger for a fraction of fields that could not be matched to the PanSTARRS, SDSS, or 2MASS reference systems. In this paper we describe the database design with emphasis on those aspects that enable the users to fully exploit the catalog while avoiding common misunderstandings and potential pitfalls. We provide usage examples to illustrate some of the science capabilities and data quality characteristics, and briefly discuss plans for future improvements to the Hubble Source Catalog.
△ Less
Submitted 15 February, 2016;
originally announced February 2016.
-
The Evolution of Planet-Disk Systems That Are Mildly Inclined to the Orbit of a Binary Companion
Authors:
Stephen H. Lubow,
Rebecca G. Martin
Abstract:
We determine the evolution of a giant planet-disk system that orbits a member of a binary star system and is mildly inclined with respect to the binary orbital plane. The planet orbit and disk are initially mutually coplanar. We analyze the evolution of the planet and the disk by analytic means and hydrodynamic simulations. We generally find that the planet and the disk do not remain coplanar unle…
▽ More
We determine the evolution of a giant planet-disk system that orbits a member of a binary star system and is mildly inclined with respect to the binary orbital plane. The planet orbit and disk are initially mutually coplanar. We analyze the evolution of the planet and the disk by analytic means and hydrodynamic simulations. We generally find that the planet and the disk do not remain coplanar unless the disk mass is very large or the gap that separates the planet from the disk is very small. The relative planet-disk tilt undergoes secular oscillations whose initial amplitudes are typically of order the initial disk tilt relative to the binary orbital plane for disk masses ~1% of the binary mass or less. The effects of a secular resonance and the disk tilt decay enhance the planet-disk misalignment. The secular resonance plays an important role for disk masses greater than the planet mass. At later times, the accretion of disk gas by the planet causes its orbit to evolve towards alignment, if the disk mass is sufficiently large. The results have several implications for the evolution of massive planets in binary systems.
△ Less
Submitted 7 December, 2015;
originally announced December 2015.
-
The Kozai-Lidov Mechanism in Hydrodynamical Disks. III. Effects of Disk Mass and Self-Gravity
Authors:
Wen Fu,
Stephen H. Lubow,
Rebecca G. Martin
Abstract:
Previously we showed that a substantially misaligned viscous accretion disk with pressure that orbits around one component of a binary system can undergo global damped Kozai-Lidov (KL) oscillations. These oscillations produce periodic exchanges of the disk eccentricity with inclination. The disk KL mechanism is quite robust and operates over a wide range of binary and disk parameters. However, the…
▽ More
Previously we showed that a substantially misaligned viscous accretion disk with pressure that orbits around one component of a binary system can undergo global damped Kozai-Lidov (KL) oscillations. These oscillations produce periodic exchanges of the disk eccentricity with inclination. The disk KL mechanism is quite robust and operates over a wide range of binary and disk parameters. However, the effects of self-gravity, which are expected to suppress the KL oscillations for sufficiently massive disks, were ignored. Here, we analyze the effects of disk self-gravity by means of hydrodynamic simulations and compare the results with the expectations of analytic theory. The disk mass required for suppression in the simulations is a few percent of the mass of the central star and this roughly agrees with an analytical estimate. The conditions for suppression of the KL oscillations in the simulations are close to requiring that the disk be gravitationally unstable. We discuss some implications of our results for the dynamics of protoplanetary disks and the related planet formation.
△ Less
Submitted 9 November, 2015; v1 submitted 3 September, 2015;
originally announced September 2015.
-
The Kozai-Lidov Mechanism in Hydrodynamical Disks. II. Effects of binary and disk parameters
Authors:
Wen Fu,
Stephen H. Lubow,
Rebecca G. Martin
Abstract:
Martin et al. (2014b) showed that a substantially misaligned accretion disk around one component of a binary system can undergo global damped Kozai-Lidov oscillations. During these oscillations, the inclination and eccentricity of the disk are periodically exchanged. However, the robustness of this mechanism and its dependence on the system parameters were unexplored. In this paper, we use three-d…
▽ More
Martin et al. (2014b) showed that a substantially misaligned accretion disk around one component of a binary system can undergo global damped Kozai-Lidov oscillations. During these oscillations, the inclination and eccentricity of the disk are periodically exchanged. However, the robustness of this mechanism and its dependence on the system parameters were unexplored. In this paper, we use three-dimensional hydrodynamical simulations to analyze how various binary and disk parameters affect the Kozai-Lidov mechanism in hydrodynamical disks. The simulations include the effect of gas pressure and viscosity, but ignore the effects of disk self-gravity. We describe results for different numerical resolutions, binary mass ratios and orbital eccentricities, initial disk sizes, initial disk surface density profiles, disk sound speeds, and disk viscosities. We show that the Kozai-Lidov mechanism can operate for a wide range of binary-disk parameters. We discuss the applications of our results to astrophysical disks in various accreting systems.
△ Less
Submitted 2 July, 2015; v1 submitted 1 April, 2015;
originally announced April 2015.
-
Resonances in retrograde circumbinary discs
Authors:
Chris Nixon,
Stephen H. Lubow
Abstract:
We analyse the interaction of an eccentric binary with a circular coplanar circumbinary disc that rotates in a retrograde sense with respect to the binary. In the circular binary case, no Lindblad resonances lie within the disc and no Lindblad resonant torques are produced, as was previously known. By analytic means, we show that when the binary orbit is eccentric, there exist components of the gr…
▽ More
We analyse the interaction of an eccentric binary with a circular coplanar circumbinary disc that rotates in a retrograde sense with respect to the binary. In the circular binary case, no Lindblad resonances lie within the disc and no Lindblad resonant torques are produced, as was previously known. By analytic means, we show that when the binary orbit is eccentric, there exist components of the gravitational potential of the binary which rotate in a retrograde sense to the binary orbit and so rotate progradely with respect to this disc, allowing a resonant interaction to occur between the binary and the disc. The resulting resonant torques distinctly alter the disc response from the circular binary case. We describe results of three-dimensional hydrodynamic simulations to explore this effect and categorise the response of the disc in terms of modes whose strengths vary as a function of binary mass ratio and eccentricity. These mode strengths are weak compared to the largest mode strengths expected in the prograde case where the binary and disc rotate in the same sense. However, for sufficiently high binary eccentricity, resonant torques open a gap in a retrograde circumbinary disc, while permitting gas inflow on to the binary via gas streams. The inflow results in a time varying accretion rate on to the binary that is modulated over the binary orbital period, as was previously found to occur in the prograde case.
△ Less
Submitted 28 January, 2015;
originally announced January 2015.
-
The Inner Disk Structure, Disk-Planet Interactions, and Temporal Evolution in the Beta Pictoris System: A Two-Epoch HST/STIS Coronagraphic Study
Authors:
Daniel Apai,
Glenn Schneider,
Carol A. Grady,
Mark C. Wyatt,
Anne-Marie Lagrange,
Marc J. Kuchner,
Christopher J. Stark,
Stephen H. Lubow
Abstract:
We present deep HST/STIS coronagraphic images of the Beta Pic debris disk obtained at two epochs separated by 15 years. The new images and the re-reduction of the 1997 data provide the most sensitive and detailed views of the disk at optical wavelengths as well as the yet smallest inner working angle optical coronagraphic image of the disk. Our observations characterize the large-scale and inner-d…
▽ More
We present deep HST/STIS coronagraphic images of the Beta Pic debris disk obtained at two epochs separated by 15 years. The new images and the re-reduction of the 1997 data provide the most sensitive and detailed views of the disk at optical wavelengths as well as the yet smallest inner working angle optical coronagraphic image of the disk. Our observations characterize the large-scale and inner-disk asymmetries and we identify multiple breaks in the disk radial surface brightness profile. We study in detail the radial and vertical disk structure and show that the disk is warped. We explore the disk at the location of the Beta Pic b super-jupiter and find that the disk surface brightness slope is continuous between 0.5 and 2.0 arcsec, arguing for no change at the separations where Beta Pic b orbits. The two epoch images constrain the disk surface brightness evolution on orbital and radiation pressure blow-out timescales. We place an upper limit of 3% on the disk surface brightness change between 3-5 arcsec, including the locations of the disk warp, and the CO and dust clumps. We discuss the new observations in the context of high-resolution multi-wavelength images and divide the disk asymmetries in two groups: axisymmetric and non-axisymmetric. The axisymmetric structures (warp, large-scale butterfly, etc.) are consistent with disk structure models that include interactions of a planetesimal belt and a non-coplanar giant planet. The non-axisymmetric features, however, require a different explanation.
△ Less
Submitted 13 January, 2015;
originally announced January 2015.
-
Tidal Torques on Misaligned Disks in Binary Systems
Authors:
Stephen H. Lubow,
Rebecca G. Martin,
Chris Nixon
Abstract:
We extend previous studies of the tidal truncation of coplanar disks in binary systems to the more general case of noncoplanar disks. As in the prograde coplanar case, Lindblad resonances play a key role in tidal truncation. We analyze the tidal torque acting on a misaligned nearly circular disk in a circular orbit binary system. We concentrate on the 2:1 inner Lindblad resonance associated with t…
▽ More
We extend previous studies of the tidal truncation of coplanar disks in binary systems to the more general case of noncoplanar disks. As in the prograde coplanar case, Lindblad resonances play a key role in tidal truncation. We analyze the tidal torque acting on a misaligned nearly circular disk in a circular orbit binary system. We concentrate on the 2:1 inner Lindblad resonance associated with the m=2 tidal forcing (for azimuthal wavenumber m) that plays a major role in the usual coplanar case. We determine the inclination dependence of this torque, which is approximately cos^8(i/2) for misalignment angle i. Compared to the prograde coplanar case (i=0), this torque decreases by a factor of about 2 for i = pi/6 and by a factor of about 20 for i=pi/2. The Lindblad torque decreases to zero for a tilt angle of pi (counter-rotation), consistent with previous investigations. The effects of higher order resonances associated with m>2 tidal forcing may contribute somewhat, but are much more limited than in the i=0 case. These results suggest that misaligned disks in binary systems can be significantly extended compared to their coplanar counterparts. In cases where a disk is sufficiently inclined and viscous, it can overrun all Lindblad resonances and overflow the Roche lobe of the disk central object.
△ Less
Submitted 24 December, 2014;
originally announced December 2014.