Engineering topics

Malik, Matej

Publications and source records attributed to Malik, Matej.

7 records

DOE OSTI2021

Spitzer Phase-curve Observations and Circulation Models of the Inflated Ultrahot Jupiter WASP-76b

The large radii of many hot Jupiters can only be matched by models that have hot interior adiabats, and recent theoretical work has shown that the interior evolution of hot Jupiters has a significant impact on their atmospheric structure. Due to its inflated radius, low gravity, and ultrahot equilibrium temperature, WASP-76b is an ideal case study for the impact of internal evolution on observable properties. Hot interiors should most strongly affect the nonirradiated side of the planet, and thus full phase-curve observations are critical to ascertain the effect of the interior on the atmospheres of hot Jupiters. In this work, we present the first Spitzer phase-curve observations of WASP-76b. We find that WASP-76b has an ultrahot dayside and relatively cold nightside with brightness temperatures of 2471 ± 27 K/1518 ± 61 K at 3.6 μm and 2699 ± 32 K/1259 ± 44 K at 4.5 μm, respectively. These results provide evidence for a dayside thermal inversion. Both channels exhibit small phase offsets of 0.68 ± 0.°48 at 3.6 μm and 0.67 ± 0.°2 at 4.5 μm. We compare our observations to a suite of general circulation models (GCMs) that consider two endmembers of interior temperature along with a broad range of frictional drag strengths. Strong frictional drag is necessary to match the small phase offsets and cold nightside temperatures observed. From our suite of cloud-free GCMs, we find that only cases with a cold interior can reproduce the cold nightsides and large phase-curve amplitude at 4.5 μm, hinting that the hot interior adiabat of WASP-76b does not significantly impact its atmospheric dynamics or that clouds blanket its nightside.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Optical Transmission Spectroscopy of the Terrestrial Exoplanet LHS 3844b from 13 Ground-based Transit Observations

Atmospheric studies of spectroscopically accessible terrestrial exoplanets lay the groundwork for comparative planetology between these worlds and the solar system terrestrial planets. LHS 3844b is a highly irradiated terrestrial exoplanet (R = 1.303 ± 0.022R {sub ⊕}) orbiting a mid-M dwarf 15 parsecs away. Work based on near-infrared Spitzer phase curves ruled out atmospheres with surface pressures ≥10 bars on this planet. We present 13 transit observations of LHS 3844b taken with the Magellan Clay telescope and the LDSS3C multi-object spectrograph covering 620–1020 nm. We analyze each of the 13 data sets individually using a Gaussian process regression, and present both white and spectroscopic light curves. In the combined white light curve we achieve an rms precision of 65 ppm when binning to 10 minutes. The mean white light-curve value of (R {sub p}/R {sub s}){sup 2} is 0.4170 ± 0.0046%. To construct the transmission spectrum, we split the white light curves into 20 spectrophotometric bands, each spanning 20 nm, and compute the mean values of (R {sub p}/R {sub s}){sup 2} in each band. We compare the transmission spectrum to two sets of atmospheric models. We disfavor a clear, solar composition atmosphere (μ = 2.34) with surface pressures ≥0.1 bar to 5.2σ confidence. We disfavor a clear, H{sub 2}O steam atmosphere (μ = 18) with surface pressures ≥0.1 bar to low confidence (2.9σ). Our observed transmission spectrum favors a flat line. For solar composition atmospheres with surface pressures ≥1 bar we rule out clouds with cloud-top pressures of 0.1 bar (5.3σ), but we cannot address high-altitude clouds at lower pressures. Our results add further evidence that LHS 3844b is devoid of an atmosphere.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

Self-luminous and Irradiated Exoplanetary Atmospheres Explored with HELIOS

