Source author record

G. Chabrier

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

36works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

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

36 published item(s)

preprint2021arXiv

Non-ideal self-gravity and cosmology: the importance of correlations in the dynamics of the large-scale structures of the Universe

Inspired by the statistical mechanics of an ensemble of interacting particles (BBGKY hierarchy), we propose to account for small-scale inhomogeneities in self-gravitating astrophysical fluids by deriving a non-ideal Virial theorem and non-ideal Navier-Stokes equations using a decomposition of the gravitational force into a near- and far-field component. These equations involve the pair radial distribution function (similar to the two-point correlation function), similarly to the interaction energy and equation of state in liquids. Small-scale correlations lead to a non-ideal amplification of the gravitational interaction energy, whose omission leads to a missing mass problem, e.g., in galaxies and galaxy clusters. We also propose an extension of the Friedmann equations in the non-ideal regime. We estimate the non-ideal amplification factor of the gravitational interaction energy of the baryons to lie between 5 and 20, potentially explaining the observed value of the Hubble parameter. Within this framework, the acceleration of the expansion emerges naturally because of the increasing number of sub-structures induced by gravitational collapse, which increases their contribution to the total gravitational energy. A simple estimate predicts a non-ideal deceleration parameter qni~-1;this is potentially the first determination of the observed value based on an intuitively physical argument. We suggest that correlations and gravitational interactions could produce a transition to a viscous regime that can lead to flat rotation curves. This transition could also explain the dichotomy between (Keplerian) LSB elliptical galaxy and (non-Keplerian) spiral galaxy rotation profiles. Overall, our results demonstrate that non-ideal effects induced by inhomogeneities must be taken into account in order to properly determine the gravitational dynamics of galaxies and the larger scale universe.

preprint2020arXiv

Crust structure and thermal evolution of neutron stars in soft X-ray transients

We study the effect of physics input on thermal evolution of neutron stars in soft X-ray transients (SXTs). In particular, we consider different modern models of the sources of deep crustal heating during accretion episodes and the effects brought about by impurities embedded in the crust during its formation. We simulate thermal structure and evolution of episodically accreting neutron stars under different assumptions on the crust composition and on the distribution of heat sources and impurities. For the nonaccreted crust, we consider the nuclear charge fluctuations that arise at crust formation. For the accreted crust, we compare different theoretical models of composition and internal heating. We also compare results of numerical simulations with observations of the crust cooling in SXT MXB 1659-29. We found that the nonaccreted part of the inner crust of a neutron star can have a layered structure, with almost pure crystalline layers interchanging with layers composed of mixtures of different nuclei. The latter layers have relatively low thermal conductivities, which affects thermal evolution of the transients. The impurity distribution in the crust strongly depends on the models of the dense matter and the crust formation scenario. The shallow heating that is needed to reach agreement between the theory and observations depends on characteristics of the crust and envelope.

preprint2016arXiv

Cloudless atmospheres for L/T dwarfs and extra-solar giant planets

The admitted, conventional scenario to explain the complex spectral evolution of brown dwarfs (BD) since their first detections twenty years ago, has always been the key role played by micron-size condensates, called "dust" or "clouds", in their atmosphere. This scenario, however, faces major problems, in particular the J-band brightening and the resurgence of FeH absorption at the L to T transition, and a physical first-principle understanding of this transition is lacking. In this paper, we propose a new, completely different explanation for BD and extrasolar giant planet (EGP) spectral evolution, without the need to invoke clouds. We show that, due to the slowness of the CO/CH4 and N2/NH3 chemical reactions, brown dwarf (L and T, respectively) and EGP atmospheres are subject to a thermo-chemical instability similar in nature to the fingering or chemical convective instability present in Earth oceans and at the Earth core/mantle boundary. The induced small-scale turbulent energy transport reduces the temperature gradient in the atmosphere, explaining the observed increase in near infrared J-H and J-K colors of L dwarfs and hot EGPs, while a warming up of the deep atmosphere along the L to T transition, as the CO/CH4 instability vanishes, naturally solves the two aforementioned puzzles, and provides a physical explanation of the L to T transition. This new picture leads to a drastic revision of our understanding of BD and EGP atmospheres and their evolution.

preprint2016arXiv

Near-Infrared Spectroscopy of the Y0 WISEP J173835.52+273258.9 and the Y1 WISE J035000.32-565830.2: the Importance of Non-Equilibrium Chemistry

We present new near-infrared spectra, obtained at Gemini Observatory, for two Y dwarfs: WISE J035000.32-565830.2 (W0350) and WISEP J173835.52+273258.9 (W1738). A FLAMINGOS-2 R=540 spectrum was obtained for W0350, covering 1.0 < lambda um < 1.7, and a cross-dispersed GNIRS R=2800 spectrum was obtained for W1738, covering 0.993-1.087 um, 1.191-1.305 um, 1.589-1.631 um, and 1.985-2.175 um, in four orders. We also present revised YJH photometry for W1738, using new NIRI Y and J imaging, and a re-analysis of the previously published NIRI H band images. We compare these data, together with previously published data for late-T and Y dwarfs, to cloud-free models of solar metallicity, calculated both in chemical equilibrium and with disequilibrium driven by vertical transport. We find that for the Y dwarfs the non-equilibrium models reproduce the near-infrared data better than the equilibrium models. The remaining discrepancies suggest that fine-tuning the CH_4/CO and NH_3/N_2 balance is needed. Improved trigonometric parallaxes would improve the analysis. Despite the uncertainties and discrepancies, the models reproduce the observed near-infrared spectra well. We find that for the Y0, W1738, T_eff = 425 +/- 25 K and log g = 4.0 +/- 0.25, and for the Y1, W0350, T_eff = 350 +/- 25 K and log g = 4.0 +/- 0.25. W1738 may be metal-rich. Based on evolutionary models, these temperatures and gravities correspond to a mass range for both Y dwarfs of 3-9 Jupiter masses, with W0350 being a cooler, slightly older, version of W1738; the age of W0350 is 0.3-3 Gyr, and the age of W1738 is 0.15-1 Gyr.

