Eileen C. Gonzales

In this study, we present the full (97.3% complete) 0.8–12.5 μ m spectral energy distribution (SED) of an L/T transition object, CWISE J210640.16+250729.0 (CW 2106), using the James Webb Space Telescope (JWST). We provide a full characterization of the host star’s elemental abundances and age. We empirically derive the bolometric luminosity (≈ −4.825 dex) of CW 2106, and obtain estimates of its mass ( M ≈ 50–62 M _Jup ), radius ( R ≈ 0.83–0.87 R _Jup ), effective temperature ( T _eff ≈ 1213 K), and surface gravity ( log g ≈ 5.28 dex). We find the near-infrared (near-IR) spectrum (0.8–2.5 μ m) is best reproduced with cloudy atmospheric models while the mid-infrared (mid-IR) spectrum (5–12.5 μ m) is best reproduced with cloudless models. This suggests a cloud layer restricted to only the deepest observable parts of the atmosphere and is qualified by the lack of a 9 μ m silicate feature. Making use of the Mg/Si ratio of the primary, alongside thermochemical models, we predict the clouds in CW 2106 to be composed primarily of enstatite (MgSiO ₃ ), removing ∼23% of the bulk oxygen out of the atmosphere. Future retrieval studies will be able to help investigate the existence and full impact of these cloud species.

Galactic archaeology relies on stellar kinematics and chemical abundances to identify various stellar populations and associations. With JWST, Galactic archaeology of ancient metal-poor brown dwarfs is now possible through high-sensitivity, medium-resolution spectroscopy. We present the first metal-poor brown dwarf radial velocity (RV) survey using JWST/NIRSpec G395H spectra ( λ /Δ λ ∼ 3000) for a sample of 23 LTY dwarfs across a wide range of temperatures and metallicities. We introduce a forward-modeling framework to measure robust RVs for brown dwarfs by focusing on the CO fundamental band at ∼4.5 μ m. By comparing to sources with existing high-resolution RV measurements, we demonstrate that G395H RVs can achieve a systematic uncertainty of 5 km s ⁻¹ , sufficient to examine their Galactic space motions. We identify kinematic members of the Milky Way’s thin disk, thick disk, and halo populations in this sample. We also confirm WISE J155349.98+693355.2 as likely associated with the Gaia–Enceladus merger remnant, and 2MASS J05325346+8246465 is likely associated with the Thamnos stream. Our study systematically demonstrates that JWST G395H spectroscopy can provide robust RVs in addition to detailed chemical abundances, necessary for Galactic archaeology studies. Future discoveries and characterization of brown dwarfs associated with Galactic halo substructures and streams will enable a stringent test of stellar/substellar formation and evolutionary models for diverse metallicities.

Brown dwarf variability directly probes atmospheric dynamics beyond the Solar System, and recent James Webb Space Telescope (JWST) time-resolved spectroscopy has opened a new window into these processes. Principal component analysis (PCA) offers a data-driven framework to identify the dominant, independent patterns of spectral variability of variable targets without relying on prior atmospheric assumptions. simp (simp) is a young, T2.5, brown dwarf at the planetary-mass boundary, making it an ideal analogue for directly imaged exoplanets. We analysed one rotation of JWST /NIRSpec PRISM time-series spectroscopy to investigate the drivers of its variability using PCA. Prior analyses of these data revealed wavelength-dependent variability consistent with multi-layer atmospheric structure and thermal variations atop patchy clouds, motivating a complementary, model-independent approach. Our PCA showed that the variability is intrinsically low dimensional. Two principal components are sufficient to reduce the residual spectra to the propagated noise floor, indicating that the detectable coherent spectroscopic variability is captured by the first two components. The leading principal component captures broadband variability consistent with temperature changes, while the second traces chromatic variability linked to vertical cloud structure. The existence of two dominant components implies that the spectra can be described as mixtures of three distinct atmospheric states. We identified these three extreme spectral states from PCA projections and mapped their relative contributions as a function of rotational phase. The observed spectra are described as evolving linear combinations of these states, indicating that the variability arises from the changing visibility of spatially distinct atmospheric regions. By projecting Sonora Diamondback forward models into the same principal component space, we found that the eigenspectra capture a large fraction of the model variance, demonstrating that the same physical processes govern much of the model grid and observed variability. These results show that changes in temperature and cloud vertical structure account for most of the variability in simp, and establish PCA as a computationally efficient, physically interpretable framework for analysing JWST time-resolved spectroscopy of substellar atmospheres.

