Disks In Binary Systems Research Articles

Sites of Planet Formation in Binary Systems. I. Evidence for Disk−Orbit Alignment in the Close Binary FO Tau

Close binary systems present challenges to planet formation. As binary separations decrease, so do the occurrence rates of protoplanetary disks in young systems and planets in mature systems. For systems that do retain disks, their disk masses and sizes are altered by the presence of the binary companion. Through the study of protoplanetary disks in binary systems with known orbital parameters, we seek to determine the properties that promote disk retention and therefore planet formation. In this work, we characterize the young binary−disk system FO Tau. We determine the first full orbital solution for the system, finding masses of 0.35−0.05+0.06M⊙ and 0.34 ± 0.05 M ⊙ for the stellar components, a semimajor axis of 22(−1+2) au, and an eccentricity of 0.21(−0.03+0.04) . With long-baseline Atacama Large Millimeter/submillimeter Array interferometry, we detect 1.3 mm continuum and 12CO (J = 2–1) line emission toward each of the binary components; no circumbinary emission is detected. The protoplanetary disks are compact, consistent with being truncated by the binary orbit. The dust disks are unresolved in the image plane, and the more extended gas disks are only marginally resolved. Fitting the continuum and CO visibilities, we determine the inclination of each disk, finding evidence for alignment of the disk and binary orbital planes. This study is the first of its kind linking the properties of circumstellar protoplanetary disks to a precisely known binary orbit. In the case of FO Tau, we find a dynamically placid environment (coplanar, low eccentricity), which may foster its potential for planet formation.

Planet migration in accretion discs in binary systems

AbstractWe model evolution of exoplanets of S-type in close binary systems at the stage when the companion starts to lose mass via a slow stellar wind. At this stage an accretion disc is formed around the planets’ host. Detailed structure of such discs is calculated in quasi-stationary and non-stationary approaches. We model migration of planets embedded in these discs.

On dust evolution in planet-forming discs in binary systems – II. Comparison with Taurus and ρ Ophiuchus (sub-)millimetre observations: discs in binaries have small dust sizes

ABSTRACT The recently discovered exoplanets in binary or higher order multiple stellar systems sparked a new interest in the study of protoplanetary discs in stellar aggregations. Here, we focus on disc solids, as they make up the reservoir out of which exoplanets are assembled and dominate (sub-)millimetre disc observations. These observations suggest that discs in binary systems are fainter and smaller than in isolated systems. In addition, disc dust sizes are consistent with tidal truncation only if they orbit very eccentric binaries. In a previous study, we showed that the presence of a stellar companion hastens the radial migration of solids, shortening disc lifetime, and challenging planet formation. In this paper, we confront our theoretical and numerical results with observations: Disc dust fluxes and sizes from our models are computed at ALMA wavelengths and compared with Taurus and ρ Ophiuchus data. A general agreement between theory and observations is found. In particular, we show that the dust disc sizes are generally smaller than the binary truncation radius due to the combined effect of grain growth and radial drift: Therefore, small disc sizes do not require implausibly high eccentricities to be explained. Furthermore, the observed binary discs are compatible within 1σ with a quadratic flux–radius correlation similar to that found for single-star discs and show a close match with the models. However, the observational sample of resolved binary discs is still small and additional data are required to draw more robust conclusions on the flux–radius correlation and how it depends on the binary properties.

On dust evolution in planet-forming discs in binary systems – I. Theoretical and numerical modelling: radial drift is faster in binary discs

ABSTRACT Many stars are in binaries or higher order multiple stellar systems. Although in recent years a large number of binaries have been proven to host exoplanets, how planet formation proceeds in multiple stellar systems has not been studied much yet from the theoretical standpoint. In this paper, we focus on the evolution of the dust grains in planet-forming discs in binaries. We take into account the dynamics of gas and dust in discs around each component of a binary system under the hypothesis that the evolution of the circumprimary and the circumsecondary discs is independent. It is known from previous studies that the secular evolution of the gas in binary discs is hastened due to the tidal interactions with their hosting stars. Here, we prove that binarity affects dust dynamics too, possibly in a more dramatic way than the gas. In particular, the presence of a stellar companion significantly reduces the amount of solids retained in binary discs because of a faster, more efficient radial drift, ultimately shortening their lifetime. We prove that how rapidly discs disperse depends both on the binary separation, with discs in wider binaries living longer, and on the disc viscosity. Although the less-viscous discs lose high amounts of solids in the earliest stages of their evolution, they are dissipated slowly, while those with higher viscosities show an opposite behaviour. The faster radial migration of dust in binary discs has a striking impact on planet formation, which seems to be inhibited in this hostile environment, unless other disc substructures halt radial drift further in. We conclude that if planetesimal formation were viable in binary discs, this process would take place on very short time-scales.

