Researcher profile

Matej Malik

Matej Malik contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

10 published item(s)

preprint2022arXiv

Confirmation of Water Absorption in the Thermal Emission Spectrum of the Hot Jupiter WASP-77Ab with HST/WFC3

Secondary eclipse observations of hot Jupiters can reveal both their compositions and thermal structures. Previous observations have shown a diversity of hot Jupiter eclipse spectra, including absorption features, emission features, and featureless blackbody-like spectra. We present a secondary eclipse spectrum of the hot Jupiter WASP-77Ab observed between $1-5$ $μ$m with the Hubble Space Telescope (HST) and the Spitzer Space Telescope. The HST observations show signs of water absorption indicative of a non-inverted thermal structure. We fit the data with both a one-dimensional free retrieval and a grid of one-dimensional self-consistent forward models to confirm this non-inverted structure. The free retrieval places a $3σ$ lower limit on the atmospheric water abundance of $\log(n_\mathrm{H_2O})>-4.78$ and can not constrain the CO abundance. The grid fit produces a slightly super-stellar metallicity and constrains the carbon-to-oxygen ratio to less than or equal to the solar value. We also compare our data to recent high-resolution observations of WASP-77Ab taken with the Gemini-South/IGRINS spectrograph and find that our observations are consistent with the best-fit model to the high-resolution data. However, the metallicity derived from the IGRINS data is significantly lower than that derived from our self-consistent model fit. We find that this difference may be due to disequilibrium chemistry, and the varying results between the models applied here demonstrate the difficulty of constraining disequilibrium chemistry with low-resolution, low wavelength coverage data alone. Future work to combine observations from IGRINS, HST, and JWST will improve our estimate of the atmospheric composition of WASP-77Ab.

preprint2022arXiv

Strong H$_2$O and CO emission features in the spectrum of KELT-20b driven by stellar UV irradiation

Know thy star, know thy planetary atmosphere. Every exoplanet with atmospheric measurements orbits around a star, and the stellar environment directly affects the planetary atmosphere. Here we present the emission spectrum of ultra-hot Jupiter KELT-20b which provides an observational link between host star properties and planet atmospheric thermal structure. It is currently the only planet with thermal emission measurements in the $T_{eq}\sim$2200K range that orbits around an early A-type star. By comparing it with other similar ultra-hot Jupiters around FGK stars, we can better understand how different host star types influence planetary atmospheres. The emission spectrum covers 0.6 to 4.5 $μm$ with data from TESS, HST WFC3/G141, and Spitzer 4.5 $μm$ channel. KELT-20b has a 1.4 $μm$ water feature strength metric of S$_{H_2O}$ = -0.097$\pm$0.02 and a blackbody brightness temperature difference of 528K between WFC3/G141 (T$_b$=2402$\pm$14K) and Spitzer 4.5 $μm$ channel (T$_b$=2930$\pm59$K). These very large H$_2$O and CO emission features combined with the A-type host star make KELT-20b a unique planet among other similar hot Jupiters. The abundant FUV, NUV, and optical radiation from its host star (T$_{eff}=8720\pm250$K) is expected to be the key that drives its strong thermal inversion and prominent emission features based on previous PHOENIX models calculations.

preprint2022arXiv

The THOR+HELIOS general circulation model: multi-wavelength radiative transfer with accurate scattering by clouds/hazes

General circulation models (GCMs) provide context for interpreting multi-wavelength, multi-phase data of the atmospheres of tidally locked exoplanets. In the current study, the non-hydrostatic THOR GCM is coupled with the HELIOS radiative transfer solver for the first time, supported by an equilibrium chemistry solver (FastChem), opacity calculator (HELIOS-K) and Mie scattering code (LX-MIE). To accurately treat the scattering of radiation by medium-sized to large aerosols/condensates, improved two-stream radiative transfer is implemented within a GCM for the first time. Multiple scattering is implemented using a Thomas algorithm formulation of the two-stream flux solutions, which decreases the computational time by about 2 orders of magnitude compared to the iterative method used in past versions of HELIOS. As a case study, we present four GCMs of the hot Jupiter WASP-43b, where we compare the temperature, velocity, entropy, and streamfunction, as well as the synthetic spectra and phase curves, of runs using regular versus improved two-stream radiative transfer and isothermal versus non-isothermal layers. While the global climate is qualitatively robust, the synthetic spectra and phase curves are sensitive to these details. A THOR+HELIOS WASP-43b GCM (horizontal resolution of about 4 degrees on the sphere and with 40 radial points) with multi-wavelength radiative transfer (30 k-table bins) running for 3000 Earth days (864,000 time steps) takes about 19-26 days to complete depending on the type of GPU.

preprint2021arXiv

A Comparative Study of Atmospheric Chemistry with VULCAN

