Search

Irradiated brown dwarfs offer a unique opportunity to bridge the gap between stellar and planetary atmospheres. We present high-quality $\mathit{HST}$/WFC3/G141 phase-resolved spectra of the white dwarf + brown dwarf binary GD 1400, covering more than one full rotation of the brown dwarf. Accounting for brightness variations caused by ZZ Ceti pulsations, we revealed weak ($\sim$1\%) phase curve amplitude modulations originating from the brown dwarf. Sub-band light curve exploration in various bands showed no significant wavelength dependence on amplitude or phase shift. Extracted day- and night-side spectra indicated chemically similar hemispheres, with slightly higher day-side temperatures, suggesting efficient heat redistribution or dominance of radiative escape over atmospheric circulation. A simple radiative and energy redistribution model reproduced observed temperatures well. Cloud-inclusive models fit the day and night spectra better than cloudless models, indicating global cloud coverage. We also begin qualitatively exploring atmospheric trends across six irradiated brown dwarfs, from the now complete "Dancing with Dwarfs" WD$-$BD sample. The trend we find in the day-side/night-side temperature and irradiation levels is consistent with efficient heat redistribution for irradiation levels less than $\sim$10$^9$ ergs/s/cm$^2$ and decreasing efficiency above that level., Comment: accepted to the Astrophysical Journal, 24 Dec 2024

The coldest Y spectral type brown dwarfs are similar in mass and temperature to cool and warm ($\sim$200 -- 400 K) giant exoplanets. We can therefore use their atmospheres as proxies for planetary atmospheres, testing our understanding of physics and chemistry for these complex, cool worlds. At these cold temperatures, their atmospheres are cold enough for water clouds to form, and chemical timescales increase, increasing the likelihood of disequilibrium chemistry compared to warmer classes of planets. JWST observations are revolutionizing the characterization of these worlds with high signal-to-noise, moderate resolution near- and mid-infrared spectra. The spectra have been used to measure the abundances of prominent species like water, methane, and ammonia; species that trace chemical reactions like carbon monoxide; and even isotopologues of carbon monoxide and ammonia. Here, we present atmospheric retrieval results using both published fixed-slit (GTO program 1230) and new averaged time series observations (GO program 2327) of the coldest known Y dwarf, WISE 0855-0714 (using NIRSpec G395M spectra), which has an effective temperature of $\sim$ 264 K. We present a detection of deuterium in an atmosphere outside of the solar system via a relative measurement of deuterated methane (CH$_{3}$D) and standard methane. From this, we infer the D/H ratio of a substellar object outside the solar system for the first time. We also present a well-constrained part-per-billion abundance of phosphine (PH$_{3}$). We discuss our interpretation of these results and the implications for brown dwarf and giant exoplanet formation and evolution., Comment: 17 pages, 7 figures, accepted to ApJ Letters

Located 10.888 pc from Earth, COCONUTS-2b is a planetary-mass companion to a young (150-800 Myr) M3 star, with a wide orbital separation (6471 au) and a low companion-to-host mass ratio ($0.021\pm0.005$). We have studied the atmospheric properties of COCONUTS-2b using newly acquired 1.0-2.5 $\mu$m spectroscopy from Gemini/Flamingos-2. The spectral type of COCONUTS-2b is refined to T$9.5 \pm 0.5$ based on comparisons with T/Y dwarf spectral templates. We have conducted an extensive forward-modeling analysis, comparing the near-infrared spectrum and mid-infrared broadband photometry with sixteen state-of-the-art atmospheric model grids developed for brown dwarfs and self-luminous exoplanets near the T/Y transition. The PH$_{3}$-free ATMO2020++, ATMO2020++, and Exo-REM models best match the specific observations of COCONUTS-2b, regardless of variations in the input spectrophotometry. This analysis suggests the presence of disequilibrium chemistry, along with a diabatic thermal structure and/or clouds, in the atmosphere of COCONUTS-2b. All models predict fainter $Y$-band fluxes than observed, highlighting uncertainties in the alkali chemistry models and opacities. We determine a bolometric luminosity of $\log{(L_{\rm bol}/L_{\odot})}=-6.18$ dex, with a 0.5 dex-wide range of $[-6.43,-5.93]$ dex that accounts for various assumptions of models. Using thermal evolution models, we derive an effective temperature of $T_{\rm eff}=483^{+44}_{-53}$ K, a surface gravity of $\log{(g)}=4.19^{+0.18}_{-0.13}$ dex, a radius of $R=1.11^{+0.03}_{-0.04}$ R$_{\rm Jup}$, and a mass of $M=8 \pm 2$ M$_{\rm Jup}$. Various atmospheric model grids consistently indicate that COCONUTS-2b's atmosphere has sub- or near-solar metallicity and C/O. These findings provide valuable insights into COCONUTS-2b's formation history and the potential outward migration to its current wide orbit., Comment: Accepted for publication in AJ. Main text: Pages 1-25, Figures 1-11, Tables 1-4; Appendix: Pages 26-43, Figures 12-15. Mostly unchanged from the previous version, except for footnotes 6-15, which were updated based on suggestions from the AJ data editor. The Gemini/F2 spectrum of COCONUTS-2b is accessible via https://doi.org/10.5281/zenodo.13975825

Due to their short orbital periods and relatively high flux ratios, irradiated brown dwarfs in binaries with white dwarfs offer better opportunities to study irradiated atmospheres than hot Jupiters, which have lower planet-to-star flux ratios. WD1032+011 is an eclipsing, tidally locked white dwarf-brown dwarf binary with a 9950 K white dwarf orbited by a 69.7 M$_{Jup}$ brown dwarf in a 0.09 day orbit. We present time-resolved Hubble Space Telescope Wide Field Camera 3 spectrophotometric data of WD1032+011. We isolate the phase-dependent spectra of WD1032+011B, finding a 210 K difference in brightness temperature between the dayside and nightside. The spectral type of the brown dwarf is identified as L1 peculiar, with atmospheric retrievals and comparison to field brown dwarfs showing evidence for a cloud-free atmosphere. The retrieved temperature of the dayside is $1748^{+66}_{-67}$ K, with a nightside temperature of $1555^{+76}_{-62}$ K, showing an irradiation-driven temperature contrast coupled with inefficient heat redistribution from the dayside to the nightside. The brown dwarf radius is inflated, likely due to the constant irradiation from the white dwarf, making it the only known inflated brown dwarf in an eclipsing white dwarf-brown dwarf binary., Comment: 20 pages, 23 figures. Accepted for Publication in MNRAS

Aerosol opacity has emerged as a critical factor controlling transmission and emission spectra. We provide a simple guideline for the effects of aerosol morphology on opacity and residence time in the atmosphere, as it pertains to transit observations, particularly those with flat spectra due to high altitude aerosols. This framework can be used for understanding complex cloud and haze particle properties before getting into detailed microphysical modeling. We consider high altitude aerosols to be composed of large fluffy particles that can have large residence times in the atmosphere and influence the deposition of stellar flux and/or the emergence of thermal emission in a different way than compact droplet particles, as generally modeled to date for extrasolar planetary atmospheres. We demonstrate the important influence of aggregate particle porosity and composition on the extent of the wavelength independent regime. We also consider how such fluffy particles reach such high altitudes and conclude that the most likely scenario is their local production at high altitudes via UV bombardment and subsequent blanketing of the atmosphere, rather than some mechanism of lofting or transport from the lower atmosphere., Comment: 29 pages, 20 figures, accepted for publication in PASP

Direct mid-infrared signatures of silicate clouds in substellar atmospheres were first detected in Spitzer observations of brown dwarfs, although their existence was previously inferred from near-infrared spectra. With JWST's Mid-Infrared Instrument (MIRI) instrument, we can now more deeply probe silicate features from 8 to 10 microns, exploring specific particle composition, size, and structure. Recent characterization efforts have led to the identification in particular of silica (silicon dioxide, SiO$_2$) cloud features in brown dwarfs and giant exoplanets. Previous modeling, motivated by chemical equilibrium, has primarily focused on magnesium silicates (forsterite, enstatite), crystalline quartz, and amorphous silica to match observations. Here, we explore the previously neglected possibility that other crystalline structures of silica, i.e. polymorphs, may be more likely to form at the pressure and temperature conditions of substellar upper atmospheres. We evaluate JWST's diagnostic potential for these polymorphs and find that existing published transmission data are only able to conclusively distinguish tridymite, but future higher signal-to-noise transmission observations, directly imaged planet observations, and brown dwarf observations may be able to disentangle all four of the silica polymorphs. We ultimately propose that accounting for the distinct opacities arising from the possible crystalline structure of cloud materials may act as a powerful, observable diagnostic tracer of atmospheric conditions, where particle crystallinity records the history of the atmospheric regions through which clouds formed and evolved. Finally, we highlight that high fidelity, accurate laboratory measurements of silica polymorphs are critically needed to draw meaningful conclusions about the identities and structures of clouds in substellar atmospheres., Comment: 17 pages, 6 figures, 3 tables. Resubmitted to ApJL after review. Updated to fix typo in Equation 5 and the resulting Figure 1. Thanks to the community for spotting this issue

