Physical Properties of 29 sdB+dM Eclipsing Binaries in Zwicky Transient Facility

The development of large-scale time-domain surveys provides an opportunity to study the physical properties as well as the evolutionary scenario of B-type subdwarfs (sdB) and M-type dwarfs (dM). Here, we obtained 33 sdB+dM eclipsing binaries based on the Zwicky Transient Facility (ZTF) light curves and {\sl Gaia} early data release 3 (EDR3) parallaxes. By using the PHOEBE code for light curve analysis, we obtain probability distributions for parameters of 29 sdB+dM. $R_1$, $R_2$, and $i$ are well determined, and the average uncertainty of mass ratio $q$ is 0.08. Our parameters are in good agreement with previous works if a typical mass of sdB is assumed. Based on parameters of 29 sdB+dM, we find that both the mass ratio $q$ and the companion's radius $R_2$ decrease with the shortening of the orbital period. For the three sdB+dMs with orbital periods less than 0.075 days, their companions are all brown dwarfs. The masses and radii of the companions satisfy the mass--radius relation for low-mass stars and brown dwarfs. Companions with radii between $0.12R_\odot$ and $0.15R_\odot$ seem to be missing in the observations. As more short-period sdB+dM eclipsing binaries are discovered and classified in the future with ZTF and {\sl Gaia}, we will have more information to constrain the evolutionary ending of sdB+dM.

Let the Great World Spin: Revealing the Stormy, Turbulent Nature of Young Giant Exoplanet Analogs with the Spitzer Space Telescope

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

Photometric Properties of L and T Dwarf Binaries

Recent studies using volume-limited samples of brown dwarfs have revealed the photometric evolution of single L and T dwarfs and unresolved binaries as they age and cool, in particular as they transition from L to T spectral types. We demonstrate that the near-infrared photometric evolution of the resolved components of L and T dwarf binaries is consistent with that of single objects using a volume-limited sample. In addition, we provide supporting evidence that the L9 dwarf WISE J185101.83+593508.6 is an unresolved binary and identify three more candidate or confirmed unresolved binaries all from being overluminous for their spectral type. We calculate a new binary fraction of 13.1% ± 1.8% for L0–T8 dwarfs.

Spectroscopy of Candidate Members of the Sco-Cen Complex*

We present spectroscopy of 285 previously identified candidate members of populations in the Sco-Cen complex, primarily Ophiuchus, Upper Sco, and Lupus. The spectra are used to measure spectral types and diagnostics of youth. We find that 269 candidates exhibit signatures of youth in our spectra or previous data, which is consistent with their membership in Sco-Cen. We have constructed compilations of candidate members of Ophiuchus, Upper Sco, and Lupus that have spectral classifications and evidence of youth, which contain a total of 2274 objects. In addition, we have used spectra from previous studies to classify three sources in Ophiuchus that have been proposed to be protostellar brown dwarfs: ISO Oph 70, 200, and 203. We measure spectral types of early M from those data, which are earlier than expected for young brown dwarfs based on evolutionary models (≳M6.5) and instead are indicative of stellar masses (∼0.6 M ⊙).

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

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–1.7 μm phase curves demonstrate rotational modulations with semi-amplitudes of 5.27% ± 0.02% and 29.1% ± 0.1%; both can be fit well 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 respective day- and nightside spectra. Fitting the irradiated brown dwarf model grids to the dayside spectra require a filling factor of ∼50%, further supporting a hot spot dominating the dayside emission. The nightside spectrum of WD 0137-349B is fit reasonably well by non-irradiated substellar models, and the one of EPIC 21223521B can be approximated by a Planck function. We find strong spectral variations in the brown dwarfs’ day/night flux and brightness temperature contrasts, highlighting 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, suggesting 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.

Mapping the Pressure-dependent Day–Night Temperature Contrast of a Strongly Irradiated Atmosphere with HST Spectroscopic Phase Curve

