Dust entrainment in photoevaporative winds: The impact of X-rays

X-ray- and EUV- (XEUV-) driven photoevaporative winds acting on protoplanetary disks around young T-Tauri stars may crucially impact disk evolution, affecting both gas and dust distributions. We investigate the dust entrainment in XEUV-driven photoevaporative winds and compare our results to existing MHD and EUV-only models. For an X-ray luminosity of $L_X = 2 \cdot 10^{30}\,\mathrm{erg/s}$ emitted by a $M_* = 0.7\,\mathrm{M}_\odot$ star, corresponding to a wind mass-loss rate of $\dot{M}_\mathrm{w} \simeq 2.6 \cdot 10^{-8} \,\mathrm{M_\odot/yr}$, we find dust entrainment for sizes $a_0 \lesssim 11\,\mu$m ($9\,\mu$m) from the inner $25\,$AU ($120\,$AU). This is an enhancement over dust entrainment in less vigorous EUV-driven winds with $\dot{M}_\mathrm{w} \simeq 10^{-10}\,\mathrm{M_\odot/yr}$. Our numerical model also shows deviations of dust grain trajectories from the gas streamlines even for $\mu$m-sized particles. In addition, we find a correlation between the size of the entrained grains and the maximum height they reach in the outflow.Comment: Accepted for publication in A&A. 12+6 pages, 10+9 figure

High-resolution [OI] line spectral mapping of TW Hya consistent with X-ray driven photoevaporation

Theoretical models indicate that photoevaporative and magnetothermal winds play a crucial role in the evolution and dispersal of protoplanetary disks and affect the formation of planetary systems. However, it is still unclear what wind-driving mechanism is dominant or if both are at work, perhaps at different stages of disk evolution. Recent spatially resolved observations by Fang et al. (2023) of the [OI] 6300 Angstrom spectral line, a common disk wind tracer, in TW Hya revealed that about 80% of the emission is confined to the inner few au of the disk. In this work, we show that state-of-the-art X-ray driven photoevaporation models can reproduce the compact emission and the line profile of the [OI] 6300 Angstrom line. Furthermore, we show that the models also simultaneously reproduce the observed line luminosities and detailed spectral profiles of both the [OI] 6300 Angstrom and the [NeII] 12.8 micron lines. While MHD wind models can also reproduce the compact radial emission of the [OI] 6300 Angstrom line, they fail to match the observed spectral profile of the [OI] 6300 Angstrom line and underestimate the luminosity of the [NeII] 12.8 micron line by a factor of three. We conclude that, while we cannot exclude the presence of an MHD wind component, the bulk of the wind structure of TW Hya is predominantly shaped by a photoevaporative flow.Comment: 7 pages, 4 figures, accepted for publication in Astrophysical Journal Letter

Influence of the circumbinary disk gravity on planetesimal accumulation in the Kepler 16 system

Recent observations from NASA's Kepler mission detected the first planets in circumbinary orbits. The question we try to answer is where these planets formed in the circumbinary disk and how far inside they migrated to reach their present location. We investigate the first and more delicate phase of planet formation when planetesimals accumulate to form planetary embryos. We use the hydrodynamical code FARGO to study the evolution of the disk and of a test population of planetesimals embedded in it. With this hybrid hydrodynamical--N--body code we can properly account for the gas drag force on the planetesimals and for the gravitational force of the disk on them. The numerical simulations show that the gravity of the eccentric disk on the planetesimal swarm excites their eccentricities to values much larger than those induced by the binary perturbations only within 10 AU from the stars. Moreover, the disk gravity prevents a full alignment of the planetesimal pericenters. Both these effects lead to large impact velocities, beyond the critical value for erosion. Planetesimals accumulation in circumbinary disks appears to be prevented close to the stellar pair by the gravitational perturbations of the circumbinary disk. The observed planets possibly formed in the outer regions of the disk and then migrated inside by tidal interaction with the disk.Comment: Accepted for publication in A&

The dispersal of protoplanetary discs - I. A new generation of X-ray photoevaporation models

