Showing 1–50 of 265 results for author: Marley, M S

JADES: Spectroscopic Confirmation and Proper Motion for a T-Dwarf at 2 Kiloparsecs

Authors: Kevin N. Hainline, Francesco D'Eugenio, Fengwu Sun, Jakob M. Helton, Brittany E. Miles, MArk S. Marley, Ben W. P. Lew, Jarron M. Leisenring, Andrew J. Bunker, Phillip A. Cargile, Stefano Carniani, Daniel J. Eisenstein, Ignas Juodzbalis, Benjamin D. Johnson, Brant Robertson, Sandro Tacchella, Christina C. Williams, Christopher N. A. Willmer

Abstract: 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 sta… ▽ More 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. △ Less

Submitted 11 July, 2024; originally announced July 2024.

Comments: 10 pages, 4 figures, submitted to AAS Journals

Precise Bolometric Luminosities and Effective Temperatures of 23 late-T and Y dwarfs Obtained with JWST

Authors: Samuel A. Beiler, Michael C. Cushing, J. Davy Kirkpatrick, Adam C. Schneider, Sagnick Mukherjee, Mark S. Marley, Federico Marocco, Richard L. Smart

Abstract: 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 $μ$m wavelength range at $λ/ Δλ\approx 100$ and broadband photometry at 15, 18, and 21 $μ$m. The spectra exhibit absorption features common to these objects including H… ▽ More 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 $μ$m wavelength range at $λ/ Δλ\approx 100$ and broadband photometry at 15, 18, and 21 $μ$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 $λ< 3$ $μ$m and $λ> 8$ $μ$m, it shows no clear trend in the 5 $μ$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. △ Less

Submitted 11 July, 2024; originally announced July 2024.

Comments: Accepted, 33 pages, 19 figures, 8 tables

Influence of Orbit and Mass Constraints on Reflected Light Characterization of Directly Imaged Rocky Exoplanets

Authors: Arnaud Salvador, Tyler D. Robinson, Jonathan J. Fortney, Mark S. Marley

Abstract: 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 pr… ▽ More 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. △ Less

Submitted 11 June, 2024; originally announced June 2024.

Comments: Accepted for publication in ApJL

Probing the Heights and Depths of Y Dwarf Atmospheres: A Retrieval Analysis of the JWST Spectral Energy Distribution of WISE J035934.06$-$540154.6

Authors: Harshil Kothari, Michael C. Cushing, Ben Burningham, Samuel A. Beiler, J. Davy Kirkpatrick, Adam C. Schneider, Sagnick Mukherjee, Mark S. Marley

Abstract: We present an atmospheric retrieval analysis of the Y0 brown dwarf WISE J035934.06$-$540154.6 using the low-resolution 0.96--12 $μ$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-… ▽ More We present an atmospheric retrieval analysis of the Y0 brown dwarf WISE J035934.06$-$540154.6 using the low-resolution 0.96--12 $μ$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. △ Less

Submitted 10 June, 2024; originally announced June 2024.

Time-resolved Hubble Space Telescope Wide Field Camera 3 Spectrophotometry Reveals Inefficient Day-to-Night Heat Redistribution in the Highly Irradiated Brown Dwarf SDSS 1557B

Authors: Rachael C. Amaro, Daniel Apai, Ben W. P. Lew, Yifan Zhou, Joshua D. Lothringer, Sarah L. Casewell, Xianyu Tan, Travis Barman, Mark S. Marley, L. C. Mayorga, Vivien Parmentier

Abstract: 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… ▽ More 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 $μ$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. △ Less

Submitted 11 April, 2024; originally announced April 2024.

Comments: 19 pages and 11 figures. Accepted to Astrophysical Journal

The Sonora Substellar Atmosphere Models. III. Diamondback: Atmospheric Properties, Spectra, and Evolution for Warm Cloudy Substellar Objects

Authors: Caroline V. Morley, Sagnick Mukherjee, Mark S. Marley, Jonathan J. Fortney, Channon Visscher, Roxana Lupu, Ehsan Gharib-Nezhad, Daniel Thorngren, Richard Freedman, Natasha Batalha 7

Abstract: 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 pressur… ▽ More 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. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: 26 pages, 18 figures, submitted for publication in ApJ

The Sonora Substellar Atmosphere Models. IV. Elf Owl: Atmospheric Mixing and Chemical Disequilibrium with Varying Metallicity and C/O Ratios

Authors: Sagnick Mukherjee, Jonathan J. Fortney, Caroline V. Morley, Natasha E. Batalha, Mark S. Marley, Theodora Karalidi, Channon Visscher, Roxana Lupu, Richard Freedman, Ehsan Gharib-Nezhad

Abstract: 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, at… ▽ More 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 $μ$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. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: Accepted for Publication in The Astrophysical Journal, 16 Figures, 3 Tables, 28 Pages

The JWST Early Release Science Program for Direct Observations of Exoplanetary Systems V: Do Self-Consistent Atmospheric Models Represent JWST Spectra? A Showcase With VHS 1256 b

Authors: Simon Petrus, Niall Whiteford, Polychronis Patapis, Beth A. Biller, Andrew Skemer, Sasha Hinkley, Genaro Suárez, Anna Lueber, Paulina Palma-Bifani, Jordan M. Stone, Johanna M. Vos, Caroline V. Morley, Pascal Tremblin, Benjamin Charnay, Christiane Helling, Brittany E. Miles, Aarynn L. Carter, Jason J. Wang, Markus Janson, Eileen C. Gonzales, Ben Sutlieff, Kielan K. W. Hoch, Mickaël Bonnefoy, Gaël Chauvin, Olivier Absil , et al. (97 additional authors not shown)

Abstract: 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. W… ▽ More 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. △ Less

Submitted 31 January, 2024; v1 submitted 6 December, 2023; originally announced December 2023.

Comments: 32 pages, 16 figures, 6 tables, 2 appendices

The \textit{JWST} Early Release Science Program for Direct Observations of Exoplanetary Systems III: Aperture Masking Interferometric Observations of the star HIP\,65426 at $\boldsymbol{3.8\,\rm{μm}}$

Authors: Shrishmoy Ray, Steph Sallum, Sasha Hinkley, Anand Sivamarakrishnan, Rachel Cooper, Jens Kammerer, Alexandra Z. Greebaum, Deepashri Thatte, Cecilia Lazzoni, Andrei Tokovinin, Matthew de Furio, Samuel Factor, Michael Meyer, Jordan M. Stone, Aarynn Carter, Beth Biller, Andrew Skemer, Genaro Suarez, Jarron M. Leisenring, Marshall D. Perrin, Adam L. Kraus, Olivier Absil, William O. Balmer, Mickael Bonnefoy, Marta L. Bryan , et al. (98 additional authors not shown)

