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

Gilles Chabrier

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

ResearcherUnclaimed source record
astro-ph.SRastro-ph.EPastro-ph.GAastro-phastro-ph.IMphysics.comp-phphysics.flu-dynastro-ph.COphysics.ao-ph

Catalog footprint

What is connected

38works
9topics
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

38 published item(s)

preprint2022arXiv

Statistical properties and correlation length in star-forming molecular clouds: I. Formalism and application to observations

The proper characterization of the general statistical behavior of these fluctuations, from a limited sample of observations or simulations, is of prime importance to understand the process of star formation. In this article, we use the ergodic theory for any random field of fluctuations, as commonly used in statistical physics, to derive rigorous statistical results. We outline how to evaluate the autocovariance function (ACF) and the characteristic correlation length of these fluctuations. We then apply this statistical approach to astrophysical systems characterized by a field of density fluctuations, notably star-forming clouds. When it is difficult to determine the correlation length from the empirical ACF, we show alternative ways to estimate the correlation length. We show that the statistics of the column-density field is hampered by biases introduced by integration effects along the line of sight and we explain how to reduce these biases. The statistics of the probability density function (PDF) ergodic estimator also yields the derivation of the proper statistical error bars. We provide a method that can be used by observers and numerical simulation specialists to determine the latter. We show that they (i) cannot be derived from simple Poisson statistics and (ii) become increasingly large for increasing density contrasts, severely hampering the accuracy of the low and high end part of the PDF because of a sample size that is too small. As templates of various stages of star formation in MCs, we then examine the case of the Polaris and Orion B clouds in detail. We calculate, from the observations, the ACF and the correlation length in these clouds and show that the latter is on the order of $\sim$1\% of the size of the cloud.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Statistical properties and correlation length in star-forming molecular clouds: II. Gravitational potential and virial parameter

In the first article of this series, we have used the ergodic theory to assess the validity of a statistical approach to characterize various properties of star-forming molecular clouds (MCs) from a limited number of observations or simulations. This allows the proper determination of confidence intervals for various volumetric averages of statistical quantities obtained form observations or numerical simulations. In this joint paper, we apply the same formalism to a different kind of (observational or numerical) study of MCs. Indeed, as observations cannot fully unravel the complexity of the inner density structure of star forming clouds, it is important to know whether global observable estimates, such as the total mass and size of the cloud, can give an accurate estimation of various key physical quantities that characterize the dynamics of the cloud. Of prime importance is the correct determination of the total gravitational (binding) energy and virial parameter of a cloud. We show that, whereas for clouds that are not in a too advanced stage of star formation, such as Polaris or Orion B, the knowledge of only their mass and size is sufficient to yield an accurate determination of the aforementioned quantities from observations (i.e. in real space). In contrast, we show that this is no longer true for numerical simulations in a periodic box. We derive a relationship for the ratio of the virial parameter in these two respective cases.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A library of ATMO forward model transmission spectra for hot Jupiter exoplanets

We present a grid of forward model transmission spectra, adopting an isothermal temperature-pressure profile, alongside corresponding equilibrium chemical abundances for 117 observationally significant hot exoplanets (Equilibrium Temperatures of 547-2710 K). This model grid has been developed using a 1D radiative-convective-chemical equilibrium model termed ATMO, with up-to-date high temperature opacities. We present an interpretation of observations of ten exoplanets, including best fit parameters and $χ^{2}$ maps. In agreement with previous works, we find a continuum from clear to hazy/cloudy atmospheres for this sample of hot Jupiters. The data for all the 10 planets are consistent with sub-solar to solar C/O ratio, 0.005 to 10 times solar metallicity and water rather than a methane dominated infrared spectra. We then explore the range of simulated atmospheric spectra for different exoplanets, based on characteristics such as temperature, metallicity, C/O-ratio, haziness and cloudiness. We find a transition value for the metallicity between 10 and 50 times solar, which leads to substantial changes in the transmission spectra. We also find a transition value of C/O ratio, from water to carbon species dominated infrared spectra, as found by previous works, revealing a temperature dependence of this transition point ranging from $\sim$0.56 to $\sim$1-1.3 for equilibrium temperatures from $\sim$900 to $\sim$2600 K. We highlight the potential of the spectral features of HCN and C$_2$H$_2$ to constrain the metallicities and C/O ratios of planets, using JWST observations. Finally, our entire grid ($\sim$460,000 simulations) is publicly available and can be used directly with the JWST simulator PandExo for planning observations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A new set of atmosphere and evolution models for cool T-Y brown dwarfs and giant exoplanets

We present a new set of solar metallicity atmosphere and evolutionary models for very cool brown dwarfs and self-luminous giant exoplanets, which we term ATMO 2020. Atmosphere models are generated with our state-of-the-art 1D radiative-convective equilibrium code ATMO, and are used as surface boundary conditions to calculate the interior structure and evolution of $0.001-0.075\,\mathrm{M_{\odot}}$ objects. Our models include several key improvements to the input physics used in previous models available in the literature. Most notably, the use of a new H-He equation of state including ab initio quantum molecular dynamics calculations has raised the mass by $\sim1-2\%$ at the stellar-substellar boundary and has altered the cooling tracks around the hydrogen and deuterium burning minimum masses. A second key improvement concerns updated molecular opacities in our atmosphere model ATMO, which now contains significantly more line transitions required to accurately capture the opacity in these hot atmospheres. This leads to warmer atmospheric temperature structures, further changing the cooling curves and predicted emission spectra of substellar objects. We present significant improvement for the treatment of the collisionally broadened potassium resonance doublet, and highlight the importance of these lines in shaping the red-optical and near-infrared spectrum of brown dwarfs. We generate three different grids of model simulations, one using equilibrium chemistry and two using non-equilibrium chemistry due to vertical mixing, all three computed self-consistently with the pressure-temperature structure of the atmosphere. We show the impact of vertical mixing on emission spectra and in colour-magnitude diagrams, highlighting how the $3.5-5.5\,\mathrm{μm}$ flux window can be used to calibrate vertical mixing in cool T-Y spectral type objects.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Protostellar collapse: regulation of the angular momentum and onset of an ionic precursor

