Search

Juno's observations of Jupiter's gravity field have revealed extremely low values for the gravitational moments that are difficult to reconcile with the high abundance of metals observed in the atmosphere by both Galileo and Juno. Recent studies chose to arbitrarily get rid of one of these two constraints in order to build models of Jupiter.In this presentation, I will detail our new Jupiter structure models reconciling Juno and Galileo observational constraints. These models confirm the need to separate Jupiter into at least 4 layers: an outer convective shell, a non-convective zone of compositional change, an inner convective shell and a diluted core representing about 60 percent of the planet in radius. Compared to other studies, these models propose a new idea with important consequences: a decrease in the quantity of metals between the outer and inner convective shells. This would imply that the atmospheric composition is not representative of the internal composition of the planet, contrary to what is regularly admitted, and would strongly impact the Jupiter formation scenarios (localization, migration, accretion).In particular, the presence of an internal non-convective zone prevents mixing between the two convective envelopes. I will detail the physical processes of this semi-convective zone (layered convection or H-He immiscibility) and explain how they may persist during the evolution of the planet.These models also impose a limit mass on the compact core, which cannot be heavier than 5 Earth masses. Such a mass, lower than the runaway gas accretion minimum mass, needs to be explained in the light of our understanding of the formation and evolution of giant planets.I will finally detail the application of our work to Saturn, and what we can expect to learn about the interior of the giant planets in the years to come.

Interior models of giant planets traditionally assume that at a given radius (i.e. pressure) the density should be larger than or equal to the one corresponding to a homogeneous, adiabatic stratification throughout the planet (referred to as the 'outer adiabat'). The observations of Jupiter's gravity field by Juno combined with the constraints on its atmospheric composition appear to be incompatible with such a profile. In this letter, we show that the above assumption stems from an incorrect understanding of the Schwarzschild-Ledoux criterion, which is only valid on a local scale. In order to fulfil the buoyancy stability condition, the density gradient with pressure in a non-adiabatic region must indeed rise more steeply than the {\it local} adiabatic density gradient. However, the density gradient can be smaller than the one corresponding to the outer adiabat at the same pressure because of the higher temperature in an inhomogeneously stratified medium. Deep enough, the density can therefore be lower than the one corresponding to the outer adiabat. We show that this is permitted only if the slope of the local adiabat becomes shallower than the slope of the outer adiabat at the same pressure, as found in recent Jupiter models due to the increase of both specific entropy and adiabatic index with depth. We examine the dynamical stability of this structure and show that it is stable against non-adiabatic perturbations. The possibility of such unconventional density profile in Jupiter complicates further our understanding of the internal structure and evolution of (extrasolar) giant planets., Accepted in ApJL

In this Letter, we study the evolution of the autocovariance function (ACF) of density field fluctuations in star-forming clouds and thus of the correlation length $l_c(\rho)$ of these fluctuations, which can be identified as the average size of the most correlated structures within the cloud. Generalizing the transport equation derived by Chandrasekhar (1951) for static, homogeneous turbulence, we show that the mass contained within these structures is an invariant, i.e. that the average mass contained in the most correlated structures remains constant during the evolution of the cloud, whatever dominates the global dynamics (gravity or turbulence). We show that the growing impact of gravity on the turbulent flow yields an increase of the variance of the density fluctuations and thus a drastic decrease of the correlation length. Theoretical relations are successfully compared to numerical simulations. This picture brings a robust support to star formation paradigms where the mass concentration in turbulent star-forming clouds evolves from initially large, weakly correlated filamentary structures to smaller, denser more correlated ones, and eventually to small, tightly correlated prestellar cores. We stress that the present results rely on a pure statistical approach of density fluctuations and do not involve any specific condition for the formation of prestellar cores. Interestingly enough, we show that, under average conditions typical of Milky Way molecular clouds, this invariant average mass is about a solar mass, providing an appealing explanation for the apparent universality of the IMF under such environments., Comment: Accepted for publications in APJ Letters

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., Comment: 17 pages, 20 figures,accepted for publication in A&A

Juno's observations of Jupiter's gravity field have revealed extremely low values for the gravitational moments that are difficult to reconcile with the high abundance of metals observed in the atmosphere by Galileo. Recent studies chose to arbitrarily get rid of one of these two constraints in order to build models of Jupiter. In this presentation, I will detail our new Jupiter structure models reconciling Juno and Galileo observational constraints. These models confirm the need to separate Jupiter into at least 4 layers: an outer convective shell, a non-convective zone of compositional change, an inner convective shell and a diluted core representing about 60 percent of the planet in radius. Compared to other studies, these models propose a new idea with important consequences: a decrease in the quantity of metals between the outer and inner convective shells. This would imply that the atmospheric composition is not representative of the internal composition of the planet, contrary to what is regularly admitted, and would strongly impact the Jupiter formation scenarios (localization, migration, accretion). In particular, the presence of an internal non-convective zone prevents mixing between the two convective envelopes. I will detail the physical processes of this semi-convective zone (layered convection or H-He immiscibility) and explain how they may persist during the evolution of the planet. These models also impose a limit mass on the compact core, which cannot be heavier than 5 Earth masses. Such a mass, lower than the runaway gas accretion minimum mass, needs to be explained in the light of our understanding of the formation and evolution of giant planets. Using these models of Jupiter, I will finally detail the application of our new understanding of the interior of this planet to giant exoplanets. At a time of direct imaging of extrasolar planets and atmospheric characterization of hot Jupiters, a good understanding of the internal processes of planets in the solar system is paramount to make the best use of all the observations.

