Search

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

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 the 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 (S/N). Because of the strong correlation between the orbit-related parameters and the planetary radius, prior information on the orbital distance and planetary phase angle yield much tighter constraints on the planetary radius: from ${R}_{{\rm{p}}}={2.95}_{-1.95}^{+2.69}\,{R}_{\oplus }$ without prior knowledge, to ${R}_{{\rm{p}}}={1.01}_{-0.19}^{+0.33}\,{R}_{\oplus }$ with prior determination of the orbit for S/N = 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.14}^{+0.17}\,{R}_{\oplus }$ ) or atmospheric characterization. Also, prior knowledge of the mass alone does not yield a tight radius constraint ( ${R}_{{\rm{p}}}={1.64}_{-0.80}^{+1.29}\,{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 S/N is more beneficial than additional prior observations.

The unprecedented medium-resolution ( R _λ ∼ 1500–3500) near- and mid-infrared (1–18 μ m) spectrum provided by JWST for the young (140 ± 20 Myr) low-mass (12–20 M _Jup ) L–T transition (L7) companion VHS 1256 b gives access to a catalog of molecular absorptions. In this study, we present a comprehensive analysis of this data set utilizing a forward-modeling approach applying our Bayesian framework, ForMoSA . We explore five distinct atmospheric models to assess their performance in estimating key atmospheric parameters: T _eff , log( g ), [M/H], C/O, γ , f _sed , 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 VHS 1256 b, 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 T _eff 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 VHS 1256 b 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.

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 ∼9–10 mag at ≳ λ / 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.

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 to 5.5, appropriate for a broad range of objects with masses between 1 and 84 M _J . 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 A. S. Ackerman & M. S. Marley cloud model, including cloud parameter f _sed values from 1 to 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. For solar-metallicity evolution models including clouds in warm objects, we find a hydrogen-burning minimum mass of 70.2 M _J , close to empirical measurements; we find a deuterium-burning minimum mass of 12.05 M _J (50% of initial D burned). 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 on Zenodo (doi: https://doi.org/10.5281/zenodo.12735103 ).

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–900 K (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] ≤ −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 ± 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.

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 low-resolution spectrometer spectra covering the ∼1–12 μ m wavelength range at λ /Δ λ ≈ 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 Hubble Space Telescope, 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 ∼1000 to 350 K with a median uncertainty of ± 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.

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

We present an atmospheric retrieval analysis of the Y0 brown dwarf WISE J035934.06−540154.6 using the low-resolution 0.96–12 μ m James Webb Space Telescope (JWST) spectrum presented in Beiler et al. 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 three to five times more precise than previous work that used Hubble Space Telescope (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) _atm , and strengths of vertical mixing as encapsulated by the eddy diffusion coefficient K _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 _zz = 10 ^9 [cm ^2 s ^−1 ]. To test the sensitivity of these results to our five-knot spline thermal profile model, we performed a second retrieval using the Madhusudhan & Seager 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 five-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 JWST have to characterize the atmospheres of cool brown dwarfs.

Brown dwarfs (BDs) in ultra-short-period orbits around white dwarfs (WDs) 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 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–1.7 μ m phase curves show rotational modulations with semiamplitudes of 10.5% ± 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 re-creates 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 nightside, whereas the featureless dayside spectrum is likely dominated by H ^− opacity and a temperature inversion, much like the other highly irradiated BD EPIC 2122B.

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 _zz diffusion parameter. The role of mixing in altering the abundances of C-N-O-bearing molecules has mostly been explored for atmospheres with a solar composition. However, atmospheric metallicity and the C/O ratio also impact atmospheric chemistry. Therefore, we present the Sonora Elf Owl grid of self-consistent cloud-free 1D radiative-convective equilibrium model atmospheres for JWST observations, which includes a variation in K _zz across several orders of magnitude and also encompasses subsolar to supersolar metallicities and C/O ratios. We find that the impact of K _zz on the T ( P ) profile and spectra is a strong function of both T _eff and metallicity. For metal-poor objects, K _zz has large impacts on the atmosphere at significantly higher T _eff than in metal-rich atmospheres, where the impact of K _zz is seen to occur at lower T _eff . We identify significant spectral degeneracies between varying K _zz and metallicity in multiple wavelength windows, in particular, at 3–5 μ m. We use the Sonora Elf Owl atmospheric grid to fit the observed spectra of a sample of nine early to late T-type objects from T _eff = 550–1150 K. We find evidence for very inefficient vertical mixing in these objects, with inferred K _zz values lying in the range between ∼10 ^1 and 10 ^4 cm ^2 s ^−1 . Using self-consistent models, we find that this slow vertical mixing is due to the observations, which probe mixing in the deep detached radiative zone in these atmospheres.

Approximate methods for radiative transfer equations that are fast, reliable, and accurate 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. (hereafter, Toon89) outlined a two-stream method for computing reflected light and thermal spectra that 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 that 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 four-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 four-term method is twice as slow as our two-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.

The presence of aerosols is intimately linked to the global energy budget and the composition of a planet’s atmosphere. Their ability to reflect incoming light prevents energy from being deposited into the atmosphere, and they shape the spectra of exoplanets. We observed five near-infrared secondary eclipses of WASP-80b with the Wide Field Camera 3 (WFC3) aboard the Hubble Space Telescope to provide constraints on the presence and properties of atmospheric aerosols. We detect a broadband eclipse depth of 34 ± 10 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.

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

We present JWST Early Release Science coronagraphic observations of the super-Jupiter exoplanet, HIP 65426b, with the Near-Infrared Camera (NIRCam) from 2 to 5 μ m, and with the Mid-Infrared Instrument (MIRI) from 11 to 16 μ m. At a separation of ∼0.″82 (87 ${}_{-31}^{+108}$ au), HIP 65426b 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 ∼1 × 10 ^−5 and ∼2 × 10 ^−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 _Jup beyond separations of ∼100 au. Together with existing ground-based near-infrared data, the JWST photometry are fit well by a BT-SETTL atmospheric model from 1 to 16 μ m, and they span ∼97% of HIP 65426b'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 and −4.14, which in turn provides a robust mass constraint of 7.1 ± 1.2 M _Jup . In totality, these observations confirm that JWST presents a powerful and exciting opportunity to characterize the population of exoplanets amenable to high-contrast imaging in greater detail.

We present the highest fidelity spectrum to date of a planetary-mass object. VHS 1256 b is a

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 ± 0.31 dex). The planet’s dynamical mass ( ${2.8}_{-0.5}^{+0.6}$ M _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\mathrm{ln}T/d\mathrm{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 _eff ≈ 800 K, $\mathrm{log}(g)\approx 3.7$ dex, and the presence of silicate clouds and disequilibrium chemistry in the atmosphere. Straddling the L/T transition, AF Lep b is thus far the coldest exoplanet with suggested evidence of silicate clouds.

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. 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 two to four, 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 four-term method (SH4) against Toon et al. and two independent and higher-fidelity methods ( CDISORT and doubling method). On average, the SH4 method provides an order-of-magnitude increase in accuracy, compared to Toon et al. Finally, we implement SH4 within PICASO and observe only a modest increase in computational time, compared to two-stream methods (20% increase).

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 ± 3 M _Jup brown dwarf, orbiting a T _eff = 7756 ± 35 K white dwarf with an ultra-short period of ∼102 minutes. We present Hubble Space Telescope/Wide Field Camera 3 spectroscopic phase curves of NLTT5306, consisting of 47 spectra from 1.1 to 1.7 μ m with an average signal-to-noise ratio ∼ 65 per wavelength. We extracted the phase-resolved spectra of the brown dwarf NLTT5306B, finding a small

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 Gemini Planet Imager monitoring of the companion’s orbit with six 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 inner radii that are close to the companion’s semimajor axis. Including a posteriori knowledge of the system’s SED further narrows the constraints on the disk’s inner radius and places 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 test bed for dynamical interaction between a debris disk and a substellar companion.

We use observations with the infrared-optimized Magellan Adaptive Optics (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 a 180° dark region (gvAPP-180). To remove residual starlight in postprocessing, we apply a time-domain principal-components-analysis-based algorithm we call PCA-Temporal, which uses eigen time series rather than eigenimages to subtract starlight. By casting the problem in terms of eigen time series, we reduce the computational cost of postprocessing the data, enabling the use of the fully sampled data set 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 postprocessed contrast limits of 1.5 × 10 ^−6 –9.8 × 10 ^−6 outside of 0.″75, which correspond to planet masses of 2.6–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 ≥9 M _J planets from 2.5 to 7 au.

Upcoming James Webb Space Telescope 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 1D and 3D 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 1Dl atmospheric radiative-convective equilibrium (RCE) model. This model has heritage from the Fortran-based code, 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 RCE 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).