Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userEffects of disc midplane evolution on CO snowline location
116
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Temperature changes in the planet forming disc midplanes carry important physico-chemical consequences, such as the effect on the locations of the condensation fronts of molecules - the snowlines. Snowlines impose major chemical gradients and possibly foster grain growth. The aim of this paper is to understand how disc midplane temperature changes with gas and dust evolution, and identify trends that may influence planet formation or allow to constrain disc evolution observationally. We calculate disc temperature, hydrostatic equilibrium and dust settling in a mutually consistent way from a grid of disc models at different stages of gas loss, grain growth and hole opening. We find that the CO snowline location depends very strongly on disc properties. The CO snowline location migrates closer to the star for increasing degrees of gas dispersal and dust growth. Around a typical A type star, the snowline can be anywhere between several tens and a few hundred au, depending on the disc properties such as gas mass and grain size. In fact, gas loss is as efficient as dust evolution in settling discs, and flat discs may be gas-poor counterparts of flared discs. Our results, in the context of different pre-main sequence evolution of the luminosity in low- and intermediate-mass stars suggests very different thermal (and hence chemical) histories in these two types of discs. Discs of T Tauri stars settle and cool down while discs of Herbig Ae stars may remain rather warm throughout the pre-main sequence.
rate research
Read More
We investigate the impact of pre-main sequence stellar luminosity evolution on the thermal and chemical properties of disc midplanes. We create template disc models exemplifying initial conditions for giant planet formation for a variety of stellar masses and ages. These models include the 2D physical structure of gas as well as 1D chemical structure in the disc midplane. The disc temperature profiles are calculated using fully physically consistent radiative transfer models for stars between 0.5 and 3 Msun and ages up to 10 Myr. The resulting temperature profiles are used to determine how the chemical conditions in the mid-plane change over time. We therefore obtain gas and ice-phase abundances of the main carbon and oxygen carrier species. While the temperature profiles produced are not markedly different for the stars of different masses at early stages (<1 Myr), they start to diverge significantly beyond 2 Myr. Discs around stars with mass >1.5 Msun become warmer over time as the stellar luminosity increases, whereas low-mass stars decrease in luminosity leading to cooler discs. This has an observable effect on the location of the CO snowline, which is located >200 au in most models for a 3 Msun star, but is always within 80 au for 0.5 Msun star. The chemical compositions calculated show that a well defined stellar mass and age range exists in which high C/O gas giants can form. In the case of the exoplanet HR8799b, our models show it must have formed before the star was 1 Myr old.
CO is thought to be the main reservoir of volatile carbon in protoplanetary disks, and thus the primary initial source of carbon in the atmospheres of forming giant planets. However, recent observations of protoplanetary disks point towards low volatile carbon abundances in many systems, including at radii interior to the CO snowline. One potential explanation is that gas phase carbon is chemically reprocessed into less volatile species, which are frozen on dust grain surfaces as ice. This mechanism has the potential to change the primordial C/H ratio in the gas. However, current observations primarily probe the upper layers of the disk. It is not clear if the low volatile carbon abundances extend to the midplane, where planets form. We have run a grid of 198 chemical models, exploring how the chemical reprocessing of CO depends on disk mass, dust grain size distribution, temperature, cosmic ray and X-ray ionization rate, and initial water abundance. Building on our previous work focusing on the warm molecular layer, here we analyze the results for our grid of models in the disk midplane at 12 au. We find that either an ISM level cosmic ray ionization rate or the presence of UV photons due to a low dust surface density are needed to chemically reduce the midplane CO gas abundance by at least an order of magnitude within 1 Myr. In the majority of our models CO does not undergo substantial reprocessing by in situ chemistry and there is little change in the gas phase C/H and C/O ratios over the lifetime of the typical disk. However, in the small sub-set of disks where the disk midplane is subject to a source of ionization or photolysis, the gas phase C/O ratio increases by up to nearly 9 orders of magnitude due to conversion of CO into volatile hydrocarbons.
The interaction between a planet located in the inner region of a disc and the warped outer region is studied. We consider the stage of evolution after the planet has cleared-out a gap, so that the planetary orbit evolves only under the gravitational potential from the disc. We develop a secular analysis and compute the evolution of the orbital elements by solving Lagranges equations valid to second order in the eccentricity. We also perform numerical simulations with the full disc potential. In general, the interaction between the disc and the planet leads to the precession of the orbit. The orbital plane therefore becomes tilted relative to the discs inner parts, with no change in the eccentricity. When the inclination approaches 90 degrees, there is an instability and the eccentricity increases. In this case, both the inclination and the eccentricity develop large variations, with the orbit becoming retrograde. As the eccentricity reaches high values, we would expect tidal capture on a short orbit of the planet by the star to occur. This instability happens when the disc is severely warped, or if there is a significant amount of mass in a ring inclined by at least 45 degrees relative to the initial orbital plane. The inclination of the orbit does not depend on the semimajor axis nor on the planets mass. However, for a significant inclination to be generated on a timescale of at most a few Myr, the planet should be beyond the snow line. The process described here would therefore produce two distinct populations of inclined planets: one with objects beyond the snow line with at most moderate eccentricities, and another with objects on short circularized orbits.
[Abridged] Snowlines in protoplanetary disks play an important role in planet formation and composition. Since the CO snowline is difficult to observe directly with CO emission, its location has been inferred in several disks from spatially resolved ALMA observations of DCO+ and N2H+. N2H+ is considered to be a good tracer of the CO snowline based on astrochemical considerations predicting an anti-correlation between N2H+ and gas-phase CO. In this work, the robustness of N2H+ as a tracer of the CO snowline is investigated. A simple chemical network is used in combination with the radiative transfer code LIME to model the N2H+ distribution and corresponding emission in the disk around TW Hya. The assumed CO and N2 abundances, corresponding binding energies, cosmic ray ionization rate, and degree of large-grain settling are varied to determine the effects on the N2H+ emission and its relation to the CO snowline. For the adopted physical structure of the TW Hya disk and molecular binding energies for pure ices, the balance between freeze-out and thermal desorption predicts a CO snowline at 19 AU, corresponding to a CO midplane freeze-out temperature of 20 K. A model with a total, i.e. gas plus ice, CO abundance of 3e-6 with respect to H2 fits the position of the emission peak observed by Qi et al. 2013 for the TW Hya disk. However, the relationship between N2H+ and the CO snowline is more complicated than generally assumed: for the investigated parameters, the N2H+ column density peaks at least 5 AU outside the CO snowline. Moreover, the N2H+ emission can peak much further out, as far as ~50 AU beyond the snowline. Hence, chemical modeling, as done here, is necessary to derive a CO snowline location from N2H+ observations.
During the process of planet formation, the planet-discs interactions might excite (or damp) the orbital eccentricity of the planet. In this paper, we present two long ($tsim 3times 10^5$ orbits) numerical simulations: (a) one (with a relatively light disc, $M_{rm d}/M_{rm p}=0.2$) where the eccentricity initially stalls before growing at later times and (b) one (with a more massive disc, $M_{rm d}/M_{rm p}=0.65$) with fast growth and a late decrease of the eccentricity. We recover the well-known result that a more massive disc promotes a faster initial growth of the planet eccentricity. However, at late times the planet eccentricity decreases in the massive disc case, but increases in the light disc case. Both simulations show periodic eccentricity oscillations superimposed on a growing/decreasing trend and a rapid transition between fast and slow pericentre precession. The peculiar and contrasting evolution of the eccentricity of both planet and disc in the two simulations can be understood by invoking a simple toy model where the disc is treated as a second point-like gravitating body, subject to secular planet-planet interaction and eccentricity pumping/damping provided by the disc. We show how the counterintuitive result that the more massive simulation produces a lower planet eccentricity at late times can be understood in terms of the different ratios of the disc-to-planet angular momentum in the two simulations. In our interpretation, at late times the planet eccentricity can increase more in low-mass discs rather than in high-mass discs, contrary to previous claims in the literature.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
