skip to main content

Title: Cloud busting: enstatite and quartz clouds in the atmosphere of 2M2224-0158
ABSTRACT 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 μ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 model 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 ± 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 towards the extremes of the local distribution is consistent with this target being an outlier in the ultracool dwarf population.  more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range / eLocation ID:
1944 to 1961
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We present a reanalysis of the K2-106 transiting planetary system, with a focus on the composition of K2-106b, an ultra-short-period, super-Mercury candidate. We globally model existing photometric and radial velocity data and derive a planetary mass and radius for K2-106b of Mp = 8.53 ± 1.02 M⊕ and = - + Rp 1.71 0.057 RÅ 0.069 , which leads to a density of r = - + 9.4 p 1.5 1.6 g cm−3 , a significantly lower value than previously reported in the literature. We use planet interior models that assume a two-layer planet comprised of a liquid, pure Fe core and an iron-free, MgSiO3 mantle, and we determine that the range of the core mass fractions are consistent with the observed mass and radius. We use existing high-resolution spectra of the host star to derive the Fe/Mg/Si abundances ([Fe/ H] = −0.03 ± 0.01, [Mg/H] = 0.04 ± 0.02, [Si/H] = 0.03 ± 0.06) to infer the composition of K2-106b. We find that K2-106b has a density and core mass fraction ( - + 44 %15 12 ) consistent with that of Earth (CMF⊕ = 32%). Furthermore, its composition is consistent with what is expected, assuming that it reflects the relative refractory abundances of its host star. K2-106b is therefore unlikely to be a super-Mercury, as has been suggested in previous literature. 
    more » « less

    We present new [${\rm O\, {\small III}}$] 88-$\mu \mathrm{{m}}$ observations of five bright z ∼ 7 Lyman-break galaxies spectroscopically confirmed by ALMA through [${\rm C\, {\small II}}$] 158 $\mu \mathrm{{m}}$, unlike recent [${\rm O\, {\small III}}$] detections where Lyman α was used. This nearly doubles the sample of Epoch of Reionization galaxies with robust (5σ) [${\rm C\, {\small II}}$] and [${\rm O\, {\small III}}$] detections. We perform a multiwavelength comparison with new deep HST images of the rest-frame UV, whose compact morphology aligns well with [${\rm O\, {\small III}}$] tracing ionized gas. In contrast, we find more spatially extended [${\rm C\, {\small II}}$] emission likely produced in neutral gas, as indicated by an [${\rm N\, {\small II}}$] 205-$\mu \mathrm{{m}}$ non-detection in one source. We find a correlation between the optical ${[{\rm O\, {\small III}}]}+ {\mathrm{H\,\beta }}$ equivalent width and [${\rm O\, {\small III}}$]/[${\rm C\, {\small II}}$], as seen in local metal-poor dwarf galaxies. cloudy models of a nebula of typical density harbouring a young stellar population with a high-ionization parameter adequately reproduce the observed lines. Surprisingly, however, our models fail to reproduce the strength of [${\rm O\, {\small III}}$] 88-$\mu \mathrm{{m}}$, unless we assume an α/Fe enhancement and near-solar nebular oxygen abundance. On spatially resolved scales, we find [${\rm O\, {\small III}}$]/[${\rm C\, {\small II}}$] shows a tentative anticorrelation with infrared excess, LIR/LUV, also seen on global scales in the local Universe. Finally, we introduce the far-infrared spectral energy distribution fitting code mercurius to show that dust-continuum measurements of one source appear to favour a low dust temperature and correspondingly high dust mass. This implies a high stellar metallicity yield and may point towards the need of dust production or grain-growth mechanisms beyond supernovae.

    more » « less

    Interstellar chemistry is important for galaxy formation, as it determines the rate at which gas can cool, and enables us to make predictions for observable spectroscopic lines from ions and molecules. We explore two central aspects of modelling the chemistry of the interstellar medium (ISM): (1) the effects of local stellar radiation, which ionizes and heats the gas, and (2) the depletion of metals on to dust grains, which reduces the abundance of metals in the gas phase. We run high-resolution (400 M⊙ per baryonic particle) simulations of isolated disc galaxies, from dwarfs to Milky Way-mass, using the fire galaxy formation models together with the chimes non-equilibrium chemistry and cooling module. In our fiducial model, we couple the chemistry to the stellar fluxes calculated from star particles using an approximate radiative transfer scheme; and we implement an empirical density-dependent prescription for metal depletion. For comparison, we also run simulations with a spatially uniform radiation field, and without metal depletion. Our fiducial model broadly reproduces observed trends in H i and H2 mass with stellar mass, and in line luminosity versus star formation rate for [C ii]$_{158 \rm {\mu m}}$, [O i]$_{63 \rm {\mu m}}$, [O iii]$_{88 \rm {\mu m}}$, [N ii]$_{122 \rm {\mu m}}$, and H α6563Å. Our simulations with a uniform radiation field predict fainter luminosities, by up to an order of magnitude for [O iii]$_{88 \rm {\mu m}}$ and H α6563Å, while ignoring metal depletion increases the luminosity of carbon and oxygen lines by a factor ≈ 2. However, the overall evolution of the galaxy is not strongly affected by local stellar fluxes or metal depletion, except in dwarf galaxies where the inclusion of local fluxes leads to weaker outflows and hence higher gas fractions.

    more » « less
  4. null (Ed.)
    ABSTRACT Advanced LIGO and Advanced Virgo are detecting a large number of binary stellar origin black hole (BH) mergers. A promising channel for accelerated BH merger lies in active galactic nucleus (AGN) discs of gas around supermasssive BHs. Here, we investigate the relative number of compact object (CO) mergers in AGN disc models, including BH, neutron stars (NS), and white dwarfs, via Monte Carlo simulations. We find the number of all merger types in the bulk disc grows ∝ t1/3 which is driven by the Hill sphere of the more massive merger component. Median mass ratios of NS–BH mergers in AGN discs are $\tilde{q}=0.07\pm 0.06(0.14\pm 0.07)$ for mass functions (MF) M−1(− 2). If a fraction fAGN of the observed rate of BH–BH mergers (RBH–BH) come from AGN, the rate of NS–BH (NS–NS) mergers in the AGN channel is ${R}_{\mathrm{ BH}\!-\!\mathrm{ NS}} \sim f_{\mathrm{ AGN}}[10,300]\, \rm {Gpc}^{-3}\, \rm {yr}^{-1},({\mathit{ R}}_{NS\!-\!NS} \le \mathit{ f}_{AGN}400\, \rm {Gpc}^{-3}\, \rm {yr}^{-1}$). Given the ratio of NS–NS/BH–BH LIGO search volumes, from preliminary O3 results the AGN channel is not the dominant contribution to observed NS–NS mergers. The number of lower mass gap events expected is a strong function of the nuclear MF and mass segregation efficiency. CO merger ratios derived from LIGO can restrict models of MF, mass segregation, and populations embedded in AGN discs. The expected number of electromagnetic (EM) counterparts to NS–BH mergers in AGN discs at z < 1 is $\sim [30,900]\, {\rm {yr}}^{-1}(f_{\mathrm{ AGN}}/0.1)$. EM searches for flaring events in large AGN surveys will complement LIGO constraints on AGN models and the embedded populations that must live in them. 
    more » « less
  5. Abstract

    We present an atmospheric retrieval analysis of a pair of highly variable, ∼200 Myr old, early T type planetary-mass exoplanet analogs SIMP J01365662+0933473 and 2MASS J21392676+0220226 using the Brewster retrieval framework. Our analysis, which makes use of archival 1–15μm spectra, finds almost identical atmospheres for both objects. For both targets, we find that the data is best described by a patchy, high-altitude forsterite (Mg2SiO4) cloud above a deeper, optically thick iron (Fe) cloud. Our model constrains the cloud properties well, including the cloud locations and cloud particle sizes. We find that the patchy forsterite slab cloud inferred from our retrieval may be responsible for the spectral behavior of the observed variability. Our retrieved cloud structure is consistent with the atmospheric structure previously inferred from spectroscopic variability measurements, but clarifies this picture significantly. We find consistent C/O ratios for both objects, which supports their formation within the same molecular cloud in the Carina-Near moving group. Finally, we note some differences in the constrained abundances of H2O and CO, which may be caused by data quality and/or astrophysical processes such as auroral activity and their differing rotation rates. The results presented in this work provide a promising preview of the detail with which we will characterize extrasolar atmospheres with JWST, which will yield higher-quality spectra across a wider wavelength range.

    more » « less