The sensitivity and spectral coverage of JWST is enabling us to test our assumptions of ultracool dwarf atmospheric chemistry, especially with regards to the abundances of phosphine (PH$_3$) and carbon dioxide (CO$_2$). In this paper, we use NIRSpec PRISM spectra ($\sim$0.8$-$5.5 $\mu$m, $R\sim$100) of four late T and Y dwarfs to show that standard substellar atmosphere models have difficulty replicating the 4.1$-$4.4 $\mu$m wavelength range as they predict an overabundance of phosphine and an underabundance of carbon dioxide. To help quantify this discrepancy, we generate a grid of models using PICASO based on the Elf Owl chemical and temperature profiles where we include the abundances of these two molecules as parameters. The fits to these PICASO models show a consistent preference for orders of magnitude higher CO$_2$ abundances and a reduction in PH$_3$ abundance as compared to the nominal models. This tendency means that the claimed phosphine detection in UNCOVER$-$BD$-$3 could instead be explained by a CO$_2$ abundance in excess of standard atmospheric model predictions; however the signal-to-noise of the spectrum is not high enough to discriminate between these cases. We discuss atmospheric mechanisms that could explain the observed underabundance of PH$_3$ and overabundance of CO$_2$, including a vertical eddy diffusion coefficient ($K_{\mathrm{zz}}$) that varies with altitude, incorrect chemical pathways, or elements condensing out in forms such as NH$_4$H$_2$PO$_4$. However, our favored explanation for the required CO$_2$ enhancement is that the quench approximation does not accurately predict the CO$_2$ abundance, as CO$_2$ remains in chemical equilibrium with CO after CO quenches., Comment: Accepted, 15 pages, 9 figures

Large area observations of extragalactic deep fields with the James Webb Space Telescope (JWST) have provided a wealth of candidate low-mass L- and T-class brown dwarfs. The existence of these sources, which are at derived distances of hundreds of parsecs to several kiloparsecs from the Sun, has strong implications for the low-mass end of the stellar initial mass function, and the link between stars and planets at low metallicities. In this letter, we present a JWST/NIRSpec PRISM spectrum of brown dwarf JADES-GS-BD-9, confirming its photometric selection from observations taken as part of the JWST Advanced Deep Extragalactic Survey (JADES) program. Fits to this spectrum indicate that the brown dwarf has an effective temperature of 800-900K (T5 - T6) at a distance of $1.8 - 2.3$kpc from the Sun, with evidence of the source being at low metallicity ([M/H] $\leq -0.5$). Finally, because of the cadence of JADES NIRCam observations of this source, we additionally uncover a proper motion between the 2022 and 2023 centroids, and we measure a proper motion of $20 \pm 4$ mas yr$^{-1}$ (a transverse velocity of 214 km s$^{-1}$ at 2.25 kpc). At this predicted metallicity, distance, and transverse velocity, it is likely that this source belongs either to the edge of the Milky Way thick disk or the galactic halo. This spectral confirmation demonstrates the efficacy of photometric selection of these important sources across deep extragalactic JWST imaging., Comment: 10 pages, 4 figures, Accepted for publication in The Astrophysical Journal, August 30 2024

We present infrared spectral energy distributions of 23 late-type T and Y dwarfs obtained with the James Webb Space Telescope. The spectral energy distributions consist of NIRSpec PRISM and MIRI LRS spectra covering the $\sim$1--12 $\mu$m wavelength range at $\lambda/ \Delta \lambda \approx 100$ and broadband photometry at 15, 18, and 21 $\mu$m. The spectra exhibit absorption features common to these objects including H$_2$O, CH$_4$, CO, CO$_2$, and NH$_3$. Interestingly, while the spectral morphology changes relatively smoothly with spectral type at $\lambda < 3$ $\mu$m and $\lambda > 8$ $\mu$m, it shows no clear trend in the 5 $\mu$m region where a large fraction of the flux emerges. The broad wavelength coverage of the data enables us to compute the first accurate measurements of the bolometric fluxes of cool brown dwarfs. Combining these bolometric fluxes with parallaxes from Spitzer and HST, we also obtain the first accurate bolometric luminosities of these cool dwarfs. We then used the Sonora Bobcat solar metallicity evolutionary models to estimate the radii of the dwarfs which results in effective temperature estimates ranging from $\sim$1000 to 350 K with a median uncertainty of $\pm$20 K which is nearly an order of magnitude improvement over previous work. We also discuss how various portions of the spectra either do or do not exhibit a clear sequence when ordered by their effective temperatures., Comment: Accepted, 33 pages, 19 figures, 8 tables, Fixed Appendix table (photometry ordering) and Fixed Table 1 (full designation of WISE2018 and WISE0410 photometry)

Survey strategies for upcoming exoplanet direct imaging missions have considered varying assumptions of prior knowledge. Precursor radial velocity surveys could have detected nearby exo-Earths and provided prior orbit and mass constraints. Alternatively, a direct imaging mission performing astrometry could yield constraints on orbit and phase angle of target planets. Understanding the impact of prior mass and orbit information on planetary characterization is crucial for efficiently recognizing habitable exoplanets. To address this question, we use a reflected-light retrieval tool to infer the atmospheric and bulk properties of directly imaged Earth-analogs while considering varying levels of prior information and signal-to-noise ratio (SNR). Because of the strong correlation between the orbit-related parameters and the planetary radius, prior information on the orbital distance and planetary phase yield tight constraints on the planetary radius: from $R_{\rm{p}}=2.95^{+2.69}_{-1.95}~R_{\oplus}$ without prior knowledge, to $R_{\rm{p}}=1.01^{+0.33}_{-0.19}~R_{\oplus}$ with prior determination of the orbit for $\rm{SNR}=20$ in the visible/near-infrared spectral range, thus allowing size determination from reflected light observations. However, additional knowledge of planet mass does not notably enhance radius ($R_{\rm{p}}=0.98^{+0.17}_{-0.14}~R_{\oplus}$) or atmospheric characterization. Also, prior knowledge of the mass alone does not yield a tight radius constraint ($R_{\rm{p}}=1.64^{+1.29}_{-0.80}~R_{\oplus}$) nor improves atmospheric composition inference. By contrast, because of its sensitivity to gas column abundance, detecting a Rayleigh scattering slope or bounding Rayleigh opacity helps to refine gas mixing ratio inferences without requiring prior mass knowledge. Overall, apart from radius determination, increasing the SNR is more beneficial than additional prior observations., Comment: Accepted for publication in ApJL

We present an atmospheric retrieval analysis of the Y0 brown dwarf WISE J035934.06$-$540154.6 using the low-resolution 0.96--12 $\mu$m JWST spectrum presented in \citet{Beiler_2023}. We obtain volume number mixing ratios of the major gas-phase absorbers (H$_2$O, CH$_4$, CO, CO$_2$, PH$_3$, and H$_2$S) that are 3--5$\times$ more precise than previous work that used HST spectra. We also find an order-of-magnitude improvement in the precision of the retrieved thermal profile, a direct result of the broad wavelength coverage of the JWST data. We used the retrieved thermal profile and surface gravity to generate a grid of chemical forward models with varying metallicity, (C/O)$_\textrm{atm}$, and strengths of vertical mixing as encapsulated by the eddy diffusion coefficient $K_\textrm{zz}$. Comparison of the retrieved abundances with this grid of models suggests that the deep atmosphere of WISE 0359$-$54 shows signs of vigorous vertical mixing with $K_\textrm{zz}=10^9$ [cm$^{2}$ s$^{-1}$]. To test the sensitivity of these results to our 5-knot spline thermal profile model, we performed a second retrieval using the \citet{Madhusudhan_2009} thermal profile model. While the results of the two retrievals generally agree well, we do find differences between the retrieved values of mass and volume number mixing ratio of H$_2$S with fractional differences of the median values of $-$0.64 and $-$0.10, respectively. In addition, the 5-knot thermal profile is consistently warmer at pressure between 1 and 70 bar. Nevertheless, our results underscore the power that the broad-wavelength infrared spectra obtainable with the James Webb Space Telescope have to characterize the atmospheres of cool brown dwarfs.