A Fast-growing Tilt Instability of Detached Circumplanetary Disks

Abstract 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 (phantom) and grid-based (athena++) hydrodynamic simulations. For a disk with outer aspect ratio H/r ≃ 0.1, initially moderate tilts double on a timescale 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.

Hydrodynamical Simulations of Misaligned Accretion Discs in Binary Systems: Companions tear discs

AbstractAccretion discs appear in many astrophysical systems. In most cases, these discs are probably not completely axisymmetric. Discs in binary systems are often found to be misaligned with respect to the binary orbit. In this case, the gravitational torque from a companion induces nodal precession in misaligned rings of gas. We first calculate whether this precession is strong enough to overcome the internal disc torques communicating angular momentum. For typical parameters, precession torque wins. To check this result, we perform numerical simulations using the Smoothed Particle Hydrodynamics code, PHANTOM, and confirm that sufficiently thin and sufficiently inclined discs can break into distinct planes that precess effectively independently. Disc tearing is widespread and severely changes the disc structure. It enhances dissipation and promotes stronger accretion onto the central object. We also perform a stability analysis on isolated warped discs to understand the physics of disc breaking and tearing observed in numerical simulations. The instability appears in the form of viscous anti-diffusion of the warp amplitude and the surface density. The discovery of disc breaking and tearing has revealed new physical processes that dramatically change the evolution of accretion discs, with obvious implications for observed systems.

Instability of non-Keplerian warped discs

ABSTRACT Many accretion discs are thought to be warped. Recent hydrodynamical simulations show that (i) discs can break into distinct planes when the amplitude of an imposed warp is sufficiently high and the viscosity sufficiently low, and that (ii) discs can tear up into discrete rings when an initially planar disc is subject to a forced precession. Previously, we investigated the local stability of isolated, Keplerian, warped discs in order to understand the physics causing an accretion disc to break into distinct planes, finding that antidiffusion of the warp amplitude is the underlying cause. Here, we explore the behaviour of this instability in disc regions where the rotation profile deviates from Keplerian. We find that at small warp amplitudes non-Keplerian rotation can stabilize the disc by increasing the critical warp amplitude for instability, while at large warp amplitudes non-Keplerian rotation can lead to an increased growth rate for discs that are unstable. Tidal effects on discs in binary systems are typically weak enough such that the disc remains close to Keplerian rotation. However, the inner regions of discs around black holes are strongly affected, with the smallest radius at which the disc can break into discrete planes being a function of the black hole spin. We suggest that interpreting observed frequencies in the power spectra of light curves from accreting compact objects as nodal and apsidal precession of discrete orbits requires an instability that can break the disc into discrete rings such as the one explored here.

ALMA 0.88 mm Survey of Disks around Planetary-mass Companions

Abstract Characterizing the physical properties and compositions of circumplanetary disks can provide important insights into the formation of giant planets and satellites. We report Atacama Large Millimeter/submillimeter Array 0.88 mm (Band 7) continuum observations of six planetary-mass (10–20 M Jup) companions: CT Cha b, 1RXS 1609 b, ROXs 12 b, ROXs 42B b, DH Tau b, and FU Tau b. No continuum sources are detected at the locations of the companions down to 3σ limits of 120–210 μJy. Given these nondetections, it is not clear whether disks around planetary-mass companions indeed follow the disk-flux–host-mass trend in the stellar regime. The faint radio brightness of these companion disks may result from a combination of fast radial drift and a lack of dust traps. Alternatively, as disks in binary systems are known to have significantly lower millimeter fluxes due to tidal interactions, these companion disks may instead follow the relationship of moderate-separation binary stars. This scenario can be tested with sensitive continuum imaging at rms levels of ≲10 μJy.

Self-gravitating disks in binary systems: an SPH approach

The study of the stability of massive gaseous disks around a star in a nonisolated context is a difficult task and becomes even more complicated for disks that are hosted by binary systems. The role of self-gravity is thought to be significant when the ratio of the disk-to-star mass is non-negligible. To solve these problems, we implemented, tested, and applied our own smoothed particle hydrodynamics (SPH) algorithm. The code (named GaSPH) passed various quality tests and shows good performances, and it can therefore be reliably applied to the study of disks around stars when self-gravity needs to be accounted for. We here introduce and describe the algorithm, including some performance and stability tests. This paper is the first part of a series of studies in which self-gravitating disks in binary systems are let evolve in larger environments such as open clusters.

