Search

White dwarfs are a class of stars with unique physical properties. They present many challenging problems whose solution requires the application of advanced theories of dense matter, state-of-the-art experimental techniques, and extensive computing efforts. New ground- and space-based observatories will soon provide an increasingly detailed view of white dwarf stars and reveal new phenomena that will challenge our models. This review is an introduction for researchers who are not in the field of white dwarf astrophysics with the intent to entice them to contribute their expertise to advance our knowledge of these exotic stars. We discuss a wide variety of currently unsolved or partially resolved problems that are broadly related to equations of state, transport processes and opacities., Comment: Review article to appear in Physics Reports. 108 pages of text, 33 figures, 372 references

White dwarfs represent the final state of evolution for most stars1–3. Certain classes of white dwarfs pulsate4,5, leading to observable brightness variations, and analysis of these variations with theoretical stellar models probes their internal structure. Modelling of these pulsating stars provides stringent tests of white dwarf models and a detailed picture of the outcome of the late stages of stellar evolution6. However, the high-energy-density states that exist in white dwarfs are extremely difficult to reach and to measure in the laboratory, so theoretical predictions are largely untested at these conditions. Here we report measurements of the relationship between pressure and density along the principal shock Hugoniot (equations describing the state of the sample material before and after the passage of the shock derived from conservation laws) of hydrocarbon to within five per cent. The observed maximum compressibility is consistent with theoretical models that include detailed electronic structure. This is relevant for the equation of state of matter at pressures ranging from 100 million to 450 million atmospheres, where the understanding of white dwarf physics is sensitive to the equation of state and where models differ considerably. The measurements test these equation-of-state relations that are used in the modelling of white dwarfs and inertial confinement fusion experiments7,8, and we predict an increase in compressibility due to ionization of the inner-core orbitals of carbon. We also find that a detailed treatment of the electronic structure and the electron degeneracy pressure is required to capture the measured shape of the pressure–density evolution for hydrocarbon before peak compression. Our results illuminate the equation of state of the white dwarf envelope (the region surrounding the stellar core that contains partially ionized and partially degenerate non-ideal plasmas), which is a weak link in the constitutive physics informing the structure and evolution of white dwarf stars9. Researchers have measured the equation of state of hydrocarbon in a high-density regime, which is necessary for accurate modelling of the oscillations of white dwarf stars.

We present the results of the first Charged-Particle Transport Coefficient Code Comparison Workshop, which was held in Albuquerque, NM October 4-6, 2016. In this first workshop, scientists from eight institutions and four countries gathered to compare calculations of transport coefficients including thermal and electrical conduction, electron-ion coupling, inter-ion diffusion, ion viscosity, and charged particle stopping powers. Here, we give general background on Coulomb coupling and computational expense, review where some transport coefficients appear in hydrodynamic equations, and present the submitted data. Large variations are found when either the relevant Coulomb coupling parameter is large or computational expense causes difficulties. Understanding the general accuracy and uncertainty associated with such transport coefficients is important for quantifying errors in hydrodynamic simulations of inertial confinement fusion and high-energy density experiments., Comment: 45 pages, 17 figures

Exoplanet and brown dwarf atmospheres commonly show signs of disequilibrium chemistry. In the James Webb Space Telescope era high resolution spectra of directly imaged exoplanets will allow the characterization of their atmospheres in more detail, and allow systematic tests for the presence of chemical species that deviate from thermochemical equilibrium in these atmospheres. Constraining the presence of disequilibrium chemistry in these atmospheres as a function of parameters such as their effective temperature and surface gravity will allow us to place better constrains in the physics governing these atmospheres. This paper is part of a series of works presenting the Sonora grid of atmosphere models (Marley et al 2021, Morley et al in prep.). In this paper we present a grid of cloud-free, solar metallicity atmospheres for brown dwarfs and wide separation giant planets with key molecular species such as CH4, H2O, CO and NH3 in disequilibrium. Our grid covers atmospheres with Teff~[500 K,1300 K], logg~[3.0,5.5] (cgs) and an eddy diffusion parameter of logKzz=2, 4 and 7 (cgs). We study the effect of different parameters within the grid on the temperature and composition profiles of our atmospheres. We discuss their effect on the near-infrared colors of our model atmospheres and the detectability of CH4, H2O, CO and NH3 using the JWST. We compare our models against existing MKO and Spitzer observations of brown dwarfs and verify the importance of disequilibrium chemistry for T dwarf atmospheres. Finally, we discuss how our models can help constrain the vertical structure and chemical composition of these atmospheres., accepted to ApJ

The effect of ionic disorder on the principal Hugoniot is investigated using multiple scattering theory to very high pressure (Gbar). Calculations using molecular dynamics to simulate ionic disorder are compared to those with a fixed crystal lattice, for both carbon and aluminum. For the range of conditions considered here we find that ionic disorder has a relatively minor influence. It is most important at the onset of shell ionization and we find that, at higher pressures, the subtle effect of the ionic environment is overwhelmed by the larger number of ionized electrons with higher thermal energies.

We present a new generation of substellar atmosphere and evolution models, appropriate for application to studies of L, T, and Y-type brown dwarfs and self-luminous extrasolar planets. The atmosphere models describe the expected temperature-pressure profiles and emergent spectra of atmospheres in radiative-convective equilibrium with effective temperatures and gravities within the ranges $200\le T_{\rm eff}\le2400\,\rm K$ and $2.5\le \log g \le 5.5$. These ranges encompass masses from about 0.5 to 85 Jupiter masses for a set of metallicities ($[{\rm M/H}] = -0.5$ to $+0.5$), C/O ratios (from 0.5 to 1.5 times that of solar), and ages. The evolution tables describe the cooling of these substellar objects through time. These models expand the diversity of model atmospheres currently available, notably to cooler effective temperatures and greater ranges in C/O. Notable improvements from past such models include updated opacities and atmospheric chemistry. Here we describe our modeling approach and present our initial tranche of models for cloudless, chemical equilibrium atmospheres. We compare the modeled spectra, photometry, and evolution to various datasets., 27 pages, 16 figures, accepted for Astrophysical Journal. Models available at https://doi.org/10.5281/zenodo.5063476