Brown dwarfs in ultra-short period orbits around white dwarfs offer a unique opportunity to study the properties of tidally-locked, fast rotating (1-3 hr), and highly-irradiated atmospheres. Here, we present phase-resolved spectrophotometry of the white dwarf-brown dwarf (WD-BD) binary SDSS 1557, which is the fifth WD-BD binary in our six-object sample. Using the Hubble Space Telescope Wide Field Camera 3 Near-infrared G141 instrument, the 1.1 to 1.7 $\mu$m phase curves show rotational modulations with semi-amplitudes of 10.5$\pm$0.1%. We observe a wavelength dependent amplitude, with longer wavelengths producing larger amplitudes, while no wavelength dependent phase shifts were identified. The phase-resolved extracted BD spectra exhibit steep slopes and are nearly featureless. A simple radiative energy redistribution atmospheric model recreates the hemisphere integrated brightness temperatures at three distinct phases and finds evidence for weak redistribution efficiency. Our model also predicts a higher inclination than previously published. We find that SDSS 1557B, the second most irradiated BD in our sample, is likely dominated by clouds on the night side, whereas the featureless day side spectrum is likely dominated by H$^-$ opacity and a temperature inversion, much like the other highly-irradiated BD EPIC2122B., Comment: 19 pages and 11 figures. Accepted to Astrophysical Journal

The launch of the James Webb Space Telescope (JWST) marks a pivotal moment for precise atmospheric characterization of Y dwarfs, the coldest brown dwarf spectral type. In this study, we leverage moderate spectral resolution observations (R $\sim$ 2700) with the G395H grating of the Near-Infrared Spectrograph (NIRSpec) onboard of JWST to characterize the nearby (9.9 pc) Y dwarf WISEPA J182831.08+265037.8 (WISE 1828). With the NIRSpec G395H 2.88-5.12 $\mathrm{\mu}$m spectrum, we measure the abundances of CO, CO$_2$, CH$_4$, H$_2$S, NH$_3$, and H$_2$O, which are the major carbon, nitrogen, oxygen, and sulfur bearing species in the atmosphere. Based on the retrieved volume mixing ratios with the atmospheric retrieval framework CHIMERA, we report that the C/O ratio is $0.45 \pm 0.01$, close to the solar C/O value of 0.55, and the metallicity to be +0.30 $\pm$ 0.02 dex. Comparison between the retrieval results with the forward modeling results suggests that the model bias for C/O and metallicity could be as high as 0.03 and 0.97 dex respectively. We also report a lower limit of the $^{12}$CO/$^{13}$CO ratio of $>40 $, being consistent with the nominal solar value of 90. Our results highlight the potential of JWST in measuring the C/O ratios down to percent-level precision and characterizing isotopologues of cold planetary atmospheres similar to WISE 1828., Comment: 18 pages + references, including 11 figures, accepted for publication in The Astronomical Journal

We present a new grid of cloudy atmosphere and evolution models for substellar objects. These models include the effect of refractory cloud species, including silicate clouds, on the spectra and evolution. We include effective temperatures from 900 to 2400 K and surface gravities from log g=3.5-5.5, appropriate for a broad range of objects with masses between 1 and 84 Jupiter masses. Model pressure-temperature structures are calculated assuming radiative-convective and chemical equilibrium. We consider the effect of both clouds and metallicity on the atmospheric structure, resulting spectra, and thermal evolution of substellar worlds. We parameterize clouds using the Ackerman & Marley (2001) cloud model, including cloud parameter fsed values from 1-8; we include three metallicities (-0.5, 0.0, and +0.5). Refractory clouds and metallicity both alter the evolution of substellar objects, changing the inferred temperature at a given age by up to 100-200 K. We compare to the observed photometry of brown dwarfs, finding broad agreement with the measured photometry. We publish the spectra, evolution, and other data products online with open access., Comment: 26 pages, 18 figures, submitted for publication in ApJ

Disequilibrium chemistry due to vertical mixing in the atmospheres of many brown dwarfs and giant exoplanets is well-established. Atmosphere models for these objects typically parameterize mixing with the highly uncertain $K_{\rm zz}$ diffusion parameter. The role of mixing in altering the abundances of C-N-O-bearing molecules has mostly been explored for solar composition atmospheres. However, atmospheric metallicity and the C/O ratio also impact atmospheric chemistry. Therefore, we present the \texttt{Sonora Elf Owl} grid of self-consistent cloud-free 1D radiative-convective equilibrium model atmospheres for JWST observations, which includes a variation of $K_{\rm zz}$ across several orders of magnitude and also encompasses sub-solar to super-solar metallicities and C/O ratios. We find that the impact of $K_{\rm zz}$ on the $T(P)$ profile and spectra is a strong function of both $T_{\rm eff}$ and metallicity. For metal-poor objects $K_{\rm zz}$ has large impacts on the atmosphere at significantly higher $T_{\rm eff}$ compared to metal-rich atmospheres where the impact of $K_{\rm zz}$ is seen to occur at lower $T_{\rm eff}$. We identify significant spectral degeneracies between varying $K_{\rm zz}$ and metallicity in multiple wavelength windows, in particular at 3-5 $\mu$m. We use the \texttt{Sonora Elf Owl} atmospheric grid to fit the observed spectra of a sample of 9 early to late T- type objects from $T_{\rm eff}=550-1150$ K. We find evidence for very inefficient vertical mixing in these objects with inferred $K_{\rm zz}$ values lying in the range between $\sim$ 10$^1$-10$^4$ cm$^2$s$^{-1}$. Using self-consistent models, we find that this slow vertical mixing is due to the observations probing mixing in the deep detached radiative zone in these atmospheres., Comment: Accepted for Publication in The Astrophysical Journal, 16 Figures, 3 Tables, 28 Pages

The unprecedented medium-resolution (R~1500-3500) near- and mid-infrared (1-18um) spectrum provided by JWST for the young (140+/-20Myr) low-mass (12-20MJup) L-T transition (L7) companion VHS1256b gives access to a catalogue of molecular absorptions. In this study, we present a comprehensive analysis of this dataset utilizing a forward modelling approach, applying our Bayesian framework, ForMoSA. We explore five distinct atmospheric models to assess their performance in estimating key atmospheric parameters: Teff, log(g), [M/H], C/O, gamma, fsed, and R. Our findings reveal that each parameter's estimate is significantly influenced by factors such as the wavelength range considered and the model chosen for the fit. This is attributed to systematic errors in the models and their challenges in accurately replicating the complex atmospheric structure of VHS1256b, notably the complexity of its clouds and dust distribution. To propagate the impact of these systematic uncertainties on our atmospheric property estimates, we introduce innovative fitting methodologies based on independent fits performed on different spectral windows. We finally derived a Teff consistent with the spectral type of the target, considering its young age, which is confirmed by our estimate of log(g). Despite the exceptional data quality, attaining robust estimates for chemical abundances [M/H] and C/O, often employed as indicators of formation history, remains challenging. Nevertheless, the pioneering case of JWST's data for VHS1256b has paved the way for future acquisitions of substellar spectra that will be systematically analyzed to directly compare the properties of these objects and correct the systematics in the models., Comment: 32 pages, 16 figures, 6 tables, 2 appendices

We present aperture masking interferometry (AMI) observations of the star HIP 65426 at $3.8\,\rm{\mu m}$ as a part of the JWST Direct Imaging Early Release Science (ERS) program obtained using the Near Infrared Imager and Slitless Spectrograph (NIRISS) instrument. This mode provides access to very small inner working angles (even separations slightly below the Michelson limit of $0.5\lambda/D$ for an interferometer), which are inaccessible with the classical inner working angles of the JWST coronagraphs. When combined with JWST's unprecedented infrared sensitivity, this mode has the potential to probe a new portion of parameter space across a wide array of astronomical observations. Using this mode, we are able to achieve a $5\sigma$ contrast of $\Delta m{\sim}7.62{\pm}0.13$ mag relative to the host star at separations ${\gtrsim}0.07{"}$, and the contrast deteriorates steeply at separations ${\lesssim}0.07{"}$. However, we detect no additional companions interior to the known companion HIP 65426 b (at separation ${\sim}0.82{"}$ or, $87^{+108}_{-31}\,\rm{au}$). Our observations thus rule out companions more massive than $10{-}12\,\rm{M_{Jup}}$ at separations ${\sim}10{-}20\,\rm{au}$ from HIP 65426, a region out of reach of ground or space-based coronagraphic imaging. These observations confirm that the AMI mode on JWST is sensitive to planetary mass companions at close-in separations (${\gtrsim}0.07{"}$), even for thousands of more distant stars at $\sim$100 pc, in addition to the stars in the nearby young moving groups as stated in previous works. This result will allow the planning and successful execution of future observations to probe the inner regions of nearby stellar systems, opening an essentially unexplored parameter space., Comment: 23 pages, 14 figures, accepted for publication in ApJL

