The interiors of giant planets remain poorly understood. Even for the planets in the Solar System, difficulties in observation lead to major uncertainties in the properties of planetary cores. Exoplanets that have undergone rare evolutionary pathways provide a new route to understanding planetary interiors. We present the discovery of TOI-849b, the remnant core of a giant planet, with a radius smaller than Neptune but an anomalously high mass $M_p=40.8^{+2.4}_{-2.5}M_{\oplus}$ and density of $5.5\pm0.8$gcm$^{-3}$, similar to the Earth. Interior structure models suggest that any gaseous envelope of pure hydrogen and helium consists of no more than $3.9^{+0.8}_{-0.9}$% of the total mass of the planet. TOI-849b transits a late G type star (T$_{\rm mag}=11.5$) with an orbital period of 18.4 hours, leading to an equilibrium temperature of 1800K. The planet's mass is larger than the theoretical threshold mass for runaway gas accretion. As such, the planet could have been a gas giant before undergoing extreme mass loss via thermal self-disruption or giant planet collisions, or it avoided substantial gas accretion, perhaps through gap opening or late formation. Photoevaporation rates cannot provide the mass loss required to reduce a Jupiter-like gas giant, but can remove a few $M_\oplus$ hydrogen and helium envelope on timescales of several Gyr, implying that any remaining atmosphere is likely to be enriched by water or other volatiles from the planetary interior. TOI-849b represents a unique case where material from the primordial core is left over from formation and available to study.

Ultra hot jupiters (UHJs), giant exoplanets with equilibrium temperatures above 2000 K, are ideal laboratories for studying metal compositions of planetary atmospheres. At these temperatures the thermal dissociation of metal-rich molecules into their constituent elements makes these atmospheres conducive for elemental characterisation. Several elements, mostly ionized metals, have been detected in UHJs recently using high resolution transit spectroscopy. Even though a number of neutral transition metals (e.g., Fe, Ti, V, Cr) are expected to be strong sources of optical/NUV opacity and, hence, influence radiative processes in the lower atmospheres of UHJs, only Fe I has been detected to date. We conduct a systematic search for atomic species in the UHJ WASP-121 b. Using theoretical models we present a metric to predict the atomic species likely to be detectable in such planets with high resolution transmission spectroscopy. We search for the predicted species in observations of WASP-121 b and report the first detections of neutral transition metals Cr I and V I in an exoplanet at 3.6 $\sigma$ and 4.5 $\sigma$, respectively. We confirm previous detections of Fe I and Fe II. Whereas Fe II was detected previously in the NUV, we detect it in the optical. We infer that the neutral elements Fe I, V I, and Cr I are present in the lower atmosphere, as predicted by thermochemical equilibrium, while Fe II is a result of photoionisation in the upper atmosphere. Our study highlights the rich chemical diversity of UHJs.

Recent detections of thermal inversions in the dayside atmospheres of some hot Jupiters are motivating new avenues to understand the interplay between their temperature structures and other atmospheric conditions. In particular, TiO has long been proposed to cause thermal inversions in hot Jupiters, depending on other factors such as stellar irradiation, C/O, and vertical mixing. TiO also has spectral features in the optical and near-infrared that have been detected. However, interpretations of TiO signatures rely on the accuracy of TiO opacity used in the models. The recently reported Toto TiO line list provides a new opportunity to investigate these dependencies, which is the goal of the present work. First, we investigate how the Toto line list affects observable transmission and emission spectra of hot Jupiters at low and high resolution. The improvement in the Toto line list compared to a previous line list results in observable differences in the model spectra, particularly in the optical at high resolution. Secondly, we explore the interplay between temperature structure, irradiation and composition with TiO as the primary source of optical opacity, using 1D self-consistent atmospheric models. Among other trends, we find that the propensity for thermal inversions due to TiO peaks at C/O$\sim$0.9, consistent with recent studies. Using these models, we further assess metrics to quantify thermal inversions due to TiO, compared to frequently-used Spitzer photometry, over a range in C/O, irradiation, metallicity, gravity and stellar type.

