Source author record

Stephen H. Lubow

Stephen H. Lubow 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

34works
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

34 published item(s)

preprint2022arXiv

Accretion onto a binary from a polar circumbinary disc

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

preprint2022arXiv

Circumbinary disk evolution in the presence of an outer companion star

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

preprint2022arXiv

Eccentric dust ring formation in Kozai-Lidov gas disks

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

preprint2022arXiv

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

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

preprint2020arXiv

Formation of the polar debris disc around 99 Herculis

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

preprint2020arXiv

Polar planets around highly eccentric binaries are the most stable

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

preprint2020arXiv

Sustained Kozai-Lidov oscillations in misaligned circumstellar gas disks

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

preprint2019arXiv

Multi-planet disc interactions in binary systems

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

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

preprint2016arXiv

Turbulence, Transport and Waves in Ohmic Dead Zones

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.

preprint2016arXiv

Version 1 of the Hubble Source Catalog

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.

preprint2015arXiv

Resonances in retrograde circumbinary discs

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.

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 Inner Disk Structure, Disk-Planet Interactions, and Temporal Evolution in the Beta Pictoris System: A Two-Epoch HST/STIS Coronagraphic Study

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.

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.

preprint2014arXiv

An Analytic Model for Buoyancy Resonances in Protoplanetary Disks

Zhu, Stone, and Rafikov (2012) found in 3D shearing box simulations a new form of planet-disk interaction that they attributed to a vertical buoyancy resonance in the disk. We describe an analytic linear model for this interaction. We adopt a simplified model involving azimuthal forcing that produces the resonance and permits an analytic description of its structure. We derive an analytic expression for the buoyancy torque and show that the vertical torque distribution agrees well with results of Athena simulations and a Fourier method for linear numerical calculations carried out with the same forcing. The buoyancy resonance differs from the classic Lindblad and corotation resonances in that the resonance lies along tilted planes. Its width depends on damping effects and is independent of the gas sound speed. The resonance does not excite propagating waves. At a given large azimuthal wavenumber k_y > 1/h (for disk thickness h), the buoyancy resonance exerts a torque over a region that lies radially closer to the corotation radius than the Lindblad resonance. Because the torque is localized to the region of excitation, it is potentially subject to the effects of nonlinear saturation. In addition, the torque can be reduced by the effects of radiative heat transfer between the resonant region and its surroundings. For each azimuthal wavenumber, the resonance establishes a large scale density wave pattern in a plane within the disk.

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

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

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.

preprint2012arXiv

Catalog Matching with Astrometric Correction and its Application to the Hubble Legacy Archive

Object cross-identification in multiple observations is often complicated by the uncertainties in their astrometric calibration. Due to the lack of standard reference objects, an image with a small field of view can have significantly larger errors in its absolute positioning than the relative precision of the detected sources within. We present a new general solution for the relative astrometry that quickly refines the World Coordinate System of overlapping fields. The efficiency is obtained through the use of infinitesimal 3-D rotations on the celestial sphere, which do not involve trigonometric functions. They also enable an analytic solution to an important step in making the astrometric corrections. In cases with many overlapping images, the correct identification of detections that match together across different images is difficult to determine. We describe a new greedy Bayesian approach for selecting the best object matches across a large number of overlapping images. The methods are developed and demonstrated on the Hubble Legacy Archive, one of the most challenging data sets today. We describe a novel catalog compiled from many Hubble Space Telescope observations, where the detections are combined into a searchable collection of matches that link the individual detections. The matches provide descriptions of astronomical objects involving multiple wavelengths and epochs. High relative positional accuracy of objects is achieved across the Hubble images, often sub-pixel precision in the order of just a few milli-arcseconds. The result is a reliable set of high-quality associations that are publicly available online.

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

Three Dimensional MHD Simulation of Circumbinary Accretion Disks: Disk Structures and Angular Momentum Transport

We present the first three-dimensional magnetohydrodynamic (MHD) simulations of a circumbinary disk surrounding an equal mass binary. The binary maintains a fixed circular orbit of separation $a$. As in previous hydrodynamical simulations, strong torques by the binary can maintain a gap of radius $\simeq 2a$. Streams curve inward from $r \simeq 2a$ toward the binary; some of their mass passes through the inner boundary, while the remainder swings back out to the disk. However, we also find that near its inner edge the disk develops both a strong $m=1$ asymmetry and growing orbital eccentricity. Because the MHD stresses introduce more matter into the gap, the total torque per unit disk mass is $\simeq 14$ times larger than found previously. The inner boundary accretion rate per unit disk mass is $\simeq 40$ times greater than found from previous hydrodynamical calculations. The implied binary shrinkage rate is determined by a balance between the rate at which the binary gains angular momentum by accretion and loses it by gravitational torque. The large accretion rate brings these two rates nearly into balance, but in net, we find that $\dot a/a < 0$, and its magnitude is about 2.7 times larger than predicted by the earlier hydrodynamic simulations. If the binary comprises two massive black holes, the accretion rate may be great enough for one or both to be AGN, with consequences for the physical state of the gas both in the disk body and in its inner gap.

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

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