We present new methodological features and physical ingredients included in the one-dimensional radiative transfer code HELIOS, improving the hemispheric two-stream formalism. We conduct a thorough intercomparison survey with several established forward models, including COOLTLUSTY and PHOENIX, and find satisfactory consistency with their results. Then, we explore the impact of (i) different groups of opacity sources, (ii) a stellar path length adjustment, and (iii) a scattering correction on self-consistently calculated atmospheric temperatures and planetary emission spectra. First, we observe that temperature–pressure (T–P) profiles are very sensitive to the opacities included, with metal oxides, hydrides, and alkali atoms (and ionized hydrogen) playing an important role in the absorption of shortwave radiation (in very hot surroundings). Moreover, if these species are sufficiently abundant, they are likely to induce nonmonotonic T–P profiles. Second, without the stellar path length adjustment, the incoming stellar flux is significantly underestimated for zenith angles above 80°, which somewhat affects the upper atmospheric temperatures and the planetary emission. Third, the scattering correction improves the accuracy of the computation of the reflected stellar light by ∼10%. We use HELIOS to calculate a grid of cloud-free atmospheres in radiative–convective equilibrium for self-luminous planets for a range of effective temperatures, surface gravities, metallicities, and C/O ratios to be used by planetary evolution studies. Furthermore, we calculate dayside temperatures and secondary eclipse spectra for a sample of exoplanets for varying chemistry and heat redistribution. These results may be used to make predictions on the feasibility of atmospheric characterizations with future observations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2018

Revisiting the Phase Curves of WASP-43b: Confronting Re-analyzed Spitzer Data with Cloudy Atmospheres

Recently acquired Hubble and Spitzer phase curves of the short-period hot Jupiter WASP-43b make it an ideal target for confronting theory with data. On the observational front, we re-analyze the 3.6 and 4.5 μm Spitzer phase curves and demonstrate that our improved analysis better removes residual red noise due to intra-pixel sensitivity, which leads to greater fluxes emanating from the nightside of WASP-43b, thus reducing the tension between theory and data. On the theoretical front, we construct cloud-free and cloudy atmospheres of WASP-43b using our Global Circulation Model (GCM), THOR, which solves the non-hydrostatic Euler equations (compared to GCMs that typically solve the hydrostatic primitive equations). The cloud-free atmosphere produces a reasonable fit to the dayside emission spectrum. The multi-phase emission spectra constrain the cloud deck to be confined to the nightside and have a finite cloud-top pressure. The multi-wavelength phase curves are naturally consistent with our cloudy atmospheres, except for the 4.5 μm phase curve, which requires the presence of enhanced carbon dioxide in the atmosphere of WASP-43b. Multi-phase emission spectra at higher spectral resolution, as may be obtained using the James Webb Space Telescope, and a reflected-light phase curve at visible wavelengths would further constrain the properties of clouds in WASP-43b.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Retrieval Analysis of the Emission Spectrum of WASP-12b: Sensitivity of Outcomes to Prior Assumptions and Implications for Formation History

We analyze the emission spectrum of the hot Jupiter WASP-12b using our HELIOS-R retrieval code and HELIOS-K opacity calculator. When interpreting Hubble and Spitzer data, the retrieval outcomes are found to be prior-dominated. When the prior distributions of the molecular abundances are assumed to be log-uniform, the volume mixing ratio of HCN is found to be implausibly high. A VULCAN chemical kinetics model of WASP-12b suggests that chemical equilibrium is a reasonable assumption even when atmospheric mixing is implausibly rigorous. Guided by (exo)planet formation theory, we set Gaussian priors on the elemental abundances of carbon, oxygen, and nitrogen with the Gaussian peaks being centered on the measured C/H, O/H, and N/H values of the star. By enforcing chemical equilibrium, we find substellar O/H and stellar to slightly superstellar C/H for the dayside atmosphere of WASP-12b. The superstellar carbon-to-oxygen ratio is just above unity, regardless of whether clouds are included in the retrieval analysis, consistent with Madhusudhan et al. Furthermore, whether a temperature inversion exists in the atmosphere depends on one’s assumption for the Gaussian width of the priors. Our retrieved posterior distributions are consistent with the formation of WASP-12b in a solar-composition protoplanetary disk, beyond the water iceline, via gravitational instability or pebble accretion (without core erosion) and migration inward to its present orbital location via a disk-free mechanism, and are inconsistent with both in situ formation and core accretion with disk migration, as predicted by Madhusudhan et al. We predict that the interpretation of James Webb Space Telescope WASP-12b data will not be prior-dominated.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

