BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

F. Hersant

F. Hersant appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.EPastro-ph.GAastro-ph.SRastro-ph.IMastro-ph.HEcond-mat.otherphysics.flu-dyn

Catalog footprint

What is connected

24works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

24 published item(s)

preprint2020arXiv

GG Tau A: gas properties and dynamics from the cavity to the outer disk

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Chemical and physical characterization of collapsing low-mass prestellar dense cores

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Importance of the H2 abundance in protoplanetary disk ices for the molecular layer chemical composition

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

2D photochemical modeling of Saturn's stratosphere. Part II: Feedback between composition and temperature

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Chemistry in Protoplanetary Disks: the gas-phase CO/H2 ratio and the Carbon reservoir

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The EChO science case

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Modeling Complex Organic Molecules in dense regions: Eley-Rideal and complex induced reaction

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

First submillimeter observation of CO in the stratosphere of Uranus

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The effect of rotation and tidal heating on the thermal lightcurves of Super Mercuries

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A KInetic Database for Astrochemistry (KIDA)

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Chemistry in Disks VIII: the CS molecule as an analytic tracer of turbulence in disks

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

CID: Chemistry In Disks VII. First detection of HC3N in protoplanetary disks

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).

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Protoplanetary migration in non-isothermal discs with turbulence driven by stochastic forcing

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The impact of atmospheric circulation on the chemistry of the hot Jupiter HD 209458b

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The potential of discs from a "mean Green function"

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Momentum transport and torque scaling in Taylor-Couette flow from an analogy with turbulent convection

We generalize an analogy between rotating and stratified shear flows. This analogy is summarized in Table 1. We use this analogy in the unstable case (centrifugally unstable flow v.s. convection) to compute the torque in Taylor-Couette configuration, as a function of the Reynolds number. At low Reynolds numbers, when most of the dissipation comes from the mean flow, we predict that the non-dimensional torque $G=T/ν^2L$, where $L$ is the cylinder length, scales with Reynolds number $R$ and gap width $η$, $G=1.46 η^{3/2} (1-η)^{-7/4}R^{3/2}$. At larger Reynolds number, velocity fluctuations become non-negligible in the dissipation. In these regimes, there is no exact power law dependence the torque versus Reynolds. Instead, we obtain logarithmic corrections to the classical ultra-hard (exponent 2) regimes: $$ G=0.50\frac{η^{2}}{(1-η)^{3/2}}\frac{R^{2}}{\ln[η^2(1-η)R^ 2/10^4]^{3/2}}.$$ These predictions are found to be in excellent agreement with available experimental data. Predictions for scaling of velocity fluctuations are also provided.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On the dynamics of resonant super-Earths in disks with turbulence driven by stochastic forcing

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 α 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}< α< 10^{-3} but the planets become locked in stronger p+1:p resonances, with p increasing as the value for α increases. For α> 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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Oxygen depletion in dense molecular clouds: a clue to a low O2 abundance?

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Sulfur chemistry: 1D modeling in massive dense cores

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Thermal phase curves of non-transiting terrestrial exoplanets 2. Characterizing airless planets

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 μm range. [...] As inclination above 60° 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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Chemistry in Disks. III. -- Photochemistry and X-ray driven chemistry probed by the ethynyl radical (CCH) in DM Tau, LkCa 15, and MWC 480

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Sensitivity analyses of dense cloud chemical models

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

On the pressure of collisionless particle fluids. The case of solids settling in disks

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Self-gravity at the scale of the polar cell

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.

0 citations0 reviewsNo code linkedOpen access