Context: The anomalously large radii of hot Jupiters has long been a mystery. However, by combining both theoretical arguments and 2D models, a recent study has suggested that the vertical advection of potential temperature leads to an adiabatic temperature profile in the deep atmosphere hotter than the profile obtained with standard 1D models. Aims: In order to confirm the viability of that scenario, we extend this investigation to three dimensional, time-dependent, models. Methods: We use a 3D GCM, DYNAMICO to perform a series of calculations designed to explore the formation and structure of the driving atmospheric circulations, and detail how it responds to changes in both upper and deep atmospheric forcing. Results: In agreement with the previous, 2D, study, we find that a hot adiabat is the natural outcome of the long-term evolution of the deep atmosphere. Integration times of order $1500$ years are needed for that adiabat to emerge from an isothermal atmosphere, explaining why it has not been found in previous hot Jupiter studies. Models initialised from a hotter deep atmosphere tend to evolve faster toward the same final state. We also find that the deep adiabat is stable against low-levels of deep heating and cooling, as long as the Newtonian cooling time-scale is longer than $\sim 3000$ years at $200$ bar. Conclusions: We conclude that the steady-state vertical advection of potential temperature by deep atmospheric circulations constitutes a robust mechanism to explain hot Jupiter inflated radii. We suggest that future studies of hot Jupiters are evolved for a longer time than currently done, and, when possible, include models initialised with a hot deep adiabat. We stress that this mechanism stems from the advection of entropy by irradiation induced mass flows and does not require (finely tuned) dissipative process, in contrast with most previously suggested scenarios., Comment: 13 pages, 12 figures. Accepted for publication in A&A

In a recent paper (Chabrier et al. 2019), we have derived a new equation of state (EOS) for dense hydrogen/helium mixtures which covers the temperature-density domain from solar-type stars to brown dwarfs and gaseous planets. This EOS is based on the so-called additive volume law and thus does not take into account the interactions between the hydrogen and helium species. In the present paper, we go beyond these calculations by taking into account H/He interactions, derived from quantum molecular dynamics simulations. These interactions, which eventually lead to H/He phase separation, become important at low temperature and high density, in the domain of brown dwarfs and giant planets. The tables of this new EOS are made publicly available., Comment: To be published in The Astrophysical Journal

The study of brown dwarfs and giant exoplanets is rapidly evolving as ever-improving instrumentation becomes sensitive to cooler objects. Accurate and reliable atmosphere and evolutionary models are important for placing mass and age constraints on these newly discovered objects, and understanding the rich chemistry and physics taking place in their atmospheres. We are expanding on the widely used COND evolutionary models by developing a grid of model atmospheres (Teff=200-2000K, log(g)=2.5-5.5) with our state-of-the-art 1D radiative-convective equilibrium code ATMO. ATMO includes the latest opacities for important molecular absorbers such as H2O, CH4and NH3, and the latest line shapes for the collisionally broadened Na and K resonance lines.These model improvements allow us to follow the evolution of 1-75MJupobjects down to the coolest effective temperatures (Teff=200K). I will present comparisons of this new model set to those previously published, illustrating how the evolutionary tracks and the substellar boundary have changed due to improved opacities and the usage of a new equation of state. We also compare our model grid to observational datasets in colour-magnitude diagrams and investigate the impact of our new Na and K line shapes in reproducing brown dwarf spectra. Our future work will involve expanding on this initial grid, to investigate the effects of metallicity, C/O ratio and radiative diabatic convection in cool brown dwarfs and giant exoplanets.

The stellar initial mass function (IMF) is playing a critical role in the history of our universe. We propose a theory that is based solely on local processes, namely the dust opacity limit, the tidal forces and the properties of the collapsing gas envelope. The idea is that the final mass of the central object is determined by the location of the nearest fragments, which accrete the gas located further away, preventing it to fall onto the central object. To estimate the relevant statistics in the neighbourhood of an accreting protostar, we perform high resolution numerical simulations. We also use these simulations to further test the idea that fragmentation in the vicinity of an existing protostar is determinant in setting the peak of the stellar mass spectrum. We develop an analytical model, which is based on a statistical counting of the turbulent density fluctuations, generated during the collapse, that are at least equal to the mass of the first hydrostatic core, and sufficiently important to supersede tidal and pressure forces to be self-gravitating. The analytical mass function presents a peak located at roughly 10 times the mass of the first hydrostatic core in good agreement with the numerical simulations. Since the physical processes involved are all local, i.e. occurs at scales of a few 100 AU or below, and do not depend on the gas distribution at large scale and global properties such as the mean Jeans mass, the mass spectrum is expected to be relatively universal., accepted for publication in ApJ