HELIOS–RETRIEVAL: An Open-source, Nested Sampling Atmospheric Retrieval Code; Application to the HR 8799 Exoplanets and Inferred Constraints for Planet Formation

We present an open-source retrieval code named HELIOS–RETRIEVAL, designed to obtain chemical abundances and temperature–pressure profiles by inverting the measured spectra of exoplanetary atmospheres. In our forward model, we use an exact solution of the radiative transfer equation, in the pure absorption limit, which allows us to analytically integrate over all of the outgoing rays. Two chemistry models are considered: unconstrained chemistry and equilibrium chemistry (enforced via analytical formulae). The nested sampling algorithm allows us to formally implement Occam’s Razor based on a comparison of the Bayesian evidence between models. We perform a retrieval analysis on the measured spectra of the four HR 8799 directly imaged exoplanets. Chemical equilibrium is disfavored for HR 8799b and c. We find supersolar C/H and O/H values for the outer HR 8799b and c exoplanets, while the inner HR 8799d and e exoplanets have a range of C/H and O/H values. The C/O values range from being superstellar for HR 8799b to being consistent with stellar for HR 8799c and being substellar for HR 8799d and e. If these retrieved properties are representative of the bulk compositions of the exoplanets, then they are inconsistent with formation via gravitational instability (without late-time accretion) and consistent with a core accretion scenario in which late-time accretion of ices occurred differently for the inner and outer exoplanets. For HR 8799e, we find that spectroscopy in the K band is crucial for constraining C/O and C/H. HELIOS–RETRIEVAL is publicly available as part of the Exoclimes Simulation Platform (http://www.exoclime.org).

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

HELIOS: AN OPEN-SOURCE, GPU-ACCELERATED RADIATIVE TRANSFER CODE FOR SELF-CONSISTENT EXOPLANETARY ATMOSPHERES

We present the open-source radiative transfer code named HELIOS, which is constructed for studying exoplanetary atmospheres. In its initial version, the model atmospheres of HELIOS are one-dimensional and plane-parallel, and the equation of radiative transfer is solved in the two-stream approximation with nonisotropic scattering. A small set of the main infrared absorbers is employed, computed with the opacity calculator HELIOS-K and combined using a correlated-k approximation. The molecular abundances originate from validated analytical formulae for equilibrium chemistry. We compare HELIOS with the work of Miller-Ricci and Fortney using a model of GJ 1214b, and perform several tests, where we find: model atmospheres with single-temperature layers struggle to converge to radiative equilibrium; k-distribution tables constructed with ≳0.01 cm{sup −1} resolution in the opacity function (≲10{sup 3} points per wavenumber bin) may result in errors ≳1%–10% in the synthetic spectra; and a diffusivity factor of 2 approximates well the exact radiative transfer solution in the limit of pure absorption. We construct “null-hypothesis” models (chemical equilibrium, radiative equilibrium, and solar elemental abundances) for six hot Jupiters. We find that the dayside emission spectra of HD 189733b and WASP-43b are consistent with the null hypothesis, while the latter consistently underpredicts the observed fluxes of WASP-8b, WASP-12b, WASP-14b, and WASP-33b. We demonstrate that our results are somewhat insensitive to the choice of stellar models (blackbody, Kurucz, or PHOENIX) and metallicity, but are strongly affected by higher carbon-to-oxygen ratios. The code is publicly available as part of the Exoclimes Simulation Platform (exoclime.net).

79 ASTRONOMY AND ASTROPHYSICS↗