preprint2016arXiv

Self-consistent evolution of accreting low-mass stars and brown dwarfs

We present self-consistent calculations coupling numerical hydrodynamics simulations of collapsing pre-stellar cores and stellar evolution models of accreting objects. We analyse the main impact of consistent accretion history on the evolution and lithium depletion of young low-mass stars and brown dwarfs. These consistent models confirm the generation of a luminosity spread in Herzsprung-Russell diagrams at ages $\sim$ 1-10 Myr. They also confirm that early accretion can produce objects with abnormal Li depletion, as found in a previous study that was based on arbitrary accretion rates. The results strengthen that objects with anomalously high level of Li depletion in young clusters should be extremely rare. We also find that early phases of burst accretion can produce coeval models of similar mass with a range of different Li surface abundances, and in particular with Li-excess compared to the predictions of non-accreting counterparts. This result is due to a subtle competition between the effect of burst accretion and its impact on the central stellar temperature, the growth of the stellar radiative core and the accretion of fresh Li from the accretion disk. Only consistent models could reveal such a subtle combination of effects. This new result could explain the recent, puzzling observations of Li-excess of fast rotators in the young cluster NGC 2264. Present self-consistent accreting models are available in electronic form.

preprint2015arXiv

Fingering convection and cloudless models for cool brown dwarf atmospheres

This work aims to improve the current understanding of the atmospheres of brown dwarfs, especially cold ones with spectral type T and Y, whose modeling is a current challenge. Silicate and iron clouds are believed to disappear at the photosphere at the L/T transition, but cloudless models fail to reproduce correctly the spectra of T dwarfs, advocating for the addition of more physics, e.g. other types of clouds or internal energy transport mechanisms. We use a one-dimensional (1D) radiative/convective equilibrium code ATMO to investigate this issue. This code includes both equilibrium and out-of-equilibrium chemistry and solves consistently the PT structure. Included opacity sources are H2-H2, H2-He, H2O, CO, CO2, CH4, NH3, K, Na, and TiO, VO if they are present in the atmosphere. We show that the spectra of Y dwarfs can be accurately reproduced with a cloudless model if vertical mixing and NH3 quenching are taken into account. T dwarf spectra still have some reddening in e.g. J - H compared to cloudless models. This reddening can be reproduced by slightly reducing the temperature gradient in the atmosphere. We propose that this reduction of the stabilizing temperature gradient in these layers, leading to cooler structures, is due to the onset of fingering convection, triggered by the destabilizing impact of condensation of very thin dust.

preprint2015arXiv

Incorporating Ambipolar and Ohmic Diffusion in the AMR MHD code RAMSES

We have implemented non-ideal Magneto-Hydrodynamics (MHD) effects in the Adaptive Mesh Refinement (AMR) code RAMSES, namely ambipolar diffusion and Ohmic dissipation, as additional source terms in the ideal MHD equations. We describe in details how we have discretized these terms using the adaptive Cartesian mesh, and how the time step is diminished with respect to the ideal case, in order to perform a stable time integration. We have performed a large suite of test runs, featuring the Barenblatt diffusion test, the Ohmic diffusion test, the C-shock test and the Alfven wave test. For the latter, we have performed a careful truncation error analysis to estimate the magnitude of the numerical diffusion induced by our Godunov scheme, allowing us to estimate the spatial resolution that is required to address non-ideal MHD effects reliably. We show that our scheme is second-order accurate, and is therefore ideally suited to study non-ideal MHD effects in the context of star formation and molecular cloud dynamics.

preprint2015arXiv

New evolutionary models for pre-main sequence and main sequence low-mass stars down to the hydrogen-burning limit

We present new models for low-mass stars down to the hydrogen-burning limit that consistently couple atmosphere and interior structures, thereby superseding the widely used BCAH98 models. The new models include updated molecular linelists and solar abundances, as well as atmospheric convection parameters calibrated on 2D/3D radiative hydrodynamics simulations. Comparison of these models with observations in various colour-magnitude diagrams for various ages shows significant improvement over previous generations of models. The new models can solve flaws that are present in the previous ones, such as the prediction of optical colours that are too blue compared to M dwarf observations. They can also reproduce the four components of the young quadruple system LkCa 3 in a colour-magnitude diagram with one single isochrone, in contrast to any presently existing model. In this paper we also highlight the need for consistency when comparing models and observations, with the necessity of using evolutionary models and colours based on the same atmospheric structures.

preprint2014arXiv

Giant planet and brown dwarf formation