We present individual dynamical masses for the nearby M9.5+T5.5 binary WISE J072003.20$-$084651.2AB, a.k.a. Scholz's star. Combining high-precision CFHT/WIRCam photocenter astrometry and Keck adaptive optics resolved imaging, we measure the first high-quality parallactic distance ($6.80_{-0.06}^{+0.05}$ pc) and orbit ($8.06_{-0.25}^{+0.24}$ yr period) for this system composed of a low-mass star and brown dwarf. We find a moderately eccentric orbit ($e = 0.240_{-0.010}^{+0.009}$), incompatible with previous work based on less data, and dynamical masses of $99\pm6$ $M_{\rm Jup}$ and $66\pm4$ $M_{\rm Jup}$ for the two components. The primary mass is marginally inconsistent (2.1$\sigma$) with the empirical mass$-$magnitude$-$metallicity relation and models of main-sequence stars. The relatively high mass of the cold ($T_{\rm eff} = 1250\pm40$ K) brown dwarf companion indicates an age older than a few Gyr, in accord with age estimates for the primary star, and is consistent with our recent estimate of $\approx$70 $M_{\rm Jup}$ for the stellar/substellar boundary among the field population. Our improved parallax and proper motion, as well as an orbit-corrected system velocity, improve the accuracy of the system's close encounter with the solar system by an order of magnitude. WISE J0720$-$0846AB passed within $68.7\pm2.0$ kAU of the Sun $80.5\pm0.7$ kyr ago, passing through the outer Oort cloud where comets can have stable orbits., Comment: accepted to AJ

We present a new equation of state (EOS) for dense hydrogen/helium mixtures which covers a range of densities from $10^{-8}$ to $10^6$ g.cm$^{-3}$, pressures from $10^{-9}$ to $10^{13}$ GPa and temperatures from $10^{2}$ to $10^{8}$ K. The calculations combine the EOS of Saumon, Chabrier & vanHorn (1995) in the low density, low temperature molecular/atomic domain, the EOS of Chabrier & Potekhin (1998) in the high-density, high-temperature fully ionized domain, the limits of which differ for H and He, and ab initio quantum molecular dynamics (QMD) calculations in the intermediate density and temperature regime, characteristic of pressure dissociation and ionization. The EOS for the H/He mixture is based on the so-called additive volume law and thus does not take into account the interactions between the two species. A major improvement of the present calculations over existing ones is that we calculate the entropy over the entire density-temperature domain, a necessary quantity for stellar or planetary evolution calculations. The EOS results are compared with existing experimental data, namely Hugoniot shock experiments for pure H and He, and with first principle numerical simulations for both the single elements and the mixture. This new EOS covers a wide range of physical and astrophysical conditions, from jovian planets to solar-type stars, and recovers the existing relativistic EOS at very high densities, in the domains of white dwarfs and neutron stars., To appear in Astrophysical Journal. Figures in the published version will be larger

By generalizing the theory of convection to any type of thermal and compositional source terms (diabatic processes), we show that thermohaline convection in Earth oceans, fingering convection in stellar atmospheres, and moist convection in Earth atmosphere are deriving from the same general diabatic convective instability. We show also that "radiative convection" triggered by CO/CH4 transition with radiative transfer in the atmospheres of brown dwarfs is analog to moist and thermohaline convection. We derive a generalization of the mixing length theory to include the effect of source terms in 1D codes. We show that CO/CH4 radiative convection could significantly reduce the temperature gradient in the atmospheres of brown dwarfs similarly to moist convection in Earth atmosphere thus possibly explaining the reddening in brown-dwarf spectra. By using idealized two-dimensional hydrodynamic simulations in the Ledoux unstable regime, we show that compositional source terms can indeed provoke a reduction of the temperature gradient. The L/T transition could be explained by a bifurcation between the adiabatic and diabatic convective transports and could be seen as a giant cooling crisis: an analog of the boiling crisis in liquid/steam-water convective flows. This mechanism with other chemical transitions could be present in many giant and earth-like exoplanets. The study of the impact of different parameters (effective temperature, compositional changes) on CO/CH4 radiative convection and the analogy with Earth moist and thermohaline convection is opening the possibility to use brown dwarfs to better understand some aspects of the physics at play in the climate of our own planet., accepted in ApJ

