Accretion disc dynamics in extreme mass ratio compact binaries

An analysis is presented of a numerical investigation of the dynamics and geometry of accretion discs in binary systems with mass ratios q < 0.1, applicable to ultra-compact X-ray binaries, AM CVn stars and very short period cataclysmic variables. The steady-state geometry of the disc in the binary reference frame is found to be quite different from that expected at higher mass ratios. For q ~ 0.1, the disc takes on the usual elliptical shape, with the major axis aligned perpendicular to the line of centres of the two stars. However, at smaller mass ratios the elliptical gaseous orbits in the outer regions of the disc are rotated in the binary plane. The angle of rotation increases with gas temperature, but is found to vary inversely with q. At q = 0.01, the major axis of these orbits is aligned almost parallel to the line of centres of the two stars. These effects may be responsible for the similar disc structure inferred from Doppler tomography of the AM CVn star GP Com (Morales-Rueda et al. 2003), which has q = 0.02. The steady-state geometry at low mass ratios is not predicted by an inviscid, restricted three-body model of gaseous orbits; it is related to the effects of tidal-viscous truncation of the disc near the Roche lobe boundary. Since the disc geometry can be inferred observationally for some systems, it is proposed that this may offer a useful diagnostic for the determination of mass ratios in ultra-compact binaries.Comment: 17 pages, 9 figures, 7 in colour. Accepted for publication in MNRAS. Plain article formatting to get round arXiv problems with mn2e.st

Secular interactions between inclined planets and a gaseous disk

In a planetary system, a secular particle resonance occurs at a location where the precession rate of a test particle (e.g. an asteroid) matches the frequency of one of the precessional modes of the planetary system. We investigate the secular interactions of a system of mutually inclined planets with a gaseous protostellar disk that may contain a secular nodal particle resonance. We determine the normal modes of some mutually inclined planet-disk systems. The planets and disk interact gravitationally, and the disk is internally subject to the effects of gas pressure, self-gravity, and turbulent viscosity. The behavior of the disk at a secular resonance is radically different from that of a particle, owing mainly to the effects of gas pressure. The resonance is typically broadened by gas pressure to the extent that global effects, including large-scale warps, dominate. The standard resonant torque formula is invalid in this regime. Secular interactions cause a decay of the inclination at a rate that depends on the disk properties, including its mass, turbulent viscosity, and sound speed. For a Jupiter-mass planet embedded within a minimum-mass solar nebula having typical parameters, dissipation within the disk is sufficient to stabilize the system against tilt growth caused by mean-motion resonances.Comment: 30 pages, 6 figures, to be published in The Astrophysical Journa

Source to Accretion Disk Tilt

Many different system types retrogradely precess, and retrograde precession could be from a tidal torque by the secondary on a misaligned accretion disk. However, a source to cause and maintain disk tilt is unknown. In this work, we show that accretion disks can tilt due to a force called lift. Lift results from differing gas stream supersonic speeds over and under an accretion disk. Because lift acts at the disk's center of pressure, a torque is applied around a rotation axis passing through the disk's center of mass. The disk responds to lift by pitching around the disk's line of nodes. If the gas stream flow ebbs, then lift also ebbs and the disk attempts to return to its original orientation. To first approximation, lift does not depend on magnetic fields or radiation sources but does depend on mass and the surface area of the disk. Also, for disk tilt to be initiated, a minimum mass transfer rate must be exceeded. For example, a $10^{-11}M_{\odot}$ disk around a 0.8$M_{\odot}$ compact central object requires a mass transfer rate greater than $\sim10^{-13}$M$_{\odot}$yr$^{-1}$, a value well below known mass transfer rates in Cataclysmic Variable Dwarf Novae systems that retrogradely precess and that exhibit negative superhumps in their light curves and a value well below mass transfer rates in protostellar forming systems

Forced oscillations in a hydrodynamical accretion disk and QPOs