Understanding the dominant brown dwarf and giant planet formation processes, and finding out whether these processes rely on completely different mechanisms or share common channels represents one of the major challenges of astronomy and remains the subject of heated debates. It is the aim of this review to summarize the latest developments in this field and to address the issue of origin by confronting different brown dwarf and giant planet formation scenarios to presently available observational constraints. As examined in the review, if objects are classified as "Brown Dwarfs" or "Giant Planets" on the basis of their formation mechanism, it has now become clear that their mass domains overlap and that there is no mass limit between these two distinct populations. Furthermore, while there is increasing observational evidence for the existence of non-deuterium burning brown dwarfs, some giant planets, characterized by a significantly metal enriched composition, might be massive enough to ignite deuterium burning in their core. Deuterium burning (or lack of) thus plays no role in either brown dwarf or giant planet formation. Consequently, we argue that the IAU definition to distinguish these two populations has no physical justification and brings scientific confusion. In contrast, brown dwarfs and giant planets might bear some imprints of their formation mechanism, notably in their mean density and in the physical properties of their atmosphere. Future direct imaging surveys will undoubtedly provide crucial information and perhaps provide some clear observational diagnostics to unambiguously distinguish these different astrophysical objects.

preprint2014arXiv

Opacities and spectra of hydrogen atmospheres of moderately magnetized neutron stars

There is observational evidence that central compact objects (CCOs) in supernova remnants have moderately strong magnetic fields $B\sim10^{11}$ G. Meanwhile, available models of partially ionized hydrogen atmospheres of neutron stars with strong magnetic fields are restricted to $B\gtrsim10^{12}$ G. We extend the equation of state and radiative opacities, presented in previous papers for $10^{12}\mbox{ G}\lesssim B \lesssim 10^{15}$ G, to weaker fields. An equation of state and radiative opacities for a partially ionized hydrogen plasma are obtained at magnetic fields $B$, temperatures $T$, and densities $ρ$ typical for atmospheres of CCOs and other isolated neutron stars with moderately strong magnetic fields. The first- and second-order thermodynamic functions, monochromatic radiative opacities, and Rosseland mean opacities are calculated and tabulated, taking account of partial ionization, for $3\times10^{10}\mbox{ G}\lesssim B\lesssim 10^{12}$ G, $10^5$ K $\lesssim T\lesssim 10^7$ K, and a wide range of densities. Atmosphere models and spectra are calculated to verify the applicability of the results and to determine the range of magnetic fields and effective temperatures where the incomplete ionization of the hydrogen plasma is important.

preprint2014arXiv

Planetary internal structures

This chapter reviews the most recent advancements on the topic of terrestrial and giant planet interiors, including Solar System and extrasolar objects. Starting from an observed mass-radius diagram for known planets in the Universe, we will discuss the various types of planets appearing in this diagram and describe internal structures for each type. The review will summarize the status of theoretical and experimental works performed in the field of equation of states (EOS) for materials relevant to planetary interiors and will address the main theoretical and experimental uncertainties and challenges. It will discuss the impact of new EOS on interior structures and bulk composition determination. We will discuss important dynamical processes which strongly impact the interior and evolutionary properties of planets (e.g plate tectonics, semiconvection) and describe non standard models recently suggested for our giant planets. We will address the case of short-period, strongly irradiated exoplanets and critically analyse some of the physical mechanisms which have been suggested to explain their anomalously large radius.

preprint2014arXiv

Variations of the stellar initial mass function in the progenitors of massive early-type galaxies and in extreme starburst environments

We examine variations of the stellar initial mass function (IMF) in extreme environments within the formalism derived by Hennebelle \& Chabrier. We focus on conditions encountered in progenitors of massive early type galaxies and starburst regions. We show that, when applying the concept of turbulent Jeans mass as the characteristic mass for fragmentation in a turbulent medium, instead of the standard thermal Jeans mass for purely gravitational fragmentation, the peak of the IMF in such environments is shifted towards smaller masses, leading to a bottom-heavy IMF, as suggested by various observations. In very dense and turbulent environments, we predict that the high-mass tail of the IMF can become even steeper than the standard Salpeter IMF, with a limit for the power law exponent $α\simeq -2.7$, in agreement with recent observational determinations. This steepening is a direct consequence of the high densities and Mach values in such regions but also of the time dependence of the fragmentation process, as incorporated in the Hennebelle-Chabrier theory. We provide analytical parametrizations of these IMFs in such environments, to be used in galaxy evolution calculations. We also calculate the star formation rates and the mass-to-light ratios expected under such extreme conditions and show that they agree well with the values inferred in starburst environments and massive high-redshift galaxies. This reinforces the paradigm of star formation as being a universal process, i.e. the direct outcome of gravitationally unstable fluctuations in a density field initially generated by large scale shock-dominated turbulence. This globally enables us to infer the variations of the stellar IMF and related properties for atypical galactic conditions.

preprint2013arXiv

Electron screening effect on stellar thermonuclear fusion

We study the impact of plasma correlation effects on nonresonant thermonuclear reactions for various stellar objects, namely in the liquid envelopes of neutron stars, and the interiors of white dwarfs, low-mass stars, and substellar objects. We examine in particular the effect of electron screening on the enhancement of thermonuclear reactions in dense plasmas within and beyond the linear mixing rule approximation as well as the corrections due to quantum effects at high density. In addition, we examine some recent unconventional (Yukawa-potential and "quantum-tail") theoretical results on stellar thermonuclear fusions and show that these scenarios do not apply to stellar conditions.