Atmospheres of a number of ultra-hot Jupiters (UHJs) with temperatures $\gtrsim$2000~K have been observed recently. Many of these planets show largely featureless thermal spectra in the near-infrared observed with the HST WFC3 spectrograph (1.1-1.7~$\mu$m) even though this spectral range contains strong H$_2$O opacity. Recent works have proposed the possibility of H- opacity masking the H$_2$O feature and/or thermal dissociation of H$_2$O causing its apparent depletion at the high temperatures of UHJs. In this work we test these hypotheses using observations of the exoplanet WASP-18b as a case study. We report detailed atmospheric retrievals of the planet using the HyDRA retrieval code, extended to include the effects of H- opacity and thermal dissociation. We report constraints on the H$_2$O, CO and H- abundances as well as the pressure-temperature profile of the dayside atmosphere for retrievals with and without H-/dissociation for each dataset. We find that the H$_2$O and H- abundances are relatively unconstrained given the featureless WFC3 spectra. We do not conclusively detect H- in the planet contrary to previous studies which used equilibrium models to infer its presence. The constraint on the CO abundance depends on the combination of WFC3 and Spitzer data, ranging from solar to super-solar CO values. We additionally see signs of a thermal inversion from two of the datasets. Our study demonstrates the potential of atmospheric retrievals of UHJs including the effects of H- and thermal dissociation of molecules.

High resolution spectroscopy (HRS) has been used to detect a number of species in the atmospheres of hot Jupiters. Key to such detections is accurately and precisely modelled spectra for cross-correlation against the R$\gtrsim$20,000 observations. There is a need for the latest generation of opacities which form the basis for high signal-to-noise detections using such spectra. In this study we present and make publicly available cross sections for six molecular species, H$_2$O, CO, HCN, CH$_4$, NH$_3$ and CO$_2$ using the latest line lists most suitable for low- and high-resolution spectroscopy. We focus on the infrared (0.95-5~$\mu$m) and between 500-1500~K where these species have strong spectral signatures. We generate these cross sections on a grid of pressures and temperatures typical for the photospheres of super Earth, warm Neptunes and hot Jupiters using the latest H$_2$ and He pressure broadening. We highlight the most prominent infrared spectral features by modelling three representative exoplanets, GJ~1214~b, GJ~3470~b and HD~189733~b, which encompass a wide range in temperature, mass and radii. In addition, we verify the line lists for H$_2$O, CO and HCN with previous high resolution observations of hot Jupiters. However, we are unable to detect CH$_4$ with our new cross sections from HRS observations of HD~102195~b. These high accuracy opacities are critical for atmospheric detections with HRS and will be continually updated as new data becomes available.

High resolution spectroscopy (R > 20,000) is currently the only known method to constrain the orbital solution and atmospheric properties of non-transiting hot Jupiters. It does so by resolving the spectral features of the planet into a forest of spectral lines and directly observing its Doppler shift while orbiting the host star. In this study, we analyse VLT/CRIRES (R = 100,000) L-band observations of the non-transiting giant planet HD 179949 b centred around 3.5 microns. We observe a weak (3.0 sigma, or S/N = 4.8) spectral signature of H2O in absorption contained within the radial velocity of the planet at superior-conjunction, with a mild dependence on the choice of line list used for the modelling. Combining this data with previous observations in the K-band, we measure a detection significance of 8.4 sigma for an atmosphere that is most consistent with a shallow lapse-rate, solar C/O ratio, and with CO and H2O being the only major sources of opacity in this wavelength range. As the two sets of data were taken three years apart, this points to the absence of strong radial-velocity anomalies due, e.g., to variability in atmospheric circulation. We measure a projected orbital velocity for the planet of KP = (145.2 +- 2.0)kms^{-1} (1 sigma) and improve the error bars on this parameter by ~70%. However, we only marginally tighten constraints on orbital inclination (66.2 +3.7 -3.1 degrees) and planet mass (0.963 +0.036 -0.031 Jupiter masses), due to the dominant uncertainties of stellar mass and semi-major axis. Follow ups of radial-velocity planets are thus crucial to fully enable their accurate characterisation via high resolution spectroscopy.