Through the magnetic braking and the launching of protostellar outflows, magnetic fields play a major role in the regulation of angular momentum in star formation, which directly impacts the formation and evolution of protoplanetary disks and binary systems. The aim of this paper is to quantify those phenomena in the presence of non-ideal magnetohydrodynamics effects, namely the Ohmic and ambipola r diffusion. We perform three-dimensional simulations of protostellar collapses varying the mass of the prestellar dense core, the thermal support (the $α$ ratio) and the dust grain size-distribu tion. The mass mostly influences the magnetic braking in the pseudo-disk, while the thermal support impacts the accretion rate and hence the properties of the disk. Removing the grains smaller than 0. 1 $μ$m in the Mathis, Rumpl, Nordsieck (MRN) distribution enhances the ambipolar diffusion coefficient. Similarly to previous studies, we find that this change in the distribution reduces the magnet ic braking with an impact on the disk. The outflow is also significantly weakened. In either case, the magnetic braking largely dominates the outflow as a process to remove the angular momentum from t he disk. Finally, we report a large ionic precursor to the outflow with velocities of several km s$^{-1}$, which may be observable.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Parallax of VHS J1256-1257 from CFHT and Pan-STARRS 1

We present new parallax measurements from the CFHT Infrared Parallax Program and the Pan-STARRS 3$π$ Steradian Survey for the young ($\approx150-300$ Myr) triple system VHS J125601.92$-$125723.9. This system is composed of a nearly equal-flux binary ("AB") and a wide, possibly planetary-mass companion ("b"). The system's published parallactic distance ($12.7\pm1.0$ pc) implies absolute magnitudes unusually faint compared to known young objects and is in tension with the spectrophotometric distance for the central binary ($17.2\pm2.6$ pc). Our CFHT and Pan-STARRS parallaxes are consistent, and the more precise CFHT result places VHS J1256-1257 at $22.2^{+1.1}_{-1.2}$ pc. Our new distance results in higher values for the companion's mass ($19\pm5$ M$_{\rm Jup}$) and temperature ($1240\pm50$ K), and also brings the absolute magnitudes of all three components into better agreement with known young objects.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A chemical solver to compute molecule and grain abundances and non-ideal MHD resistivities in prestellar core collapse calculations

We develop a detailed chemical network relevant to the conditions characteristic of prestellar core collapse. We solve the system of time-dependent differential equations to calculate the equilibrium abundances of molecules and dust grains, with a size distribution given by size-bins for these latter. These abundances are used to compute the different non-ideal magneto-hydrodynamics resistivities (ambipolar, Ohmic and Hall), needed to carry out simulations of protostellar collapse. For the first time in this context, we take into account the evaporation of the grains, the thermal ionisation of Potassium, Sodium and Hydrogen at high temperature, and the thermionic emission of grains in the chemical network, and we explore the impact of various cosmic ray ionisation rates. All these processes significantly affect the non-ideal magneto-hydrodynamics resistivities, which will modify the dynamics of the collapse. Ambipolar diffusion and Hall effect dominate at low densities, up to n_H = 10^12 cm^-3, after which Ohmic diffusion takes over. We find that the time-scale needed to reach chemical equilibrium is always shorter than the typical dynamical (free fall) one. This allows us to build a large, multi-dimensional multi-species equilibrium abundance table over a large temperature, density and ionisation rate ranges. This table, which we make accessible to the community, is used during first and second prestellar core collapse calculations to compute the non-ideal magneto-hydrodynamics resistivities, yielding a consistent dynamical-chemical description of this process.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Ambipolar diffusion in low-mass star formation. I. General comparison with the ideal MHD case

In this paper, we provide a more accurate description of the evolution of the magnetic flux redistribution during prestellar core collapse by including resistive terms in the magnetohydrodynamics (MHD) equations. We focus more particularly on the impact of ambipolar diffusion. We use the adaptive mesh refinement code RAMSES to carry out such calculations. The resistivities required to calculate the ambipolar diffusion terms were computed using a reduced chemical network of charged, neutral and grain species. The inclusion of ambipolar diffusion leads to the formation of a magnetic diffusion barrier in the vicinity of the core, preventing accumulation of magnetic flux in and around the core and amplification of the field above 0.1G. The mass and radius of the first Larson core remain similar between ideal and non-ideal MHD models. This diffusion plateau has crucial consequences on magnetic braking processes, allowing the formation of disk structures. Magnetically supported outflows launched in ideal MHD models are weakened when using non-ideal MHD. Contrary to ideal MHD misalignment between the initial rotation axis and the magnetic field direction does not significantly affect the results for a given mu, showing that the physical dissipation truly dominate over numerical diffusion. We demonstrate severe limits of the ideal MHD formalism, which yield unphysical behaviours in the long-term evolution of the system. This includes counter rotation inside the outflow, interchange instabilities, and flux redistribution triggered by numerical diffusion, none observed in non-ideal MHD. Disks with Keplerian velocity profiles form in all our non-ideal MHD simulations, with final mass and size which depend on the initial magnetisation. This ranges from a few 0.01 solar masses and 20-30 au for the most magnetised case (mu=2) to 0.2 solar masses and 40-80 au for a lower magnetisation (mu=5).

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Early Results from VLT-SPHERE: Long-Slit Spectroscopy of 2MASS 0122-2439B, a Young Companion Near the Deuterium Burning Limit