preprint2013arXiv

Simulations of protostellar collapse using multigroup radiation hydrodynamics. II. The second collapse

Star formation begins with the gravitational collapse of a dense core inside a molecular cloud. As the collapse progresses, the centre of the core begins to heat up as it becomes optically thick. The temperature and density in the centre eventually reach high enough values where fusion reactions can ignite; the protostar is born. This sequence of events entail many physical processes, of which radiative transfer is of paramount importance. Many simulations of protostellar collapse make use of a grey treatment of radiative transfer coupled to the hydrodynamics. However, interstellar gas and dust opacities present large variations as a function of frequency. In this paper, we follow-up on a previous paper on the collapse and formation of Larson's first core using multigroup radiation hydrodynamics (Paper I) by extending the calculations to the second phase of the collapse and the formation of Larson's second core. We have made the use of a non-ideal gas equation of state as well as an extensive set of spectral opacities in a spherically symmetric fully implicit Godunov code to model all the phases of the collapse of a 0.1, 1 and 10 solar mass cloud cores. We find that, for a same central density, there are only small differences between the grey and multigroup simulations. The first core accretion shock remains supercritical while the shock at the second core border is found to be strongly subcritical with all the accreted energy being transfered to the core. The size of the first core was found to vary somewhat in the different simulations (more unstable clouds form smaller first cores) while the size, mass and temperature of the second cores are independent of initial cloud mass, size and temperature. Our simulations support the idea of a standard (universal) initial second core size of 0.003 AU and mass 0.0014 solar masses.

preprint2013arXiv

The effect of episodic accretion on the phase transition of CO and CO_2 in low-mass star formation

We study the evaporation and condensation of CO and CO_2 during the embedded stages of low-mass star formation by using numerical simulations. We focus on the effect of luminosity bursts, similar in magnitude to FUors and EXors, on the gas-phase abundance of CO and CO_2 in the protostellar disk and infalling envelope. The evolution of a young protostar and its environment is followed based on hydrodynamical models using the thin-disk approximation, coupled with a stellar evolution code and phase transformations of CO and CO_2. The accretion and associated luminosity bursts in our model are caused by disk gravitational fragmentation followed by quick migration of the fragments onto the forming protostar. We found that bursts with luminosity on the order of 100-200 L_sun can evaporate CO ices in part of the envelope. The typical freeze-out time of the gas-phase CO onto dust grains in the envelope (a few kyr) is much longer than the burst duration (100-200 yr). This results in an increased abundance of the gas-phase CO in the envelope long after the system has returned into a quiescent stage. In contrast, luminosity bursts can evaporate CO_2 ices only in the disk, where the freeze-out time of the gas-phase CO_2 is comparable to the burst duration. We thus confirm that luminosity bursts can leave long-lasting traces in the abundance of gas-phase CO in the infalling envelope, enabling the detection of recent bursts as suggested by previous semi-analytical studies.

preprint2013arXiv

The influence of frequency-dependent radiative transfer on the structures of radiative shocks

Radiative shocks are shocks in a gas where the radiative energy and flux coming from the very hot post-shock material are non-negligible in the shock's total energy budget, and are often large enough to heat the material ahead of the shock. Many simulations of radiative shocks, both in the contexts of astrophysics and laboratory experiments, use a grey treatment of radiative transfer coupled to the hydrodynamics. However, the opacities of the gas show large variations as a function of frequency and this needs to be taken into account if one wishes to reproduce the relevant physics. We have performed radiation hydrodynamics simulations of radiative shocks in Ar using multigroup (frequency dependent) radiative transfer with the HERACLES code. The opacities were taken from the ODALISC database. We show the influence of the number of frequency groups used on the dynamics and morphologies of subcritical and supercritical radiative shocks in Ar gas, and in particular on the extent of the radiative precursor. We find that simulations with even a low number of groups show significant differences compared to single-group (grey) simulations, and that in order to correctly model such shocks, a minimum number of groups is required. Results appear to eventually converge as the number of groups increases above 50. We were also able to resolve in our simulations of supercritical shocks the adaptation zones which connect the cooling layer to the final post-shock state and the precursor. Inside these adaptation zones, we find that the radiative flux just ahead of the shock in one or several high-opacity groups can heat the gas to a temperature higher than the post-shock temperature. Through the use of Hugoniot curves, we have checked the consistency of our radiation hydrodynamics scheme by showing that conservation of mass, momentum and energy (including radiative flux) holds. ABRIDGED.

preprint2013arXiv

The structure of warm dense matter modeled with an average atom model with ion-ion correlations

We present a new model of warm dense matter that represents an intermediate approach between the relative simplicity of ''one-ion'' average atom models and the more realistic but computationally expensive ab initio simulation methods. Physical realism is achieved primarily by including the correlations in the plasma that surrounds a central ion. The plasma is described with the Ornstein-Zernike integral equations theory of fluids, which is coupled to an average atom model for the central ion. In this contribution we emphasize the key elements and approximations and how they relate to and expand upon a typical average atom model. Besides being relatively inexpensive computationally, this approach offers several advantages over ab initio simulations but also has a number of limitations. The model is validated by comparisons with numerical solutions for the pair distribution function of the ions from ab initio simulations for several elements and a wide range of plasma conditions. Simulations results are reproduced remarkably well and simpler limiting theories are recovered as well. This model has many potential applications to calculation of the properties of warm dense matter such as the equation of state and conductivities for a wide range of temperatures and densities.