We present the most detailed data-driven exploration of cloud opacity in a substellar object to-date. We have tested over 60 combinations of cloud composition and structure, particle size distribution, scattering model, and gas phase composition assumptions against archival $1-15 {\rm \mu m}$ spectroscopy for the unusually red L4.5~dwarf 2MASSW~J2224438-015852 using the Brewster retrieval framework. We find that, within our framework, a model that includes enstatite and quartz cloud layers at shallow pressures, combined with a deep iron cloud deck fits the data best. This models assumes a Hansen distribution for particle sizes for each cloud, and Mie scattering. We retrieved particle effective radii of $\log_{10} a {\rm (\mu m)} = -1.41^{+0.18}_{-0.17}$ for enstatite, $-0.44^{+0.04}_{-0.20}$ for quartz, and $-0.77^{+0.05}_{-0.06}$ for iron. Our inferred cloud column densities suggest ${\rm (Mg/Si)} = 0.69^{+0.06}_{-0.08}$ if there are no other sinks for magnesium or silicon. Models that include forsterite alongside, or in place of, these cloud species are strongly rejected in favour of the above combination. We estimate a radius of $0.75 \pm 0.02$ Rjup, which is considerably smaller than predicted by evolutionary models for a field age object with the luminosity of 2M2224-0158. Models which assume vertically constant gas fractions are consistently preferred over models that assume thermochemical equilibrium. From our retrieved gas fractions we infer ${\rm [M/H]} = +0.38^{+0.07}_{-0.06}$ and ${\rm C/O} = 0.83^{+0.06}_{-0.07}$. Both these values are towards the upper end of the stellar distribution in the Solar neighbourhood, and are mutually consistent in this context. A composition toward the extremes of the local distribution is consistent with this target being an outlier in the ultracool dwarf population., Comment: 20 pages, 10 figures, MNRAS accepted

Accurate modeling of the electronic structure of warm dense matter is a challenging problem whose solution would allow a better understanding of material properties like equation of state, opacity, and conductivity, with resulting applications from astrophysics to fusion energy research. Here we explore the real-space Green’s function method as a technique for solving the Kohn–Sham density functional theory equations under warm dense matter conditions. We find the method to be tractable and accurate throughout the density and temperature range of interest, in contrast to other approaches. Good agreement on equation of state is found when comparing to other methods, where they are thought to be accurate.

We present new near-infrared Gemini Planet Imager (GPI) spectroscopy of HD 206893 B, a substellar companion orbiting within the debris disk of its F5V star. The $J$, $H$, $K1$, and $K2$ spectra from GPI demonstrate the extraordinarily red colors of the object, confirming it as the reddest substellar object observed to date. The significant flux increase throughout the infrared presents a challenging atmosphere to model with existing grids. Best-fit values vary from 1200 K to 1800 K for effective temperature and from 3.0 to 5.0 for log($g$), depending on which individual wavelength band is fit and which model suite is applied. The extreme redness of the companion can be partially reconciled by invoking a high-altitude layer of sub-micron dust particles, similar to dereddening approaches applied to the peculiar red field L-dwarf population. However, reconciling the HD 206893 B spectra with even those of the reddest low-gravity L-dwarf spectra still requires the contribution of additional atmospheric dust, potentially due to the debris disk environment in which the companion resides. Orbit fitting from four years of astrometric monitoring is consistent with a $\sim$30-year period, orbital inclination of 147$^{\circ}$, and semimajor axis of 10 au, well within the estimated disk inner radius of $\sim$50 au. As one of very few substellar companions imaged interior to a circumstellar disk, the properties of this system offer important dynamical constraints on companion-disk interaction and provide a benchmark for substellar and planetary atmospheric study., Accepted for publication in AJ (October 15, 2020). 35 pages, 20 figures, with tables and appendices presented in their entirety

Thanks to their continuous cooling and relative simplicity, white dwarf stars are routinely used to measure the ages of stellar populations. The usefulness of white dwarfs as cosmochronometers depends on the availability of accurate cooling models. A key ingredient of those models are the conductive opacities, which largely govern the cooling rate. In this work, we present improved conductive opacities for the regime of moderate coupling and moderate degeneracy that characterizes an important portion of the envelopes of DA and DB white dwarfs. We find differences of up to a factor 3 between our calculations and the commonly used opacities of Cassisi et al. (2007), which we attribute to an improved account of electron-electron scattering. The cooling models are strongly affected by those changes in the conductive opacities: the age of a 4000 K white dwarf can be reduced by as much as 2 Gyr. We provide analytical fits to our new opacities to facilitate the implementation of this important effect in white dwarf evolution codes., Comment: 11 pages, 6 figures, 3 tables. Accepted for publication in The Astrophysical Journal