We present here a minor modification of our numerical implementation of the Hall effect for the 2D Riemann solver used in Constrained Transport schemes, as described in Marchand et al. (2018). In the previous work, the tests showed that the angular momentum was not conserved during protostellar collapse simulations, with significant impact. By removing the whistler waves speed from the characteristic speeds of non-magnetic variables in the 1D Riemann solver, we are able to improve the angular momentum conservation in our test-case by one order of magnitude, while keeping the second-order numerical convergence of the scheme. We also reproduce the simulations of Tsukamoto et al. (2015) with consistent resistivities, the three non-ideal MHD effects and initial rotation, and agree with their results. In this case, the violation of angular momentum conservation is negligible in regard to the total angular momentum and the angular momentum of the disk., Comment: 6 pages, 6 figures, published in A&A

International audience; Observations of Jupiter’s gravity field by Juno have revealed surprisingly low values for the high-order gravitational moments, considering the abundances of heavy elements measured by Galileo 20 years ago. The derivation of recent equations of state for hydrogen and helium, which are much denser in the megabar region, exacerbates the conflict between these two observations. In order to circumvent this puzzle, current Jupiter model studies either ignore the constraint from Galileo or invoke an ad hoc modification of the equations of state. In this paper, we derive Jupiter models that satisfy constraints of both Juno and Galileo. We confirm that Jupiter’s structure must encompass at least four different regions: an outer convective envelope, a region of compositional and thus entropy change, an inner convective envelope, an extended diluted core enriched in heavy elements, and potentially a central compact core. We show that in order to reproduce Juno and Galileo observations, one needs a significant entropy increase between the outer and inner envelopes and a lower density than for an isentropic profile, which is associated with some external differential rotation. The best way to fulfill this latter condition is an inward-decreasing abundance of heavy elements in this region. We examine in detail the three physical mechanisms that can yield such a change of entropy and composition: a first-order molecular-metallic hydrogen transition, immiscibility between hydrogen and helium, or a region of layered convection. Given our present knowledge of hydrogen pressure ionization, a combination of the two latter mechanisms seems to be the most favored solution.

We study long-term thermal evolution of neutron stars in soft X-ray transients (SXTs), taking the deep crustal heating into account consistently with the changes of the composition of the crust. We collect observational estimates of average accretion rates and thermal luminosities of such neutron stars and compare the theory with observations. We perform simulations of thermal evolution of accreting neutron stars, considering the gradual replacement of the original nonaccreted crust by the reprocessed accreted matter, the neutrino and photon energy losses, and the deep crustal heating due to nuclear reactions in the accreted crust. We test and compare results for different modern theoretical models. We update a compilation of the observational estimates of the thermal luminosities in quiescence and average accretion rates in the SXTs and compare the observational estimates with the theoretical results. Long-term thermal evolution of transiently accreting neutron stars is nonmonotonic. The quasi-equilibrium temperature in quiescence reaches a minimum and then increases toward the final steady state. The quasi-equilibrium thermal luminosity of a neutron star in an SXT can be substantially lower at the minimum than in the final state. This enlarges the range of possibilities for theoretical interpretation of observations of such neutron stars. The updates of the theory and observations leave unchanged the previous conclusions that the direct Urca process operates in relatively cold neutron stars and that an accreted heat-blanketing envelope is likely present in relatively hot neutron stars in the SXTs in quiescence. The results of the comparison of theory with observations favor suppression of the triplet pairing type of nucleon superfluidity in the neutron-star matter., 16 pages, 2 tables, 14 figures. In v4, Eq.(18) is corrected (results are not affected)

As a first step toward a multi-phase equation of state for dense water, we develop a temperature-dependent equation of state for dense water covering the liquid and plasma regimes and extending to the super-ionic and gas regimes. This equation of state covers the complete range of conditions encountered in planetary modeling. We use first principles quantum molecular dynamics simulations and its Thomas-Fermi extension to reach the highest pressures encountered in giant planets several times the size of Jupiter. Using these results, as well as the data available at lower pressures, we obtain a parametrization of the Helmholtz free energy adjusted over this extended temperature and pressure domain. The parametrization ignores the entropy and density jumps at phase boundaries but we show that it is sufficiently accurate to model interior properties of most planets and exoplanets. We produce an equation of state given in analytical form that is readily usable in planetary modeling codes and dynamical simulations (a fortran implementation can be found at http://www.ioffe.ru/astro/H2O/). The EOS produced is valid for the entire density range relevant to planetary modeling, for densities where quantum effects for the ions can be neglected, and for temperatures below 50,000K. We use this equation of state to calculate the mass-radius relationship of exoplanets up to 5,000M_Earth, explore temperature effects in ocean and wet Earth-like planets, and quantify the influence of the water EOS for the core on the gravitational moments of Jupiter., Comment: 13 pages, 18 figures. In v.2, Eq.(13) and the value of S_0 are corrected

In spite of decades of theoretical efforts, the physical origin of the stellar initial mass function (IMF) is still debated. Particularly crucial is the question of what sets the peak of the distribution. To investigate this issue we perform high resolution numerical simulations with radiative feedback exploring in particular the role of the stellar and accretion luminosities. We also perform simulations with a simple effective equation of state (eos) and we investigate 1000 solar mass clumps having respectively 0.1 and 0.4 pc of initial radii. We found that most runs, both with radiative transfer or an eos, present similar mass spectra with a peak broadly located around 0.3-0.5 M$_\odot$ and a powerlaw-like mass distribution at higher masses. However, when accretion luminosity is accounted for, the resulting mass spectrum of the most compact clump tends to be moderately top-heavy. The effect remains limited for the less compact one, which overall remains colder. Our results support the idea that rather than the radiative stellar feedback, this is the transition from the isothermal to the adiabatic regime, which occurs at a gas density of about 10$^{10}$ cm$^{-3}$, that is responsible for setting the peak of the initial mass function. This stems for the fact that $i)$ extremely compact clumps for which the accretion luminosity has a significant influence are very rare and $ii)$ because of the luminosity problem, which indicates that the effective accretion luminosity is likely weaker than expected., Comment: accepted for publication in ApJ