preprint2012arXiv

Thermonuclear fusion in dense stars: Electron screening, conductive cooling, and magnetic field effects

We study the plasma correlation effects on nonresonant thermonuclear reactions of carbon and oxygen in the interiors of white dwarfs and liquid envelopes of neutron stars. We examine the effects of electron screening on thermodynamic enhancement of thermonuclear reactions in dense plasmas beyond the linear mixing rule. Using these improved enhancement factors, we calculate carbon and oxygen ignition curves in white dwarfs and neutron stars. The energy balance and ignition conditions in neutron star envelopes are evaluated, taking their detailed thermal structure into account. The result is compared to the simplified "one-zone model," which is routinely used in the literature. We also consider the effect of strong magnetic fields on the ignition curves in the ocean of magnetars.

preprint2011arXiv

Equation of state for partially ionized carbon and oxygen mixtures at high temperatures

The equation of state (EOS) for partially ionized carbon, oxygen, and carbon-oxygen mixtures at temperatures 3\times10^5 K <~ T <~ 3\times10^6 K is calculated over a wide range of densities, using the method of free energy minimization in the framework of the chemical picture of plasmas. The free energy model is an improved extension of our model previously developed for pure carbon (Phys. Rev. E, 72, 046402; arXiv:physics/0510006). The internal partition functions of bound species are calculated by a self-consistent treatment of each ionization stage in the plasma environment taking into account pressure ionization. The long-range Coulomb interactions between ions and screening of the ions by free electrons are included using our previously published analytical model, recently improved, in particular for the case of mixtures. We also propose a simple but accurate method of calculation of the EOS of partially ionized binary mixtures based on detailed ionization balance calculations for pure substances.

preprint2010arXiv

The effect of magnetic activity on low-mass stars in eclipsing binaries

In recent years, analyses of eclipsing binary systems have unveiled differences between the observed fundamental properties of low-mass stars and those predicted by stellar structure models. Particularly, radius and effective temperatures computed from models are ~ 5-10% lower and ~ 3-5% higher than observed, respectively. These discrepancies have been attributed to different factors, notably to the high levels of magnetic activity present on these stars. In this paper, we test the effect of magnetic activity both on models and on the observational analysis of eclipsing binaries using a sample of such systems with accurate fundamental properties. Regarding stellar models, we have found that unrealistically high spot coverages need to be assumed to reproduce the observations. Tests considering metallicity effects and missing opacities on models indicate that these are not able to explain the radius discrepancies observed. With respect to the observations, we have tested the effect of several spot distributions on the light curve analysis. Our results show that spots cause systematic deviations on the stellar radii derived from light curve analysis when distributed mainly over the stellar poles. Assuming the existence of polar spots, overall agreement between models and observations is reached when ~ 35% spot coverage is considered on stellar models. Such spot coverage induces a systematic deviation in the radius determination from the light curve analysis of ~ 3% and is also compatible with the modulations observed on the light curves of these systems. Finally, we have found that the effect of activity or rotation on convective transport in partially radiative stars may also contribute to explain the differences seen in some of the systems with shorter orbital periods.

preprint2010arXiv

The physical properties of extrasolar planets

Tremendous progress in the science of extrasolar planets has been achieved since the discovery of a Jupiter orbiting the nearby Sun-like star 51 Pegasi in 1995. Theoretical models have now reached enough maturity to predict the characteristic properties of these new worlds, mass, radius, atmospheric signatures, and can be confronted with available observations. We review our current knowledge of the physical properties of exoplanets, internal structure and composition, atmospheric signatures, including expected biosignatures for exo-Earth planets, evolution, and the impact of tidal interaction and stellar irradiation on these properties for the short-period planets. We discuss the most recent theoretical achievements in the field and the still pending questions. We critically analyse the different solutions suggested to explain abnormally large radii of a significant fraction of transiting exoplanets. Special attention is devoted to the recently discovered transiting objects in the overlapping mass range between massive planets and low-mass brown dwarfs, stressing the ambiguous nature of these bodies, and we discuss the possible observable diagnostics to identify these two distinct populations. We also review our present understanding of planet formation and critically examine the different suggested formation mechanisms. We expect the present review to provide the basic theoretical background to capture the essential of the physics of exoplanet formation, structure and evolution, and the related observable signatures.

preprint2010arXiv

Thermodynamic functions of dense plasmas: analytic approximations for astrophysical applications

We briefly review analytic approximations of thermodynamic functions of fully ionized nonideal electron-ion plasmas, applicable in a wide range of plasma parameters, including the domains of nondegenerate and degenerate, nonrelativistic and relativistic electrons, weakly and strongly coupled Coulomb liquids, classical and quantum Coulomb crystals. We present improvements to previously published approximations. Our code for calculation of thermodynamic functions based on the reviewed approximations is made publicly available.

preprint2010arXiv

Understanding exoplanet formation, structure and evolution in 2010