We present a performance analysis for the aperture masking interferometry (AMI) mode on board the James Webb Space Telescope Near Infrared Imager and Slitless Spectrograph (JWST/NIRISS). Thanks to self-calibrating observables, AMI accesses inner working angles down to and even within the classical diffraction limit. The scientific potential of this mode has recently been demonstrated by the Early Release Science (ERS) 1386 program with a deep search for close-in companions in the HIP 65426 exoplanetary system. As part of ERS 1386, we use the same data set to explore the random, static, and calibration errors of NIRISS AMI observables. We compare the observed noise properties and achievable contrast to theoretical predictions. We explore possible sources of calibration errors and show that differences in charge migration between the observations of HIP 65426 and point-spread function calibration stars can account for the achieved contrast curves. Lastly, we use self-calibration tests to demonstrate that with adequate calibration NIRISS F380M AMI can reach contrast levels of $\sim9-10$ mag at $\gtrsim \lambda/D$. These tests lead us to observation planning recommendations and strongly motivate future studies aimed at producing sophisticated calibration strategies taking these systematic effects into account. This will unlock the unprecedented capabilities of JWST/NIRISS AMI, with sensitivity to significantly colder, lower-mass exoplanets than lower-contrast ground-based AMI setups, at orbital separations inaccessible to JWST coronagraphy., Comment: 20 pages, 12 figures, accepted to Astrophysical Journal Letters

We present model constraints on the atmospheric structure of HD 106906 b, a planetary-mass companion orbiting at a ~700 AU projected separation around a 15 Myr-old stellar binary, using the APOLLO retrieval code on spectral data spanning 1.1-2.5 $\mu$m. C/O ratios can provide evidence for companion formation pathways, as such pathways are ambiguous both at wide separations and at star-to-companion mass ratios in the overlap between the distributions of planets and brown dwarfs. We benchmark our code against an existing retrieval of the field L dwarf 2M2224-0158, returning a C/O ratio consistent with previous fits to the same JHKs data, but disagreeing in the thermal structure, cloud properties, and atmospheric scale height. For HD 106906 b, we retrieve C/O $=0.53^{+0.15}_{-0.25}$, consistent with the C/O ratios expected for HD 106906's stellar association and therefore consistent with a stellar-like formation for the companion. We find abundances of H$_2$O and CO near chemical equilibrium values for a solar metallicity, but a surface gravity lower than expected, as well as a thermal profile with sharp transitions in the temperature gradient. Despite high signal-to-noise and spectral resolution, more accurate constraints necessitate data across a broader wavelength range. This work serves as preparation for subsequent retrievals in the era of JWST, as JWST's spectral range provides a promising opportunity to resolve difficulties in fitting low-gravity L dwarfs, and also underscores the need for simultaneous comparative retrievals on L dwarf companions with multiple retrieval codes., Comment: 39 pages, 21 figures. Accepted for publication in the Astronomical Journal

AF Lep A+b is a remarkable planetary system hosting a gas-giant planet that has the lowest dynamical mass among directly imaged exoplanets. We present an in-depth analysis of the atmospheric composition of the star and planet to probe the planet's formation pathway. Based on new high-resolution spectroscopy of AF Lep A, we measure a uniform set of stellar parameters and elemental abundances (e.g., [Fe/H] = $-0.27 \pm 0.31$ dex). The planet's dynamical mass ($2.8^{+0.6}_{-0.5}$ M$_{\rm Jup}$) and orbit are also refined using published radial velocities, relative astrometry, and absolute astrometry. We use petitRADTRANS to perform chemically-consistent atmospheric retrievals for AF Lep b. The radiative-convective equilibrium temperature profiles are incorporated as parameterized priors on the planet's thermal structure, leading to a robust characterization for cloudy self-luminous atmospheres. This novel approach is enabled by constraining the temperature-pressure profiles via the temperature gradient $(d\ln{T}/d\ln{P})$, a departure from previous studies that solely modeled the temperature. Through multiple retrievals performed on different portions of the $0.9-4.2$ $\mu$m spectrophotometry, along with different priors on the planet's mass and radius, we infer that AF Lep b likely possesses a metal-enriched atmosphere ([Fe/H] $> 1.0$ dex). AF Lep b's potential metal enrichment may be due to planetesimal accretion, giant impacts, and/or core erosion. The first process coincides with the debris disk in the system, which could be dynamically excited by AF Lep b and lead to planetesimal bombardment. Our analysis also determines $T_{\rm eff} \approx 800$ K, $\log{(g)} \approx 3.7$ dex, and the presence of silicate clouds and dis-equilibrium chemistry in the atmosphere. Straddling the L/T transition, AF Lep b is thus far the coldest exoplanet with suggested evidence of silicate clouds., Comment: AJ, in press. Main text: Pages 1-32, Figures 1-15, Tables 1-6. All figures and tables after References belong to the Appendix (Pages 32-58, Figures 16-20, Table 7). For supplementary materials, please refer to the Zenodo repository https://doi.org/10.5281/zenodo.8267466

The presence of aerosols is intimately linked to the global energy budget and the composition of a planet's atmospheres. Their ability to reflect incoming light prevents energy from being deposited into the atmosphere, and they shape spectra of exoplanets. We observed five near-infrared secondary eclipses of WASP-80b with the Wide Field Camera 3 (WFC3) aboard the \textit{Hubble Space Telescope} to provide constraints on the presence and properties of atmospheric aerosols. We detect a broadband eclipse depth of $34\pm10$\,ppm for WASP-80b. We detect a higher planetary flux than expected from thermal emission alone at $1.6\sigma$, which hints toward the presence of reflecting aerosols on this planet's dayside, indicating a geometric albedo of $A_g<0.33$ at 3$\sigma$. We paired the WFC3 data with Spitzer data and explored multiple atmospheric models with and without aerosols to interpret this spectrum. Albeit consistent with a clear dayside atmosphere, we found a slight preference for near-solar metallicities and for dayside clouds over hazes. We exclude soot haze formation rates higher than $10^{-10.7}$ g cm$^{-2}$s$^{-1}$ and tholin formation rates higher than $10^{-12.0}$ g cm$^{-2}$s$^{-1}$ at $3\sigma$. We applied the same atmospheric models to a previously published WFC3/Spitzer transmission spectrum for this planet and found weak haze formation. A single soot haze formation rate best fits both the dayside and the transmission spectra simultaneously. However, we emphasize that no models provide satisfactory fits in terms of the chi-square of both spectra simultaneously, indicating longitudinal dissimilarity in the atmosphere's aerosol composition., Comment: Published in ApJ Letters (20 Oct 2023)

Understanding the driving forces behind spiral arms in protoplanetary disks remains a challenge due to the faintness of young giant planets. MWC 758 hosts such a protoplanetary disk with a two-armed spiral pattern that is suggested to be driven by an external giant planet. We present new thermal infrared observations that are uniquely sensitive to redder (i.e., colder or more attenuated) planets than past observations at shorter wavelengths. We detect a giant protoplanet, MWC 758c, at a projected separation of ~100 au from the star. The spectrum of MWC 758c is distinct from the rest of the disk and consistent with emission from a planetary atmosphere with Teff = 500 +/- 100 K for a low level of extinction (AV<30), or a hotter object with a higher level of extinction. Both scenarios are commensurate with the predicted properties of the companion responsible for driving the spiral arms. MWC 758c provides evidence that spiral arms in protoplanetary disks can be caused by cold giant planets or by those whose optical emission is highly attenuated. MWC 758c stands out both as one of the youngest giant planets known, and also as one of the coldest and/or most attenuated. Furthermore, MWC 758c is among the first planets to be observed within a system hosting a protoplanetary disk., Comment: Published in Nature Astronomy