The class of ultra-hot Jupiters comprises giant exoplanets undergoing intense irradiation from their host stars. They have proved to be a particularly interesting population for their orbital and atmospheric properties. One such planet, WASP-121~b, is in a highly misaligned orbit close to its Roche limit, and its atmosphere exhibits a thermal inversion. These properties make WASP-121~b an interesting target for additional atmospheric characterization. In this paper, we present analysis of archival high-resolution optical spectra obtained during transits of WASP-121~b. Artifacts from the Rossiter-McLaughlin effect and Center-to-Limb Variation are deemed negligible. However, we discuss scenarios where these effects warrant more careful treatment by modeling the WASP-121 system and varying its properties. We report a new detection of atmospheric absorption from H$\alpha$ in the planet with a transit depth of $1.87\pm0.11\%$. We further confirm a previous detection of the Na I doublet, and report a new detection of Fe I via cross-correlation with a model template. We attribute the H$\alpha$ absorption to an extended Hydrogen atmosphere, potentially undergoing escape, and the Fe I to equilibrium chemistry at the planetary photosphere. These detections help to constrain the composition and chemical processes in the atmosphere of WASP-121~b.

High-resolution Doppler spectroscopy has been used to detect several chemical species in exoplanetary atmospheres. Such detections have traditionally relied on cross correlation of observed spectra against spectral model templates, an approach that is successful for detecting chemical species but not optimised for constraining abundances. Recent work has explored ways to perform atmospheric retrievals on high-resolution spectra (HRS) and combine them with retrievals routinely performed for low-resolution spectra (LRS) by developing a mapping from the cross correlation function to a likelihood metric. We build upon previous studies and report HyDRA-H, a hybrid retrieval code for simultaneous analysis of low- and high- resolution thermal emission spectra of exoplanets in a fully Bayesian approach. We demonstrate HyDRA-H on the hot Jupiter HD 209458b as a case study. We validate our HRS retrieval capability by confirming previous results and report a simultaneous hybrid retrieval using both HRS and LRS data. The LRS data span the HST WFC3 (1.1-1.7 $\mu$m) and Spitzer photometry (IRAC 3.6-8 $\mu$m) bands, while the HRS data were obtained with CRIRES on VLT at 2.3 $\mu$m. The constraints on the composition and temperature profiles for the hybrid retrieval are more stringent than retrievals with either LRS or HRS datasets individually. We retrieve abundances of $\log(\mathrm{H_2O)} = -4.11^{+0.91}_{-0.30}$ and $\log(\mathrm{CO}) = {-2.16}^{+0.99}_{-0.47}$, and $\mathrm{C/O} = 0.99^{+0.01}_{-0.02}$, consistent with previous works. We constrain the photospheric temperature to be $1498^{+216}_{-57}$ K, consistent with the equilibrium temperature. Our results demonstrate the significant advantages of hybrid retrievals by combining strengths of both HRS and LRS observations which probe complementary aspects of exoplanetary atmospheres.

Thermal emission spectra of hot Jupiters have led to key constraints on thermal inversions (or 'stratospheres') in their atmospheres with important implications for their atmospheric processes. Canonically, thermal inversions in hot Jupiters have been suggested to be caused by species such as TiO and VO which have strong visible opacity to absorb incident starlight. We explore two new avenues for thermal inversions in hot Jupiters, exploring both the visible and infrared opacities in their atmospheres. Firstly, by exploring a range of metal-rich species we find that four species (AlO, CaO, NaH and MgH) provide visible opacities comparable to TiO/VO and can cause strong inversions with reasonable abundances. Secondly, we show that a low infrared opacity caused by a low H$_2$O abundance, e.g. through a C/O~1, can also lead to strong thermal inversions even with low abundances of the visible absorbers mentioned above. We find that increasing the C/O ratio towards unity requires almost 2 orders of magnitude lower abundances for the visible absorbers in order for an inversion to form. Finally, we explore the thermal inversion in WASP-121b and find that it can be explained by all the visible absorbers listed above for different C/O ratios. Our study demonstrates the importance of both the refractory and volatile species in governing the physicochemical processes in hot Jupiter atmospheres. Spectroscopic observations in the visible have the potential to detect the newly proposed refractory species that can cause thermal inversions in addition to TiO and VO.