Photoevaporation of planet-forming discs by high-energy radiation from the central star is potentially a crucial mechanism for disc evolution and it may play an important role in the formation and evolution of planetary systems. We present here a new generation of X-ray photoevaporation models for solar-type stars, based on hydrodynamical simulations, which account for stellar irradiation via a significantly improved parametrization of gas temperatures, based on detailed photoionization and radiation transfer calculations. This is the first of a series of papers aiming at providing a library of models which cover the observed parameter space in stellar and disc mass, metallicity, and stellar X-ray properties. We focus here on solar-type stars (0.7 M⊙) with relatively low-mass discs (1 per cent of the stellar mass) and explore the dependence of the wind mass-loss rates on stellar X-ray luminosity. We model primordial discs and transition discs at various stages of evolution. Our two-dimensional hydrodynamical models are then used to derive simple recipes for the mass-loss rates that are suitable for one-dimensional disc evolution and/or planet formation models typically employed for population synthesis studies. Line profiles from typical wind diagnostics ([O I]6300 Å and [Ne II]12.8 μm) are also calculated for our models and found to be roughly in agreement with previous studies. Finally, we perform a population study of transition discs by means of one-dimensional viscous evolution models including our new photoevaporation prescription and find that roughly a half of observed transition discs cavities and accretion rates could be reproduced by our models

Detectability of embedded protoplanets from hydrodynamical simulations

We predict magnitudes for young planets embedded in transition discs, still affected by extinction due to material in the disc. We focus on Jupiter-size planets at a late stage of their formation, when the planet has carved a deep gap in the gas and dust distributions and the disc starts being transparent to the planet flux in the infrared (IR). Column densities are estimated by means of three-dimensional hydrodynamical models, performed for several planet masses. Expected magnitudes are obtained by using typical extinction properties of the disc material and evolutionary models of giant planets. For the simulated cases located at $5.2$ AU in a disc with local unperturbed surface density of $127$ $\mathrm{g} \cdot \mathrm{cm}^{-2}$, a $1$ $M_J$ planet is highly extincted in J-, H- and K-bands, with predicted absolute magnitudes $\ge 50$ mag. In L- and M-bands extinction decreases, with planet magnitudes between $25$ and $35$ mag. In the N-band, due to the silicate feature on the dust opacities, the expected magnitude increases to $40$ mag. For a $2$ $M_J$ planet, the magnitudes in J-, H- and K-bands are above $22$ mag, while for L-, M- and N-bands the planet magnitudes are between $15$ and $20$ mag. For the $5$ $M_J$ planet, extinction does not play a role in any IR band, due to its ability to open deep gaps. Contrast curves are derived for the transition discs in CQ Tau, PDS70, HL Tau, TW Hya and HD163296. Planet mass upper-limits are estimated for the known gaps in the last two systems.Comment: Accepted for publication on January 8, 2020 in MNRAS. 15 pages of main text with 14 figures, and 5 pages of appendices A and B with 4 figure

Interpreting molecular hydrogen and atomic oxygen line emission of T Tauri disks with photoevaporative disk-wind models

Context. Winds in protoplanetary disks play an important role in their evolution and dispersal. However, the physical process that is actually driving the winds is still unclear (i.e. magnetically versus thermally driven), and can only be understood by directly confronting theoretical models with observational data. Aims. We aim to interpret observational data for molecular hydrogen and atomic oxygen lines that show kinematic disk-wind signatures in order to investigate whether or not purely thermally driven winds are consistent with the data. Methods. We use hydrodynamic photoevaporative disk-wind models and post-process them with a thermochemical model to produce synthetic observables for the spectral lines o-H2 1-0 S(1) at 2.12 μm and [OI] 1D2-3P2 at 0.63 μm and directly compare the results to a sample of observations. Results. We find that our photoevaporative disk-wind model is consistent with the observed signatures of the blueshifted narrow low-velocity component (NLVC) -which is usually associated with slow disk winds -for both tracers. Only for one out of seven targets that show blueshifted NLVCs does the photoevaporative model fail to explain the observed line kinematics. Our results also indicate that interpreting spectral line profiles using simple methods, such as the thin-disk approximation, to determine the line emitting region is not appropriate for the majority of cases and can yield misleading conclusions. This is due to the complexity of the line excitation, wind dynamics, and the impact of the actual physical location of the line-emitting regions on the line profiles. Conclusions. The photoevaporative disk-wind models are largely consistent with the studied observational data set, but it is not possible to clearly discriminate between different wind-driving mechanisms. Further improvements to the models are necessary, such as consistent modelling of the dynamics and chemistry, and detailed modelling of individual targets (i.e. disk structure) would be beneficial. Furthermore, a direct comparison of magnetically driven disk-wind models to the observational data set is necessary in order to determine whether or not spatially unresolved observations of multiple wind tracers are sufficient to discriminate between theoretical models