Planet migration in wind-fed accretion discs in binaries

Planet migration originally refers to protoplanetary disks, which are more massive and dense than typical accretion disks in binary systems. We study planet migration in an accretion disk in a binary system consisting of a solar-like star hosting a planet and a red giant donor star. The accretion disk is fed by a stellar wind. %, disk self-gravity is neglected. We use the $\alpha$-disk model and consider that the stellar wind is time-dependent. Assuming the disk is quasi-stationary we calculate its temperature and surface density profiles. In addition to the standard disk model, when matter is captured by the disk at its outer edge, we study the situation when the stellar wind delivers matter on the whole disc surface inside the accretion radius with the rate depending on distance from the central star. Implying that a planet experiences classical type I/II migration we calculate migration time for a planet on a circular orbit coplanar with the disk. Potentially, rapid inward planet migration can result in a planet-star merger which can be accompanied by an optical or/and UV/X-ray transient. We calculate timescale of migration for different parameters of planets and binaries. Our results demonstrate that planets can fall on their host stars within the lifetime of the late-type donor for realistic sets of parameters.

Simulating the synchrotron emission of AGN with grid based relativistic hydrodynamics

Grid based numerical hydrodynamic codes have been used to simulate the evolution of a variety of astrophysical environments, ranging from stellar winds to accretion disks in binary systems. These simulations can be used as valuable tools to correlate the theoretical models and observational data. In this study a grid based relativistic hydrodynamic simulation is created, using the freely available numerical code PLUTO ver 4.2, to study possible causes of variability in AGN. The code uses shock-capturing Godunov-type methods to solve the fluid dynamic conservation equations on a structured mesh grid. To simulate the production of a jet similar to those observed in AGN a three dimensional static grid was created containing 256 × 256 × 512 computational cells. A uniform background medium was assigned to this grid, while a nozzle was defined on the z = 0 boundary to inject a pressure matched relativistic jet, with Γ = 10 (v = 0:995c), into the medium. The simulation showed the formation of a central collimated relativistic beam along with a surrounding cocoon region. In addition, a post-processing emission modelling code is being developed, in order to determine the emission produced by the relativistic jet. The emission model utilizes the data produced by the hydrodynamic code to calculate intensity maps due to synchrotron self absorption. The emission code takes into account the effects produced by the relativistic motion of the emitting regions (excluding light travel time) which can lead to effects such as the Doppler boosting observed in blazars. Synchrotron based emission models can be applied to other astrophysical objects, such as X-ray binary systems, which means that the application of the emission code can be adapted to model a variety of astrophysical sources.

Be discs in binary systems – II. Misaligned orbits

We use a smoothed particle hydrodynamics (SPH) code to examine the effects of misaligned binary companions on Be star discs. We systematically vary the degree of misalignment between the disc and the binary orbit, as well as the disc viscosity and orbital period to study their effects on the density in the inner and outer parts of the disc. We find that varying the degree of misalignment, the viscosity, and the orbital period affects both the truncation radius and the density structure of the outer disc, while the inner disc remains mostly unaffected. We also investigate the tilting of the disc in the innermost part of the disc and find the tilt increases with radius until reaching a maximum around 5 stellar radii. The direction of the line of nodes, with respect to the equator of the central star, is found to be offset compared to the orbital line of nodes, and to vary periodically in time, with a period of half a orbital phase. We also compare the scale height of our discs with the analytical scale height of an isothermal disc, which increases with radius as $r^{1.5}$. We find that this formula reproduces the scale height well for both aligned and misaligned systems but underestimates the scale height in regions of the disc where density enhancements develop.

Excitation of turbulence in accretion disks of binary stars by non-linear perturbations

Accretion disks in binary systems can experience hydrodynamic impact at inner as well as outer edges. The first case is typical for protoplanetary disks around young T Tau stars. The second one is typical for circumstellar disks in close binaries. As a result of such an impact, perturbations with different scales and amplitudes are excited in the disk. We investigated the nonlinear evolution of perturbations of a finite, but small amplitude, at the background of sub-Keplerian flow. Nonlinear effects at the front of perturbations lead to the formation of a shock wave, namely the discontinuity of the density and radial velocity. At this, the tangential flow in the neighborhood of the shock becomes equivalent to the flow in in the boundary layer. Instability of the tangential flow further leads to turbulization of the disk. Characteristics of the turbulence depend on perturbation parameters, but alpha-parameter of Shakura-Sunyaev does not exceed ~0.1.