Many brown dwarfs are on ultrashort-period and tidally locked orbits around white dwarf hosts. Because of these small orbital separations, the brown dwarfs are irradiated at levels similar to hot Jupiters. Yet, they are easier to observe than hot Jupiters because white dwarfs are fainter than main-sequence stars at near-infrared wavelengths. Irradiated brown dwarfs are, therefore, ideal hot Jupiter analogs for studying the atmospheric response under strong irradiation and fast rotation. We present the 1.1–1.67 μm spectroscopic phase curve of the irradiated brown dwarf (SDSS1411-B) in the SDSS J141126.20 + 200911.1 brown dwarf–white dwarf binary with the near-infrared G141 grism of the Hubble Space Telescope Wide Field Camera 3. SDSS1411-B is a 50M Jup brown dwarf with an irradiation temperature of 1300 K and has an orbital period of 2.02864 hr. Our best-fit model suggests a phase-curve amplitude of 1.4% and places an upper limit of 11° for the phase offset from the secondary eclipse. After fitting the white dwarf spectrum, we extract the phase-resolved brown dwarf emission spectra. We report a highly wavelength-dependent day–night spectral variation, with a water-band flux variation of about 360% ± 70% and a comparatively small J-band flux variation of 37% ± 2%. By combining the atmospheric modeling results and the day–night brightness temperature variations, we derive a pressure-dependent temperature contrast. We discuss the difference in the spectral features of SDSS1411-B and hot Jupiter WASP-43b, as well as the lower-than-predicted day–night temperature contrast of J4111-BD. Our study provides the high-precision observational constraints on the atmospheric structures of an irradiated brown dwarf at different orbital phases.

The Sonora Substellar Atmosphere Models. II. Cholla: A Grid of Cloud-free, Solar Metallicity Models in Chemical Disequilibrium for the JWST Era

Exoplanet and brown dwarf atmospheres commonly show signs of disequilibrium chemistry. In the James Webb Space Telescope (JWST) era, high-resolution spectra of directly imaged exoplanets will allow the characterization of their atmospheres in more detail, and allow systematic tests for the presence of chemical species that deviate from thermochemical equilibrium in these atmospheres. Constraining the presence of disequilibrium chemistry in these atmospheres as a function of parameters such as their effective temperature and surface gravity will allow us to place better constraints on the physics governing these atmospheres. This paper is part of a series of works presenting the Sonora grid of atmosphere models. In this paper, we present a grid of cloud-free, solar metallicity atmospheres for brown dwarfs and wide-separation giant planets with key molecular species such as CH4, H2O, CO, and NH3 in disequilibrium. Our grid covers atmospheres with T eff ∈ [500 K, 1300 K], log g ∈ [3.0, 5.5] (cgs) and an eddy diffusion parameter of log K zz = 2 , 4 and 7 (cgs). We study the effect of different parameters within the grid on the temperature and composition profiles of our atmospheres. We discuss their effect on the near-infrared colors of our model atmospheres and the detectability of CH4, H2O, CO, and NH3 using the JWST. We compare our models against existing MKO and Spitzer observations of brown dwarfs and verify the importance of disequilibrium chemistry for T dwarf atmospheres. Finally, we discuss how our models can help constrain the vertical structure and chemical composition of these atmospheres.

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

The detection of disk-integrated polarization from Luhman 16 A and B in the H band, and subsequent modeling, has been interpreted in the framework of zonal cloud bands on these bodies. Recently, Tan and Showman investigated the 3D atmospheric circulation and cloud structures of brown dwarfs with general circulation models (GCMs), and their simulations yielded 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 16 A and B, along with the 3D 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 objects with cloud particle sizes between 0.5 and 1 μm for Luhman 16 A and of 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, overpredict 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.

Improved Dynamical Masses for Six Brown Dwarf Companions Using Hipparcos and Gaia EDR3

We present comprehensive orbital analyses and dynamical masses for the substellar companions Gl 229 B, Gl 758 B, HD 13724 B, HD 19467 B, HD 33632 Ab, and HD 72946 B. Our dynamical fits incorporate radial velocities, relative astrometry, and, most importantly, calibrated Hipparcos-Gaia EDR3 accelerations. For HD 33632 A and HD 72946 we perform three-body fits that account for their outer stellar companions. We present new relative astrometry of Gl 229 B with Keck/NIRC2, extending its observed baseline to 25 yr. We obtain a <1% mass measurement of 71.4 ± 0.6 M Jup for the first T dwarf Gl 229 B and a 1.2% mass measurement of its host star (0.579 ± 0.007 M ⊙) that agrees with the high-mass end of the M-dwarf mass–luminosity relation. We perform a homogeneous analysis of the host stars’ ages and use them, along with the companions’ measured masses and luminosities, to test substellar evolutionary models. Gl 229 B is the most discrepant, as models predict that an object this massive cannot cool to such a low luminosity within a Hubble time, implying that it may be an unresolved binary. The other companions are generally consistent with models, except for HD 13724 B, which has a host star activity age 3.8σ older than its substellar cooling age. Examining our results in context with other mass–age–luminosity benchmarks, we find no trend with spectral type but instead note that younger or lower-mass brown dwarfs are overluminous compared to models, while older or higher-mass brown dwarfs are underluminous. The presented mass measurements for some companions are so precise that the stellar host ages, not the masses, limit the analysis.