The evolution of a supermassive binary caused by an accretion disc

The interaction of a massive binary and a non-self-gravitating circumbinary accretion disc is considered. The shape of the stationary twisted disc produced by the binary is calculated. It is shown that the inner part of the disc must lie in the binary orbital plane for any value of viscosity. When the inner disc midplane is aligned with the binary orbital plane on the scales of interest and it rotates in the same sense as the binary, the modification of the disc structure and the rate of decay of the binary orbit, assumed circular, due to tidal exchange of angular momentum with the disc, are calculated. It is shown that the modified disc structure is well described by a self-similar solution of the non-linear diffusion equation governing the evolution of the disc surface density. The calculated time scale for decay of the binary orbit is always smaller than the "accretion" time $t_{acc}=m/{\dot M}$ ($m$ is the mass of the secondary component, and $\dot M$ is the disc accretion rate), and is determined by ratio of secondary mass $m$, assumed to be much smaller than the primary mass, the disc mass inside the initial binary orbit, and the form of viscosity in the disc.Comment: to be published in MNRA

Numerical simulations of the type III migration:I. Disc model and convergence tests

We investigate the fast (type III) migration regime of high-mass protoplanets orbiting in protoplanetary disks. This type of migration is dominated by corotational torques. We study the details of flow structure in the planet's vicinity, the dependence of migration rate on the adopted disc model, and the numerical convergence of models (independence of certain numerical parameters such as gravitational softening). We use two-dimensional hydrodynamical simulations with adaptive mesh refinement,based on the FLASH code with improved time-stepping scheme. We perform global disk simulations with sufficient resolution close to the planet, which is allowed to freely move throughout the grid. We employ a new type of equation of state in which the gas temperature depends on both the distance to the star and planet, and a simplified correction for self-gravity of the circumplanetary gas. We find that the migration rate in the type III migration regime depends strongly on the gas dynamics inside the Hill sphere (Roche lobe of the planet) which, in turn, is sensitive to the aspect ratio of the circumplanetary disc. Furthermore, corrections due to the gas self-gravity are necessary to reduce numerical artifacts that act against rapid planet migration. Reliable numerical studies of Type III migration thus require consideration of both the thermal andthe self-gravity corrections, as well as a sufficient spatial resolution and the calculation of disk-planet attraction both inside and outside the Hill sphere. With this proviso, we find Type III migration to be a robust mode of migration, astrophysically promising because of a speed much faster than in the previously studied modes of migration.Comment: 17 pages, 15 figures, submitted to MNRAS. Comments welcom

Quasi-Periodic Formaldehyde Maser Flares in the Massive Protostellar Object IRAS18566+0408

We report results of an extensive observational campaign of the 6 cm formaldehyde maser in the young massive stellar object IRAS18566+0408 (G37.55+0.20) conducted from 2002 to 2009. Using Arecibo, VLA, and GBT, we discovered quasi-periodic formaldehyde flares (P ~ 237 days). Based on Arecibo observations, we also discovered correlated variability between formaldehyde (H2CO) and methanol (CH3OH) masers. The H2CO and CH3OH masers are not spatially coincident, as demonstrated by different line velocities and high angular resolution MERLIN observations. The flares could be caused by variations in the infrared radiation field, possibly modulated by periodic accretion onto a young binary system.Comment: 16 pages, 3 figures, accepted for publication in the Astrophysical Journal Letter

Pumping of a Planetesimal Disc by a Rapidly Migrating Planet

We examine the effect of a rapidly migrating protoplanet on a ring of planetesimals. The eccentricities of the planetesimals are usually increased by $\Delta e \in (0.01, 0.1)$, with the exact increase being proportional to the protoplanet's mass, and inversely proportional to its migration rate. The eccentricity distribution is also substantially changed from a Rayleigh distribution. We discuss the possible implications for further planet formation, and suggest that the rapid passage of a protoplanet may not prevent the planetesimal disc from forming further planets.Comment: Five pages, two figures, accepted by MNRA

Limiting eccentricity of sub-parsec massive black hole binaries surrounded by self-gravitating gas discs