Revealing the structure of the outer disks of Be stars

Context. The structure of the inner parts of Be star disks (20 stellar radii) is well explained by the viscous decretion disk (VDD) model, which is able to reproduce the observable properties of most of the objects studied so far. The outer parts, on the ther hand, are not observationally well-explored, as they are observable only at radio wavelengths. A steepening of the spectral slope somewhere between infrared and radio wavelengths was reported for several Be stars that were previously detected in the radio, but a convincing physical explanation for this trend has not yet been provided. Aims. We test the VDD model predictions for the extended parts of a sample of six Be disks that have been observed in the radio to address the question of whether the observed turndown in the spectral energy distribution (SED) can be explained in the framework of the VDD model, including recent theoretical development for truncated Be disks in binary systems. Methods. We combine new multi-wavelength radio observations from the Karl. G. Jansky Very Large Array (JVLA) and Atacama Pathfinder Experiment (APEX) with previously published radio data and archival SED measurements at ultraviolet, visual, and infrared wavelengths. The density structure of the disks, including their outer parts, is constrained by radiative transfer modeling of the observed spectrum using VDD model predictions. In the VDD model we include the presumed effects of possible tidal influence from faint binary companions. Results. For 5 out of 6 studied stars, the observed SED shows strong signs of SED turndown between far-IR and radio wavelengths. A VDD model that extends to large distances closely reproduces the observed SEDs up to far IR wavelengths, but fails to reproduce the radio SED. ... (abstract continues but did not fit here)

Kozai–Lidov disc instability

Recent results by Martin et al. (2014) showed in 3D SPH simulations that tilted discs in binary systems can be unstable to the development of global, damped Kozai-Lidov (KL) oscillations in which the discs exchange tilt for eccentricity. We investigate the linear stability of KL modes for tilted inviscid discs under the approximations that the disc eccentricity is small and the disc remains flat. By using 1D equations, we are able to probe regimes of large ratios of outer to inner disc edge radii that are realistic for binary systems of hundreds of AU separations and are not easily probed by multi-dimensional simulations. For order unity binary mass ratios, KL instability is possible for a window of disc aspect ratios H/r in the outer parts of a disc that roughly scale as (n_b/n)^2 < H/r < n_b/n, for binary orbital frequency n_b and orbital frequency n at the disc outer edge. We present a framework for understanding the zones of instability based on the determination of branches of marginally unstable modes. In general, multiple growing eccentric KL modes can be present in a disc. Coplanar apsidal-nodal precession resonances delineate instability branches. We determine the range of tilt angles for unstable modes as a function of disc aspect ratio. Unlike the KL instability for free particles that involves a critical (minimum) tilt angle, disc instability is possible for any nonzero tilt angle depending on the disc aspect ratio.

On the relationship between circumstellar disc size and X-ray outbursts in Be/X-ray binaries

We present long term H$\alpha$ monitoring results of five Be/X-ray binaries to study the Be disc size variations and their influence on Type II (giant) X-ray outbursts. The work is done in the context of the viscous decretion disc model which predicts that Be discs in binary systems are truncated by resonant torques induced by the neutron star in its orbit. Our observations show that type II outbursts are not correlated(nor anti-correlated) with the disc size, as they are seen to occur both at relatively small and large Be disc radii. We discuss these observations in context of alternate interpretation of Be disc behaviour, such as precession, elongation and density effects, and with cognisance of the limitations of our disc size estimates.

Be discs in binary systems – I. Coplanar orbits

Be stars are surrounded by outflowing circumstellar matter structured in the form of decretion discs. They are often members of binary systems, where it is expected that the decretion disc interacts both radiatively and gravitationally with the companion. In this work we study how various orbital (period, mass ratio and eccentricity) and disc (viscosity) parameters affect the disc structure in coplanar systems. We simulate such binaries with the use of a smoothed particle hydrodynamics code. The main effects of the secondary on the disc are its truncation and the accumulation of material inwards of truncation. We find two cases with respect to the effects of eccentricity: (i) In circular or nearly circular prograde orbits, the disc maintains a rotating, constant in shape, configuration, which is locked to the orbital phase. The disc is smaller in size, more elongated and more massive for low viscosity parameter, small orbital separation and/or high mass ratio. (ii) Highly eccentric orbits are more complex, with the disc structure and total mass strongly dependent on the orbital phase and the distance to the secondary. We also study the effects of binarity in the disc continuum emission. Since the infrared and radio SED are sensitive to the disc size and density slope, the truncation and matter accumulation result in considerable modifications in the emergent spectrum. We conclude that binarity can serve as an explanation for the variability exhibited in observations of Be stars, and that our model can be used to detect invisible companions.