We present the first JWST spectral energy distribution of a Y dwarf. This spectral energy distribution of the Y0 dwarf WISE J035934.06$-$540154.6 consists of low-resolution ($\lambda$/$\Delta\lambda$ $\sim$ 100) spectroscopy from 1$-$12 $\mu$m and three photometric points at 15, 18, and 21 $\mu$m. The spectrum exhibits numerous fundamental, overtone, and combination rotational-vibrational bands of H$_2$O, CH$_4$, CO, CO$_2$, and NH$_3$, including the previously unidentified $\nu_3$ band of NH$_3$ at 3 $\mu$m. Using a Rayleigh-Jeans tail to account for the flux emerging at wavelengths greater than 21 $\mu$m, we measure a bolometric luminosity of $1.523\pm0.090\times10^{20}$ W. We determine a semi-empirical effective temperature estimate of $467^{+16}_{-18}$ K using the bolometric luminosity and evolutionary models to estimate a radius. Finally, we compare the spectrum and photometry to a grid of atmospheric models and find reasonably good agreement with a model having $T_{\mathrm{eff}}$=450 K, log $g$=3.25 [cm s$^{-2}$], [M/H]=$-0.3$. However, the low surface gravity implies an extremely low mass of 1 $M_{\rm{Jup}}$ and a very young age of 20 Myr, the latter of which is inconsistent with simulations of volume-limited samples of cool brown dwarfs.

Understanding differences between sub-stellar spectral data and models has proven to be a major challenge, especially for self-consistent model grids that are necessary for a thorough investigation of brown dwarf atmospheres. Using the supervised machine learning method of the random forest, we study the information content of 14 previously published model grids of brown dwarfs (from 1997 to 2021). The random forest method allows us to analyze the predictive power of these model grids, as well as interpret data within the framework of Approximate Bayesian Computation (ABC). Our curated dataset includes 3 benchmark brown dwarfs (Gl 570D, {\epsilon} Indi Ba and Bb) as well as a sample of 19 L and T dwarfs; this sample was previously analyzed in Lueber et al. (2022) using traditional Bayesian methods (nested sampling). We find that the effective temperature of a brown dwarf can be robustly predicted independent of the model grid chosen for the interpretation. However, inference of the surface gravity is model-dependent. Specifically, the BT-Settl, Sonora Bobcat and Sonora Cholla model grids tend to predict logg ~3-4 (cgs units) even after data blueward of 1.2 {\mu}m have been disregarded to mitigate for our incomplete knowledge of the shapes of alkali lines. Two major, longstanding challenges associated with understanding the influence of clouds in brown dwarf atmospheres remain: our inability to model them from first principles and also to robustly validate these models., Comment: Accepted for publication in The Astrophysical Journal

Approximate methods to estimate solutions to the radiative transfer equation are essential for the understanding of atmospheres of exoplanets and brown dwarfs. The simplest and most popular choice is the "two-stream method" which is often used to produce simple yet effective models for radiative transfer in scattering and absorbing media. Toon et al. (1989) (Toon89) outlined a two-stream method for computing reflected light and thermal spectra and was later implemented in the open-source radiative transfer model PICASO. In Part~I of this series, we developed an analytical spherical harmonics method for solving the radiative transfer equation for reflected solar radiation (Rooney et al. 2023), which was implemented in PICASO to increase the accuracy of the code by offering a higher-order approximation. This work is an extension of this spherical harmonics derivation to study thermal emission spectroscopy. We highlight the model differences in the approach for thermal emission and benchmark the 4-term method (SH4) against Toon89 and a high-stream discrete-ordinates method, CDISORT. By comparing the spectra produced by each model we demonstrate that the SH4 method provides a significant increase in accuracy, compared to Toon89, which can be attributed to the increased order of approximation and to the choice of phase function. We also explore the trade-off between computational time and model accuracy. We find that our 4-term method is twice as slow as our 2-term method, but is up to five times more accurate, when compared with CDISORT. Therefore, SH4 provides excellent improvement in model accuracy with minimal sacrifice in numerical expense., Comment: Submitted ApJ; 17 pages; 7 figures; Code available at https://github.com/natashabatalha/picaso; Zenodo release at https://zenodo.org/record/7765171#.ZC3G7uzMI8Y; Tutorials/figure reproducibility at https://natashabatalha.github.io/picaso/notebooks/10c_AnalyzingApproximationsThermal.html

A significant challenge in radiative transfer theory for atmospheres of exoplanets and brown dwarfs is the derivation of computationally efficient methods that have adequate fidelity to more precise, numerically demanding solutions. In this work, we extend the capability of the first open-source radiative transfer model for computing the reflected light of exoplanets at any phase geometry, PICASO: Planetary Intensity Code for Atmospheric Spectroscopy Observations. Until now, PICASO has implemented two-stream approaches to the solving the radiative transfer equation for reflected light, in particular following the derivations of Toon et al. (1989) (Toon89). In order to improve the model accuracy, we have considered higher-order approximations of the phase functions, namely, we have increased the order of approximation from 2 to 4, using spherical harmonics. The spherical harmonics approximation decouples spatial and directional dependencies by expanding the intensity and phase function into a series of spherical harmonics, or Legendre polynomials, allowing for analytical solutions for low-order approximations to optimize computational efficiency. We rigorously derive the spherical harmonics method for reflected light and benchmark the 4-term method (SH4) against Toon89 and two independent and higher-fidelity methods (CDISORT & doubling-method). On average, the SH4 method provides an order of magnitude increase in accuracy, compared to Toon89. Lastly, we implement SH4 within PICASO and observe only modest increase in computational time, compared to two-stream methods (20% increase)., Comment: Accepted ApJ; 27 pages; 5 figures; Code available at https://github.com/natashabatalha/picaso; Zenodo release at https://zenodo.org/record/7765171#.ZC3G7uzMI8Y; Tutorials/figure reproducibility at https://natashabatalha.github.io/picaso/notebooks/10b_AnalyzingApproximationsReflectedLightSH.html

We reanalyze near-infrared spectra of the young extrasolar giant planet 51 Eridani b which was originally presented in (Macintosh et al. 2015) and (Rajan et al. 2017) using modern atmospheric models which include a self-consistent treatment of disequilibrium chemistry due to turbulent vertical mixing. In addition, we investigate the possibility that significant opacity from micrometeors or other impactors in the planet's atmosphere may be responsible for shaping the observed spectral energy distribution (SED). We find that disequilibrium chemistry is useful for describing the mid-infrared colors of the planet's spectra, especially in regards to photometric data at M band around 4.5 $\mu$m which is the result of super-equilibrium abundances of carbon monoxide, while the micrometeors are unlikely to play a pivotal role in shaping the SED. The best-fitting, micrometeroid-dust-free, disequilibrium chemistry, patchy cloud model has the following parameters: effective temperature $T_\textrm{eff} = 681$ K with clouds (or without clouds, i.e. the grid temperature $T_\textrm{grid}$ = 900 K), surface gravity $g$ = 1000 m/s$^2$, sedimentation efficiency $f_\textrm{sed}$ = 10, vertical eddy diffusion coefficient $K_\textrm{zz}$ = 10$^3$ cm$^2$/s, cloud hole fraction $f_\textrm{hole}$ = 0.2, and planet radius $R_\textrm{planet}$ = 1.0 R$_\textrm{Jup}$., Comment: 22 pages, 14 figures, Accepted to AJ

At the lowest masses, the distinction between brown dwarfs and giant exoplanets is often blurred and literature classifications rarely reflect the deuterium burning boundary. Atmospheric characterisation may reveal the extent to which planetary formation pathways contribute to the population of very-low mass brown dwarfs, by revealing if their abundance distributions differ from those of the local field population or, in the case of companions, their primary stars. The T8 dwarf Ross 458c is a possible planetary mass companion to a pair of M dwarfs, and previous work suggests that it is cloudy. We here present the results of the retrieval analysis of Ross 458c, using archival spectroscopic data in the 1.0 to 2.4 micron range. We test a cloud free model as well as a variety of cloudy models and find that the atmosphere of Ross 458c is best described by a cloudy model (strongly preferred). The CH4/H2O is higher than expected at 1.97 +0.13 -0.14. This value is challenging to understand in terms of equilibrium chemistry and plausible C/O ratios. Comparisons to thermochemical grid models suggest a C/O of ~ 1.35, if CH4 and H2O are quenched at 2000 K, requiring vigorous mixing. We find a [C/H] ratio of +0.18, which matches the metallicity of the primary system, suggesting that oxygen is missing from the atmosphere. Even with extreme mixing, the implied C/O is well beyond the typical stellar regime, suggesting a either non-stellar formation pathway, or the sequestration of substantial quantities of oxygen via hitherto unmodeled chemistry or condensation processes., Comment: 16 pages, 9 figures