We present 0.95-1.80 $μ$m spectroscopy of the $\sim$12-27 $M_{\rm Jup}$ companion orbiting the faint ($R$$\sim$13.6), young ($\sim$120 Myr) M-dwarf 2MASS J01225093--2439505 ("2M0122--2439 B") at 1.5 arcsecond separation (50 AU). Our coronagraphic long-slit spectroscopy was obtained with the new high contrast imaging platform VLT-SPHERE during Science Verification. The unique long-slit capability of SPHERE enables spectral resolution an order of magnitude higher than other extreme AO exoplanet imaging instruments. With a low mass, cool temperature, and very red colors, 2M0122-2439 B occupies a particularly important region of the substellar color-magnitude diagram by bridging the warm directly imaged hot planets with late-M/early-L spectral types (e.g. $β$ Pic b and ROXs 42Bb) and the cooler, dusty objects near the L/T transition (e.g. HR 8799bcde and 2MASS 1207b). We fit BT-Settl atmospheric models to our $R$$\approx$350 spectrum and find $T_{\rm eff}$=1600$\pm$100 K and $\log(g)$=4.5$\pm$0.5 dex. Visual analysis of our 2M0122-2439 B spectrum suggests a spectral type L3-L4, and we resolve shallow $J$-band alkali lines, confirming its low gravity and youth. Specifically, we use the Allers & Liu (2013) spectral indices to quantitatively measure the strength of the FeH, VO, KI, spectral features, as well as the overall $H$-band shape. Using these indices, along with the visual spectral type analysis, we classify 2M0122-2439 B as an intermediate gravity (INT-G) object with spectral type L3.7$\pm$1.0.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Theoretical limits on magnetic field strengths in low-mass stars

Observations have suggested that some low-mass stars have larger radii than predicted by 1-D structure models. Some theoretical models have invoked very strong interior magnetic fields (of order 1 MG or more) as a possible cause of such large radii. Whether fields of that strength could in principle by generated by dynamo action in these objects is unclear, and we do not address the matter directly. Instead, we examine whether such fields could remain in the interior of a low mass object for a significant time, and whether they would have any other obvious signatures. First, we estimate timescales for the loss of strong fields by magnetic buoyancy instabilities. We consider a range of field strengths and simple morphologies, including both idealized flux tubes and smooth layers of field. We confirm some of our analytical estimates using thin flux tube magnetohydrodynamic (MHD) simulations of the rise of buoyant fields in a fully-convective M-dwarf. Separately, we consider the Ohmic dissipation of such fields. We find that dissipation provides a complementary constraint to buoyancy: while small-scale, fibril fields might be regenerated faster than they rise, the dissipative heating associated with such fields would in some cases greatly exceed the luminosity of the star. We show how these constraints combine to yield limits on the internal field strength and morphology in low-mass stars. In particular, we find that for stars of 0.3 solar masses, no fields in flux tubes stronger than about 800 kG are simultaneously consistent with both constraints.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the role of the H2 ortho:para ratio in gravitational collapse during star formation

Hydrogen molecules (H2) come in two forms in the interstellar medium, ortho- and para-hydrogen, corresponding to the two different spin configurations of the two hydrogen atoms. The relative abundances of the two flavours in the interstellar medium are still very uncertain, and this abundance ratio has a significant impact on the thermal properties of the gas. In the context of star formation, theoretical studies have recently adopted two different strategies when considering the ortho:para ratio (OPR) of H2 molecules; the first considers the OPR to be frozen at 3:1 while the second assumes that the species are in thermal equilibrium. As the OPR potentially affects the protostellar cores which form as a result of the gravitational collapse of a dense molecular cloud, the aim of this paper is to quantify precisely what role the choice of OPR plays in the properties and evolution of the cores. We used two different ideal gas equations of state for a hydrogen and helium mix in a radiation hydrodynamics code to simulate the collapse of a dense cloud and the formation of the first and second Larson cores; the first equation of state uses a fixed OPR of 3:1 while the second assumes thermal equilibrium. Simulations using an equilibrium ratio collapse faster at early times and show noticeable oscillations around hydrostatic equilibrium, to the point where the core expands for a short time right after its formation before resuming its contraction. In the case of a fixed 3:1 OPR, the core's evolution is a lot smoother. The OPR was however found to have little impact on the size, mass and radius of the two Larson cores. We conclude that if one is solely interested in the final properties of the cores when they are formed, it does not matter which OPR is used. On the other hand, if one's focus lies primarily in the evolution of the first core, the choice of OPR becomes important.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Mass-Dependence of Angular Momentum Evolution in Sun-Like Stars

To better understand the observed distributions of rotation rate and magnetic activity of sun-like and low-mass stars, we derive a physically motivated scaling for the dependence of the stellar-wind torque on Rossby number. The torque also contains an empirically-derived scaling with stellar mass (and radius), which provides new insight into the mass-dependence of stellar magnetic and wind properties. We demonstrate that this new formulation explains why the lowest mass stars are observed to maintain rapid rotation for much longer than solar-mass stars, and simultaneously, why older populations exhibit a sequence of slowly rotating stars, in which the low-mass stars rotate more slowly than solar-mass stars. The model also reproduces some previously unexplained features in the period-mass diagram for the Kepler field, notably: the particular shape of the "upper envelope" of the distribution, suggesting that ~95% of Kepler field stars with measured rotation periods are younger than ~4 Gyr; and the shape of the "lower envelope," corresponding to the location where stars transition between magnetically saturated and unsaturated regimes.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Analytical theory for the initial mass function: III time dependence and star formation rate

