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.

Extending our recent studies of two-dimensional stellar convection to 3D, we compare three-dimensional hydrodynamic simulations to identically set-up two-dimensional simulations, for a realistic pre-main sequence star. We compare statistical quantities related to convective flows including: average velocity, vorticity, local enstrophy, and penetration depth beneath a convection zone. These statistics are produced during stationary, steady-state compressible convection in the star's convection zone. Our simulations with the MUSIC code confirm the common result that two-dimensional simulations of stellar convection have a higher magnitude of velocity on average than three-dimensional simulations. Boundary conditions and the extent of the spherical shell can affect the magnitude and variability of convective velocities. The difference between 2D and 3D velocities is dependent on these background points; in our simulations this can have an effect as large as the difference resulting from the dimensionality of the simulation. Nevertheless, radial velocities near the convective boundary are comparable in our 2D and 3D simulations. The average local enstrophy of the flow is lower for two-dimensional simulations than for three-dimensional simulations, indicating a different shape and structuring of 3D stellar convection. We perform a statistical analysis of the depth of convective penetration below the convection zone, using the model proposed in our recent study (Pratt et al. 2017). Here we analyze the convective penetration in three dimensional simulations, and compare the results to identically set-up 2D simulations. In 3D the penetration depth is as large as the penetration depth calculated from 2D simulations.

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.

Context. Atmospheric superrotating flows at the equator are an almost ubiquitous result of simulations of hot Jupiters, and a theory explaining how this zonally coherent flow reaches an equilibrium has been developed in the literature. However, this understanding relies on the existence of either an initial superrotating or a sheared flow, coupled with a slow evolution such that a linear steady state can be reached. Aims. A consistent physical understanding of superrotation is needed for arbitrary drag and radiative timescales, and the relevance of considering linear steady states needs to be assessed. Methods. We obtain an analytical expression for the structure, frequency and decay rate of propagating waves in hot Jupiter atmospheres around a state at rest in the 2D shallow water beta plane limit. We solve this expression numerically and confirm the robustness of our results with a 3D linear wave algorithm. We then compare with 3D simulations of hot Jupiter atmospheres and study the non linear momentum fluxes. Results. We show that under strong day night heating the dynamics does not transit through a linear steady state when starting from an initial atmosphere in solid body rotation. We further show that non linear effects favour the initial spin up of superrotation and that the acceleration due to the vertical component of the eddy momentum flux is critical to the initial development of superrotation. Conclusions. Overall, we describe the initial phases of the acceleration of superrotation, including consideration of differing radiative and drag timescales, and conclude that eddy-momentum driven superrotating equatorial jets are robust, physical phenomena in simulations of hot Jupiter atmospheres.

Context. The study of linear waves and instabilities is necessary to understand the physical evolution of an atmosphere, and can provide physical interpretation of the complex flows found in simulations performed using Global Circulation Models (GCM). In particular, the acceleration of superrotating flow at the equator of hot Jupiters has mostly been studied under several simplifying assumptions, the relaxing of which may impact final results. Aims. We develop and benchmark a publicly available algorithm to identify the eigenmodes of an atmosphere around any initial steady state. We also solve for linear steady states. Methods. We linearise the hydrodynamical equations of a planetary atmosphere in a steady state with arbitrary velocities and thermal profile. We then discretise the linearised equations on an appropriate staggered grid, and solve for eigenvectors and linear steady solutions with the use of a parallel library for linear algebra: ScaLAPACK. We also implement a posteriori calculation of an energy equation in order to obtain more information on the underlying physics of the mode. Results. Our code is benchmarked against classical wave and instability test cases in multiple geometries. The steady linear circulation calculations also reproduce expected results for the atmosphere of hot Jupiters. We finally show the robustness of our energy equation, and its power to obtain physical insight into the modes. Conclusions. We have developed and benchmarked a code for the study of linear processes in planetary atmospheres, with an arbitrary steady state. The calculation of an a posteriori energy equation provides both increased robustness and physical meaning to the obtained eigenmodes. This code can be applied to various problems, and notably to further study the initial spin up of superrotation of GCM simulations of hot Jupiter.

Characterizing stellar convection in multiple dimensions is a topic at the forefront of stellar astrophysics. Numerical simulations are an essential tool for this task. We present an extension of the existing numerical tool-kit A-MaZe that enables such simulations of stratified flows in a gravitational field. The finite-volume based, cell-centered, and time-explicit hydrodynamics solver of A-MaZe was extended such that the scheme is now well-balanced in both momentum and energy. The algorithm maintains an initially static balance between gravity and pressure to machine precision. Quasi-stationary convection in slab-geometry preserves gas energy (internal plus kinetic) on average despite strong local up- and down-drafts. By contrast, a more standard numerical scheme is demonstrated to result in substantial gains of energy within a short time on purely numerical grounds. The test is further used to point out the role of dimensionality, viscosity, and Rayleigh number for compressible convection. Applications to a young sun in 2D and 3D, covering a part of the inner radiative zone as well as the outer convective zone, demonstrate that the scheme meets its initial design goal. Comparison with results obtained for a physically identical setup with a time-implicit code show qualitative agreement.

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.

In this white paper, we recommend the European Space Agency plays a proactive role in developing a global collaborative effort to construct a large high-contrast imaging space telescope, e.g. as currently under study by NASA. Such a mission will be needed to characterize a sizable sample of temperate Earth-like planets in the habitable zones of nearby Sun-like stars and to search for extraterrestrial biological activity. We provide an overview of relevant European expertise, and advocate ESA to start a technology development program towards detecting life outside the Solar system.