Rotational modulations of emission spectra in brown dwarf and exoplanet atmospheres show that clouds are often distributed non-uniformly in these ultracool atmospheres. The spatial heterogeneity in cloud distribution demonstrates the impact of atmospheric dynamics on cloud formation and evolution. In this study, we update the Hubble Space Telescope (HST) time-series data analysis of the previously reported rotational modulations of WISEP J004701+680352 -- an unusually red late-L brown dwarf with a spectrum similar to that of the directly imaged planet HR8799e. We construct a self-consistent spatially heterogeneous cloud model to explain the Hubble Space Telescope and the Spitzer time-series observations, as well as the time-averaged spectra of WISE0047. In the heterogeneous cloud model, a cloud thickness variation of around one pressure scale height explains the wavelength dependence in the HST near-IR spectral variability. By including disequilibrium CO/$CH_4$ chemistry, our models also reproduce the redder $J-K_{\rm s}$ color of WISE0047 compared to that of field brown dwarfs. We discuss the impact of vertical cloud structure on atmospheric profile and estimate the minimum eddy diffusivity coefficient for other objects with redder colors. Our data analysis and forward modeling results demonstrate that time-series spectrophotometry with a broad wavelength coverage is a powerful tool for constraining heterogeneous atmospheric structure., Comment: accepted for publication in The Astrophysical Journal

We have developed an efficient and versatile model to describe warm and hot dense matter that couples an average atom model with the integral fluid equations for the ion correlations. This model provides all bound and free electronic states and wave functions, the interaction potentials and all correlation functions without adjustable parameters. The electrons can be described quantum mechanically (Schrodinger equation) or semi-classically with the Thomas-Fermi model. The ion-ion pair potential can be used in a classical molecular dynamics simulation to yield equation of state and dynamic properties of dense plasmas, including mixtures. This "pseudo-atom molecular dynamics" model, or PAMD, is more approximate than ab initio methods but presents distinct computational advantages. Extensions of the model allow the calculation of diffusion coefficients, viscosity, X-ray Thomson Scattering spectra, DC conductivities and opacities. LA-UR-19-27814

We present $Spitzer$ IRS 5--14 $\mu$m spectra and 16 $\mu$m and 22 $\mu$m photometry of the T2.5 companion to the $\sim$300 Myr-old G0V star HN Peg. We incorporate previous 0.8--5 $\mu$m observations to obtain the most comprehensive spectral energy distribution of an intermediate-gravity L/T-transition dwarf which, together with an accurate Gaia EDR3 parallax of the primary, enable us to derive precise fundamental parameters. We find that young ($\approx$0.1--0.3 Gyr) early-T dwarfs on average have $\approx$140 K lower effective temperatures, $\approx$20% larger radii, and similar bolometric luminosities compared to $\gtrsim$1 Gyr-old field dwarfs with similar spectral types. Our accurate infrared spectrophotometry offers new detail at wavelengths where the dominant carbon-bearing molecules have their strongest transitions: at 3.4 $\mu$m for methane and at 4.6 $\mu$m for carbon monoxide. We assess the performance of various widely available photospheric models and find that models with condensates and/or clouds better reproduce the full SED of this moderately young early-T dwarf. However, cloud-free models incorporating a more general convective instability treatment reproduce at least the low-resolution near-IR spectrum similarly well. Our analysis of $R\approx2300$ $J$-band spectra shows that the near-infrared potassium absorption lines in HN Peg B have similar strengths to those seen in both younger and older T2-T3 dwarfs. We conclude that while alkali lines are well-established as surface gravity indicators for L-type or warmer stars, they are insensitive to surface gravity in early-T dwarfs, Comment: Accepted for publication in ApJ

We present a statistical analysis of the first 300 stars observed by the Gemini Planet Imager Exoplanet Survey (GPIES). This subsample includes six detected planets and three brown dwarfs; from these detections and our contrast curves we infer the underlying distributions of substellar companions with respect to their mass, semi-major axis, and host stellar mass. We uncover a strong correlation between planet occurrence rate and host star mass, with stars M $>$ 1.5 $M_\odot$ more likely to host planets with masses between 2-13 M$_{\rm Jup}$ and semi-major axes of 3-100 au at 99.92% confidence. We fit a double power-law model in planet mass (m) and semi-major axis (a) for planet populations around high-mass stars (M $>$ 1.5M$_\odot$) of the form $\frac{d^2 N}{dm da} \propto m^\alpha a^\beta$, finding $\alpha$ = -2.4 $\pm$ 0.8 and $\beta$ = -2.0 $\pm$ 0.5, and an integrated occurrence rate of $9^{+5}_{-4}$% between 5-13 M$_{\rm Jup}$ and 10-100 au. A significantly lower occurrence rate is obtained for brown dwarfs around all stars, with 0.8$^{+0.8}_{-0.5}$% of stars hosting a brown dwarf companion between 13-80 M$_{\rm Jup}$ and 10-100 au. Brown dwarfs also appear to be distributed differently in mass and semi-major axis compared to giant planets; whereas giant planets follow a bottom-heavy mass distribution and favor smaller semi-major axes, brown dwarfs exhibit just the opposite behaviors. Comparing to studies of short-period giant planets from the RV method, our results are consistent with a peak in occurrence of giant planets between ~1-10 au. We discuss how these trends, including the preference of giant planets for high-mass host stars, point to formation of giant planets by core/pebble accretion, and formation of brown dwarfs by gravitational instability., Comment: 52 pages, 18 figures. AJ in press

The precise astrometric measurements of the Gaia Data Release 2 have opened the door to detailed tests of the predictions of white dwarf cooling models. Significant discrepancies between theory and observations have been identified, the most striking affecting ultramassive white dwarfs. Cheng et al. (2019) found that a small fraction of white dwarfs on the so-called Q branch must experience an extra cooling delay of $\sim 8\,$Gyr not predicted by current models. $^{22}$Ne phase separation in a crystallizing C/O white dwarf can lead to a distillation process that efficiently transports $^{22}$Ne toward its center, thereby releasing a considerable amount of gravitational energy. Using state-of-the-art Monte Carlo simulations, we show that this mechanism can largely resolve the ultramassive cooling anomaly if the delayed population consists of white dwarfs with moderately above-average $^{22}$Ne abundances. We also argue that $^{22}$Ne phase separation can account for the smaller cooling delay currently missing for models of white dwarfs with more standard compositions., 8 pages, 3 figures, 1 table. Accepted for publication in ApJL on 2021-03-23