There is a dichotomy in the Milky Way in the $[\alpha/$Fe]-[Fe/H] plane, in which stars fall into high-$\alpha$, and low-$\alpha$ sequences. The high-$\alpha$ sequence comprises mostly old stars, and the low-$\alpha$ sequence comprises primarily young stars. The origin of this dichotomy is uncertain. To better understand how the high- and low-$\alpha$ stars are affiliated, we examine if the high- and low-$\alpha$ sequences have distinct orbits at all ages, or if age sets the orbital properties of stars irrespective of their $\alpha$-enhancement. Orbital actions $J_R$, $J_z$, and $J_\phi$ (or $L_z$) are our labels of stellar dynamics. We use ages for 58,278 LAMOST stars (measured to a precision of 40\%) within $\leq$2kpc of the Sun and we calculate orbital actions from proper motions and parallaxes given by Gaia's DR2. We find that \emph{at all ages}, the high- and low-$\alpha$ sequences are dynamically distinct. This implies separate formation and evolutionary histories for the two sequences; a star's membership in the high- or low-$\alpha$ sequence indicates its dynamical properties at a given time. We use action space to make an efficient selection of halo stars and subsequently report a group of old, low-$\alpha$ stars in the halo, which may be a discrete population from an infall event.

High-resolution doppler spectroscopy provides a powerful means for chemical detections in exoplanetary atmospheres. This approach involves monitoring hundreds of molecular lines in the planetary spectrum doppler shifted by the orbital motion of the planet. The molecules are detected by cross-correlating the observed spectrum of the system with a model planetary spectrum. The method has led to molecular detections of H2O, CO, and TiO in hot Jupiters using large ground-based telescopes. Critical to this method, however, is the accurate removal of the stellar and telluric features from the observed spectrum, also known as detrending. Previous molecular detections have relied on specific choices of detrending methods and parameters. However, the robustness of molecular detections across the different choices has not been investigated in detail. We conduct a systematic investigation of the effect of detrending algorithms, parameters, and optimizations on chemical detections using high-resolution spectroscopy. As a case study, we consider the hot Jupiter HD 189733 b. Using multiple methods, we confirm high-significance detections of H2O (4.8$\sigma$) and CO (4.7$\sigma$). Additionally, we report evidence for HCN at high significance (5.0$\sigma$). On the other hand, our results highlight the need for improved metrics and extended observations for robust confirmations of such detections. In particular, we show that detection significances of $\gtrsim$ 4$\sigma$ can be obtained by optimizing detrending at incorrect locations in the planetary velocity space; such false positives can occur in nearly 30\% of cases. We discuss approaches to help distinguish molecular detections from spurious noise.

Transmission spectroscopy of exoplanets has the potential to provide precise measurements of atmospheric chemical abundances, in particular of hot Jupiters whose large sizes and high temperatures make them conducive to such observations. To date, several transmission spectra of hot Jupiters have revealed low amplitude features of water vapour compared to expectations from cloud-free atmospheres of solar metallicity. The low spectral amplitudes in such atmospheres could either be due to the presence of aerosols that obscure part of the atmosphere or due to inherently low abundances of H$_2$O in the atmospheres. A recent survey of transmission spectra of ten hot Jupiters used empirical metrics to suggest atmospheres with a range of cloud/haze properties but with no evidence for H$_2$O depletion. Here, we conduct a detailed and homogeneous atmospheric retrieval analysis of the entire sample and report the H$_2$O abundances, cloud properties, terminator temperature profiles, and detection significances of the chemical species. Our present study finds that the majority of hot Jupiters have atmospheres consistent with sub-solar H$_2$O abundances at their day-night terminators. The best constrained abundances range from $\mathrm{log(H_2O)}$ of $-5.04^{+0.46}_{-0.30}$ to $-3.16^{+0.66}_{-0.69}$, which compared to expectations from solar-abundance equilibrium chemistry correspond to $0.018^{+0.035}_{-0.009}\times$ solar to $1.40^{+4.97}_{-1.11}\times$ solar. Besides H$_2$O we report statistical constraints on other chemical species and cloud/haze properties, including cloud/haze coverage fractions which range from $0.18^{+0.26}_{-0.12}$ to $0.76^{+0.13}_{-0.15}$. The retrieved H$_2$O abundances suggest sub-solar oxygen and/or super-solar C/O ratios, and can provide important constraints on the formation and migration pathways of hot giant exoplanets.