The detection of brown dwarfs in globular star clusters will allow us to break the degeneracies in age, mass and composition that affect our current models, and therefore to constrain the physics of their atmospheres and interiors. Furthermore, detecting brown dwarfs will help us constrain the properties of the clusters themselves, as they carry information about the clusters' age and dynamics. Great advancements in this field are to be expected in the next ten years, thanks to the extraordinary sensitivity in the infrared of upcoming telescopes like JWST and the ELTs.

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 significant differences in the simulated atmospheric flow for warm, tidally-locked small Neptunes and super Earths (based on a nominal GJ 1214b) when solving the simplified, and commonly used, primitive dynamical equations or the full Navier-Stokes equations. The dominant prograde, superrotating zonal jet is markedly different between the simulations which are performed using practically identical numerical setups, within the same model. The differences arise due to the breakdown of the so-called `shallow-fluid' and traditional approximations, which worsens when rotation rates are slowed, and day-night temperature contrasts are increased. The changes in the zonal advection between simulations solving the full and simplified equations, give rise to significant differences in the atmospheric redistribution of heat, altering the position of the hottest part of the atmosphere and temperature contrast between the day and night sides. The implications for the atmospheric chemistry and, therefore, observations need to be studied with a model including a more detailed treatment of the radiative transfer and chemistry. Small Neptunes and super Earths are extremely abundant and important, potentially bridging the structural properties (mass, radius, composition) of terrestrial and gas giant planets. Our results indicate care is required when interpreting the output of models solving the primitive equations of motion for such planets.

In this paper we present three-dimensional atmospheric simulations of the hot Jupiter HD~189733b under two different scenarios: local chemical equilibrium and including advection of the chemistry by the resolved wind. Our model consistently couples the treatment of dynamics, radiative transfer and chemistry, completing the feedback cycle between these three important processes. The effect of wind--driven advection on the chemical composition is qualitatively similar to our previous results for the warmer atmosphere of HD~209458b, found using the same model. However, we find more significant alterations to both the thermal and dynamical structure for the cooler atmosphere of HD~189733b, with changes in both the temperature and wind velocities reaching $\sim10\%$. We also present the contribution function, diagnosed from our simulations, and show that wind--driven chemistry has a significant impact on its three--dimensional structure, particularly for regions where methane is an important absorber. Finally, we present emission phase curves from our simulations and show the significant effect of wind--driven chemistry on the thermal emission, particularly within the 3.6 \textmu m Spitzer/IRAC channel.

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.

The precise measurement of the masses and radii of stars in eclipsing binary systems provides a window into uncertain processes in stellar evolution, especially mixing at convective boundaries. Recently, these data have been used to calibrate models of convective overshooting in the cores of main sequence stars. In this study we have used a small representative sample of eclipsing binary stars with $1.25 \leq M/\text{M}_\odot < 4.2$ to test how precisely this method can constrain the overshooting and whether the data support a universal stellar mass--overshooting relation. We do not recover the previously reported stellar mass dependence for the extent of overshooting and in each case we find there is a substantial amount of uncertainty, that is, the same binary pair can be matched by models with different amounts of overshooting. Models with a moderate overshooting parameter $0.013 \leq f_\text{os} \leq 0.014$ (using the scheme from Herwig et al. 1997) are consistent with all eight systems studied. Generally, a much larger range of $f_\text{os}$ is suitable for individual systems. In the case of main sequence and early post-main sequence stars, large changes in the amount of overshooting have little effect on the radius and effective temperature, and therefore the method is of extremely limited utility.

Atmospheric composition provides essential markers of the most fundamental properties of giant exoplanets, such as their formation mechanism or internal structure. New-generation exoplanet imagers, like VLT/SPHERE or Gemini/GPI, have been designed to achieve very high contrast (>15 mag) at small angular separations ($<$0.5\as) for the detection of young giant planets in the near-infrared, but they only provide very low spectral resolutions ($R<100$) for their characterization. High-dispersion spectroscopy at resolutions up to $10^5$ is one of the most promising pathways for the detailed characterization of exoplanets, but it is currently out of reach for most directly imaged exoplanets because current high-dispersion spectrographs in the near-infrared lack coronagraphs to attenuate the stellar signal and the spatial resolution necessary to resolve the planet. Project HiRISE (High-Resolution Imaging and Spectroscopy of Exoplanets) ambitions to develop a demonstrator that will combine the capabilities of two flagship instruments installed on the ESO Very Large Telescope, the high-contrast exoplanet imager SPHERE and the high-resolution spectrograph CRIRES+, with the goal of answering fundamental questions on the formation, composition and evolution of young planets. In this work, we will present the project, the first set of realistic simulations and the preliminary design of the fiber injection unit that will be implemented in SPHERE.

We present a study of the effect of wind-driven advection on the chemical composition of hot Jupiter atmospheres using a fully-consistent 3D hydrodynamics, chemistry and radiative transfer code, the Met Office Unified Model (UM). Chemical modelling of exoplanet atmospheres has primarily been restricted to 1D models that cannot account for 3D dynamical processes. In this work we couple a chemical relaxation scheme to the UM to account for the chemical interconversion of methane and carbon monoxide. This is done consistently with the radiative transfer meaning that departures from chemical equilibrium are included in the heating rates (and emission) and hence complete the feedback between the dynamics, thermal structure and chemical composition. In this letter we simulate the well studied atmosphere of HD~209458b. We find that the combined effect of horizontal and vertical advection leads to an increase in the methane abundance by several orders of magnitude; directly opposite to the trend found in previous works. Our results demonstrate the need to include 3D effects when considering the chemistry of hot Jupiter atmospheres. We calculate transmission and emission spectra, as well as the emission phase curve, from our simulations. We conclude that gas-phase non-equilibrium chemistry is unlikely to explain the model-observation discrepancy in the 4.5\,{\textmu m} {\it Spitzer}/IRAC channel. However, we highlight other spectral regions, observable with the James Webb Space Telescope, where signatures of wind-driven chemistry are more prominant.

In this work we have performed a series of simulations of the atmosphere of GJ~1214b assuming different metallicities using the Met Office Unified Model (UM). The UM is a general circulation model (GCM) that solves the deep, non-hydrostatic equations of motion and uses a flexible and accurate radiative transfer scheme, based on the two-stream and correlated-$k$ approximations, to calculate the heating rates. In this work we consistently couple a well-tested Gibbs energy minimisation scheme to solve for the chemical equilibrium abundances locally in each grid cell for a general set of elemental abundances, further improving the flexibility and accuracy of the model. As the metallicity of the atmosphere is increased we find significant changes in the dynamical and thermal structure, with subsequent implications for the simulated phase curve. The trends that we find are qualitatively consistent with previous works, though with quantitative differences. We investigate in detail the effect of increasing the metallicity by splitting the mechanism into constituents, involving the mean molecular weight, the heat capacity and the opacities. We find the opacity effect to be the dominant mechanism in altering the circulation and thermal structure. This result highlights the importance of accurately computing the opacities and radiative transfer in 3D GCMs.

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.

We study lithium depletion in low-mass and solar-like stars as a function of time, using a new diffusion coefficient describing extra-mixing taking place at the bottom of a convective envelope. This new form is motivated by multi-dimensional fully compressible, time implicit hydrodynamic simulations performed with the MUSIC code. Intermittent convective mixing at the convective boundary in a star can be modeled using extreme value theory, a statistical analysis frequently used for finance, meteorology, and environmental science. In this letter, we implement this statistical diffusion coefficient in a one-dimensional stellar evolution code, using parameters calibrated from multi-dimensional hydrodynamic simulations of a young low-mass star. We propose a new scenario that can explain observations of the surface abundance of lithium in the Sun and in clusters covering a wide range of ages, from $\sim$ 50 Myr to $\sim$ 4 Gyr. Because it relies on our physical model of convective penetration, this scenario has a limited number of assumptions. It can explain the observed trend between rotation and depletion, based on a single additional assumption, namely that rotation affects the mixing efficiency at the convective boundary. We suggest the existence of a threshold in stellar rotation rate above which rotation strongly prevents the vertical penetration of plumes and below which rotation has small effects. In addition to providing a possible explanation for the long standing problem of lithium depletion in pre-main sequence and main sequence stars, the strength of our scenario is that its basic assumptions can be tested by future hydrodynamic simulations.

We examine a penetration layer formed between a central radiative zone and a large convection zone in the deep interior of a young low-mass star. Using the Multidimensional Stellar Implicit Code (MUSIC) to simulate two-dimensional compressible stellar convection in a spherical geometry over long times, we produce statistics that characterize the extent and impact of convective penetration in this layer. We apply extreme value theory to the maximal extent of convective penetration at any time. We compare statistical results from simulations which treat non-local convection, throughout a large portion of the stellar radius, with simulations designed to treat local convection in a small region surrounding the penetration layer. For each of these situations, we compare simulations of different resolution, which have different velocity magnitudes. We also compare statistical results between simulations that radiate energy at a constant rate to those that allow energy to radiate from the stellar surface according to the local surface temperature. Based on the frequency and depth of penetrating convective structures, we observe two distinct layers that form between the convection zone and the stable radiative zone. We show that the probability density function of the maximal depth of convective penetration at any time corresponds closely in space with the radial position where internal waves are excited. We find that the maximal penetration depth can be modeled by a Weibull distribution with a small shape parameter. Using these results, and building on established scalings for diffusion enhanced by large-scale convective motions, we propose a new form for the diffusion coefficient that may be used for one-dimensional stellar evolution calculations in the large P\'eclet number regime. These results should contribute to the 321D link.

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 present highlights from a large set of simulations of a hot Jupiter atmosphere, nominally based on HD 209458b, aimed at exploring both the evolution of the deep atmosphere, and the acceleration of the zonal flow or jet. We find the occurrence of a super-rotating equatorial jet is robust to changes in various parameters, and over long timescales, even in the absence of strong inner or bottom boundary drag. This jet is diminished in one simulation only, where we strongly force the deep atmosphere equator-to-pole temperature gradient over long timescales. Finally, although the eddy momentum fluxes in our atmosphere show similarities with the proposed mechanism for accelerating jets on tidally-locked planets, the picture appears more complex. We present tentative evidence for a jet driven by a combination of eddy momentum transport and mean flow.

We present the results of a numerical benchmark study for the MUlti-dimensional Stellar Implicit Code (MUSIC) based on widely applicable two- and three-dimensional compressible hydrodynamics problems relevant to stellar interiors. MUSIC is an implicit large eddy simulation code that uses implicit time integration, implemented as a Jacobian-free Newton Krylov method. A physics based preconditioning technique which can be adjusted to target varying physics is used to improve the performance of the solver. The problems used for this benchmark study include the Rayleigh-Taylor and Kelvin-Helmholtz instabilities, and the decay of the Taylor-Green vortex. Additionally we show a test of hydrostatic equilibrium, in a stellar environment which is dominated by radiative effects. In this setting the flexibility of the preconditioning technique is demonstrated. This work aims to bridge the gap between the hydrodynamic test problems typically used during development of numerical methods and the complex flows of stellar interiors. A series of multi-dimensional tests are performed and analysed. Each of these test cases is analysed with a simple, scalar diagnostic, with the aim of enabling direct code comparisons. As the tests performed do not have analytic solutions we verify MUSIC by comparing to established codes including ATHENA and the PENCIL code. MUSIC is able to both reproduce behaviour from established and widely-used codes as well as results expected from theoretical predictions. This benchmarking study concludes a series of papers describing the development of the MUSIC code and provides confidence in the future applications.

The correlated-k method is frequently used to speed up radiation calculations in both one-dimensional and three-dimensional atmosphere models. An inherent difficulty with this method is how to treat overlapping absorption, i.e. absorption by more than one gas in a given spectral region. We have evaluated the applicability of three different methods in hot Jupiter and brown dwarf atmosphere models, all of which have been previously applied within models in the literature: (i) Random overlap, both with and without resorting and rebinning, (ii) equivalent extinction and (iii) pre-mixing of opacities, where (i) and (ii) combine k-coefficients for different gases to obtain k-coefficients for a mixture of gases, while (iii) calculates k-coefficients for a given mixture from the corresponding mixed line-by-line opacities. We find that the random overlap method is the most accurate and flexible of these treatments, and is fast enough to be used in one-dimensional models with resorting and rebinning. In three-dimensional models such as GCMs it is too slow, however, and equivalent extinction can provide a speed-up of at least a factor of three with only a minor loss of accuracy while at the same time retaining the flexibility gained by combining k-coefficients computed for each gas individually. Pre-mixed opacities are significantly less flexible, and we also find that particular care must be taken when using this method in order to to adequately resolve steep variations in composition at important chemical equilibrium boundaries. We use the random overlap method with resorting and rebinning in our one-dimensional atmosphere model and equivalent extinction in our GCM, which allows us to e.g. consistently treat the feedback of non-equilibrium chemistry on the total opacity and therefore the calculated P-T profiles in our models.

To study the complexity of hot Jupiter atmospheres revealed by observations of increasing quality, we have adapted the UK Met Office global circulation model (GCM), the Unified Model (UM), to these exoplanets. The UM solves the full 3D Navier-Stokes equations with a height-varying gravity, avoiding the simplifications used in most GCMs currently applied to exoplanets. In this work we present the coupling of the UM dynamical core to an accurate radiation scheme based on the two-stream approximation and correlated-k method with state-of-the-art opacities from ExoMol. Our first application of this model is devoted to the extensively studied hot Jupiter HD 209458b. We have derived synthetic emission spectra and phase curves, and compare them to both previous models also based on state-of-the-art radiative transfer, and to observations. We find a reasonable a agreement between observations and both our day side emission and hot spot offset, however, our night side emission is too large. Overall our results are qualitatively similar to those found by Showman et al. (2009) with the SPARC/MITgcm, however, we note several quantitative differences: Our simulations show significant variation in the position of the hottest part of the atmosphere with pressure, as expected from simple timescale arguments, and in contrast to the "vertical coherency" found by Showman et al. (2009). We also see significant quantitative differences in calculated synthetic observations. Our comparisons strengthen the need for detailed intercomparisons of dynamical cores, radiation schemes and post-processing tools to understand these differences. This effort is necessary in order to make robust conclusions about these atmospheres based on GCM results.

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.

We present a description of a new instrument development, HARPS3, planned to be installed on an upgraded and roboticized Isaac Newton Telescope by end-2018. HARPS3 will be a high resolution (R = 115,000) echelle spectrograph with a wavelength range from 380-690 nm. It is being built as part of the Terra Hunting Experiment - a future 10 year radial velocity measurement programme to discover Earth-like exoplanets. The instrument design is based on the successful HARPS spectrograph on the 3.6m ESO telescope and HARPS-N on the TNG telescope. The main changes to the design in HARPS3 will be: a customised fibre adapter at the Cassegrain focus providing a stabilised beam feed and on-sky fibre diameter ~ 1.4 arcsec, the implementation of a new continuous flow cryostat to keep the CCD temperature very stable, detailed characterisation of the HARPS3 CCD to map the effective pixel positions and thus provide an improved accuracy wavelength solution, an optimised integrated polarimeter and the instrument integrated into a robotic operation. The robotic operation will optimise our programme which requires our target stars to be measured on a nightly basis. We present an overview of the entire project, including a description of our anticipated robotic operation.

In this work we investigate the impact of calculating non-equilibrium chemical abundances consistently with the temperature structure for the atmospheres of highly-irradiated, close-in gas giant exoplanets. Chemical kinetics models have been widely used in the literature to investigate the chemical compositions of hot Jupiter atmospheres which are expected to be driven away from chemical equilibrium via processes such as vertical mixing and photochemistry. All of these models have so far used pressure--temperature (P-T) profiles as fixed model input. This results in a decoupling of the chemistry from the radiative and thermal properties of the atmosphere, despite the fact that in nature they are intricately linked. We use a one-dimensional radiative-convective equilibrium model, ATMO, which includes a sophisticated chemistry scheme to calculate P-T profiles which are fully consistent with non-equilibrium chemical abundances, including vertical mixing and photochemistry. Our primary conclusion is that, in cases of strong chemical disequilibrium, consistent calculations can lead to differences in the P-T profile of up to 100 K compared to the P-T profile derived assuming chemical equilibrium. This temperature change can, in turn, have important consequences for the chemical abundances themselves as well as for the simulated emission spectra. In particular, we find that performing the chemical kinetics calculation consistently can reduce the overall impact of non-equilibrium chemistry on the observable emission spectrum of hot Jupiters. Simulated observations derived from non-consistent models could thus yield the wrong interpretation. We show that this behaviour is due to the non-consistent models violating the energy budget balance of the atmosphere.

Context: We study the impact of two-dimensional spherical shells on compressible convection. Realistic profiles for density and temperature from a one-dimensional stellar evolution code are used to produce a model of a large stellar convection zone representative of a young low-mass star. Methods: We perform hydrodynamic implicit large-eddy simulations of compressible convection using the MUltidimensional Stellar Implicit Code (MUSIC). Because MUSIC has been designed to use realistic stellar models produced from one-dimensional stellar evolution calculations, MUSIC simulations are capable of seamlessly modeling a whole star. Simulations in two-dimensional spherical shells that have different radial extents are performed over hundreds of convective turnover times, permitting the collection of well-converged statistics. Results: We evaluate basic statistics of the convective turnover time, the convective velocity, and the overshooting layer. These quantities are selected for their relevance to one-dimensional stellar evolution calculations, so that our results are focused toward the 321D link. The inclusion in the spherical shell of the boundary between the radiative and convection zones decreases the amplitude of convective velocities in the convection zone. The inclusion of near-surface layers in the spherical shell can increase the amplitude of convective velocities, although the radial structure of the velocity profile established by deep convection is unchanged. The impact from including the near-surface layers depends on the speed and structure of small-scale convection in the near-surface layers. Larger convective velocities in the convection zone result in a commensurate increase in the overshooting layer width and decrease in the convective turnover time. These results provide support for non-local aspects of convection.

HD 142527A is one of the most studied Herbig Ae/Be stars with a transitional disk, as it has the largest imaged gap in any protoplanetary disk: the gas is cleared from 30 to 90 AU. The HD142527 system is also unique in that it has a stellar companion with a small mass compared to the mass of the primary star. This factor of ~20 in mass ratio between the two objects makes this binary system different from any other YSO. The HD142527 system could therefore provide a valuable test bed. This low-mass stellar object may be responsible for both the gap and dust trapping observed by ALMA at longer distances. We observed this system with the NACO and GPI instruments using the aperture masking technique. Aperture masking is ideal for providing high dynamic range even at very small angular separations. We present the spectral energy distribution for HD142527A and B. Brightness of the companion is now known from the R band up to the M' band. We also followed the orbital motion of HD 142527B over a period of more than two years. The SED of the companion is compatible with a T=3000+/-100K object in addition to a 1700K blackbody environment (likely a circus-secondary disk). From evolution models, we find that it is compatible with an object of mass 0.13+/-0.03Msun, radius 0.90+/-0.15Rsun, and age $1.0^{+1.0}_{-0.75}$Myr. This age is significantly younger than the age previously estimated for HD142527A. Computations to constrain the orbital parameters found a semi major axis of $140^{+120}_{-70}$mas, an eccentricity of 0.5+/-0.2, an inclination of 125+/-15 degrees, and a position angle of the right ascending node of -5+/-40 degrees. Inclination and position angle of the ascending node are in agreement with an orbit coplanar with the inner disk, not coplanar with the outer disk. Despite its high eccentricity, it is unlikely that HD142527B is responsible for truncating the inner edge of the outer disk.

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.

This work is the first attempt to describe the multi-dimensional structure of accreting young stars based on fully compressible time implicit multi-dimensional hydrodynamics simulations. One major motivation is to analyse the validity of accretion treatment used in previous 1D stellar evolution studies. We analyse the effect of accretion on the structure of a realistic stellar model of the young Sun. Our work is inspired by the numerical work of Kley \& Lin (1996, ApJ, 461, 933) devoted to the structure of the boundary layer in accretion disks. We analyse the redistribution of accreted material with a range of values of specific entropy relative to the bulk specific entropy of the material in the accreting object's convective envelope. A primary goal is to understand whether and how accreted energy deposited onto a stellar surface is redistributed in the interior. This study focusses on the high accretion rates characteristic of FU Ori systems. We find that the highest entropy cases produce a distinctive behaviour in the mass redistribution, rms velocities, and enthalpy flux in the convective envelope. This change in behaviour is characterised by the formation of a hot layer on the surface of the accreting object, which tends to suppress convection in the envelope. We analyse the long-term effect of such a hot buffer zone on the structure and evolution of the accreting object with 1D stellar evolution calculations. We study the relevance of the assumption of redistribution of accreted energy into the stellar interior used in the literature. One conclusion is that, for a given amount of accreted energy transferred to the accreting object, a treatment assuming accretion energy redistribution throughout the stellar interior could significantly overestimate the effects on the stellar structure, in particular, on the resulting expansion.

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.

This work is a continuation of our efforts to develop an efficient implicit solver for multidimensional hydrodynamics for the purpose of studying important physical processes in stellar interiors, such as turbulent convection and overshooting. We present an implicit solver that results from the combination of a Jacobian-Free Newton-Krylov method and a preconditioning technique tailored to the inviscid, compressible equations of stellar hydrodynamics. We assess the accuracy and performance of the solver for both 2D and 3D problems for Mach numbers down to $10^{-6}$. Although our applications concern flows in stellar interiors, the method can be applied to general advection and/or diffusion-dominated flows. The method presented in this paper opens up new avenues in 3D modeling of realistic stellar interiors allowing the study of important problems in stellar structure and evolution.

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.

We present high-speed, three-colour photometry of the eclipsing dwarf nova PHL 1445, which, with an orbital period of 76.3 min, lies just below the period minimum of ~82 min for cataclysmic variable stars. Averaging four eclipses reveals resolved eclipses of the white dwarf and bright spot. We determined the system parameters by fitting a parameterised eclipse model to the averaged lightcurve. We obtain a mass ratio of q = 0.087 +- 0.006 and inclination i = 85.2 +- 0.9 degrees. The primary and donor masses were found to be Mw = 0.73 +- 0.03 Msun and Md = 0.064 +- 0.005 Msun, respectively. Through multicolour photometry a temperature of the white dwarf of Tw = 13200 +- 700 K and a distance of 220 +- 50 pc were determined. The evolutionary state of PHL 1445 is uncertain. We are able to rule out a significantly evolved donor, but not one that is slightly evolved. Formation with a brown dwarf donor is plausible; though the brown dwarf would need to be no older than 600 Myrs at the start of mass transfer, requiring an extremely low mass ratio (q = 0.025) progenitor system. PHL 1445 joins SDSS 1433 as a sub-period minimum CV with a substellar donor. These existence of two such systems raises an alternative possibility; that current estimates for the intrinsic scatter and/or position of the period minimum may be in error.

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.

The treatment of radiation transport in global circulation models (GCMs) is crucial to correctly describe Earth and exoplanet atmospheric dynamics processes. The two-stream approximation and correlated-k method are currently state-of-the-art approximations applied in both Earth and hot Jupiter GCM radiation schemes to facilitate rapid calculation of fluxes and heating rates. Their accuracy have been tested extensively for Earth-like conditions, but verification of the methods' applicability to hot Jupiter-like conditions is lacking in the literature. We are adapting the UK Met Office GCM, the Unified Model (UM), for the study of hot Jupiters, and present in this work the adaptation of the Edwards-Slingo radiation scheme based on the two-stream approximation and the correlated-k method. We discuss the calculation of absorption coefficients from high temperature line lists and highlight the large uncertainty in the pressure-broadened line widths. We compare fluxes and heating rates obtained with our adapted scheme to more accurate discrete ordinate (DO) line-by-line (LbL) calculations ignoring scattering effects. We find that, in most cases, errors stay below 10 % for both heating rates and fluxes using ~ 10 k-coefficients in each band and a diffusivity factor D = 1.66. The two-stream approximation and the correlated-k method both contribute non-negligibly to the total error. We also find that using band-averaged absorption coefficients, which have previously been used in radiative-hydrodynamical simulations of a hot Jupiter, may yield errors of ~ 100 %, and should thus be used with caution.

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.

The physics of brown dwarfs has continuously improved since the discovery of these astrophysical bodies. The first important developments were devoted to the description of their mechanical structure, with the derivation of an appropriate equation of state, and the modelling of their atmosphere characterised by strong molecular absorption. New challenges are arising with progress in observational techniques which provide data of unprecedented accuracy. The goal of this chapter is to describe some of the current challenges for the theory of brown dwarfs. Those challenges concerns atmospheric dust and cloud, non-equilibrium atmospheric chemistry, the effect of rotation and magnetic fields on internal structure and the very early phases of evolution characterised by accretion processes. The field remains lively as more and more high quality observational data become available and because of increasing discoveries of exoplanets. Indeed, many physical properties of giant exoplanets can be described by the same theory as brown dwarfs, as described in this chapter.

We are adapting the Global Circulation Model (GCM) of the UK Met Office, the so--called Unified Model (UM), for the study of hot Jupiters. In this work we demonstrate the successful adaptation of the most sophisticated dynamical core, the component of the GCM which solves the equations of motion for the atmosphere, available within the UM, ENDGame (Even Newer Dynamics for General atmospheric modelling of the environment). Within the same numerical scheme ENDGame supports solution to the dynamical equations under varying degrees of simplification. We present results from a simple, shallow (in atmospheric domain) hot Jupiter model (SHJ), and a more realistic (with a deeper atmosphere) HD 209458b test case. For both test cases we find that the large--scale, time--averaged (over the 1200 days prescribed test period), dynamical state of the atmosphere is relatively insensitive to the level of simplification of the dynamical equations. However, problems exist when attempting to reproduce the results for these test cases derived from other models. For the SHJ case the lower (and upper) boundary intersects the dominant dynamical features of the atmosphere meaning the results are heavily dependent on the boundary conditions. For the HD 209458b test case, when using the more complete dynamical models, the atmosphere is still clearly evolving after 1200 days, and in a transient state. Solving the complete (deep atmosphere and non--hydrostatic) dynamical equations allows exchange between the vertical and horizontal momentum of the atmosphere, via Coriolis and metric terms. Subsequently, interaction between the upper atmosphere and the deeper more slowly evolving (radiatively inactive) atmosphere significantly alters the results, and acts over timescales longer than 1200 days.

We demonstrate that both the current (New Dynamics), and next generation (ENDGame) dynamical cores of the UK Met Office global circulation model, the UM, reproduce consistently, the long-term, large-scale flows found in several published idealised tests. The cases presented are the Held-Suarez test, a simplified model of Earth (including a stratosphere), and a hypothetical tidally locked Earth. Furthermore, we show that using simplifications to the dynamical equations, which are expected to be justified for the physical domains and flow regimes we have studied, and which are supported by the ENDGame dynamical core, also produces matching long-term, large-scale flows. Finally, we present evidence for differences in the detail of the planetary flows and circulations resulting from improvements in the ENDGame formulation over New Dynamics.

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.

Time-implicit schemes are attractive since they allow numerical time steps that are much larger than those permitted by the Courant-Friedrich-Lewy criterion characterizing time-explicit methods. This advantage comes, however, with a cost: the solution of a system of nonlinear equations is required at each time step. In this work, the nonlinear system results from the discretization of the hydrodynamical equations with the Crank-Nicholson scheme. We compare the cost of different methods, based on Newton-Raphson iterations, to solve this nonlinear system, and benchmark their performances against time-explicit schemes. Since our general scientific objective is to model stellar interiors, we use as test cases two realistic models for the convective envelope of a red giant and a young Sun. Focusing on 2D simulations, we show that the best performances are obtained with the quasi-Newton method proposed by Broyden. Another important concern is the accuracy of implicit calculations. Based on the study of an idealized problem, namely the advection of a single vortex by a uniform flow, we show that there are two aspects: i) the nonlinear solver has to be accurate enough to resolve the truncation error of the numerical discretization, and ii) the time step has be small enough to resolve the advection of eddies. We show that with these two conditions fulfilled, our implicit methods exhibit similar accuracy to time-explicit schemes, which have lower values for the time step and higher computational costs. Finally, we discuss in the conclusion the applicability of these methods to fully implicit 3D calculations.

Studies of beryllium abundance in large samples of solar-type stars show a small fraction of extremely beryllium-deficient stars, which challenges our current understanding of light element depletion in these stars. We suggest two possible scenarios that may explain this high level of Be depletion: early accretion and rotational mixing. We show that in both cases, the conditions required to reach the observed level of Be depletion are quite extreme, which explains the very small fraction of detected Be outliers. We suggest that substantial Be depletion can be obtained in stars if they were fast rotators in the past, with high initial rotational velocities and short disc lifetimes. Our analysis suggests that rotational mixing may not be efficient enough to deplete Be in less than 10 Myr. Consequently, the detection of strongly Be-deficient stars in clusters younger than $\sim$ 10 Myr may provide a genuine signature of accretion process and the proof that some protostars may undergo many extreme bursts of accretion during their embedded phases of evolution.

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]

This paper describes the first steps of development of a new multidimensional time implicit code devoted to the study of hydrodynamical processes in stellar interiors. The code solves the hydrodynamical equations in spherical geometry and is based on the finite volume method. Radiation transport is taken into account within the diffusion approximation. Realistic equation of state and opacities are implemented, allowing the study of a wide range of problems characteristic of stellar interiors. We describe in details the numerical method and various standard tests performed to validate the method. We present preliminary results devoted to the description of stellar convection. We first perform a local simulation of convection in the surface layers of a A-type star model. This simulation is used to test the ability of the code to address stellar conditions and to validate our results, since they can be compared to similar previous simulations based on explicit codes. We then present a global simulation of turbulent convective motions in a cold giant envelope, covering 80% in radius of the stellar structure. Although our implicit scheme is unconditionally stable, we show that in practice there is a limitation on the time step which prevent the flow to move over several cells during a time step. Nevertheless, in the cold giant model we reach a hydro CFL number of 100. We also show that we are able to address flows with a wide range of Mach numbers (10^-3 < Ms< 0.5), which is impossible with an anelastic approach. Our first developments are meant to demonstrate that the use of an implicit scheme applied to a stellar evolution context is perfectly thinkable and to provide useful guidelines to optimise the development of an implicit multi-D hydrodynamical code.

We reconstruct the evolutionary path followed by cataclysmic variables (CVs) from the observed mass-radius relationship of their donor stars. Along the way, we update the semi-empirical CV donor sequence of Knigge (2006) and present a comprehensive review of the link between CV evolution and donor physics. After calibrating state-of-the art stellar models for use in the CV setting, we fit self-consistent theoretical evolution sequences to the observed donor masses and radii. In the standard model of CV evolution, AML below the period gap is assumed to be driven solely by gravitational radiation (GR), while AML above the gap is usually described by a magnetic braking prescription due to Rappaport, Verbunt & Joss (1983). We find that simple scaled versions of these recipes match the data quite well. However, the optimal scaling factors turn out to be f_GR = 2.47 +/- 0.22 below the gap and f_MB = 0.66 +/- 0.05 above. The implications and applications of our results include: (1) The revised evolution sequence yields correct locations for the CV minimum period and the upper edge of the period gap; the standard sequence does not. (2) A comparison of predicted and observed WD temperatures suggests an even higher value for f_GR, but this is sensitive to the assumed WD mass. (3) The absolute donor magnitudes predicted by our sequences can be used to set firm lower limits on the distances toward CVs. (4) Both standard and revised sequences predict that short-period CVs should be susceptible to dwarf nova (DN) eruptions, consistent with observations. However, both sequences also predict that the DNe fraction among long-period CVs should decline with P_orb. Observations suggest the opposite behaviour. (5) The ratio of long-period CVs to short-period, pre-bounce CV is about 3x higher for the revised sequence than the standard one. This may resolve a long-standing problem in CV evolution. [abridged]

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.

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.

We report on the development of an implicit multi-D hydrodynamic code for stellar evolution. We present two test-cases relevant for the first scientific goal of the code: the simulation of convection in pulsating stars. First results on a realistic stellar model are also presented.

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.

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.

We discuss the interior structure and composition of giant planets, and how this structure changes as these planets cool and contract over time. Here we define giant planets as those that have an observable hydrogen-helium envelope, which includes Jupiter-like planets, which are predominantly H/He gas, and Neptune-like planets which are predominantly composed of elements heavier than H/He. We describe the equations of state of planetary materials and the construction of static structural models and thermal evolution models. We apply these models to transiting planets close to their parent stars, as well as directly imaged planets far from their parent stars. Mechanisms that have been postulated to inflate the radii of close-in transiting planets are discussed. We also review knowledge gained from the study of the solar system's giant planets. The frontiers of giant planet physics are discussed with an eye towards future planetary discoveries.

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.

In the context of the space-based mission CoRoT, devoted to asteroseismology and search for planet transits, we analyse the accuracy of fundamental stellar parameters (mass, radius, luminosity) that can be obtained from asteroseismological data.} Our work is motivated by the large uncertainties on planetary radius determination of transiting planets which are mainly due to uncertainties on the stellar parameters. Our goal is to analyse uncertainties of fundamental stellar parameters for a given accuracy of oscillation frequency determination. We generate grids of equilibrium models of stars and compute their pulsation spectra based on a linear nonadiabatic stability analysis. Using differents methods of comparison of oscillation mode spectra, we derive uncertainties on fundamental stellar parameters and analyse the effect of varying the number of considered modes.} The limits obtained depend strongly on the adapted method to compare spectra. We find a degeneracy in the stellar parameter solutions, up to a few % in mass (from less than 1% to more than 7% depending on the method used and the number of considered modes), luminosity (from 2% to more than 10%) or radius (from less than 1% to 3%), for a given pulsation spectrum.

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 use deep Hubble Space Telescope photometry of the rich, young (~20-45 Myr-old) star cluster NGC 1818 in the Large Magellanic Cloud to derive its stellar mass function (MF) down to ~0.15 Msun. This represents the deepest robust MF thus far obtained for a stellar system in an extragalactic, low-metallicity ([Fe/H]~-0.4 dex) environment. Combining our results with the published MF for masses above 1.0 Msun, we obtain a complete present-day MF. This is a good representation of the cluster's initial MF (IMF), particularly at low masses, because our observations are centred on the cluster's uncrowded half-mass radius. Therefore, stellar and dynamical evolution of the cluster will not have affected the low-mass stars significantly. The NGC 1818 IMF is well described by both a lognormal and a broken power-law distribution with slopes of Gamma=0.46+/-0.10 and Gamma~-1.35 (Salpeter-like) for masses in the range from 0.15 to 0.8 Msun and greater than 0.8 Msun, respectively. Within the uncertainties, the NGC 1818 IMF is fully consistent with both the Kroupa solar-neighbourhood and the Chabrier lognormal mass distributions.

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.

Despite the revolution in our knowledge resulting from the detection of planets around mature stars, we know almost nothing about planets orbiting young stars because rapid rotation and active photospheres preclude detection by radial velocities or transits and because direct imaging has barely penetrated the requisite range of high contrast and angular resolution. Of the techniques presently under consideration for the coming decade, only space-based astrometry offers the prospect of discovering gas giants (100 to >> 300 Mearth), lower mass systems such as icy giants (10 to 100 Mearth), and even a few rocky, super-Earths 300 Mearth) orbiting stars ranging in age from 1 to 100 Myr. Astrometry will complement high contrast imaging which should be able to detect gas giants (1~10 MJup) in orbits from a few to a few hundred AU. An astrometric survey in combination with imaging data for a subsample of objects will allow a detailed physical understanding of the formation and evolution of young gas giant planets impossible to achieve by any one technique.

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

(Abridged) We present high-speed, three-colour photometry of seven short period (Porb <= 95 mins) eclipsing CVs from the Sloan Digital Sky Survey. We determine the system parameters via a parametrized model of the eclipse fitted to the observed lightcurve by chi2 minimization. Three out of seven of the systems possess brown dwarf donor stars and are believed to have evolved past the orbital period minimum. This is in line with predictions that 40-70 per cent of CVs should have evolved past the orbital period minimum. Therefore, the main result of our study is that the missing population of post-period minimum CVs has finally been identified. The donor star masses and radii are, however, inconsistent with model predictions; the donor stars are approximately 10 per cent larger than expected across the mass range studied here. One explanation for the discrepancy is enhanced angular momentum loss (e.g. from circumbinary discs), however the mass-transfer rates, as deduced from white dwarf effective temperatures, are not consistent with enhanced angular momentum loss. We show it is possible to explain the large donor radii without invoking enhanced angular momentum loss by a combination of geometrical deformation and the effects of starspots due to strong rotation and expected magnetic activity. Choosing unambiguously between these different solutions will require independent estimates of the mass-transfer rates in short period CVs.

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

We conducted a search for brown dwarfs (BDs) and very low mass (VLM) stars in the 625 Myr-old Hyades cluster in order to derive the cluster's mass function across the stellar-substellar boundary. We performed a deep (I=23, z=22.5) photometric survey over 16 sq.deg. around the cluster center, followed up with K-band photometry to measure the proper motion of candidate members, and optical and near-IR spectroscopy of probable BD and VLM members. We report the discovery of the first 2 brown dwarfs in the Hyades cluster. The 2 objects have a spectral type early-T and their optical and near-IR photometry as well as their proper motion are consistent with them being cluster members. According to models, their mass is 50 Jupiter masses at an age of 625 Myr. We also report the discovery of 3 new very low mass stellar members of the cluster, and confirm the membership of 16 others. We combine these results with a list of previously known cluster members to build the present-day mass function (PDMF) of the Hyades cluster from 50 Jupiter masses to 3Mo. We find the Hyades PDMF to be strongly deficient in very low mass objects and brown dwarfs compared to the IMF of younger open clusters such as the Pleiades. We interpret this deficiency as the result of dynamical evolution over the past few 100 Myr, i.e., the preferential evaporation of low mass cluster members due to weak gravitational encounters. We thus estimate that the Hyades cluster currently hosts about 10-15 brown dwarfs, while its initial substellar population may have amounted up to 150-200 members.

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 present high-speed, three-colour photometry of the eclipsing cataclysmic variable SDSS J150722.30+523039.8 (hereafter SDSS J1507). This system has an orbital period of 66.61 minutes, placing it below the observed ``period minimum'' for cataclysmic variables. We determine the system parameters via a parameterised model of the eclipse fitted to the observed lightcurve by chi^2 minimisation. We obtain a mass ratio of q = 0.0623 +/- 0.0007 and an orbital inclination i = 83.63 +/- 0.05 degrees. The primary mass is M_w = 0.90 +/- 0.01 M_sun. The secondary mass and radius are found to be M_r =0.056 +/- 0.001 M_sun and R_r = 0.096 +/- 0.001 R_sun respectively. We find a distance to the system of 160 +/- 10 pc. The secondary star in SDSS J1507 has a mass substantially below the hydrogen burning limit, making it the second confirmed sub-stellar donor in a cataclysmic variable. The very short orbital period of SDSS J1507 is readily explained if the secondary star is nuclearly evolved, or if SDSS J1507 formed directly from a detached white dwarf/brown dwarf binary. Given the lack of any visible contribution from the secondary star, the very low secondary mass and the low HeI(6678AA)/Halpha emission line ratio, we argue that SDSS J1507 probably formed directly from a detached white dwarf/brown dwarf binary. If confirmed, SDSS J1507 will be the first such system identified. The implications for binary star evolution, the brown-dwarf desert and the common envelope phase are discussed.

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.

We analyse the presence of nonradial oscillations in Cepheids, a problem which has not been theoretically revised since the work of Dziembowsky (1977) and Osaki (1977). Our analysis is motivated by a work of Moskalik et al. (2004) which reports the detec tion of low amplitude periodicities in a few Cepheids of the large Magellanic cloud. These newly discovered periodicities were interpreted as nonradial modes.} {Based on linear nonadiabatic stability analysis, our goal is to reanalyse the presence and stability of nonradial modes, taking into account improvement in the main input phys ics required for the modelling of Cepheids.} {We compare the results obtained from two different numerical methods used to solve the set of differential equations: a matrix method and the Ricatti method.} {We show the limitation of the matrix method to find low order p-modes ($l<6$), because of their dual character in evolved stars such as Cepheids. For higher order p-modes, we find an excellent agreement between the two methods.} {No nonradial instability is found below $l=5$, whereas many unstable nonradial modes exist for higher orders. We also find that nonradial modes remain unstable, even at hotter effective temperatures than the blue edge of the Cepheid instability strip, where no radial pulsations are expected.

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

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.

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 present the results of a high-resolution imaging survey for brown dwarf binaries in the Pleiades open cluster. The observations were carried out with the Advance Camera for Surveys onboard the Hubble Space Telescope. Our sample consists of 15 bona-fide brown dwarfs. We confirm 2 binaries and detect their orbital motion, but we did not resolve any new binary candidates in the separation range between 5.4AU and 1700AU and masses in the range 0.035--0.065~Msun. Together with the results of our previous study (Martin et al., 2003), we can derive a visual binary frequency of 13.3$^{+13.7}_{-4.3}$\% for separations greater than 7~AU masses between 0.055--0.065~M$_{\sun}$ and mass ratios between 0.45--0.9$<q<$1.0. The other observed properties of Pleiades brown dwarf binaries (distributions of separation and mass ratio) appear to be similar to their older counterparts in the field.

We have used 3.5 to 8 micron data from the Cores to Disks (c2d) Legacy survey and our own deep IJHKs images of a 0.5 square degree portion of the c2d fields in Ophiuchus to produce a sample of candidate young objects with probable masses between 1 and 10 Jupiter masses. The availability of photometry over whole range where these objects emit allows us to discriminate between young, extremely low-mass candidates and more massive foreground and background objects and means our survey will have fewer false positives than existing near-IR surveys. The sensitive inventory of a star forming cloud from the red to the mid-IR will allow us to constrain the IMF for these non-clustered star formation regions to well below the deuterium burning limit. For stars with fluxes in the broad gap between the 2MASS limits and our limits, our data will provide information about the photospheres. We will use the Spitzer results in combination with current disk models to learn about the presence and nature of circumstellar disks around young brown dwarfs.

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.

Very low-mass stars and brown dwarfs can undergo pulsational instability excited by central deuterium burning during the initial phases of their evolution. We present the results of evolutionary and nonadiabatic linear stability models that show the presence of unstable fundamental modes. The pulsation periods vary bewteen ~5 hr for a 0.1 Msun star and ~1 hr for a 0.02 M$_\odot$ brown dwarf. The results are rather insensitive to variations in the input physics of the models. We show the location of the instability strip in the HR and c-m diagrams and discuss the observational searches for young pulsators in nearby star forming regions.

We present the results of a nonadiabatic, linear stability analysis of models of very low-mass stars (VLMSs) and brown dwarfs (BDs) during the deuterium burning phase in the center. We find unstable fundamental modes with periods varying between ~5 hr for a 0.1 Msun star and ~1 hr for a 0.02 Msun BD. The growth time of the instability decreases with decreasing mass and remains well below the deuterium burning time scale in the mass range considered (0.1--0.02 Msun). These results are robust against variations of the relevant input physics in the evolutionary models. We identify possible candidates for pulsational variability among known VLMSs and BDs in nearby star forming regions whose location in the HR diagram falls within or close to the boundary of the instability strip. Finally, we discuss the possibility that the variability observed in a few objects with periods of ~1 hr can be interpreted in terms of pulsation.

We analyze the possible existence of a pulsational instability excited by the $\epsilon$-mechanism during the last few centuries of evolution of a Chandrasekhar mass white dwarf prior to its explosion as a Type Ia supernova. Our analysis is motivated by the temperature sensitivity of the nuclear energy generation rate ($\sim T^{23}$) in a white dwarf whose structural adiabatic index is near 4/3. Based upon a linear stability analysis, we find that the fundamental mode and higher order radial modes are indeed unstable and that the fundamental mode has the shortest growth time scale. However, the growth time scale for such instability never becomes shorter than the evolutionary timescale. Therefore, even though the star \emph{is} pulsationally unstable, we do not expect these radial modes to have time to grow and to affect the structure and explosion properties of Type Ia supernovae.

We present here the results of astrometric, photometric and spectroscopic observations leading to the determination of the orbit and dynamical masses of the binary L dwarf 2MASSW J0746425+2000321. High angular resolution observations spread over almost 4 years and obtained with the Hubble Space Telescope (HST), the ESO Very Large Telescope (VLT), and a the W. M. Keck Observatory (Keck) allow us to cover 36% of the period, corresponding to 60% of the orbit, and, for the first time, to derive a precise estimate of the total and individual masses of such a late-type object. We find an orbital period of 3850.9$^{+904}_{-767}$ days. The corresponding total mass is 0.146$^{+0.016}_{-0.006}$ M$_{\sun}$, with uncertainties depending on the distance. Spatially resolved low resolution optical (550--1025 nm) spectra have been obtained with HST/STIS, allowing us to measure the spectral types of the two components (L0$\pm$0.5 for the primary and L1.5$\pm$0.5 for the secondary). We also present precise photometry of the individual components measured on the high angular resolution images obtained with HST/ACS and WFPC2 (visible), VLT/NACO (J, H and Ks bands) and Keck I (Ks) band). These spectral and photometric measurements enable us to estimate their effective temperatures and mass ratio, and to place the object accurately in a H--R diagram. The binary system is most likely formed by a primary with a mass of 0.085$\pm$0.010 M$_{\sun}$ and a secondary with a mass of 0.066$\pm$0.006 M$_{\sun}$, thus clearly substellar, for an age of approximately 300$\pm$150 Myr. H$\alpha$ variability indicates chromospheric and/or magnetic activity.

We present imaging and spectroscopic observations with HST (WFPC2, ACS/HRC and STIS), VLT (FORS2) and Keck (HIRES) of the dM8 ultra-cool dwarf DENIS-P J185950.9-370632, located in the R-CrA region. The presence of lithium absorption at 670.8 nm and the strong Halpha emission indicate a young age and a sub-stellar mass. Our diffraction-limited images resolve a companion at the separation limit of HST/ACS (~0.06"). The 2.1 mJy flux in the LW2 filter (5.0-8.5 microns) of the Infrared Space Observatory (see 1999A&A...350..883) likely corresponds to an infrared excess, suggesting the presence of circumstellar material. Proper motion and photometric measurements, as well as the Halpha activity, confirm membership in the R-CrA star forming region. If confirmed by further observations, DENIS-P J185950.9-370632 would be the first accreting sub-stellar multiple system observed to date.

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 present the results of a high-resolution imaging survey for brown dwarf (BD) binaries in two open clusters. The observations were carried out with WFPC2 onboard HST. Our sample consists of 8 BD candidates in the alpha Per cluster and 25 BD candidates in the Pleiades. We have resolved 4 binaries in the Pleiades with separations in the range 0".094--0".058, corresponding to projected separations between 11.7~AU and 7.2~AU. No binaries were found among the alpha Per targets. Three of the binaries have proper motions consistent with cluster membership in the Pleiades cluster, and for one of them we report the detection of Halpha in emission and LiI absorption obtained from Keck~II/ESI spectroscopy. One of the binaries does not have a proper motion consistent with Pleiades membership. We estimate that BD binaries wider than 12~AU are less frequent than 9% in the alphaPer and Pleiades clusters. This is consistent with an extension to substellar masses of a trend observed among stellar binaries: the maximum semimajor axis of binary systems decreases with decreasing primary mass. We find a binary frequency of 2 binaries over 13 BDs with confirmed proper motion membership in the Pleiades, corresponding to a binary fraction of 15%(1 sigma error bar +15%/-5%). These binaries are limited to the separation range 7-12~AU and their mass ratios are larger than 0.7. The relatively high binary frequency (>10%), the bias to separations smaller than about 15 AU and the trend to high mass ratios (q>0.7) are fundamental properties of BDs. Current theories of BD formation do not appear to provide a good description of all these properties.

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 present spectroscopy and time-series photometry of the dwarf nova QZ Ser. The spectrum shows a rich absorption line spectrum of type K4 +- 2. K-type secondary stars are generally seen in dwarf novae with orbital periods P-orb around 6 h, but in QZ Ser the absorption radial velocities show an obvious modulation (semi-amplitude 207(5) km/s) at P-orb = 119.752(2) min, much shorter than typical for such a relatively warm and prominent secondary spectrum. The H-alpha emission-line velocity is modulated at the same period and roughly opposite phase. Time-series photometry shows flickering superposed on a modulation with two humps per orbit, consistent with ellipsoidal variation of the secondary's light. QZ Ser is a second example of a relatively short-period dwarf nova with a surprisingly warm secondary. Model calculations suggest that the secondary is strongly enhanced in helium, and had already undergone significant nuclear evolution when mass transfer began. Several sodium absorption features in the secondary spectrum are unusually strong, which may indicate that the present-day surface was the site of CNO-cycle hydrogen burning in the past.

Based on SPH simulations, we quantify the geometrical distortion effect due to tidal and rotational forces on polytropic secondaries in semi-detached binaries. The main effect is an expansion of the polytropic star, with an effect on the radius of $\sim$ 5%-12%, depending on the polytropic index and the mass ratio. We apply such distortion effects to the secular evolution of secondaries in cataclysmic variable systems. We focus on systems below the 2-3h period gap and that approach the minimum period. We find a significant increase of the predicted minimum period ($\sim$ 4% if changes in the secondary's thermal relaxation are approximately taken into account). Though an improvement, the effect is not big enough to solve the mismatch between predicted and observed minimum period at 80 min.

We present spectroscopy and time-series photometry of the newly discovered dwarf nova 1RXS J232953.9+062814. Photometry in superoutburst reveals a superhump with a period of 66.06(6) minutes. The low state spectrum shows Balmer and HeI emission on a blue continuum, and in addition shows a rich absorption spectrum of type K4 +- 2. The absorption velocity is modulated sinusoidally at P_orb = 64.176(5) min, with semi-amplitude K = 348(4) km/s. The low-state light curve is double-humped at this period, and phased as expected for ellipsoidal variations. The absorption strength does not vary appreciably around the orbit. The orbital period is shorter than any other cataclysmic variable save for a handful of helium-star systems and V485 Centauri (59 minutes). The secondary is much hotter than main sequence stars of similar mass, but is well-matched by helium-enriched models, indicating that the secondary evolved from a more massive progenitor. A preliminary calculation in which a 1.2 solar-mass star begins mass transfer near the end of H burning matches this system's characteristics remarkably well.

While the modern stellar IMF shows a rapid decline with increasing mass, theoretical investigations suggest that very massive stars (>100 solar masses) may have been abundant in the early universe. Other calculations also indicate that, lacking metals, these same stars reach their late evolutionary stages without appreciable mass loss. After central helium burning, they encounter the electron-positron pair instability, collapse, and burn oxygen and silicon explosively. If sufficient energy is released by the burning, these stars explode as brilliant supernovae with energies up to 100 times that of an ordinary core collapse supernova. They also eject up to 50 solar masses of radioactive Ni56. Stars less massive than 140 solar masses or more massive than 260 solar masses should collapse into black holes instead of exploding, thus bounding the pair-creation supernovae with regions of stellar mass that are nucleosynthetically sterile. Pair-instability supernovae might be detectable in the near infrared out to redshifts of 20 or more and their ashes should leave a distinctive nucleosynthetic pattern.

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.

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 theoretical period - magnitude relationships for Cepheids in different filters for fundamental and first overtone pulsators, completing the work by Alibert et al. (1999). The results are provided for different metallicities characteristic of the Magellanic Clouds and the Milky Way. In contrast to the fundamental mode, we find a small metallicity effect on the period - luminosity relationship for the first overtone, due to the sensitivity of the period ratio $P_1$/$P_0$ with metallicity. Comparison is made with observations from OGLE and EROS in the Small and Large Magellanic Clouds. We emphasize the constraint on theoretical predictions provided by the combination of both fundamental and first overtone observed sequences. We obtain excellent agreement between models and data in a $\log P$ - $\wi$ (Wesenheit index) diagram for a distance modulus for the LMC $\mu_0$ = 18.60 - 18.70. We analyse the uncertainties of the fundamental period - magnitude relationships and the consequences on distance determination. We show that an arbitrary shift of the instability strip by 350 K in $\te$ yields up to 0.45 mag effect on $\mv$ at a given period, whereas the effect is less than 0.1 mag in the $K$-band. Using recent near-IR observations in the Large Magellanic Cloud and our $P$ - $\mk$ relationship, we derive a distance modulus for the LMC in agreement with the value based on $\wi$ data.

The stability of metal-free very massive stars ($Z$ = 0; $M = 120 - 500 \msol$) is analyzed and compared with metal-enriched stars. Such zero-metal stars are unstable to nuclear-powered radial pulsations on the main sequence, but the growth time scale for these instabilities is much longer than for their metal-rich counterparts. Since they stabilize quickly after evolving off the ZAMS, the pulsation may not have sufficient time to drive appreciable mass loss in Z = 0 stars. For reasonable assumptions regarding the efficiency of converting pulsational energy into mass loss, we find that, even for the larger masses considered, the star may die without losing a large fraction of its mass. We find a transition between the $\epsilon$- and $\kappa$-mechanisms for pulsational instability at $Z\sim 2\E{-4} - 2\E{-3}$. For the most metal-rich stars, the $\kappa$-mechanism yields much shorter $e$-folding times, indicating the presence of a strong instability. We thus stress the fundamental difference of the stability and late stages of evolution between very massive stars born in the early universe and those that might be born today.

We investigate the evolution, final fate, and nucleosynthetic yields of rotating and non-rotating very massive stars (VMS) of zero metallicity. First we address the issue of mass loss during hydrogen burning due to vibrational instabilities. We find that these objects are much more stable than what was found in previous studies of VMS of solar composition, and expect only negligible mass loss driven by the pulsations. As these stars thus reach the end of their evolution with massive helium cores, they encounter the pair-creation instability. We find that for helium core masses of ~64...133 solar masses these stars are completely disrupted with explosion energies of up to ~1E53 erg and eject up to ~60 solar masses of Ni56 Stars with more massive helium cores collapse into black holes. We present the first calculations that follow the collapse of such a massive rotating star and predict that X-ray burst and significant gravitational wave emission could result.

Using stellar models we investigate the relation between mass and the spectral type of the secondary star in low-mass short-period compact binaries such as cataclysmic variables and soft X-ray transients. Allowing for different mass transfer rates and different system ages prior to mass transfer we find that the secondaries should populate a band in the spectral type vs. mass plane. The mass M_ms of a ZAMS star with the same spectral type as the donor is effectively an upper limit to the donor mass. The lower mass limit for a given spectral type depends on the mixing length parameter. If this is large, there is no lower limit if the spectral type is later than K6. The band width decreases from 0.4 M_sun at K6 to less than 0.2 M_sun at K0. If the mixing length parameter is small, there is no lower mass limit if the spectral type is later than M2, and the band width decreases from 0.2 M_sun at M2 to less than 0.1 M_sun for types earlier than K0. We also point out an error in the method suggested by Beekman et al. (1997) to estimate the primary mass in a soft X-ray transient.

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

We investigate why the spectral type of most cataclysmic variable (CV) secondaries is significantly later than that of a ZAMS star with the same mean density. Using improved stellar input physics, tested against observations of low-mass stars at the bottom of the main sequence, we calculate the secular evolution of CVs with low-mass donors. We consider sequences with different mass transfer rates and with a different degree of nuclear evolution of the donor prior to mass transfer. Systems near the upper edge of the gap ($P \sim 3 - 6$ h) can be reproduced by models with a wide range of mass transfer rates from $1.5 \times 10^{-9} \msolyr$ to $10^{-8} \msolyr$. Evolutionary sequences with a small transfer rate and donors that are substantially evolved off the ZAMS (central hydrogen content $0.05-0.5$) reproduce CVs with late spectral types above $P \simgr$ 6 h. Systems with the most discrepant (late) spectral type should have the smallest donor mass at any given $P$. Consistency with the period gap suggests that the mass transfer rate increases with decreasing donor mass for evolved sequences above the period gap. In this case, a single-parameter family of sequences with varying $\xc$ and increasing mass transfer rate reproduces the full range of observed spectral types. This would imply that CVs with such evolved secondaries dominate the CV population.

We apply the new generation of theoretical models of low-mass stars to secondaries in CVs, focussing on systems above the period gap. The models confirm that the spectral type should be a good indicator of the donor mass. The orbital period-spectral type diagram can potentially constrain the long-term mean mass transfer rate. A transfer rate that increases with decreasing period is most easily reconciled with the observational data.

Using improved, up-to-date stellar input physics tested against observations of low-mass stars and brown dwarfs we calculate the secular evolution of low-mass donor cataclysmic variables (CVs), including those which form with a brown dwarf donor. Our models confirm the mismatch between the calculated minimum period (Pmin ~ 70 min) and the observed short-period cut-off (~ 80 min) in the CV period histogram. We find that tidal and rotational corrections applied to the one-dimensional stellar structure equations have no significant effect on the period minimum. Theoretical period distributions synthesized from our model sequences always show an accumulation of systems at the minimum period, a feature absent from the observed distribution. We suggest that non-magnetic CVs become unobservable as they are effectively trapped in permanent quiescence before they reach Pmin, and that small-number statistics may hide the period spike for magnetic CVs.

The visual surface brightness of K/M giants and dwarfs with near-solar metallicity differ slightly in agreement with the gravity effects predicted by recent theoretical models. We show that M-dwarfs display also a metallicity dependence of the surface brightness in the infrared $K$-band in agreement with theory. Based on these results, we present improved Barnes-Evans type relations and estimate the radii of 60 single or presumed M and K-dwarfs

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

Based on consistent evolutionary and pulsation calculations, we analyse the effect of metallicity and of different convection treatments in the stellar models on period - magnitude, - color and - radius relationships. In order to perform an accurate comparison with observations, we have computed grids of atmosphere models and synthetic spectra for different metallicities, covering the range of effective temperatures and gravities relevant for Cepheids. The models are compared to recent observations of galactic and Magellanic Clouds Cepheids. Unprecedented level of agreement is found between models and observations. We show that within the range of metallicity for the Galaxy and the Magellanic Clouds, a change of slope in the period - luminosity (PL) relationship is predicted at low periods, due to the reduction of the blue loop during core He burning. The minimum mass undergoing a blue loop and consequently the critical period at which this change of slope occurs depend on the metallicity Z and on the convection treatment in the stellar models. However, besides this change of slope, we do not find any significant effect of metallicity on period - magnitude relationships from V to K bands, and on period - color relationships in IR colors. We only find a detectable effect of $Z$ on $(B-V)$ colors. These results are not affected by uncertainties inherent to current stellar models, mainly due to convection treatment.

Hubble Space Telescope NICMOS observations are presented of six brown dwarf candidates in the Pleiades open cluster. One of them, namely CFHT-Pl-18, is clearly resolved as a binary with an angular separation of 0".33. The very low density of contaminating background stars in our images and the photometry of the components support that this system is a physical binary rather than a chance projection. All the available photometric and spectroscopic data indicate that the CFHT-Pl-18 system is likely a member of the Pleiades cluster, but a final confirmation will have to wait until lithium can be detected. Assuming cluster membership, we compare our NICMOS photometry with evolutionary models, and find that the inclusion of the effects of dust grains is necessary for fitting the data. We estimate that the masses of the components are about 0.045 Msol and 0.035 Msol. The binary system has a projected separation of 42 AU (for a distance of 125 pc) that is common among stellar binaries.

Using improved, up-to-date stellar input physics tested against observations of low-mass stars and brown dwarfs we calculate the secular evolution of low-mass donor CVs, including those which form with a brown dwarf donor star. Our models confirm the mismatch between the calculated minimum period (P_min ~ 70min) and the observed short-period cut-off (~80min) in the CV period histogram. Theoretical period distributions synthesized from our model sequences always show an accumulation of systems at the minimum period, a feature absent in the observed distribution. We suggest that non-magnetic CVs become unobservable as they are effectively trapped in permanent quiescence before they reach P_min, and that small-number statistics may hide the period spike for magnetic CVs. We calculate the minimum period for high mass transfer rate sequences and discuss the relevance of these for explaining the location of CV secondaries in the orbital period - spectral type diagram. We also show that a recently suggested revised mass-radius relation for low-mass main-sequence stars cannot explain the CV period gap.

We show that the secondaries in short-period cataclysmic variables with orbital periods $P < 3$ hr are close to the solar-abundance main sequence defined by single field stars. In cataclysmic variables with $P > 3$ hr, the earliest spectral types at a given period correspond to main sequence stars, while the majority of secondaries have later spectral types. Possible causes are nuclear evolution prior to mass transfer and lack of thermal equilibrium due to mass transfer. A comparison with evolutionary sequences obtained with up--to--date stellar models implies unusually high transfer rates and a large fraction of systems with evolved donors. There is no evidence for a secondary of low metallicity in any of the well-studied cataclysmic variables.

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.

We investigate the pulsational properties of RSG models --- which we evolve from ZAMS masses in the range 10 to $20 \Msun$ --- by means of linear and non-linear calculations. We find period and growth rate of the dominant fundamental mode to increase with increasing luminosity-to-mass ratio $L/M$. Our models obtain relatively large $L/M$ values due to the inclusion of rotation in the evolutionary calculations; however, the largest values are obtained at and beyond central He-exhaustion due to major internal rearrangements of the nuclear burning regions. Our non-linear calculations as well as the behavior of the linear period and growth rate of the pulsations for periods approaching the Kelvin-Helmholtz time scale of the H-rich stellar envelope point towards the possibility of large amplitude pulsations. Such properties are similar to that found in AGB stars and suggest the possibility of a ``superwind'' to occur before the RSGs explode as supernovae. We conclude that changes in global stellar properties during the last few $10^4 \yr$ before core collapse may lead to drastic changes in the pulsational and wind properties of pre-supernova stars, with marked consequences for the immediate pre-supernova structure of the star and the circumstellar medium. We compare our results with observations of long-period OH/IR variables and discuss observational evidence for our scenario from observed supernova light curves, spectra and remnants.

Over the past decade considerable effort, both observational and theoretical, has been directed towards a more accurate determination of the stellar lower main sequence and of the sub-stellar domain covered by Brown Dwarfs and Planets. A review of the recent improvement in the theory of cool model atmospheres is presented.

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