A new density functional theory, average atom based model for the electrical conductivity of dense plasmas with a mixture of ion species, containing no adjustable parameters, is presented. The model takes the temperature, mass density and relative abundances of the species as input. It takes into account partial ionization, ionic structure, and core-valence orthogonality, and uses quantum mechanical calculations of cross sections. Comparison to an existing high fidelity but computationally expensive method reveals good agreement. The new model is computationally efficient and can reach high temperatures. A new mixing rule is also presented that gives reasonably accurate conductivities for high temperature plasma mixtures.

The continuous cooling of a white dwarf is punctuated by events that affect its cooling rate. Probably the most significant of those is the crystallization of its core, a phase transition that occurs once the C/O interior has cooled down below a critical temperature. This transition releases latent heat as well as gravitational energy due to the redistribution of the C and O ions during solidification, thereby slowing down the evolution of the white dwarf. The unambiguous observational signature of core crystallization - a pile-up of objects in the cooling sequence - was recently reported. However, existing evolution models struggle to quantitatively reproduce this signature, casting doubt on their accuracy when used to measure the ages of stellar populations. The timing and amount of the energy released during crystallization depend on the exact form of the C/O phase diagram. Using the advanced Gibbs-Duhem integration method and state-of-the-art Monte Carlo simulations of the solid and liquid phases, we have obtained a very accurate version of this phase diagram, allowing a precise modeling of the phase transition. Despite this improvement, the magnitude of the crystallization pile-up remains underestimated by current evolution models. We conclude that latent heat release and O sedimentation alone are not sufficient to explain the observations and that other unaccounted physical mechanisms, possibly $^{22}$Ne phase separation, play an important role., 4 pages, 2 figures. Accepted for publication in A&A Letters

International audience; Material equation-of-state (EOS) models, generally providing the pressure and internal energy for a given density and temperature, are required to close the equations of hydrodynamics. As a result they are an essential piece of physics used to simulate inertial confinement fusion (ICF) implosions. Historically, EOS models based on different physical/chemical pictures of matter have been developed for ICF relevant materials such as the deuterium (D2) or deuterium-tritium (DT) fuel, as well as candidate ablator materials such as polystyrene (CH), glow-discharge polymer (GDP), beryllium (Be), carbon (C), and boron carbide (B4C). The accuracy of these EOS models can directly affect the reliability of ICF target design and understanding, as shock timing and material compressibility are essentially determined by what EOS models are used in ICF simulations. Systematic comparisons of current EOS models, benchmarking with experiments, not only help us to understand what the model differences are and why they occur, but also to identify the state-of-the-art EOS models for ICF target designers to use. For this purpose, the first Equation-of-State Workshop, supported by the US Department of Energy’s ICF program, was held at the Laboratory for Laser Energetics (LLE), University of Rochester on 31 May - 2nd June, 2017. This paper presents a detailed review on the findings from this workshop: (1) 5-10% model-model variations exist throughout the relevant parameter space, and can be much larger in regions where ionization and dissociation are occurring, (2) the D2 EOS is particularly uncertain, with no single model able to match the available experimental data, and this drives similar uncertainties in the CH EOS, and (3) new experimental capabilities such as Hugoniot measurements around 100 Mbar and high-quality temperature measurements are essential to reducing EOS uncertainty.

The binary brown dwarf WISE J104915.57$-$531906.1 (also Luhman 16AB), composed of a late L and early T dwarf, is a prototypical L/T transition flux reversal binary located at only 2 pc distance. Luhman 16B is a known variable whose light curves evolve rapidly. We present spatially resolved spectroscopic time-series of Luhman 16A and B covering 6.5 h using HST/WFC3 at 1.1 to 1.66 $\mu$m. The small, count-dependent variability of Luhman 16A at the beginning of the observations likely stems from instrumental systematics; Luhman 16A appears non-variable above $\approx$0.4%. Its spectrum is well fit by a single cloud layer with intermediate cloud thickness (f_sed=2, Teff=1200 K). Luhman 16B varies at all wavelengths with peak-to-valley amplitudes of 7-11%. The amplitude and light curve shape changes over only one rotation period. The lowest relative amplitude is found in the deep water absorption band at 1.4 $\mu$m, otherwise it mostly decreases gradually from the blue to the red edge of the spectrum. This is very similar to the other two known highly variable early T dwarfs. A two-component cloud model accounts for most of the variability, although small deviations are seen in the water absorption band. We fit the mean spectrum and relative amplitudes with a linear combination of two models of a warm, thinner cloud (Teff=1300 K, fsed=3) and a cooler, thicker cloud (Teff=1000-1100 K, f_sed=1), assuming out-of-equilibrium atmospheric chemistry. A cloud as for Luhman 16A but with holes cannot reproduce the variability of Luhman 16B, indicating more complex cloud evolution through the L/T transition. The projected separation of the binary has decreased by $\approx$0.3'' in 8 months., Comment: 14 pages, 12 figures. Accepted for publication in ApJ