We present results from a spectroscopic survey of 32 FGK stars hosting brown dwarfs, using high-resolution optical spectra ( R = 130,000 and 50,000) obtained with the Potsdam Echelle Polarimetric and Spectroscopic Instrument spectrograph on the Large Binocular Telescope. The primary goal of this survey is to determine precise stellar parameters and abundances for 11 elements (C, O, Mg, Si, Ca, Al, Ti, Fe, Y, S, and N) in these systems. We employ spectral synthesis within the BACCHUS framework to derive precise stellar properties and elemental abundance ratios. For our data with average signal-to-noise ratio > 200, we achieve a typical error of 42 K in T _eff and ∼0.03 dex for [Fe/H]. We observe a significant dispersion from a solar C/O ratio among the sample of brown dwarf host stars that host primarily wide-orbit brown dwarfs. Using established theoretical chemical frameworks, we discuss the implications of the observed Mg/Si and Ca/Al ratios for cloud properties in the brown dwarf companions. Finally, we evaluate the applicability of the [Y/Mg] stellar clock for our sample and discuss the broader implications of our results. This work provides a timely and uniform abundance analysis of host stars, supporting extended wavelength brown dwarf observations in the era of JWST.

Studies suggest Jupiter’s aurorae are supplied with plasma from volcanic outflows on the planet’s innermost moon, Io. Repeating bursts of radio emission thought to trace massively scaled-up analogs of Jupiter’s aurorae have been detected around nearly a dozen isolated substellar worlds, yet the source of the electrons fueling the aurorae remains unknown. Volcanism from tidally heated exosatellites may provide the plasma that fuels the aurorae on these worlds. We assess whether transit observations provide a viable means of detecting exosatellites around aurorally active substellar worlds, thereby enabling future tests of this hypothesis. Specifically, we analyze JWST near- and mid-infrared light curves of SIMP J013656.5+093347.3 (SIMP 0136+0933), a 12.7 M _Jup “super-Jupiter,” known to exhibit auroral emission. We demonstrate the capability to detect exosatellites in the SIMP 0136+0933 system with satellite-to-host mass ratios comparable to those of Jupiter’s Galilean moons, achieving detection success rates of 66% for Io-to-Jupiter mass-ratio satellites and 93% for Ganymede-to-Jupiter mass-ratio satellites. Although the existing light curve is sufficient to demonstrate that this technique is capable of detecting transiting exosatellites, the available archival data are too short in duration to place meaningful constraints on the presence of a transiting satellite in this system. We conclude that JWST light curves spanning ∼1.5 days for 4–12 known aurorally active super-Jupiters would be sufficient to yield evidence for or against this hypothesis. A small target sample may suffice, as short satellite periods boost transit probabilities and aurorally active worlds may be preferentially observed near edge-on inclinations.

We report the detection of photometric variability in the directly imaged super-Jupiter β Pictoris b ( β Pic b). Using JWST NIRCam dual-band coronagraphic imaging, we conducted a 16 hr continuous photometric monitoring campaign in the F210M and F410M filters. We developed and validated a time-series photometry framework that combines point-spread function subtraction, principal component analysis for systematic noise removal, and injection-and-recovery tests to confirm signal fidelity. Both light curves show consistent sinusoidal variability at ∼5 σ and ≫5 σ significance in the F210M and F410M bands, respectively. A joint sinusoidal fit yields a rotation period of P _rot = 9.00 ± 0.13 hr and variability amplitudes of 0.85% ± 0.07% and 0.89% ± 0.04% in F210M and F410M, respectively. The near-identical amplitudes and periods in both bands confirm a common astrophysical origin in a heterogeneous atmosphere. Combining P _rot with the previously measured projected rotational velocity, we constrain the line-of-sight spin axis inclination of β Pic b. The result favors an equator-on viewing geometry, consistent with line-of-sight spin–orbit alignment: the planetary spin axis, orbital plane, debris disk, and stellar equator are all mutually aligned. This stands in sharp contrast to the large obliquities of wide-separation companions that are likely formed via gravitational fragmentation. Together with the system’s young age, this observation provides independent dynamical evidence that β Pic b formed via core accretion. This result constitutes the first detection of rotational modulation in a close-in, high-contrast exoplanet that likely formed via core accretion, demonstrating that time-series coronagraphic imaging with JWST opens a powerful new window onto the rotation, atmospheric dynamics, and spin–orbit architecture of this population.

We report the direct imaging discovery of a third exoplanet in the β Pictoris ( β Pic) system. We detect β Pictoris d ( β Pic d) in noncoronagraphic observations obtained with the Very Large Telescope (VLT) Enhanced Resolution Imager and Spectrograph (ERIS), as well as multi-epoch archival datasets from the JWST Near Infrared Camera (NIRCam) and VLT/SPHERE. Astrometric measurements over an 11 yr baseline demonstrate that it is consistent with a gravitationally bound source with orbital motion. Joint multi-planet orbit fits of all three planets in the system yield a semimajor axis of 26 . 0 − 6.1 + 2.2 au and inclination 89 . 0 − 0.6 + 0.7 deg for planet d. β Pic d has a larger orbital semimajor axis than the other known planets in the system, but is coplanar with the inner two planets, and its orbit is consistent with sculpting the inner edge of the debris disk. β Pic d has a contrast of Δ L ′ = 12.11 ± 0.15 mag, with colors and luminosity that closely match those of 51 Eri b, another exoplanet in the β Pic moving group. Its VLT/ERIS and JWST/NIRCam colors are distinct from those of free-floating planetary-mass objects of a similar age and temperature. Its red F410M − F444W color indicates strong CO ₂ absorption in its atmosphere and suggests significant enhancement in metals compared to free-floating objects. From the ATMO hot-start evolutionary models, we estimate an effective temperature of 60 0 − 60 + 45 K and mass of 2.4 ± 0.6 M _Jup , which also closely matches similar estimates for 51 Eri b. β Pic d is among the lowest-mass exoplanets imaged from the ground. This discovery highlights the deep sensitivity achievable with ground-based imaging in the mid-infrared and the discovery potential of future high-contrast observations with the Extremely Large Telescope.