Molecular species in planetary atmospheres provide key insights into their atmospheric processes and formation conditions. In recent years, high-resolution Doppler spectroscopy in the near-infrared has allowed detections of H$_2$O and CO in the atmospheres of several hot Jupiters. This method involves monitoring the spectral lines of the planetary thermal emission Doppler-shifted due to the radial velocity of the planet over its orbit. However, aside from CO and H$_2$O, which are the primary oxygen- and carbon-bearing species in hot H$_2$-rich atmospheres, little else is known about molecular compositions of hot Jupiters. Several recent studies have suggested the importance and detectability of nitrogen-bearing species in such atmospheres. In this Letter, we confirm potential detections of CO and H$_2$O in the hot Jupiter HD 209458b using high-resolution spectroscopy. We also report a cross-correlation peak with a signal-to-noise ratio of $4.7$ from a search for HCN. The results are obtained using high-resolution phase-resolved spectroscopy with the Very Large telescope CRyogenic high-resolution InfraRed Echelle Spectrograph (VLT CRIRES) and standard analysis methods reported in the literature. A more robust treatment of telluric contamination and other residuals would improve confidence and enable unambiguous molecular detections. The presence of HCN could provide constraints on the C/O ratio of HD~209458b and its potential origins.

Exoplanets with relatively clear atmospheres are prime targets for detailed studies of chemical compositions and abundances in their atmospheres. Alkali metals have long been suggested to exhibit broad wings due to pressure broadening, but most of the alkali detections only show very narrow absorption cores, probably because of the presence of clouds. We report the strong detection of the pressure-broadened spectral profiles of Na, K, and Li absorption in the atmosphere of the super-Neptune WASP-127b, at 4.1$\sigma$, 5.0$\sigma$, and 3.4$\sigma$, respectively. We performed a spectral retrieval modeling on the high-quality optical transmission spectrum newly acquired with the 10.4 m Gran Telescopio Canarias (GTC), in combination with the re-analyzed optical transmission spectrum obtained with the 2.5 m Nordic Optical Telescope (NOT). By assuming a patchy cloudy model, we retrieved the abundances of Na, K, and Li, which are super-solar at 3.7$\sigma$ for K and 5.1$\sigma$ for Li (and only 1.8$\sigma$ for Na). We constrained the presence of haze coverage to be around 52%. We also found a hint of water absorption, but cannot constrain it with the global retrieval owing to larger uncertainties in the probed wavelengths. WASP-127b will be extremely valuable for atmospheric characterization in the era of James Webb Space Telescope.

We find evidence for a strong thermal inversion in the dayside atmosphere of the highly irradiated hot Jupiter WASP-18b (T$_{eq}=2411K$, $M=10.3M_{J}$) based on emission spectroscopy from Hubble Space Telescope secondary eclipse observations and Spitzer eclipse photometry. We demonstrate a lack of water vapor in either absorption or emission at 1.4$\mu$m. However, we infer emission at 4.5$\mu$m and absorption at 1.6$\mu$m that we attribute to CO, as well as a non-detection of all other relevant species (e.g., TiO, VO). The most probable atmospheric retrieval solution indicates a C/O ratio of 1 and a high metallicity (C/H=$283^{+395}_{-138}\times$ solar). The derived composition and T/P profile suggest that WASP-18b is the first example of both a planet with a non-oxide driven thermal inversion and a planet with an atmospheric metallicity inconsistent with that predicted for Jupiter-mass planets at $>2\sigma$. Future observations are necessary to confirm the unusual planetary properties implied by these results.

Thermal emission spectra of exoplanets provide constraints on the chemical compositions, pressure-temperature (P-T) profiles, and energy transport in exoplanetary atmospheres. Accurate inferences of these properties rely on the robustness of the atmospheric retrieval methods employed. While extant retrieval codes have provided significant constraints on molecular abundances and temperature profiles in several exoplanetary atmospheres, the constraints on their deviations from thermal and chemical equilibria have yet to be fully explored. Our present work is a step in this direction. We report HyDRA, a disequilibrium retrieval framework for thermal emission spectra of exoplanetary atmospheres. The retrieval code uses the standard architecture of a parametric atmospheric model coupled with Bayesian statistical inference using the Nested Sampling algorithm. For a given dataset, the retrieved compositions and P-T profiles are used in tandem with the GENESIS self-consistent atmospheric model to constrain layer-by-layer deviations from chemical and radiative-convective equilibrium in the observable atmosphere. We demonstrate HyDRA on the Hot Jupiter WASP-43b with a high-precision emission spectrum. We retrieve an H2O mixing ratio of log(H2O) = -3.54^{+0.82}_{-0.52}, consistent with previous studies. We detect H2O and a combined CO/CO2 at 8-sigma significance. We find the dayside P-T profile to be consistent with radiative-convective equilibrium within the 1-sigma limits and with low day-night redistribution, consistent with previous studies. The derived compositions are also consistent with thermochemical equilibrium for the corresponding distribution of P-T profiles. In the era of high precision and high resolution emission spectroscopy, HyDRA provides a path to retrieve disequilibrium phenomena in exoplanetary atmospheres.