We present spectro-photometry spanning 1-5 $\mu$m of 51 Eridani b, a 2-10 M$_\text{Jup}$ planet discovered by the Gemini Planet Imager Exoplanet Survey. In this study, we present new $K1$ (1.90-2.19 $\mu$m) and $K2$ (2.10-2.40 $\mu$m) spectra taken with the Gemini Planet Imager as well as an updated $L_P$ (3.76 $\mu$m) and new $M_S$ (4.67 $\mu$m) photometry from the NIRC2 Narrow camera. The new data were combined with $J$ (1.13-1.35 $\mu$m) and $H$ (1.50-1.80 $\mu$m) spectra from the discovery epoch with the goal of better characterizing the planet properties. 51 Eri b photometry is redder than field brown dwarfs as well as known young T-dwarfs with similar spectral type (between T4-T8) and we propose that 51 Eri b might be in the process of undergoing the transition from L-type to T-type. We used two complementary atmosphere model grids including either deep iron/silicate clouds or sulfide/salt clouds in the photosphere, spanning a range of cloud properties, including fully cloudy, cloud free and patchy/intermediate opacity clouds. Model fits suggest that 51 Eri b has an effective temperature ranging between 605-737 K, a solar metallicity, a surface gravity of $\log$(g) = 3.5-4.0 dex, and the atmosphere requires a patchy cloud atmosphere to model the SED. From the model atmospheres, we infer a luminosity for the planet of -5.83 to -5.93 ($\log L/L_{\odot}$), leaving 51 Eri b in the unique position as being one of the only directly imaged planet consistent with having formed via cold-start scenario. Comparisons of the planet SED against warm-start models indicates that the planet luminosity is best reproduced by a planet formed via core accretion with a core mass between 15 and 127 M$_{\oplus}$., Comment: 27 pages, 19 figures, Accepted for publication in The Astronomical Journal

We present the first results from applying the spectral inversion technique in the cloudy L dwarf regime. Our new framework provides a flexible approach to modelling cloud opacity which can be built incrementally as the data requires, and improves upon previous retrieval experiments in the brown dwarf regime by allowing for scattering in two stream radiative transfer. Our first application of the tool to two mid-L dwarfs is able to reproduce their near-infrared spectra far more closely than grid models. Our retrieved thermal, chemical, and cloud profiles allow us to estimate $T_{\rm eff} = 1796^{+23}_{-25}$ K and $\log g = 5.21^{+0.05}_{-0.08}$ for 2MASS J05002100+0330501 and for 2MASSW J2224438-015852 we find $T_{\rm eff} = 1723^{+18}_{-19}$ K and $\log g = 5.31^{+0.04}_{-0.08}$, in close agreement with previous empirical estimates. Our best model for both objects includes an optically thick cloud deck which passes $\tau_{cloud} \geq 1$ (looking down) at a pressure of around 5 bar. The temperature at this pressure is too high for silicate species to condense, and we argue that corundum and/or iron clouds are responsible for this cloud opacity. Our retrieved profiles are cooler at depth, and warmer at altitude than the forward grid models that we compare, and we argue that some form of heating mechanism may be at work in the upper atmospheres of these L dwarfs. We also identify anomalously high CO abundance in both targets, which does not correlate with the warmth of our upper atmospheres or our choice of cloud model, and find similarly anomalous alkali abundance for one of our targets. These anomalies may reflect unrecognised shortcomings in our retrieval model, or inaccuracies in our gas phase opacities., Comment: Accepted for publication in MNRAS. Changes in review: additional validation against T dwarf test case, posteriors for cloud free retrievals and contribution functions added to Appendix. Also, entire investigation re-run with new UCL H2O opacities, with minimal impact on results

A model for dense homo-nuclear plasmas that couples an average atom model for the calculation of the electronic structure to the quantum Ornstein–Zernike equations describing the ionic structure is summarized and described pedagogically. The model is applied to the calculation of ion–ion pair distribution functions g II ( r ) for tungsten in the warm and hot dense matter regimes. These results are compared to orbital-free molecular dynamics simulations and excellent agreement is found. Calculations of g II ( r ) with a simple version of the model (which we call the ion-sphere model) are in remarkable agreement with those of the full model. This ion-sphere model provides a simple and efficient method of calculating accurate g II ( r ) for warm and hot dense matter for many applications involving low- to high- Z elements with a modest investment of effort.

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

We present an approximation for calculating the equation of state (EOS) of warm and hot dense matter that is built on the previously published pseudoatom molecular dynamics (PAMD) model of dense plasmas [Starrett et al., Phys. Rev. E 91, 013104 (2015)PLEEE81539-375510.1103/PhysRevE.91.013104]. While the EOS calculation with PAMD was previously limited to orbital-free density functional theory (DFT), the new approximation presented here allows a Kohn-Sham DFT treatment of the electrons. The resulting EOS thus includes a quantum mechanical treatment of the electrons with a self-consistent model of the ionic structure, while remaining tractable at high temperatures. The method is validated by comparisons with pressures from ab initio simulations of Be, Al, Si, and Fe. The EOS in the Thomas-Fermi approximation shows remarkable thermodynamic consistency over a wide range of temperatures for aluminum. We calculate the principal Hugoniots of aluminum and silicon up to 500 eV. We find that the ionic structure of the plasma has a modest effect that peaks at temperatures of a few eV and that the features arising from the electronic structure agree well with ab initio simulations.

Modeling warm dense matter, where a combination of partial ionization, partial electron degeneracy, and strong ion–ion and ion–electron coupling occur, is a frontier of equation of state research. We present the quantum hypernetted chain model which can be applied to studies of liquid metals, warm dense matter, and plasmas. This is an all-electron model that considers a mixture of a classical fluid of ions (with bound electrons) and a quantum electron fluid. The model describes self-consistently the structure of the ion fluid as well as the bound states of the ions and the non-linear response of the electron fluid. We present our initial results and compare them with experimental and ab initio results for liquid metals and low-temperature plasmas.