With infrared flux contrasts larger than typically seen in hot Jupiter, tidally-locked white dwarf-brown dwarf binaries offer a superior opportunity to investigate atmospheric processes in irradiated atmospheres. NLTT5306 is such a system, with a M$_{BD}$ = 52 $\pm$ 3 M$_{Jup}$ brown dwarf orbiting a T$_{eff}$ = 7756 $\pm$ 35 K white dwarf with an ultra-short period of $\sim$102 min. We present HST/WFC3 spectroscopic phase curves of NLTT5306, consisting of 47 spectra from 1.1 to 1.7 microns with an average S/N=65 per wavelength. We extracted the phase-resolved spectra of the brown dwarf NLTT5306B, finding a small (<100 K) day/night temperature difference ($\sim$5 percent of the average day-side temperature). Our best-fit model phase curves revealed a complex wavelength-dependence on amplitudes and relative phase offsets, suggesting longitudinal-vertical atmospheric structure. The night-side spectrum was well-fit by a cloudy, non-irradiated atmospheric model while the day-side was best-matched by a cloudy, weakly irradiated model. Additionally, we created a simple radiative energy redistribution model of the atmosphere and found evidence for efficient day-to-night heat redistribution and a moderately high Bond albedo. We also discovered an internal heat flux much higher than expected given the published system age, leading to an age reassessment that resulted in NLTT5306B most likely being much younger. We find that NLTT5306B is the only known significantly irradiated brown dwarf where the global temperature structure is not dominated by external irradiation, but rather its own internal heat. Our study provides an essential insight into the drivers of global circulation and day-to-night heat transport as a function of irradiation, rotation rate, and internal heat., Comment: 31 pages, 25 figures, Accepted to Astrophysical Journal with video abstract. Video abstract available at https://www.rachaelamaro.com/research-1

We use observations with the infrared-optimized MagAO system and Clio camera in 3.9 $\mu$m light to place stringent mass constraints on possible undetected companions to Sirius A. We suppress the light from Sirius A by imaging it through a grating vector-apodizing phase plate coronagraph with 180-degree dark region (gvAPP-180). To remove residual starlight in post-processing, we apply a time-domain principal-components-analysis-based algorithm we call PCA-Temporal (PCAT), which uses eigen-time-series rather than eigen-images to subtract starlight. By casting the problem in terms of eigen-time-series, we reduce the computational cost of post-processing the data, enabling the use of the fully sampled dataset for improved contrast at small separations. We also discuss the impact of retaining fine temporal sampling of the data on final contrast limits. We achieve post-processed contrast limits of $1.5 \times 10^{-6}$ to $9.8 \times 10^{-6}$ outside of 0.75 arcsec which correspond to planet masses of 2.6 to 8.0 $M_J$. These are combined with values from the recent literature of high-contrast imaging observations of Sirius to synthesize an overall completeness fraction as a function of mass and separation. After synthesizing these recent studies and our results, the final completeness analysis rules out 99% of $\ge 9 \ M_J$ planets from 2.5-7 AU., Comment: 19 pages, 22 figures, accepted to AJ

The HR 2562 system is a rare case where a brown dwarf companion resides in a cleared inner hole of a debris disk, offering invaluable opportunities to study the dynamical interaction between a substellar companion and a dusty disk. We present the first ALMA observation of the system as well as the continued GPI monitoring of the companion's orbit with 6 new epochs from 2016 to 2018. We update the orbital fit and, in combination with absolute astrometry from GAIA, place a 3$\sigma$ upper limit of 18.5 $M_J$ on the companion's mass. To interpret the ALMA observations, we used radiative transfer modeling to determine the disk properties. We find that the disk is well resolved and nearly edge on. While the misalignment angle between the disk and the orbit is weakly constrained due to the short orbital arc available, the data strongly support a (near) coplanar geometry for the system. Furthermore, we find that the models that describe the ALMA data best have an inner radius that is close to the companion's semi-major axis. Including a posteriori knowledge of the system's SED further narrows the constraints on the disk's inner radius and place it at a location that is in reasonable agreement with, possibly interior to, predictions from existing dynamical models of disk truncation by an interior substellar companion. HR\,2562 has the potential over the next few years to become a new testbed for dynamical interaction between a debris disk and a substellar companion., Comment: 22 pages, 17 figures

JWST has begun its scientific mission, which includes the atmospheric characterization of transiting exoplanets. Some of the first exoplanets to be observed by JWST have equilibrium temperatures below 1000 K, which is a regime where photochemical hazes are expected to form. The optical properties of these hazes, which controls how they interact with light, are critical for interpreting exoplanet observations, but relevant data are not available. Here we measure the optical properties of organic haze analogues generated in water-rich exoplanet atmosphere experiments. We report optical constants (0.4 to 28.6 micron) of organic hazes for current and future observational and modeling efforts covering the entire wavelength range of JWST instrumentation and a large part of Hubble. We use these optical constants to generate hazy model atmospheric spectra. The synthetic spectra show that differences in haze optical constants have a detectable effect on the spectra, impacting our interpretation of exoplanet observations. This study emphasizes the need to investigate the optical properties of hazes formed in different exoplanet atmospheres, and establishes a practical procedure to determine such properties., Comment: 4 figures, 1 Table, Paper under review in Nature Astronomy

We present an atmospheric retrieval analysis of a pair of highly variable, $\sim200~$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~\mu$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 (Mg$_2$SiO$_4$) 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 H$_2$O 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., Comment: Accepted for publication in ApJ

We present a $L$-band (2.98--3.96$\mu$m) spectroscopic study of 8 young L dwarfs with spectral types ranging from L2 to L7. Our spectra (${\lambda}/{\Delta \lambda}\approx$ 250 to 600) were collected using the Gemini Near-InfraRed Spectrograph. We first examine the young $L$-band spectral sequence, most notably analyzing the evolution of the $Q$-branch of methane absorption feature at 3.3 $\mu$m. We find the $Q$-branch feature first appears between L3 and L6, as previously seen in older field dwarfs. Secondly, we analyze how well various atmospheric models reproduce the $L$-band and published near-IR (0.7--2.5 $\mu$m) spectra of our objects by fitting five different grids of model spectra to the data. Best-fit parameters for the combined near-IR and $L$-band data are compared to best-fit parameters for just the near-IR data, isolating the impact that the addition of the $L$-band has on the results. This addition notably causes a $\sim$100 K drop in the best-fit effective temperature. Also, when clouds and a vertical mixing rate ($K_{\mathrm{zz}}$) are included in the models, thick clouds and higher $K_{\mathrm{zz}}$ values are preferred. Five of our objects also have previously published effective temperatures and surface gravities derived using evolutionary models, age estimates, and bolometric luminosities. Comparing model spectra matching these parameters to our spectra, we find disequilibrium chemistry and clouds are needed to match these published effective temperatures. Three of these objects are members of AB Dor, allowing us to show the temperature dependence of the $Q$-branch of methane., Comment: 26 pages, 19 figures, Accepted for publication in MNRAS

Comparisons of atmospheric retrievals can reveal powerful insights on the strengths and limitations of our data and modeling tools. In this paper, we examine a sample of 5 similar effective temperature (Teff) or spectral type L dwarfs to compare their pressure-temperature (P-T) profiles. Additionally, we explore the impact of an object's metallicity and the observations' signal-to-noise (SNR) on the parameters we can retrieve. We present the first atmospheric retrievals: 2MASS J15261405$+$2043414, 2MASS J05395200$-$0059019, 2MASS J15394189$-$0520428, and GD 165B increasing the small but growing number of L-dwarfs retrieved. When compared to atmospheric retrievals of SDSS J141624.08+134826.7, a low-metallicity d/sdL7 primary in a wide L+T binary, we find similar Teff sources have similar P-T profiles with metallicity differences impacting the relative offset between their P-T profiles in the photosphere. We also find that for near-infrared spectra, when the SNR is $\gtrsim80$ we are in a regime where model uncertainties dominate over data measurement uncertainties. As such, SNR does not play a role in the retrieval's ability to distinguish between a cloud-free and cloudless model, but may impact the confidence of the retrieved parameters. Lastly, we also discuss how to break cloud model degeneracies and the impact of extraneous gases in a retrieval model., Comment: 19 pages, 3 figures, 7 tables

