skip to main content

Title: Modeling Polarization Signals from Cloudy Brown Dwarfs Luhman 16 A and B in Three Dimensions
Abstract 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.  more » « less
Award ID(s):
1828315 2103241
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
The Astrophysical Journal
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The photometric and spectral variability of brown dwarfs probes heterogeneous temperature and cloud distributions and traces the atmospheric circulation patterns. We present a new 42 hr Hubble Space Telescope (HST) Wide Field Camera 3 G141 spectral time series of VHS 1256-1257 b, a late L-type planetary-mass companion that has been shown to have one of the highest variability amplitudes among substellar objects. The light curve is rapidly evolving and best fit by a combination of three sine waves with different periods and a linear trend. The amplitudes of the sine waves and the linear slope vary with the wavelength, and the corresponding spectral variability patterns match the predictions by models invoking either heterogeneous clouds or thermal profile anomalies. Combining these observations with previous HST monitoring data, we find that the peak-to-valley flux difference is 33% ± 2% with an even higher amplitude reaching 38% in the J band, the highest amplitude ever observed in a substellar object. The observed light curve can be explained by maps that are composed of zonal waves, spots, or a mixture of the two. Distinguishing the origin of rapid light curve evolution requires additional long-term monitoring. Our findings underscore the essential role of atmospheric dynamics in shaping brown-dwarf atmospheres and highlight VHS 1256-1257 b as one of the most favorable targets for studying the atmospheres, clouds, and atmospheric circulation of planets and brown dwarfs. 
    more » « less
  2. null (Ed.)
    Context. Young giant planets and brown dwarf companions emit near-infrared radiation that can be linearly polarized up to several percent. This polarization can reveal the presence of an (unresolved) circumsubstellar accretion disk, rotation-induced oblateness of the atmosphere, or an inhomogeneous distribution of atmospheric dust clouds. Aims. We aim to measure the near-infrared linear polarization of 20 known directly imaged exoplanets and brown dwarf companions. Methods. We observed the companions with the high-contrast imaging polarimeter SPHERE-IRDIS at the Very Large Telescope. We reduced the data using the IRDAP pipeline to correct for the instrumental polarization and crosstalk of the optical system with an absolute polarimetric accuracy <0.1% in the degree of polarization. We employed aperture photometry, angular differential imaging, and point-spread-function fitting to retrieve the polarization of the companions. Results. We report the first detection of polarization originating from substellar companions, with a polarization of several tenths of a percent for DH Tau B and GSC 6214-210 B in H -band. By comparing the measured polarization with that of nearby stars, we find that the polarization is unlikely to be caused by interstellar dust. Because the companions have previously measured hydrogen emission lines and red colors, the polarization most likely originates from circumsubstellar disks. Through radiative transfer modeling, we constrain the position angles of the disks and find that the disks must have high inclinations. For the 18 other companions, we do not detect significant polarization and place subpercent upper limits on their degree of polarization. We also present images of the circumstellar disks of DH Tau, GQ Lup, PDS 70, β Pic, and HD 106906. We detect a highly asymmetric disk around GQ Lup and find evidence for multiple scattering in the disk of PDS 70. Both disks show spiral-like features that are potentially induced by GQ Lup B and PDS 70 b, respectively. Conclusions. The presence of the disks around DH Tau B and GSC 6214-210 B as well as the misalignment of the disk of DH Tau B with the disk around its primary star suggest in situ formation of the companions. The non-detections of polarization for the other companions may indicate the absence of circumsubstellar disks, a slow rotation rate of young companions, the upper atmospheres containing primarily submicron-sized dust grains, and/or limited cloud inhomogeneity. 
    more » « less
  3. Abstract Direct imaging studies have mainly used low-resolution spectroscopy ( R ∼ 20–100) to study the atmospheres of giant exoplanets and brown dwarf companions, but the presence of clouds has often led to degeneracies in the retrieved atmospheric abundances (e.g., carbon-to-oxygen ratio, metallicity). This precludes clear insights into the formation mechanisms of these companions. The Keck Planet Imager and Characterizer (KPIC) uses adaptive optics and single-mode fibers to transport light into NIRSPEC ( R ∼ 35,000 in the K band), and aims to address these challenges with high-resolution spectroscopy. Using an atmospheric retrieval framework based on petitRADTRANS , we analyze the KPIC high-resolution spectrum (2.29–2.49 μ m) and the archival low-resolution spectrum (1–2.2 μ m) of the benchmark brown dwarf HD 4747 B ( m = 67.2 ± 1.8 M Jup , a = 10.0 ± 0.2 au, T eff ≈ 1400 K). We find that our measured C/O and metallicity for the companion from the KPIC high-resolution spectrum agree with those of its host star within 1 σ –2 σ . The retrieved parameters from the K -band high-resolution spectrum are also independent of our choice of cloud model. In contrast, the retrieved parameters from the low-resolution spectrum are highly sensitive to our chosen cloud model. Finally, we detect CO, H 2 O, and CH 4 (volume-mixing ratio of log(CH 4 ) = −4.82 ± 0.23) in this L/T transition companion with the KPIC data. The relative molecular abundances allow us to constrain the degree of chemical disequilibrium in the atmosphere of HD 4747 B, and infer a vertical diffusion coefficient that is at the upper limit predicted from mixing length theory. 
    more » « less
  4. Abstract Extremely large telescopes (ELTs) provide an opportunity to observe surface inhomogeneities for ultracool objects including M dwarfs, brown dwarfs (BDs), and gas giant planets via Doppler imaging and spectrophotometry techniques. These inhomogeneities can be caused by star spots, clouds, and vortices. Star spots and associated stellar flares play a significant role in habitability, either stifling life or catalyzing abiogenesis depending on the emission frequency, magnitude, and orientation. Clouds and vortices may be the source of spectral and photometric variability observed at the L/T transition of BDs and are expected in gas giant exoplanets. We develop a versatile analytical framework to model and infer surface inhomogeneities that can be applied to both spectroscopic and photometric data. This model is validated against a slew of numerical simulations. Using archival spectroscopic and photometric data, we infer starspot parameters (location, size, and contrast) and generate global surface maps for Luhman 16B (an early T dwarf and one of our solar system’s nearest neighbors at a distance of ≈2 pc). We confirm previous findings that Luhman 16B’s atmosphere is inhomogeneous with time-varying features. In addition, we provide tentative evidence of longer timescale atmospheric structures such as dark equatorial and bright midlatitude to polar spots. These findings are discussed in the context of atmospheric circulation and dynamics for ultracool dwarfs. Our analytical model will be valuable in assessing the feasibility of using ELTs to study surface inhomogeneities of gas giant exoplanets and other ultracool objects. 
    more » « less
  5. Abstract

    The geographical variability, frequency content, and vertical structure of near‐surface oceanic kinetic energy (KE) are important for air‐sea interaction, marine ecosystems, operational oceanography, pollutant tracking, and interpreting remotely sensed velocity measurements. Here, KE in high‐resolution global simulations (HYbrid Coordinate Ocean Model; HYCOM, and Massachusetts Institute of Technology general circulation model; MITgcm), at the sea surface (0 m) and at 15 m, are compared with KE from undrogued and drogued surface drifters, respectively. Global maps and zonal averages are computed for low‐frequency (<0.5 cpd), near‐inertial, diurnal, and semidiurnal bands. Both models exhibit low‐frequency equatorial KE that is low relative to drifter values. HYCOM near‐inertial KE is higher than in MITgcm, and closer to drifter values, probably due to more frequently updated atmospheric forcing. HYCOM semidiurnal KE is lower than in MITgcm, and closer to drifter values, likely due to inclusion of a parameterized topographic internal wave drag. A concurrent tidal harmonic analysis in the diurnal band demonstrates that much of the diurnal flow is nontidal. We compute simple proxies of near‐surface vertical structure—the ratio 0 m KE/(0 m KE + 15 m KE) in model outputs, and the ratio undrogued KE/(undrogued KE + drogued KE) in drifter observations. Over most latitudes and frequency bands, model ratios track the drifter ratios to within error bars. Values of this ratio demonstrate significant vertical structure in all frequency bands except the semidiurnal band. Latitudinal dependence in the ratio is greatest in diurnal and low‐frequency bands.

    more » « less