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.

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

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

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.

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.

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.

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.

The lowest-mass stars, brown dwarfs and giant exoplanets span a minimum in the mass-radius relationship that probes the fundamental physics of extreme states of matter, magnetism, and fusion. This White Paper outlines scientific opportunities and the necessary resources for modeling and measuring the mass-radius relationship in this regime.

In the original manuscript arXiv:1710.10269, we presented a grid of forward model transmission spectra for hot Jupiter exoplanets. However, we recently identified an error in the treatment of rainout in our 1D atmosphere model ATMO. We explain the error, the correction, its validation and the changed conclusions in this erratum. The correction of this error led to changes in the equilibrium chemical abundances using rainout condensation and thereby the transmission spectra. We note that this error only affects the online library that includes rainout condensation, the library with local condensation (without rainout) is unaffected. The online library was updated with the correction in July 2018.

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.

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.

Observations of Jupiter's gravity field by Juno have revealed surprisingly small 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, much denser in the Mbar region, worsen 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 which satisfy both Juno and Galileo constraints. We confirm that Jupiter's structure must encompass at least four different regions: an outer convective envelope, a region of compositional, thus entropy change, an inner convective envelope and 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 smaller density than for an isentropic profile, 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 details the three physical mechanisms able to 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, combination of the two latter mechanisms seems to be the most favoured solution.

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 {\bf (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, {\bf 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.

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.

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.

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.

We develop an analytical equation of state (EOS) for magnetized fully ionized plasmas, which cover a wide range of temperatures and densities, from low-density classical plasmas to relativistic, quantum plasma conditions. This EOS directly applies to calculations of structure and evolution of strongly magnetized white dwarfs and neutron stars. We review available analytical and numerical results for thermodynamic functions of the nonmagnetized and magnetized Coulomb gases, liquids, and solids. We propose a new analytical expression for the free energy of solid Coulomb mixtures. Based on recent numerical results, we construct analytical approximations for the thermodynamic functions of harmonic Coulomb crystals in quantizing magnetic fields. The analytical description ensures a consistent evaluation of all astrophysically important thermodynamic functions based on the first, second, and mixed derivatives of the free energy. Our numerical code for calculation of thermodynamic functions based on these approximations is made publicly available. Using this code, we calculate and discuss the effects of electron screening and magnetic quantization on the position of the melting point in a range of densities and magnetic fields relevant to white dwarfs and outer envelopes of neutron stars. We consider also the thermal and mechanical structure of a magnetar envelope and argue that it can have a frozen surface which covers the liquid ocean above the solid crust.

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.

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.

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.

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 gravo-turbulent fragmentation theory of Hennebelle & Chabrier 08, 09, 13 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 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 odd 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: while 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 afore-described two-step process.

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.

We review the current status of the theory of thermal emission from the surface layers of neutron stars with strong magnetic fields $B\sim 10^{10}-10^{15}$ 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\sim 10^{10}-10^{11}$ 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 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.

The formation of protoplanetary discs during the collapse of molecular dense cores is significantly influenced by angular momentum transport, notably by the magnetic torque. In turn, the evolution of the magnetic field is determined by dynamical processes and non-ideal MHD effects such as ambipolar diffusion. Considering simple relations between various timescales characteristic of the magnetized collapse, we derive an expression for the early disc radius, $ r \simeq 18 \, {\rm AU} \, \left({\eta_{\rm AD} / 0.1 \, {\rm s}} \right)^{2/9} \left({B_z / 0.1\, {\rm G}} \right) ^{-4/9} \left({M / 0.1 \msol} \right) ^{1/3},$ where $M$ is the total disc plus protostar mass, $\eta_\mathrm{AD}$ is the ambipolar diffusion coefficient and $B_z$ is the magnetic field in the inner part of the core. This is about significantly smaller than the discs that would form if angular momentum was conserved. The analytical predictions are confronted against a large sample of 3D, non-ideal MHD collapse calculations covering variations of a factor 100 in core mass, a factor 10 in the level of turbulence, a factor 5 in rotation, and magnetic mass-to-flux over critical mass-to-flux ratios 2 and 5. The disc radius estimates are found to agree with the numerical simulations within less than a factor 2. A striking prediction of our analysis is the weak dependence of circumstellar disc radii upon the various relevant quantities, suggesting weak variations among class-0 disc sizes. In some cases, we note the onset of large spiral arms beyond this radius.

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.

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.

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.

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

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.

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

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.

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.

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.

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

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 $\alpha\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.

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 $\rho$ 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.

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.

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.

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

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

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.

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.

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.

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

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

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.

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.

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

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

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

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.

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

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

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

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

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

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

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

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.

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.

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

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

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

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

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.

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

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.

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.

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.

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

The determination of Saturn's atmospheric noble gas abundances are critical to understanding the formation and evolution of Saturn, and giant planets in general. These measurements can only be performed with an entry probe. A Saturn probe will address whether enhancement in heavy noble gases, as was found in Jupiter, are a general feature of giant planets, and their ratios will be a powerful constraint on how they form. The helium abundance will show the extent to which helium has phase separated from hydrogen in the planet's deep interior. Jupiter's striking neon depletion may also be tied to its helium depletion, and must be confirmed or refuted in Saturn. Together with Jupiter's measured atmospheric helium abundance, a consistent evolutionary theory for both planets, including "helium rain" will be possible. We will then be able to calibrate the theory of the evolution of all giant planets, including exoplanets. In addition, high pressure H/He mixtures under giant planet conditions are an important area of condensed matter physics that are beyond the realm of experiment.

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.

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.

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

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)

We present an equation of state and radiative opacities for a strongly magnetized hydrogen plasma at magnetic fields B, temperatures T, and densities \rho typical for atmospheres of isolated neutron stars. The first- and second-order thermodynamic functions, monochromatic radiative opacities, and Rosseland mean opacities are calculated and tabulated, taking account of partial ionization, for 8x10^{11} G < B < 3x10^{13} G, 2x10^5 K < T < 10^7 K, and a wide range of \rho. We show that bound-bound and bound-free transitions give an important contribution to the opacities at T < (1-5)x10^6 K in the considered range of B in the outer neutron-star atmosphere layers, which may substantially modify the X-ray spectrum of a typical magnetized neutron star. In addition, we re-evaluate opacities due to free-free transitions, taking into account the motion of both interacting particles, electron and proton, in a strong magnetic field. Compared to the previous neutron-star atmosphere models, the free-free absorption is strongly suppressed at photon frequencies below the proton cyclotron frequency. The latter result holds for any field strength, which prompts a revision of existing theoretical models of X-ray spectra of magnetar atmospheres. Note added in v4: for free-free photoabsorption (Appendix B), equivalent but much simpler expressions are available in Appendix B of Potekhin & Lai, arXiv:astro-ph/0701285.

This White Paper, submitted to the National Academy of Sciences' Astro2010 Decadal Review Committee, focuses on 2 central themes in the study of young brown dwarfs -- their formation mechanism and disk characteristics -- which are of direct relevance to fundamental questions of stellar and planetary origins and properties.

We study the thermal structure and evolution of magnetars as cooling neutron stars with a phenomenological heat source in an internal layer. We focus on the effect of magnetized (B > 10^{14} G) non-accreted and accreted outermost envelopes composed of different elements, from iron to hydrogen or helium. We discuss a combined effect of thermal conduction and neutrino emission in the outer neutron star crust and calculate the cooling of magnetars with a dipole magnetic field for various locations of the heat layer, heat rates and magnetic field strengths. Combined effects of strong magnetic fields and light-element composition simplify the interpretation of magnetars in our model: these effects allow one to interpret observations assuming less extreme (therefore, more realistic) heating. Massive magnetars, with fast neutrino cooling in their cores, can have higher thermal surface luminosity.

We revisit the tidal stability of extrasolar systems harboring a transiting planet and demonstrate that, independently of any tidal model, none but one (HAT-P-2b) of these planets has a tidal equilibrium state, which implies ultimately a collision of these objects with their host star. Consequently, conventional circularization and synchronization timescales cannot be defined because the corresponding states do not represent the endpoint of the tidal evolution. Using numerical simulations of the coupled tidal equations for the spin and orbital parameters of each transiting planetary system, we confirm these predictions and show that the orbital eccentricity and the stellar obliquity do not follow the usually assumed exponential relaxation but instead decrease significantly, reaching eventually a zero value, only during the final runaway merging of the planet with the star. The only characteristic evolution timescale of {\it all} rotational and orbital parameters is the lifetime of the system, which crucially depends on the magnitude of tidal dissipation within the star. These results imply that the nearly circular orbits of transiting planets and the alignment between the stellar spin axis and the planetary orbit are unlikely to be due to tidal dissipation. Other dissipative mechanisms, for instance interactions with the protoplanetary disk, must be invoked to explain these properties.

We develop analytic approximations of thermodynamic functions of fully ionized nonideal electron-ion plasma mixtures. In the regime of strong Coulomb coupling, we use our previously developed analytic approximations for the free energy of one-component plasmas with rigid and polarizable electron background and apply the linear mixing rule (LMR). Other thermodynamic functions are obtained through analytic derivation of this free energy. In order to obtain an analytic approximation for the intermediate coupling and transition to the Debye-Hueckel limit, we perform hypernetted-chain calculations of the free energy, internal energy, and pressure for mixtures of different ion species and introduce a correction to the LMR, which allows a smooth transition from strong to weak Coulomb coupling in agreement with the numerical results.

In this review, we summarize our present knowledge of the behaviour of the mass-radius relationship from solar-type stars down to terrestrial planets, across the regime of substellar objects, brown dwarfs and giant planets. Particular attention is paid to the identification of the main physical properties or mechanisms responsible for this behaviour. Indeed, understanding the mechanical structure of an object provides valuable information about its internal structure, composition and heat content as well as its formation history. Although the general description of these properties is reasonably well mastered, disagreement between theory and observation in certain cases points to some missing physics in our present modelling of at least some of these objects. The mass-radius relationship in the overlaping domain between giant planets and low-mass brown dwarfs is shown to represent a powerful diagnostic to distinguish between these two different populations and shows once again that the present IAU distinction between these two populations at a given mass has no valid foundation.

In the poster presented in Cool Star 15, we analyzed the effect of disk accretion on the evolution of very low mass pre-main sequence stars and young brown dwarfs and the resulting uncertainties on the determination of masses and ages. We use the Lyon evolutionary 1-D code assuming a magnetospheric accretion process, i.e., the material falls covering a small area of the radiative surface, and we take into account the internal energy added from the accreted material as a free parameter $\epsilon$. Even if the approach to this problem is phenomenological, our formalism provides important hints about characteristics of disk accretion, which are useful for improved stellar interior calculations. Using the accretion rates derived from observations our results show that accretion does not affect considerably the position of theoretical isochrones as well as the luminosity compared with standard non-accreting models. See more discussions in a forthcoming paper by Gallardo, Baraffe and Chabrier (2008).

We present the discovery of fifteen new T2.5-T7.5 dwarfs (with estimated distances between ~24-93pc, identified in the first three main data releases of the UKIRT Infrared Deep Sky Survey. This brings the total number of T dwarfs discovered in the Large Area Survey (to date) to 28. These discoveries are confirmed by near infrared spectroscopy, from which we derive spectral types on the unified scheme of Burgasser et al. (2006). Seven of the new T dwarfs have spectral types of T2.5-T4.5, five have spectral types of T5-T5.5, one is a T6.5p, and two are T7-7.5. We assess spectral morphology and colours to identify T dwarfs in our sample that may have non-typical physical properties (by comparison to solar neighbourhood populations). The colours of the full sample of LAS T dwarfs show a possible trend to bluer Y-J with decreasing effective temperature beyond T8. By accounting for the main sources of incompleteness (selection, follow-up and spatial) as well as the effects of unresolved binarity and Malmquist bias, we estimate that there are 17+-4 >=T4 dwarfs in the J<=19 volume of the LAS second data release. Comparing this to theoretical predictions is most consistent with a sub-stellar mass function exponent alpha between -1.0 and 0. This is consistent with the latest 2MASS/SDSS constraint (which is based on lower number statistics), and is significantly lower than the alpha~1.0 suggested by L dwarf field populations, possibly a result of the lower mass range probed by the T dwarf class.

We derive an analytical theory of the prestellar core initial mass function based on an extension of the Press-Schechter statistical formalism. With the same formalism, we also obtain the mass spectrum for the non self-gravitating clumps produced in supersonic flows. The mass spectrum of the self-gravitating cores reproduces very well the observed initial mass function and identifies the different mechanisms responsible for its behaviour. The theory predicts that the shape of the IMF results from two competing contributions, namely a power-law at large scales and an exponential cut-off (lognormal form) centered around the characteristic mass for gravitational collapse. The cut-off exists already in the case of pure thermal collapse, provided that the underlying density field has a lognormal distribution. Whereas pure thermal collapse produces a power-law tail steeper than the Salpeter value, dN/dlog M\propto M^{-x}, with x=1.35, this latter is recovered exactly for the (3D) value of the spectral index of the velocity power spectrum, n\simeq 3.8, found in observations and in numerical simulations of isothermal supersonic turbulence. Indeed, the theory predicts that x=(n+1)/(2n-4) for self-gravitating structures and x=2-n'/3 for non self-gravitating structures, where n' is the power spectrum index of log(rho). We show that, whereas supersonic turbulence promotes the formation of both massive stars and brown dwarfs, it has an overall negative impact on star formation, decreasing the star formation efficiency. This theory provides a novel theoretical foundation to understand the origin of the IMF and to infer its behaviour in different environments. It also provides a complementary approach and useful guidance to numerical simulations exploring star formation, while making testable predictions.

We construct partially ionized hydrogen atmosphere models for magnetized neutron stars in radiative equilibrium with fixed surface fields between B=10^12 and 2x10^13 G and effective temperatures logT_eff=5.5-6.8, as well as with surface B and T_eff distributions around these values. The models are based on the latest equation of state and opacity results for magnetized, partially ionized hydrogen plasmas. The atmospheres directly determine the characteristics of thermal emission from the surface of neutron stars. We also incorporate these model spectra into XSPEC, under the model name NSMAX, thus allowing them to be used by the community to fit X-ray observations of neutron stars.

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.

The development of parallel supercomputers allows today the detailed study of the collapse and the fragmentation of prestellar cores with increasingly accurate numerical simulations. Thanks to the advances in sub-millimeter observations, a wide range of observed initial conditions enable us to study the different modes of low-mass star formation. The challenge for the simulations is to reproduce the observational results. Two main numerical methods, namely AMR and SPH, are widely used to simulate the collapse and the fragmentation of prestellar cores. We compare thoroughly these two methods within their standard framework. We use the AMR code RAMSES and the SPH code DRAGON. Our physical model is as simple as possible and consists of an isothermal sphere rotating around the z-axis. We first study the conservation of angular momentum as a function of the resolution. Then, we explore a wide range of simulation parameters to study the fragmentation of prestellar cores. There seems to be a convergence between the two methods, provided resolution in each case is sufficient. Resolution criteria adapted to our physical cases, in terms of resolution per Jeans mass, for an accurate description of the formation of protostellar cores are deduced from the present study. This convergence is encouraging for future work in simulations of low-mass star formation, providing the aforementioned criteria are fulfilled. Higher resolution figures can be downloaded at http://www-dapnia.cea.fr/Projets/COAST/paper_amrvssph.pdf

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.

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.

The development of parallel supercomputers allows today the detailed study of the collapse and the fragmentation of prestellar cores with increasingly accurate numerical simulations. Thanks to the advances in sub-millimeter observations, a wide range of observed initial conditions enable us to study the different modes of low-mass star formation. The challenge for the simulations is to reproduce the observational results. Two main numerical methods, namely AMR and SPH, are widely used to simulate the collapse and the fragmentation of prestellar cores. We compare here thoroughly these two methods with numerical resolution requirements deduced from previous studies. Our physical model is as simple as possible, and consists of an isothermal sphere rotating around the z-axis. We first study the conservation of angular momentum as a function of the resolution. Then, we explore a wide range of simulation parameters to study the fragmentation of prestellar cores. There seems to be a convergence between the two methods, provided resolution in each case is sufficient. Resolution criteria adapted to our physical cases, in terms of resolution per Jeans mass, for an accurate description of the formation of protostellar cores are deduced from the present study. This convergence is encouraging for future work in simulations of low-mass star formation, providing the aforementioned criteria are fulfilled.

Strong (B >> 10^9 G) and superstrong (B > 10^{14} G) magnetic fields profoundly affect many thermodynamic and kinetic characteristics of dense plasmas in neutron star envelopes. In particular, they produce strongly anisotropic thermal conductivity in the neutron star crust and modify the equation of state and radiative opacities in the atmosphere, which are major ingredients of the cooling theory and spectral atmosphere models. As a result, both the radiation spectrum and the thermal luminosity of a neutron star can be affected by the magnetic field. We briefly review these effects and demonstrate the influence of magnetic field strength on the thermal structure of an isolated neutron star, putting emphasis on the differences brought about by the superstrong fields and high temperatures of magnetars. For the latter objects, it is important to take proper account of a combined effect of the magnetic field on thermal conduction and neutrino emission at densities \rho > 10^{10} g cm^{-3}. We show that the neutrino emission puts a B-dependent upper limit on the effective surface temperature of a cooling neutron star.

(Abridged) 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 the orbital period 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 orbital periods of approximately 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 $\sim 10^4$-$10^6$ yrs, significantly longer than prior estimates. We predict the donor's luminosity and effective temperature. During the adiabatic expansion phase (ignoring irradiation effects), the luminosity is approximately $10^{-6}$--$10^{-4} L_\odot$ and the effective temperature is approximately 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 orbital periods and alters the donor's observational characteristics. Irradiated donors during the adiabatic phase can attain a surface luminosity of up to $\approx10^{-2} L_\odot$. We argue that the turn-on and cooling phases both will leave significant imprints on the AM CVn population's orbital period distribution.

We examine the evolution of low-mass star and brown dwarf eclipsing binaries. These objects are rapid rotators and are believed to shelter large magnetic fields. We suggest that reduced convective efficiency, due to fast rotation and large field strengths, and/or to magnetic spot coverage of the radiating surface significantly affect their evolution, leading to a reduced heat flux and thus larger radii and cooler effective temperatures than for regular objects. We have considered such processes in our evolutionary calculations, using a phenomenological approach. This yields mass-radius and effective temperature-radius relationships in agreement with the observations. We also reproduce the effective temperature ratio and the radii of the two components of the recently discovered puzzling eclipsing brown dwarf system. These calculations show that fast rotation and/or magnetic activity may significantly affect the evolution of eclipsing binaries and that the mechanical and thermal properties of these objects depart from the ones of non-active low-mass objects. We find that, for internal field strengths compatible with the observed surface value of a few kiloGauss, convection can be severely inhibited. The onset of a central radiative zone for rapidly rotating active low-mass stars might thus occur below the usual $\sim 0.35 \msol$ limit.

We explore the possibility that large-scale convection be inhibited over some regions of giant planet interiors, as a consequence of a gradient of composition inherited either from their formation history or from particular events like giant impacts or core erosion during their evolution. Under appropriate circumstances, the redistribution of the gradient of molecular weight can lead to double diffusive layered or overstable convection. This leads to much less efficient heat transport and compositional mixing than large-scale adiabatic convection. We show that this process can explain the abnormally large radius of the transit planet HD209458b and similar objects, and may be at play in some giant planets, with short-period planets offering the most favorable conditions. Observational signatures of this transport mechanism are a large radius and a reduced heat flux output compared with uniformly mixed objects. If our suggestion is correct, it bears major consequences on our understanding of giant planet formation, structure and evolution, including possibly our own jovian planets.

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.

Eccentricity or obliquity tides have been proposed as the missing energy source that may explain the anomalously large radius of some transiting ``hot Jupiters''. To maintain a non-zero and large obliquity, it was argued that the planets can be locked in a Cassini state, i.e. a resonance between spin and orbital precessions. We compute the tidal heating within ``inflated'' close-in giant planets with a non-zero eccentricity or obliquity. We further inspect whether the spin of a ``hot Jupiter'' could have been trapped and maintained in a Cassini state during its early despinning and migration. We estimate the capture probability in a spin-orbit resonance between $\sim$ 0.5 AU (a distance where tidal effects become significant) and 0.05 AU for a wide range of secular orbital frequencies and amplitudes of gravitational perturbations. Numerical simulations of the spin evolution are performed to explore the influence of tidal despinning and migration processes on the resonance stability. We find that tidal heating within a non-synchronous giant planet is about twice larger than previous estimates based on the hypothesis of synchronization. Chances of capture in a spin-orbit resonance are very good around 0.5 AU but they decrease dramatically with the semi-major axis. Furthermore, even if captured, both tidal despinning and migration processes cause the tidal torque to become large enough that the obliquity ultimately leaves the resonance and switches to near $0^{\circ}$. Locking a ``hot Jupiter'' in an isolated spin-orbit resonance is unlikely at 0.05 AU but could be possible at larger distances. Another mechanism is then required to maintain a large obliquity and create internal heating through obliquity tides

We discuss results of our study on AM CVn binaries formed with donors that never ignited He before contact. For the first time, we treat the donor's in these systems in the context of a full stellar structure evolution theory and find that the binary's evolution can described in terms of 3 phases: contact, adiabatic donor expansion, and late-time donor cooling. Details of the first and third phase are new results from this study and we focus on generally characterizing these two phases. Finally, we present our predictions for the donor's light in these systems.

Since the discovery of the first giant planet outside the solar system in 1995 (Mayor & Queloz 1995), more than 180 extrasolar planets have been discovered. With improving detection capabilities, a new class of planets with masses 5-20 times larger than the Earth, at close distance from their parent star is rapidly emerging. Recently, the first system of three Neptune-mass planets has been discovered around the solar type star HD69830 (Lovis et al. 2006). Here, we present and discuss a possible formation scenario for this planetary system based on a consistent coupling between the extended core accretion model and evolutionary models (Alibert et al. 2005a, Baraffe et al. 2004,2006). We show that the innermost planet formed from an embryo having started inside the iceline is composed essentially of a rocky core surrounded by a tiny gaseous envelope. The two outermost planets started their formation beyond the iceline and, as a consequence, accrete a substantial amount of water ice during their formation. We calculate the present day thermodynamical conditions inside these two latter planets and show that they are made of a rocky core surrounded by a shell of fluid water and a gaseous envelope.

We review recent theoretical progress aimed at understanding the formation and the early stages of evolution of giant planets, low-mass stars and brown dwarfs. Calculations coupling giant planet formation, within a modern version of the core accretion model, and subsequent evolution yield consistent determinations of the planet structure and evolution. Because of the uncertainties in the initial conditions, however, it is not possible to say whether young planets are faint or bright compared with low-mass young brown dwarfs. We review the effects of irradiation and evaporation on the evolution of short period planets and argue that substantial mass loss may have occurred for these objects. Concerning star formation, geometrical effects in protostar core collapse are examined by comparing 1D and 3D calculations. Spherical collapse is shown to overestimate the core inner density and temperature and thus to yield incorrect initial conditions for PMS or young brown dwarf evolution. Accretion is also shown to occur over a very limited fraction of the protostar surface. Accretion affects the evolution of young brown dwarfs and yields more compact structures for a given mass and age, thus fainter luminosities. This can lead to severe misinterpretations of the mass and/or age of young accreting objects from their location in the HR diagram. We argue that newborn stars and brown dwarfs should appear rapidly over an extended area in the HR diagram, depending on their accretion history, rather than on a well defined birth line. Finally, we suggest that the distinction between planets and brown dwarfs be based on an observational diagnostic, reflecting the different formation mechanisms between these two distinct populations, rather than on an arbitrary, confusing definition.

This paper presents a consistent description of the formation and the subsequent evolution of gaseous planets, with special attention to short-period, low-mass hot-Neptune planets characteristic of $\mu$ Ara-like systems. We show that core accretion including migration and disk evolution and subsequent evolution taking into account irradiation and evaporation provide a viable formation mechanism for this type of strongly irradiated light planets. At an orbital distance $a \simeq$ 0.1 AU, this revised core accretion model leads to the formation of planets with total masses ranging from $\sim$ 14 $\mearth$ (0.044 $\mjup$) to $\sim$ 400 $\mearth$ (1.25 $\mjup$). The newly born planets have a dense core of $\sim$ 6 $\mearth$, independent of the total mass, and heavy element enrichments in the envelope, $M_{\rm Z,env}/M_{\rm env} $, varying from 10% to 80% from the largest to the smallest planets. We examine the dependence of the evolution of the born planet on the evaporation rate due to the incident XUV stellar flux. In order to reach a $\mu$ Ara-like mass ($\sim$ 14 $\mearth$) after $\sim $ 1 Gyr, the initial planet mass must range from 166 $\mearth$ ($\sim$ 0.52 $\mjup$) to about 20 $\mearth$, for evaporation rates varying by 2 orders of magnitude, corresponding to 90% to 20% mass loss during evolution. The presence of a core and heavy elements in the envelope affects appreciably the structure and the evolution of the planet and yields $\sim 8%-9%$ difference in radius compared to coreless objects of solar composition for Saturn-mass planets. These combinations of evaporation rates and internal compositions translate into different detection probabilities, and thus different statistical distributions for hot-Neptunes and hot-Jupiters.

We study the equation of state, polarization and radiation properties for nonideal, strongly magnetized plasmas which compose outer envelopes of magnetic neutron stars. Detailed calculations are performed for partially ionized hydrogen atmospheres and for condensed hydrogen or iron surfaces of these stars. This is a companion paper to astro-ph/0511803

We briefly examine the properties of dense plasmas characteristic of the interior of giant planets and the atmospheres of neutron stars. Special attention is devoted to the equation of state of hydrogen and helium at high density and to the effect of magnetic fields on the properties of dense matter.

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

Equation of state for partially ionized carbon at temperatures T > ~ 10^5 K is calculated in 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 includes the internal partition functions of bound species. The latter 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.

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.

We study the relation between the mean effective surface temperature T_s and the internal temperature T_b for magnetic neutron stars, assuming that the magnetic field near the surface has a presumably small-scale structure. The heavy-element (iron) and light-element (accreted) heat-blanketing envelopes are considered, and the results are compared with the case of a dipole magnetic field. We argue that the difference in the T_b(T_s)-relation for different magnetic configurations is always much smaller than a possible difference caused by variations of the chemical composition in the envelope.

Source blending in microlensing experiments is known to modify the Einstein time of the observed events. In this paper, we have conducted Monte-Carlo calculations, using the analytical relationships derived by Han (1999) to quantify the effect of blending on the observed event time distribution and optical depth. We show that short-time events are affected significantly by source blending and that, for moderately blended sources, the optical depth $\tau$ is globally overestimated, because of an underestimation of the exposure. For high blending situations, on the opposite, blending leads to an {\it under}estimation of the optical depth. Our results are in agreement with the most recent optical depth determinations toward the Galactic Center of the MACHO collaboration (Popowski et al. 2004) and the OGLE-II collaboration (Sumi et al. 2005) that use clump giants (less affected by the blending effect) as sources. The blending-corrected, lower optical depth toward the Galactic Bulge is now in good agreement with the value inferred from galactic models, reconciling theoretical and observational determinations.

We compare evolutionary models for close-in exoplanets coupling irradiation and evaporation due respectively to the thermal and high energy flux of the parent star with observations of recently discovered new transiting planets. The models provide an overall good agreement with observations, although at the very limit of the quoted error bars of OGLE-TR-10, depending on its age. Using the same general theory, we show that the three recently detected hot-Neptune planets (GJ436, $\rho$ Cancri, $\mu$ Ara) may originate from more massive gas giants which have undergone significant evaporation. We thus suggest that hot-Neptunes and hot-Jupiters may share the same origin and evolution history. Our scenario provides testable predictions in terms of the mass-radius relationships of these hot-Neptunes.

We construct hydrogen atmosphere models for strongly magnetized neutron stars in thermodynamic equilibrium, taking into account partial ionization. The presence of bound states affects the equation of state, absorption coefficients, and polarizability tensor of a strongly magnetized plasma. Therefore the partial ionization influences the polarization vectors and opacities of normal electromagnetic waves, and thus the spectra of outgoing radiation. Here we review a model suitable for the most typical neutron-star atmospheres and focus on the problems that remain to be solved for its extension to other atmospheric parameters.

We present a semi-analytical free-energy model aimed at characterizing the thermodynamic properties of dense fluid helium, from the low-density atomic phase to the high-density fully ionized regime. The model is based on a free-energy minimization method and includes various different contributions representative of the correlations between atomic and ionic species and electrons. This model allows the computation of the thermodynamic properties of dense helium over an extended range of density and temperature and leads to the computation of the phase diagram of dense fluid helium, with its various temperature and pressure ionization contours. One of the predictions of the model is that pressure ionization occurs abruptly at $\rho \simgr 10$ g cm$^{-3}$, {\it i.e.} $P\simgr 20$ Mbar, from atomic helium He to fully ionized helium He$^{2+}$, or at least to a strongly ionized state, without He$^{+}$ stage, except at high enough temperature for temperature ionization to become dominant. These predictions and this phase diagram provide a guide for future dynamical experiments or numerical first-principle calculations aimed at studying the properties of helium at very high density, in particular its metallization. Indeed, the characterization of the helium phase diagram bears important consequences for the thermodynamic, magnetic and transport properties of cool and dense astrophysical objects, among which the solar and the numerous recently discovered extrasolar giant planets.

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

We derive a new formulation to calculate the excess chemical potential of a fraction of $N_1$ particles interacting with $N_2$ particles of a different species. The excess chemical potential is calculated numerically from first principles by coupling molecular dynamics and Thomas-Fermi density functional theory to take into account the contribution arising from the quantum electrons on the forces acting on the ions. The choice of this simple functional is motivated by the fact that the present paper is devoted to the derivation and the validation of the method but more complicated functionals can and will be implemented in the future. This new method is applied in the microcanonical ensemble, the most natural ensemble for molecular dynamics simulations. This avoids the introduction of a thermostat in the simulation, and thus uncontroled modifications of the trajectories calculated from the forces between particles. The calculations are conducted for three values of the input thermodynamic quantities, energy and density, and for different total numbers of particles in order to examine the uncertainties due to finite size effects. This method and these calculations lie the basic foundation to study the thermodynamic stability of dense mixtures, without any {\it a priori} assumption on the degree of ionization of the different species.

Polarizability tensor of a strongly magnetized plasma and the polarization vectors and opacities of normal electromagnetic waves are studied for the conditions typical of neutron star atmospheres, taking account of partial ionization effects. Vacuum polarization is also included using a new set of fitting formulae that are accurate for wide range of field strengths. The full account of the coupling of the quantum mechanical structure of the atoms to their center-of-mass motion across the magnetic field is shown to be crucial for the correct evaluation of the polarization properties and opacities of the plasma. The self-consistent treatment of the polarizability and absorption coefficients proves to be necessary if the ionization degree of the plasma is low, which can occur in the atmospheres of cool or ultramagnetized neutron stars. Atmosphere models and spectra are presented to illustrate the importance of such self-consistent treatment.

We briefly review recent theoretical studies of the effects of accreted envelopes, magnetic fields and crustal superfluidity on the cooling of neutron stars. These effects are especially important for slowly cooling low-mass neutron stars, where direct Urca process of neutrino emission is forbidden. The effects are useful for interpretation of observations of several isolated middle-aged neutron stars.

We determine the expected surface brightness and photometric signature of a white dwarf remnant population, issued from primordial low-mass stars formed at high redshifts, in today galactic halos. We examine the radial dependence of such a contribution as well as its redshift dependence. Such a halo diffuse radiation is below the detection limit of present large field ground-based surveys, but should be observable with the HST and with the future JWST project. Since the surface brightness does not depend on the distance, the integration of several galactic dark halos along the line of sight will raise appreciably the chances of detection. Both the detection or the non-detection of such a remnant diffuse radiation within relevant detection limits offer valuable information on the minimum mass for star formation in the early universe and on the evolution of the stellar initial mass function.

We include the effect of evaporation in our evolutionary calculations of close-in giant planets, based on a recent model for thermal evaporation taking into account the XUV flux of the parent star (Lammer et al. 2003). Our analysis leads to the existence of a critical mass for a given orbital distance $m_{\rm crit}(a)$ below which the evaporation timescale becomes shorter than the thermal timescale of the planet. For planets with initial masses below $m_{\rm crit}$, evaporation leads to a rapid expansion of the outer layers and of the total planetary radius, speeding up the evaporation process. Consequently, the planet does not survive as long as estimated by a simple application of mass loss rates without following consistently its evolution. We find out that the transit planet HD 209458b might be in such a dramatic phase, although with an extremely small probability. As a consequence, we predict that, after a certain time, only planets above a value $m_{\rm crit}(a)$ should be present at an orbital distance $a$ of a star. For planets with initial masses above $m_{\rm crit}$, evaporation does not affect the evolution of the radius with time.

We review our present understanding of the physical properties of substellar objects, brown dwarfs and irradiated or non-irradiated gaseous exoplanets. This includes a description of their internal properties, mechanical structure and heat content, their atmospheric properties, thermal profile and emergent spectrum, and their evolution, in particular as irradiated companions of a close parent star. The general theory can be used to make predictions in term of detectability for the future observational projects. Special attention is devoted to the evolution of the two presently detected transit planets, HD209458B and OGLE-TR-56B. For this latter, we present a consistent evolution for its recently revised mass and show that we reproduce the observed radius within its error bars. We briefly discuss differences between brown dwarfs and gaseous planets, both in terms of mass function and formation process. We outline several arguments to show that the minimum mass for deuterium burning, recently adopted officially as the limit to distinguish the two types of objects, is unlikely to play any specific role in star formation, so that such a limit is of purely semantic nature and is not supported by a physical justification.

Extending the theory we derived recently for HD209458b to different cases of strongly irradiated gaseous exoplanets, we have calculated the consistent evolution of the new transiting planet, OGLE-TR-56b, for its recently revised mass determination. The theory is shown to successfully reproduce the observed radius, for the proper age of the system. We also examine the dissipation of kinetic energy at the planet's internal adiabat due to atmospheric winds, and place constraints on the efficiency of this process. We show that a fraction $\sim 0.1$ -- 0.5% of the incident flux transformed into thermal energy deposited at the adiabatic level can accommodate the observed radii of both OGLE-TR-56b and HD209458b. The present theory yields quantitative predictions on the evolution of the emergent spectrum and fundamental properties of hot-jupiters. The predictions for radius, luminosity, temperature as a function of the planet's mass and orbital distance can be used as benchmarks for future detections of transit planets.

We calculate the thermal structure and quiescent thermal luminosity of accreting neutron stars (warmed by deep crustal heating in accreted matter) in soft X-ray transients (SXTs). We consider neutron stars with nucleon and hyperon cores and with accreted envelopes. It is assumed that an envelope has an outer helium layer (of variable depth) and deeper layers of heavier elements, either with iron or with much heavier nuclei (of atomic weight A > 100) on the top (Haensel & Zdunik 1990, 2003, astro-ph/0305220). The relation between the internal and surface stellar temperatures is obtained and fitted. The quiescent luminosity of the hottest (low-mass) and coldest (high-mass) neutron stars is calculated, together with the ranges of its possible variations due to variable thickness of the helium layer. The results are compared with observations of SXTs, particularly, containing the coldest (SAX J1808.4-3658) and the hottest (Aql X-1) neutron stars. The observations of SAX J1808.4-3658 in a quiescent state on March 24, 2001 (Campana et al. 2002, astro-ph/0206376) can be explained only if this SXT contains a massive neutron star with a nucleon/hyperon core; a hyperon core with a not too low fraction of electrons is preferable. Future observations may discriminate between the various models of hyperon/nucleon dense matter. The thermal emission of SAX J1808.4-3658 is also sensitive to the models of plasma ionization in the outermost surface layers and can serve for testing such models.

We study the thermal structure of neutron stars with magnetized envelopes composed of accreted material, using updated thermal conductivities of plasmas in quantizing magnetic fields, as well as equation of state and radiative opacities for partially ionized hydrogen in strong magnetic fields. The relation between the internal and local surface temperatures is calculated and fitted by an analytic function of the internal temperature, magnetic field strength, angle between the field lines and the normal to the surface, surface gravity, and the mass of the accreted material. The luminosity of a neutron star with a dipole magnetic field is calculated for various values of the accreted mass, internal temperature, and magnetic field strength. Using these results, we simulate cooling of superfluid neutron stars with magnetized accreted envelopes. We consider slow and fast cooling regimes, paying special attention to very slow cooling of low-mass superfluid neutron stars. In the latter case, the cooling is strongly affected by the combined effect of magnetized accreted envelopes and neutron superfluidity in the stellar crust. Our results are important for interpretation of observations of isolated neutron stars hottest for their age, such as RX J0822-43 and PSR B1055-52.

The equation of state and radiative opacities of partially ionized, strongly magnetized hydrogen plasmas, presented in a previous paper [ApJ 585, 955 (2003), astro-ph/0212062] for the magnetic field strengths 8.e11 G < B < 3.e13 G, are extended to the field strengths 3.e13 G < B < 1.e15 G, relevant for magnetars. The first- and second-order thermodynamic functions and radiative opacities are calculated and tabulated for 5.e5 < T < 4.e7 K in a wide range of densities. We show that bound-free transitions give an important contribution to the opacities in the considered range of B in the outer neutron-star atmosphere layers. Unlike the case of weaker fields, bound-bound transitions are unimportant.

We construct partially ionized hydrogen atmosphere models for magnetized neutron stars in radiative equilibrium with surface fields B=10^12-5 \times 10^14 G and effective temperatures T_eff \sim a few \times 10^5-10^6 K. These models are based on the latest equation of state and opacity results for magnetized, partially ionized hydrogen plasmas that take into account various magnetic and dense medium effects. The atmospheres directly determine the characteristics of thermal emission from isolated neutron stars. For the models with B=10^12-10^13 G, the spectral features due to neutral atoms lie at extreme UV and very soft X-ray energy bands and therefore are difficult to observe. However, the continuum flux is also different from the fully ionized case, especially at lower energies. For the superstrong field models (B\ga 10^14 G), we show that the vacuum polarization effect not only suppresses the proton cyclotron line as shown previously, but also suppresses spectral features due to bound species; therefore spectral lines or features in thermal radiation are more difficult to observe when the neutron star magnetic field is \ga 10^14 G.

We review recent determinations of the present day and initial mass functions in various components of the Galaxy, disk, spheroid, young and globular clusters. As a general feature, the IMF is well described by a power-law form for $m\ga 1 \msol$ and a lognormal form below. The extension of the disk IMF into the brown dwarf (BD) regime is in good agreement with observations and yields a disk BD number-density comparable to the stellar one $\sim 0.1 \pc3$. The IMF of young clusters is found to be consistent with the disk field IMF, providing the same correction for unresolved binaries. The spheroid IMF relies on much less robust grounds. Within all the uncertainties, it is found to be similar to the one derived for globular clusters, and is well represented also by a lognormal form with a characteristic mass slightly larger than for the disk. The IMF characteristic of early star formation remains undetermined, but different observational constraints suggest that it does not extend below $\sim 1 \msol$. These IMFs allow a reasonably robust determination of the Galactic present-day and initial stellar and brown dwarf contents. They also have important galactic implications in yielding more accurate mass-to-light ratio determinations. The M/L ratios obtained with the disk and the spheroid IMF yield values 1.8 and 1.4 smaller than a Salpeter IMF, respectively. This general IMF determination is examined in the context of star formation theory. (shortened)

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

We present evolutionary models for cool brown dwarfs and extra-solar giant planets. The models reproduce the main trends of observed methane dwarfs in near-IR color-magnitude diagrams. We also present evolutionary models for irradiated planets, coupling for the first time irradiated atmosphere profiles and inner structures. We focus on HD 209458-like systems and show that irradiation effects can substantially affect the radius of sub-jovian mass giant planets. Irradiation effects, however, cannot alone explain the large observed radius of HD 209458b. Adopting assumptions which optimise irradiation effects and taking into account the extension of the outer atmospheric layers, we still find $\sim$ 20% discrepancy between observed and theoretical radii. An extra source of energy seems to be required to explain the observed value of the first transit planet.

We construct hydrogen atmosphere models for magnetized neutron stars in radiative equilibrium with surface fields B=10^12-5x10^14 G and effective temperatures T_eff a few x 10^5-10^6 K by solving the full radiative transfer equations for both polarization modes in the magnetized hydrogen plasma. The atmospheres directly determine the characteristics of thermal emission from isolated neutron stars. We study the effects of vacuum polarization and bound atoms on the atmosphere structure and spectra. For the lower magnetic field models (B 10^12 G), the spectral features due to neutral atoms lie at extreme UV and very soft X-ray energies and therefore are not likely to be observed. However, the continuum flux is also different from the fully ionized case, especially at lower energies. For the higher magnetic field models, we find that vacuum polarization softens the high energy tail of the thermal spectrum. We show that this depression of continuum flux strongly suppresses not only the proton cyclotron line but also spectral features due to bound species; therefore spectral lines or features in thermal radiation are more difficult to observe when the neutron star magnetic field is > 10^14 G.

We briefly examine the properties of dense plasmas characteristic of the atmospheres of neutron stars and of the interior of massive white dwarfs. These astrophysical bodies are natural laboratories to study respectively the problem of pressure ionization of hydrogen in a strong magnetic field and the crystallization of the quantum one-component-plasma at finite temperature.

We study the thermodynamic properties of a partially ionized hydrogen plasma in strong magnetic fields, B ~ 10^{12}-10^{13} G, typical of neutron stars. The properties of the plasma depend significantly on the quantum-mechanical sizes and binding energies of the atoms, which are strongly modified by thermal motion across the field. We use new fitting formulas for the atomic binding energies and sizes, based on accurate numerical calculations and valid for any state of motion of the atom. In particular, we take into account decentered atomic states, neglected in previous studies of thermodynamics of magnetized plasmas. We also employ analytic fits for the thermodynamic functions of nonideal fully ionized electron-ion Coulomb plasmas. This enables us to construct an analytic model of the free energy. An ionization equilibrium equation is derived, taking into account the strong magnetic field effects and the nonideality effects. This equation is solved by an iteration technique. Ionization degrees, occupancies, and the equation of state are calculated.

The analytic equation of state of nonideal Coulomb plasmas consisting of pointlike ions immersed in a polarizable electron background (physics/9807042) is improved, and its applicability range is considerably extended. First, the fit of the electron screening contribution in the free energy of the Coulomb liquid is refined at high densities where the electrons are relativistic. Second, we calculate the screening contribution for the Coulomb solid (bcc and fcc) and derive an analytic fitting expression. Third, we propose a simple approximation to the internal and free energy of the liquid one-component plasma of ions, accurate within the numerical errors of the most recent Monte Carlo simulations. We obtain an updated value of the coupling parameter at the solid-liquid phase transition for the one-component plasma: Gamma_m = 175.0 (+/- 0.4).

Recent observations show that chromospheric activity in late-M and L dwarfs is much lower than in the earlier M types, in spite of comparatively rapid rotation. We investigate the possibility that this drop-off in activity results from the very high electrical resistivities in the dense, cool and predominantly neutral atmospheres of late-M and L dwarfs. We calculate magnetic field diffusivities in the atmospheres of objects with effective temperatures in the range 3000-1500 (mid-M to L), using the atmospheric structure models of Allard and Hauschildt. We find that the combination of very low ionization fraction and high density in these atmospheres results in very large resistivities due to neutral-charged particle collisions, and efficient field diffusion. The resistivities are found to increase with both decreasing optical depth, and decreasing effective temperature. As a result, any existing magnetic fields are increasingly decoupled from atmospheric motions as one moves from mid-M to L; we quantify this through a simple Reynolds number calculation. This, coupled with the difficulty in transporting magnetic stresses through the highly resistive atmosphere, can account for the observed drop in activity from mid-M to L, assuming activity in these objects is magnetically driven. We also examine the issue of acoustic heating, and find that this appears inadequate to explain the observed H-alpha fluxes in mid-M to L dwarfs. Consequently, magnetic heating does seem to be the most viable mechanism for generating activity in these objects. Finally, we speculate on a possible flare mechanism in these cool dwarfs.

The low-frequency electric microfield distribution in a Coulomb plasma is calculated for various plasma parameters, from weak to strong Coulomb coupling and from zero to strong electron screening. Two methods of numerical calculations are employed: the adjustable-parameter exponential approximation and the Monte Carlo simulation. The results are represented by analytic fitting formulas suitable for applications.

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. $\rho_{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 $\rho_\odot \approx 5.0\times 10^{-2} \mvol$.

We analyse pre-Main Sequence evolutionary tracks for low mass stars with masses $m \le 1.4 \msol$ based on the Baraffe et al. (1998) input physics. We also extend the recent Chabrier et al. (2000) evolutionary models based on dusty atmosphere to young brown dwarfs down to one mass of Jupiter. We analyse current theoretical uncertainties due to molecular line lists, convection and initial conditions. Simple tests on initial conditions show the high uncertainties of models at ages $\simle$ 1 Myr. We find a significant sensitivity of atmosphere profiles to the treatment of convection at low gravity and $\te < 4000$ K, whereas it vanishes as gravity increases. This effect adds another source of uncertainty on evolutionary tracks at very early phases. We show that at low surface gravity ($\log g \simle 3.5$,) the common picture of vertical Hayashi lines with constant $\te$ is oversimplified. The effect of a variation of initial deuterium abundance is studied. We compare our models with evolutionary tracks available in the literature and discuss the main differences. We finally analyse to which extent current observations of young systems provide a good test for pre-Main Sequence tracks.

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.

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.

We present new spectroscopic and photometric measurements of the white dwarfs LHS 3250 and WD 0346+246. Along with F351-50, these white dwarfs are the coolest ones known, all with effective temperatures below 4000 K. Their membership in the Galactic halo population is discussed, and detailed comparisons of all three objects with new atmosphere models are presented. The new models consider the effects of mixed H/He atmospheres and indicate that WD 0346+246 and F351-50 have predominantly helium atmospheres with only traces of hydrogen. LHS 3250 may be a double degenerate whose average radiative temperature is between 2000 and 4000 K, but the new models fail to explain this object.

Outer envelopes of neutron stars consist mostly of fully ionized, strongly coupled Coulomb plasmas characterized by typical densities about 10^4-10^{11} g/cc and temperatures about 10^4-10^9 K. Many neutron stars possess magnetic fields about 10^{11}-10^{14} G. Here we briefly review recent theoretical advances which allow one to calculate thermodynamic functions and electron transport coefficients for such plasmas with an accuracy required for theoretical interpretation of observations.

We consider the depletion of primordial deuterium in the interior of substellar objects as a function of mass, age and absolute magnitude in several photometric passbands. We characterize potential spectroscopic signatures of deuterium in the lines of deuterated water HDO. These results will serve as a useful, independent diagnostic to characterize the mass and/or the age of young substellar objects, and to provide an independent age determination of very young clusters. These results can serve to identify objects at the deuterium-burning limit and to confront the theoretical prediction that D-burning is a necessary condition to form star-like objects.

We present evolutionary models for low mass stars and brown dwarfs ($m \le 1.2 \msol$) based on recent improvement of the theory: equation of state, atmosphere models, ... We concentrate on early evolutionary phases from the initial deuterium burning phase to the zero-age Main Sequence. Evolutionary models for young brown dwarfs are also presented. We discuss the uncertainties of the present models. We analyse the difficulties arising when comparing models with observations for very young objects, in particular concerning the problem of reddening.

Since the discovery of the first bona-fide brown dwarfs and extra-solar planets in 1995, the field of low mass stars and substellar objects has considerably progressed, both from theoretical and observational viewpoints.Recent developments in the physics entering the modeling of these objects have led to significant improvements in the theory and to a better understanding of their mechanical and thermal properties. This theory can now be confronted with observations directly in various observational diagrams (color-color, color-magnitude, mass-magnitude, mass-spectral type), a stringent and unavoidable constraint which became possible only recently, with the generation of synthetic spectra. In this paper, we present the current state-of-the-art general theory of low-mass stars and sub-stellar objects, from one solar mass to one Jupiter mass, regarding primarily their interior structure and evolution. This review is a natural complement to the previous review on the atmosphere of low-mass stars and brown dwarfs (Allard et al 1997). Special attention is devoted to the comparison of the theory with various available observations. The contribution of low-mass stellar and sub-stellar objects to the Galactic mass budget is also analysed.

We present new cooling sequences, color-magnitude diagrams, and color-color diagrams for cool white dwarfs with pure hydrogen atmospheres down to an effective temperature $\te=1500$ K. We include a more detailed treatment of the physics of the fully-ionized interior, particularly an improved discussion of the thermodynamics of the temperature-dependent ion-ion and ion-electron contributions of the quantum, relativistic electron-ion plasma. The present calculations also incorporate accurate boundary conditions between the degenerate core and the outermost layers as well as updated atmosphere models including the H$_2$-H$_2$ induced-dipole absorption. We examine the differences on the cooling time of the star arising from uncertainties in the initial carbon-oxygen profile and the core-envelope $L$-$T_c$ relation. The maximum time delay due to crystallization-induced chemical fractionation remains substantial, from $\sim 1.0$ Gyr for 0.5 and 1.2 $\msol$ white dwarfs to $\sim 1.5$ Gyr for 0.6 to 0.8 $\msol$ white dwarfs, even with initial stratified composition profiles, and cannot be ignored in detailed white dwarf cooling calculations. These cooling sequences provide theoretical support to search for or identify old disk or halo hydrogen-rich white dwarfs by characterizing their mass and age from their observational signatures.

We present evolutionary calculations for very-low-mass stars and brown dwarfs based on synthetic spectra and non-grey atmosphere models which include dust formation and opacity, i.e. objects with $\te\simle 2800$ K. The interior of the most massive brown dwarfs is shown to develop a conductive core after $\sim 2$ Gyr which slows down their cooling. Comparison is made in optical and infrared color-magnitude diagrams with recent late-M and L-dwarf observations. The saturation in optical colors and the very red near-infrared colors of these objects are well explained by the onset of dust formation in the atmosphere. Comparison of the faintest presently observed L-dwarfs with these dusty evolutionary models suggests that dynamical processes such as turbulent diffusion and gravitational settling are taking place near the photosphere. As the effective temperature decreases below $\te\approx 1300-1400$ K, the colors of these objects move to very blue near-infrared colors, a consequence of the ongoing methane absorption in the infrared. We suggest the possibility ofa brown dwarf dearth in $J,H,K$ color-magnitude diagrams around this temperature.

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.

The evolution of white dwarfs is a cooling process that depends on the energy stored in the core and on the way in which it is transferred through the envelope. In this paper we show that despite some (erroneous) claims, the redistribution of chemical elements ensuing the crystallization of C/O white dwarfs provides between the 10% and the 20% of the total energy released during the crystallization process, depending on the internal chemical composition, which is not negligible at all, given the present state of the art of the white dwarf cooling theory.

We briefly summarize our present knowledge of the theory of low-mass stars and substellar objects and their contribution to the Galactic population.

We briefly outline the physics underlying the mechanical and thermal properties of brown dwarfs, which characterizes their interiors and their atmospheres. We mention the most recent improvements realyzed in the theory of brown dwarfs and the connection with experimental and observational tests of this theory.

We calculate the expected white dwarf luminosity functions and discovery functions in photometric passbands, if these stellar remnants provide a substantial fraction of the sought Galactic dark matter, as suggested on various observational grounds. We demonstrate the extremely rapid variation of the photometic signature of halo white dwarfs with time and thus the powerful diagnostic of white dwarf colors to determine the age of the Galactic halo. We also consider the various indirect constraints implied by a white dwarf dominated halo. These calculations will guide present and future observational projects at faint magnitudes. This will enable us to determine not only the nature of the Galactic dark matter but also the age and the initial conditions of the Galaxy formation.

We study hydrogen plasmas at magnetic fields B ~ 10^{12}-10^{13} Gauss, densities ~ 10^{-3}-10^3 g/cc and temperatures T ~ 10^{5.5}-10^{6.5} K, typical of photospheres of middle-aged cooling neutron stars. We construct an analytical free energy model of the partially ionized plasma, including into consideration the decentred atomic states, which arise due to the thermal motion across the strong field. We show that these states, neglected in previous studies, may contribute appreciably into thermodynamics of the outer atmospheric layers at density below 1 g/cc and typical B and T. We take into account Coulomb non-ideality of the ionized component of the plasma affected by intense magnetic field. Ionization degree, occupancies and equation of state are calculated, and their dependences on the temperature, density and magnetic field are studied.

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.

We present evolutionary models for low mass stars from 0.075 to 1 $\msol$ for solar-type metallicities [M/H]= 0 and -0.5. The calculations include the most recent interior physics and the latest generation of non-grey atmosphere models. We provide mass-age-color-magnitude relationships for both metallicities. The mass-M$_V$ and mass-M$_K$ relations are in excellent agreement with the empirical relations derived observationally. The theoretical color-magnitude diagrams are compared with the sequences of globular clusters (47 Tucanae) and open clusters (NGC2420 and NGC2477) observed with the Hubble Space Telescope. Comparison is also made with field star sequences in $M_V$-$(V-I)$, $M_K$-$(I-K)$ and $M_K$-$(J-K)$ diagrams. These comparisons show that the most recent improvements performed in low-mass star atmosphere models yield now reliable stellar models in the near-infrared. These models can be used for metallicity, mass, temperature and luminosity calibrations. Uncertainties still remain, however, in the optical spectral region below $T_{eff} \sim 3700K$, where predicted (V-I) colors are too blue by 0.5 mag for a given magnitude. The possible origins for such a discrepancy, most likely a missing source of opacity in the optical and the onset of grain formation are examined in detail.

We have computed stellar evolutionary models for stars in a mass range characteristic of Cepheid variables ($3<m/\Msol<12$) for different metallicities representative of the Galaxy and the Magellanic Clouds populations. The stellar evolution calculations are coupled to a linear non adiabatic stability analysis to get self-consistent mass-period-luminosity relations. The period - luminosity relation as a function of metallicity is analysed and compared to the recent EROS observations in the Magellanic Clouds. The models reproduce the observed width of the instability strips for the SMC and LMC. We determine a statistical P-L relationship, taking into account the evolutionary timescales and a mass distribution given by a Salpeter mass function. Excellent agreement is found with the SMC PL relationship determined by Sasselov et al. (1997). The models reproduce the change of slope in the P-L relationship near $P\sim 2.5$ days discovered recently by the EROS collaboration (Bauer 1997; Bauer et al. 1998) and thus explain this feature in term of stellar evolution. Some discrepancy, however, remains for the LMC Cepheids. The models are also in good agreement with Beat Cepheids observed by the MACHO and EROS collaborations. We show that most of the 1H/2H Beat Cepheids have not yet ignited central helium burning; they are just evolving off the Main Sequence toward the red giant branch.

In dense stars the nuclear reaction rates are influenced by screening effects arising from both ions and electrons. In this paper we calculate the enhancement factors due to electron polarization in the high-density, degenerate and relativistic regime, for non-resonant nuclear reaction rates. In an earlier analysis, Sahrling had proposed the possibility that the polarized electrons would lower the reaction rate instead of enhancing it. This analysis was based on Monte Carlo simulations with only one choice of density, temperature and charge. Here we extend the analysis to a wider range of densities, temperatures and charges and calculate analytical expressions for the enhancement factors. We concentrate on carbon and oxygen ions and show that at very high-densities, high-order quantum effects will be important and act to reduce the zeroth order, classical value for the enhancement factor. We show that in any case, the total electron contribution remains weak, namely an enhancement in the reaction of about a factor 2, contrarily to what had been claimed by some authors in previous calculations. We examine the astrophysical implications of these results on the final stages of massive white dwarfs, near the carbon-ignition curve.

We use the calculations derived in a previous paper (M\'era, Chabrier and Schaeffer, 1997), based on observational constraints arising from star counts, microlensing experiments and kinematic properties, to determine the amount of dark matter under the form of stellar and sub-stellar objects in the different parts of the Galaxy. This yields the derivation of different mass-models for the Galaxy. In the light of all the afore-mentioned constraints, we discuss two models that correspond to different conclusions about the nature and the location of the Galactic dark matter. In the first model there is a small amount of dark matter in the disk, and a large fraction of the dark matter in the halo is still undetected and likely to be non-baryonic. The second, less conventional model is consistent with entirely, or at least predominantly baryonic dark matter, under the form of brown dwarfs in the disk and white dwarfs in the dark halo. We derive observational predictions for these two models which should be verifiable by near future infrared and microlensing observations.

We examine the most recent observational constraints arising from i) small-scale and large-scale Galactic dynamical properties, ii) star counts at faint magnitude and iii) microlensing experiments. From these constraints, we determine the halo and disk stellar mass functions and stellar content down to the bottom of the main sequence, which yields the normalization of the halo/disk total stellar population, and we infer the contributions of sub-stellar objects to the mass budget of the various Galactic regions. The consistent analysis of star counts and of the overall microlensing observations in the Bulge are compatible with a small contribution of brown dwarfs to the Galactic mass budget $\rho_{BD}/\rho_* \leq 0.2 $. However the separate bulge/disk analysis based on the bulge clump giants is compatible with a substantial population of disk brown dwarfs, $\Sigma_{BD}/\Sigma_*\leq 1 $. More statistics of microlensing events towards the Galactic center and a better determination of the velocity dispersions in the bulge should break this degeneracy of solutions. For the halo, we show that a steep mass-function in the dark halo is excluded and that low-mass stars and brown dwarfs represent a negligible fraction of the halo dark matter, and thus of the observed events towards the LMC. The nature of these events remains a puzzle and halo white dwarfs remain the least unlikely candidates.

The relationships between the effective surface (T_e) and internal temperatures of neutron stars (NSs) with and without accreted envelopes are calculated for T_e > 5\times10^4 K using new data on the equation of state and opacities in the outer NS layers. We examine various models of accreted layers (H, He, C, O shells produced by nuclear transformations in accreted matter). We employ new Opacity Library (OPAL) radiative opacities for H, He, and Fe. In the outermost NS layers, we implement the modern OPAL equation of state for Fe, and the Saumon-Chabrier equation of state for H and He. The updated thermal conductivities of degenerate electrons include the Debye-Waller factor for the electron-phonon scattering in solidified matter, while in liquid matter they include the contributions from electron-ion collisions (evaluated with non-Born corrections and with the ion structure factors in responsive electron background) and from the electron-electron collisions. For T_e < 10^{5.5} K, the electron conduction in non-degenerate layers of the envelope becomes important, reducing noticeably the temperature gradient. The accreted matter further decreases this gradient at T_e > 10^5 K. Even a small amount of accreted matter (with mass > 10^{-16} Msun) affects appreciably the NS cooling, leading to higher T_e at the neutrino cooling stage and to lower T_e at the subsequent photon stage.

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

We present a review of the most recent improvements in the theory of low-mass star-like objects, in the stellar and sub-stellar domains. Comparison is made with various recent observations.

We summarize recent improvements in model atmosphere and internal structure of faint white dwarfs. We derive an analytical cooling theory which illustrates the effects of various physical processes on the energy budget and the cooling history of cool white dwarfs. We consider in particular the effect of chemical fractionation at crystallization. This process, although it liberates a small amount of energy, affects significantly the cooling time of white dwarfs at their low-luminosity stages. This modifies the age of the faintest white dwarf ever observed, ESO 439-26, and bears important consequences for the Galactic implications of white dwarf evolution, age of the Galactic disk and expected contribution to the halo mass. We briefly consider the remaining major uncertainties in white dwarf cooling theory.

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

We present extensive calculations of the structure and the evolution of low-massstars in the range 0.07-0.8 $\msol$, for metallicities $-2.0\le \mh \le 0.0$. These calculations are based on the most recent description of the microphysics characteristic of these dense and cool objects and on the lattest generation of grainless non-grey atmosphere models. We examine the evolution of the different mechanical and thermal properties of these objects as a function of mass and metallicity. We also demonstrate the inaccuracy of grey models and $T(\tau)$ relationships under these conditions. We provide detailed tables of the mass-radius-luminosity-effective temperature relations for various ages and metallicities, aimed at calibrating existing or future observations of low-mass stars and massive brown dwarfs. We derive new hydrogen-burning minimum masses, within the afore-mentioned metallicity range. These minimum masses are found to be smaller than previous estimates, a direct consequence of non-grey effects. At last, we examine the evolution of the abundance of light elements, $Li, Be$ and $B$, as a function of age, mass and metallicity.

We examine the recent results of the MACHO collaboration towards the Large Magellanic Cloud (Alcock et al. 1996) in terms of a halo brown dwarf or white dwarf population. The possibility for most of the microlensing events to be due to brown dwarfs is totally excluded by large-scale kinematic properties. The white dwarf scenario is examined in details in the context of the most recent white dwarf cooling theory (Segretain et al. 1994) which includes explicitely the extra source of energy due to carbon-oxygen differentiation at crystallization, and the subsequent Debye cooling. We show that the observational constraints arising from the luminosity function of high-velocity white dwarfs in the solar neighborhood and from the recent HST deep field counts are consistent with a white dwarf contribution to the halo missing mass as large as 50 %, provided i) an IMF strongly peaked around 1.7 Msol and ii) a halo age older than 18 Gyr.