The present paper extends our previous theory of the stellar initial mass function (IMF) by including the time-dependence, and by including the impact of magnetic field. The predicted mass spectra are similar to the time independent ones with slightly shallower slopes at large masses and peak locations shifted toward smaller masses by a factor of a few. Assuming that star-forming clumps follow Larson type relations, we obtain core mass functions in good agreement with the observationally derived IMF, in particular when taking into account the thermodynamics of the gas. The time-dependent theory directly yields an analytical expression for the star formation rate (SFR) at cloud scales. The SFR values agree well with the observational determinations of various Galactic molecular clouds. Furthermore, we show that the SFR does not simply depend linearly on density, as sometimes claimed in the literature, but depends also strongly on the clump mass/size, which yields the observed scatter. We stress, however, that {\it any} SFR theory depends, explicitly or implicitly, on very uncertain assumptions like clump boundaries or the mass of the most massive stars that can form in a given clump, making the final determinations uncertain by a factor of a few. Finally, we derive a fully time-dependent model for the IMF by considering a clump, or a distribution of clumps accreting at a constant rate and thus whose physical properties evolve with time. In spite of its simplicity, this model reproduces reasonably well various features observed in numerical simulations of converging flows. Based on this general theory, we present a paradigm for star formation and the IMF.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Layered convection as the origin of Saturn's luminosity anomaly

As they keep cooling and contracting, Solar System giant planets radiate more energy than they receive from the Sun. Applying the first and second principles of thermodynamics, one can determine their cooling rate, luminosity, and temperature at a given age. Measurements of Saturn's infrared intrinsic luminosity, however, reveal that this planet is significantly brighter than predicted for its age. This excess luminosity is usually attributed to the immiscibility of helium in the hydrogen-rich envelope, leading to "rains" of helium-rich droplets. Existing evolution calculations, however, suggest that the energy released by this sedimentation process may not be sufficient to resolve the puzzle. Here, we demonstrate using planetary evolution models that the presence of layered convection in Saturn's interior, generated, like in some parts of Earth oceans, by the presence of a compositional gradient, significantly reduces its cooling. It can explain the planet's present luminosity for a wide range of configurations without invoking any additional source of energy. This suggests a revision of the conventional homogeneous adiabatic interior paradigm for giant planets, and questions our ability to assess their heavy element content. This reinforces the possibility for layered convection to help explaining the anomalously large observed radii of extrasolar giant planets.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Star Formation Rate of Molecular Clouds

We review recent advances in the analytical and numerical modeling of the star formation rate in molecular clouds and discuss the available observational constraints. We focus on molecular clouds as the fundamental star formation sites, rather than on the larger-scale processes that form the clouds and set their properties. Molecular clouds are shaped into a complex filamentary structure by supersonic turbulence, with only a small fraction of the cloud mass channeled into collapsing protostars over a free-fall time of the system. In recent years, the physics of supersonic turbulence has been widely explored with computer simulations, leading to statistical models of this fragmentation process, and to the prediction of the star formation rate as a function of fundamental physical parameters of molecular clouds, such as the virial parameter, the rms Mach number, the compressive fraction of the turbulence driver, and the ratio of gas to magnetic pressure. Infrared space telescopes, as well as ground-based observatories have provided unprecedented probes of the filamentary structure of molecular clouds and the location of forming stars within them.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A new vision on giant planet interiors: the impact of double diffusive convection

While conventional interior models for Jupiter and Saturn are based on the simplistic assumption of a solid core surrounded by a homogeneous gaseous envelope, we derive new models with an inhomogeneous distribution of heavy elements, i.e. a gradient of composition, within these planets. Such a compositional stratification hampers large scale convection which turns into double-diffusive convection, yielding an inner thermal profile which departs from the traditionally assumed adiabatic interior, affecting these planet heat content and cooling history. To address this problem, we develop an analytical approach of layered double-diffusive convection and apply this formalism to Solar System gaseous giant planet interiors. These models satisfy all observational constraints and yield a metal enrichment for our gaseous giants up to 30 to 60% larger than previously thought. The models also constrain the size of the convective layers within the planets. As the heavy elements tend to be redistributed within the gaseous envelope, the models predict smaller than usual central cores inside Saturn and Jupiter, with possibly no core for this latter. These models open a new window and raise new challenges on our understanding of the internal structure of giant (solar and extrasolar) planets, in particular on the determination of their heavy material content, a key diagnostic for planet formation theories.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Observed luminosity spread in young clusters and Fu Ori stars: a unified picture

The idea that non steady accretion during the embedded phase of protostar evolution can produce the observed luminosity spread in the Herzsprung-Russell diagram (HRD) of young clusters has recently been called into question. Observations of Fu Ori, for instance, suggest an expansion of the star during strong accretion events whereas the luminosity spread implies a contraction of the accreting objects, decreasing their radiating surface. In this paper, we present a global scenario based on calculations coupling episodic accretion histories derived from numerical simulations of collapsing cloud prestellar cores of various masses and subsequent protostar evolution. Our calculations show that, assuming an initial protostar mass $\mi \sim 1\,\mjup$, typical of the second Larson's core, both the luminosity spread in the HRD and the inferred properties of Fu Ori events (mass, radius, accretion rate) can be explained by this scenario, providing two conditions. First, there must be some variation within the fraction of accretion energy absorbed by the protostar during the accretion process. Second the range of this variation should increase with increasing accretion burst intensity, and thus with the initial core mass and final star mass. The numerical hydrodynamics simulations of collapsing cloud prestellar cores indeed show that the intensity of the accretion bursts correlates with the mass and initial angular momentum of the prestellar core. Massive prestellar cores with high initial angular momentum are found to produce intense bursts characteristic of Fu Ori like events. Our results thus suggest a link between the burst intensities and the fraction of accretion energy absorbed by the protostar, with some threshold in the accretion rate, of the order of $10^{-5}\msolyr$, delimitating the transition from "cold" to "hot" accretion. [Abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Simulations of protostellar collapse using multigroup radiation hydrodynamics. I. The first collapse