We derive an analytical theory of the PDF of density fluctuations in supersonic turbulence in the presence of gravity in star-forming clouds. The theory is based on a rigorous derivation of a combination of the Navier-Stokes continuity equations for the fluid motions and the Poisson equation for the gravity. It extends upon previous approaches first by including gravity, second by considering the PDF as a dynamical system, not a stationary one. We derive the transport equations of the density PDF, characterize its evolution and determine the density threshold above which gravity strongly affects and eventually dominates the dynamics of turbulence. We demonstrate the occurence of {\it two} power law tails in the PDF, with two characteristic exponents, corresponding to two different stages in the balance between turbulence and gravity. Another important result of this study is to provide a procedure to relate the observed {\it column density} PDFs to the corresponding {\it volume density} PDFs. This allows to infer, from the observation of column-densities, various physical parameters characterizing molecular clouds, notably the virial parameter. Furthermore, the theory offers the possibility to date the clouds in units of ${t}_{\rm coll}$, the time since a statistically significant fraction of the cloud started to collapse. The theoretical results and diagnostics reproduce very well numerical simulations and observations of star-forming clouds. The theory provides a sound theoretical foundation and quantitative diagnostics to analyze observations or numerical simulations of star-forming regions and to characterize the evolution of the density PDF in various regions of molecular clouds., Comment: Accepted for publications in APJ Letters

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 $\alpha$ 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 $\mu$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., Comment: 20 pages, 16 figures, Accepted in ApJ

This is the author accepted manuscript. The final version is available from the American Astronomical Society via the DOI in this record

We present an implementation of the Hall term in the non-ideal magnetohydrodynamics equations into the adaptive-mesh-refinement code RAMSES to study its impact on star formation. Recent works show that the Hall effect heavily influences the regulation of the angular momentum in collapsing dense cores, strengthening or weakening the magnetic braking. Our method consists of a modification of the two-dimensional constrained transport scheme. Our scheme shows convergence of second-order in space and the frequency of the propagation of whistler waves is accurate. We confirm previous results, namely that during the collapse, the Hall effect generates a rotation of the fluid with a direction in the mid-plane that depends on the sign of the Hall resistivity, while counter-rotating envelopes develop on each side of the mid-plane. However, we find that the predictability of our numerical results is severely limited. The angular momentum is not conserved in any of our dense core-collapse simulations with the Hall effect: a large amount of angular momentum is generated within the first Larson core, a few hundred years after its formation, without compensation by the surrounding gas. This issue is not mentioned in previous studies and may be correlated to the formation of the accretion shock on the Larson core. We expect that this numerical effect could be a serious issue in star formation simulations., 17 pages, 21 figures. Accepted for publication in Astronomy & Astrophysics