We study the dynamics of supermassive black hole binaries embedded in circumbinary gaseous discs, with the SPH code Gadget-2. The sub-parsec binary (of total mass M and mass ratio q=1/3) has excavated a gap and transfers its angular momentum to the self--gravitating disc (M_disc=0.2 M). We explore the changes of the binary eccentricity e, by simulating a sequence of binary models that differ in the initial eccentricity e_0, only. In initially low-eccentric binaries, the eccentricity increases with time, while in high-eccentric binaries e declines, indicating the existence of a limiting eccentricity e_crit that is found to fall in the interval [0.6,0.8]. We also present an analytical interpretation for this saturation limit. An important consequence of the existence of e_crit is the detectability of a significant residual eccentricity e_LISA} by the proposed gravitational wave detector LISA. It is found that at the moment of entering the LISA frequency domain e_LISA ~ 10^{-3}-10^{-2}; a signature of its earlier coupling with the massive circumbinary disc. We also observe large periodic inflows across the gap, occurring on the binary and disc dynamical time scales rather than on the viscous time. These periodic changes in the accretion rate (with amplitudes up to ~100%, depending on the binary eccentricity) can be considered a fingerprint of eccentric sub-parsec binaries migrating inside a circumbinary disc.Comment: 10 pages, 7 figures, accepted for publication in MNRA

Modelling Circumbinary Gas Flows in Close T Tauri Binaries

Young close binaries open central gaps in the surrounding circumbinary accretion disc, but the stellar components may still gain mass from gas crossing through the gap. It is not well understood how this process operates and how the stellar components are affected by such inflows. Our main goal is to investigate how gas accretion takes place and evolves in close T Tauri binary systems. In particular, we model the accretion flows around two close T Tauri binaries, V4046 Sgr and DQ Tau, both showing periodic changes in emission lines, although their orbital characteristics are very different. In order to derive the density and velocity maps of the circumbinary material, we employ two-dimensional hydrodynamic simulations with a locally isothermal equation of state. The flow patterns become quasi-stable after a few orbits in the frame co-rotating with the system. Gas flows across the circumbinary gap through the co-rotating Lagrangian points, and local circumstellar discs develop around both components. Spiral density patterns develop in the circumbinary disc that transport angular momentum efficiently. Mass is preferentially channelled towards the primary and its circumstellar disc is more massive than the disc around the secondary. We also compare the derived density distribution to observed line profile variability. The line profile variability tracing the gas flows in the central cavity shows clear similarities with the corresponding observed line profile variability in V4046 Sgr, but only when the local circumstellar disc emission was excluded. Closer to the stars normal magnetospheric accretion may dominate while further out the dynamic accretion process outlined here dominates. Periodic changes in the accretion rates onto the stars can explain the outbursts of line emission observed in eccentric systems such as DQ Tau.Comment: Accepted for publication in MNRA

Eccentricity Evolution of Extrasolar Multiple Planetary Systems due to the Depletion of Nascent Protostellar Disks

Most extrasolar planets are observed to have eccentricities much larger than those in the solar system. Some of these planets have sibling planets, with comparable masses, orbiting around the same host stars. In these multiple planetary systems, eccentricity is modulated by the planets' mutual secular interaction as a consequence of angular momentum exchange between them. For mature planets, the eigenfrequencies of this modulation are determined by their mass and semi-major axis ratios. But, prior to the disk depletion, self gravity of the planets' nascent disks dominates the precession eigenfrequencies. We examine here the initial evolution of young planets' eccentricity due to the apsidal libration or circulation induced by both the secular interaction between them and the self gravity of their nascent disks. We show that as the latter effect declines adiabatically with disk depletion, the modulation amplitude of the planets' relative phase of periapse is approximately invariant despite the time-asymmetrical exchange of angular momentum between planets. However, as the young planets' orbits pass through a state of secular resonance, their mean eccentricities undergo systematic quantitative changes. For applications, we analyze the eccentricity evolution of planets around Upsilon Andromedae and HD168443 during the epoch of protostellar disk depletion. We find that the disk depletion can change the planets' eccentricity ratio. However, the relatively large amplitude of the planets' eccentricity cannot be excited if all the planets had small initial eccentricities.Comment: 50 pages including 11 figures, submitted to Ap

Global m=1 modes and migration of protoplanetary cores in eccentric protoplanetary discs