This is the second of a series of papers aimed to look for an explanation on the generation of high frequency quasi-periodic oscillations (QPOs) in accretion disks around neutron star, black hole, and white dwarf binaries. The model is inspired by the general idea of a resonance mechanism in the accretion disk oscillations as was already pointed out by Abramowicz & Klu{\'z}niak (\cite{Abramowicz2001}). In a first paper (P\'etri \cite{Petri2005a}, paper I), we showed that a rotating misaligned magnetic field of a neutron star gives rise to some resonances close to the inner edge of the accretion disk. In this second paper, we suggest that this process does also exist for an asymmetry in the gravitational potential of the compact object. We prove that the same physics applies, at least in the linear stage of the response to the disturbance in the system. This kind of asymmetry is well suited for neutron stars or white dwarfs possessing an inhomogeneous interior allowing for a deviation from a perfectly spherically symmetric gravitational field. We show by a linear analysis that the disk initially in a cylindrically symmetric stationary state is subject to three kinds of resonances: a corotation resonance, a Lindblad resonance due to a driven force and a parametric sonance. The highest kHz QPOs are then interpreted as the orbital frequency of the disk at locations where the response to the resonances are maximal. It is also found that strong gravity is not required to excite the resonances.Comment: Accepte

Hydrodynamic simulations of irradiated secondaries in dwarf novae

We investigate numerically the surface flow on the secondary star during outbursts. We use a simple model for the irradiation and the geometry of the secondary star: the irradiation temperature is treated as a free parameter and the secondary is replaced by a spherical star with a space-dependent Coriolis force that mimics the effect of the Roche geometry. The Euler equations are solved in spherical coordinates with the TVD-MacCormack scheme. We show that the Coriolis force leads to the formation of a circulation flow from high latitude region to the close vicinity of the $L_1$ point. However no heat can be efficiently transported to the $L_1$ region due to the rapid radiative cooling of the hot material as it enters the equatorial belt shadowed from irradiation. Under the assumption of hydrostatic equilibrium, the Coriolis force could lead to a moderate increase of the mass transfer rate by pushing the gas in the vertical direction at the $L_1$ point, but only during the initial phases of the outburst (about 15 -- 20 orbital periods). We conclude that the Coriolis force does not prevent a flow from the heated regions of the secondary towards the $L_1$ region, at least during the initial phase of an outburst, but the resulting increase of the mass transfer rate is moderate, and it is unlikely to be able to account for the duration of long outbursts.Comment: 11 pages, 11 figures, accepted for publication in A&

A warped m=2 water maser disc in V778 Cyg?

The silicate carbon star V778 Cyg is a source of 22 GHz water maser emission which was recently resolved by MERLIN. Observations revealed an elongated S-like structure along which the velocities of the maser features show a linear dependence on the impact parameter. This is consistent with a doubly-warped m=2 disc observed edge-on. Water masers and silicate dust emission (detected by IRAS and ISO) have a common origin in O-rich material and are likely to be co-located in the disc. We propose a detailed self-consistent model of a masing gas-dust disc around a companion to the carbon star in a binary system, which allows us to estimate the companion mass of 1.7 +- 0.1 M_sun, the disc radius of 40 +-3 AU and the distance between companions of about 80 AU. Using a dust-gas coupling model for water masing, we calculate the maser power self-consistently, accounting for both the gas and the dust energy balances. Comparing the simulation results with the observational data, we deduce the main physical parameters of the masing disc, such as the gas and dust temperatures and their densities. We also present an analysis of the stability of the disc.Comment: 7 pages, 5 figures. This paper is accepted for publication in MNRA

Accretion Discs with an Inner Spiral Density Wave

In Montgomery (2009a), we show that accretion discs in binary systems could retrogradely precess by tidal torques like the Moon and the Sun on a tilted, spinning, non-spherical Earth. In addition, we show that the state of matter and the geometrical shape of the celestial object could significantly affect the precessional value. For example, a Cataclysmic Variable (CV) Dwarf Novae (DN) non-magnetic system that shows negative superhumps in its light curve can be described by a retrogradely precessing, differentially rotating, tilted disc. Because the disc is a fluid and because the gas stream overflows the tilted disc and particles can migrate into inner disc annuli, coupled to the disc could be a retrogradely precessing inner ring that is located near the innermost annuli of the disc. However, numerical simulations by Bisikalo et al. (2003, 2004) and this work show that an inner spiral density wave can be generated instead of an inner ring. Therefore, we show that retrograde precession in non-magnetic, spinning, tilted CV DN systems can equally be described by a retrogradely precessing and differentially rotating disc with an attached retrogradely precessing inner spiral density wave so long as the wave appears at the same radius as the ring and within the plane of the tilted disc. We find that the theoretical results generated in this work agree well with the theoretical results presented in Montgomery (2009a) and thus with the numerical simulations and select CV DN systems in Montgomery (2009b) that may have a main sequence secondary. Therefore, pressure effects do need to be considered in CV DN systems that exhibit negative superhumps if the accretion discs are tilted and have an inner spiral density wave that is in the plane of the disc

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.Comment: 20 pages 16 figures, 3 tables. To appear in The Astrophysical Journal vol.652 (December 1, 2006 issue

Three-dimensional Calculations of High and Low-mass Planets Embedded in Protoplanetary Discs

We analyse the non-linear, three-dimensional response of a gaseous, viscous protoplanetary disc to the presence of a planet of mass ranging from one Earth mass (1 M$_e$) to one Jupiter mass (1 M$_J$) by using the ZEUS hydrodynamics code. We determine the gas flow pattern, and the accretion and migration rates of the planet. The planet is assumed to be in a fixed circular orbit about the central star. It is also assumed to be able to accrete gas without expansion on the scale of its Roche radius. Only planets with masses M \gsim 0.1 M$_J$ produce significant perturbations in the disc's surface density. The flow within the Roche lobe of the planet is fully three-dimensional. Gas streams generally enter the Roche lobe close to the disc midplane, but produce much weaker shocks than the streams in two-dimensional models. The streams supply material to a circumplanetary disc that rotates in the same sense as the planet's orbit. Much of the mass supply to the circumplanetary disc comes from non-coplanar flow. The accretion rate peaks with a planet mass of approximately 0.1 M$_J$ and is highly efficient, occurring at the local viscous rate. The migration timescales for planets of mass less than 0.1 M$_J$, based on torques from disc material outside the planets' Roche lobes, are in excellent agreement with the linear theory of Type I (non-gap) migration for three-dimensional discs. The transition from Type I to Type II (gap) migration is smooth, with changes in migration times of about a factor of 2. Starting with a core which can undergo runaway growth, a planet can gain up to a few M$_J$ with little migration. Planets with final masses of order 10 M$_J$ would undergo large migration, which makes formation and survival difficult.Comment: Accepted by MNRAS, 18 pages, 13 figures (6 degraded resolution). Paper with high-resolution figures available at http://www.astro.ex.ac.uk/people/mbate

Tidal Decay of Close Planetary Orbits

The 4.2-day orbit of the newly discovered planet around 51~Pegasi is formally unstable to tidal dissipation. However, the orbital decay time in this system is longer than the main-sequence lifetime of the central star. Given our best current understanding of tidal interactions, a planet of Jupiter's mass around a solar-like star could have dynamically survived in an orbit with a period as short as $\sim10\,$hr. Since radial velocities increase with decreasing period, we would expect to find those planets close to the tidal limit first and, unless this is a very unusual system, we would expect to find many more. We also consider the tidal stability of planets around more evolved stars and we re-examine in particular the question of whether the Earth can dynamically survive the red-giant phase in the evolution of the Sun.Comment: AAS LaTeX macros v.4, 14 pages, 2 postscript figures, also available from http://ensor.mit.edu/~rasio/, to appear in Ap