We study the relative importance of several recent updates of microphysics input to the neutron star cooling theory and the effects brought about by superstrong magnetic fields of magnetars, including the effects of the Landau quantization in their crusts. We use a finite-difference code for simulation of neutron-star thermal evolution on timescales from hours to megayears with an updated microphysics input. The consideration of short timescales ($\lesssim1$ yr) is made possible by a treatment of the heat-blanketing envelope without the quasistationary approximation inherent to its treatment in traditional neutron-star cooling codes. For the strongly magnetized neutron stars, we take into account the effects of Landau quantization on thermodynamic functions and thermal conductivities. We simulate cooling of ordinary neutron stars and magnetars with non-accreted and accreted crusts and compare the results with observations. Suppression of radiative and conductive opacities in strongly quantizing magnetic fields and formation of a condensed radiating surface substantially enhance the photon luminosity at early ages, making the life of magnetars brighter but shorter. These effects together with the effect of strong proton superfluidity, which slows down the cooling of kiloyear-aged neutron stars, can explain thermal luminosities of about a half of magnetars without invoking heating mechanisms. Observed thermal luminosities of other magnetars are still higher than theoretical predictions, which implies heating, but the effects of quantizing magnetic fields and baryon superfluidity help to reduce the discrepancy., 18 pages, 14 figures. In v3, a typo in a graphics code has been fixed, which caused a vertical shift by 0.3 of cooling curves in the figures (except the inset in Fig.2). In v2 and v4, typos were fixed in the bibliography list and in Eqs.(1)-(3)

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 $\chi^{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., Comment: No changes in replacement, main paper now made visible on arXiv as earlier the published Erratum associated with the paper was instead visible. 34 pages, 28 figures. Accepted for Publication in MNRAS. Full grid of model transmission spectra and chemical abundances are available here, https://drive.google.com/drive/folders/1tvlgSWyEAA0cX_yxnA5URD4Zo3Zj4GcH

A few years ago, Hubbard (2012, 2013) presented an elegant, non-perturbative method, called concentric MacLaurin spheroid (CMS), to calculate with very high accuracy the gravitational moments of a rotating fluid body following a barotropic pressure-density relationship. Having such an accurate method is of great importance for taking full advantage of the Juno mission, and its extremely precise determination of Jupiter gravitational moments, to better constrain the internal structure of the planet. Recently, several authors have applied this method to the Juno mission with 512 spheroids linearly spaced in altitude. We demonstrate in this paper that such calculations lead to errors larger than Juno's error bars, invalidating the aforederived Jupiter models at the level required by Juno's precision. We show that, in order to fulfill Juno's observational constraints, at least 1500 spheroids must be used with a cubic, square or exponential repartition, the most reliable solutions. When using a realistic equation of state instead of a polytrope, we highlight the necessity to properly describe the outermost layers to derive an accurate boundary condition, excluding in particular a zero pressure outer condition. Providing all these constraints are fulfilled, the CMS method can indeed be used to derive Jupiter models within Juno's present observational constraints. However, we show that the treatment of the outermost layers leads to irreducible errors in the calculation of the gravitational moments and thus on the inferred physical quantities for the planet. We have quantified these errors and evaluated the maximum precision that can be reached with the CMS method in the present and future exploitation of Juno's data., Comment: 21 pages, 11 figures

Recently, Jao et al. (2018) discovered a gap in the mid-M dwarfs main sequence revealed by the analysis of Gaia data Release 2. They suggested the feature is linked to the onset of full convection in M dwarfs. Following the announcement of this discovery, MacDonald & Gizis (2018) proposed an explanation based on standard stellar evolution models. In this paper we re-examine the explanation suggested by MacDonald & Gizis (2018). We confirm that nuclear burning and mixing process of $^3$He provide the best explanation for the observed feature. We also find that a change in the energy transport from convection to radiation does not induce structural changes that could be visible. Regarding the very details of the process, however, we disagree with MacDonald & Gizis (2018) and propose a different explanation., 5 pages, 6 figures, accepted for publication in Astronomy & Astrophysics

The transport of angular momentum is capital during the formation of low-mass stars; too little removal and rotation ensures stellar densities are never reached, too much and the absence of rotation means no protoplanetary disks can form. Magnetic diffusion is seen as a pathway to resolving this long-standing problem. We investigate the impact of including resistive MHD in simulations of the gravitational collapse of a 1 solar mass gas sphere, from molecular cloud densities to the formation of the protostellar seed; the second Larson core. We used the AMR code RAMSES to perform two 3D simulations of collapsing magnetised gas spheres, including self-gravity, radiative transfer, and a non-ideal gas equation of state to describe H2 dissociation which leads to the second collapse. The first run was carried out under the ideal MHD approximation, while ambipolar and ohmic diffusion was incorporated in the second calculation. In the ideal MHD simulation, the magnetic field dominates the energy budget everywhere inside and around the first core, fueling interchange instabilities and driving a low-velocity outflow. High magnetic braking removes essentially all angular momentum from the second core. On the other hand, ambipolar and ohmic diffusion create a barrier which prevents amplification of the magnetic field beyond 0.1 G in the first Larson core which is now fully thermally supported. A significant amount of rotation is preserved and a small Keplerian-like disk forms around the second core. When studying the radiative efficiency of the first and second core accretion shocks, we found that it can vary by several orders of magnitude over the 3D surface of the cores. Magnetic diffusion is a pre-requisite to star-formation; it enables the formation of protoplanetary disks in which planets will eventually form, and also plays a determinant role in the formation of the protostar itself., 18 pages, 11 figures, accepted for publication in Astronomy & Astrophysics

Atmospheric modeling of low-gravity (VL-G) young brown dwarfs remains a challenge. The presence of very thick clouds has been suggested because of their extremely red near-infrared (NIR) spectra, but no cloud models provide a good fit to the data with a radius compatible with evolutionary models for these objects. We show that cloudless atmospheres assuming a temperature gradient reduction caused by fingering convection provides a very good model to match the observed VL-G NIR spectra. The sequence of extremely red colors in the NIR for atmospheres with effective temperature from ~2000 K down to ~1200 K is very well reproduced with predicted radii typical of young low-gravity objects. Future observations with NIRSPEC and MIRI on the James Webb Space Telescope (JWST) will provide more constrains in the mid-infrared, helping to confirm/refute whether or not the NIR reddening is caused by fingering convection. We suggest that the presence/absence of clouds will be directly determined by the silicate absorption features that can be observed with MIRI. JWST will therefore be able to better characterize the atmosphere of these hot young brown dwarfs and their low-gravity exoplanet analogues., Comment: Accepted in ApJ

International audience; Observations suggest that star formation in filamentary molecular clouds occurs in a two-step process, with the formation of filaments preceding that of prestellar cores and stars. Here, we apply the gravoturbulent fragmentation theory of Hennebelle & Chabrier to a filamentary environment, taking into account magnetic support. We discuss the induced geometrical effect on the cores, with a transition from 3D geometry at small scales to 1D at large ones. The model predicts the fragmentation behavior of a filament for a given mass per unit length (MpL) and level of magnetization. This core mass function (CMF) for individual filaments is then convolved with the distribution of filaments to obtain the final system CMF. The model yields two major results. (I) The filamentary geometry naturally induces a hierarchical fragmentation process, first into groups of cores, separated by a length equal to a few filament Jeans lengths, I.e., a few times the filament width. These groups then fragment into individual cores. (II) Non-magnetized filaments with high MpL are found to fragment excessively, at odds with observations. This is resolved by taking into account the magnetic field (treated simply as additional pressure support). The present theory suggests two complementary modes of star formation: although small (spherical or filamentary) structures will collapse directly into prestellar cores, according to the standard Hennebelle-Chabrier theory, the large (filamentary) ones, the dominant population according to observations, will follow the aforedescribed two-step process.

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., 4 pages, 2 figures ; to appear in "IAU Symposium 270: Computational Star formation", Eds. J. Alves, B. Elmegreen, J. Girart, V. Trimble

The anomalously large radii of strongly irradiated exoplanets have remained a major puzzle in astronomy. Based on a 2D steady state atmospheric circulation model, the validity of which is assessed by comparison to 3D calculations, we reveal a new mechanism, namely the advection of the potential temperature due to mass and longitudinal momentum conservation, a process occuring in the Earth's atmosphere or oceans. At depth, the vanishing heating flux forces the atmospheric structure to converge to a hotter adiabat than the one obtained with 1D calculations, implying a larger radius for the planet. Not only do the calculations reproduce the observed radius of HD209458b, but also the observed correlation between radius inflation and irradiation for transiting planets. Vertical advection of potential temperature induced by non uniform atmospheric heating thus provides a robust mechanism explaining the inflated radii of irradiated hot Jupiters., Comment: accepted in ApJ

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., Comment: 14 pages, 8 figures, accepted for publication in Astronomy and Astrophysics

This is the final version. Available from the American Astronomical Society via the DOI in this record

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., Comment: 19 pages, 12 figures, accepted for publication in JQSRT

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., 9 pages, 4 figures; Proceedings of the 14th International Conference on Physics of Nonideal Plasmas "PNP-14" (9-14 September 2012, Rostock, Germany)

We review the current status of the theory of thermal emission from the surface layers of neutron stars with strong magnetic fields B ~ 1010 - 1015 G, including formation of the spectrum in a partially ionized atmosphere and at a condensed surface. In particular, we describe recent progress in modeling partially ionized atmospheres of central compact objects in supernova remnants, which may have moderately strong fields B ~ 1010 - 1011 G. Special attention is given to polarization of thermal radiation emitted by a neutron star surface. Finally, we briefly describe applications of the theory to observations of thermally emitting isolated neutron stars

We report on the first ultracool dwarf discoveries from the UKIRT Infrared Deep Sky Survey (UKIDSS) Large Area Survey Early Data Release (LAS EDR), in particular the discovery of T dwarfs which are fainter and more distant than those found using the 2MASS and SDSS surveys. We aim to show that our methodologies for searching the ~27 sq degs of the LAS EDR are successful for finding both L and T dwarfs $via$ cross-correlation with the Sloan Digital Sky Survey (SDSS) DR4 release. While the area searched so far is small, the numbers of objects found shows great promise for near-future releases of the LAS and great potential for finding large numbers of such dwarfs. Ultracool dwarfs are selected by combinations of their YJH(K) UKIDSS colours and SDSS DR4 z-J and i-z colours, or, lower limits on these red optical/infrared colours in the case of DR4 dropouts. After passing visual inspection tests, candidates have been followed up by methane imaging and spectroscopy at 4m and 8m-class facilities. Our main result is the discovery following CH4 imaging and spectroscopy of a T4.5 dwarf, ULASJ 1452+0655, lying ~80pc distant. A further T dwarf candidate, ULASJ 1301+0023, has very similar CH4 colours but has not yet been confirmed spectroscopically. We also report on the identification of a brighter L0 dwarf, and on the selection of a list of LAS objects designed to probe for T-like dwarfs to the survey J-band limit. Our findings indicate that the combination of the UKIDSS LAS and SDSS surveys provide an excellent tool for identifying L and T dwarfs down to much fainter limits than previously possible. Our discovery of one confirmed and one probable T dwarf in the EDR is consistent with expectations from the previously measured T dwarf density on the sky., 7 pages, 5 figures, A&A in press

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., Accepted for publication in A&A, 14 pages, 26 figures

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., Accepted on March 30 2016 for publication in ApJ

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., Comment: Accepted in ApJL, comments welcome

Thanks to the combination of transit photometry and radial velocity doppler measurements, we are now able to constrain theoretical models of the structure and evolution of objects in the whole mass range between icy giants and stars, including the giant planet/brown dwarf overlapping mass regime (Leconte et al. 2009). In the giant planet mass range, the significant fraction of planets showing a larger radius than predicted by the models suggests that a missing physical mechanism which is either injecting energy in the deep convective zone or reducing the net outward thermal flux is taking place in these objects. Several possibilities have been suggested for such a mechanism: •downward transport of kinetic energy originating from strong winds generated at the planet's surface (Showman & Guillot 2002),•enhanced opacity sources in hot-Jupiter atmospheres (Burrows et al. 2007),•ohmic dissipation in the ionized atmosphere (Batygin & Stevenson 2010),•(inefficient) layered or oscillatory convection in the planet's interior (Chabrier & Baraffe 2007),•Tidal heating due to circularization of the orbit, as originally suggested by Bodenheimer, Lin & Mardling (2001). Here we first review the differences between current models of tidal evolution and their uncertainties. We then revisit the viability of the tidal heating hypothesis using a tidal model which treats properly the highly eccentric and misaligned orbits commonly encountered in exoplanetary systems. We stress again that the low order expansions in eccentricity often used in constant phase lag tidal models (i.e. constant Q) necessarily yields incorrect results as soon as the (present or initial) eccentricity exceeds ~ 0.2, as can be rigorously demonstrated from Kepler's equations.

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, 32 pages, 11 figures, 2 tables. Published in the Astrophysical Journal

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., 6 pages, 2 figures, Contributions to Plasma Physics, 2010, vol.50 (in press)

We calculate the full stellar-structural evolution of donors in AM CVn systems formed through the WD channel coupled to the binary’s evolution. Contrary to assumptions made in prior modelling, these donors are not fully convective over much of the AM CVn phase and do not evolve adiabatically under mass loss indefinitely. Instead, we identify three distinct phases of evolution: a mass transfer turn-on phase (during which Porb continues to decrease even after contact, the donor contracts, and the mass transfer rate accelerates to its maximum), a phase in which the donor expands adiabatically in response to mass loss, and a cooling phase beginning at Porb � 45–55 minutes during which the donor contracts. The physics that determines the behaviour in the first and third phases, both of which are new outcomes of this study, are discussed in some detail. We find the overall duration of the turn-on phase to be between � 10 4 -10 6 yrs, significantly longer than prior estimates. We predict the donor’s luminosity, L, and effective temperature, Teff. During the adiabatic expansion phase (ignoring irradiation effects), L � 10 −6 –10 −4 L⊙ and Teff � 1000–1800 K. However, the flux generated in the accretion flow dominates the donor’s intrinsic light at all times. The impact of irradiation on the donor extends the phase of adiabatic expansion to longer Porb, slows the contraction during the cooling phase, and alters the donor’s observational characteristics. Irradiated donors during the adiabatic phase can attain surface luminosities up to � 10 −2 L⊙. We argue that the turn-on and cooling phases both will leave significant imprints on the AM CVn population’s Porb-distribution. Finally, we show that the eclipsing AM CVn system SDSS J0926+3624 provides evidence that WD-channel systems with non-zero entropy donors contribute to the AM CVn population, and we discuss the observational signature of the donor in this system.

We briefly summarize the present status of theoretical and experimental investigations aimed at describing the thermodynamic properties of hydrogen and helium at high density. We confront various theoretical models to presently available experiments and we consider the astrophysical implications for giant planet interiors.

Jupiter-sized brown dwarfs are found in the solar neighborhood with effective temperature Teff as low as 250 K [1]. Iron, silicates, chlorides and sulfides condense in the atmospheres of the Teff ≈ 2000 K L-type and Teff ≈ 1000 K T-type dwarfs [2]. At the T-/Y-type boundary, Teff ≈ 500 K and atmospheres are clear [3]. The next species to condense are H2O at Teff ≈ 350 K and NH3 at Teff ≈ 200 K [4].

We present 0.95-1.80 $\mu$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. $\beta$ 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., Comment: Accepted to ApJ Letters, 8 pages, 4 figures, some minor typographical issues were fixed

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., 7 pages, 8 figures, Astronomy & Astrophysics in press