Radiative transfer plays a major role in the process of star formation. Many simulations of gravitational collapse of a cold gas cloud followed by the formation of a protostellar core use a grey treatment of radiative transfer coupled to the hydrodynamics. However, dust opacities which dominate extinction show large variations as a function of frequency. In this paper, we used frequency-dependent radiative transfer to investigate the influence of the opacity variations on the properties of Larson's first core. We used a multigroup M1 moment model in a 1D radiation hydrodynamics code to simulate the spherically symmetric collapse of a 1 solar mass cloud core. Monochromatic dust opacities for five different temperature ranges were used to compute Planck and Rosseland means inside each frequency group. The results are very consistent with previous studies and only small differences were observed between the grey and multigroup simulations. For a same central density, the multigroup simulations tend to produce first cores with a slightly higher radius and central temperature. We also performed simulations of the collapse of a 10 and 0.1 solar mass cloud, which showed the properties of the first core to be independent of the initial cloud mass, with again no major differences between grey and multigroup models. For Larson's first collapse, where temperatures remain below 2000 K, the vast majority of the radiation energy lies in the IR regime and the system is optically thick. In this regime, the grey approximation does a good job reproducing the correct opacities, as long as there are no large opacity variations on scales much smaller than the width of the Planck function. The multigroup method is however expected to yield more important differences in the later stages of the collapse when high energy (UV and X-ray) radiation is present and matter and radiation are strongly decoupled.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Analytical star formation rate from gravoturbulent fragmentation

We present an analytical determination of the star formation rate (SFR) in molecular clouds, based on a time-dependent extension of our analytical theory of the stellar initial mass function (IMF). The theory yields SFR's in good agreement with observations, suggesting that turbulence {\it is} the dominant, initial process responsible for star formation. In contrast to previous SFR theories, the present one does not invoke an ad-hoc density threshold for star formation; instead, the SFR {\it continuously} increases with gas density, naturally yielding two different characteristic regimes, thus two different slopes in the SFR vs gas density relationship, in agreement with observational determinations. Besides the complete SFR derivation, we also provide a simplified expression, which reproduces reasonably well the complete calculations and can easily be used for quick determinations of SFR's in cloud environments. A key property at the heart of both our complete and simplified theory is that the SFR involves a {\it density-dependent dynamical time}, characteristic of each collapsing (prestellar) overdense region in the cloud, instead of one single mean or critical freefall timescale. Unfortunately, the SFR also depends on some ill determined parameters, such as the core-to-star mass conversion efficiency and the crossing timescale. Although we provide estimates for these parameters, their uncertainty hampers a precise quantitative determination of the SFR, within less than a factor of a few.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dimensional argument for the impact of turbulent support on the stellar initial mass function

We present a simple dimensional argument to illustrate the impact of nonthermal support from turbulent velocity dispersion on the shape of the prestellar core mass function (CMF), precursor of the stellar initial mass function (IMF). The argument demonstrates the need to invoke such support to recover the Salpeter slope in the high-mass part of the CMF/IMF, whereas pure thermal support leads to a much steeper slope. This simple dimensional argument clearly highlights the results obtained in the complete Hennebelle-Chabrier theory of the IMF.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Distorted, non-spherical transiting planets: impact on the transit depth and on the radius determination

We quantify the systematic impact of the non-spherical shape of transiting planets and brown dwarfs, due to tidal forces and rotation, on the observed transit depth. Such a departure from sphericity leads to a bias in the derivation of the transit radius from the light curve and affects the comparison with planet structure and evolution models which assume spherical symmetry. As the tidally deformed planet projects its smallest cross section area during the transit, the measured effective radius is smaller than the one of the unperturbed spherical planet. This effect can be corrected by calculating the theoretical shape of the observed planet. We derive simple analytical expressions for the ellipsoidal shape of a fluid object (star or planet) accounting for both tidal and rotational deformations and calibratre it with fully numerical evolution models in the 0.3Mjup-75Mjup mass range. Our calculations yield a 20% effect on the transit depth, i.e. a 10% decrease of the measured radius, for the extreme case of a 1Mjup planet orbiting a Sun-like star at 0.01AU. For the closest planets detected so far (< 0.05 AU), the effect on the radius is of the order of 1 to 10%, by no means a negligible effect, enhancing the puzzling problem of the anomalously large bloated planets. These corrections must thus be taken into account for a correct determination of the radius from the transit light curve. Our analytical expressions can be easily used to calculate these corrections, due to the non-spherical shape of the planet, on the observed transit depth and thus to derive the planet's real equilibrium radius. They can also be used to model ellipsoidal variations of the stellar flux now detected in the CoRoT and Kepler light curves. We also derive directly usable analytical expressions for the moment of inertia, oblateness and Love number (k_2) of a fluid planet as a function of its mass.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Mach Number Dependence of Turbulent Magnetic Field Amplification: Solenoidal versus Compressive Flows

We study the growth rate and saturation level of the turbulent dynamo in magnetohydrodynamical simulations of turbulence, driven with solenoidal (divergence-free) or compressive (curl-free) forcing. For models with Mach numbers ranging from 0.02 to 20, we find significantly different magnetic field geometries, amplification rates, and saturation levels, decreasing strongly at the transition from subsonic to supersonic flows, due to the development of shocks. Both extreme types of turbulent forcing drive the dynamo, but solenoidal forcing is more efficient, because it produces more vorticity.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Physical and radiative properties of the first core accretion shock