Abstract: We present aperture masking interferometry (AMI) observations of the star HIP 65426 at $3.8\,\rm{μm}$ as a part of the \textit{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λ/D$ f… ▽ More We present aperture masking interferometry (AMI) observations of the star HIP 65426 at $3.8\,\rm{μm}$ as a part of the \textit{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λ/D$ for an interferometer), which are inaccessible with the classical inner working angles of the \textit{JWST} coronagraphs. When combined with \textit{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 contrast of $Δm_{F380M}{\sim }7.8$\,mag relative to the host star at a separation of ${\sim}0.07\arcsec$ but detect no additional companions interior to the known companion HIP\,65426\,b. Our observations thus rule out companions more massive than $10{-}12\,\rm{M\textsubscript{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 \textit{JWST} is sensitive to planetary mass companions orbiting at the water frost line, even for more distant stars at $\sim$100\,pc. This result will allow the planning and successful execution of future observations to probe the inner regions of nearby stellar systems, opening essentially unexplored parameter space. △ Less

Submitted 17 October, 2023; originally announced October 2023.

Comments: 15 pages, 9 figures, submitted to ApJ Letters

The JWST Early Release Science Program for Direct Observations of Exoplanetary Systems IV: NIRISS Aperture Masking Interferometry Performance and Lessons Learned

Authors: Steph Sallum, Shrishmoy Ray, Jens Kammerer, Anand Sivaramakrishnan, Rachel Cooper, Alexandra Z. Greebaum, Deepashri Thatte, Matthew de Furio, Samuel Factor, Michael Meyer, Jordan M. Stone, Aarynn Carter, Beth Biller, Sasha Hinkley, Andrew Skemer, Genaro Suarez, Jarron M. Leisenring, Marshall D. Perrin, Adam L. Kraus, Olivier Absil, William O. Balmer, Mickael Bonnefoy, Marta L. Bryan, Sarah K. Betti, Anthony Boccaletti , et al. (98 additional authors not shown)

Abstract: 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… ▽ More 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 λ/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. △ Less

Submitted 11 March, 2024; v1 submitted 17 October, 2023; originally announced October 2023.

Comments: 20 pages, 12 figures, accepted to Astrophysical Journal Letters

ELemental abundances of Planets and brown dwarfs Imaged around Stars (ELPIS): I. Potential Metal Enrichment of the Exoplanet AF Lep b and a Novel Retrieval Approach for Cloudy Self-luminous Atmospheres

Authors: Zhoujian Zhang, Paul Mollière, Keith Hawkins, Catherine Manea, Jonathan J. Fortney, Caroline V. Morley, Andrew Skemer, Mark S. Marley, Brendan P. Bowler, Aarynn L. Carter, Kyle Franson, Zachary G. Maas, Christopher Sneden

Abstract: 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.,… ▽ More 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$ $μ$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. △ Less

Submitted 5 September, 2023; originally announced September 2023.

Comments: 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

Probing reflection from aerosols with the near-infrared dayside spectrum of WASP-80b

Authors: Bob Jacobs, Jean-Michel Désert, Peter Gao, Caroline V. Morley, Jacob Arcangeli, Saugata Barat, Mark S. Marley, Julianne I. Moses, Jonathan J. Fortney, Jacob L. Bean, Kevin B. Stevenson, Vatsal Panwar

Abstract: 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… ▽ More 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σ$, which hints toward the presence of reflecting aerosols on this planet's dayside, indicating a geometric albedo of $A_g<0.33$ at 3$σ$. 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σ$. 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. △ Less

Submitted 26 October, 2023; v1 submitted 26 July, 2023; originally announced July 2023.

Comments: Published in ApJ Letters (20 Oct 2023)

Journal ref: ApJL 2023 Volume 956, Number 2, page L43

Saturn's Interior After the Cassini Grand Finale

Authors: J. J. Fortney, B. Militzer, C. R. Mankovich, R. Helled, S. M. Wahl, N. Nettelmann, W. B. Hubbard, D. J. Stevenson, L. Iess, M. S. Marley, N. Movshovitz

Abstract: We present a review of Saturn's interior structure and thermal evolution, with a particular focus on work in the past 5 years. Data from the Cassini mission, including a precise determination of the gravity field from the Grand Finale orbits, and the still ongoing identification of ring wave features in Saturn's C-ring tied to seismic modes in the planet, have led to dramatic advances in our under… ▽ More We present a review of Saturn's interior structure and thermal evolution, with a particular focus on work in the past 5 years. Data from the Cassini mission, including a precise determination of the gravity field from the Grand Finale orbits, and the still ongoing identification of ring wave features in Saturn's C-ring tied to seismic modes in the planet, have led to dramatic advances in our understanding of Saturn's structure. Models that match the gravity field suggest that differential rotation, as seen in the visible atmosphere, extends down to at least a depth of 10,000 km (1/6$^{\rm th}$ the planet's radius). At greater depths, a variety of different investigations all now point to a deep Saturn rotation rate of 10 hours and 33 minutes. There is very compelling evidence for a central heavy element concentration (``core''), that in most recent models is 12-20 Earth masses. Ring seismology strongly suggests that the core is not entirely compact, but is dilute (mixed in with the overlying H/He), and has a substantial radial extent, perhaps out to around one-half of the planet's radius. A wide range of thermal evolution scenarios can match the planet's current luminosity, with progress on better quantifying the helium rain scenario hampered by Saturn's poorly known atmospheric helium abundance. We discuss the relevance of magnetic field data on understanding the planet's current interior structure. We point towards additional future work that combines seismology and gravity within a framework that includes differential rotation, and the utility of a Saturn entry probe. △ Less

Submitted 18 April, 2023; originally announced April 2023.

Comments: Invited review. Accepted for publication in "Saturn: The Grand Finale", K. H. Baines et al., eds., Cambridge University Press. All-new follow-up to previous 2016 (pre-Grand Finale) review chapter here: arXiv:1609.06324

Spherical Harmonics for the 1D Radiative Transfer Equation II: Thermal Emission

Authors: Caoimhe M. Rooney, Natasha E. Batalha, Mark S. Marley

Abstract: 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 fo… ▽ More 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. △ Less

Submitted 10 April, 2023; originally announced April 2023.

Comments: 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;

Spherical Harmonics for the 1D Radiative Transfer Equation I: Reflected Light

Authors: Caoimhe M. Rooney, Natasha E. Batalha, Mark S. Marley

Abstract: 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:… ▽ More 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). △ Less

Submitted 10 April, 2023; originally announced April 2023.

Comments: 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

The puzzle of the formation of T8 dwarf Ross 458c

Authors: Josefine Gaarn, Ben Burningham, Jacqueline K. Faherty, Channon Visscher, Mark S. Marley, Eileen C. Gonzales, Emily Calamari, Daniella Bardalez Gagliuffi, Roxana Lupu, Richard Freedman

Abstract: 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… ▽ More 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. △ Less

Submitted 3 April, 2023; v1 submitted 29 March, 2023; originally announced March 2023.

Comments: 16 pages, 9 figures

Hotter Than Expected: HST/WFC3 Phase-resolved Spectroscopy of a Rare Irradiated Brown Dwarf with Strong Internal Heat Flux

Authors: Rachael C. Amaro, Daniel Apai, Yifan Zhou, Ben W. P. Lew, Sarah L. Casewell, L. Mayorga, Mark S. Marley, Xianyu Tan, Joshua D. Lothringer, Vivien Parmentier, Travis Barman

Abstract: 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 HS… ▽ More 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. △ Less

Submitted 13 March, 2023; originally announced March 2023.

Comments: 31 pages, 25 figures, Accepted to Astrophysical Journal with video abstract. Video abstract available at https://www.rachaelamaro.com/research-1

Improved companion mass limits for Sirius A with thermal infrared coronagraphy using a vector-apodizing phase plate and time-domain starlight-subtraction techniques

Authors: Joseph D. Long, Jared R. Males, Sebastiaan Y. Haffert, Logan Pearce, Mark S. Marley, Katie M. Morzinski, Laird M. Close, Gilles P. P. L. Otten, Frans Snik, Matthew A. Kenworthy, Christoph U. Keller, Philip Hinz, John D. Monnier, Alycia Weinberger, Volker Tolls

Abstract: We use observations with the infrared-optimized MagAO system and Clio camera in 3.9 $μ$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… ▽ More We use observations with the infrared-optimized MagAO system and Clio camera in 3.9 $μ$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. △ Less

Submitted 9 March, 2023; originally announced March 2023.

Comments: 19 pages, 22 figures, accepted to AJ

Testing the Interaction Between a Substellar Companion and a Debris Disk in the HR 2562 System

Authors: Stella Yimiao Zhang, Gaspard Duchêne, Robert J. De Rosa, Megan Ansdell, Quinn Konopacky, Thomas Esposito, Eugene Chiang, Malena Rice, Brenda Matthews, Paul Kalas, Bruce Macintosh, Franck Marchis, Stan Metchev, Jenny Patience, Julien Rameau, Kimberly Ward-Duong, Schuyler Wolff, Michael P. Fitzgerald, Vanessa P. Bailey, Travis S. Barman, Joanna Bulger, Christine H. Chen, Jeffrey K. Chilcotte, Tara Cotten, René Doyon , et al. (29 additional authors not shown)

Abstract: 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… ▽ More 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$σ$ 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. △ Less

Submitted 9 February, 2023; originally announced February 2023.

Comments: 22 pages, 17 figures

Optical Properties of Organic Haze Analogues in Water-rich Exoplanet Atmospheres Observable with JWST

Authors: Chao He, Michael Radke, Sarah E. Moran, Sarah M. Horst, Nikole K. Lewis, Julianne I. Moses, Mark S. Marley, Natasha E. Batalha, Eliza M. -R. Kempton, Caroline V. Morley, Jeff A. Valenti, Veronique Vuitton

Abstract: 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 o… ▽ More 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 experimental data are not available. Here we measure the density and optical properties of organic haze analogues generated in water-rich exoplanet atmosphere experiments. We report optical constants (0.4 to 28.6 μm) of organic haze analogues 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. △ Less

Submitted 28 November, 2023; v1 submitted 6 January, 2023; originally announced January 2023.

Comments: 4 figures, 1 Table, Published in Nature Astronomy

Journal ref: Nature Astronomy (2023)

The JWST Early Release Science Program for Direct Observations of Exoplanetary Systems II: A 1 to 20 Micron Spectrum of the Planetary-Mass Companion VHS 1256-1257 b

Authors: Brittany E. Miles, Beth A. Biller, Polychronis Patapis, Kadin Worthen, Emily Rickman, Kielan K. W. Hoch, Andrew Skemer, Marshall D. Perrin, Niall Whiteford, Christine H. Chen, B. Sargent, Sagnick Mukherjee, Caroline V. Morley, Sarah E. Moran, Mickael Bonnefoy, Simon Petrus, Aarynn L. Carter, Elodie Choquet, Sasha Hinkley, Kimberly Ward-Duong, Jarron M. Leisenring, Maxwell A. Millar-Blanchaer, Laurent Pueyo, Shrishmoy Ray, Karl R. Stapelfeldt , et al. (79 additional authors not shown)

Abstract: 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… ▽ More 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 $μ$m to 20 $μ$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. △ Less

Submitted 4 July, 2024; v1 submitted 1 September, 2022; originally announced September 2022.

Comments: 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

The JWST Early Release Science Program for Direct Observations of Exoplanetary Systems I: High Contrast Imaging of the Exoplanet HIP 65426 b from 2-16 $μ$m

Authors: Aarynn L. Carter, Sasha Hinkley, Jens Kammerer, Andrew Skemer, Beth A. Biller, Jarron M. Leisenring, Maxwell A. Millar-Blanchaer, Simon Petrus, Jordan M. Stone, Kimberly Ward-Duong, Jason J. Wang, Julien H. Girard, Dean C. Hines, Marshall D. Perrin, Laurent Pueyo, William O. Balmer, Mariangela Bonavita, Mickael Bonnefoy, Gael Chauvin, Elodie Choquet, Valentin Christiaens, Camilla Danielski, Grant M. Kennedy, Elisabeth C. Matthews, Brittany E. Miles , et al. (86 additional authors not shown)

Abstract: 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 $μ$m, and with the Mid-Infrared Instrument (MIRI) from 11-16 $μ$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 exo… ▽ More 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 $μ$m, and with the Mid-Infrared Instrument (MIRI) from 11-16 $μ$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 $μ$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$σ$ contrast limits of $\sim$1$\times10^{-5}$ and $\sim$2$\times10^{-4}$ at 1" for NIRCam at 4.4 $μ$m and MIRI at 11.3 $μ$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 $μ$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. △ Less

Submitted 3 May, 2023; v1 submitted 31 August, 2022; originally announced August 2022.

Comments: 35 pages, 16 figures, 4 tables, 1 wonderful telescope; Submitted to AAS Journals

Probing the Extent of Vertical Mixing in Brown Dwarf Atmospheres with Disequilibrium Chemistry

Authors: Sagnick Mukherjee, Jonathan J. Fortney, Natasha E. Batalha, Theodora Karalidi, Mark S. Marley, Channon Visscher, Brittany E. Miles, Andrew J. I. Skemer

Abstract: 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 convecti… ▽ More 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. △ Less

Submitted 30 August, 2022; originally announced August 2022.

Comments: Accepted for Publication in The Astrophysical Journal, 31 Pages, 22 Figures. The model grid will be released via Zenodo

Identification of carbon dioxide in an exoplanet atmosphere

Authors: The JWST Transiting Exoplanet Community Early Release Science Team, Eva-Maria Ahrer, Lili Alderson, Natalie M. Batalha, Natasha E. Batalha, Jacob L. Bean, Thomas G. Beatty, Taylor J. Bell, Björn Benneke, Zachory K. Berta-Thompson, Aarynn L. Carter, Ian J. M. Crossfield, Néstor Espinoza, Adina D. Feinstein, Jonathan J. Fortney, Neale P. Gibson, Jayesh M. Goyal, Eliza M. -R. Kempton, James Kirk, Laura Kreidberg, Mercedes López-Morales, Michael R. Line, Joshua D. Lothringer, Sarah E. Moran, Sagnick Mukherjee , et al. (107 additional authors not shown)

Abstract: 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… ▽ More 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 μm in wavelength and show a prominent CO2 absorption feature at 4.3 μm (26σ 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 μm that is not reproduced by these models. △ Less

Submitted 24 August, 2022; originally announced August 2022.

Comments: 27 pages, 6 figures, Accepted for publication in Nature, data and models available at https://doi.10.5281/zenodo.6959427

PICASO 3.0: A One-Dimensional Climate Model for Giant Planets and Brown Dwarfs

Authors: Sagnick Mukherjee, Natasha E. Batalha, Jonathan J. Fortney, Mark S. Marley

Abstract: 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,… ▽ More 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). △ Less

Submitted 16 August, 2022; originally announced August 2022.

Comments: 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

Ring Seismology of the Ice Giants Uranus and Neptune

Authors: Joseph A. A'Hearn, Matthew M. Hedman, Christopher R. Mankovich, Hima Aramona, Mark S. Marley

Abstract: 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 re… ▽ More 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, $α$, and $β$) of Uranus, while f-mode resonances with $2 \leq m \leq 12$ fall in the tenuous $ζ$ 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. △ Less

Submitted 10 June, 2022; originally announced June 2022.

Comments: 27 pages, 8 figures, accepted for publication in The Planetary Science Journal

A Uniform Retrieval Analysis of Ultra-cool Dwarfs. IV. A Statistical Census from 50 Late-T Dwarfs

Authors: Joseph A Zalesky, Kezman Saboi, Michael R. Line, Zhoujian Zhang, Adam C Schneider, Michael C Liu, William M J Best, Mark S Marley

Abstract: 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 vertica… ▽ More 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. △ Less

Submitted 2 June, 2022; originally announced June 2022.

Comments: 21 pages, 11 figures

Diurnal variations in the stratosphere of the ultrahot giant exoplanet WASP-121b

Authors: Thomas Mikal-Evans, David K. Sing, Joanna K. Barstow, Tiffany Kataria, Jayesh Goyal, Nikole Lewis, Jake Taylor, Nathan. J. Mayne, Tansu Daylan, Hannah R. Wakeford, Mark S. Marley, Jessica J. Spake

Abstract: 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 spe… ▽ More 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. △ Less

Submitted 20 February, 2022; originally announced February 2022.

Comments: Published in Nature Astronomy (publisher version is Open Access)

Polarization of Rotationally Oblate Self-Luminous Exoplanets with Anisotropic Atmospheres

Authors: Aritra Chakrabarty, Sujan Sengupta, Mark S. Marley

Abstract: 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.… ▽ More 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 $β$ 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. △ Less

Submitted 19 January, 2022; originally announced January 2022.

Comments: 18 pages, 13 figures, Accepted for publication in The Astrophysical journal

HST/WFC3 Complete Phase-resolved Spectroscopy of White Dwarf-Brown Dwarf Binaries WD 0137 and EPIC 2122

Authors: Yifan Zhou, Dániel Apai, Xianyu Tan, Joshua D. Lothringer, Ben W. P. Lew, Sarah L. Casewell, Vivien Parmentier, Mark S. Marley, Siyi Xu, L. C. Mayorga

Abstract: Brown dwarfs in close-in orbits around white dwarfs offer an excellent opportunity to investigate properties of fast-rotating, tidally-locked, and highly-irradiated atmospheres. We present Hubble Space Telescope Wide Field Camera 3 G141 phase-resolved observations of two brown dwarf-white dwarf binaries: WD 0137-349 and EPIC 212235321. Their 1.1 to 1.7 $μ$m phase curves demonstrate rotational modu… ▽ More Brown dwarfs in close-in orbits around white dwarfs offer an excellent opportunity to investigate properties of fast-rotating, tidally-locked, and highly-irradiated atmospheres. We present Hubble Space Telescope Wide Field Camera 3 G141 phase-resolved observations of two brown dwarf-white dwarf binaries: WD 0137-349 and EPIC 212235321. Their 1.1 to 1.7 $μ$m phase curves demonstrate rotational modulations with semi-amplitudes of $5.27\pm0.02$% and $29.1\pm0.1$%; both can be well fit by multi-order Fourier series models. The high-order Fourier components have the same phase as the first order and are likely caused by hot spots located at the substellar points, suggesting inefficient day/night heat transfer. Both brown dwarfs' phase-resolved spectra can be accurately represented by linear combinations of their day- and night-side spectra. Fitting the irradiated brown dwarf model grids to the day-side spectra require a filling factor of ~50%, further supporting a hot spot dominating the emission of the day-sides. The night-side spectrum of WD 0137-349B is reasonably well fit by non-irradiated substellar models and the one of EPIC 212235321B can be approximated by a Planck function. We find strong spectral variations in the brown dwarfs' day/night flux and brightness temperature contrasts, which highlights the limitations of band-integrated measurements in probing heat transfer in irradiated objects. On the color-magnitude diagram, WD 0137-349B evolves along a cloudless model track connecting the early-L and mid-T spectral types, demonstrating that clouds and disequilibrium chemistry have a negligible effect on this object. A full interpretation of these high-quality phase-resolved spectra calls for new models that couple atmospheric circulation and radiative transfer under high-irradiation conditions. △ Less

Submitted 19 October, 2021; originally announced October 2021.

Comments: Accepted for publication in the AAS Journals. Please also see a companion paper led by Ben W. P. Lew, in which we discuss results from another set of observations of the same program. A video abstract of the paper can be viewed at https://youtu.be/WNz4iBKO99M

A New Sedimentation Model for Greater Cloud Diversity in Giant Exoplanets and Brown Dwarfs

Authors: Caoimhe M. Rooney, Natasha E. Batalha, Peter Gao, Mark S. Marley

Abstract: The observed atmospheric spectrum of exoplanets and brown dwarfs depends critically on the presence and distribution of atmospheric condensates. The Ackerman & Marley (2001) methodology for predicting the vertical distribution of condensate particles is widely used to study cloudy atmospheres and has recently been implemented in an open-source python package virga. The model relies upon input para… ▽ More The observed atmospheric spectrum of exoplanets and brown dwarfs depends critically on the presence and distribution of atmospheric condensates. The Ackerman & Marley (2001) methodology for predicting the vertical distribution of condensate particles is widely used to study cloudy atmospheres and has recently been implemented in an open-source python package virga. The model relies upon input parameter $f_{\text{sed}}$, the sedimentation efficiency, which until now has been held constant. The relative simplicity of this model renders it useful for retrieval studies due to its rapidly attainable solutions. However, comparisons with more complex microphysical models such as CARMA have highlighted inconsistencies between the two approaches, namely that the cloud parameters needed for radiative transfer produced by virga are dissimilar to those produced by CARMA. To address these discrepancies, we have extended the original Ackerman and Marley methodology in virga to allow for non-constant $f_{\text{sed}}$ values, in particular those that vary with altitude. We discuss one such parameterization and compare the cloud mass mixing ratio produced by virga with constant and variable $f_{\text{sed}}$ profiles to that produced by CARMA. We find that the variable $f_{\text{sed}}$ formulation better captures the profile produced by CARMA with heterogeneous nucleation, yet performs comparatively to constant $f_{\text{sed}}$ for homogeneous nucleation. In general, virga has the capacity to handle any $f_{\text{sed}}$ with an explicit anti-derivative, permitting a plethora of alternative cloud profiles that are otherwise unattainable by constant $f_{\text{sed}}$ values. The ensuing flexibility has the potential to better agree with increasingly complex models and observed data. △ Less

Submitted 26 October, 2021; v1 submitted 12 October, 2021; originally announced October 2021.

Comments: Accepted ApJ; GitHub: http://github.com/natashabatalha/virga Docs: http://natashabatalha.github.io/virga/# Virga Tutorials: http://natashabatalha.github.io/virga/tutorials.html

Modeling Polarization Signals from Cloudy Brown Dwarfs: Luhman 16 A and B in Three Dimensions

Authors: Sagnick Mukherjee, Jonathan J. Fortney, Rebecca Jensen-Clem, Xianyu Tan, Mark S. Marley, Natasha E. Batalha

Abstract: The detection of disk-integrated polarization from Luhman 16A and B in H-band, and subsequent modeling, has been interpreted in the framework of zonal cloud bands on these bodies. Recently, Tan and Showman (2021) investigated three-dimensional atmospheric circulation and cloud structures of brown dwarfs with general circulation models (GCMs), and their simulations yield complex cloud distributions… ▽ More The detection of disk-integrated polarization from Luhman 16A and B in H-band, and subsequent modeling, has been interpreted in the framework of zonal cloud bands on these bodies. Recently, Tan and Showman (2021) investigated three-dimensional atmospheric circulation and cloud structures of brown dwarfs with general circulation models (GCMs), and their simulations yield complex cloud distributions showing some aspects of zonal jets, but also complex vortices that cannot be captured by a simple model. Here we use these 3D GCMs specific to Luhman 16A and B, along with the three-dimensional Monte Carlo radiative transfer code ARTES, to calculate their polarization signals. We adopt the 3D temperature-pressure and cloud profiles from the GCMs as our input atmospheric structures. Our polarization calculations at 1.6 $μ$m agree well with the measured degree of linear polarization from both Luhman 16 A and B. Our calculations reproduce the measured polarization for both the objects with cloud particle sizes between 0.5-1 \,$μ$m for Luhman 16 A and 5 \,$μ$m for Luhman 16 B. We find that the degree of linear polarization can vary on hour-long timescales over the course of a rotation period. We also show that models with azimuthally symmetric band-like cloud geometries, typically used for interpreting polarimetry observations of brown dwarfs, over-predict the polarization signal if the cloud patterns do not include complex vortices within these bands. This exploratory work shows that GCMs are promising for modeling and interpreting polarization signals of brown dwarfs. △ Less

Submitted 12 October, 2021; originally announced October 2021.

Comments: 13 pages, 6 figures, Accepted for publication in The Astrophysical journal

Ultracool Dwarfs Observed with the Spitzer Infrared Spectrograph. I. An Accurate Look at the L-to-T Transition at $\sim$300 Myr from Optical through Mid-infrared Spectrophotometry

Authors: Genaro Suárez, Stanimir Metchev, Sandy K. Leggett, Didier Saumon, Mark S. Marley

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

Submitted 29 July, 2021; originally announced July 2021.

Comments: Accepted for publication in ApJ

The Sonora Brown Dwarf Atmosphere and Evolution Models I. Model Description and Application to Cloudless Atmospheres in Rainout Chemical Equilibrium

Authors: Mark S. Marley, Didier Saumon, Channon Visscher, Roxana Lupu, Richard Freedman, Caroline Morley, Jonathan J. Fortney, Christopher Seay, Adam J. R. W. Smith, D. J. Teal, Ruoyan Wang

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

Submitted 15 July, 2021; originally announced July 2021.

Comments: 27 pages, 16 figures, accepted for Astrophysical Journal. Models available at https://doi.org/10.5281/zenodo.5063476

Detection and Bulk Properties of the HR 8799 Planets with High Resolution Spectroscopy

Authors: Jason J. Wang, Jean-Baptiste Ruffio, Evan Morris, Jacques-Robert Delorme, Nemanja Jovanovic, Jacklyn Pezzato, Daniel Echeverri, Luke Finnerty, Callie Hood, J. J. Zanazzi, Marta L. Bryan, Charlotte Z. Bond, Sylvain Cetre, Emily C. Martin, Dimitri Mawet, Andy Skemer, Ashley Baker, Jerry W. Xuan, J. Kent Wallace, Ji Wang, Randall Bartos, Geoffrey A. Blake, Andy Boden, Cam Buzard, Benjamin Calvin , et al. (27 additional authors not shown)

Abstract: Using the Keck Planet Imager and Characterizer (KPIC), we obtained high-resolution (R$\sim$35,000) $K$-band spectra of the four planets orbiting HR 8799. We clearly detected \water{} and CO in the atmospheres of HR 8799 c, d, and e, and tentatively detected a combination of CO and \water{} in b. These are the most challenging directly imaged exoplanets that have been observed at high spectral reso… ▽ More Using the Keck Planet Imager and Characterizer (KPIC), we obtained high-resolution (R$\sim$35,000) $K$-band spectra of the four planets orbiting HR 8799. We clearly detected \water{} and CO in the atmospheres of HR 8799 c, d, and e, and tentatively detected a combination of CO and \water{} in b. These are the most challenging directly imaged exoplanets that have been observed at high spectral resolution to date when considering both their angular separations and flux ratios. We developed a forward modeling framework that allows us to jointly fit the spectra of the planets and the diffracted starlight simultaneously in a likelihood-based approach and obtained posterior probabilities on their effective temperatures, surface gravities, radial velocities, and spins. We measured $v\sin(i)$ values of $10.1^{+2.8}_{-2.7}$~km/s for HR 8799 d and $15.0^{+2.3}_{-2.6}$~km/s for HR 8799 e, and placed an upper limit of $< 14$~km/s of HR 8799 c. Under two different assumptions of their obliquities, we found tentative evidence that rotation velocity is anti-correlated with companion mass, which could indicate that magnetic braking with a circumplanetary disk at early times is less efficient at spinning down lower mass planets. △ Less

Submitted 14 July, 2021; originally announced July 2021.

Comments: 31 pages, 12 figures, Accepted to AJ

An Improved Near-Infrared Spectrum of the Archetype Y Dwarf WISEP J182831.08+265037.8

Authors: Michael C. Cushing, Adam C. Schneider, J. Davy Kirkpatrick, Caroline V. Morley, Mark S. Marley, Christopher R. Gelino, Gregory N. Mace, Edward L. Wright, Peter R. Eisenhardt, Michael F. Skrutskie, Kenneth A. Marsh

Abstract: We present a Hubble Space Telescope/Wide-Field Camera 3 near infrared spectrum of the archetype Y dwarf WISEP 182831.08+265037.8. The spectrum covers the 0.9-1.7 um wavelength range at a resolving power of lambda/Delta lambda ~180 and is a significant improvement over the previously published spectrum because it covers a broader wavelength range and is uncontaminated by light from a background sta… ▽ More We present a Hubble Space Telescope/Wide-Field Camera 3 near infrared spectrum of the archetype Y dwarf WISEP 182831.08+265037.8. The spectrum covers the 0.9-1.7 um wavelength range at a resolving power of lambda/Delta lambda ~180 and is a significant improvement over the previously published spectrum because it covers a broader wavelength range and is uncontaminated by light from a background star. The spectrum is unique for a cool brown dwarf in that the flux peaks in the Y, J, and H band are of near equal intensity in units of f_lambda. We fail to detect any absorption bands of NH_3 in the spectrum, in contrast to the predictions of chemical equilibrium models, but tentatively identify CH_4 as the carrier of an unknown absorption feature centered at 1.015 um. Using previously published ground- and spaced-based photometry, and using a Rayleigh Jeans tail to account for flux emerging longward of 4.5 um, we compute a bolometric luminosity of log (L_bol/L_sun)=-6.50+-0.02 which is significantly lower than previously published results. Finally, we compare the spectrum and photometry to two sets of atmospheric models and find that best overall match to the observed properties of WISEP 182831.08+265037.8 is a ~1 Gyr old binary composed of two T_eff~325 K, ~5 M_Jup brown dwarfs with subsolar [C/O] ratios. △ Less

Submitted 1 July, 2021; originally announced July 2021.

Comments: Accepted for publication in the Astrophysical Journal

Following the Lithium: Tracing Li-bearing Molecules Across Age, Mass, and Gravity in Brown Dwarfs

Authors: Ehsan Gharib-Nezhad, Mark S. Marley, Natasha E. Batalha, Channon Visscher, Richard S. Freedman, Roxana E. Lupu

Abstract: Lithium is an important element for the understanding of ultracool dwarfs because it is lost to fusion at masses above $\sim 68\, M_{\rm J}$. Hence, the presence or absence of atomic Li has served as an indicator of the nearby H-burning boundary at about $75\,M_{\rm J}$ between brown-dwarfs and very low-mass stars. Historically the "Lithium test", a search for the presence and strength of the Li l… ▽ More Lithium is an important element for the understanding of ultracool dwarfs because it is lost to fusion at masses above $\sim 68\, M_{\rm J}$. Hence, the presence or absence of atomic Li has served as an indicator of the nearby H-burning boundary at about $75\,M_{\rm J}$ between brown-dwarfs and very low-mass stars. Historically the "Lithium test", a search for the presence and strength of the Li line at 670.8 nm, has been a marker if an object has a substellar mass with stellar-like spectral energy distribution (e.g., a late-type M dwarf). While the Li test could in principle also be used to distinguish masses of later-type L-T dwarfs, Li is predominantly no longer found as an atomic gas, but rather a molecular species such as LiH, LiF, LiOH, and LiCl in their cooler atmospheres. L- and T-type brown dwarfs are also quite faint at 670 nm and thus challenging targets for high resolution spectroscopy. But only recently have experimental molecular line lists become available for the molecular Li species, allowing molecular Li mass discrimination. In this study, we generated the latest opacity of each of these Li-bearing molecules and performed thermochemical equilibrium atmospheric composition calculation of the abundance of these molecules. Finally, we computed thermal emission spectra for a series of radiative-convective equilibrium models of cloudy and cloudless brown dwarf atmospheres (with $T_{\rm eff}=$ 500--2400~K, and $\log g$=4.0, 4.5, 5.0) to understand where the presence or absence of atmospheric lithium-bearing species is most easily detected as a function of brown dwarf mass and age. After atomic Li, the best spectral signatures were found to be LiF at $10.5-12.5$~\micron and LiCl at $14.5-18.5$ $\micron$. LiH also shows a narrow feature at $\sim 9.38$ $\micron$. △ Less

Submitted 1 June, 2021; originally announced June 2021.

Comments: 13 pages, 7 figures, 1 table, In final revision at ApJ. Comments welcome

Variable Irradiation on 1D Cloudless Eccentric Exoplanet Atmospheres

Authors: L. C. Mayorga, Tyler D. Robinson, Mark S. Marley, E. M. May, Kevin B. Stevenson

Abstract: Exoplanets on eccentric orbits experience an incident stellar flux that can be markedly larger at periastron versus apoastron. This variation in instellation can lead to dramatic changes in atmospheric structure in regions of the atmosphere where the radiative and advective heating/cooling timescales are shorter than the orbital timescale. To explore this phenomenon, we develop a sophisticated one… ▽ More Exoplanets on eccentric orbits experience an incident stellar flux that can be markedly larger at periastron versus apoastron. This variation in instellation can lead to dramatic changes in atmospheric structure in regions of the atmosphere where the radiative and advective heating/cooling timescales are shorter than the orbital timescale. To explore this phenomenon, we develop a sophisticated one-dimensional (vertical) time-stepping atmospheric structure code, EGP+, capable of simulating the dynamic response of atmospheric thermal and chemical structure to time-dependent perturbations. Critically, EGP+ can efficiently simulate multiple orbits of a planet, thereby providing new opportunities for exoplanet modeling without the need for more computationally-expensive models. We make the simplifying assumption of cloud-free atmospheres, and apply our model to HAT-P-2b, HD~17156b, and HD~80606b, which are known to be on higher-eccentricity orbits. We find that for those planets which have Spitzer observations, our planet-to-star ratio predictions are roughly consistent with observations. However, we are unable to reproduce the observed peak offsets from periastron passage. Finally, we discuss promising pathways forward for adding new model complexity that would enable more detailed studies of clear and cloudy eccentric planets as well as worlds orbiting active host stars. △ Less

Submitted 17 May, 2021; originally announced May 2021.

Comments: 13 pages, 6 figures, 2 tables, accepted to ApJ

Impacts of Water Latent Heat on the Thermal Structure of Ultra-Cool Objects: Brown Dwarfs and Free-Floating Planets

Authors: Shih-Yun Tang, Tyler D. Robinson, Mark S. Marley, Natasha E. Batalha, Roxana Lupu, L. Prato

Abstract: Brown dwarfs are essential targets for understanding planetary and sub-stellar atmospheres across a wide range of thermal and chemical conditions. As surveys continue to probe ever deeper, and as observing capabilities continue to improve, the number of known Y dwarfs -- the coldest class of sub-stellar objects, with effective temperatures below about 600 K -- is rapidly growing. Critically, this… ▽ More Brown dwarfs are essential targets for understanding planetary and sub-stellar atmospheres across a wide range of thermal and chemical conditions. As surveys continue to probe ever deeper, and as observing capabilities continue to improve, the number of known Y dwarfs -- the coldest class of sub-stellar objects, with effective temperatures below about 600 K -- is rapidly growing. Critically, this class of ultra-cool objects has atmospheric conditions that overlap with Solar System worlds and, as a result, tools and ideas developed from studying Earth, Jupiter, Saturn and other nearby worlds are well-suited for application to sub-stellar atmospheres. To that end, we developed a one-dimensional (vertical) atmospheric structure model for ultra-cool objects that includes moist adiabatic convection, as this is an important process for many Solar System planets. Application of this model across a range of effective temperatures (350, 300, 250, 200 K), metallicities ([M/H] of 0.0, 0.5, 0.7, 1.5), and gravities (log $g$ of 4.0, 4.5, 4.7, 5.0) demonstrates strong impacts of water latent heat release on simulated temperature-pressure profiles. At the highest metallicities, water vapor mixing ratios reach an Earth-like 3%, with associated major alterations to the thermal structure in the atmospheric regions where water condenses. Spectroscopic and photometric signatures of metallicity and moist convection should be readily detectable at near- and mid-infrared wavelengths, especially with James Webb Space Telescope observations, and can help indicate the formation history of an object. △ Less

Submitted 16 August, 2021; v1 submitted 14 May, 2021; originally announced May 2021.

Comments: 16 pages, 10 figures, ApJ accepted version, data public available on zenodo

Uniform Forward-Modeling Analysis of Ultracool Dwarfs. II. Atmospheric Properties of 55 Late-T Dwarfs

Authors: Zhoujian Zhang, Michael C. Liu, Mark S. Marley, Michael R. Line, William M. J. Best

Abstract: We present a large forward-modeling analysis for 55 late-T (T7-T9) dwarfs, using low-resolution ($R\approx150$) near-infrared spectra and cloudless Sonora-Bobcat model atmospheres. We derive the objects' effective temperatures, surface gravities, metallicities, radii, masses, and luminosities using our newly developed Bayesian framework, and use the resulting population properties to test the mode… ▽ More We present a large forward-modeling analysis for 55 late-T (T7-T9) dwarfs, using low-resolution ($R\approx150$) near-infrared spectra and cloudless Sonora-Bobcat model atmospheres. We derive the objects' effective temperatures, surface gravities, metallicities, radii, masses, and luminosities using our newly developed Bayesian framework, and use the resulting population properties to test the model atmospheres. We find (1) our objects' fitted metallicities are 0.3-0.4 dex lower than those of nearby stars; (2) their ages derived from spectroscopic parameters are implausibly young; (3) their fitted temperatures show a similar spread as empirical temperature scales at a given spectral type but are $\sim100$ K hotter for $\geqslant$T8 dwarfs; and (4) their spectroscopically inferred masses are unphysically small. These results suggest the Sonora-Bobcat assumptions of cloudless and chemical-equilibrium atmospheres do not adequately reproduce late-T dwarf spectra. We also find a gravity- and a metallicity-dependence of temperatures. Combining the resulting parameter posteriors of our sample, we quantify the degeneracy between surface gravity and metallicity such that an increase in $Z$ combined with a $3.4\times$ increase in $\log{g}$ results in a spectrum that has similar fitted parameters. We note the systematic difference between our 1.0-2.5 $μ$m spectra and the Sonora-Bobcat models is $\approx$2-4% of the objects' peak $J$-band fluxes, implying modeling systematics will exceed measurement uncertainties when analyzing data with $J$-band S/N $\gtrsim50$. Using our large sample, we examine the fitting residuals as a function of wavelength and atmospheric properties to discern how to improve the models. Our work constitutes the largest analysis of brown dwarf spectra using multi-metallicity models and the most systematic examination of ultracool model atmospheres to date. △ Less

Submitted 11 May, 2021; originally announced May 2021.

Comments: Accepted by ApJ. 52 pages including 21 figures and 9 tables

Cloud busting: enstatite and quartz clouds in the atmosphere of 2M2224-0158

Authors: Ben Burningham, Jacqueline K. Faherty, Eileen C. Gonzales, Mark S. Marley, Channon Visscher, Roxana Lupu, Josefine Gaarn, Michelle Fabienne Bieger, Richard Freedman, Didier Saumon

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

Submitted 10 May, 2021; originally announced May 2021.

Comments: 20 pages, 10 figures, MNRAS accepted

EXOPLINES: Molecular Absorption Cross-Section Database for Brown Dwarf and Giant Exoplanet Atmospheres

Authors: Ehsan Gharib-Nezhad, Aishwarya R. Iyer, Michael R. Line, Richard S. Freedman, Mark S. Marley, Natasha E. Batalha

Abstract: Stellar, substellar, and planetary atmosphere models are all highly sensitive to the input opacities. Generational differences between various state-of-the-art stellar/planetary models are primarily because of incomplete and outdated atomic/molecular line-lists. Here we present a database of pre-computed absorption cross-sections for all isotopologues of key atmospheric molecules relevant to late-… ▽ More Stellar, substellar, and planetary atmosphere models are all highly sensitive to the input opacities. Generational differences between various state-of-the-art stellar/planetary models are primarily because of incomplete and outdated atomic/molecular line-lists. Here we present a database of pre-computed absorption cross-sections for all isotopologues of key atmospheric molecules relevant to late-type stellar, brown dwarf, and planetary atmospheres: MgH, AlH, CaH, TiH, CrH, FeH, SiO, TiO, VO, and H2O. The pressure and temperature ranges of the computed opacities are between 10$^{-6}$--3000~bar and 75--4000~K, and their spectral ranges are 0.25--330~$μ$m for many cases where possible. For cases with no pressure-broadening data, we use collision theory to bridge the gap. We also probe the effect of absorption cross-sections calculated from different line lists in the context of Ultra-Hot Jupiter and M-dwarf atmospheres. Using 1-D self-consistent radiative-convective thermochemical equilibrium models, we report significant variations in the theoretical spectra and thermal profiles of substellar atmospheres. With a 2000 K representative Ultra-Hot Jupiter, we report variations of up to 320 and 80 ppm in transmission and thermal emission spectra, respectively. For a 3000 K M-dwarf, we find differences of up to 125$\%$ in the spectra. We find that the most significant differences arise due to the choice of TiO line-lists, primarily below 1$μ$m. In sum, we present (1) a database of pre-computed molecular absorption cross-sections, and (2) quantify biases that arise when characterizing substellar/exoplanet atmospheres due to line list differences, therefore highlighting the importance of correct and complete opacities for eventual applications to high precision spectroscopy and photometry. △ Less

Submitted 1 April, 2021; originally announced April 2021.

Comments: accepted in ApJS. 44 pages, 17 figures, 7 tables. Comments are welcome

Weather on Other Worlds. V. The Three Most Rapidly Rotating Ultra-Cool Dwarfs

Authors: Megan E. Tannock, Stanimir Metchev, Aren Heinze, Paulo A. Miles-Páez, Jonathan Gagné, Adam Burgasser, Mark S. Marley, Dániel Apai, Genaro Suárez, Peter Plavchan

Abstract: We present the discovery of rapid photometric variability in three ultra-cool dwarfs from long-duration monitoring with the Spitzer Space Telescope. The T7, L3.5, and L8 dwarfs have the shortest photometric periods known to date: ${1.080}^{+0.004}_{-0.005}$ h, ${1.14}^{+0.03}_{-0.01}$ h, and ${1.23}^{+0.01}_{-0.01}$ h, respectively. We confirm the rapid rotation through moderate-resolution infrare… ▽ More We present the discovery of rapid photometric variability in three ultra-cool dwarfs from long-duration monitoring with the Spitzer Space Telescope. The T7, L3.5, and L8 dwarfs have the shortest photometric periods known to date: ${1.080}^{+0.004}_{-0.005}$ h, ${1.14}^{+0.03}_{-0.01}$ h, and ${1.23}^{+0.01}_{-0.01}$ h, respectively. We confirm the rapid rotation through moderate-resolution infrared spectroscopy, which reveals projected rotational velocities between 79 and 104 km s$^{-1}$. We compare the near-infrared spectra to photospheric models to determine the objects' fundamental parameters and radial velocities. We find that the equatorial rotational velocities for all three objects are $\gtrsim$100 km s$^{-1}$. The three L and T dwarfs reported here are the most rapidly spinning and likely the most oblate field ultra-cool dwarfs known to date. Correspondingly, all three are excellent candidates for seeking auroral radio emission and net optical/infrared polarization. As of this writing, 78 L-, T-, and Y-dwarf rotation periods have now been measured. The clustering of the shortest rotation periods near 1 h suggests that brown dwarfs are unlikely to spin much faster. △ Less

Submitted 25 March, 2021; v1 submitted 2 March, 2021; originally announced March 2021.

Comments: 35 pages, 13 figures, 6 tables, accepted for publication in the Astronomical Journal, typos corrected in updated version

Cloud Parameterizations and their Effect on Retrievals of Exoplanet Reflection Spectroscopy

Authors: Sagnick Mukherjee, Natasha E. Batalha, Mark S. Marley

Abstract: Future space-based direct imaging missions will perform low-resolution (R$<$100) optical (0.3-1~$μ$m) spectroscopy of planets, thus enabling reflected spectroscopy of cool giants. Reflected light spectroscopy is encoded with rich information about the scattering and absorbing properties of planet atmospheres. Given the diversity of clouds and hazes expected in exoplanets, it is imperative we solid… ▽ More Future space-based direct imaging missions will perform low-resolution (R$<$100) optical (0.3-1~$μ$m) spectroscopy of planets, thus enabling reflected spectroscopy of cool giants. Reflected light spectroscopy is encoded with rich information about the scattering and absorbing properties of planet atmospheres. Given the diversity of clouds and hazes expected in exoplanets, it is imperative we solidify the methodology to accurately and precisely retrieve these scattering and absorbing properties that are agnostic to cloud species. In particular, we focus on determining how different cloud parameterizations affect resultant inferences of both cloud and atmospheric composition. We simulate mock observations of the reflected spectra from three top priority direct imaging cool giant targets with different effective temperatures, ranging from 135 K to 533 K. We perform retrievals of cloud structure and molecular abundances on these three planets using four different parameterizations, each with increasing levels of cloud complexity. We find that the retrieved atmospheric and scattering properties strongly depend on the choice of cloud parameterization. For example, parameterizations that are too simplistic tend to overestimate the abundances. Overall, we are unable to retrieve precise/accurate gravity beyond $\pm$50\%. Lastly, we find that even low SNR=5, low R=40 reflected light spectroscopy gives cursory zeroth order insights into cloud deck position relative to molecular and Rayleigh optical depth level. △ Less

Submitted 10 February, 2021; originally announced February 2021.

Comments: 21 Pages, 18 Figures, Accepted for publication in The Astrophysical Journal (ApJ)

Keys of a Mission to Uranus or Neptune, the Closest Ice Giants

Authors: Tristan Guillot, Jonathan Fortney, Emily Rauscher, Mark S. Marley, Vivien Parmentier, Mike Line, Hannah Wakeford, Yohai Kaspi, Ravit Helled, Masahiro Ikoma, Heather Knutson, Kristen Menou, Diana Valencia, Daniele Durante, Shigeru Ida, Scott J. Bolton, Cheng Li, Kevin B. Stevenson, Jacob Bean, Nicolas B. Cowan, Mark D. Hofstadter, Ricardo Hueso, Jeremy Leconte, Liming Li, Christoph Mordasini , et al. (4 additional authors not shown)

Abstract: Uranus and Neptune are the archetypes of "ice giants", a class of planets that may be among the most common in the Galaxy. They hold the keys to understand the atmospheric dynamics and structure of planets with hydrogen atmospheres inside and outside the solar system; however, they are also the last unexplored planets of the Solar System. Their atmospheres are active and storms are believed to be… ▽ More Uranus and Neptune are the archetypes of "ice giants", a class of planets that may be among the most common in the Galaxy. They hold the keys to understand the atmospheric dynamics and structure of planets with hydrogen atmospheres inside and outside the solar system; however, they are also the last unexplored planets of the Solar System. Their atmospheres are active and storms are believed to be fueled by methane condensation which is both extremely abundant and occurs at low optical depth. This means that mapping temperature and methane abundance as a function of position and depth will inform us on how convection organizes in an atmosphere with no surface and condensates that are heavier than the surrounding air, a general feature of giant planets. Owing to the spatial and temporal variability of these atmospheres, an orbiter is required. A probe would provide a reference atmospheric profile to lift ambiguities inherent to remote observations. It would also measure the abundances of noble gases which can be used to reconstruct the history of planet formation in the Solar System. Finally, mapping the planets' gravity and magnetic fields will be essential to constrain their global composition, atmospheric dynamics, structure and evolution. An exploration of Uranus or Neptune will be essential to understand these planets and will also be key to constrain and analyze data obtained at Jupiter, Saturn, and for numerous exoplanets with hydrogen atmospheres. △ Less

Submitted 17 December, 2020; originally announced December 2020.

Comments: arXiv admin note: substantial text overlap with arXiv:1908.02092

Uniform Forward-Modeling Analysis of Ultracool Dwarfs. I. Methodology and Benchmarking

Authors: Zhoujian Zhang, Michael C. Liu, Mark S. Marley, Michael R. Line, William M. J. Best

Abstract: We present a forward-modeling framework using the Bayesian inference tool Starfish and cloudless Sonora-Bobcat model atmospheres to analyze low-resolution ($R\approx80-250$) near-infrared ($1.0-2.5$ $μ$m) spectra of T dwarfs. Our approach infers effective temperatures, surface gravities, metallicities, radii, and masses, and by accounting for uncertainties from model interpolation and correlated r… ▽ More We present a forward-modeling framework using the Bayesian inference tool Starfish and cloudless Sonora-Bobcat model atmospheres to analyze low-resolution ($R\approx80-250$) near-infrared ($1.0-2.5$ $μ$m) spectra of T dwarfs. Our approach infers effective temperatures, surface gravities, metallicities, radii, and masses, and by accounting for uncertainties from model interpolation and correlated residuals due to instrumental effects and modeling systematics, produces more realistic parameter posteriors than traditional ($χ^2$-based) spectral-fitting analyses. We validate our framework by fitting the model atmospheres themselves and finding negligible offsets between derived and input parameters. We apply our methodology to three well-known benchmark late-T dwarfs, HD 3651B, GJ 570D, and Ross 458C, using both solar and non-solar metallicity atmospheric models. We also derive these benchmarks' physical properties using their bolometric luminosities, their primary stars' ages and metallicities, and Sonora-Bobcat evolutionary models. Assuming the evolutionary-based parameters are more robust, we find our atmospheric-based, forward-modeling analysis produces two outcomes. For HD 3615B and GJ 570D, spectral fits provide accurate $T_{\rm eff}$ and $R$ but underestimated $\log{g}$ (by $\approx1.2$ dex) and $Z$ (by $\approx0.35$ dex), likely due to the systematics from modeling the potassium line profiles. For Ross 458C, spectral fits provide accurate $\log{g}$ and $Z$ but overestimated $T_{\rm eff}$ (by $\approx120$ K) and underestimated $R$ (by $\approx1.6\times$), likely because our model atmospheres lack clouds, reduced vertical temperature gradients, or disequilibrium processes. Finally, the spectroscopically inferred masses of these benchmarks are all considerably underestimated. △ Less

Submitted 11 May, 2021; v1 submitted 24 November, 2020; originally announced November 2020.

Comments: Accepted by ApJ. Minor updates since v1: (1) calculation of (very small) effect by different photometric calibrations (Sec 4.1.6), (2) expanded literature review about possible systematics in evolutionary models (Sec 6), (3) examination of the (very small) difference between Sonora vs ATMO models for luminosity calculations (footnote 10), (4) slight revision to calculation plotted in Figure 2

Into the UV: The Atmosphere of the Hot Jupiter HAT-P-41b Revealed

Authors: Nikole K. Lewis, Hannah R. Wakeford, Ryan J. MacDonald, Jayesh M. Goyal, David K. Sing, Joanna Barstow, Diana Powell, Tiffany Kataria, Ishan Mishra, Mark S. Marley, Natasha E. Batalha, Julie I. Moses, Peter Gao, Tom J. Wilson, Katy L. Chubb, Thomas Mikal-Evans, Nikolay Nikolov, Nor Pirzkal, Jessica J. Spake, Kevin B. Stevenson, Jeff Valenti, Xi Zhang

Abstract: For solar-system objects, ultraviolet spectroscopy has been critical in identifying sources for stratospheric heating and measuring the abundances of a variety of hydrocarbon and sulfur-bearing species, produced via photochemical mechanisms, as well as oxygen and ozone. To date, less than 20 exoplanets have been probed in this critical wavelength range (0.2-0.4 um). Here we use data from Hubble's… ▽ More For solar-system objects, ultraviolet spectroscopy has been critical in identifying sources for stratospheric heating and measuring the abundances of a variety of hydrocarbon and sulfur-bearing species, produced via photochemical mechanisms, as well as oxygen and ozone. To date, less than 20 exoplanets have been probed in this critical wavelength range (0.2-0.4 um). Here we use data from Hubble's newly implemented WFC3 UVIS G280 grism to probe the atmosphere of the hot Jupiter HAT-P-41b in the ultraviolet through optical in combination with observations at infrared wavelengths. We analyze and interpret HAT-P-41b's 0.2-5.0 um transmission spectrum using a broad range of methodologies including multiple treatments of data systematics as well as comparisons with atmospheric forward, cloud microphysical, and multiple atmospheric retrieval models. Although some analysis and interpretation methods favor the presence of clouds or potentially a combination of Na, VO, AlO, and CrH to explain the ultraviolet through optical portions of HAT-P-41b's transmission spectrum, we find that the presence of a significant H- opacity provides the most robust explanation. We obtain a constraint for the abundance of H-, log(H-) = -8.65 +/- 0.62 in HAT-P-41b's atmosphere, which is several orders of magnitude larger than predictions from equilibrium chemistry for a 1700 - 1950 K hot Jupiter. We show that a combination of photochemical and collisional processes on hot hydrogen-dominated exoplanets can readily supply the necessary amount of H- and suggest that such processes are at work in HAT-P-41b and many other hot Jupiter atmospheres. △ Less

Submitted 16 October, 2020; originally announced October 2020.

Comments: 17 pages, 7 figures. Published in ApJL (October 8th 2020)

Journal ref: ApJL (2020): 902, L19

NLTT5306B: an inflated, weakly irradiated brown dwarf

Authors: S. L. Casewell, J. Debes, I. P. Braker, M. C. Cushing, G. Mace, M. S. Marley, J. Davy Kirkpatrick

Abstract: We present Spitzer observations at 3.6 and 4.5 microns and a near-infrared IRTF SpeX spectrum of the irradiated brown dwarf NLTT5306B. We determine that the brown dwarf has a spectral type of L5 and is likely inflated, despite the low effective temperature of the white dwarf primary star. We calculate brightness temperatures in the Spitzer wavebands for both the model radius, and Roche Lobe radius… ▽ More We present Spitzer observations at 3.6 and 4.5 microns and a near-infrared IRTF SpeX spectrum of the irradiated brown dwarf NLTT5306B. We determine that the brown dwarf has a spectral type of L5 and is likely inflated, despite the low effective temperature of the white dwarf primary star. We calculate brightness temperatures in the Spitzer wavebands for both the model radius, and Roche Lobe radius of the brown dwarf, and conclude that there is very little day-night side temperature difference. We discuss various mechanisms by which NLTT5306B may be inflated, and determine that while low mass brown dwarfs (M<35 MJup) are easily inflated by irradiation from their host star, very few higher mass brown dwarfs are inflated. The higher mass brown dwarfs that are inflated may be inflated by magnetic interactions or may have thicker clouds. △ Less

Submitted 14 October, 2020; originally announced October 2020.

Comments: 7 pages, 4 figures. Accepted for publication in MNRAS

Beyond Equilibrium Temperature: How the Atmosphere/Interior Connection Affects the Onset of Methane, Ammonia, and Clouds in Warm Transiting Giant Planets

Authors: Jonathan J. Fortney, Channon Visscher, Mark S. Marley, Callie E. Hood, Michael R. Line, Daniel P. Thorngren, Richard S. Freedman, Roxana Lupu

Abstract: The atmospheric pressure-temperature profiles for transiting giant planets cross a range of chemical transitions. Here we show that the particular shape of these irradiated profiles for warm giant planets below 1300 K lead to striking differences in the behavior of non-equilibrium chemistry compared to brown dwarfs of similar temperatures. Our particular focus is H$_2$O, CO, CH$_4$, CO$_2$, and NH… ▽ More The atmospheric pressure-temperature profiles for transiting giant planets cross a range of chemical transitions. Here we show that the particular shape of these irradiated profiles for warm giant planets below 1300 K lead to striking differences in the behavior of non-equilibrium chemistry compared to brown dwarfs of similar temperatures. Our particular focus is H$_2$O, CO, CH$_4$, CO$_2$, and NH$_3$ in Jupiter- and Neptune-class planets. We show the cooling history of a planet, which depends most significantly on planetary mass and age, can have a dominant effect on abundances in the visible atmosphere, often swamping trends one might expect based on Teq alone. The onset of detectable CH$_4$ in spectra can be delayed to lower Teq for some planets compared to equilibrium, or pushed to higher Teq. The detectability of NH$_3$ is typically enhanced compared to equilibrium expectations, which is opposite to the brown dwarf case. We find that both CH$_4$ and NH$_3$ can become detectable at around the same Teq (at Teq values that vary with mass and metallicity) whereas these "onset" temperatures are widely spaced for brown dwarfs. We suggest observational strategies to search for atmospheric trends and stress that non-equilibrium chemistry and clouds can serve as probes of atmospheric physics. As examples of atmospheric complexity, we assess three Neptune-class planets GJ 436b, GJ 3470b, and WASP-107, all around Teq=700 K. Tidal heating due to eccentricity damping in all three planets heats the deep atmosphere by thousands of degrees, and may explain the absence of CH$_4$ in these cool atmospheres. Atmospheric abundances must be interpreted in the context of physical characteristics of the planet. △ Less

Submitted 30 October, 2020; v1 submitted 30 September, 2020; originally announced October 2020.

Comments: Accepted to AJ. No additional significant changes

Haze Formation in Warm H2-rich Exoplanet Atmospheres

Authors: Chao He, Sarah M. Horst, Nikole K. Lewis, Xinting Yu, Julianne I. Moses, Patricia McGuiggan, Mark S. Marley, Eliza M. -R. Kempton, Caroline V. Morley, Jeff A. Valenti, Veronique Vuitton

Abstract: New observing capabilities coming online over the next few years will provide opportunities for characterization of exoplanet atmospheres. However, clouds/hazes could be present in the atmospheres of many exoplanets, muting the amplitude of spectral features. We use laboratory simulations to explore photochemical haze formation in H2-rich exoplanet atmospheres at 800 K with metallicity either 100… ▽ More New observing capabilities coming online over the next few years will provide opportunities for characterization of exoplanet atmospheres. However, clouds/hazes could be present in the atmospheres of many exoplanets, muting the amplitude of spectral features. We use laboratory simulations to explore photochemical haze formation in H2-rich exoplanet atmospheres at 800 K with metallicity either 100 and 1000 times solar. We find that haze particles are produced in both simulated atmospheres with small particle size (20 to 140 nm) and relative low production rate (2.4 x 10-5 to 9.7 x 10-5 mg cm-3 h-1), but the particle size and production rate is dependent on the initial gas mixtures and the energy sources used in the simulation experiments. The gas phase mass spectra show that complex chemical processes happen in these atmospheres and generate new gas products that can further react to form larger molecules and solid haze particles. Two H2-rich atmospheres with similar C/O ratios (~0.5) yield different haze particles size, haze production rate, and gas products, suggesting both the elemental abundances and their bonding environments in an atmosphere can significantly affect the photochemistry. There is no methane (CH4) in our initial gas mixtures, although CH4 is often believed to be required to generate organic hazes. However, haze production rates from our experiments with different initial gas mixtures indicate that CH4 is neither required to generate organic hazes nor necessary to promote the organic haze formation. The variety and relative yield of the gas products indicate that CO and N2 enrich chemical reactions in H2-rich atmospheres. △ Less

Submitted 21 August, 2020; originally announced August 2020.

Comments: Accepted in The Planetary Science Journal; 21 pages, 5 figures, 1 table