A summary of a variational average atom model which is used to close the quantum Ornstein–Zernike relations is presented. The first numerical results are presented from a code developed to solve these equations with two simplifying assumptions, described in the text. The inputs are the nuclear charge of the plasma species, the average material density and the temperature. Results generated include the spatially dependent electronic and nuclear densities, which are related to the electron-nucleus and nucleus–nucleus pair distribution functions. These in turn are simply related to the static structure factors. Numerical results presented are in the form of average ionization and nucleus–nucleus pair distribution functions. Average ionizations for carbon are compared to experiment and other models, showing good agreement.

We present near-infrared photometry and spectroscopy, and warm-Spitzer IRAC photometry of the young very cool T dwarf Ross 458C, which we have typed as T8.5p. By applying the fiducial age constraints (< 1Gyr) imposed by the properties of the active M dwarf Ross 458A, we have used these data to determine that Ross 458C has Teff = 695 \pm 60K, log g = 4.0 - 4.7 and an inferred mass of 5-20 MJup. We have compared fits of the near-infrared spectrum and IRAC photometry to the BT Settl and Saumon & Marley model grids, and have found that both sets provide best fits that are consistent with our derived properties, whilst the former provide a marginally closer match to the data for all scenarios explored here. The main difference between the model grids arises in the 4.5{\mu}m region, where the BT Settl models are able to better predict the flux through the IRAC filter, suggesting that non-equilibrium effects on the CO-CO2 ratio are important for shaping the mid-infrared spectra of very cool T dwarfs. We have also revisited the issue of dust opacity in the spectra of Ross 458C that was raised by Burgasser et al (2010). We have found that the BT Settl models which also incorporate a condensate cloud model, provide a better match to the near-infrared spectrum of this target than the Saumon & Marley model with fsed = 2, and we briefly discuss the influence of condensate clouds on T dwarf spectra.

We present an analysis of the recently discovered blue L dwarf SDSS J141624.08+134826.7. We extend the spectral coverage of its published spectrum to ~4 microns by obtaining a low-resolution L band spectrum with SpeX on the NASA IRTF. The spectrum exhibits a tentative weak CH4 absorption feature at 3.3 microns but is otherwise featureless. We derive the atmospheric parameters of SDSS J141624.08+134826.7 by comparing its 0.7-4.0 micron spectrum to the atmospheric models of Marley and Saumon which include the effects of both condensate cloud formation and non-equilibrium chemistry due to vertical mixing and find the best fitting model has Teff=1700 K, log g=5.5 [cm s-2], fsed=4, and Kzz=10^4 cm2 s-1. The derived effective temperature is significantly cooler than previously estimated but we confirm the suggestion by Bowler et al. that the peculiar spectrum of SDSS J141624.08+134826.7 is primarily a result of thin condensate clouds. In addition, we find strong evidence of vertical mixing in the atmosphere of SDSS J141624.08+134826.7 based on the absence of the deep 3.3 micron CH4 absorption band predicted by models computed in chemical equilibrium. This result suggests that observations of blue L dwarfs are an appealing way to quantitatively estimate the vigor of mixing in the atmospheres of L dwarfs because of the dramatic impact such mixing has on the strength of the 3.3 micron CH4 band in the emergent spectra of L dwarfs with thin condensate clouds., Accepted for publication in the Astronomical Journal

We present 7.5-14.2um low-resolution spectroscopy, obtained with the Spitzer Infrared Spectrograph, of the T8.5 dwarf Wolf 940 B, which is a companion to an M4 dwarf with a projected separation of 400 AU. We combine these data with previously published near-infrared spectroscopy and mid-infrared photometry, to produce the spectral energy distribution for the very low-temperature T dwarf. We use atmospheric models to derive the bolometric correction and obtain a luminosity of log L/Lsun = -6.01 +/- 0.05. Evolutionary models are used with the luminosity to constrain the values of effective temperature (T_eff) and surface gravity, and hence mass and age for the T dwarf. We further restrict the allowed range of T_eff and gravity using age constraints implied by the M dwarf primary, and refine the physical properties of the T dwarf by comparison of the observed and modelled spectroscopy and photometry. This comparison indicates that Wolf 940 B has a metallicity within 0.2 dex of solar, as more extreme values give poor fits to the data - lower metallicity produces a poor fit at lambda > 2um while higher metallicity produces a poor fit at lambda < 2um. This is consistent with the independently derived value of [m/H] = +0.24 +/- 0.09 for the primary star, using the Johnson & Apps (2008) M_K:V-K relationship. We find that the T dwarf atmosphere is undergoing vigorous mixing, with an eddy diffusion coefficient K_zz of 10^4 to 10^6 cm^2 s^-1. We derive an effective temperature of 585 K to 625 K, and surface gravity log g = 4.83 to 5.22 (cm s^-2), for an age range of 3 Gyr to 10 Gyr, as implied by the kinematic and H alpha properties of the M dwarf primary. The lower gravity corresponds to the lower temperature and younger age for the system, and the higher value to the higher temperature and older age. The mass of the T dwarf is 24 M_Jupiter to 45 M_Jupiter for the younger to older age limit., 24 pages which include 5 Figures and 3 Tables. Accepted for publication in the Astrophysical Journal July 2 2010

We briefly review the close connection between research on the equation of state of dense hydrogen and helium and models of the structure of Jupiter, Saturn and exoplanets and how they inform the process of planet formation. Recent progress in equation of state experiments and theory is high lighted.