We present preliminary non-grey model atmospheres and interiors for cool brown dwarfs. The resulting synthetic spectra are compared to available spectroscopic and photometric observations of the coolest brown dwarf yet discovered, Gl229B (Nakajima \etal, 1995). Despite recognized shortcomings of current methane opacities, we find that the model spectra provide an excellent fit to most of the spectral features of the brown dwarf. We confirm the presence of methane absorption and the substellar nature of Gl229B. These preliminary models set an upper limit for the effective temperature of 1000~K. We also compute the evolution of brown dwarfs with solar composition and masses from 0.02 to 0.065 $\msol$. While uncertainties in the age of the system yield some undetermination for the mass of Gl229B, the most likely solution is $m\approx 0.04-0.055 \msol$. In any case, we can set an upper limit $m= 0.065\msol$ for a very unlikely age $t=10$ Gyr.

We present a broad suite of models of extrasolar giant planets (EGP's), ranging in mass from 0.3 to 15 Jupiter masses. The models predict luminosity (both reflected and emitted) as a function of age, mass, deuterium abundance and distance from parent stars of various spectral type. We also explore the effects of helium mass fraction, rotation rate and the presence of a rock-ice core. The models incorporate the most accurate available equation of state for the interior, including a new theory for the enhancement of deuterium fusion by electron screening which is potentially important in these low mass objects. The results of our calculations reveal the enormous sensitivity of EGP's to the presence of the parent star, particularly for G and earlier spectral types. They also show a strong sensitivity of the flux contrast in the mid-infrared between parent star and EGP to the mass and age of the EGP's. We interpret our results in terms of search strategies for ground- and space-based observatories in place or anticipated in the near future.