WOBBLING AND PRECESSING JETS FROM WARPED DISKS IN BINARY SYSTEMS

We present results of the first ever three-dimensional (3D) magnetohydrodynamic (MHD) simulations of the accretion-ejection structure. We investigate the 3D evolution of jets launched symmetrically from single stars but also jets from warped disks in binary systems. We have applied various model setups and tested them by simulating a stable and bipolar symmetric 3D structure from a single star-disk-jet system. Our reference simulation maintains a good axial symmetry and also a bipolar symmetry for more than 600 rotations of the inner disk confirming the quality of our model setup. We have then implemented a 3D gravitational potential (Roche potential) due to a companion star and run a variety of simulations with different binary separations and mass ratios. These simulations show typical 3D deviations from axial symmetry, such as jet inclination outside the Roche lobe or spiral arms forming in the accretion disk. In order to find indication for precession effects, we have also run an exemplary parameter setup, essentially governed by a small binary separation of only $\simeq{200}$ inner disk radii. This simulation shows strong indication that we observe the onset of a jet precession caused by the wobbling of the jet launching disk. We estimate the opening angle of the precession cone defined by the lateral motion of the the jet axis of about 4 degree after about 5000 dynamical time steps.

Circumstellar discs in Galactic centre clusters: Disc-bearing B-type stars in the Quintuplet and Arches clusters

We investigate the circumstellar disc fraction as determined from L-band excess observations of the young, massive Arches and Quintuplet clusters residing in the central molecular zone of the Milky Way. The Quintuplet cluster was searched for L-band excess sources for the first time. We find a total of 26 excess sources in the Quintuplet cluster and 21 in the Arches cluster, of which 13 are new detections. With the aid of proper motion membership samples, the disc fraction of the Quintuplet cluster was derived for the first time to be 4.0 +/- 0.7%. There is no evidence for a radially varying disc fraction in this cluster. In the case of the Arches cluster, a disc fraction of 9.2 +/- 1.2% approximately out to the cluster's predicted tidal radius, r < 1.5 pc, is observed. This excess fraction is consistent with our previously found disc fraction in the cluster in the radial range 0.3 < r < 0.8 pc. In both clusters, the host star mass range covers late A- to early B-type stars, 2 < M < 15 Msun, as derived from J-band photospheric magnitudes. We discuss the unexpected finding of dusty circumstellar discs in these UV intense environments in the context of primordial disc survival and formation scenarios of secondary discs. We consider the possibility that the L-band excess sources in the Arches and Quintuplet clusters could be the high-mass counterparts to T Tauri pre-transitional discs. As such a scenario requires a long pre-transitional disc lifetime in a UV intense environment, we suggest that mass transfer discs in binary systems are a likely formation mechanism for the B-star discs observed in these starburst clusters.

Stellar multiplicity and debris discs: an unbiased sample

Circumstellar dust disks have been observed around many nearby stars. However, many stars are part of binary or multiple stellar systems. A natural question arises regarding the presence and properties of such disks in systems with more than one star. To address this, we consider a sample of 449 systems (spectral types A-M) observed with the Herschel Space Observatory as part of the DEBRIS program. We have examined the stellar multiplicity of this sample by gathering information from the literature and performing an adaptive optics imaging survey at Lick Observatory. Five new companions were revealed with our program. In total, we identify 188 (42%) binary or multiple star systems. The multiplicity of the sample is examined with regards to the detection of circumstellar disks for stars of spectral types AFGK. In general, disks are less commonly detected around binaries than single stars, though the disk frequency is comparable among A stars regardless of multiplicity. However, this sample reveals the period distribution of disk-bearing binaries is consistent with that of non-disk binaries and with comparison field samples. We find that the properties of disks in binary systems are not statistically different from those around single stars. Although the frequency of disk-bearing FGK binaries may be lower than in single star systems, the processes behind disk formation and the characteristics of these disks are comparable among both populations.