Planet formation in Binaries

Spurred by the discovery of numerous exoplanets in multiple systems, binaries have become in recent years one of the main topics in planet formation research. Numerous studies have investigated to what extent the presence of a stellar companion can affect the planet formation process. Such studies have implications that can reach beyond the sole context of binaries, as they allow to test certain aspects of the planet formation scenario by submitting them to extreme environments. We review here the current understanding on this complex problem. We show in particular how each of the different stages of the planet-formation process is affected differently by binary perturbations. We focus especially on the intermediate stage of kilometre-sized planetesimal accretion, which has proven to be the most sensitive to binarity and for which the presence of some exoplanets observed in tight binaries is difficult to explain by in-situ formation following the "standard" planet-formation scenario. Some tentative solutions to this apparent paradox are presented. The last part of our review presents a thorough description of the problem of planet habitability, for which the binary environment creates a complex situation because of the presence of two irradation sources of varying distance.Comment: Review chapter to appear in "Planetary Exploration and Science: Recent Advances and Applications", eds. S. Jin, N. Haghighipour, W.-H. Ip, Springer (v2, numerous typos corrected

Circumstellar discs: What will be next?

This prospective chapter gives our view on the evolution of the study of circumstellar discs within the next 20 years from both observational and theoretical sides. We first present the expected improvements in our knowledge of protoplanetary discs as for their masses, sizes, chemistry, the presence of planets as well as the evolutionary processes shaping these discs. We then explore the older debris disc stage and explain what will be learnt concerning their birth, the intrinsic links between these discs and planets, the hot dust and the gas detected around main sequence stars as well as discs around white dwarfs.Comment: invited review; comments welcome (32 pages

Disks in close binary stars: γ-Cephei revisited

Context. Close binaries (abin < 20 au) are known to harbor planets, yet planet formation is unlikely to succeed in such systems. Studying the dynamics of disks in close binaries can help to understand how those planets could have formed. Aims. We study the impact that numerical and physical parameters have on the dynamics of disks in close binaries. We use the γ-Cephei system as an example and focus on disk quantities such as disk eccentricity and the precession rate as indicators for the dynamical state of the disks. Methods. We simulate disks in close binaries by performing two-dimensional radiative hydrodynamical simulations using a modified version of the FARGO code. First, we perform a parameter study for different numerical parameters to confirm that our results are robust. In the second part, we study the effects of different masses and different viscosities on the disks' dynamics. Results. Previous studies on radiative disks in close binaries used too low resolutions and too small simulation domains, which impacted the disk's dynamics. We find that radiative disks in close binaries, after an initialization phase, become eccentric with mean eccentricities between 0:06 and 0:27 and display a slow retrograde precession with periods ranging from 4-40Tbin which depends quadratically on the disk's mean aspect ratio. In general, the disks show a coherent, rigid precession which can be broken, however, by changes in the opacity law reducing the overall eccentricity of the disk

Radiation-hydrodynamical models of X-ray photoevaporation in carbon-depleted circumstellar discs

The so-called transition discs provide an important tool to probe various mechanisms that might influence the evolution of protoplanetary discs and therefore the formation of planetary systems. One of these mechanisms is photoevaporation due to energetic radiation from the central star, which can in principal explain the occurrence of discs with inner cavities like transition discs. Current models, however, fail to reproduce a subset of the observed transition discs, namely objects with large measured cavities and vigorous accretion. For these objects the presence of (multiple) giant planets is often invoked to explain the observations. In our work, we explore the possibility of X-ray photoevaporation operating in discs with different gas-phase depletion of carbon and show that the influence of photoevaporation can be extended in such low-metallicity discs. As carbon is one of the main contributors to the X-ray opacity, its depletion leads to larger penetration depths of X-rays in the disc and results in higher gas temperatures and stronger photoevaporative winds. We present radiation-hydrodynamical models of discs irradiated by internal X-ray + EUV radiation assuming carbon gas-phase depletions by factors of three, 10, and 100 and derive realistic mass-loss rates and profiles. Our analysis yields robust temperature prescriptions as well as photoevaporative mass-loss rates and profiles which may be able to explain a larger fraction of the observed diversity of transition discs