Determinations of the luminosity and atmospheric properties of the T8 brown dwarf 2MASS J09393548 - 2448279 are presented, based on Spitzer IRAC and IRS observations and ground-based astrometry. We find log_(10)(L_(bol)/L⊙) = -5.69 ± 0.03 for this source, comparable to the current low-luminosity record holder 2MASS J04151954 - 0935066. However, modeling of near- and mid-infrared spectral data indicates an effective temperature of 600 ± 35 K, roughly 100 K cooler than 2M0415. These parameters require a highly inflated radius for 2M0939 (R ≈ 0.13 R⊙) which cannot be reconciled with brown dwarf structure models. However, if this source is an unresolved, equal-mass binary, then the reduced luminosity of each component (L_(bol) ≈ 10^-6 L⊙) can be brought into agreement with the inferred atmospheric parameters for an age of 0.4-12 Gyr and component masses of 0.01-0.05 M⊙. This hypothesis can be tested through future high-resolution imaging and/or spectroscopic observations.

We present an analysis of the 0.95-14.5 micron spectral energy distributions of nine field ultracool dwarfs with spectral types ranging from L1 to T4.5. Effective temperatures, gravities, and condensate cloud sedimentation efficiencies are derived by comparing the data to synthetic spectra computed from atmospheric models that self-consistently include the formation of condensate clouds. Derived effective temperatures decrease steadily through the L1 to T4.5 spectral types and we confirm that the effective temperatures of ultracool dwarfs at the L/T transition are nearly constant, decreasing by only ~200 K from spectral types L7.5 to T4.5. The two objects in our sample with very red J-Ks colors are best fitted with synthetic spectra that have thick clouds which hints at a possible correlation between the near-infrared colors of L dwarfs and the condensate cloud properties. The fits to the two T dwarfs in our sample (T2 and T4.5) also suggest that the clouds become thinner in this spectral class, in agreement with previous studies. Restricting the fits to narrower wavelength ranges (i.e., individual photometric bands) almost always yields excellent agreement between the data and models. Limitations in our knowledge of the opacities of key absorbers such as FeH, VO, and CH4 at certain wavelengths remain obvious, however. The effective temperatures obtained by fitting the narrower wavelength ranges can show a large scatter compared to the values derived by fitting the full spectral energy distributions; deviations are typically ~200 K and in the worst cases, up to 700 K.

We present 10 - 19 um moderate resolution spectra of ten M dwarfs, one L dwarf, and two T dwarf systems obtained with the Infrared Spectrograph (IRS) onboard the Spitzer Space Telescope. The IRS allows us to examine molecular spectroscopic features/lines at moderate spectral resolution in a heretofore untapped wavelength regime. These R~600 spectra allow for a more detailed examination of clouds, non-equilibrium chemistry, as well as the molecular features of H2O, NH3, and other trace molecular species that are the hallmarks of these objects. A cloud-free model best fits our mid-infrared spectrum of the T1 dwarf epsilon Indi Ba, and we find that the NH3 feature in epsilon Indi Bb is best explained by a non-equilibrium abundance due to vertical transport in its atmosphere. We examined a set of objects (mostly M dwarfs) in multiple systems to look for evidence of emission features, which might indicate an atmospheric temperature inversion, as well as trace molecular species; however, we found no evidence of either., Comment: 19 pages, 7 figures, accepted ApJ 1/12/07

We present new L' (3.75 um) photometry of six L and T dwarfs, and M' (4.70 um) photometry of ten L and T dwarfs, observed at Gemini Observatory, and new 3.55, 4.49, 5.73 and 7.87 um photometry of nine L and T dwarfs, obtained with the Spitzer Space Telescope. The sample includes unusually blue and red dwarfs from our near-infrared studies. The data are combined with published L', M' and Spitzer photometry of L and T dwarfs, and trends of colors with spectral type and other colors are examined. Model atmospheres by Marley and Saumon are used to generate synthetic colors for ranges of effective temperature, gravity, grain sedimentation efficiency, metallicity and vertical mixing efficiency. We explore how these parameters affect the mid-infrared colors of L and T dwarfs and find that the data are modelled satisfactorily only if substantial vertical mixing occurs in both L- and T-dwarf atmospheres. The location and range of the L and T dwarf sequences in IRAC color-color and color - magnitude diagrams is also only reproduced if this mixing occurs, with a range of efficiency described by K_zz ~ 10^2-10^6 cm^2/s. The colors of the unusually red dwarfs are best reproduced by non-equilibrium models with low sedimentation efficiency, i.e. thick cloud decks, and those of the unusually blue dwarfs by non-equilibrium models with high sedimentation efficiency, i.e. thin cloud decks. The K-L' and Spitzer [3.55]-[4.49] colors can be used as indicators of effective temperature for L and T dwarfs, but care must be taken to include gravity and metallicity effects for late-T dwarfs and vertical mixing for both late-L and T dwarfs., Accepted for publication in the Astrophysical Journal; 21 pages, 8 Figures, 3 Tables

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

We present a theoretical model that allows a fast and accurate evaluation of ionic transport properties of realistic plasmas spanning from warm and dense to hot and dilute conditions, including mixtures. This is achieved by combining a recent kinetic theory based on effective interaction potentials with a model for the equilibrium radial density distribution based on an average atom model and the integral equations theory of fluids. The model should find broad use in applications where nonideal plasma conditions are traversed, including inertial confinement fusion, compact astrophysical objects, solar and extrasolar planets, and numerous present-day high energy density laboratory experiments.