We calculate global $m=1$ modes with low pattern speed corresponding to introducing a finite eccentricity into a protoplanetary disc. We consider disc models which are either isolated or contain one or two protoplanets orbiting in an inner cavity. Global modes that are strongly coupled to inner protoplanets are found to have disc orbits which tend to have apsidal lines antialigned with respect to those of the inner protoplanets. Other modes corresponding to free disc modes may be global over a large range of length scales and accordingly be long lived. We consider the motion of a protoplanet in the earth mass range embedded in an eccentric disc and determine the equilibrium orbits which maintain fixed apsidal alignment with respect to the disc gas orbits. Equilibrium eccentricities are found to be comparable or possibly exceed the disc eccentricity. We then approximately calculate the tidal interaction with the disc in order to estimate the orbital migration rate. Results are found to deviate from the case of axisymmetric disc with near circular protoplanet orbit once eccentricities of protoplanet and disc orbits become comparable to the disc aspect ratio in magnitude. Aligned protoplanet orbits with very similar eccentricity to that of the gas disc are found to undergo litle eccentricity change while undergoing inward migration in general. However, for significantly larger orbital eccentricities, migration may be significantly reduced or even reverse from inwards to outwards. Thus the existence of global non circular motions in discs with radial excursions comparable to the semi-thickness may have important consequences for the migration and survival of protoplanetary cores in the earth mass range.Comment: Accepted for publication by A &

Narrow-Angle Astrometry with the Space Interferometry Mission: The Search for Extra-Solar Planets. II. Detection and Characterization of Planetary Systems

(Abridged) The probability of detecting additional companions is essentially unchanged with respect to the single-planet configurations, but after fitting and subtraction of orbits with astrometric signal-to-noise ratio $\alpha/\sigma_d\to 1$ the false detection rates can be enhanced by up to a factor 2; the periodogram approach results in robust multiple-planet detection for systems with periods shorter than the SIM mission length, even at low values of $\alpha/\sigma_d$, while the least squares technique combined with Fourier series expansions is arguably preferable in the long-period regime. The accuracy on multiple-planet orbit reconstruction and mass determination suffers a typical degradation of 30-40% with respect to single-planet solutions; mass and orbital inclination can be measured to better than 10% for periods as short as 0.1 yr, and for $\alpha/\sigma_d$ as low as $\sim 5$, while $\alpha/\sigma_d\simeq 100$ is required in order to measure with similar accuracy systems harboring objects with periods as long as three times the mission duration. For systems with all components producing $\alpha/\sigma_d\simeq 10$ or greater, quasi-coplanarity can be reliably established with uncertainties of a few degrees, for periods in the range $0.1\leq T\leq 15$ yr; in systems where at least one component has $\alpha/\sigma_d\to 1$, coplanarity measurements are compromised, with typical uncertainties on the mutual inclinations of order of $30^\circ-40^\circ$. Our findings are illustrative of the importance of the contribution SIM will make to the fields of formation and evolution of planetary systems.Comment: 61 pages, 14 figures, 5 tables, to appear in the September 2003 Issue of the Publications of the Astronomical Society of the Pacifi

Twins Among the Low Mass Spectroscopic Binaries

We report an analysis of twins of spectral types F or later in the 9th Catalog of Spectroscopic Binaries (SB9). Twins, the components of binaries with mass ratio within 2% of 1.0, are found among the binaries with primaries of F and G spectral type. They are most prominent among the binaries with periods less than 43 days, a cutoff first identified by Lucy. Within the subsample of binaries with P<43 days, the twins do not differ from the other binaries in their distributions of periods (median P~7d), masses, or orbital eccentricities. Combining the mass ratio distribution in the SB9 in the mass range 0.6 to 0.85 Msun with that measured by Mazeh et al. for binaries in the Carney-Latham high proper motion survey, we estimate that the frequency of twins in a large sample of spectroscopic binaries is about 3%. Current theoretical understanding indicates that accretion of high specific angular momentum material by a protobinary tends to equalize its masses. We speculate that the excess of twins is produced in those star forming regions where the accretion processes were able to proceed to completion for a minority of protobinaries. This predicts that the components of a young twin may appear to differ in age and that, in a sample of spectroscopic binaries in a star formation region, the twins are, on average, older than the binaries with mass ratios much smaller than 1.Comment: Accepted by the Astronomical Journa