Eccentricity Growth Rates of Tidally Distorted Discs

We consider discs that orbit a central object and are tidally perturbed by a circular orbit companion. Such discs are sometimes subject to an eccentric instability due to the effects of certain resonances. Eccentric instabilities may be present in planetary rings perturbed by satellites, protostellar discs perturbed by planets, and discs in binary star systems. Although the basic mechanism for eccentric instability is well understood, the detailed response of a gaseous disc to such an instability is not understood. We apply a linear eccentricity evolution equation developed by Goodchild and Ogilvie. We explore how the eccentricity is distributed in such a disc and how the distribution in turn affects the instability growth rate for a range of disc properties. We identify a disc mode, termed the superhump mode, that is likely at work in the superhump binary star case. The mode results from the excitation of the fundamental free precession mode. We determine an analytic expression for the fundamental free mode precession rate that is applicable to a sufficiently cool disc. Depending on the disc sound speed and disc edge location, other eccentric modes can grow faster than the superhump mode and dominate.

preprint2010arXiv

Planet Migration

Planet migration is the process by which a planet's orbital radius changes in time. The main agent for causing gas giant planet migration is the gravitational interaction of the young planet with the gaseous disk from which it forms. We describe the migration rates resulting from these interactions based on a simple model for disk properties. These migration rates are higher than is reasonable for planet survival. We discuss some proposed models for which the migration rates are lower. There are major uncertainties in migration rates due to a lack of knowledge about the detailed physical properties of disks. We also describe some additional forms of migration.

preprint2010arXiv

Three-dimensional Disk-Planet Torques in a Locally Isothermal Disk

We determine an expression for the Type I planet migration torque involving a locally isothermal disk, with moderate turbulent viscosity (~0.0005 < alpha < ~0.05), based on three-dimensional nonlinear hydrodynamical simulations. The radial gradients (in a dimensionless logarithmic form) of density and temperature are assumed to be constant near the planet. We find that the torque is roughly equally sensitive to the surface density and temperature radial gradients. Both gradients contribute to inward migration when they are negative. Our results indicate that two-dimensional calculations with a smoothed planet potential, used to account for the effects of the third dimension, do not accurately determine the effects of density and temperature gradients on the three-dimensional torque. The results suggest that substantially slowing or stopping planet migration by means of changes in disk opacity or shadowing is difficult and appears unlikely for a disk that is locally isothermal. The scalings of the torque and torque density with planet mass and gas sound speed follow the expectations of linear theory. We also determine an improved formula for the torque density distribution that can be used in one-dimensional long-term evolution studies of planets embedded in locally isothermal disks. This formula can be also applied in the presence of mildly varying radial gradients and of planets that open gaps. We illustrate its use in the case of migrating super-Earths and determine some conditions sufficient for survival.

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.

preprint2006arXiv

Evolution of Giant Planets in Eccentric Disks

We investigate the interaction between a giant planet and a viscous circumstellar disk by means of high-resolution, two-dimensional hydrodynamical simulations. We consider planet masses that range from 1 to 3 Jupiter masses (Mjup) and initial orbital eccentricities that range from 0 to 0.4. We find that a planet can cause eccentricity growth in a disk region adjacent to the planet's orbit, even if the planet's orbit is circular. Disk-planet interactions lead to growth in a planet's orbital eccentricity. The orbital eccentricities of a 2 Mjup and a 3 Mjup planet increase from 0 to 0.11 within about 3000 orbits. Over a similar time period, the orbital eccentricity of a 1 Mjup planet grows from 0 to 0.02. For a case of a 1 Mjup planet with an initial eccentricity of 0.01, the orbital eccentricity grows to 0.09 over 4000 orbits. Radial migration is directed inwards, but slows considerably as a planet's orbit becomes eccentric. If a planet's orbital eccentricity becomes sufficiently large, e > ~0.2, migration can reverse and so be directed outwards. The accretion rate towards a planet depends on both the disk and the planet orbital eccentricity and is pulsed over the orbital period. Planet mass growth rates increase with planet orbital eccentricity. For e~0.2 the mass growth rate of a planet increases by approximately 30% above the value for e=0. For e > ~0.1, most of the accretion within the planet's Roche lobe occurs when the planet is near the apocenter. Similar accretion modulation occurs for flow at the inner disk boundary which represents accretion toward the star.