We study the influence of transport processes on the chemical evolution of DM Tau-like protoplanetary disks. Turbulent transport of gases and ices is implicitly modeled in full two dimensions (2D), ...using the mixing-length approximation, along with the time-dependent chemistry. We find that turbulent transport enhances abundances and column densities of many gas-phase species and ices, particularly, complex ones. The influence of turbulent mixing on disk chemistry is more pronounced in the inner, planet-forming disk region where gradients of temperature and high-energy radiation intensities are steeper than in the outer region. The molecules that are unresponsive to transport include, e.g., C2H, C+, CH4, CN, CO, HCN, HNC, H2CO, OH, as well as water and ammonia ice. Their column densities computed with the laminar and 2D mixing model differ by a factor of 2-5. Molecules whose vertical column densities in the laminar and dynamical models differ by up to two orders of magnitude include, e.g., C2H2, some carbon chains, CS, H2CS, H2O, HCO+, HCOOH, HNCO, N2H+, NH3, CO ice, H2CO ice, CH3OH ice, and electrons. Molecules whose column densities are altered by diffusion by more than two orders of magnitude include, e.g., C2S, C3S, C6H6, CO2, O2, SiO, SO, SO2, long carbon chain ices, CH3CHO ice, HCOOH ice, O2 ice, and OCN ice. We indicate several observable or potentially detectable tracers of transport processes in protoplanetary disks and the solar nebula, such as heavy hydrocarbon ices, complex organics, CO2, O2, SO, SO2, C2S, C3S compared to CO and water ice.
Context.
Clouds are ubiquitous in exoplanet atmospheres and they represent a challenge for the model interpretation of their spectra. When generating a large number of model spectra, complex cloud ...models often prove too costly numerically, whereas more efficient models may be overly simplified.
Aims.
We aim to constrain the atmospheric properties of the directly imaged planet HR 8799e with a free retrieval approach.
Methods.
We used our radiative transfer code petitRADTRANS for generating the spectra, which we coupled to the PyMultiNest tool. We added the effect of multiple scattering which is important for treating clouds. Two cloud model parameterizations are tested: the first incorporates the mixing and settling of condensates, the second simply parameterizes the functional form of the opacity.
Results.
In mock retrievals, using an inadequate cloud model may result in atmospheres that are more isothermal and less cloudy than the input. Applying our framework on observations of HR 8799e made with the GPI, SPHERE, and GRAVITY, we find a cloudy atmosphere governed by disequilibrium chemistry, confirming previous analyses. We retrieve that C/O = 0.60
−0.08
+0.07
. Other models have not yet produced a well constrained C/O value for this planet. The retrieved C/O values of both cloud models are consistent, while leading to different atmospheric structures: either cloudy or more isothermal and less cloudy. Fitting the observations with the self-consistent Exo-REM model leads to comparable results, without constraining C/O.
Conclusions.
With data from the most sensitive instruments, retrieval analyses of directly imaged planets are possible. The inferred C/O ratio of HR 8799e is independent of the cloud model and thus appears to be a robust. This C/O is consistent with stellar, which could indicate that the HR 8799e formed outside the CO
2
or CO iceline. As it is the innermost planet of the system, this constraint could apply to all HR 8799 planets.
Formation and evolution of water in the solar system and the origin of water on Earth constitute one of the most interesting questions in astronomy. In this paper, the isotopic and chemical evolution ...of water during the early history of the solar nebula, before the onset of planetesimal formation, is studied. A gas-grain chemical model that includes multiply deuterated species and nuclear spin-states is combined with a steady-state solar nebula model. To calculate initial abundances, we simulated 1 Myr of evolution of a cold and dark TMC-1-like prestellar core. Two time-dependent chemical models of the solar nebula are calculated over 1 Myr: (1) a laminar model and (2) a model with two-dimensional (2D) turbulent mixing. We find that with regards to the water isotopic composition and the origin of the comets, the mixing model seems to be favored over the laminar model.
Aims. We aim to obtain a spatially resolved measurement of velocity dispersions in the disk of TW Hya. Methods. We obtained images with high spatial and spectral resolution of the CO J = 2–1, CN N = ...2–1 and CS J = 5–4 emission with ALMA in Cycle 2. The radial distribution of the turbulent broadening was derived with two direct methods and one modelling approach. The first method requires a single transition and derives Tex directly from the line profile, yielding a vturb. The second method assumes that two different molecules are co-spatial, which allows using their relative line widths for calculating Tkin and vturb. Finally we fitted a parametric disk model in which the physical properties of the disk are described by power laws, to compare our direct methods with previous values. Results. The two direct methods were limited to the outer r > 40 au disk because of beam smear. The direct method found vturb to range from ≈130 m s-1 at 40 au, and to drop to ≈50 m s-1 in the outer disk, which is qualitatively recovered with the parametric model fitting. This corresponds to roughly 0.2−0.4 cs. CN was found to exhibit strong non-local thermal equilibrium effects outside r ≈ 140 au, so that vturb was limited to within this radius. The assumption that CN and CS are co-spatial is consistent with observed line widths only within r ≲ 100 au, within which vturb was found to drop from 100 m s-1 (≈0.4 cs) to zero at 100 au. The parametric model yielded a nearly constant 50 m s-1 for CS (0.2−0.4 cs). We demonstrate that absolute flux calibration is and will be the limiting factor in all studies of turbulence using a single molecule. Conclusions. The magnitude of the dispersion is comparable with or below that predicted by the magneto-rotational instability theory. A more precise comparison would require reaching an absolute calibration precision of about 3%, or finding a suitable combination of light and heavy molecules that are co-located in the disk.
The cross-disciplinary field of astrochemistry exists to understand the formation, destruction, and survival of molecules in astrophysical environments. Molecules in space are synthesized via a large ...variety of gas-phase reactions, and reactions on dust-grain surfaces, where the surface acts as a catalyst. A broad consensus has been reached in the astrochemistry community on how to suitably treat gas-phase processes in models, and also on how to present the necessary reaction data in databases; however, no such consensus has yet been reached for grain-surface processes. A team of
∼
25
experts covering observational, laboratory and theoretical (astro)chemistry met in summer of 2014 at the Lorentz Center in Leiden with the aim to provide solutions for this problem and to review the current state-of-the-art of grain surface models, both in terms of technical implementation into models as well as the most up-to-date information available from experiments and chemical computations. This review builds on the results of this workshop and gives an outlook for future directions.
Understanding the chemical evolution of young (high-mass) star-forming regions is a central topic in star formation research. Chemistry is employed as a unique tool 1) to investigate the underlying ...physical processes and 2) to characterize the evolution of the chemical composition. With these aims in mind, we observed a sample of 59 high-mass star-forming regions at different evolutionary stages varying from the early starless phase of infrared dark clouds to high-mass protostellar objects to hot molecular cores and, finally, ultra-compact Hii regions at 1 mm and 3 mm with the IRAM 30 m telescope. We determined their large-scale chemical abundances and found that the chemical composition evolves along with the evolutionary stages. On average, the molecular abundances increase with time. We modeled the chemical evolution, using a 1D physical model where density and temperature vary from stage to stage coupled with an advanced gas-grain chemical model and derived the best-fit χ2 values of all relevant parameters. A satisfying overall agreement between observed and modeled column densities for most of the molecules was obtained. With the best-fit model we also derived a chemical age for each stage, which gives the timescales for the transformation between two consecutive stages. The best-fit chemical ages are ~10 000 years for the IRDC stage, ~60 000 years for the HMPO stage, ~40 000 years for the HMC stage, and ~10 000 years for the UCHii stage. Thus, the total chemical timescale for the entire evolutionary sequence of the high-mass star formation process is on the order of 105 years, which is consistent with theoretical estimates. Furthermore, based on the approach of a multiple-line survey of unresolved data, we were able to constrain an intuitive and reasonable physical and chemical model. The results of this study can be used as chemical templates for the different evolutionary stages in high-mass star formation.
We present ALMA Cycle 2 observations at 0 5 resolution of TW Hya of CS emission. The radial profile of the integrated line emission displays oscillatory features outward of 1 5 ( au). A dip-like ...feature at 1 6 is coincident in location, depth, and width with features observed in dust scattered light at near-infrared wavelengths. Using a thermochemical model indicative of TW Hya, gas-grain chemical modeling, and non-LTE radiative transfer, we demonstrate that such a feature can be reproduced with a surface density depression, consistent with the modeling performed for scattered-light observations of TW Hya. We further demonstrate that a gap in the dust distribution and dust opacity only cannot reproduce the observed CS feature. The outer enhancement at 3 1 is identified as a region of intensified desorption due to enhanced penetration of the interstellar far-UV radiation at the exponential edge of the disk surface density, which intensifies the photochemical processing of gas and ices.
Context. Hydrodynamic, non-magnetic instabilities can provide turbulent stress in the regions of protoplanetary discs, where the magneto-rotational instability can not develop. The induced motions ...influence the grain growth, from which formation of planetesimals begins. Thermal relaxation of the gas constrains origins of the identified hydrodynamic sources of turbulence in discs. Aims. We aim to estimate the radiative relaxation timescale of temperature perturbations in protoplanetary discs. We study the dependence of the thermal relaxation on the perturbation wavelength, the location within the disc, the disc mass, and the dust-to-gas mass ratio. We then apply thermal relaxation criteria to localise modes of the convective overstability, the vertical shear instability, and the zombie vortex instability. Methods. For a given temperature perturbation, we estimated two timescales: the radiative diffusion timescale tthick and the optically thin emission timescale tthin. The longest of these timescales governs the relaxation: trelax = max (tthick, tthin). We additionally accounted for the collisional coupling to the emitting species. Our calculations employed the latest tabulated dust and gas mean opacities. Results. The relaxation criterion defines the bulk of a typical T Tauri disc as unstable to the development of linear hydrodynamic instabilities. The midplane is unstable to the convective overstability from at most 2au and up to 40au, as well as beyond 140au. The vertical shear instability can develop between 15au and 180au. The successive generation of (zombie) vortices from a seeded noise can work within the inner 0.8au. Conclusions. A map of relaxation timescale constrains the origins of the identified hydrodynamic turbulence-driving mechanisms in protoplanetary discs. Dynamic disc modelling with the evolution of dust and gas opacities is required to clearly localise the hydrodynamic turbulence, and especially its non-linear phase.
Context. Determining the gas density and temperature structures of protoplanetary disks is a fundamental task in order to constrain planet formation theories. This is a challenging procedure and most ...determinations are based on model-dependent assumptions. Aims. We attempt a direct determination of the radial and vertical temperature structure of the Flying Saucer disk, thanks to its favorable inclination of 90 degrees. Methods. We present a method based on the tomographic study of an edge-on disk. Using ALMA, we observe at 0.5″ resolution the Flying Saucer in CO J = 2–1 and CS J = 5–4. This edge-on disk appears in silhouette against the CO J = 2–1 emission from background molecular clouds in ρ Oph. The combination of velocity gradients due to the Keplerian rotation of the disk and intensity variations in the CO background as a function of velocity provide a direct measure of the gas temperature as a function of radius and height above the disk mid-plane. Results. The overall thermal structure is consistent with model predictions, with a cold (<12−15 K) CO-depleted mid-plane and a warmer disk atmosphere. However, we find evidence for CO gas along the mid-plane beyond a radius of about 200 au, coincident with a change of grain properties. Such behavior is expected in the case of efficient rise of UV penetration re-heating the disk and thus allowing CO thermal desorption or favoring direct CO photo-desorption. CO is also detected at up to 3–4 scale heights, while CS is confined to around 1 scale height above the mid-plane. The limits of the method due to finite spatial and spectral resolutions are also discussed. Conclusions. This method appears to be a very promising way to determine the gas structure of planet-forming disks, provided that the molecular data have an angular resolution which is high enough, on the order of 0.3−0.1″ at the distance of the nearest star-forming regions.
We present the results of a search for dark matter weakly interacting massive particles (WIMPs) in the mass range below 20 GeV/c^{2} using a target of low-radioactivity argon with a 6786.0 kg d ...exposure. The data were obtained using the DarkSide-50 apparatus at Laboratori Nazionali del Gran Sasso. The analysis is based on the ionization signal, for which the DarkSide-50 time projection chamber is fully efficient at 0.1 keVee. The observed rate in the detector at 0.5 keVee is about 1.5 event/keVee/kg/d and is almost entirely accounted for by known background sources. We obtain a 90% C.L. exclusion limit above 1.8 GeV/c^{2} for the spin-independent cross section of dark matter WIMPs on nucleons, extending the exclusion region for dark matter below previous limits in the range 1.8-6 GeV/c^{2}.