We present an update of the open-source photochemical kinetics code VULCAN (Tsai et al. 2017; https://github.com/exoclime/VULCAN) to include C-H-N-O-S networks and photochemistry. Additional new features are advection transport, condensation, various boundary conditions, and temperature-dependent UV cross-sections. First, we validate our photochemical model for hot Jupiter atmospheres by performing an intercomparison of HD 189733b models between Moses et al. (2011), Venot et al. (2012), and VULCAN, to diagnose possible sources of discrepancy. Second, we set up a model of Jupiter extending from the deep troposphere to upper stratosphere to verify the kinetics for low temperature. Our model reproduces hydrocarbons consistent with observations, and the condensation scheme successfully predicts the locations of water and ammonia ice clouds. We show that vertical advection can regulate the local ammonia distribution in the deep atmosphere. Third, we validate the model for oxidizing atmospheres by simulating Earth and find agreement with observations. Last, VULCAN is applied to four representative cases of extrasolar giant planets: WASP-33b, HD 189733b, GJ 436b, and 51 Eridani b. We look into the effects of the C/O ratio and chemistry of titanium/vanadium species for WASP-33b; we revisit HD 189733b for the effects of sulfur and carbon condensation; the effects of internal heating and vertical mixing ($K_{\textrm{zz}}$) are explored for GJ 436b; we test updated planetary properties for 51 Eridani b with S$_8$ condensates. We find sulfur can couple to carbon or nitrogen and impact other species such as hydrogen, methane, and ammonia. The observable features of the synthetic spectra and trends in the photochemical haze precursors are discussed for each case.

preprint2021arXiv

EOS: Atmospheric Radiative Transfer in Habitable Worlds with HELIOS

We present EOS, a procedure for determining the Outgoing Longwave Radiation (OLR) and top-of-atmosphere (TOA) albedo for a wide range of conditions expected to be present in the atmospheres of rocky planets with temperate conditions. EOS is based on HELIOS and HELIOS-K, which are novel and publicly available atmospheric radiative transfer (RT) codes optimized for fast calculations with GPU processors. These codes were originally developed for the study of giant planets. In this paper we present an adaptation for applications to terrestrial-type, habitable planets, adding specific physical recipes for the gas opacity and vertical structure of the atmosphere. To test the reliability of the procedure we assessed the impact of changing line opacity profile, continuum opacity model, atmospheric lapse rate and tropopause position prescriptions on the OLR and the TOA albedo. The results obtained with EOS are in line with those of other RT codes running on traditional CPU processors, while being at least one order of magnitude faster. The adoption of OLR and TOA albedo data generated with EOS in a zonal and seasonal climate model correctly reproduce the fluxes of the present-day Earth measured by the CERES spacecraft. The results of this study disclose the possibility to incorporate fast RT calculations in climate models aimed at characterizing the atmospheres of habitable exoplanets.

preprint2021arXiv

No umbrella needed: Confronting the hypothesis of iron rain on WASP-76b with post-processed general circulation models

High-resolution spectra are unique indicators of three-dimensional processes in exoplanetary atmospheres. For instance, in 2020, Ehrenreich et al. reported transmission spectra from the ESPRESSO spectrograph yielding an anomalously large Doppler blueshift from the ultra-hot Jupiter WASP-76b. Interpretations of these observations invoke toy model depictions of gas-phase iron condensation in lower-temperature regions of the planet's atmosphere. In this work, we forward model the atmosphere of WASP-76b with double-gray general circulation models (GCMs) and ray-striking radiative transfer to diagnose the planet's high-resolution transmission spectrum. We confirm that a physical mechanism driving strong east-west asymmetries across the terminator must exist to reproduce large Doppler blueshifts in WASP-76b's transmission spectrum. We identify low atmospheric drag and a deep radiative-convective boundary as necessary components of our GCM to produce this asymmetry (the latter is consistent with existing Spitzer phase curves). However, we cannot reproduce either the magnitude or the time-dependence of the WASP-76b Doppler signature with gas-phase iron condensation alone. Instead, we find that high-altitude, optically thick clouds composed of $\rm Al_2O_3$, Fe, or $\rm Mg_2SiO_4$ provide reasonable fits to the Ehrenreich et al. observations -- with marginal contributions from condensation. This fit is further improved by allowing a small orbital eccentricity ($e \approx 0.017$), consistent with prior WASP-76b orbital constraints. We additionally validate our forward-modeled spectra by reproducing lines of nearly all species detected in WASP-76b by Tabernero et al. 2021. Our procedure's success in diagnosing phase-resolved Doppler shifts demonstrates the benefits of physical, self-consistent, three-dimensional simulations in modeling high-resolution spectra of exoplanet atmospheres.

preprint2020arXiv

Evidence for H$_{2}$ Dissociation and Recombination Heat Transport in the Atmosphere of KELT-9b

Phase curve observations provide an opportunity to study the full energy budgets of exoplanets by quantifying the amount of heat redistributed from their daysides to their nightsides. Theories explaining the properties of phase curves for hot Jupiters have focused on the balance between radiation and dynamics as the primary parameter controlling heat redistribution. However, recent phase curves have shown deviations from the trends that emerge from this theory, which has led to work on additional processes that may affect hot Jupiter energy budgets. One such process, molecular hydrogen dissociation and recombination, can enhance energy redistribution on ultra-hot Jupiters with temperatures above $\sim2000$ K. In order to study the impact of H$_{2}$ dissociation on ultra-hot Jupiters, we present a phase curve of KELT-9b observed with the Spitzer Space Telescope at 4.5 $μ$m. KELT-9b is the hottest known transiting planet, with a 4.5-$μ$m dayside brightness temperature of $4566^{+140}_{-136}$ K and a nightside temperature of $2556^{+101}_{-97}$ K. We observe a phase curve amplitude of $0.609 \pm 0.020$ and a hot spot offset of $18.7^{+2.1}_{-2.3}$ degrees. The observed amplitude is too small to be explained by a simple balance between radiation and advection. General circulation models (GCMs) and an energy balance model that include the effects of H$_{2}$ dissociation and recombination provide a better match to the data. The GCMs, however, predict a maximum hot spot offset of $5$ degrees, which disagrees with our observations at $>5σ$ confidence. This discrepancy may be due to magnetic effects in the planet's highly ionized atmosphere.

preprint2019arXiv

Identifying Atmospheres on Rocky Exoplanets Through Inferred High Albedo

The upcoming launch of the James Webb Space Telescope (JWST) means that we will soon have the capability to characterize the atmospheres of rocky exoplanets. However, it is still unknown whether such planets orbiting close to M dwarf stars can retain their atmospheres, or whether high-energy irradiation from the star will strip the gaseous envelopes from these objects. We present a new method to detect an atmosphere on a synchronously rotating rocky exoplanet around a K/M dwarf, by using thermal emission during secondary eclipse to infer a high dayside albedo that could only be explained by bright clouds. Based on calculations for plausible surface conditions, we conclude that a high albedo could be unambiguously interpreted as a signal of an atmosphere for planets with substellar temperatures of $T_{sub}=$ 410-1250 K. This range corresponds to equilibrium temperatures of $T_{eq}=$ 300-880 K. We compare the inferred albedos of eight possible planet surface compositions to cloud albedo calculations. We determine that a layer of clouds with optical depths greater than $τ=0.5$-$7$, would have high enough albedos to be distinguishable from a bare rock surface. This method of detecting an atmosphere on a rocky planet is complementary to existing methods for detecting atmospheres, because it provides a way to detect atmospheres with pressures below 1 bar (e.g. Mars), which are too tenuous to transport significant heat but thick enough to host high-albedo clouds.

preprint2019arXiv

Identifying Candidate Atmospheres on Rocky M dwarf Planets via Eclipse Photometry

Most rocky planets in the galaxy orbit a cool host star, and there is large uncertainty among theoretical models whether these planets can retain an atmosphere. The James Webb Space Telescope (JWST) might be able to settle this question empirically, but most proposals for doing so require large observational effort because they are based on spectroscopy. Here we show that infrared photometry of secondary eclipses could quickly identify "candidate" atmospheres, by searching for rocky planets with atmospheres thick enough that atmospheric heat transport noticeably reduces their dayside thermal emission compared to that of a bare rock. For a planet amenable to atmospheric follow-up, we find that JWST should be able to confidently detect the heat redistribution signal of an O(1) bar atmosphere with one to two eclipses. One to two eclipses is generally much less than the effort needed to infer an atmosphere via transmission or emission spectroscopy. Candidate atmospheres can be further validated via follow-up spectroscopy or phase curves. In addition, because this technique is fast it could enable a first atmospheric survey of rocky exoplanets with JWST. We estimate that the TESS mission will find ~100 planets that are too hot to be habitable but that can be quickly probed via eclipse photometry. Knowing whether hot, rocky planets around M dwarfs have atmospheres is important not only for understanding the evolution of uninhabitable worlds: if atmospheres are common on hot planets, then cooler, potentially habitable planets around M dwarfs are also likely to have atmospheres.

preprint2019arXiv

Supervised Machine Learning for Intercomparison of Model Grids of Brown Dwarfs: Application to GJ 570D and the Epsilon Indi B Binary System

Self-consistent model grids of brown dwarfs involve complex physics and chemistry, and are often computed using proprietary computer codes, making it challenging to identify the reasons for discrepancies between model and data as well as between the models produced by different research groups. In the current study, we demonstrate a novel method for analyzing brown dwarf spectra, which combines the use of the Sonora, AMES-Cond and HELIOS model grids with the supervised machine learning method of the random forest. Besides performing atmospheric retrieval, the random forest enables information content analysis of the three model grids as a natural outcome of the method, both individually on each grid and by comparing the grids against one another, via computing large suites of mock retrievals. Our analysis reveals that the different choices made in modelling the alkali line shapes hinder the use of the alkali lines as gravity indicators. Nevertheless, the spectrum longward of 1.2 micron encodes enough information on the surface gravity to allow its inference from retrieval. Temperature may be accurately and precisely inferred independent of the choice of model grid, but not the surface gravity. We apply random forest retrieval to three objects: the benchmark T7.5 brown dwarf GJ 570D; and Epsilon Indi Ba (T1.5 brown dwarf) and Bb (T6 brown dwarf), which are part of a binary system and have measured dynamical masses. For GJ 570D, the inferred effective temperature and surface gravity are consistent with previous studies. For Epsilon Indi Ba and Bb, the inferred surface gravities are broadly consistent with the values informed by the dynamical masses.