BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Franck Hersant

Franck Hersant contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.EPastro-ph.GAastro-ph.SRastro-ph.IMastro-ph.COphysics.space-ph

Trust snapshot

Quick read

Trust 21 - Emerging
23works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

23 published item(s)

preprint2016arXiv

Gas and grain chemical composition in cold cores as predicted by the Nautilus 3-phase model

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

2D photochemical modeling of Saturn's stratosphere. Part I: Seasonal variation of atmospheric composition without meridional transport

Saturn's axial tilt of 26.7° 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° in both hemispheres

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Accretion of Uranus and Neptune from inward-migrating planetary embryos blocked by Jupiter and Saturn

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Evidence for DCO+ as a probe of ionization in the warm disk surface

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Gas giant planets as dynamical barriers to inward-migrating super-Earths

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Mercury-T: A new code to study tidally evolving multi-planet systems. Applications to Kepler-62

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Formation, tidal evolution and habitability of the Kepler-186 system

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Grain Surface Reactions in Molecular Clouds: The Effect of Cosmic Rays and Quantum Tunneling

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Hot super-Earths and giant planet cores from different migration histories

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Physical and chemical structure of planet-forming disks probed by millimeter observations and modeling

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Pseudo 2D chemical model of hot Jupiter atmospheres: application to HD 209458b and HD 189733b

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Self-gravity in curved mesh elements

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Self-gravity in thin discs and edge effects: an extension of Paczynski's approximation

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 πG Σ(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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dynamical and collisional constraints on a stochastic late veneer on the terrestrial planets

Given their tendency to be incorporated into the core during differentiation, the highly-siderophile elements (HSEs) in Earth&#39;s mantle are thought to have been accreted as a `late veneer&#39; 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&#39; 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&#39;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&#39;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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Survival of interstellar molecules to prestellar dense core collapse and early phases of disk formation

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Tidal dissipation and eccentricity pumping: Implications for the depth of the secondary eclipse of 55 Cnc e

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&#39;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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Water in Protoplanetary Disks: Deuteration and Turbulent Mixing

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A chemical model for the atmosphere of hot Jupiters

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$μ$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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Chemistry in the First Hydrostatic Core Stage By Adopting Three-Dimensional Radiation Hydrodynamic Simulations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

From Prestellar to Protostellar Cores II. Time Dependence and Deuterium Fractionation

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

AGN disks and black holes on the weighting scales

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

CID: Chemistry in disks VI.sulfur-bearing molecules in the protoplanetary disks surrounding LkCa15, MWC480, DM Tau, and GO Tau

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Newtonian potential of thin disks

The one-dimensional, ordinary differential equation (ODE) by Huré & 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 &#34;softening length&#34; $λ$. The solution of this &#34;softened ODE&#34; 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 $λ\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.

0 citations0 reviewsNo code linkedOpen access