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 the 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 the Rossby number. The torque also contains an empirically derived scaling with stellar mass (and radius), which provides new insight into the mass-dependence of stellar magnetic and wind properties. We demonstrate that this new formulation explains why the lowest mass stars are observed to maintain rapid rotation for much longer than solar-mass stars, and simultaneously why older populations exhibit a sequence of slowly rotating stars, in which the low-mass stars rotate more slowly than solar-mass stars. The model also reproduces some previously unexplained features in the period-mass diagram for the Kepler field, notably: the particular shape of the "upper envelope" of the distribution, suggesting that ~95% of Kepler field stars with measured rotation periods are younger than ~4 Gyr; and the shape of the "lower envelope," corresponding to the location where stars transition between magnetically saturated and unsaturated regimes.
We explore the possibility that large-scale convection is inhibited over some regions of giant planet interiors, as a consequence of a gradient of composition inherited from either their formation ...history or 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 HD 209458b 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.
In this paper we present 3D 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 ∼10%. We also present the contribution function, diagnosed from our simulations, and show that wind-driven chemistry has a significant impact on its 3D 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 m Spitzer/IRAC channel.
The treatment of radiation transport in global circulation models (GCMs) is crucial for correctly describing 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 the 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.
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 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.
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 {approx}5%-10% lower and {approx}3%-5% higher than observed, respectively. These discrepancies have been attributed to different factors, notably 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 mainly distributed over the stellar poles. Assuming the existence of polar spots, overall agreement between models and observations is reached when {approx}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 {approx}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 the explanation of the differences seen in some of the systems with shorter orbital periods.
Recently, an analysis of Gaia Data Release 2 revealed a gap in the mid-M dwarf main sequence. The authors suggested the feature is linked to the onset of full convection in M dwarfs. Following the ...announcement of this discovery, an explanation has been proposed based on standard stellar evolution models. In this paper we re-examine this explanation. We confirm that nuclear burning and mixing processes of 3He 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 the details of the published explanation and propose an alternative.
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 ≤ M/M⊙ < 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 ≤ fos ≤ 0.014 (using the scheme from Herwig et al. 1997, A&A, 324, L81) are consistent with all eight systems studied. Generally, a much larger range of f 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.