We present the highest fidelity spectrum to date of a planetary-mass object. VHS 1256 b is a $<$20 M$_\mathrm{Jup}$ widely separated ($\sim$8\arcsec, a = 150 au), young, planetary-mass companion that shares photometric colors and spectroscopic features with the directly imaged exoplanets HR 8799 c, d, and e. As an L-to-T transition object, VHS 1256 b exists along the region of the color-magnitude diagram where substellar atmospheres transition from cloudy to clear. We observed VHS 1256~b with \textit{JWST}'s NIRSpec IFU and MIRI MRS modes for coverage from 1 $\mu$m to 20 $\mu$m at resolutions of $\sim$1,000 - 3,700. Water, methane, carbon monoxide, carbon dioxide, sodium, and potassium are observed in several portions of the \textit{JWST} spectrum based on comparisons from template brown dwarf spectra, molecular opacities, and atmospheric models. The spectral shape of VHS 1256 b is influenced by disequilibrium chemistry and clouds. We directly detect silicate clouds, the first such detection reported for a planetary-mass companion., Comment: Accepted ApJL. Iterations of spectra reduced by the ERS team are hosted at this link: https://github.com/bemiles/JWST_VHS1256b_Reduction/tree/main/reduced_spectra

We present JWST Early Release Science (ERS) coronagraphic observations of the super-Jupiter exoplanet, HIP 65426 b, with the Near-Infrared Camera (NIRCam) from 2-5 $\mu$m, and with the Mid-Infrared Instrument (MIRI) from 11-16 $\mu$m. At a separation of $\sim$0.82" (86$^{+116}_{-31}$ au), HIP 65426 b is clearly detected in all seven of our observational filters, representing the first images of an exoplanet to be obtained by JWST, and the first ever direct detection of an exoplanet beyond 5 $\mu$m. These observations demonstrate that JWST is exceeding its nominal predicted performance by up to a factor of 10, depending on separation and subtraction method, with measured 5$\sigma$ contrast limits of $\sim$1$\times10^{-5}$ and $\sim$2$\times10^{-4}$ at 1" for NIRCam at 4.4 $\mu$m and MIRI at 11.3 $\mu$m, respectively. These contrast limits provide sensitivity to sub-Jupiter companions with masses as low as 0.3$M_\mathrm{Jup}$ beyond separations of $\sim$100 au. Together with existing ground-based near-infrared data, the JWST photometry are well fit by a BT-SETTL atmospheric model from 1-16 $\mu$m, and span $\sim$97% of HIP 65426 b's luminous range. Independent of the choice of model atmosphere we measure an empirical bolometric luminosity that is tightly constrained between $\mathrm{log}\!\left(L_\mathrm{bol}/L_{\odot}\right)$=-4.31 to $-$4.14, which in turn provides a robust mass constraint of 7.1$\pm$1.2 $M_\mathrm{Jup}$. In totality, these observations confirm that JWST presents a powerful and exciting opportunity to characterise the population of exoplanets amenable to high-contrast imaging in greater detail., Comment: 35 pages, 16 figures, 4 tables, 1 wonderful telescope; Submitted to AAS Journals

Evidence of disequilibrium chemistry due to vertical mixing in the atmospheres of many T and Y-dwarfs has been inferred due to enhanced mixing ratios of CO and reduced NH$_3$. Atmospheric models of planets and brown dwarfs typically parameterize this vertical mixing phenomenon with the vertical eddy diffusion coefficient, $K_{\rm zz}$. While $K_{\rm zz}$ can perhaps be approximated in the convective regions in the atmosphere with mixing length theory, in radiative regions the strength of vertical mixing is uncertain by many orders of magnitude. With a new grid of self-consistent 1D model atmospheres from $T_{\rm eff}$ of 400 - 1000 K, computed with a new radiative-convective equilibrium python code PICASO 3.0, we aim to assess how molecular abundances and corresponding spectra can be used as a probe of depth-dependent $K_{\rm zz}$. At a given surface gravity, we find non-monotonic behavior in the CO abundance as a function of $T_{\rm eff}$, as chemical abundances are sometimes quenched in either of two potential atmospheric convective zones, or quenched in either of two possible radiative zones. The temperature structure and chemical quenching behavior also changes with gravity. We compare our models with available near-infrared and M-band spectroscopy of several T and Y-dwarfs and assess their atmospheric vertical mixing profiles. We also compare to color-magnitude diagrams and make predictions for JWST spectra. This work yields new constraints, and points the way to significant future gains, in determining $K_{\rm zz}$, a fundamental atmospheric parameter in substellar atmospheres, with significant implications for chemistry and cloud modeling., Comment: Accepted for Publication in The Astrophysical Journal, 31 Pages, 22 Figures. The model grid will be released via Zenodo

Direct imaging characterization of extrasolar planets is often done at low spectral resolution. We model the spectrograph for the Gemini Planet Imager upgrade (GPI 2.0) and assess the instrument's potential for allowing observers to constrain exoplanet properties through analysis of near-infrared spectra. We simulated noisy observations followed by calculations of posterior distributions from maximum likelihood comparison with the Sonora 2018 model grid. Preliminary results suggest that GPI 2.0 should allow observers to constrain temperature with sufficient accuracy, but gravity remains largely uncertain. We also explore the effects of incorporating convolution with the instrument line spread function into our simulation and compare the results with our preliminary findings., Comment: 11 pages, 10 figures, to be published in proceedings for SPIE Astronomical Telescopes + Instrumentation 2022

Carbon dioxide (CO2) is a key chemical species that is found in a wide range of planetary atmospheres. In the context of exoplanets, CO2 is an indicator of the metal enrichment (i.e., elements heavier than helium, also called "metallicity"), and thus formation processes of the primary atmospheres of hot gas giants. It is also one of the most promising species to detect in the secondary atmospheres of terrestrial exoplanets. Previous photometric measurements of transiting planets with the Spitzer Space Telescope have given hints of the presence of CO2 but have not yielded definitive detections due to the lack of unambiguous spectroscopic identification. Here we present the detection of CO2 in the atmosphere of the gas giant exoplanet WASP-39b from transmission spectroscopy observations obtained with JWST as part of the Early Release Science Program (ERS). The data used in this study span 3.0 to 5.5 {\mu}m in wavelength and show a prominent CO2 absorption feature at 4.3 {\mu}m (26{\sigma} significance). The overall spectrum is well matched by one-dimensional, 10x solar metallicity models that assume radiative-convective-thermochemical equilibrium and have moderate cloud opacity. These models predict that the atmosphere should have water, carbon monoxide, and hydrogen sulfide in addition to CO2, but little methane. Furthermore, we also tentatively detect a small absorption feature near 4.0 {\mu}m that is not reproduced by these models., Comment: 27 pages, 6 figures, Accepted for publication in Nature, data and models available at https://doi.10.5281/zenodo.6959427

Upcoming James Webb Space Telescope (JWST) observations will allow us to study exoplanet and brown dwarf atmospheres in great detail. The physical interpretation of these upcoming high signal-to-noise observations requires precise atmospheric models of exoplanets and brown dwarfs. While several one-dimensional and three-dimensional atmospheric models have been developed in the past three decades, these models have often relied on simplified assumptions like chemical equilibrium and are also often not open-source, which limits their usage and development by the wider community. We present a python-based one-dimensional atmospheric radiative-convective equilibrium model. This model has heritage from the Fortran-based code (Marley et al.,1996} which has been widely used to model the atmospheres of Solar System objects, brown dwarfs, and exoplanets. In short, the basic capability of the original model is to compute the atmospheric state of the object under radiative-convective equilibrium given its effective or internal temperature, gravity, and host--star properties (if relevant). In the new model, which has been included within the well-utilized code-base PICASO, we have added these original features as well as the new capability of self-consistently treating disequilibrium chemistry. This code is widely applicable to Hydrogen-dominated atmospheres (e.g., brown dwarfs and giant planets)., Comment: Under review in ApJ; posting after incorporating the first referee response; official code release will coincide with paper acceptance; comments welcome; 29 Pages; 15 Figures; 3 Tables

We assess the prospect of using ring seismology to probe the interiors of the ice giants Uranus and Neptune. We do this by calculating normal mode spectra for different interior models of Uranus and Neptune using the stellar oscillation code GYRE. These spectra provide predictions of where in these planets' ring systems the effects of interior oscillations might be detected. We find that f-mode resonances with azimuthal order $m=2$ or $7 \leq m \leq 19$ fall among the inner rings (6, 5, 4, $\alpha$, and $\beta$) of Uranus, while f-mode resonances with $2 \leq m \leq 12$ fall in the tenuous $\zeta$ ring region. In addition, f-mode resonances with $m=2$ or $6 \leq m \leq 13$ may give azimuthal structure to Neptune's tenuous Galle ring. We also find that g-mode resonances may fall in the middle to outer rings of these planets. Although an orbiter is most likely required to confirm the association between any waves in the rings and planetary normal modes, the diversity of normal mode spectra implies that identification of just one or two modes in the rings of Uranus or Neptune would eliminate a variety of interior models, and thus aid in the interpretation of Voyager observations and future spacecraft measurements., Comment: 27 pages, 8 figures, accepted for publication in The Planetary Science Journal