In this short review, we summarize our present understanding (and non-understanding) of exoplanet formation, structure and evolution, in the light of the most recent discoveries. Recent observations of transiting massive brown dwarfs seem to remarkably confirm the predicted theoretical mass-radius relationship in this domain. This mass-radius relationship provides, in some cases, a powerful diagnostic to distinguish planets from brown dwarfs of same mass, as for instance for Hat-P-20b. If confirmed, this latter observation shows that planet formation takes place up to at least 8 Jupiter masses. Conversely, observations of brown dwarfs down to a few Jupiter masses in young, low-extinction clusters strongly suggest an overlapping mass domain between (massive) planets and (low-mass) brown dwarfs, i.e. no mass edge between these two distinct (in terms of formation mechanism) populations. At last, the large fraction of heavy material inferred for many of the transiting planets confirms the core-accretion scenario as been the dominant one for planet formation.

preprint2009arXiv

Addendum to "Equation of state of classical Coulomb plasma mixtures"

Recently developed analytic approximation for the equation of state of fully ionized nonideal electron-ion plasma mixtures [Potekhin et al., Phys. Rev. E, 79, 016411 (2009); arXiv:0812.4344], which covers the transition between the weak and strong Coulomb coupling regimes and reproduces numerical results obtained in the hypernetted chain (HNC) approximation, is modified in order to fit the small deviations from the linear mixing in the strong coupling regime, revealed by recent Monte Carlo simulations. In addition, a mixing rule is proposed for the regime of weak coupling, which generalizes post-Debye density corrections to the case of mixtures and numerically agrees with the HNC approximation in that regime.

preprint2009arXiv

Episodic accretion at early stages of evolution of low mass stars and brown dwarfs: a solution for the observed luminosity spread in HR diagrams?

We present evolutionary models for young low mass stars and brown dwarfs taking into account episodic phases of accretion at early stages of the evolution, a scenario supported by recent large surveys of embedded protostars. An evolution including short episodes of vigorous accretion ($\mdot \ge 10^{-4} \msolyr$) followed by longer quiescent phases ($\mdot < 10^{-6} \msolyr$) can explain the observed luminosity spread in HR diagrams of star forming regions at ages of a few Myr, for objects ranging from a few Jupiter masses to a few tenths of a solar mass. The gravitational contraction of these accreting objects strongly departs from the standard Hayashi track at constant $\te$. The best agreement with the observed luminosity scatter is obtained if most of the accretion shock energy is radiated away. The obtained luminosity spread at 1 Myr in the HR diagram is equivalent to what can be misinterpreted as a $\sim$ 10 Myr age spread for non-accreting objects. We also predict a significant spread in radius at a given $\te$, as suggested by recent observations. These calculations bear important consequences on our understanding of star formation and early stages of evolution and on the determination of the IMF for young ($\le$ a few Myr) clusters. Our results also show that the concept of a stellar birthline for low-mass objects has no valid support.

preprint2009arXiv

Protostellar collapse: radiative and magnetic feedbacks on small scale fragmentation

It is established that both radiative transfer and magnetic field have a strong impact on the collapse and the fragmentation of prestellar dense cores, but no consistent calculation exists yet at such scales. We present original AMR calculations including magnetic field (in the ideal MHD limit) and radiative transfer, within the Flux Limited Diffusion approximation, of the collapse of a 1 solar mass dense core. We compare the results with calculations performed with a barotropic EOS. We show that radiative transfer has an important impact on the collapse and the fragmentation, through the cooling or heating of the gas, and is complementary of the magnetic field. A larger field yields a stronger magnetic braking, increasing the accretion rate and thus the effect of the radiative feedback. Even for a strongly magnetized core, where the dynamics of the collapse is dominated by the magnetic field, radiative transfer is crucial to determine the temperature and optical depth distributions, two potentially accessible observational diagnostics. A barotropic EOS cannot account for realistic fragmentation. The diffusivity of the numerical scheme, however, is found to strongly affect the output of the collapse, leading eventually to spurious fragmentation. Both radiative transfer and magnetic field must be included in numerical calculations of star formation to obtain realistic collapse configurations and observable signatures. Nevertheless, the numerical resolution and the robustness of the solver are of prime importance to obtain reliable results. When using an accurate solver, the fragmentation is found to always remain inhibited by the magnetic field, at least in the ideal MHD limit, even when radiative transfer is included.

preprint2009arXiv

Structure and evolution of the first CoRoT exoplanets: Probing the Brown Dwarf/Planet overlapping mass regime

We present detailed structure and evolution calculations for the first transiting extrasolar planets discovered by the space-based CoRoT mission. Comparisons between theoretical and observed radii provide information on the internal composition of the CoRoT objects. We distinguish three different categories of planets emerging from these discoveries and from previous ground-based surveys: (i) planets explained by standard planetary models including irradiation, (ii) abnormally bloated planets and (iii) massive objects belonging to the overlapping mass regime between planets and brown dwarfs. For the second category, we show that tidal heating can explain the relevant CoRoT objects, providing non-zero eccentricities. We stress that the usual assumption of a quick circularization of the orbit by tides, as usually done in transit light curve analysis, is not justified a priori, as suggested recently by Levrard et al. (2009), and that eccentricity analysis should be carefully redone for some observations. Finally, special attention is devoted to CoRoT-3b and to the identification of its very nature: giant planet or brown dwarf ? The radius determination of this object confirms the theoretical mass-radius predictions for gaseous bodies in the substellar regime but, given the present observational uncertainties, does not allow an unambiguous identification of its very nature. This opens the avenue, however, to an observational identification of these two distinct astrophysical populations, brown dwarfs and giant planets, in their overlapping mass range, as done for the case of the 8 Jupiter-mass object Hat-P-2b. (abridged)

