Source author record

Rebecca G. Martin

Rebecca G. Martin appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

47works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

47 published item(s)

preprint2022arXiv

A radial limit on polar circumbinary orbits from general relativity

A particle orbiting a misaligned eccentric orbit binary undergoes nodal precession either around the binary angular momentum vector (a circulating orbit) or around a stationary inclination (a librating orbit). In the absence of general relativity, the stationary inclination is inclined by 90 degrees to the binary angular momentum vector (aligned with the binary eccentricity vector) and does not depend on the particle semi-major axis. General relativity causes apsidal precession of the binary orbit. Close to the binary, the behaviour of the particle is not significantly affected, a librating particle precesses with the binary. However, we find that the stationary inclination and the minimum inclination required for libration both increase with the particle semi-major axis. There is a critical radius beyond which there are no librating orbits, only circulating orbits, and therefore there is a maximum orbital radius for a stationary polar orbiting body. The critical radius is within planet forming regions around binaries with semi-major axis <= 1 au. This has implications for the search for misaligned circumbinary planets and the radial extent of polar circumbinary disks.

preprint2022arXiv

Accretion onto a binary from a polar circumbinary disc

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

preprint2022arXiv

Asteroids and life: How special is the solar system?

Asteroid impacts with the Earth may have played an essential role in the emergence of life on Earth through their creation of favorable niches for life, changes to the atmosphere and delivery of water. Consequently, we suggest two potential requirements for life in an exoplanetary system: first, that the system has an asteroid belt, and second, that there is a mechanism to drive asteroids to impact the terrestrial habitable planet. Since in the solar system, the $ν_6$ secular resonance has been shown to have been important in driving these impacts, we explore how the masses and locations of two giant planets determine the location and strength of this secular resonance. Examining observed exoplanetary systems with two giant planets, we find that a secular resonance within the asteroid belt region may not be uncommon. Hence the solar system is somewhat special, but the degree of fine-tuning that may be necessary for the emergence of life is not excessive. Finally, with $n$-body simulations, we show that when the two giant planets are close to the 2:1 mean motion resonance, the asteroid belt is unstable but this does not lead to increased asteroid delivery.

preprint2022arXiv

Circumbinary disk evolution in the presence of an outer companion star

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

preprint2022arXiv

Eccentric dust ring formation in Kozai-Lidov gas disks

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

preprint2022arXiv

Fast nodal precession of the disc around Pleione requires a broken disc

Pleione is a Be star that is in a 218 day orbit with a low-mass binary companion. Recent numerical simulations have shown that a Be star disc can be subject to breaking when material is actively being fed into the inner parts of the disc. After breaking, the disc is composed of two rings: an inner ring that is anchored to the stellar equator and an outer ring that is free to nodally precess. A double ring disc may explain some of the observed variability in Pleione. We model the nodal precession of the outer disc ring that is driven by the companion on an observed timescale of $80.5\,\rm yr$. We find that the outer ring of a broken disc in a binary with an eccentricity of $e_{\rm b}= 0.6$ can precess on the observed timescale and have an outer radius that is in rough agreement with the observed disc size. An unbroken disc model cannot fit both the observed precession rate and disc size. Suppression of Kozai-Lidov driven disc eccentricity is more likely for a high binary eccentricity if the disc extends to the tidal truncation radius.

preprint2022arXiv

Misaligned circumbinary disks as efficient progenitors of interstellar asteroids