The spectra of brown dwarfs are key to exploring the chemistry and physics that take place in their atmospheres. Late-T dwarf spectra are particularly diagnostic due to their relatively cloud-free atmospheres and deep molecular bands. With the use of powerful atmospheric retrieval tools applied to the spectra of these objects, direct constraints on molecular/atomic abundances, gravity, and vertical thermal profiles can be obtained enabling a broad exploration of the chemical/physical mechanisms operating in their atmospheres. We present a uniform retrieval analysis on low-resolution IRTF SpeX near-IR spectra of a sample of 50 T dwarfs, including new observations as part of a recent volume-limited survey. This analysis more than quadruples the sample of T dwarfs with retrieved temperature profiles and abundances (H$_2$O, CH$_4$, NH$_3$, K and subsequent C/O and metallicities). We are generally able to constrain effective temperatures to within 50K, volume mixing ratios for major species to within 0.25dex, atmospheric metallicities [M/H] to within 0.2, and C/O ratios to within 0.2. We compare our retrieved constraints on the thermal structure, chemistry, and gravities of these objects with predictions from self-consistent radiative-convective equilibrium models and find, in general though with substantial scatter, consistency with solar composition chemistry and thermal profiles of the neighboring stellar FGK population. Objects with notable discrepancies between the two modeling techniques and potential mechanisms for their differences, be they related to modeling approach or physically motivated, are discussed more thoroughly in the text., Comment: 21 pages, 11 figures

The direct characterization of exoplanetary systems with high contrast imaging is among the highest priorities for the broader exoplanet community. As large space missions will be necessary for detecting and characterizing exo-Earth twins, developing the techniques and technology for direct imaging of exoplanets is a driving focus for the community. For the first time, JWST will directly observe extrasolar planets at mid-infrared wavelengths beyond 5$\mu$m, deliver detailed spectroscopy revealing much more precise chemical abundances and atmospheric conditions, and provide sensitivity to analogs of our solar system ice-giant planets at wide orbital separations, an entirely new class of exoplanet. However, in order to maximise the scientific output over the lifetime of the mission, an exquisite understanding of the instrumental performance of JWST is needed as early in the mission as possible. In this paper, we describe our 55-hour Early Release Science Program that will utilize all four JWST instruments to extend the characterisation of planetary mass companions to $\sim$15$\mu$m as well as image a circumstellar disk in the mid-infrared with unprecedented sensitivity. Our program will also assess the performance of the observatory in the key modes expected to be commonly used for exoplanet direct imaging and spectroscopy, optimize data calibration and processing, and generate representative datasets that will enable a broad user base to effectively plan for general observing programs in future cycles., Comment: 21 pages, 7 figures. Accepted for Publication in PASP

The temperature profile of a planetary atmosphere is a key diagnostic of radiative and dynamical processes governing the absorption, redistribution, and emission of energy. Observations have revealed dayside stratospheres that either cool or warm with altitude for a small number of gas giant exoplanets, while other dayside stratospheres are consistent with constant temperatures. Here we report spectroscopic phase curve measurements for the gas giant WASP-121b, which constrain stratospheric temperatures throughout the diurnal cycle. Variations measured for a water vapour spectral feature reveal a temperature profile that transitions from warming with altitude on the dayside hemisphere to cooling with altitude on the nightside hemisphere. The data are well explained by models assuming chemical equilibrium, with water molecules thermally dissociating at low pressures on the dayside and recombining on the nightside. Nightside temperatures are low enough for perovskite (CaTiO3) to condense, which could deplete titanium from the gas phase and explain recent non-detections at the day-night terminator. Nightside temperatures are also consistent with the condensation of refractory species such as magnesium, iron, and vanadium. Detections of these metals at the day-night terminator suggest, however, that if they do form nightside clouds, cold trapping does not efficiently remove them from the upper atmosphere. Horizontal winds and vertical mixing could keep these refractory condensates aloft in the upper atmosphere of the nightside hemisphere until they are recirculated to the hotter dayside hemisphere and vaporised., Comment: Published in Nature Astronomy (publisher version is Open Access)

Young self-luminous giant exoplanets are expected to be oblate in shape owing to the high rotational speeds observed for some objects. Similar to the case of brown dwarfs, the thermal emission from these planets should be polarized by scatterings of molecules and condensate cloud particles, and the rotation-induced asymmetry of the planet's disk would yield to net non-zero detectable polarization. Considering an anisotropic atmosphere, we present here a three-dimensional approach to estimate the disk-averaged polarization that arises due to the oblateness of the planets. We solve the multiple-scattering vector radiative transfer equations at each location on the planet's disk and calculate the local Stokes vectors and then calculate the disk-integrated flux and linear polarization. For a cloud-free atmosphere, the polarization signal is observable only in the visible wavelength region. However, the presence of clouds in the planetary atmospheres leads to a detectable amount of polarization in the infrared wavelength region where the planetary thermal emission peaks. Considering different broad-band filters of the SPHERE-IRDIS instrument of the Very Large Telescope, we present generic models for the polarization at different wavelength bands as a function of their rotation period. We also present polarization models for the Exoplanets $\beta$ Pic b and ROXs 42B b as two representative cases which can guide future observations. Our insights on the polarization of young giant planets presented here would be useful for the upcoming polarimetric observations of the directly imaged planets., Comment: 18 pages, 13 figures, Accepted for publication in The Astrophysical journal

We present a survey for photometric variability in young, low-mass brown dwarfs with the Spitzer Space Telescope. The 23 objects in our sample show robust signatures of youth and share properties with directly-imaged exoplanets. We present three new young objects: 2MASS J03492367$+$0635078, 2MASS J09512690 $-$8023553 and 2MASS J07180871$-$6415310. We detect variability in 13 young objects, and find that young brown dwarfs are highly likely to display variability across the L2--T4 spectral type range. In contrast, the field dwarf variability occurrence rate drops for spectral types $>$L9. We examine the variability amplitudes of young objects and find an enhancement in maximum amplitudes compared to field dwarfs. We speculate that the observed range of amplitudes within a spectral type may be influenced by secondary effects such as viewing inclination and/or rotation period. We combine our new rotation periods with the literature to investigate the effects of mass on angular momentum evolution. While high mass brown dwarfs ($>30 M_{\mathrm{Jup}}$) spin up over time, the same trend is not apparent for lower mass objects ($<30 M_{\mathrm{Jup}}$), likely due to the small number of measured periods for old, low-mass objects. The rotation periods of companion brown dwarfs and planetary-mass objects are consistent with those of isolated objects with similar ages and masses, suggesting similar angular momentum histories. Within the AB Doradus group, we find a high variability occurrence rate and evidence for common angular momentum evolution. The results are encouraging for future variability searches in directly-imaged exoplanets with facilities such as the James Webb Space Telescope and 30-meter telescopes.

We present a near-infrared transmission spectrum of the long period (P=542 days), temperate ($T_{eq}$=294 K) giant planet HIP 41378 f obtained with the Wide-Field Camera 3 (WFC3) instrument aboard the Hubble Space Telescope (HST). With a measured mass of 12 $\pm$ 3 $M_{\oplus}$ and a radius of 9.2 $\pm$ 0.1 $R_{\oplus}$, HIP 41378 f has an extremely low bulk density (0.09 $\pm$ 0.02 g/cm$^{3}$). We measure the transit depth with a median precision of 84 ppm in 30 spectrophotometric channels with uniformly-sized widths of 0.018 microns. Within this level of precision, the spectrum shows no evidence of absorption from gaseous molecular features between 1.1-1.7 microns. Comparing the observed transmission spectrum to a suite of 1D radiative-convective-thermochemical-equilibrium forward models, we rule out clear, low-metallicity atmospheres and find that the data prefer high-metallicity atmospheres or models with an additional opacity source such as high-altitude hazes and/or circumplanetary rings. We explore the ringed scenario for this planet further by jointly fitting the K2 and HST light curves to constrain the properties of putative rings. We also assess the possibility of distinguishing between hazy, ringed, and high-metallicity scenarios at longer wavelengths with JWST. HIP 41378 f provides a rare opportunity to probe the atmospheric composition of a cool giant planet spanning the gap between the Solar System giants, directly imaged planets, and the highly-irradiated hot Jupiters traditionally studied via transit spectroscopy., Comment: 12 pages, 4 figures, accepted to ApJL. Please also see companion paper Ohno & Fortney (2022) for further details on the ring modeling methodology (arxiv DOI: 2201.02794)