preprint2008arXiv

Structure and evolution of super-Earth to super-Jupiter exoplanets: I. heavy element enrichment in the interior

We examine the uncertainties in current planetary models and we quantify their impact on the planet cooling histories and mass-radius relationships. These uncertainties include (i) the differences between the various equations of state used to characterize the heavy material thermodynamical properties, (ii) the distribution of heavy elements within planetary interiors, (iii) their chemical composition and (iv) their thermal contribution to the planet evolution. Our models, which include a gaseous H/He envelope, are compared with models of solid, gasless Earth-like planets in order to examine the impact of a gaseous envelope on the cooling and the resulting radius. We find that for a fraction of heavy material larger than 20% of the planet mass, the distribution of the heavy elements in the planet's interior affects substantially the evolution and thus the radius at a given age. For planets with large core mass fractions ($\simgr$ 50%), such as the Neptune-mass transiting planet GJ436b, the contribution of the gravitational and thermal energy from the core to the planet cooling history is not negligible, yielding a $\sim$ 10% effect on the radius after 1 Gyr. We show that the present mass and radius determinations of the massive planet Hat-P-2b require at least 200 $\mearth$ of heavy material in the interior, at the edge of what is currently predicted by the core-accretion model for planet formation. We show that if planets as massive as $\sim$ 25 $\mjup$ can form, as predicted by improved core-accretion models, deuterium is able to burn in the H/He layers above the core, even for core masses as large as $\sim$ 100 $\mearth$. We provide extensive grids of planetary evolution models from 10 $\mearth$ to 10 M$_{\rm Jup}$, with various fractions of heavy elements.

preprint2007arXiv

Fundamental properties of low-mass stars

Numerous investigations on the fundamental properties of low-mass stars using eclipsing binaries indicate a strong discrepancy between theory and observations that is still awaiting explanation. Current models seem to predict radii for stars less massive than the Sun that are some 10% smaller than observed, while their effective temperatures are some 5% larger. Here we discuss recent new observational data that are relevant to this issue and review the progress made in understanding the origin of the important differences with theoretical calculations. Notably, we provide evidence that stellar activity may be responsible for the mismatch between observations and theory through two different channels: inhibition of convection or effects of a significant starspot coverage. The activity hypothesis is put to a test with observational diagnostics and some of the consequences of the large starspot coverage are evaluated. We conclude that stellar activity likely plays a key role in defining the properties of active low-mass stars and that this should be properly taken into account when investigating young, active stars in clusters or star-forming regions.

preprint2007arXiv

The stable magnetic field of the fullly-convective star V374 Peg

We report in this paper phase-resolved spectropolarimetric observations of the rapidly-rotating fully-convective M4 dwarf V374 Peg, on which a strong, mainly axisymmetric, large-scale poloidal magnetic field was recently detected. In addition to the original data set secured in 2005 August, we present here new data collected in 2005 September and 2006 August. From the rotational modulation of unpolarised line profiles, we conclude that starspots are present at the surface of the star, but their contrast and fractional coverage are much lower than those of non-fully convective active stars with similar rotation rate. Applying tomographic imaging on each set of circularly polarised profiles separately, we find that the large-scale magnetic topology is remarkably stable on a timescale of 1 yr; repeating the analysis on the complete data set suggests that the magnetic configuration is sheared by very weak differential rotation (about 1/10th of the solar surface shear) and only slightly distorted by intrinsic variability. This result is at odds with various theoretical predictions, suggesting that dynamo fields of fully-convective stars should be mostly non-axisymmetric unless they succeed at triggering significant differential rotation.

preprint2005arXiv

Review on low-mass stars and brown dwarfs

In this review, we examine the successes and weaknesses of modern low-mass star and brown dwarf theory.(1) We first focus on the mechanical (equation of state) and thermal (atmosphere) properties and on the evolution. We then examine the current shortcomings of the theory and we discuss recent observational analysis which have suggested discrepancies between models and observations.(2) We then examine the stellar and brown dwarf IMF and suggest that a power-law above the average thermal Jeans mass (about 1 Msol) rolling over a lognormal form below this limit adequately reproduces the observations of field and young cluster stellar and brown dwarf distributions. This yields a reasonably accurate estimate of the stellar and brown dwarf Galactic census. Finally (3) we argue that the combination of turbulence driven fragmentation at large scale and gravity at small scales provides an appealing solution for the general star and brown dwarf formation mechanism. It also provides a physical ground for the aforementioned power-law + lognormal form for the IMF, whereas a series of different power laws lacks such a physical motivation. At last, we argue that the deuterium-burning limit as the distinction between stars and planets has no physical foundation in this modern star formation scheme. Opacity limited fragmentation extending down to a few (< 10) jupiter masses, due to shocks, anisotropy or magnetic fields, provides a much more robust limit, even though difficult to determine accurately. Therefore, the various "direct" detections of exoplanets claimed recently in the literature are most likely regular low-mass brown dwarfs and the direct detection of an extrasolar planet remains for now elusive.

preprint2001arXiv

The Galactic disk mass-budget : I. stellar mass-function and density