Studying cold brown dwarfs is key to understanding the diverse characteristics of cold giant exoplanets atmospheres. COCONUTS-2 is a wide binary system composed of a T9 brown dwarf and an M3 star, which presents a unique opportunity to characterize a cold benchmark brown dwarf. As part of a JWST program to study the range of physical and atmospheric properties of the coldest brown dwarfs, we obtained NIRSpec G395H spectra ( R ∼ 2700, 2.87−5.13 μ m) and MIRI F1000W, F1280W, and F1800W photometry for COCONUTS-2 b. In this work, we find a 99% probability of the system belonging to the Corona of Ursa Major moving group (414 ± 23 Myr) using BANYAN Σ and its full kinematics. We also reestimate the astrometry of COCONUTS-2b using the MIRI data. We support this membership with a comparison of the rotation period, metallicity, and C/O ratio of the group with those of the COCONUTS-2 system. We also calculate its bolometric luminosity, which, combined with our age estimation, allows us to derive its mass, effective temperature, surface gravity, and radius with high precision. As a result of our analysis, we support the conclusion that COCONUTS-2 b is a planetary-mass object (7.5 ± 0.4 M _Jup ), which was likely formed via the same mechanism as stars. In addition we compare the JWST spectrum to another object in the sample, J082507.35+280548.5 (0825+2805), a Y0.5 brown dwarf, which is a candidate member of the same moving group, but has a lower mass (3.7 ± 0.2 M _Jup ). We identify absorption feature differences, which could indicate that 0825+2805 has stronger vertical mixing.

Ross 458 C is a widely separated planetary mass companion at a distance of 1100 au from its host binary, Ross 458 AB. It is a member of a class of very low-mass companions at distances of hundreds to thousands of astronomical units from their host stars. We aim to constrain Ross 458 C’s formation history by fitting its near-IR spectrum with models to constrain its composition. If its composition is similar to its host star, we infer that it likely formed through turbulent fragmentation of the same molecular cloud that formed the host. If its composition is enhanced in heavy elements relative to the host, this lends evidence to formation in the disk and subsequent migration to its current separation. Here, we present high-resolution ( R ∼ 2700) emission spectra of Ross 458 C obtained with JWST NIRSpec’s fixed slit in the F070LP, F100LP, and F170LP filters from 0.8 to 3.1 μ m. We fit these spectra using both grids of forward models (Sonora Bobcat, Sonora Elf Owl, and ExoREM) and atmospheric retrievals ( POSEIDON ). We also constrain the composition of Ross 458 AB by fitting an archival SpeX spectrum with PHOENIX forward models. The forward model grids prefer an enhanced atmospheric metallicity for Ross 458 C relative to the host, but our retrievals return a metallicity consistent with the host within 1 σ . Our results offer new insights into the formation history of Ross 458 C, as well as the efficacy of fitting forward model grids versus retrievals to derive atmospheric properties of directly imaged companions.

Studying cold brown dwarfs is key to understanding the diverse characteristics of cold giant exoplanets atmospheres. COCONUTS-2, is a wide binary system composed of a T9 brown dwarf and an M3 star, which presents a unique opportunity to characterize a cold benchmark brown dwarf. As part of a JWST program to study the range of physical and atmospheric properties of the coldest brown dwarfs, we obtained NIRSpec G395H spectra (R~2700, 2.87-5.13 um) and MIRI F1000W, F1280W, and F1800W photometry for COCONUTS-2b. In this work, we find a 99% probability of the system belonging to the Corona of Ursa Major (414+/-23 Myr) using BANYAN Sigma and its full kinematics. We also re-estimate the astrometry of COCONUTS-2b using the MIRI data. We support the membership with a comparison of rotation period, metallicity and C/O ratio of the group with those of the COCONUTS-2 system. We also calculate its bolometric luminosity, which combined with our age estimation, allows us to derive its mass, effective temperature, surface gravity, and radius with high precision. As a result of our analysis, we support the conclusion that COCONUTS-2b is a planetary mass object (7.5+/-0.4 MJup) which was likely formed via the same mechanism as stars. In addition we compare the JWST spectrum to another object in the sample, J082507.35+280548.5 (0825+2805), a Y0.5 brown dwarf, which is a candidate member of the same moving group, but has a lower mass (3.7+/-0.2 MJup). We identify absorption feature differences which could indicate that 0825+2805 has stronger vertical mixing.