Gaseous circumbinary disks (CBDs) that are highly inclined to the binary orbit are commonly observed in nature. These disks harbor particles that can reach large mutual inclinations as a result of nodal precession once the gas disk has dissipated. With n-body simulations that include fragmentation we demonstrate that misaligned disks of particles can be efficient progenitors of interstellar asteroids (ISAs). Collisions that take place between particles with large mutual inclinations have large impact velocities which can result in mass ejection, with a wide range of fragment sizes and ejection velocities. We explore the binary parameters for which the majority of the terrestrial planet forming material is ejected rather than accreted into planets. The misalignment required to eject significant material decreases with binary eccentricity. If the distribution of binary eccentricity is uniform and the initial particle CBD orientation relative to the binary orbit is isotropic, about 59% of binaries are more likely to eject the majority of their CBD terrestrial planet disk mass through high velocity body-body collisions rather than retain this material and build terrestrial planets. However, binary--disk interactions during the gas disk phase with non-zero disk viscosity will reduce this fraction. The composition, small size, highly elongated shape, and tumbling motion of `Oumuamua is consistent with ISAs generated by misaligned CBDs.

preprint2022arXiv

Misalignment of terrestrial circumbinary planets as an indicator of their formation mechanism

Circumbinary gas disks are often observed to be misaligned to the binary orbit suggesting that planet formation may proceed in a misaligned disk. With N-body simulations we consider the formation of circumbinary terrestrial planets from a particle disk that is initially misaligned. We find that if terrestrial planets form in this way, in the absence of gas, they can only form close to coplanar or close to polar to the binary orbit. Planets around a circular binary form coplanar while planets around an eccentric binary can form coplanar or polar depending on the initial disk misalignment and the binary eccentricity. The more massive a terrestrial planet is, the more aligned it is (to coplanar or polar) since it has undergone more mergers that lead on average to smaller misalignment angles. Nodal precession of particle disks with very large initial inclinations lead to high mutual inclinations between the particles. This produces high relative velocities between particles that leads to mass ejections that can completely inhibit planet formation. Misaligned terrestrial circumbinary planets may be able to form in the presence of a misaligned circumbinary gas disk that may help to nodally align the particle orbits and maintain the inclination of the planets during their formation.

preprint2021arXiv

Kozai-Lidov oscillations triggered by a tilt instability of detached circumplanetary discs

Circumplanetary discs can be linearly unstable to the growth of disc tilt in the tidal potential of the star-planet system. We use three-dimensional hydrodynamical simulations to characterize the disc conditions needed for instability, together with its long term evolution. Tilt growth occurs for disc aspect ratios, evaluated near the disc outer edge, of $H/r\gtrsim 0.05$, with a weak dependence on viscosity in the wave-like regime of warp propagation. Lower mass giant planets are more likely to have circumplanetary discs that satisfy the conditions for instability. We show that the tilt instability can excite the inclination to above the threshold where the circumplanetary disc becomes unstable to Kozai--Lidov (KL) oscillations. Dissipation in the Kozai--Lidov unstable regime caps further tilt growth, but the disc experiences large oscillations in both inclination and eccentricity. Planetary accretion occurs in episodic accretion events. We discuss implications of the joint tilt--KL instability for the detectability of circumplanetary discs, for the obliquity evolution of forming giant planets, and for the formation of satellite systems.

preprint2020arXiv

A fast-growing tilt instability of detached circumplanetary disks

Accretion disks in binary systems can exhibit a tilt instability, arising from the interaction between components of the tidal potential and dissipation. Using a linear analysis, we show that the aspect ratios and outer radii of circumplanetary disks provide favorable conditions for tilt growth. We quantify the growth rate of the instability using particle-based ({\sc phantom}) and grid-based ({\sc athena++}) hydrodynamic simulations. For a disk with outer aspect ratio $H/r \simeq 0.1$, initially moderate tilts double on a time scale of about 15-30 binary orbits. Our results imply that detached circumplanetary disks, whose evolution is not entirely controlled by accretion from the circumstellar disk, may commonly be misaligned to the planetary orbital plane. We discuss implications for planetary spin evolution, and possible interactions between the tilt instability and Kozai-Lidov dynamics.

preprint2020arXiv

Asteroid belt survival through stellar evolution: dependence on the stellar mass

Polluted white dwarfs are generally accreting terrestrial-like material that may originate from a debris belt like the asteroid belt in the solar system. The fraction of white dwarfs that are polluted drops off significantly for white dwarfs with masses $M_{\rm WD}\gtrsim 0.8\,\rm M_\odot$. This implies that asteroid belts and planetary systems around main-sequence stars with mass $M_{\rm MS}\gtrsim 3\,\rm M_\odot$ may not form because of the intense radiation from the star. This is in agreement with current debris disc and exoplanet observations. The fraction of white dwarfs that show pollution also drops off significantly for low mass white dwarfs $(M_{\rm WD}\lesssim 0.55\,\rm M_\odot)$. However, the low-mass white dwarfs that do show pollution are not currently accreting but have accreted in the past. We suggest that asteroid belts around main sequence stars with masses $M_{\rm MS}\lesssim 2\,\rm M_\odot$ are not likely to survive the stellar evolution process. The destruction likely occurs during the AGB phase and could be the result of interactions of the asteroids with the stellar wind, the high radiation or, for the lowest mass stars that have an unusually close-in asteroid belt, scattering during the tidal orbital decay of the inner planetary system.

preprint2020arXiv

Evolution of $α$ Centauri B's protoplanetary disc

With hydrodynamical simulations we examine the evolution of a protoplanetary disc around $α$ Centauri B including the effect of the eccentric orbit binary companion $α$ Centauri A. The initially circular orbit disc undergoes two types of eccentricity growth. First, the eccentricity oscillates on the orbital period of the binary, $P_{\rm orb}$, due to the eccentricity of the binary orbit. Secondly, for a sufficiently small disc aspect ratio, the disc undergoes global forced eccentricity oscillations on a time-scale of around $20\,P_{\rm orb}$. These oscillations damp out through viscous dissipation leaving a quasi-steady eccentricity profile for the disc that oscillates only on the binary orbital period. The time-averaged global eccentricity is in the range 0.05-0.1, with no precession in the steady state. The periastrons of the gas particles are aligned to one another. The higher the disc viscosity, the higher the disc eccentricity. With $N$-body simulations we examine the evolution of a disc of planetesimals that forms with the orbital properties of the quasi-steady protoplanetary disc. We find that the average magnitude of the eccentricity of particles increases and their periastrons become misaligned to each other once they decouple from the gas disc. The low planetesimal collision velocity required for planet formation suggests that for planet formation to have occurred in a disc of planetesimals formed from a protoplanetary disc around $α$ Centauri B, said disc's viscosity must be have been small and planet formation must have occurred at orbital radii smaller than about $2.5\,\rm au$. Planet formation may be easier with the presence of gas.

preprint2020arXiv

Formation of the polar debris disc around 99 Herculis

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

preprint2020arXiv

GW Ori: Interactions Between a Triple-star System and its Circumtriple Disk in Action

GW Ori is a hierarchical triple system which has a rare circumtriple disk. We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of 1.3 mm dust continuum and 12CO J=2-1 molecular gas emission of the disk. For the first time, we identify three dust rings in the disk at ~46, 188, and 338 AU, with estimated dust mass of ~70-250 Earth masses, respectively. To our knowledge, the outer ring in GW Ori is the largest dust ring ever found in protoplanetary disks. We use visibility modelling of dust continuum to show that the disk has misaligned parts and the innermost dust ring is eccentric. The disk misalignment is also suggested by the CO kinematics modelling. We interpret these substructures as evidence of ongoing dynamical interactions between the triple stars and the circumtriple disk.

preprint2020arXiv

Polar planets around highly eccentric binaries are the most stable

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

preprint2020arXiv

Sustained Kozai-Lidov oscillations in misaligned circumstellar gas disks

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

preprint2019arXiv

Multi-planet disc interactions in binary systems

We investigate the evolution of a multi--planet--disc system orbiting one component of a binary star system. The planet--disc system is initially coplanar but misaligned to the binary orbital plane. The planets are assumed to be giants that open gaps in the disc. We first study the role of the disc in shaping the mutual evolution of the two planets using a secular model for low initial tilt. In general we find that the planets and the disc do not remain coplanar, in agreement with Lubow & Martin (2016) for the single planet case. Instead, the planets and the disc undergo tilt oscillations. A high mass disc between the two planets causes the planets and the disc to nodally precess at the same average rate but they are generally misaligned. The amplitude of the tilt oscillations between the planets is larger while the disc is present. We then consider higher initial tilts using hydrodynamical simulations and explore the possibility of the formation of eccentric Kozai-Lidov (KL) planets. We find that the inner planet'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.

preprint2016arXiv

On the Formation of Super-Earths with Implications for the Solar System

We first consider how the level of turbulence in a protoplanetary disk affects the formation locations for the observed close-in super-Earths in exosolar systems. We find that a protoplanetary disk that includes a dead zone (a region of low turbulence) has substantially more material in the inner parts of the disk, possibly allowing for in situ formation. For the dead zone to last the entire lifetime of the disk requires the active layer surface density to be sufficiently small, <100 g/cm^2. Migration through a dead zone may be very slow and thus super-Earth formation followed by migration towards the star through the dead zone is less likely. For fully turbulent disks, there is not enough material for in situ formation. However, in this case, super-Earths can form farther out in the disk and migrate inwards on a reasonable timescale. We suggest that both of these formation mechanisms operate in different planetary systems. This can help to explain the observed large range in densities of super-Earths because the formation location determines the composition. Furthermore, we speculate that super-Earths could have formed in the inner parts of our solar system and cleared the material in the region inside of Mercury's orbit. The super-Earths could migrate through the gas disk and fall into the Sun if the disk was sufficiently cool during the final gas disk accretion process. While it is definitely possible to meet all of these requirements, we don't expect them to occur in all systems, which may explain why the solar system is somewhat special in its lack of super-Earths.

preprint2016arXiv

Planet-disc evolution and the formation of Kozai-Lidov planets

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).

preprint2015arXiv

Nebular dead zone effects on the D/H ratio in chondrites and comets

Comets and chondrites show non-monotonic behaviour of their Deuterium to Hydrogen (D/H) ratio as a function of their formation location from the Sun. This is difficult to explain with a classical protoplanetary disk model that has a decreasing temperature structure with radius from the Sun. We want to understand if a protoplanetary disc with a dead zone, a region of zero or low turbulence, can explain the measured D/H values in comets and chondrites. We use time snapshots of a vertically layered disk model with turbulent surface layers and a dead zone at the midplane. The disc has a non-monotonic temperature structure due to increased heating from self-gravity in the outer parts of the dead zone. We couple this to a D/H ratio evolution model in order to quantify the effect of such thermal profiles on D/H enrichment in the nebula. We find that the local temperature peak in the disk can explain the diversity in the D/H ratios of different chondritic families. This disk temperature profile leads to a non-monotonic D/H enrichment evolution, allowing these families to acquire their different D/H values while forming in close proximity. The formation order we infer for these families is compatible with that inferred from their water abundances. However, we find that even for very young disks, the thermal profile reversal is too close to the Sun to be relevant for comets.

preprint2015arXiv

The Evolution of Planet-Disk Systems That Are Mildly Inclined to the Orbit of a Binary Companion

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.

preprint2015arXiv

The Kozai-Lidov Mechanism in Hydrodynamical Disks. II. Effects of binary and disk parameters

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.

preprint2015arXiv

The Kozai-Lidov Mechanism in Hydrodynamical Disks. III. Effects of Disk Mass and Self-Gravity

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.

preprint2015arXiv

The Solar System as an Exoplanetary System

With the availability of considerably more data, we revisit the question of how special our Solar System is, compared to observed exoplanetary systems. To this goal, we employ a mathematical transformation that allows for a meaningful, statistical comparison. We find that the masses and densities of the giant planets in our Solar System are very typical, as is the age of the Solar System. While the orbital location of Jupiter is somewhat of an outlier, this is most likely due to strong selection effects towards short-period planets. The eccentricities of the planets in our Solar System are relatively small compared to those in observed exosolar systems, but still consistent with the expectations for an 8-planet system (and could, in addition, reflect a selection bias towards high-eccentricity planets). The two characteristics of the Solar System that we find to be most special are the lack of super-Earths with orbital periods of days to months and the general lack of planets inside of the orbital radius of Mercury. Overall, we conclude that in terms of its broad characteristics our Solar System is not expected to be extremely rare, allowing for a level of optimism in the search for extrasolar life.

preprint2014arXiv

Giant Outbursts in Be/X-ray Binaries

Be/X-ray binary systems exhibit both periodic (Type I) X-ray outbursts and giant (Type II) outbursts, whose origin has remained elusive. We suggest that Type II X-ray outbursts occur when a highly misaligned decretion disk around the Be star becomes eccentric, allowing the compact object companion to capture a large amount of material at periastron. Using 3D smoothed particle hydrodynamics simulations we model the long term evolution of a representative Be/X-ray binary system. We find that periodic (Type I) X-ray outbursts occur when the neutron star is close to periastron for all disk inclinations. Type II outbursts occur for large misalignment angles and are associated with eccentricity growth that occurs on a timescale of about 10 orbital periods. Mass capture from the eccentric decretion disk results in an accretion disk around the neutron star whose estimated viscous time is long enough to explain the extended duration of Type II outbursts. Previous studies suggested that the outbursts are caused by a warped disk but our results suggest that this is not sufficient, the disk must be both highly misaligned and eccentric to initiate a Type II accretion event.

preprint2014arXiv

On the evolution of the CO snow line in protoplanetary disks

CO is thought to be a vital building block for prebiotic molecules that are necessary for life. Thus, understanding where CO existed in a solid phase within the solar nebula is important for understanding the origin of life. We model the evolution of the CO snow line in a protoplanetary disk. We find that the current observed location of the CO snow line in our solar system, and in the solar system analogue TW Hydra, cannot be explained by a fully turbulent disk model. With time-dependent disk models we find that the inclusion of a dead zone (a region of low turbulence) can resolve this problem. Furthermore, we obtain a fully analytic solution for the CO snow line radius for late disk evolutionary times. This will be useful for future observational attempts to characterize the demographics and predict the composition and habitability of exoplanets.

preprint2014arXiv

The Kozai-Lidov Mechanism in Hydrodynamical Disks

We use three dimensional hydrodynamical simulations to show that a highly misaligned accretion disk around one component of a binary system can exhibit global Kozai-Lidov cycles, where the inclination and eccentricity of the disk are interchanged periodically. This has important implications for accreting systems on all scales, for example, the formation of planets and satellites in circumstellar and circumplanetary disks, outbursts in X-ray binary systems and accretion on to supermassive black holes.

preprint2014arXiv

Tidal Torques on Misaligned Disks in Binary Systems

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.

preprint2013arXiv

Formation of circumbinary planets in a dead zone

Circumbinary planets have been observed at orbital radii where binary perturbations may have significant effects on the gas disk structure, on planetesimal velocity dispersion, and on the coupling between turbulence and planetesimals. Here, we note that the impact of all of these effects on planet formation is qualitatively altered if the circumbinary disk structure is layered, with a non-turbulent midplane layer (dead zone) and strongly turbulent surface layers. For close binaries, we find that the dead zone typically extends from a radius close to the inner disk edge up to a radius of around 10-20 au from the centre of mass of the binary. The peak in the surface density occurs within the dead zone, far from the inner disk edge, close to the snow line, and may act as a trap for aerodynamically coupled solids. We suggest that circumbinary planet formation may be easier near this preferential location than for disks around single stars. However, dead zones around wide binaries are less likely and hence planet formation may be more difficult there.

preprint2013arXiv

On the evolution of the snow line in protoplanetary discs II: Analytic approximations

We examine the evolution of the snow line in a protoplanetary disc that contains a dead zone (a region of zero or low turbulence). The snow line is within a self-gravitating part of the dead zone, and we obtain a fully analytic solution for its radius. Our formula could prove useful for future observational attempts to characterise the demographics of planets outside the snow line. External sources such as comic rays or X-rays from the central star can ionise the disc surface layers and allow the magneto-rotational instability to drive turbulence there. We show that provided that the surface density in this layer is less than about 50 g/cm^2, the dead zone solution exists, after an initial outbursting phase, until the disc is dispersed by photoevaporation. We demonstrate that the snow line radius is significantly larger than that predicted by a fully turbulent disc model, and that in our own solar system it remains outside of the orbital radius of the Earth. Thus, the inclusion of a dead zone into a protoplanetary disc model explains how our Earth formed with very little water.

preprint2013arXiv

Propagation of the Gravo-Magneto Disc Instability

Discs that contain dead zones are subject to the Gravo-Magneto (GM) instability that arises when the turbulence shifts from gravitational to magnetic. We have previously described this instability through a local analysis at some radius in the disc in terms of a limit cycle. A disc may be locally unstable over a radial interval. In this paper, we consider how the local instability model can describe global disc outbursts. The outburst is triggered near the middle of the range of locally unstable radii. The sudden increase in turbulence within high surface density material causes a snow plough of density that propagates both inwards and outwards. All radii inside of the trigger radius become unstable, as well as locally unstable radii outside the trigger radius. In addition, a locally stable region outside of the trigger radius may also become unstable as the gravitational instability is enhanced by the snow plough. For the circumstellar disc model we consider, we find that a quarter of the disc mass is accreted on to the central object during the outburst. The radius out to which the disc is globally unstable is twice that for which it is locally unstable.

preprint2013arXiv

The gravo-magneto disc instability with a viscous dead zone

We consider the evolution of accretion discs that contain some turbulence within a disc dead zone, a region about the disc midplane of a disc that is not sufficiently ionised for the magneto-rotational instability (MRI) to drive turbulence. In particular, we determine whether additional sources of turbulence within a dead zone are capable of suppressing gravo-magneto (GM) disc outbursts that arise from a rapid transition from gravitationally produced to MRI produced turbulence. With viscous $α$ disc models we consider two mechanisms that may drive turbulence within the dead zone. First, we examine a constant $α$ parameter within the dead zone. This may be applicable to hydrodynamical instability, such as baroclinic instability, where the turbulence level is independent of the MRI active surface layer properties. In this case, we find that the disc will not become stable to the outbursts unless the dead zone turbulent viscosity is comparable to that in the MRI active surface layers. Under such conditions, the accretion rate through the dead zone must be larger than that through the MRI active layers. In a second model, we assume that the accretion flow though the dead zone is a constant fraction (less than unity) of that through the active layers. This may be applicable to turbulence driven by hydrodynamic waves that are excited by the MRI active surface layers. We find that the instability is hardly affected by the viscous dead zone. In both cases however, we find that the triggering of the MRI during the outburst may be due to the heating from the viscosity in the dead zone, rather than self-gravity. Thus, neither mechanism for generating turbulence within the dead zone can significantly stabilise a disc or the resulting outburst behaviour.

preprint2013arXiv

Two timescale dispersal of magnetized protoplanetary disks

Protoplanetary disks are likely to be threaded by a weak net flux of vertical magnetic field that is a remnant of the much larger fluxes present in molecular cloud cores. If this flux is approximately conserved its dynamical importance will increase as mass is accreted, initially by stimulating magnetorotational disk turbulence and subsequently by enabling wind angular momentum loss. We use fits to numerical simulations of ambipolar dominated disk turbulence to construct simplified one dimensional evolution models for weakly magnetized protoplanetary disks. We show that the late onset of significant angular momentum loss in a wind can give rise to "two timescale" disk evolution in which a long phase of viscous evolution precedes rapid dispersal as the wind becomes dominant. The wide dispersion in disk lifetimes could therefore be due to varying initial levels of net flux. Magnetohydrodynamic (MHD) wind triggered dispersal differs from photoevaporative dispersal in predicting mass loss from small (less that 1 AU) scales, where thermal winds are suppressed. Our specific models are based on a limited set of simulations that remain uncertain, but qualitatively similar evolution appears likely if mass is lost from disks more quickly than flux, and if MHD winds become important as the plasma beta decreases.

preprint2012arXiv

Dead Zones around Young Stellar Objects: FU Orionis Outbursts and Transition Discs

We perform global time-dependent simulations of an accretion disc around a young stellar object with a dead zone (a region where the magneto-rotational instability cannot drive turbulence because the material is not sufficiently ionised). For infall accretion rates on to the disc of around 10^-7 Msun/yr, dead zones occur if the critical magnetic Reynolds number is larger than about 10^4. We model the collapse of a molecular gas cloud. At early times when the infall accretion rate is high, the disc is thermally ionised and fully turbulent. However, as the infall accretion rate drops, a dead zone may form if the critical magnetic Reynolds number is sufficiently large, otherwise the disc remains fully turbulent. With a dead zone the disc can become unstable to the gravo-magneto instability. The mass of the star grows in large accretion outbursts that may explain FU Orionis events. At late times there is not sufficient mass in the disc for outbursts to occur but the dead zone becomes even more prominent as the disc cools. Large inner dead zones in the later stages of disc evolution may help to explain observations of transition discs with an inner hole.

preprint2012arXiv

Dead Zones In Circumplanetary Discs as Formation Sites For Regular Satellites

Regular satellites in the solar system are thought to form within circumplanetary discs. We consider a model of a layered circumplanetary disc that consists of a nonturbulent midplane layer and and strongly turbulent disc surface layers. The dead zone provides a favorable site for satellite formation. It is a quiescent environment that permits the growth of solid bodies. Viscous torques within the disc cause it to expand to a substantial fraction of its Hill radius (~0.4 R_H) where tidal torques from the central star remove its angular momentum. For certain parameters, the dead zone develops into a high density substructure well inside the Hill sphere. The radial extent of the dead zone may explain the compactness of the regular satellites orbits for Jupiter and Saturn. The disc temperatures can be low enough to be consistent with the high ice fractions of Ganymede and Callisto.

preprint2012arXiv

On the Evolution of the Snow Line in Protoplanetary Discs

We model the evolution of the snow line in a protoplanetary disc. If the magneto-rotational instability (MRI) drives turbulence throughout the disc, there is a unique snow line outside of which the disc is icy. The snow line moves closer to the star as the infall accretion rate drops. Because the snow line moves inside the radius of the Earth's orbit, the formation of our water-devoid planet is difficult with this model. However, protoplanetary discs are not likely to be sufficiently ionised to be fully turbulent. A dead zone at the mid-plane slows the flow of material through the disc and a steady state cannot be achieved. We therefore model the evolution of the snow line also in a time-dependent disc with a dead zone. As the mass is accumulating, the outer parts of the dead zone become self gravitating, heat the massive disc and thus the outer snow line does not come inside the radius of the Earth's orbit, contrary to the fully turbulent disc model. There is a second, inner icy region, within the dead zone, that moves inwards of the Earth's orbit after a time of about 10^6 yr. With this model there is sufficient time and mass in the disc for the Earth to form from water-devoid planetesimals at a radius of 1 AU. Furthermore, the additional inner icy region predicted by this model may allow for the formation of giant planets close to their host star without the need for much migration.

preprint2012arXiv

On the formation and evolution of asteroid belts and their potential significance for life

Suggestions have been made that asteroid belts may be important both for the existence of life and perhaps even for the evolution of complex life on a planet. Using numerical models for protoplanetary discs we calculate the location of the snow line, and we propose that asteroid belts are most likely to form in its vicinity. We then show that observations of warm dust in exo-solar systems, thought to be produced by collisions between asteroids in a belt, indicate that asteroid belts (when they exist), indeed coincide with the radial location and the temperature of the snow line. Giant planets form outside the snow line and prevent planet formation just inside of their orbit creating an asteroid belt there. However, the migration of giant planets through the asteroid belt likely disperses the compact formation. We examine existing observations of giant exo-planets and find that less than 4% are at radial locations outside of the snow line. This definitely may be the consequence of observational selection effects. However, with this caveat in mind, we point out that the dearth of giant planets outside the snow line may also suggest that compact asteroid belts are not common, and more speculatively that complex life may not expected in most of the currently observed systems.

preprint2011arXiv

Dead Zones around Young Stellar Objects: Dependence on Physical Parameters

Angular momentum is transported outwards through an accretion disc by magnetohydrodynamical (MHD) turbulence thus allowing material to accrete on to the central object. The magneto-rotational instability (MRI) requires a minimum ionisation fraction to drive turbulence in a disc. The inner parts of the disc around a young stellar object are sufficiently hot to be thermally ionised. Further out, cosmic rays ionise the surface layers and a dead zone forms at the mid-plane where the disc is too cool for the MRI to operate. The surface density in the turbulent active layer is often assumed to be constant with radius because the cosmic rays penetrate a constant layer. However, if a critical magnetic Reynolds number, Re_{M,crit}, is used to determine the extent of the dead zone, the surface density in the layer generally increases with radius. For small critical magnetic Reynolds number of order 1, the constant layer approximation may be a reasonable fit. However, MHD simulations suggest the critical magnetic Reynolds number may be much larger, of order 10^4. Analytical fits for the surface density in the magnetic active layer show that Σ_m \propto Re_{M,crit}^{-2} R^{9/2} T^{\,2}$, at temperature T and radius R, are a good fit for higher critical magnetic Reynolds number. For the metallicity variation between our galaxy, the LMC and the SMC, there should be no significant difference in the extent of the dead zone. Observations suggest an increase in the lifetime of the disc with decreasing metallicity that cannot be explained by the dead zone structure (ignoring possible differences in dust abundances).

preprint2011arXiv

On Orbital Period Changes in Nova Outbursts

We propose a new mechanism that produces an orbital period change during a nova outburst. When the ejected material carries away the specific angular momentum of the white dwarf, the orbital period increases. A magnetic field on the surface of the secondary star forces a fraction of the ejected material to corotate with the star, and hence the binary system. The ejected material thus takes angular momentum from the binary orbit and the orbital period decreases. We show that for sufficiently strong magnetic fields on the surface of the secondary star, the total change to the orbital period could even be negative during a nova outburst, contrary to previous expectations. Accurate determinations of pre- and post-outburst orbital periods of recurrent nova systems could test the new mechanism, in addition to providing meaningful constraints on otherwise difficult to measure physical quantities. We apply our mechanism to outbursts of the recurrent nova U Sco.

preprint2011arXiv

The Gravo-Magneto Limit Cycle in Accretion Disks

Previous theoretical studies have found that repeating outbursts can occur in certain regions of an accretion disk, due to sudden transitions in time from gravitationally produced turbulence to magnetically produced turbulence. We analyze the disk evolution in a state diagram that plots the mass accretion rate versus disk surface density. We determine steady state accretion branches that involve gravitational and magnetic sources of turbulence. Using time-dependent numerical disk simulations, we show that cases having outbursts track along a nonsteady 'dead zone' branch and some steady state accretion branches. The outburst is the result of a rapid inter-branch transition. The gravo-magneto outbursts are then explained on this diagram as a limit cycle that is analogous to the well-known S-curve that has been applied to dwarf nova outbursts. The diagram and limit cycle provide a conceptual framework for understanding the nature of the outbursts that may occur in accretion disks of all scales, from circumplanetary to protoplanetary to AGN accretion disks.

preprint2011arXiv

Tidal Warping of Be Star Decretion Discs

Rapidly rotating Be stars are observed as shell stars when the decretion disc is viewed edge on. Transitions between the two implies that the discs may be warped and precessing. Type II X-ray outbursts are thought to occur when the warped disc interacts with the fast stellar wind. We suggest that tides from a misaligned companion neutron star can cause the observed effects. We make numerical models of a Be star decretion disc in which the spin of the Be star is misaligned with the orbital axis of a neutron star companion. Tidal torques from the neutron star truncate the disc at a radius small enough that the neutron star orbit does not intersect the disc unless the eccentricity or misalignment is very large. A magnetic torque from the Be star that is largest at the equator, where the rotation is fastest, is approximated by an inner boundary condition. There are large oscillations in the mass and inclination of the disc as it moves towards a steady state. These large variations may explain the observed changes from Be star to Be shell star and vice-versa and also the Type II X-ray outbursts. We find the tidal timescale on which the disc warps, precesses and reaches a steady state to be around a year up to a few hundred years. If present, the oscillations in mass and disc inclination occur on a fraction of this timescale depending on the orbital parameters of the binary. The timescales associated with the tidal torque for observed Be star binaries suggest that these effects are important in all but the longest period binaries.

preprint2010arXiv

On the Orbital Period Change in the Recurrent Nova U Scorpii

The orbital period of the recurrent nova U Sco has been observed to decrease during the 1999 outburst. In an outburst mass is ejected from the surface of the white dwarf. The separation of the binary system widens and the orbital period increases. We find that magnetic braking between outbursts, mass transfer to the companion, and frictional angular momentum losses during outbursts are all too small to account for this unexpected change. We find, however, that if the secondary has a sufficiently strong magnetic field, B=8x10^3 G, then the ejected material can couple to it and corrotate with the system. The ejected material gains angular momentum while the binary system loses it and the period decreases. If such a strong magnetic field is indeed present, then we predict that a period decrease should be observed also during the current 2010 outburst. If, however, the presence of such a field can be ruled out observationally, then the cause for the period decrease (if confirmed) remains unknown.

preprint2010arXiv

Tidal Truncation of Circumplanetary Discs

We analyse some properties of circumplanetary discs. Flow through such discs may provide most of the mass to gas giant planets, and such discs are likely sites for the formation of regular satellites. We model these discs as accretion discs subject to the tidal forces of the central star. The tidal torques from the star remove the disc angular momentum near the disc outer edge and permit the accreting disc gas to lose angular momentum at the rate appropriate for steady accretion. Circumplanetary discs are truncated near the radius where periodic ballistic orbits cross, where tidal forces on the disc are strong. This radius occurs at approximately 0.4 r_H for the planet Hill radius r_H. During the T Tauri stage of disc accretion, the disc is fairly thick with aspect ratio H/r > 0.2 and the disc edge tapering occurs over a radial scale ~ H ~ 0.1 r_H. For a circular or slightly eccentric orbit planet, no significant resonances lie within the main body of the disc. Tidally driven waves involving resonances nonetheless play an important role in truncating the disc, especially when it is fairly thick. We model the disc structure using one dimensional time-dependent and steady-state models and also two dimensional SPH simulations. The circumplanetary disc structure depends on the variation of the disc turbulent viscosity with radius and is insensitive to the angular distribution of the accreting gas. Dead zones may occur within the circumplanetary disc and result in density structures. If the disc is turbulent throughout, the predicted disc structure near the location of the regular Jovian and Saturnian satellites is smooth with no obvious feature that would favor formation at their current locations.

preprint2009arXiv

Supernova Kicks and Misaligned Be Star Binaries

Be stars are rapidly spinning B stars surrounded by an outflowing disc of gas in Keplerian rotation. Be star/X-ray binary systems contain a Be star and a neutron star. They are found to have non-zero eccentricities and there is evidence that some systems have a misalignment between the spin axis of the star and the spin axis of the binary orbit. The eccentricities in these systems are thought to be caused by a kick to the neutron star during the supernova that formed it. Such kicks would also give rise to misalignments. In this paper we investigate the extent to which the same kick distribution can give rise to both the observed eccentricity distribution and the observed misalignments. We find that a Maxwellian distribution of velocity kicks with a low velocity dispersion, $σ_k \approx 15\rm km s^{-1}$, is consistent with the observed eccentricity distribution but is hard to reconcile with the observed misalignments, typically $i \ge 25^\circ$. Alternatively a higher velocity kick distribution, $σ_k = 265 \rm km s^{-1}$, is consistent with the observed misalignments but not with the observed eccentricities, unless post-supernova circularisation of the binary orbits has taken place. We discuss briefly how this might be achieved.

preprint2009arXiv

Supernova Kicks and Misaligned Microquasars

The low-mass X-ray binary microquasar GRO J1655-40 is observed to have a misalignment between the jets and the binary orbital plane. Since the current black hole spin axis is likely to be parallel to the jets, this implies a misalignment between the spin axis of the black hole and the binary orbital plane. It is likely the black holes formed with an asymmetric supernova which caused the orbital axis to misalign with the spin of the stars. We ask whether the null hypothesis that the supernova explosion did not affect the spin axis of the black hole can be ruled out by what can be deduced about the properties of the explosion from the known system parameters. We find that this null hypothesis cannot be disproved but we find that the most likely requirements to form the system include a small natal black hole kick (of a few tens of km/s) and a relatively wide pre-supernova binary. In such cases the observed close binary system could have formed by tidal circularisation without a common envelope phase.

preprint2009arXiv

The Shape of an Accretion Disc in a Misaligned Black Hole Binary

We model the overall shape of an accretion disc in a semi-detached binary system in which mass is transfered on to a spinning black hole the spin axis of which is misaligned with the orbital rotation axis. We assume the disc is in a steady state. Its outer regions are subject to differential precession caused by tidal torques of the companion star. These tend to align the outer parts of the disc with the orbital plane. Its inner regions are subject to differential precession caused by the Lense-Thirring effect. These tend to align the inner parts of the disc with the spin of the black hole. We give full numerical solutions for the shape of the disc for some particular disc parameters. We then show how an analytic approximation to these solutions can be obtained for the case when the disc surface density varies as a power law with radius. These analytic solutions for the shape of the disc are reasonably accurate even for large misalignments and can be simply applied for general disc parameters. They are particularly useful when the numerical solutions would be slow.

preprint2006arXiv

Critical Mass Transfer in Double-Degenerate Type Ia Supernovae

Doubly-degenerate binary systems consisting of two white dwarfs both composed of carbon and oxygen and close enough that mass is transferred from the less massive to the more massive are possible progenitors of type Ia supernovae. If the mass transfer rate is slow enough that the accreting white dwarf can reach 1.38 solar masses then it can ignite carbon degenerately at its centre. This can lead to a thermonuclear runaway and thence a supernova explosion. However if the accretion rate is too high the outer layers of the white dwarf heat up too much and carbon ignites there non-degenerately. A series of mild carbon flashes can then propagate inwards and convert the carbon to neon relatively gently. There is no thermonuclear runaway and no supernova. We have determined the critical rate for accretion and find it to be generally about two-fifths of the Eddington rate. In the case of actual bianry stars the mass transfer rate falls off as mass is transferred so we find that this rate can be very much exceeded at the onset of mass transfer. We determine combinations of white dwarf masses that can lead to central degenerate carbon ignition.