Radiative shocks play a dominant role in star formation. The accretion shocks on the first and second Larson's cores involve radiative processes and are thus characteristic of radiative shocks. In this study, we explore the formation of the first Larson's core and characterize the radiative and dynamical properties of the accretion shock, using both analytical and numerical approaches. We develop both numerical RHD calculations and a semi-analytical model that characterize radiative shocks in various physical conditions, for radiating or barotropic fluids. Then, we perform 1D spherical collapse calculations of the first Larson's core, using a grey approximation for the opacity of the material. We consider three different models for radiative transfer, namely: the barotropic approximation, the FLD approximation and the more complete M1 model. We investigate the characteristic properties of the collapse and of the first core formation. Comparison between the numerical results and our semi-analytical model shows that this latter reproduces quite well the core properties obtained with the numerical calculations. The accretion shock on the first Larson core is found to be supercritical, implying that all the accretion shock energy on the core is radiated away. The FLD approximation is found to agree quite well with the results based on the M1 model, and is thus appropriate to study the star formation process. In contrast, the barotropic approximation does not correctly describe the thermal properties of the gas during the collapse. We show that a consistent treatment of radiation and hydrodynamics is mandatory to correctly handle the cooling of the gas during the core formation and thus to obtain the correct mechanical and thermal properties for this latter.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Radiation hydrodynamics with Adaptive Mesh Refinement and application to prestellar core collapse. I Methods

Radiative transfer has a strong impact on the collapse and the fragmentation of prestellar dense cores. We present the radiation-hydrodynamics solver we designed for the RAMSES code. The method is designed for astrophysical purposes, and in particular for protostellar collapse. We present the solver, using the co-moving frame to evaluate the radiative quantities. We use the popular flux limited diffusion approximation, under the grey approximation (one group of photon). The solver is based on the second-order Godunov scheme of RAMSES for its hyperbolic part, and on an implicit scheme for the radiation diffusion and the coupling between radiation and matter. We report in details our methodology to integrate the RHD solver into RAMSES. We test successfully the method against several conventional tests. For validation in 3D, we perform calculations of the collapse of an isolated 1 M_sun prestellar dense core, without rotation. We compare successfully the results with previous studies using different models for radiation and hydrodynamics. We have developed a full radiation hydrodynamics solver in the RAMSES code, that handles adaptive mesh refinement grids. The method is a combination of an explicit scheme and an implicit scheme, accurate to the second-order in space. Our method is well suited for star formation purposes. Results of multidimensional dense core collapse calculations with rotation are presented in a companion paper.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Theories of the initial mass function

We review the various theories which have been proposed along the years to explain the origin of the stellar initial mass function. We pay particular attention to four models, namely the competitive accretion and the theories based respectively on stopped accretion, MHD shocks and turbulent dispersion. In each case, we derive the main assumptions and calculations that support each theory and stress their respective successes and failures or difficulties.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Effect of episodic accretion on the structure and the lithium depletion of low-mass stars and planet-hosting stars

Following up our recent analysis devoted to the impact of non steady accretion on the location of young low-mass stars or brown dwarfs in the Herzsprung-Russell diagram, we perform a detailed analysis devoted to the effect of burst accretion on the internal structure of low-mass and solar type stars. We find that episodic accretion can produce objects with significantly higher central temperatures than the ones of the non accreting counterparts of same mass and age. As a consequence, lithium depletion can be severely enhanced in these objects. This provides a natural explanation for the unexpected level of lithium depletion observed in young objects for the inferred age of their parent cluster. These results confirm the limited reliability of lithium abundance as a criterion for assessing or rejecting cluster membership. They also show that lithium is not a reliable age indicator, because its fate strongly depends on the past accretion history of the star. Under the assumption that giant planets primarily form in massive disks prone to gravitational instability and thus to accretion burst episodes, the same analysis also explains the higher Li depletion observed in planet hosting stars. At last, we show that, depending on the burst rate and intensity, accretion outbursts can produce solar mass stars with lower convective envelope masses, at ages less than a few tens of Myr, than predicted by standard (non or slowly accreting) pre-main sequence models. This result has interesting, although speculative, implications for the recently discovered depletion of refractory elements in the Sun.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Is tidal heating sufficient to explain bloated exoplanets? Consistent calculations accounting for finite initial eccentricity

In this paper, we present the consistent evolution of short-period exoplanets coupling the tidal and gravothermal evolution of the planet. Contrarily to previous similar studies, our calculations are based on the complete tidal evolution equations of the Hut model, valid at any order in eccentricity, obliquity and spin. We demonstrate, both analytically and numerically, that, except if the system was formed with a nearly circular orbit (e<0.2), solving consistently the complete tidal equations is mandatory to derive correct tidal evolution histories. We show that calculations based on tidal models truncated at second order in eccentricity, as done in all previous studies, lead to erroneous tidal evolutions. As a consequence, tidal energy dissipation rates are severely underestimated in all these calculations and the characteristic timescales for the various orbital parameters evolutions can be wrong by up to three orders in magnitude. Based on these complete, consistent calculations, we revisit the viability of the tidal heating hypothesis to explain the anomalously large radius of transiting giant planets. We show that, even though tidal dissipation does provide a substantial contribution to the planet's heat budget and can explain some of the moderately bloated hot-Jupiters, this mechanism can not explain alone the properties of the most inflated objects, including HD 209458b. Indeed, solving the complete tidal equations shows that enhanced tidal dissipation and thus orbit circularization occur too early during the planet's evolution to provide enough extra energy at the present epoch. In that case another mechanisms, such as stellar irradiation induced surface winds dissipating in the planet's tidal bulges, or inefficient convection in the planet's interior must be invoked, together with tidal dissipation, to provide all the pieces of the abnormally large exoplanet puzzle.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Radiative, magnetic and numerical feedbacks on small-scale fragmentation

Radiative feedback and magnetic field are understood to have a strong impact on the protostellar collapse. We present high resolution numerical calculations of the collapse of a 1 solar mass dense core in solid body rotation, including both radiative transfer and magnetic field. Using typical parameters for low-mass cores, we study thoroughly the effect of radiative transfer and magnetic field on the first core formation and fragmentation. We show that including the two aforementioned physical processes does not correspond to the simple picture of adding them separately. The interplay between the two is extremely strong, via the magnetic braking and the radiation from the accretion shock.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Star formation: statistical measure of the correlation between the prestellar core mass function and the stellar initial mass function

