We used new ALMA $^{13}$CO and C$^{18}$O(3-2) observations obtained at high angular resolution ($\sim$0.2") together with previous CO(3-2) and (6-5) ALMA data and continuum maps at 1.3 and 0.8 mm in order to determine the gas properties (temperature, density, and kinematics) in the cavity and to a lesser extent in the outer disk of GG Tau A, the prototype of a young triple T Tauri star that is surrounded by a massive and extended Keplerian outer disk. By deprojecting, we studied the radial and azimuthal gas distribution and its kinematics. We also applied a new method to improve the deconvolution of the CO data and in particular better quantify the emission from gas inside the cavity. We perform local and nonlocal thermodynamic equilibrium studies in order to determine the excitation conditions and relevant physical parameters inside the ring and in the central cavity. Residual emission after removing a smooth-disk model indicates unresolved structures at our angular resolution, probably in the form of irregular rings or spirals. The outer disk is cold, with a temperature $<20$ K beyond 250 au that drops quickly (r$^{-1}$). The kinematics of the gas inside the cavity reveals infall motions at about 10% of the Keplerian speed. We derive the amount of gas in the cavity, and find that the brightest clumps, which contain about 10% of this mass, have kinetic temperatures 40$-$80 K, CO column densities of a few 10$^{17}$ cm$^{-2}$, and H$_2$ densities around 10$^7$ cm$^{-3}$. Although the gas in the cavity is only a small fraction of the disk mass, the mass accretion rate throughout the cavity is comparable to or higher than the stellar accretion rate. It is accordingly sufficient to sustain the circumstellar disks on a long timescale.

The formation of asteroids, comets and planets occurs in the interior of protoplanetary disks during the early phase of star formation. Consequently, the chemical composition of the disk might shape the properties of the emerging planetary system. In this context, it is crucial to understand whether and what organic molecules are synthesized in the disk. In this Letter, we report the first detection of formic acid (HCOOH) towards the TW Hydrae protoplanetary disk. The observations of the trans-HCOOH 6$_{(1,6)-5(1,5)}$ transition were carried out at 129~GHz with ALMA. We measured a disk-averaged gas-phase t-HCOOH column density of $\sim$ (2-4)$\times$10$^{12}$~cm$^{-2}$, namely as large as that of methanol. HCOOH is the first organic molecules containing two oxygen atoms detected in a protoplanetary disk, a proof that organic chemistry is very active even though difficult to observe in these objects. Specifically, this simplest acid stands as the basis for synthesis of more complex carboxylic acids used by life on Earth.

Context. Several sulfur-bearing molecules are observed in the interstellar medium and in comets, in strong contrast to protoplanetary disks where only CS, H$_2$CS and SO have been detected so far. Aims. We combine observations and chemical models to constrain the sulfur abundances and their sensitivity to physical and chemical conditions in the DM Tau protoplanetary disk. Methods. We obtained $0.5^{"}$ ALMA observations of DM Tau in Bands 4 and 6 in lines of CS, SO, SO$_2$, OCS, CCS, H$_2$CS and H$_2$S, achieving a $\sim 5$ mJy sensitivity. Using the non-LTE radiative transfer code RADEX and the forward-modeling tool DiskFit, disk-averaged CS column densities and upper limits for the other species were derived. Results. Only CS was detected with a derived column density of $\sim 2-6 \times 10^{12}$ cm$^{-2}$. We report a first tentative detection of SO$_2$ in DM Tau. The upper limits range between $\sim 10^{11}$ and $10^{14}$ cm$^{-2}$ for the other S-bearing species. The best-fit chemical model matching these values requires a gas-phase C/O ratio of > 1 at $r \sim 50-100$ au. With chemical modeling we demonstrate that sulfur-bearing species could be robust tracers of the gas-phase C/O ratio, surface reaction rates, grain size and UV intensities. Conclusions. The lack of detections of a variety of sulfur-bearing molecules in DM Tau other than CS implies a dearth of reactive sulfur in the gas phase, either through efficient freeze-out or because most of the elemental sulfur is in other large species, as found in comets. The inferred high CS/SO and CS/SO$_2$ ratios require a non-solar C/O gas-phase ratio of > 1, consistent with the recent observations of hydrocarbon rings in DM Tau. The stronger depletion of oxygen-bearing S-species compared to CS is likely linked to the low observed abundances of gaseous water in DM Tau and points to a removal mechanism of oxygen from the gas.

We aim at constraining the diffusive and advective transport processes in Jupiter's stratosphere. The IR spectrum recorded by Cassini/CIRS during the Jupiter flyby contains the fingerprints of several atmospheric compounds and allows probing its atmospheric composition. C2H2 and C2H6, the main compounds produced by methane photochemistry, were retrieved as a function of latitude at certain pressure levels. CIRS observations suggest a different meridional distribution for these two species, difficult to reconcile with their photochemical histories, which are thought to be coupled to the methane photolysis. While the C2H2 abundance decreases with latitude, C2H6 becomes more abundant at high latitudes. A new 2D (latitude-altitude) seasonal photochemical model is developed to study whether the addition of stratospheric transport processes, such as meridional diffusion and advection, can explain the latitudinal behavior of C2H2 and C2H6. C2H2 observations are fairly well reproduced without meridional diffusion. Adding meridional diffusion to the model provides an improved agreement with the C2H6 observations by flattening its meridional distribution, but degrades the fit to the C2H2 distribution. Meridional diffusion alone cannot produce the observed increase with latitude of C2H6. When adding advective transport between 30 mbar and 0.01 mbar, with upwelling winds at the equator and downwelling winds at high latitudes, we can reproduce the C2H6 abundance increase with latitude. However, the fit to the C2H2 distribution is degraded. The strength of the advective winds needed to reproduce the C2H6 abundances is very sensitive to the value of the meridional eddy diffusion coefficient. The coupled fate of C2H2 and C2H6 suggests that an additional process is missing in the model. Ion-neutral chemistry was not accounted for in this work and might be a good candidate to solve this issue.

We investigate the impact of an external pressure on the structure of self-gravitating polytropes for axially symmetric ellipsoids and rings. The confinement of the fluid by photons is accounted for through a boundary condition on the enthalpy $H$. Equilibrium configurations are determined numerically from a generalised "Self-Consistent-Field"-method. The new algorithm incorporates an intra-loop re-scaling operator ${\cal R}(H)$, which is essential for both convergence and getting self-normalised solutions. The main control parameter is the external-to-core enthalpy ratio. In the case of uniform rotation rate and uniform surrounding pressure, we compute the mass, the volume, the rotation rate and the maximum enthalpy. This is repeated for a few polytropic indices $n$. For a given axis ratio, over-pressurization globally increases all output quantities, and this is more pronounced for large $n$. Density profiles are flatter than in the absence of an external pressure. When the control parameter asymptotically tends to unity, the fluid converges toward the incompressible solution, whatever the index, but becomes geometrically singular. Equilibrium sequences, obtained by varying the axis ratio, are built. States of critical rotation are greatly exceeded or even disappear. The same trends are observed with differential rotation. Finally, the typical response to a photon point source is presented. Strong irradiation favours sharp edges. Applications concern star forming regions and matter orbiting young stars and black holes.

"Hot super-Earths" (or "Mini-Neptunes") between 1 and 4 times Earth's size with period shorter than 100 days orbit 30-50\% of Sun-like type stars. Their orbital configuration -- measured as the period ratio distribution of adjacent planets in multi-planet systems -- is a strong constraint for formation models. Here we use N-body simulations with synthetic forces from an underlying evolving gaseous disk to model the formation and long-term dynamical evolution of super-Earth systems. While the gas disk is present, planetary embryos grow and migrate inward to form a resonant chain anchored at the inner edge of the disk. These resonant chains are far more compact than the observed super-Earth systems. Once the gas dissipates resonant chains may become dynamically unstable. They undergo a phase of giant impacts that spreads the systems out. Disk turbulence has no measurable effect on the outcome. Our simulations match observations if a small fraction of resonant chains remain stable, while most super-Earths undergo a late dynamical instability. Our statistical analysis restricts the contribution of stable systems to less than $25\%$. Our results also suggest that the large fraction of observed single planet systems does not necessarily imply any dichotomy in the architecture of planetary systems. Finally, we use the low abundance of resonances in Kepler data to argue that, in reality, the survival of resonant chains happens likely only in $\sim 5\%$ of the cases. This leads to a mystery: in our simulations only 50-60\% of resonant chains became unstable whereas at least 75\% (and probably 90-95\%) must be unstable to match observations.

Thermochemical models have been used in the past to constrain the deep oxygen abundance in the gas and ice giant planets from tropospheric CO spectroscopic measurements. Knowing the oxygen abundance of these planets is a key to better understand their formation. These models have widely used dry and/or moist adiabats to extrapolate temperatures from the measured values in the upper troposphere down to the level where the thermochemical equilibrium between H$_2$O and CO is established. The mean molecular mass gradient produced by the condensation of H$_2$O stabilizes the atmosphere against convection and results in a vertical thermal profile and H$_2$O distribution that departs significantly from previous estimates. We revisit O/H estimates using an atmospheric structure that accounts for the inhibition of the convection by condensation. We use a thermochemical network and the latest observations of CO in Uranus and Neptune to calculate the internal oxygen enrichment required to satisfy both these new estimates of the thermal profile and the observations. We also present the current limitations of such modeling.

Deuterated species are unique and powerful tools in astronomy since they can probe the physical conditions, chemistry, and ionization level of various astrophysical media. Recent observations of several deuterated species along with some of their spin isomeric forms have rekindled the interest for more accurate studies on deuterium fractionation. This paper presents the first publicly available chemical network of multiply deuterated species along with spin chemistry implemented on the latest state-of-the-art gas-grain chemical code `NAUTILUS'. D/H ratios for all deuterated species observed at different positions of TMC-1 are compared with the results of our model, which considers multiply deuterated species along with the spin chemistry of light hydrogen bearing species H2, H2+, H3+ and their isotopologues. We also show the differences in the modeled abundances of non-deuterated species after the inclusion of deuteration and spin chemistry in the model. Finally, we present a list of potentially observable deuterated species in TMC-1 awaiting detection.

We compute the structure of a self-gravitating torus with polytropic equation-of-state (EOS) rotating in an imposed centrifugal potential. The Poisson-solver is based on isotropic multigrid with optimal covering factor (fluid section-to-grid area ratio). We work at $2$nd-order in the grid resolution for both finite difference and quadrature schemes. For soft EOS (i.e. polytropic index $n \ge 1$), the underlying $2$nd-order is naturally recovered for Boundary Values (BVs) and any other integrated quantity sensitive to the mass density (mass, angular momentum, volume, Virial Parameter, etc.), i.e. errors vary with the number $N$ of nodes per direction as $\sim 1/N^2$. This is, however, not observed for purely geometrical quantities (surface area, meridional section area, volume), unless a subgrid approach is considered (i.e. boundary detection). Equilibrium sequences are also much better described, especially close to critical rotation. Yet another technical effort is required for hard EOS ($n < 1$), due to infinite mass density gradients at the fluid surface. We fix the problem by using kernel splitting. Finally, we propose an accelerated version of the SCF-algorithm based on a node-by-node pre-conditionning of the mass density at each step. The computing time is reduced by a factor $2$ typically, regardless of the polytropic index. There is a priori no obstacle to applying these results and techniques to ellipsoidal configurations and even to $3$D-configurations.

In an atmosphere, a cloud condensation region is characterized by a strong vertical gradient in the abundance of the related condensing species. On Earth, the ensuing gradient of mean molecular weight has relatively few dynamical consequences because N$_2$ is heavier than water vapor, so that only the release of latent heat significantly impacts convection. On the contrary, in an hydrogen dominated atmosphere (e.g. giant planets), all condensing species are significantly heavier than the background gas. This can stabilize the atmosphere against convection near a cloud deck if the enrichment in the given species exceeds a critical threshold. This raises two questions. What is transporting energy in such a stabilized layer, and how affected can the thermal profile of giant planets be? To answer these questions, we first carry out a linear analysis of the convective and double-diffusive instabilities in a condensable medium showing that an efficient condensation can suppress double-diffusive convection. This suggests that a stable radiative layer can form near a cloud condensation level, leading to an increase in the temperature of the deep adiabat. Then, we investigate the impact of the condensation of the most abundant species---water---with a steady-state atmosphere model. Compared to standard models, the temperature increase can reach several hundred degrees at the quenching depth of key chemical tracers. Overall, this effect could have many implications for our understanding of the dynamical and chemical state of the atmosphere of giant planets, for their future observations (with Juno for example), and for their internal evolution.

Protoplanetary disks are the target of many chemical studies (both observational and theoretical) as they contain the building material for planets. Their large vertical and radial gradients in density and temperature make them challenging objects for chemical models. In the outer part of these disks, the large densities and low temperatures provide a particular environment where the binding of species onto the dust grains can be very efficient and can affect the gas-phase chemical composition. We attempt to quantify to what extent the vertical abundance profiles and the integrated column densities of molecules predicted by a detailed gas-grain code are affected by the treatment of the molecular hydrogen physisorption at the surface of the grains. We performed three different models using the Nautilus gas-grain code. One model uses a H2 binding energy on the surface of water (440 K) and produces strong sticking of H2. Another model uses a small binding energy of 23 K (as if there were already a monolayer of H2), and the sticking of H$_2$ is almost negligible. Finally, the remaining model is an intermediate solution known as the encounter desorption mechanism. We show that the efficiency of molecular hydrogen binding (and thus its abundance at the surface of the grains) can have a quantitative effect on the predicted column densities in the gas phase of major species such as CO, CS, CN, and HCN.

We present an extended version of the 2-phase gas-grain code NAUTILUS to the 3-phase modelling of gas and grain chemistry of cold cores. In this model, both the mantle and the surface are considered as chemically active. We also take into account the competition among reaction, diffusion and evaporation. The model predictions are confronted to ice observations in the envelope of low-mass and massive young stellar objects as well as toward background stars. Modelled gas-phase abundances are compared to species observed toward TMC-1 (CP) and L134N dark clouds. We find that our model successfully reproduces the observed ice species. It is found that the reaction-diffusion competition strongly enhances reactions with barriers and more specifically reactions with H2, which is abundant on grains. This finding highlights the importance to have a good approach to determine the abundance of H2 on grains. Consequently, it is found that the major N-bearing species on grains go from NH3 to N2 and HCN when the reaction-diffusion competition is accounted. In the gas-phase and before few 10^5 yrs, we find that the 3-phase model does not have a strong impact on the observed species compared to the 2-phase model. After this time, the computed abundances dramatically decrease due to the strong accretion on dust, which is not counterbalanced by the desorption less efficient than in the 2-phase model. This strongly constrains the chemical-age of cold cores to be of the order of few 10^5 yrs.

The first hydrostatic core, also called the first Larson core, is one of the first steps in low-mass star formation, as predicted by theory. With recent and future high performance telescopes, details of these first phases become accessible, and observations may confirm theory and even bring new challenges for theoreticians. In this context, we study from a theoretical point of view the chemical and physical evolution of the collapse of prestellar cores until the formation of the first Larson core, in order to better characterize this early phase in the star formation process. We couple a state-of-the-art hydrodynamical model with full gas-grain chemistry, using different assumptions on the magnetic field strength and orientation. We extract the different components of each collapsing core (i.e., the central core, the outflow, the disk, the pseudodisk, and the envelope) to highlight their specific physical and chemical characteristics. Each component often presents a specific physical history, as well as a specific chemical evolution. From some species, the components can clearly be differentiated. The different core models can also be chemically differentiated. Our simulation suggests some chemical species as tracers of the different components of a collapsing prestellar dense core, and as tracers of the magnetic field characteristics of the core. From this result, we pinpoint promising key chemical species to be observed.

Saturn's axial tilt produces seasons in a similar way as on Earth. Both the stratospheric temperature and composition are affected by this latitudinally varying insolation along the seasons. The thermal structure is controlled and regulated by the amount of hydrocarbons in the stratosphere, which act as absorbers and coolants from the UV to the far-IR spectral range, and this structure influences the amount of hydrocarbons. We study here the feedback between the chemical composition and the thermal structure by coupling a latitudinal and seasonal photochemical model with a radiative seasonal model. Our results show that the seasonal temperature peak in the higher stratosphere, associated with the seasonal increase of insolation, is shifted earlier than the maximum insolation peak. This shift is increased with increasing latitudes and is caused by the low amount of stratospheric coolants in the spring season. At 80$^{\circ}$ in both hemispheres, the temperature peak at 1d-2mbar is seen to occur half a season earlier than was previously predicted by radiative seasonal models that assumed spatially and temporally uniform distribution of coolants. This shift progressively decreases with increasing pressure, up to around the 0.5mbar pressure level where it vanishes. However, the thermal field has a small feedback on the abundance distributions. This feedback modifies the predicted equator-to-pole temperature gradient. The meridional gradients of temperature at the mbar pressure levels are better reproduced when this feedback is accounted for. At lower pressure levels, the thermal structure seems to depart from pure radiative seasonal equilibrium as previously suggested by Guerlet et al. (2014). Although the agreement with the absolute value of the stratospheric temperature observed by Cassini is moderate, it is a mandatory step toward a fully coupled GCM-photochemical model.

A large proportion of observed planetary systems contain several planets in a compact orbital configuration, and often harbor at least one close-in object. These systems are then most likely tidally evolving. We investigate how the effects of planet-planet interactions influence the tidal evolution of planets. We introduce for that purpose a new open-source addition to the Mercury N-body code, Mercury-T, which takes into account tides, general relativity and the effect of rotation-induced flattening in order to simulate the dynamical and tidal evolution of multi-planet systems. It uses a standard equilibrium tidal model, the constant time lag model. Besides, the evolution of the radius of several host bodies has been implemented (brown dwarfs, M-dwarfs of mass $0.1~M_\odot$, Sun-like stars, Jupiter). We validate the new code by comparing its output for one-planet systems to the secular equations results. We find that this code does respect the conservation of total angular momentum. We applied this new tool to the planetary system Kepler-62. We find that tides influence the stability of the system in some cases. We also show that while the four inner planets of the systems are likely to have slow rotation rates and small obliquities, the fifth planet could have a fast rotation rate and a high obliquity. This means that the two habitable zone planets of this system, Kepler-62e ad f are likely to have very different climate features, and this of course would influence their potential at hosting surface liquid water.

Reproducing Uranus and Neptune remains a challenge for simulations of solar system formation. The ice giants' peculiar obliquities suggest that they both suffered giant collisions during their formation. Thus, there must have been an epoch of accretion dominated by collisions among large planetary embryos in the primordial outer solar system. We test this idea using N-body numerical simulations including the effects of a gaseous protoplanetary disk. One strong constraint is that the masses of the ice giants are very similar -- the Neptune/Uranus mass ratio is $\sim1.18$. We show that similar-size ice giants do indeed form by collisions between planetary embryos beyond Saturn. The fraction of successful simulations varies depending on the initial number of planetary embryos in the system, their individual and total masses. Similar-sized ice giants are consistently reproduced in simulations starting with 5-10 planetary embryos with initial masses of $\sim$3-6 ${\rm M_\oplus}$. We conclude that accretion from a population of planetary embryos is a plausible scenario for the origin of Uranus and Neptune.

The gas mass of protoplanetary disks, and the gas-to-dust ratio, are two key elements driving the evolution of these disks and the formation of planetary system. We explore here to what extent CO (or its isotopologues) can be used as a tracer of gas mass. We use a detailed gas-grain chemical model and study the evolution of the disk composition, starting from a dense pre-stellar core composition. We explore a range of disk temperature profiles, cosmic rays ionization rates, and disk ages for a disk model representative of T Tauri stars. At the high densities that prevail in disks, we find that, due to fast reactions on grain surfaces, CO can be converted to less volatile forms (principally s-CO$_2$, and to a lesser extent s-CH$_4$) instead of being evaporated over a wide range of temperature. The canonical gas-phase abundance of 10$^{-4}$ is only reached above about 30-35 K. The dominant Carbon bearing entity depends on the temperature structure and age of the disk. The chemical evolution of CO is also sensitive to the cosmic rays ionization rate. Larger gas phase CO abundances are found in younger disks. Initial conditions, such as parent cloud age and density, have a limited impact. This study reveals that CO gas-phase abundance is heavily dependent on grain surface processes, which remain very incompletely understood so far. The strong dependence on dust temperature profile makes CO a poor tracer of the gas-phase content of disks.

Saturn's axial tilt of 26.7{\deg} produces seasons in a similar way as on Earth. Both the stratospheric temperature and composition are affected by this latitudinally varying insolation along Saturn's orbital path. A new time dependent 2D photochemical model is presented to study the seasonal evolution of Saturn's stratospheric composition. This study focuses on the impact of the seasonally variable thermal field on the main stratospheric C2 hydrocarbon chemistry (C2H2 and C2H6) using a realistic radiative climate model. Meridional mixing and advective processes are implemented in the model but turned off in the present study for the sake of simplicity. The results are compared to a simple study case where a latitudinally and temporally steady thermal field is assumed. Our simulations suggest that, when the seasonally variable thermal field is accounted for, the downward diffusion of the seasonally produced hydrocarbons is faster due to the seasonal compression of the atmospheric column during winter. This effect increases with increasing latitudes which experience the most important thermal changes in the course of the seasons. The seasonal variability of C2H2 and C2H6 therefore persists at higher-pressure levels with a seasonally-variable thermal field. Cassini limb-observations of C2H2 and C2H6 (Guerlet et al., 2009) are reasonably well-reproduced from the equator to 40{\deg} in both hemispheres

In this Letter we model the chemistry of DCO$^{+}$ in protoplanetary disks. We find that the overall distribution of the DCO$^{+}$ abundance is qualitatively similar to that of CO but is dominated by thin layer located at the inner disk surface. To understand its distribution, we investigate the different key gas-phase deuteration pathways that can lead to the formation of DCO$^{+}$. Our analysis shows that the recent update in the exothermicity of the reaction involving CH$_2$D$^{+}$ as a parent molecule of DCO$^{+}$ favors deuterium fractionation in warmer conditions. As a result the formation of DCO$^{+}$ is enhanced in the inner warm surface layers of the disk where X-ray ionization occurs. Our analysis points out that DCO$^{+}$ is not a reliable tracer of the CO snow line as previously suggested. We thus predict that DCO$^{+}$ is a tracer of active deuterium and in particular X-ray ionization of the inner disk.

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

Planets of 1-4 times Earth's size on orbits shorter than 100 days exist around 30-50% of all Sun-like stars. In fact, the Solar System is particularly outstanding in its lack of "hot super-Earths" (or "mini-Neptunes"). These planets -- or their building blocks -- may have formed on wider orbits and migrated inward due to interactions with the gaseous protoplanetary disk. Here, we use a suite of dynamical simulations to show that gas giant planets act as barriers to the inward migration of super-Earths initially placed on more distant orbits. Jupiter's early formation may have prevented Uranus and Neptune (and perhaps Saturn's core) from becoming hot super-Earths. Our model predicts that the populations of hot super-Earth systems and Jupiter-like planets should be anti-correlated: gas giants (especially if they form early) should be rare in systems with many hot super-Earths. Testing this prediction will constitute a crucial assessment of the validity of the migration hypothesis for the origin of close-in super-Earths.

Recent observations have revealed the existence of Complex Organic Molecules (COMs) in cold dense cores and prestellar cores. The presence of these molecules in such cold conditions is not well understood and remains a matter of debate since the previously proposed "warm- up" scenario cannot explain these observations. In this article, we study the effect of Eley- Rideal and complex induced reaction mechanisms of gas-phase carbon atoms with the main ice components of dust grains on the formation of COMs in cold and dense regions. Based on recent experiments we use a low value for the chemical desorption efficiency (which was previously invoked to explain the observed COM abundances). We show that our introduced mechanisms are efficient enough to produce a large amount of complex organic molecules in the gas-phase at temperatures as low as 10K.

Planetary embryos embedded in gaseous protoplanetary disks undergo Type I orbital migration. Migration can be inward or outward depending on the local disk properties but, in general, only planets more massive than several $M_\oplus$ can migrate outward. Here we propose that an embryo's migration history determines whether it becomes a hot super-Earth or the core of a giant planet. Systems of hot super-Earths (or mini-Neptunes) form when embryos migrate inward and pile up at the inner edge of the disk. Giant planet cores form when inward-migrating embryos become massive enough to switch direction and migrate outward. We present simulations of this process using a modified N-body code, starting from a swarm of planetary embryos. Systems of hot super-Earths form in resonant chains with the innermost planet at or interior to the disk inner edge. Resonant chains are disrupted by late dynamical instabilities triggered by the dispersal of the gaseous disk. Giant planet cores migrate outward toward zero-torque zones, which move inward and eventually disappear as the disk disperses. Giant planet cores migrate inward with these zones and are stranded at ~1-5 AU. Our model reproduces several properties of the observed extra-solar planet populations. The frequency of giant planet cores increases strongly when the mass in solids is increased, consistent with the observed giant exoplanet - stellar metallicity correlation. The frequency of hot super-Earths is not a function of stellar metallicity, also in agreement with observations. Our simulations can reproduce the broad characteristics of the observed super-Earth population.

The Kepler-186 system consists of five planets orbiting an early-M dwarf. The planets have physical radii of 1.0-1.50 R$_\oplus$ and orbital periods of 4 to 130 days. The $1.1~$R$_\oplus$ Kepler-186f with a period of 130 days is of particular interest. Its insolation of roughly $0.32~S_\odot$places it within the liquid water habitable zone. We present a multi-faceted study of the Kepler-186 system. First, we show that the distribution of planet masses can be roughly reproduced if the planets accreted from a high-surface density disk presumably sculpted by an earlier phase of migration. However, our simulations predict the existence of 1-2 undetected planets between planets e and f. Next, we present a dynamical analysis of the system including the effect of tides. The timescale for tidal evolution is short enough that the four inner planets must have small obliquities and near-synchronous rotation rates. Tidal evolution of Kepler-186f is slow enough that its current spin state depends on a combination of its dissipation rate and the stellar age. Finally, we study the habitability of Kepler-186f with a 1-D climate model. The planet's surface temperature can be raised above 273 K with 0.5-5 bars of CO$_2$, depending on the amount of N$_2$ present. Kepler-186f represents a case study of an Earth-sized planet in the cooler regions of the habitable zone of a cool star.

Grain-surface reactions play an essential role in interstellar chemistry, since dust grain catalyses reactions at its surface allowing for the formation of molecules. We used a chemical model in which both gas-phase and grain-surface reactions occur and studied particularly the diffusion mechanisms on the surface of the grains. Surface reactions can occur via thermal hopping when species cross over a potential barrier or via quantum tunneling when species cross through this barrier. We show that the thermal diffusion (hopping) can be much more efficient after a cosmic ray particle collides with a dust grain, heating it to a peak temperature of 70K. We present here the results of numerical simulations after including the quantum tunneling mechanism for species H, H2 and O and considering the effect of cosmic ray particle collision on the surface reactions. As a consequence, the gas-phase and grain-surface abundances are affected and we show that more complex molecules can be formed in molecular clouds.

We have developed a pseudo two-dimensional model of a planetary atmosphere, which takes into account thermochemical kinetics, photochemistry, vertical mixing, and horizontal transport, the latter being modeled as a uniform zonal wind. We have applied the model to the atmospheres of the hot Jupiters HD 209458b and HD 189733b. The adopted eddy diffusion coefficients are calculated by following the behaviour of passive tracers in three-dimensional general circulation models, which results in eddy values significantly below previous estimates. We find that the distribution of molecules with altitude and longitude in the atmospheres of these two hot Jupiters is complex because of the interplay of the various physical and chemical processes at work. Much of the distribution of molecules is driven by the strong zonal wind and the limited extent of vertical transport, resulting in an important homogenisation of the chemical composition with longitude. In general, molecular abundances are quenched horizontally to values typical of the hottest dayside regions, and thus the composition in the cooler nightside regions is highly contaminated by that of warmer dayside regions. As a consequence, the abundance of methane remains low, even below the predictions of previous one-dimensional models, which is likely to be in conflict with the high CH4 content inferred from observations of the dayside of HD 209458b. Another consequence of the important longitudinal homogenisation of the abundances is that the variability of the chemical composition has little effect on the way the emission spectrum is modified with phase and on the changes in the transmission spectrum from the transit ingress to the egress, these variations in the spectra being mainly due to changes in the temperature rather than in the composition between the different sides of the planet.

Protoplanetary disks composed of dust and gas are ubiquitous around young stars and are commonly recognized as nurseries of planetary systems. Their lifetime, appearance, and structure are determined by an interplay between stellar radiation, gravity, thermal pressure, magnetic field, gas viscosity, turbulence, and rotation. Molecules and dust serve as major heating and cooling agents in disks. Dust grains dominate the disk opacities, reprocess most of the stellar radiation, and shield molecules from ionizing UV/X-ray photons. Disks also dynamically evolve by building up planetary systems which drastically change their gas and dust density structures. Over the past decade significant progress has been achieved in our understanding of disk chemical composition thanks to the upgrade or advent of new millimeter/Infrared facilities (SMA, PdBI, CARMA, Herschel, e-VLA, ALMA). Some major breakthroughs in our comprehension of the disk physics and chemistry have been done since PPV. This review will present and discuss the impact of such improvements on our understanding of the disk physical structure and chemical composition.

The local character of self-gravity along with the number of spatial dimensions are critical issues when computing the potential and forces inside massive systems like stars and disks. This appears from the discretisation scale where each cell of the numerical grid is a self-interacting body in itself. There is apparently no closed-form expression yet giving the potential of a three-dimensional homogeneous cylindrical or spherical cell, in contrast with the Cartesian case. By using Green's theorem, we show that the potential integral for such polar-type 3D sectors -- initially, a volume integral with singular kernel -- can be converted into a regular line-integral running over the lateral contour, thereby generalising a formula already known under axial symmetry. It therefore is a step towards the obtention of another potential/density pair. The new kernel is a finite function of the cell's shape (with the simplest form in cylindrical geometry), and mixes incomplete elliptic integrals, inverse trigonometric and hyperbolic functions. The contour integral is easy to compute; it is valid in the whole physical space, exterior and interior to the sector itself and works in fact for a wide variety of shapes of astrophysical interest (e.g. sectors of tori or flared discs). This result is suited to easily providing reference solutions, and to reconstructing potential and forces in inhomogeneous systems by superposition. The contour integrals for the 3 components of the acceleration vector are explicitely given.

As hydrostatic equilibrium of gaseous discs is partly governed by the gravity field, we have estimated the component caused by a vertically homogeneous disc, with a special attention for the outer regions where self-gravity classically appears. The accuracy of the integral formula is better than 1%, whatever the disc thickness, radial extension and radial density profile. At order zero, the field is even algebraic for thin discs and writes $- 4 \pi G \Sigma(R) f_{edge} (R)$ at disc surface, thereby correcting Paczynski's formula by a multiplying factor $f_{edge} \gtrsim 1/2$, which depends on the relative distance to the edges and the local disc thickness. For very centrally condensed discs however, this local contribution can be surpassed by action of mass stored in the inner regions, possibly resulting in $f_{edge} \gg 1$. A criterion setting the limit between these two regimes is derived. These result are robust in the sense that the details of vertical stratification are not critical. We briefly discuss how hydrostatic equilibrium is impacted. In particular, the disc flaring should not reverse in the self-gravitating region, which contradicts what is usually obtained from Paczynski's formula. This suggests that i) these outer regions are probably not fully shadowed by the inner ones (important when illuminated by a central star), and ii) the flared shape of discs does not firmly prove the absence or weakness of self-gravity.

Context. Carbon monoxide (CO) has been detected in all Giant Planets and its origin is both internal and external in Jupiter and Neptune. Despite its first detection in Uranus a decade ago, the magnitude of its internal and external sources remains unconstrained. Aims. We targeted CO lines in Uranus in the submillimeter range to constrain its origin. Methods. We recorded disk-averaged spectra of Uranus with a very high spectral resolution at the frequencies of CO rotational lines in the submillimeter range in 2011-2012. We used empirical and diffusion models of the atmosphere of Uranus to constrain the origin of CO. We also used a thermochemical model of its troposphere to derive an upper limit on the O/H ratio in the deep atmosphere of Uranus. Results. We have detected the CO(8-7) rotational line for the first time, with Herschel-HIFI. Both empirical and diffusion model results show that CO has an external origin. An empirical profile in which CO is constant above the 100mbar level with a mole fraction of 7.1-9.0E-9 (depending on the adopted stratospheric thermal structure) reproduces the data. Sporadic and steady source models cannot be differentiated with our data. Taking Teanby & Irwin's internal source model upper limit of a mole fraction of 2.1E-9 [Teanby & Irwin 2013. ApJ, 775, L49], the deep O/H ratio of Uranus is lower than 500 times solar according to our thermochemical computations. Conclusions. Our work shows that the average mole fraction of CO decreases from the stratosphere to the troposphere and thus strongly advocates for an external source of CO in Uranus. Photochemical modeling of oxygen species in the atmosphere of Uranus and more sensitive observations are needed to reveal the nature of the external source.

We investigate water and deuterated water chemistry in turbulent protoplanetary disks. Chemical rate equations are solved with the diffusion term, mimicking turbulent mixing in vertical direction. Water near the midplane is transported to the disk atmosphere by turbulence and destroyed by photoreactions to produce atomic oxygen, while the atomic oxygen is transported to the midplane and reforms water and/or other molecules. We find that this cycle significantly decreases column densities of water ice at r < 30 AU, where dust temperatures are too high to reform water ice effectively. The radial extent of such region depends on the desorption energy of atomic hydrogen. Our model indicates that water ice could be deficient even outside the sublimation radius. Outside this radius, the cycle decreases the D/H ratio of water ice from 2x10^-2, which is set by the collapsing core model, to 10^-4-10^-2 in 10^6 yr, without significantly decreasing the water ice column density. The resultant D/H ratios depend on the strength of mixing and the radial distance from the central star. Our finding suggests that the D/H ratio of cometary water (10^-3-10^-4) could be established (i.e. cometary water could be formed) in the solar nebula, even if the D/H ratio of water ice delivered to the disk was very high (10^-2).

An outstanding question of astrobiology is the link between the chemical composition of planets, comets, and other solar system bodies and the molecules formed in the interstellar medium. Understanding the chemical and physical evolution of the matter leading to the formation of protoplanetary disks is an important step for this. We provide some new clues to this long-standing problem using three-dimensional chemical simulations of the early phases of disk formation: we interfaced the full gas-grain chemical model Nautilus with the radiation-magnetohydrodynamic model RAMSES, for different configurations and intensities of magnetic field. Our results show that the chemical content (gas and ices) is globally conserved during the collapsing process, from the parent molecular cloud to the young disk surrounding the first Larson core. A qualitative comparison with cometary composition suggests that comets are constituted of different phases, some molecules being direct tracers of interstellar chemistry, while others, including complex molecules, seem to have been formed in disks, where higher densities and temperatures allow for an active grain surface chemistry. The latter phase, and its connection with the formation of the first Larson core, remains to be modeled.

Given their tendency to be incorporated into the core during differentiation, the highly-siderophile elements (HSEs) in Earth's mantle are thought to have been accreted as a `late veneer' after the end of the giant impact phase. Bottke et al (2010) proposed that the large Earth-to-Moon HSE abundance ratio can be explained if the late veneer was characterized by large (D = 1000-4000km) impactors. Here we simulate the evolution of the terrestrial planets during a stochastic late veneer phase from the end of accretion until the start of the late heavy bombardment ~500 Myr later. We show that a late veneer population of 0.05 Earth masses dominated by large (D > 1000km) bodies naturally delivers a ~0.01 Earth mass veneer to Earth, consistent with constraints. The eccentricities and inclinations of the terrestrial planets are excited by close encounters with the largest late veneer bodies. We find the best agreement with their post-veneer orbits if either a) the terrestrial planets' pre-veneer angular momentum deficit AMD_0 was less than half of the current one AMD_now, or b) AMD_0 <= AMD_now and the veneer was limited to smaller (D_max <= 2000km) bodies. Veneer impacts on Venus, Earth and Mars were mostly accretionary but on Mercury and the Moon they were mostly erosive. In ~20% of simulations an energetic impact occurred that could have removed >25% of Mercury's mass, thereby increasing its iron mass fraction. We show that, due to the erosive nature of larger impacts, the Moon cannot accrete any material from objects larger than 500-1000km. The large Earth-to-Moon HSE abundance ratio is explained if the late veneer included large impactors (D >= 500-1000km) regardless of their size distribution, as long as most of Earth's veneer came from large bodies. The spin angular momenta imparted by stochastic late veneer impacts was far in excess of the current ones for Mercury and Venus.

Short period (<50 days) low-mass (<10Mearth) exoplanets are abundant and the few of them whose radius and mass have been measured already reveal a diversity in composition. Some of these exoplanets are found on eccentric orbits and are subjected to strong tides affecting their rotation and resulting in significant tidal heating. Within this population, some planets are likely to be depleted in volatiles and have no atmosphere. We model the thermal emission of these "Super Mercuries" to study the signatures of rotation and tidal dissipation on their infrared light curve. We compute the time-dependent temperature map at the surface and in the subsurface of the planet and the resulting disk-integrated emission spectrum received by a distant observer for any observation geometry. We calculate the illumination of the planetary surface for any Keplerian orbit and rotation. We include the internal tidal heat flow, vertical heat diffusion in the subsurface and generate synthetic light curves. We show that the different rotation periods predicted by tidal models (spin-orbit resonances, pseudo-synchronization) produce different photometric signatures, which are observable provided that the thermal inertia of the surface is high, like that of solid or melted rocks (but not regolith). Tidal dissipation can also directly affect the light curves and make the inference of the rotation more difficult or easier depending on the existence of hot spots on the surface. Infrared light curve measurement with the James Webb Space Telescope and EChO can be used to infer exoplanets' rotation periods and dissipation rates and thus to test tidal models. This data will also constrain the nature of the (sub)surface by constraining the thermal inertia.

We use the super Earth 55 Cnc e as a case study to address an observable effect of tidal heating. We investigate whether planet-planet interactions can force the eccentricity of this planet to a level affecting the eclipse depth observed with Spitzer. Using the constant time lag tidal model, we first calculate the observed planet flux as a function of albedo and eccentricity, for different tidal dissipation constants and for two extreme cases: a planet with no heat redistribution and a planet with full heat redistribution. We derive the values of albedo and eccentricity that match the observed transit depth. We then perform N-body simulations of the planetary system including tides and General Relativity to follow the evolution of the eccentricity of planet e. We compare the range of eccentricities given by the simulations with the eccentricities required to alter the eclipse depth. We find that the eccentricity of planet e can be large enough to contribute at a measurable level to the thermal emission measured with Spitzer. This affects the constraints on the albedo of the planet, which can be as high as 0.9 (instead of 0.55 when ignoring tidal heating). We also derive a maximum value for the eccentricity of planet e directly from the eclipse depth: e<0.015 assuming Earth's dissipation constant. Transiting exoplanets in multiple planet systems are exceptional targets for testing tidal models because their tidal luminosity may be observable. Future multi-wavelengths observations of eclipse depth and phase curves (for instance with EChO and JWST) should allow us to better resolve the temperature map of these planets and break the degeneracy between albedo and tidal heating that remains for single band observations. In addition, an accurate determination of the eccentricity will make it possible to constrain the dissipation rate of the planet and to probe its internal structure.

Turbulence is thought to be a key driver of the evolution of protoplanetary disks, regulating the mass accretion process, the transport of angular momentum, and the growth of dust particles. We intend to determine the magnitude of the turbulent motions in the outer parts (> 100 AU) of the disk surrounding DM Tau. Turbulent motions can be constrained by measuring the nonthermal broadening of line emission from heavy molecules. We used the IRAM Plateau de Bure interferometer to study emission from the CS molecule in the disk of DM Tau. High spatial (1.4 x 1 ") and spectral resolution (0.126 km/s) CS J=3-2 images provide constraints on the molecule distribution and velocity structure of the disk. A low sensitivity CS J=5-4 image was used in conjunction to evaluate the excitation conditions. We analyzed the data in terms of two parametric disk models, and compared the results with detailed time-dependent chemical simulations. The CS data confirm the relatively low temperature suggested by observations of other simple molecules. The intrinsic linewidth derived from the CS J=3-2 data is much larger than expected from pure thermal broadening. The magnitude of the derived nonthermal component depends only weakly on assumptions about the location of the CS molecules with respect to the disk plane. Our results indicate turbulence with a Mach number around 0.4 - 0.5 in the molecular layer. Geometrical constraints suggest that this layer is located near one scale height, in reasonable agreement with chemical model predictions.

We investigate the effects of atmospheric circulation on the chemistry of the hot Jupiter HD 209458b. We use a simplified dynamical model and a robust chemical network, as opposed to previous studies which have used a three dimensional circulation model coupled to a simple chemical kinetics scheme. The temperature structure and distribution of the main atmospheric constituents are calculated in the limit of an atmosphere that rotates as a solid body with an equatorial rotation rate of 1 km/s. Such motion mimics a uniform zonal wind which resembles the equatorial superrotation structure found by three dimensional circulation models. The uneven heating of this tidally locked planet causes, even in the presence of such a strong zonal wind, large temperature contrasts between the dayside and nightside, of up to 800 K. This would result in important longitudinal variations of some molecular abundances if the atmosphere were at chemical equilibrium. The zonal wind, however, acts as a powerful disequilibrium process. We identify the existence of a pressure level of transition between two regimes, which may be located between 100 and 0.1 mbar depending on the molecule. Below this transition layer, chemical equilibrium holds, while above it, the zonal wind tends to homogenize the chemical composition of the atmosphere, bringing molecular abundances in the limb and nightside regions close to chemical equilibrium values characteristic of the dayside, i.e. producing an horizontal quenching effect in the abundances. Reasoning based on timescales arguments indicates that horizontal and vertical mixing are likely to compete in HD 209458b's atmosphere, producing a complex distribution where molecular abundances are quenched horizontally to dayside values and vertically to chemical equilibrium values characteristic of deep layers.

We investigate the molecular evolution and D/H abundance ratios that develop as star formation proceeds from a dense-cloud core to a protostellar core, by solving a gas-grain reaction network applied to a 1-D radiative hydrodynamic model with infalling fluid parcels. Spatial distributions of gas and ice-mantle species are calculated at the first-core stage, and at times after the birth of a protostar. Gas-phase methanol and methane are more abundant than CO at radii $r\lesssim 100$ AU in the first-core stage, but gradually decrease with time, while abundances of larger organic species increase. The warm-up phase, when complex organic molecules are efficiently formed, is longer-lived for those fluid parcels in-falling at later stages. The formation of unsaturated carbon chains (warm carbon-chain chemistry) is also more effective in later stages; C$^+$, which reacts with CH$_4$ to form carbon chains, increases in abundance as the envelope density decreases. The large organic molecules and carbon chains are strongly deuterated, mainly due to high D/H ratios in the parent molecules, determined in the cold phase. We also extend our model to simulate simply the chemistry in circumstellar disks, by suspending the 1-D infall of a fluid parcel at constant disk radii. The species CH$_3$OCH$_3$ and HCOOCH$_3$ increase in abundance in $10^4-10^5$ yr at the fixed warm temperature; both also have high D/H ratios.

We investigate molecular evolution from a molecular cloud core to a first hydrostatic core in three spatial dimensions. We perform a radiation hydrodynamic simulation in order to trace fluid parcels, in which molecular evolution is investigated, using a gas-phase and grain-surface chemical reaction network. We derive spatial distributions of molecular abundances and column densities in the core harboring the first core. We find that the total of gas and ice abundances of many species in a cold era (10 K) remain unaltered until the temperature reaches ~500 K. The gas abundances in the warm envelope and the outer layer of the first core (T < 500 K) are mainly determined via the sublimation of ice-mantle species. Above 500 K, the abundant molecules, such as H2CO, start to be destroyed, and simple molecules, such as CO, H2O and N2 are reformed. On the other hand, some molecules are effectively formed at high temperature; carbon-chains, such as C2H2 and cyanopolyynes, are formed at the temperature of >700 K. We also find that large organic molecules, such as CH3OH and HCOOCH3, are associated with the first core (r < 10 AU). Although the abundances of these molecules in the first core stage are comparable or less than in the protostellar stage (hot corino), reflecting the lower luminosity of the central object, their column densities in our model are comparable to the observed values toward the prototypical hot corino, IRAS 16293-2422. We propose that these large organic molecules can be good tracers of the first cores.

Low-mass objects embedded in isothermal protoplanetary discs are known to suffer rapid inward Type I migration. In non-isothermal discs, recent work has shown that a decreasing radial profile of the disc entropy can lead to a strong positive corotation torque which can significantly slow down or reverse Type I migration in laminar viscous disc models. It is not clear however how this picture changes in turbulent disc models. The aim of this study is to examine the impact of turbulence on the torque experienced by a low-mass planet embedded in a non-isothermal protoplanetary disc. We particularly focus on the role of turbulence on the corotation torque whose amplitude depends on the efficiency of diffusion processes in the planet's horseshoe region. We performed 2D numerical simulations using a grid-based hydrodynamical code in which turbulence is modelled as stochastic forcing. In order to provide estimations for the viscous and thermal diffusion coefficients as a function of the amplitude of turbulence, we first set up non-isothermal disc models for different values of the amplitude of the turbulent forcing. We then include a low-mass planet and determine the evolution of its running time-averaged torque. We show that in non-isothermal discs, the entropy-related corotation torque can indeed remain unsaturated in the presence of turbulence. For turbulence amplitudes that do not strongly affect the disc temperature profile, we find that the running time-averaged torque experienced by an embedded protoplanet is in fairly good agreement with laminar disc models with appropriate values for the thermal and viscous diffusion coefficients. In discs with turbulence driven by stochastic forcing, the corotation torque therefore behaves similarly as in laminar viscous discs and can be responsible for significantly slowing down or reversing Type I migration.

Our purpose is to release a chemical network, and the associated rate coefficients, developed for the temperature and pressure range relevant to hot Jupiters atmospheres. Using this network, we study the vertical atmospheric composition of the two hot Jupiters (HD209458b, HD189733b) with a model that includes photolyses and vertical mixing and we produce synthetic spectra. The chemical scheme is derived from applied combustion models that have been methodically validated over a range of temperatures and pressures typical of the atmospheric layers influencing the observations of hot Jupiters. We compare the predictions obtained from this scheme with equilibrium calculations, with different schemes available in the literature that contain N-bearing species and with previously published photochemical models. Compared to other chemical schemes that were not subjected to the same systematic validation, we find significant differences whenever non-equilibrium processes take place. The deviations from the equilibrium, and thus the sensitivity to the network, are more important for HD189733b, as we assume a cooler atmosphere than for HD209458b. We found that the abundances of NH3 and HCN can vary by two orders of magnitude depending on the network, demonstrating the importance of comprehensive experimental validation. A spectral feature of NH3 at 10.5$\mu$m is sensitive to these abundance variations and thus to the chemical scheme. Due to the influence of the kinetics, we recommend the use of a validated scheme to model the chemistry of exoplanet atmospheres. Our network is robust for temperatures within 300-2500K and pressures from 10mbar up to a few hundreds of bars, for species made of C,H,O,N. It is validated for species up to 2 carbon atoms and for the main nitrogen species.

Molecular line emission from protoplanetary disks is a powerful tool to constrain their physical and chemical structure. Nevertheless, only a few molecules have been detected in disks so far. We take advantage of the enhanced capabilities of the IRAM 30m telescope by using the new broad band correlator (FTS) to search for so far undetected molecules in the protoplanetary disks surrounding the TTauri stars DM Tau, GO Tau, LkCa 15 and the Herbig Ae star MWC 480. We report the first detection of HC3N at 5 sigma in the GO Tau and MWC 480 disks with the IRAM 30-m, and in the LkCa 15 disk (5 sigma), using the IRAM array, with derived column densities of the order of 10^{12}cm^{-2}. We also obtain stringent upper limits on CCS (N < 1.5 x 10^{12} cm^{-3}). We discuss the observational results by comparing them to column densities derived from existing chemical disk models (computed using the chemical code Nautilus) and based on previous nitrogen and sulfur-bearing molecule observations. The observed column densities of HC3N are typically two orders of magnitude lower than the existing predictions and appear to be lower in the presence of strong UV flux, suggesting that the molecular chemistry is sensitive to the UV penetration through the disk. The CCS upper limits reinforce our model with low elemental abundance of sulfur derived from other sulfur-bearing molecules (CS, H2S and SO).

By using various properties of the complete elliptic integrals, we have derived an alternative expression for the gravitational potential of axially symmetric bodies, which is free of singular kernel in contrast with the classical form. This is mainly a radial integral of the local surface density weighted by a regular "mean Green function" which depends explicitly on the body's vertical thickness. Rigorously, this result stands for a wide variety of configurations, as soon as the density structure is vertically homogeneous. Nevertheless, the sensitivity to vertical stratification | the Gaussian profile has been considered | appears weak provided that the surface density is conserved. For bodies with small aspect ratio (i.e. geometrically thin discs), a first-order Taylor expansion furnishes an excellent approximation for this mean Green function, the absolute error being of the fourth order in the aspect ratio. This formula is therefore well suited to studying the structure of self-gravitating discs and rings in the spirit of the "standard model of thin discs" where the vertical structure is often ignored, but it remains accurate for discs and tori of finite thickness. This approximation which perfectly saves the properties of Newton's law everywhere (in particular at large separations), is also very useful for dynamical studies where the body is just a source of gravity acting on external test particles.

We present a novel chemical database for gas-phase astrochemistry. Named the KInetic Database for Astrochemistry (KIDA), this database consists of gas-phase reactions with rate coefficients and uncertainties that will be vetted to the greatest extent possible. Submissions of measured and calculated rate coefficients are welcome, and will be studied by experts before inclusion into the database. Besides providing kinetic information for the interstellar medium, KIDA is planned to contain such data for planetary atmospheres and for circumstellar envelopes. Each year, a subset of the reactions in the database (kida.uva) will be provided as a network for the simulation of the chemistry of dense interstellar clouds with temperatures between 10 K and 300 K. We also provide a code, named Nahoon, to study the time-dependent gas-phase chemistry of 0D and 1D interstellar sources.

Context. The photometric signal we receive from a star hosting a planet is modulated by the variation of the planet signal with its orbital phase. Such phase variations are observed for transiting hot Jupiters with current instrumentation, and have also been measured for one transiting terrestrial planet (Kepler 10 b) and one non-transiting gas giant (Ups A b). Future telescopes (JWST and EChO) will have the capability to measure thermal phase curves of exoplanets including hot rocky planets in transiting and non-transiting configurations, and at different wavelengths. Short-period planets with a mass below 10 R_EARTH are indeed frequent and nearby targets (within 10 pc) are already known and more are to be found. Aims. To test the possibility to use multi-wavelengths infrared phase curves to constrain the radius, the albedo and the orbital inclination of a non-transiting planet with no atmosphere and on a 1:1 spin orbit resonance. Methods. We model the thermal emission of a synchronous rocky planet with no atmosphere and its apparent variation with the orbital phase for a given orbital inclination. We assume that the planet is detected by radial velocity so its orbital period and minimum mass are known. We simulate observed noisy phase curves and then apply a procedure to retrieve the radius and albedo of the planet and the inclination of the orbit. Results. Airless planets can be distinguished from planets having a dense atmosphere and their radius, albedo and inclination (and therefore true mass) can be retrieved from multiband observations with MIRI-JWST and EChO in the 5-15 {\mu}m range. [...] As inclination above 60{\deg} represents half of randomly oriented orbits, the growing population of short-period terrestrial-sized planets detected by radial velocity surveys and transits should offer several nearby promising targets for this method [...]. See full abstract in the paper.

We study the content in S-bearing molecules of protoplanetary disks around low-mass stars. We used the new IRAM 30-m receiver EMIR to perform simultaneous observations of the $1_{10}-1_{01}$ line of H$_2$S at 168.8 GHz and $2_{23}-1_{12}$ line of SO at 99.3 GHz. We compared the observational results with predictions coming from the astrochemical code NAUTILUS, which has been adapted to protoplanetary disks. The data were analyzed together with existing CS J=3-2 observations. We fail to detect the SO and H$_2$S lines, although CS is detected in LkCa15, DM\,Tau, and GO\,Tau but not in MWC\,480. However, our new upper limits are significantly better than previous ones and allow us to put some interesting constraints on the sulfur chemistry. Our best modeling of disks is obtained for a C/O ratio of 1.2, starting from initial cloud conditions of H density of $2\times 10^5$ cm$^{-3}$ and age of $10^6$ yr. The results agree with the CS data and are compatible with the SO upper limits, but fail to reproduce the H$_2$S upper limits. The predicted H$_2$S column densities are too high by at least one order of magnitude. H$_2$S may remain locked onto grain surfaces and react with other species, thereby preventing the desorption of H$_2$S.

We exploit our formula for the gravitational potential of finite size, power-law disks to derive a general expression linking the mass of the black hole in active galactic nuclei (AGN), the mass of the surrounding disk, its surface density profile (through the power index s), and the differential rotation law. We find that the global rotation curve v(R) of the disk in centrifugal balance does not obey a power law of the cylindrical radius R (except in the confusing case s = -2 that mimics a Keplerian motion), and discuss the local velocity index. This formula can help to understand how, from position-velocity diagrams, mass is shared between the disk and the black hole. To this purpose, we have checked the idea by generating a sample of synthetic data with different levels of Gaussian noise, added in radius. It turns out that, when observations are spread over a large radial domain and exhibit low dispersion (standard deviation less than 10% typically), the disk properties (mass and s-parameter) and black hole mass can be deduced from a non linear fit of kinematic data plotted on a (R, Rv 2)-diagram. For a deviation higher than 10%, masses are estimated fairly well from a linear regression (corresponding to the zeroth-order treatment of the formula), but the power index s is no longer accessible. We have applied the model to 7 AGN disks whose rotation has already been probed through water maser emission. For NGC3393 and UGC3789, the masses seem well constrained through the linear approach. For IC1481, the power-law exponent s can even be deduced. Because the model is scale-free, it applies to any kind of star/disk system. Extension to disks around young stars showing deviation from Keplerian motion is thus straightforward.

The one-dimensional, ordinary differential equation (ODE) by Hur\'e & Hersant (2007) that satisfies the midplane gravitational potential of truncated, flat power-law disks is extended to the whole physical space. It is shown that thickness effects (i.e. non-flatness) can be easily accounted for by implementing an appropriate "softening length" $\lambda$. The solution of this "softened ODE" has the following properties: i) it is regular at the edges (finite radial accelerations), ii) it possesses the correct long-range properties, iii) it matches the Newtonian potential of a geometrically thin disk very well, and iv) it tends continuously to the flat disk solution in the limit $\lambda \rightarrow 0$. As illustrated by many examples, the ODE, subject to exact Dirichlet conditions, can be solved numerically with efficiency for any given colatitude at second-order from center to infinity using radial mapping. This approach is therefore particularly well-suited to generating grids of gravitational forces in order to study particles moving under the field of a gravitating disk as found in various contexts (active nuclei, stellar systems, young stellar objects). Extension to non-power-law surface density profiles is straightforward through superposition. Grids can be produced upon request.

Context: Dark cloud chemical models usually predict large amounts of O2, often above observational limits. Aims: We investigate the reason for this discrepancy from a theoretical point of view, inspired by the studies of Jenkins and Whittet on oxygen depletion. Methods: We use the gas-grain code Nautilus with an up-to-date gas-phase network to study the sensitivity of the molecular oxygen abundance to the oxygen elemental abundance. We use the rate coefficient for the reaction O + OH at 10 K recommended by the KIDA (KInetic Database for Astrochemistry) experts. Results: The updates of rate coefficients and branching ratios of the reactions of our gas-phase chemical network, especially N + CN and H3+ + O, have changed the model sensitivity to the oxygen elemental abundance. In addition, the gas-phase abundances calculated with our gas-grain model are less sensitive to the elemental C/O ratio than those computed with a pure gas-phase model. The grain surface chemistry plays the role of a buffer absorbing most of the extra carbon. Finally, to reproduce the low abundance of molecular oxygen observed in dark clouds at all times, we need an oxygen elemental abundance smaller than 1.6E-4. Conclusion: The chemistry of molecular oxygen in dense clouds is quite sensitive to model parameters that are not necessarily well constrained. That O2 abundance may be sensitive to nitrogen chemistry is an indication of the complexity of interstellar chemistry.

We investigate the evolution of a system of two super-Earths with masses < 4 Earth masses embedded in a turbulent protoplanetary disk. The aim is to examine whether or not resonant trapping can occur and be maintained in presence of turbulence and how this depends on the amplitude of the stochastic density fluctuations in the disk. We have performed 2D numerical simulations using a grid-based hydrodynamical code in which turbulence is modelled as stochastic forcing. We assume that the outermost planet is initially located just outside the position of the 3:2 mean motion resonance (MMR) with the inner one and we study the dependance of the resonance stability with the amplitude of the stochastic forcing. For systems of two equal-mass planets we find that in disk models with an effective viscous stress parameter {\alpha} 10^{-3}, damping effects due to type I migration can counteract the effects of diffusion of the resonant angles, in such a way that the 3:2 MMR can possibly remain stable over the disk lifetime. For systems of super-Earths with mass ratio q=m_i/m_o < 1/2, where m_i (m_o) is the mass of the innermost (outermost) planet, the 3:2 MMR is broken in turbulent disks with effective viscous stresses 2x10^{-4}< {\alpha}< 10^{-3} but the planets become locked in stronger p+1:p resonances, with p increasing as the value for {\alpha} increases. For {\alpha}> 2x10^{-3}, the evolution can eventually involve temporary capture in a 8:7 commensurability but no stable MMR is formed. Our results suggest that for values of the viscous stress parameter typical to those generated by MHD turbulence, MMRs between two super-Earths are likely to be disrupted by stochastic density fluctuations. For lower levels of turbulence however, as is the case in presence of a dead-zone, resonant trapping can be maintained in systems with moderate values of the planet mass ratio.

The main sulfur-bearing molecules OCS, H2S, SO, SO2, and CS have been observed in four high mass dense cores (W43-MM1, IRAS 18264, IRAS 05358, and IRAS 18162). Our goal is to put some constraints on the relative evolutionary stage of these sources by comparing these observations with time-dependent chemical modeling. We used the chemical model Nahoon, which computes the gas-phase chemistry and gas-grain interactions of depletion and evaporation. Mixing of the different chemical compositions shells in a 1D structure through protostellar envelope has been included since observed lines suggest nonthermal supersonic broadening. Observed radial profiles of the temperature and density are used to compute the chemistry as a function of time. With our model, we underproduce CS by several orders of magnitude compared to the other S-bearing molecules, which seems to contradict observations, although some uncertainties in the CS abundance observed at high temperature remain. The OCS/SO2, SO/SO2, and H2S/SO2 abundance ratios could in theory be used to trace the age of these massive protostars since they show a strong dependence with time, but the sources are too close in age compared to the accuracy of chemical models and observations. Our comparison between observations and modeling may, however, indicate that W43-MM1 could be chemically younger than the three other sources. Turbulent diffusivity through the protostellar envelopes has to be less efficient than 2e14 cm2s-1. Otherwise, it would have smoothed out the abundance profiles, and this would have been observed. The sulfur chemistry depends strongly on the 1D physical conditions. In our case, no conclusion can be given on the relative age of IRAS 18264, IRAS 18162 and IRAS 05358 except that they are very close. W43-MM1 seems younger than the other sources.

Abridged: We detail and benchmark two sophisticated chemical models developed by the Heidelberg and Bordeaux astrochemistry groups. The main goal of this study is to elaborate on a few well-described tests for state-of-the-art astrochemical codes covering a range of physical conditions and chemical processes, in particular those aimed at constraining current and future interferometric observations of protoplanetary disks. We consider three physical models: a cold molecular cloud core, a hot core, and an outer region of a T Tauri disk. Our chemical network (for both models) is based on the original gas-phase osu_03_2008 ratefile and includes gas-grain interactions and a set of surface reactions for the H-, O-, C-, S-, and N-bearing molecules. The benchmarking is performed with the increasing complexity of the considered processes: (1) the pure gas-phase chemistry, (2) the gas-phase chemistry with accretion and desorption, and (3) the full gas-grain model with surface reactions. Using atomic initial abundances with heavily depleted metals and hydrogen in its molecular form, the chemical evolution is modeled within 10^9 years. The time-dependent abundances calculated with the two chemical models are essentially the same for all considered physical cases and for all species, including the most complex polyatomic ions and organic molecules. This result however required a lot of efforts to make all necessary details consistent through the model runs, e.g. definition of the gas particle density, density of grain surface sites, the strength and shape of the UV radiation field, etc. The reference models and the benchmark setup, along with the two chemical codes and resulting time-dependent abundances are made publicly available in the Internet: http://www.mpia.de/homes/semenov/Chemistry_benchmark/home.html

Because of new telescopes that will dramatically improve our knowledge of the interstellar medium, chemical models will have to be used to simulate the chemistry of many regions with diverse properties. To make these models more robust, it is important to understand their sensitivity to a variety of parameters. In this article, we report a study of the sensitivity of a chemical model of a cold dense core, with homogeneous and time-independent physical conditions, to variations in the following parameters: initial chemical inventory, gas temperature and density, cosmic-ray ionization rate, chemical reaction rate coefficients, and elemental abundances. From the results of the parameter variations, we can quantify the sensitivity of the model to each parameter as a function of time. Our results can be used in principle with observations to constrain some parameters for different cold clouds. We also attempted to use the Monte Carlo approach with all parameters varied collectively. Within the parameter ranges studied, the most critical parameters turn out to be the reaction rate coefficients at times up to 4e5 yr and elemental abundances at later times. At typical times of best agreement with observation, models are sensitive to both of these parameters. The models are less sensitive to other parameters such as the gas density and temperature. The improvement of models will require that the uncertainties in rate coefficients of important reactions be reduced. As the chemistry becomes better understood and more robust, it should be possible to use model sensitivities concerning other parameters, such as the elemental abundances and the cosmic ray ionization rate, to yield detailed information on cloud properties and history. Nevertheless, at the current stage, we cannot determine the best values of all the parameters simultaneously based on purely observational constraints.

We studied several representative circumstellar disks surrounding the Herbig Ae star MWC 480 and the T Tauri stars LkCa 15 and DM Tau at (sub-)millimeter wavelengths in lines of CCH. Our aim is to characterize photochemistry in the heavily UV-irradiated MWC 480 disk and compare the results to the disks around cooler T Tauri stars. We detected and mapped CCH in these disks with the IRAM Plateau de Bure Interferome- ter in the C- and D-configurations in the (1-0) and (2-1) transitions. Using an iterative minimization technique, the CCH column densities and excitation conditions are con- strained. Very low excitation temperatures are derived for the T Tauri stars. These values are compared with the results of advanced chemical modeling, which is based on a steady-state flared disk structure with a vertical temperature gradient, and a gas- grain chemical network with surface reactions. Both model and observations suggest that CCH is a sensitive tracer of the X-ray and UV irradiation. The predicted radial dependency and source to source variations of CCH column densities qualitatively agree with the observed values, but the predicted column densities are too low by a factor of several. The chemical model fails to reproduce high concentrations of CCH in very cold disk midplane as derived from the observed low excitation condition for both the (1-0) and (2-1) transitions.

The determination of Saturn's atmospheric noble gas abundances are critical to understanding the formation and evolution of Saturn, and giant planets in general. These measurements can only be performed with an entry probe. A Saturn probe will address whether enhancement in heavy noble gases, as was found in Jupiter, are a general feature of giant planets, and their ratios will be a powerful constraint on how they form. The helium abundance will show the extent to which helium has phase separated from hydrogen in the planet's deep interior. Jupiter's striking neon depletion may also be tied to its helium depletion, and must be confirmed or refuted in Saturn. Together with Jupiter's measured atmospheric helium abundance, a consistent evolutionary theory for both planets, including "helium rain" will be possible. We will then be able to calibrate the theory of the evolution of all giant planets, including exoplanets. In addition, high pressure H/He mixtures under giant planet conditions are an important area of condensed matter physics that are beyond the realm of experiment.

Aims. Collections of dust, grains, and planetesimals are often treated as a pressureless fluid. We study the validity of neglecting the pressure of such a fluid by computing it exactly for the case of particles settling in a disk. Methods. We solve a modified collisionless Boltzmann equation for the particles and compute the corresponding moments of the phase space distribution: density, momentum, and pressure. Results. We find that whenever the Stokes number, defined as the ratio of the gas drag timescale to the orbital timescale, is more than 1/2, the particle fluid cannot be considered as pressureless. While we show it only in the simple case of particles settling in a laminar disk, this property is likely to remain true for most flows, including turbulent flows.

We present the exact calculus of the gravitational potential and acceleration along the symmetry axis of a plane, homogeneous, polar cell as a function of mean radius a, radial extension e, and opening angle f. Accurate approximations are derived in the limit of high numerical resolution at the geometrical mean <a> of the inner and outer radii (a key-position in current FFT-based Poisson solvers). Our results are the full extension of the approximate formula given in the textbook of Binney & Tremaine to all resolutions. We also clarify definitely the question about the existence (or not) of self-forces in polar cells. We find that there is always a self-force at radius <a> except if the shape factor a.f/e reaches ~ 3.531, asymptotically. Such cells are therefore well suited to build a polar mesh for high resolution simulations of self-gravitating media in two dimensions. A by-product of this study is a newly discovered indefinite integral involving complete elliptic integral of the first kind over modulus.

Aims. We attempt to understand the presence of gas phase CO below its freezing temperature in circumstellar disks. We study two promising mechanisms to explain this phenomenon: turbulent mixing and photodesorption. Methods. We compute the chemical evolution of circumstellar disks including grain surface reactions with and without turbulent mixing and CO photodesorption. Results. We show that photodesorption significantly enhances the gas phase CO abundance, by extracting CO from the grains when the visual extinction remains below about 5 magnitudes. However the resulting dependence of column density on radial distance is not consistent with observations so far. We propose that this inconsistency could be the result of grain growth. On the other hand, the influence of turbulent mixing is not found to be straightforward. The efficiency of turbulent mixing depends upon a variety of parameters, including the disk structure. For the set of parameters we chose, turbulent mixing is not found to have any significant influence on the CO column density.

We study the molecular content and chemistry of a circumstellar disk surrounding the Herbig Ae star AB Aur at (sub-)millimeter wavelengths. Our aim is to reconstruct the chemical history and composition of the AB Aur disk and to compare it with disks around low-mass, cooler T Tauri stars. We observe the AB Aur disk with the IRAM Plateau de Bure Interferometer in the C- and D- configurations in rotational lines of CS, HCN, C2H, CH3OH, HCO+, and CO isotopes. Using an iterative minimization technique, observed columns densities and abundances are derived. These values are further compared with results of an advanced chemical model that is based on a steady-state flared disk structure with a vertical temperature gradient, and gas-grain chemical network with surface reactions. We firmly detect HCO+ in the 1--0 transition, tentatively detect HCN, and do not detect CS, C2H, and CH3OH. The observed HCO+ and 13CO column densities as well as the upper limits to the column densities of HCN, CS, C2H, and CH3OH are in good agreement with modeling results and those from previous studies. The AB Aur disk possesses more CO, but is less abundant in other molecular species compared to the DM Tau disk. This is primarily caused by intense UV irradiation from the central Herbig A0 star, which results in a hotter disk where CO freeze out does not occur and thus surface formation of complex CO-bearing molecules might be inhibited.

The first-order ordinary differential equation (ODE) that describes the mid-plane gravitational potential in flat finite size discs in which the surface density is a power-law function of the radius R with exponent s (Hur\'e & Hersant 2007) is solved exactly in terms of infinite series. The formal solution of the ODE is derived and then converted into a series representation by expanding the elliptic integral of the first kind over its modulus before analytical integration. Inside the disc, the gravitational potential consists of three terms: a power law of radius R with index 1+s, and two infinite series of the variables R and 1/R. The convergence of the series can be accelerated, enabling the construction of reliable approximations. At the lowest-order, the potential inside large astrophysical discs (s ~ -1.5 +/- 1) is described by a very simple formula whose accuracy (a few percent typically) is easily increased by considering successive orders through a recurrence. A basic algorithm is given. Applications concern all theoretical models and numerical simulations where the influence of disc gravity must be checked and/or reliably taken into account.

We show that the gravitational potential in the plane of an axisymmetrical flat disk where the surface density varies as a power of the radius obeys an inhomogeneous first-order Ordinary Differential Equation (ODE) solvable by standard techniques. The potential being exactly known at the origin for any power index (and at infinity as well), the search for solutions consists of a Two-point Boundary Value Problem (TBVP) with Dirichlet conditions. The computating time is then linear with the number of grid points, instead of quadratic from direct summation methods. Complex mass distributions which can be decomposed into a mixture of power law surface density profiles are easily accessible through the superposition principle. This ODE definitively takes the place of the untractable bidimensional Poisson equation for planar calculations. It opens new horizons to investigate various aspects related to self-gravity in astrophysical disks (force calculations, stability analysis, etc.).

We discuss the possibility that astrophysical accretion disks are dynamically unstable to non-axisymmetric disturbances with characteristic scales much smaller than the vertical scale height. The instability is studied using three methods: one based on the energy integral, which allows the determination of a sufficient condition of stability, one using a WKB approach, which allows the determination of the necessary and sufficient condition for instability and a last one by numerical solution. This linear instability occurs in any inviscid stably stratified differential rotating fluid for rigid, stress-free or periodic boundary conditions, provided the angular velocity $\Omega$ decreases outwards with radius $r$. At not too small stratification, its growth rate is a fraction of $\Omega$. The influence of viscous dissipation and thermal diffusivity on the instability is studied numerically, with emphasis on the case when $d \ln \Omega / d \ln r =-3/2$ (Keplerian case). Strong stratification and large diffusivity are found to have a stabilizing effect. The corresponding critical stratification and Reynolds number for the onset of the instability in a typical disk are derived. We propose that the spontaneous generation of these linear modes is the source of turbulence in disks, especially in weakly ionized disks.

We investigate the analogy between circumstellar disks and the Taylor-Couette flow. Using the Reynolds similarity principle, the analogy results in a number of parameter-free predictions about stability of the disks, and their turbulent transport properties, provided the disk structure is available. We discuss how the latter can be deduced from interferometric observations of circumstellar material. We use the resulting disk structure to compute the molecular transport coefficients, including the effect of ionization by the central object. The resulting control parameter indicates that the disk is well into the turbulent regime. The analogy is also used to compute the effective accretion rate, as a function of the disk characteristic parameters (orbiting velocity, temperature and density). These values are in very good agreement with experimental, parameter-free predictions derived from the analogy. The turbulent viscosity is also computed and found to correspond to an $\alpha$-parameter $2\times 10^{-4}<\alpha<2\times 10^{-2}$. Predictions regarding fluctuations are also checked: luminosity fluctuations in disks do obey the same universal distribution as energy fluctuations observed in a laboratory turbulent flow. Radial velocity dispersion in the outer part of the disk is predicted to be of the order of 0.1 km/s, in agreement with available observations. All these issues provide a proof of the turbulent character of the circumstellar disks, as well as a parameter-free theoretical estimate of effective accretion rates.

In this paper, we show that astrophysical accretion disks are dynamically unstable to non-axisymmetric disturbances. This instability is present in any stably stratified anticyclonically sheared flow as soon as the angular velocity decreases outwards. In the large Froude number limit, the maximal growth rate is proportional to the angular rotation velocity, and is independent of the stratification. In the low Froude number limit, it decreases like the inverse of the Froude number, thereby vanishing for unstratified, centrigugally stable flows. The instability is not sensitive to disk boundaries. We discuss the possible significance of our result, and its implications on the turbulent state achieved by the disks. We conclude that this linear instability is one of the best candidates for the source of turbulence in geometrically thin disks, and that magnetic fields can be safely ignored when studying their turbulent state. The relevance of the instability for thick disks or nearly neutrally stratified disks remains to be explored.