In this paper, we use the general theory worked out within the past few years for the structure and the evolution of low-mass stars to derive the stellar mass-function in the Galactic disk down to the vicinity of the hydrogen-burning limit, from the observed nearby luminosity functions. The accuracy of the mass-magnitude relationships derived from the afore-mentioned theory is examined by comparison with recent, accurate observational relationships in the M-dwarf domain. The mass function is shown to flatten out below $\sim 1 \msol$ but to keep rising down to the bottom of the main sequence. Combining the present determination below 1 $\msol$ and Scalo's (1986) mass function for larger masses, we show that the mass function is well described over the entire stellar mass range, from $\sim 100 \msol$ to $\sim 0.1 \msol$, by three functional forms, namely a two-segment power-law, a log-normal form or an exponential form, all normalized to the Hipparcos sample at 0.8 $\msol$. Integration of this mass function yields a reasonably accurate census of the entire stellar population in the Galactic disk, and its volume and surface mass-density.

preprint2001arXiv

The Galactic disk mass-budget : II. Brown dwarf mass-function and density

In this paper, we extend the calculations conducted previously in the stellar regime to determine the brown dwarf IMF in the Galactic disk. We perform Monte Carlo calculations taking into account the brown dwarf formation rate, spatial distribution and binary fraction. Comparison with existing surveys seems to exclude a power-law MF as steep as the one determined in the stellar regime below 1 $\msol$ and tends to favor a more flatish behaviour. Comparison with methane-dwarf detections tends to favor an eventually decreasing form like the lognormal or the more general exponential distributions determined in the previous paper. We calculate predicting brown dwarf counts in near-infrared color diagrams and brown dwarf discovery functions. These calculations yield the presently most accurate determination of the brown dwarf census in the Galactic disk. The brown dwarf number density is comparable to the stellar one, $n_{BD}\simeq n_\star\simeq 0.1$ pc$^{-3}$. The corresponding brown dwarf mass density, however, represents only about 10% of the stellar contribution, i.e. $ρ_{BD}\simle 5.0\times 10^{-3} \mvol$. Adding up the local stellar density determined previously yields the density of star-like objects, stars and brown dwarfs, in the solar neighborhood $ρ_\odot \approx 5.0\times 10^{-2} \mvol$.

preprint2001arXiv

The Substellar Mass Function in sigma Orionis

We combine results from imaging searches for substellar objects in the sigma Orionis cluster and follow-up photometric and spectroscopic observations to derive a census of the brown dwarf population in a region of 847 arcmin^2. We identify 64 very low-mass cluster member candidates in this region. We have available three color (IZJ) photometry for all of them, spectra for 9 objects, and K photometry for 27% of our sample. These data provide a well defined sequence in the I vs I-J, I-K color magnitude diagrams, and indicate that the cluster is affected by little reddening despite its young age (~5 Myr). Using state-of-the-art evolutionary models, we derive a mass function from the low-mass stars (0.2 Msol) across the complete brown dwarf domain (0.075 Msol to 0.013 Msol), and into the realm of free-floating planetary-mass objects (<0.013 Msol). We find that the mass spectrum (dN/dm ~ m^{-alpha}) increases toward lower masses with an exponent alpha = 0.8+/-0.4. Our results suggest that planetary-mass isolated objects could be as common as brown dwarfs; both kinds of objects together would be as numerous as stars in the cluster. If the distribution of stellar and substellar masses in sigma Orionis is representative of the Galactic disk, older and much lower luminosity free-floating planetary-mass objects with masses down to about 0.005 Msol should be abundant in the solar vicinity, with a density similar to M-type stars.

preprint1999arXiv

Modeling Pressure-Ionization of Hydrogen in the Context of Astrophysics

The recent development of techniques for laser-driven shock compression of hydrogen has opened the door to the experimental determination of its behavior under conditions characteristic of stellar and planetary interiors. The new data probe the equation of state (EOS) of dense hydrogen in the complex regime of pressure ionization. The structure and evolution of dense astrophysical bodies depend on whether the pressure ionization of hydrogen occurs continuously or through a ``plasma phase transition'' (PPT) between a molecular state and a plasma state. For the first time, the new experiments constrain predictions for the PPT. We show here that the EOS model developed by Saumon and Chabrier can successfully account for the data, and we propose an experiment that should provide a definitive test of the predicted PPT of hydrogen. The usefulness of the chemical picture for computing astrophysical EOS and in modeling pressure ionization is discussed.

preprint1998arXiv

Equation of state of fully ionized electron-ion plasmas

Thermodynamic quantities of Coulomb plasmas consisting of point-like ions immersed in a compressible, polarizable electron background are calculated for ion charges Z=1 to 26 and for a wide domain of plasma parameters ranging from the Debye-Hueckel limit to the crystallization point and from the region of nondegenerate to fully degenerate nonrelativistic or relativistic electrons. The calculations are based on the linear-response theory for the electron-ion interaction, including the local-field corrections in the electronic dielectric function. The thermodynamic quantities are calculated in the framework of the N-body hypernetted-chain equations and fitted by analytic expressions. We present also accurate analytic approximations for the free energy of the ideal electron gas at arbitrary degeneracy and relativity and for the excess free energy of the one-component plasma of ions (OCP) derived from Monte Carlo simulations. The extension to multi-ionic mixtures is discussed within the framework of the linear mixing rule. These formulae provide a completely analytic, accurate description of the thermodynamic quantities of fully ionized electron-ion Coulomb plasmas, a useful tool for various applications from liquid state theory to dense stellar matter.