We present a simple statistical analysis of recent numerical simulations exploring the correlation between the core mass function obtained from the fragmentation of a molecular cloud and the stellar mass function which forms from these collapsing cores. Our analysis shows that the distributions of bound cores and sink particles obtained in the simulations are consistent with the sinks being formed predominantly from their parent core mass reservoir, with a statistical dispersion of the order of one third of the core mass. Such a characteristic dispersion suggests that the stellar initial mass function is relatively tightly correlated to the parent core mass function, leading to two similar distributions, as observed. This in turn argues in favor of the IMF being essentially determined at the early stages of core formation and being only weakly affected by the various environmental factors beyond the initial core mass reservoir, at least in the mass range explored in the present study. Accordingly, the final IMF of a star forming region should be determined reasonably accurately, statistically speaking, from the initial core mass function, provided some uniform efficiency factor. The calculations also show that these statistical fluctuations, due e.g. to variations among the core properties, broaden the low-mass tail of the IMF compared with the parent CMF, providing an explanation for the fact that this latter appears to underestimate the number of "pre brown dwarf" cores compared with the observationally-derived brown dwarf IMF.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Lyot Project Direct Imaging Survey of Substellar Companions: Statistical Analysis and Information from Nondetections

The Lyot project used an optimized Lyot coronagraph with Extreme Adaptive Optics at the 3.63m Advanced Electro-Optical System telescope (AEOS) to observe 86 stars from 2004 to 2007. In this paper we give an overview of the survey results and a statistical analysis of the observed nondetections around 58 of our targets to place constraints on the population of substellar companions to nearby stars. The observations did not detect any companion in the substellar regime. Since null results can be as important as detections, we analyzed each observation to determine the characteristics of the companions that can be ruled out. For this purpose we use a Monte Carlo approach to produce artificial companions, and determine their detectability by comparison with the sensitivity curve for each star. All the non-detection results are combined using a Bayesian approach and we provide upper limits on the population of giant exoplanets and brown dwarfs for this sample of stars. Our nondetections confirm the rarity of brown dwarfs around solar-like stars and we constrain the frequency of massive substellar companions (M>40Mjup) at orbital separation between and 10 and 50 AU to be <20%.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The radius anomaly in the planet/brown dwarf overlapping mass regime

The recent detection of the transit of very massive substellar companions (CoRoT-3b, Deleuil et al. 2008; CoRoT-15b, Bouchy et al. 2010; WASP-30b, Anderson et al. 2010; Hat-P-20b, Bakos et al. 2010) provides a strong constraint to planet and brown dwarf formation and migration mechanisms. Whether these objects are brown dwarfs originating from the gravitational collapse of a dense molecular cloud that, at the same time, gave birth to the more massive stellar companion, or whether they are planets that formed through core accretion of solids in the protoplanetary disk can not always been determined unambiguously and the mechanisms responsible for their short orbital distances are not yet fully understood. In this contribution, we examine the possibility to constrain the nature of a massive substellar object from the various observables provided by the combination of Radial Velocity and Photometry measurements (e.g. M_p, R_p, M_s, Age, a, e...). In a second part, developments in the modeling of tidal evolution at high eccentricity and inclination - as measured for HD 80 606 with e=0.9337 (Naef et al. 2001), XO-3 with a stellar obliquity >37.3+-3.7 deg (Hébrard et al. 2008; Winn et al. 2009) and several other exoplanets - are discussed along with their implication in the understanding of the radius anomaly problem of extrasolar giant planets.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Analytical theory for the initial mass function: II. Properties of the flow

Recently, Hennebelle and Chabrier (2008) derived an analytical theory for the mass spectrum of non self-gravitating clumps associated with overdensities in molecular clouds and for the initial mass function of gravitationally bound prestellar cores, as produced by the turbulent collapse of the cloud. In this companion paper, we examine the effects of the non-isothermality of the flow, of the turbulence forcing and of local fluctuation of the velocity dispersion, on the mass function. In particular, we investigate the influence of a polytropic equation of state and of the effective adiabatic exponent $γ$ and find that it has a drastic influence on the low mass part of the IMF. We also consider a barotropic equation of state (i.e. a piecewise polytropic eos) that mimics the thermal behaviour of the molecular gas and compare the prediction of our theory with the results of numerical simulations and with the observationally-derived IMF, for cloud parameters which satisfy Larson's type relations. We find that for clouds whose density is, at all scales, almost an order of magnitude larger than the density inferred for the CO clumps in the Galaxy, a good agreement is obtained between the theory and the observed IMF, suggesting that star formation preferentially occurs in high density environments. We derive an analytical expression for the IMF which generalizes the expression previously obtained for the isothermal case. This easy-to-implement analytical IMF should serve as a template to compare observational or numerical results with the theory.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Plasma physics and planetary astrophysics

In this review, I briefly summarize the present status of experimental and theoretical investigations of the properties of matter under conditions characteristic of planetary interiors, from terrestrial to jovian planets. I first focus on the two lightest elements, hydrogen and helium, and discuss recent theoretical and experimental investigations of their properties at high pressure and temperature. Then, I discuss the impact of these properties, as well as of the equation of state of heavier elements, on planetary interiors. Finally, I highlight the importance of exoplanet transit observations and of the inferred mass-radius relationships to determine the planetary interior compositions.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Large-scale alpha^2-dynamo in low-mass stars and brown dwarfs