We are entering the era of high-precision and high-resolution spectroscopy of exoplanets. Such observations herald the need for robust self-consistent spectral models of exoplanetary atmospheres to investigate intricate atmospheric processes and to make observable predictions. Spectral models of plane-parallel exoplanetary atmospheres exist, mostly adapted from other astrophysical applications, with different levels of sophistication and accuracy. There is a growing need for a new generation of models custom-built for exoplanets and incorporating state-of-the-art numerical methods and opacities. The present work is a step in this direction. Here we introduce GENESIS, a plane-parallel, self-consistent, line-by-line exoplanetary atmospheric modelling code which includes (a) formal solution of radiative transfer using the Feautrier method, (b) radiative-convective equilibrium with temperature correction based on the Rybicki linearisation scheme, (c) latest absorption cross sections, and (d) internal flux and external irradiation, under the assumptions of hydrostatic equilibrium, local thermodynamic equilibrium and thermochemical equilibrium. We demonstrate the code here with cloud-free models of giant exoplanetary atmospheres over a range of equilibrium temperatures, metallicities, C/O ratios and spanning non-irradiated and irradiated planets, with and without thermal inversions. We provide the community with theoretical emergent spectra and pressure-temperature profiles over this range, along with those for several known hot Jupiters. The code can generate self-consistent spectra at high resolution and has the potential to be integrated into general circulation and non-equilibrium chemistry models as it is optimised for efficiency and convergence. GENESIS paves the way for high-fidelity remote sensing of exoplanetary atmospheres at high resolution with current and upcoming observations.

As an exoplanet transits its host star, some of the light from the star is absorbed by the atoms and molecules in the planet's atmosphere, causing the planet to seem bigger; plotting the planet's observed size as a function of the wavelength of the light produces a transmission spectrum. Measuring the tiny variations in the transmission spectrum, together with atmospheric modelling, then gives clues to the properties of the exoplanet's atmosphere. Chemical species composed of light elements$-$such as hydrogen, oxygen, carbon, sodium and potassium$-$have in this way been detected in the atmospheres of several hot giant exoplanets, but molecules composed of heavier elements have thus far proved elusive. Nonetheless, it has been predicted that metal oxides such as titanium oxide (TiO) and vanadium oxide occur in the observable regions of the very hottest exoplanetary atmospheres, causing thermal inversions on the dayside. Here we report the detection of TiO in the atmosphere of the hot-Jupiter planet WASP-19b. Our combined spectrum, with its wide spectral coverage, reveals the presence of TiO (to a confidence level of 7.7{\sigma}), a strongly scattering haze (7.4{\sigma}) and sodium (3.4{\sigma}), and confirms the presence of water (7.9{\sigma}) in the atmosphere.

The formation mechanisms and chemical compositions of brown dwarfs likely span a diverse range. If they are formed predominantly in isolation by gravitational collapse like dwarf stars then their compositions might follow those of dwarf stars, spanning a relatively narrow range in metallicities (~ -0.4 to +0.4 dex) and predominantly oxygen-rich (C/O <~ 0.7). On the other hand, giant planets in the solar system are all super-solar in metallicity ([C/H] > 0.5~dex), which is thought to be a consequence of formation by core-accretion in a circumstellar disk. In this study we determine the atmospheric chemical compositions of three brown dwarfs and find them to be neither characteristic of dwarf stars nor giant planets. We derive high-precision atmospheric chemical abundances using high-S/N HST near-infrared spectra of three mid T Dwarfs which, together with two previously studied objects, display distinctly sub-solar metallicities and span C/O ratios of 0.4-1.0. These abundance patterns indicate either an old sub-stellar population and a multitude of formation environments or objects formed in circumstellar disks followed by ejection later. We show that for a simple synthetic model population of brown dwarfs with the same spectral-type as our sample the predicted age distribution is predominantly older than the Sun and, hence, more metal poor, consistent with our derived abundances. The diverse C/O ratios we find are indicative of either different chemical reservoirs in their formation environments or different formation pathways. Our results open the possibility that cool brown dwarfs may provide important probes of early galactic chemical evolution and inhomogeneity.