No Arabic abstract
Context. More than 60 planets have been discovered so far in systems that harbour two stars, some of which have binary semi-major axes as small as 20 au. It is well known that the formation of planets in such systems is strongly influenced by the stellar components, since the protoplanetary disc and the particles within are exposed to the gravitational influence of the binary. However, the question on how self-gravitating protoplanetary bodies affect the evolution of a radiative, circumprimary disc is still open. Aims. We present our 2D hydrodynamical GPU-CPU code and study the interaction of several thousands of self-gravitating particles with a viscous and radiative circumprimary disc within a binary star system. To our knowledge this program is the only one at the moment that is capable to handle this many particles and to calculate their influence on each other and on the disc. Methods. We performed hydrodynamical simulations of a circumstellar disc assuming the binary system to be coplanar. Our gridbased staggered mesh code relies on ideas from ZEUS-2D, where we implemented the FARGO algorithm and an additional energy equation for the radiative cooling according to opacity tables. To treat particle motion we used a parallelised version of the precise Bulirsch - Stoer algorithm. Four models in total where computed taking into account (i) only N-body interaction, (ii) N-body and disc interaction, (iii) the influence of computational parameters (especially smoothing) on N-body interaction, and (iv) the influence of a quiet low-eccentricity disc while running model (ii). The impact velocities where measured at two different time intervals and were compared. Results. We show that the combination of disc- and N-body self-gravity can have a significant influence on the orbit evolution of roughly Moon sized protoplanets.
Discs in binaries have a complex behavior because of the perturbations of the companion star. Planet formation in binary-star systems both depend on the companion star parameters and on the properties of the circumstellar disc. An eccentric disc may increase the impact velocity of planetesimals and therefore jeopardize the accumulation process. We model the evolution of discs in close binaries including the effects of self-gravity and adopting different prescriptions to model the discs radiative properties. We focus on the dynamical properties and evolutionary tracks of the discs. We use the hydrodynamical code FARGO and we include in the energy equation heating and cooling effects. Radiative discs have a lower disc eccentricity compared to locally isothermal discs with same temperature profile. As a consequence, we do not observe the formation of an internal elliptical low density region as in locally isothermal disc models. However, the disc eccentricity depends on the disc mass through the opacities. Akin to locally isothermal disc models, self-gravity forces the discs longitude of pericenter to librate about a fixed orientation with respect to the binary apsidal line ($pi$). The discs radiative properties play an important role in the evolution of discs in binaries. A radiative disc has an overall shape and internal structure that are significantly different compared to a locally isothermal disc with same temperature profile. This is an important finding both for describing the evolutionary track of the disc during its progressive mass loss, and for planet formation since the internal structure of the disc is relevant for planetesimals growth in binary systems. The non-symmetrical distribution of mass in these discs causes large eccentricities for planetesimals that may affect their growth.
Horseshoe-shaped brightness asymmetries of several transitional discs are thought to be caused by large-scale vortices. Anticyclonic vortices are efficiently collect dust particles, therefore they can play a major role in planet formation. Former studies suggest that the disc self-gravity weakens vortices formed at the edge of the gap opened by a massive planet in discs whose masses are in the range of 0.01<=M_disc/M_*<=0.1. Here we present an investigation on the long-term evolution of the large-scale vortices formed at the viscosity transition of the discs dead zone outer edge by means of two-dimensional hydrodynamic simulations taking disc self-gravity into account. We perform a numerical study of low mass, 0.001<=M_disc/M_*<=0.01, discs, for which cases disc self-gravity was previously neglected. The large-scale vortices are found to be stretched due to disc self-gravity even for low-mass discs with M_disc/M_*>=0.005 where initially the Toomre Q-parameter was <=50 at the vortex distance. As a result of stretching, the vortex aspect ratio increases and a weaker azimuthal density contrast develops. The strength of the vortex stretching is proportional to the disc mass. The vortex stretching can be explained by a combined action of a non-vanishing gravitational torque caused by the vortex, and the Keplerian shear of the disc. Self-gravitating vortices are subject to significantly faster decay than non-self-gravitating ones. We found that vortices developed at sharp viscosity transitions of self-gravitating discs can be described by a GNG model as long as the disc viscosity is low, i.e. alpha_dz<=10^-5.
We investigate the role of secular and mean motion resonances on the water transport from a belt of icy asteroids onto planets or embryos orbiting inside the circumprimary habitable zone (HZ) of a binary star system. In addition, the host-star has an accompanying gas giant planet. For a comparison, we perform two case studies where a secular resonance (SR) is located either inside the HZ close to 1.0 au (causing eccentric motion of a planet or embryos therein) or in the asteroid belt, beyond the snow line. In the latter case, a higher flux of icy objects moving towards the HZ is expected. Collisions between asteroids and objects in the HZ are treated analytically. Our purely dynamical study shows that the SR in the HZ boosts the water transport; however, collisions can occur at very high impact speeds. In this paper, we treat for the first time, realistic collisions using a GPU 3D-SPH code to assess the water loss in the projectile. Including the water loss into the dynamical results, we get more realistic values for the water mass fraction of the asteroid during an impact. We highlight that collisions occurring at high velocities greatly reduce the water content of the projectile and thus the amount of water transported to planets or embryos orbiting inside the HZ. Moreover, we discuss other effects that could modify our results, namely the asteroids surface rate recession due to ice sublimation and the atmospheric drag contribution on the asteroids mass loss.
We investigate how the detectability of signatures of self-gravity in a protoplanetary disc depends on its temporal evolution. We run a one-dimensional model for secular timescales to follow the disc mass as a function of time. We then combine this with three-dimensional global hydrodynamics simulations that employ a hybrid radiative transfer method to approximate realistic heating and cooling. We simulate ALMA continuum observations of these systems, and find that structures induced by the gravitational instability (GI) are readily detectable when $q=M_mathrm{disc}/M_*gtrsim 0.25$ and $R_mathrm{outer}lesssim 100$ au. The high accretion rate generated by gravito-turbulence in such a massive disc drains its mass to below the detection threshold in $sim10^4$ years, or approximately 1 % of the typical disc lifetime. Therefore, discs with spiral arms detected in ALMA dust observations, if generated by self-gravity, must either be still receiving infall to maintain a high $q$ value, or have just emerged from their natal envelope. Detection of substructure in systems with lower $q$ is possible, but would require a specialist integration with the most extended configuration over several days. This disfavours the possibility of GI-caused spiral structure in systems with $q<0.25$ being detected in relatively short integration times, such as those found in the DSHARP ALMA survey (Andrews et al. 2018; Huang et al. 2018). We find no temporal dependence of detectability on dynamical timescales
We perform a population synthesis of protoplanetary discs including infall with a total of $50,000$ simulations using a 1D vertically integrated viscous evolution code, studying a large parameter space in final stellar mass. Initial conditions and infall locations are chosen based on the results from a radiation-hydrodynamic population synthesis of circumstellar discs. We also consider a different infall prescription based on a magnetohydrodynamic (MHD) collapse simulation in order to assess the influence of magnetic fields on disc formation. The duration of the infall phase is chosen to produce a stellar mass distribution in agreement with the observationally determined stellar initial mass function. We find that protoplanetary discs are very massive early in their lives. When averaged over the entire stellar population, the discs have masses of $sim 0.3$ and $0.1,mathrm{M_odot}$ for systems based on hydrodynamic or MHD initial conditions, respectively. In systems with final stellar mass $sim 1,mathrm{M_odot}$, we find disc masses of $sim 0.7,mathrm{M_odot}$ for the `hydro case and $sim 0.2,mathrm{M_odot}$ for the `MHD case at the end of the infall phase. Furthermore, the inferred total disc lifetimes are long, $approx 5-7,mathrm{Myr}$ on average, despite our choice of a high value of $10^{-2}$ for the background viscosity $alpha$-parameter. In addition, fragmentation is common in systems that are simulated using hydrodynamic cloud collapse, with more fragments of larger mass formed in more massive systems. In contrast, if disc formation is limited by magnetic fields, fragmentation is suppressed entirely.