Directly detecting thermal emission from young extrasolar planets allows measurement of their atmospheric composition and luminosity, which is influenced by their formation mechanism. Using the Gemini Planet Imager, we discovered a planet orbiting the \$sim$20 Myr-old star 51 Eridani at a projected separation of 13 astronomical units. Near-infrared observations show a spectrum with strong methane and water vapor absorption. Modeling of the spectra and photometry yields a luminosity of L/LS=1.6-4.0 x 10-6 and an effective temperature of 600-750 K. For this age and luminosity, "hot-start" formation models indicate a mass twice that of Jupiter. This planet also has a sufficiently low luminosity to be consistent with the "cold- start" core accretion process that may have formed Jupiter., Comment: 29 pages, 3 figures, 2 tables, and Supplementary Materials. published in Science Express on Aug 13 2015. List of authors and the magnitudes of the star were correted

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

The elastic feature of x-ray scattering from warm dense aluminum has recently been measured by Fletcher et al. [Nature Photonics 9, 274 (2015)]10.1038/nphoton.2015.41 with much higher accuracy than had hitherto been possible. This measurement is a direct test of the ionic structure predicted by models of warm dense matter. We use the method of pseudoatom molecular dynamics to predict this elastic feature for warm dense aluminum with temperatures of 1-100 eV and densities of 2.7-8.1g/cm^{3}. We compare these predictions to experiments, finding good agreement with Fletcher et al. and corroborating the discrepancy found in analyses of an earlier experiment of Ma et al. [Phys. Rev. Lett. 110, 065001 (2013)]PRLTAO0031-900710.1103/PhysRevLett.110.065001. We also evaluate the validity of the Thomas-Fermi model of the electrons and of the hypernetted chain approximation in computing the elastic feature and find them both wanting in the regime currently probed by experiments.

The ion-ion dynamic structure factor of warm dense matter is determined using the recently developed pseudoatom molecular dynamics method [Starrett et al., Phys. Rev. E 91, 013104 (2015)]. The method uses density functional theory to determine ion-ion pair interaction potentials that have no free parameters. These potentials are used in classical molecular dynamics simulations. This constitutes a computationally efficient and realistic model of dense plasmas. Comparison with recently published simulations of the ion-ion dynamic structure factor and sound speed of warm dense aluminum finds good to reasonable agreement. Using this method, we make predictions of the ion-ion dynamical structure factor and sound speed of a warm dense mixture-equimolar carbon-hydrogen. This material is commonly used as an ablator in inertial confinement fusion capsules, and our results are amenable to direct experimental measurement.

We present the first mid-infrared spectra of brown dwarfs, together with observations of a low-mass star. Our targets are the M3.5 dwarf GJ 1001A, the L8 dwarf DENIS-P J0255-4700, and the T1/T6 binary system epsilon Indi Ba/Bb. As expected, the mid-infrared spectral morphology of these objects changes rapidly with spectral class due to the changes in atmospheric chemistry resulting from their differing effective temperatures and atmospheric structures. By taking advantage of the unprecedented sensitivity of the Infrared Spectrograph on the Spitzer Space Telescope we have detected the 7.8 micron methane and 10 micron ammonia bands for the first time in brown dwarf spectra.

We present the first analysis of an optical to mid-infrared spectrum of the T7.5 dwarf Gliese 570D with model atmospheres, synthetic spectra, and brown dwarf evolution sequences. We obtain precise values for the basic parameters of Gl 570D: Teff=800 - 820K, log g (cm/s^2)=5.09 - 5.23, and log L/Lsun= -5.525 to -5.551. The Spitzer IRS spectrum shows prominent features of ammonia (NH3) that can only be fitted by reducing the abundance of NH3 by about one order of magnitude from the value obtained with chemical equilibrium models. We model departures from chemical equilibrium in the atmosphere of Gl 570D by considering the kinetics of nitrogen and carbon chemistry in the presence of vertical mixing. The resulting model spectrum reproduces the data very well., Comment: Accepted for publication in the ApJ. 10 pages, including 3 figures

We investigate the atmosphere and interior of the new transiting planet HD 149026b, which appears to be very rich in heavy elements. We first compute model atmospheres at metallicities ranging from solar to ten times solar, and show how for cases with high metallicity or inefficient redistribution of energy from the day side, the planet may develop a hot stratosphere due to absorption of stellar flux by TiO and VO. The spectra predicted by these models are very different than cooler atmosphere models without stratospheres. The spectral effects are potentially detectable with the Spitzer Space Telescope. In addition the models with hot stratospheres lead to a large limb brightening, rather than darkening. We compare the atmosphere of HD 149026b to other well-known transiting planets, including the recently discovered HD 189733b, which we show have planet-to-star flux ratios twice that of HD 209458 and TrES-1. The methane abundance in the atmosphere of HD 189733b is a sensitive indicator of atmospheric temperature and metallicity and can be constrained with Spitzer IRAC observations. We then turn to interior studies of HD 149026b and use a grid of self-consistent model atmospheres and high-pressure equations of state for all components to compute thermal evolution models of the planet. We estimate that the mass of heavy elements within the planet is in the range of 60 to 93 M_earth. Finally, we discuss trends in the radii of transiting planets with metallicity in light of this new member of the class., Comment: Accepted to the Astrophysical Journal. 18 pages, including 10 figures. New section on the atmosphere of planet HD 189733b. Enhanced discussion of atmospheric Ti chemistry and core mass for HD 149026b

We briefly examine the properties of dense plasmas characteristic of the interior of giant planets and the atmospheres of neutron stars. Special attention is devoted to the equation of state of hydrogen and helium at high density and to the effect of magnetic fields on the properties of dense matter., Invited Review, Strongly Coupled Coulomb Systems, Moscow June 2005; to appear in Journal of Physics A