We develop a model based on 3D mean-field MHD for the generation of large scale magnetic fields in fully convective objects like low-mass stars, brown dwarfs and possibly gaseous planets. The dynamo process is of alpha^2 type and thus differs from the shell-dynamo at work in more massive stars.The alpha^2 dynamo becomes supercritical for Rossby numbers Ro\la 10. It generates a large-scale, non-axisymmetric, steady field that is symmetric with respect to the equatorial plane. Saturation of the alpha^2-generated field at the equipartition yields strengths of several kGauss, in agreement with observations of active M dwarfs, and provides a qualitative explanation for the observed activity saturation in late M stars. For brown dwarfs with a conductive core, as occurs at the center of the most massive and oldest of these objects, we have also studied an alpha^2-Omega dynamo, i.e. the effect of differential rotation. In this case the field is predominantly toroidal, axisymmetric and oscillatory, like the solar field. The topology of the field in the fully convective objects exhibits a high order multipole character that differs from the aligned dipole field generated by the alpha-Omega dynamo. The strong reduction of the dipolar component due to the field non-axisymmetry should considerably reduce the Alfven radius and thus the efficiency of magnetic braking, providing an appealing explanation for the decreasing angular momentum loss rate observed in low-mass stars and brown dwarfs. This may have also implications for cataclysmic variables below the period gap. In spite of this large-scale field, the decreasing conductivity in the dominantly neutral atmosphere of these objects prevents the current generation necessary to support a chromosphere and thus activity. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

The initial mass function : from Salpeter 1955 to 2005

Fifty years after Ed Salpeter's seminal paper, tremendous progress both on the observational and theoretical sides allow a fairly accurate determination of the Galactic IMF not only down to the hydrogen-burning limit but into the brown dwarf domain. The present review includes the most recent observations of low-mass stars and brown dwarfs to determine this IMF and the related Galactic mass budget. The IMF definitely exhibits a similar behaviour in various environments, disk, young and globular clusters, spheroid. Small scale dissipation of large scale compressible MHD turbulence seems to be the underlying triggering mechanism for star formation. Modern simulations of compressible MHD turbulence yield an IMF consistent with the one derived from observations.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

The Galactic disk mass function: reconciliation of the HST and nearby determinations

We derive and parametrize the Galactic mass function (MF) below 1 $\msol$ characteristic of both single objects and binary systems. We resolve the long standing discrepancy between the MFs derived from the HST and from the nearby luminosity functions, respectively. We show that this discrepancy stemmed from {\it two} cumulative effects, namely (i) incorrect color-magnitude determined distances, due a substantial fraction of M dwarfs in the HST sample belonging to the metal-depleted, thick-disk population, as corrected recently by Zheng et al. and (ii) unresolved binaries. We show that both the nearby and HST MF for unresolved systems are consistent with a fraction $\sim $50% of M-dwarf binaries, with the mass of both the primaries and the companions originating from the same underlying single MF. This implies that $\sim$30% of M dwarfs should have an M dwarf companion and $\sim$20% should have a brown dwarf companion, in agreement with recent determinations. The present calculations show that the so-called "brown-dwarf desert" should be reinterpreted as a lack of high mass-ratio $(m_2/m_1\la 0.1$) systems, and does not preclude a substantial fraction of brown dwarfs as companions of M dwarfs or for other brown dwarfs.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Determination of the globular cluster and halo stellar mass functions and stellar and brown dwarf densities

We use recent low-mass star models, which reproduce accurately the observed sequences of various globular clusters, to convert the observed luminosity functions into bolometric luminosity functions and mass functions down to the bottom of the main sequence. These mass functions are well describedby a slowly rising power-law $dN/dm\propto m^{-α}$, with $0.5\wig < α\wig < 1.5$, down to $\sim 0.1 \msol$, independently of the metallicity, suggesting a rather universal behaviour of the cluster initial mass functions. We predict luminosity functions in the NICMOS filters in the stellar and in the brown dwarf domains for different mass functions and metallicities. We apply these calculations to the determination, slope and normalization, of the mass function of the Galactic halo (spheroid and dark halo). The spheroid mass function is well described by the afore-mentioned power-law function with $α\sim 1.7\pm 0.2$ down to 0.1 $\msol$, although a slowly decreasing mass function below $\sim 0.15 \msol$ can not be excluded with the data presently available. Comparison with the Hubble Deep Field star counts is consistent with such a mass function and excludes a significant stellar population in the dark halo. Consistent analysis with recent microlensing experiments towards the LMC shows that the spheroid and the dark-halo stellar+brown dwarf populations represent at most $\sim$ 1% of the Galactic dark matter density. This clearly excludes brown dwarfs and low-mass stars as significant dark matter candidates.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Evolutionary models for metal-poor low-mass stars. Lower main sequence of globular clusters and halo field stars

We have performed evolutionary calculations of very-low-mass stars from 0.08 to 0.8 $\msol$ for different metallicites from [M/H]= -2.0 to -1.0 and we have tabulated the mechanical, thermal and photometric characteristics of these models. The calculations include the most recent interior physics and improved non-grey atmosphere models. The models reproduce the entire main sequences of the globular clusters observed with the Hubble Space Telescope over the afore-mentioned range of metallicity. Comparisons are made in the WFPC2 Flight system including the F555, F606 and F814 filters, and in the standard Johnson-Cousins system. We examine the effects of different physical parameters, mixing-length, $α$-enriched elements, helium fraction, as well as the accuracy of the photometric transformations of the HST data into standard systems. We derive mass-effective temperature and mass-magnitude relationships and we compare the results with the ones obtained with different grey-like approximations. These latter are shown to yield inaccurate relations, in particular near the hydrogen-burning limit. We derive new hydrogen-burning minimum masses, and the corresponding absolute magnitudes, for the different metallicities. We predict color-magnitude diagrams in the infrared NICMOS filters, to be used for the next generation of the HST observations, providing mass-magnitudes relationships in these colors down to the brown-dwarf limit. We show that the expected signature of the stellar to substellar transition in color-magnitude diagrams is a severe blueshift in the infrared colors, due to the increasing collision-induced absorption of molecular hydrogen with increasing density and decreasing temperature.

0 citations0 reviewsNo code linkedOpen access