Showing 1–32 of 32 results for author: Mendonca, J M

Enhancing 3D Planetary Atmosphere Simulations with a Surrogate Radiative Transfer Model

Authors: Tara P. A. Tahseen, João M. Mendonça, Kai Hou Yip, Ingo P. Waldmann

Abstract: This work introduces an approach to enhancing the computational efficiency of 3D atmospheric simulations by integrating a machine-learned surrogate model into the OASIS global circulation model (GCM). Traditional GCMs, which are based on repeatedly numerically integrating physical equations governing atmospheric processes across a series of time-steps, are time-intensive, leading to compromises in… ▽ More This work introduces an approach to enhancing the computational efficiency of 3D atmospheric simulations by integrating a machine-learned surrogate model into the OASIS global circulation model (GCM). Traditional GCMs, which are based on repeatedly numerically integrating physical equations governing atmospheric processes across a series of time-steps, are time-intensive, leading to compromises in spatial and temporal resolution of simulations. This research improves upon this limitation, enabling higher resolution simulations within practical timeframes. Speeding up 3D simulations holds significant implications in multiple domains. Firstly, it facilitates the integration of 3D models into exoplanet inference pipelines, allowing for robust characterisation of exoplanets from a previously unseen wealth of data anticipated from JWST and post-JWST instruments. Secondly, acceleration of 3D models will enable higher resolution atmospheric simulations of Earth and Solar System planets, enabling more detailed insights into their atmospheric physics and chemistry. Our method replaces the radiative transfer module in OASIS with a recurrent neural network-based model trained on simulation inputs and outputs. Radiative transfer is typically one of the slowest components of a GCM, thus providing the largest scope for overall model speed-up. The surrogate model was trained and tested on the specific test case of the Venusian atmosphere, to benchmark the utility of this approach in the case of non-terrestrial atmospheres. This approach yields promising results, with the surrogate-integrated GCM demonstrating above 99.0% accuracy and 101 factor GPU speed-up of the entire simulation compared to using the matched original GCM under Venus-like conditions. △ Less

Submitted 11 July, 2024; originally announced July 2024.

Comments: 17 pages, 11 figures

High-energy spectra of LTT 1445A and GJ 486 reveal flares and activity

Authors: H. Diamond-Lowe, G. W. King, A. Youngblood, A. Brown, W. S. Howard, J. G. Winters, D. J. Wilson, K. France, J. M. Mendonça, L. A. Buchhave, L. Corrales, L. Kreidberg, A. A. Medina, J. L. Bean, Z. K. Berta-Thompson, T. M. Evans-Soma, C. Froning, G. M. Duvvuri, E. M. -R. Kempton, Y. Miguel, J. S. Pineda, C. Schneider

Abstract: The high-energy radiative output, from the X-ray to the ultraviolet, of exoplanet host stars drives photochemical reactions and mass loss in the upper regions of planetary atmospheres. In order to place constraints on the atmospheric properties of the three closest terrestrial exoplanets transiting M dwarfs, we observe the high-energy spectra of the host stars LTT1445A and GJ486 in the X-ray with… ▽ More The high-energy radiative output, from the X-ray to the ultraviolet, of exoplanet host stars drives photochemical reactions and mass loss in the upper regions of planetary atmospheres. In order to place constraints on the atmospheric properties of the three closest terrestrial exoplanets transiting M dwarfs, we observe the high-energy spectra of the host stars LTT1445A and GJ486 in the X-ray with XMM-Newton and Chandra and in the ultraviolet with HST/COS and STIS. We combine these observations with estimates of extreme ultraviolet flux, reconstructions of the Ly-a lines, and stellar models at optical and infrared wavelengths to produce panchromatic spectra from 1A--20um for each star. While LTT1445Ab, LTT1445Ac, and GJ486b do not possess primordial hydrogen-dominated atmospheres, we calculate that they are able to retain pure CO2 atmospheres if starting with 10, 15, and 50% of Earth's total CO2 budget, respectively, in the presence of their host stars' stellar wind. We use age-activity relationships to place lower limits of 2.2 and 6.6 Gyr on the ages of the host stars LTT1445A and GJ486. Despite both LTT1445A and GJ486 appearing inactive at optical wavelengths, we detect flares at ultraviolet and X-ray wavelengths for both stars. In particular, GJ486 exhibits two flares with absolute energies of 10^29.5 and 10^30.1 erg (equivalent durations of 4357+/-96 and 19724+/-169 s) occurring three hours apart, captured with HST/COS G130M. Based on the timing of the observations, we suggest that these high-energy flares are related and indicative of heightened flaring activity that lasts for a period of days, but our interpretations are limited by sparse time-sampling. Consistent high-energy monitoring is needed to determine the duration and extent of high-energy activity on individual M dwarfs, as well as the population as a whole. △ Less

Submitted 28 June, 2024; originally announced July 2024.

Comments: 22 pages, accepted to A&A

Simultaneous retrieval of orbital phase resolved JWST/MIRI emission spectra of the hot Jupiter WASP-43b: evidence of water, ammonia and carbon monoxide

Authors: Jingxuan Yang, Mark Hammond, Anjali A. A. Piette, Jasmina Blecic, Taylor J. Bell, Patrick G. J. Irwin, Vivien Parmentier, Shang-Min Tsai, Joanna K. Barstow, Nicolas Crouzet, Laura Kreidberg, João M. Mendonça, Jake Taylor, Robin Baeyens, Kazumasa Ohno, Lucas Teinturier, Matthew C. Nixon

Abstract: Spectroscopic phase curves of hot Jupiters measure their emission spectra at multiple orbital phases, thus enabling detailed characterisation of their atmospheres. Precise constraints on the atmospheric composition of these exoplanets offer insights into their formation and evolution. We analyse four phase-resolved emission spectra of the hot Jupiter WASP-43b, generated from a phase curve observed… ▽ More Spectroscopic phase curves of hot Jupiters measure their emission spectra at multiple orbital phases, thus enabling detailed characterisation of their atmospheres. Precise constraints on the atmospheric composition of these exoplanets offer insights into their formation and evolution. We analyse four phase-resolved emission spectra of the hot Jupiter WASP-43b, generated from a phase curve observed with the MIRI/LRS onboard the JWST, to retrieve its atmospheric properties. Using a parametric 2D temperature model and assuming a chemically homogeneous atmosphere within the observed pressure region, we simultaneously fit the four spectra to constrain the abundances of atmospheric constituents, thereby yielding more precise constraints than previous work that analysed each spectrum independently. Our analysis reveals statistically significant evidence of NH3 (4$σ$) in a hot Jupiter's emission spectra for the first time, along with evidence of H2O (6.5$σ$), CO (3.1$σ$), and a non-detection of CH4. With our abundance constraints, we tentatively estimate the metallicity of WASP-43b at 0.6-6.5$\times$solar and its C/O ratio at 0.6-0.9. Our findings offer vital insights into the atmospheric conditions and formation history of WASP-43b by simultaneously constraining the abundances of carbon, oxygen, and nitrogen-bearing species. △ Less

Submitted 5 June, 2024; originally announced June 2024.

Comments: 17 pages, 14 figures, accepted for publication in MNRAS. Comments welcome!

Identifying and Fitting Eclipse Maps of Exoplanets with Cross-Validation

Authors: Mark Hammond, Neil T. Lewis, Sasha Boone, Xueqing Chen, João M. Mendonça, Vivien Parmentier, Jake Taylor, Taylor Bell, Leonardo dos Santos, Nicolas Crouzet, Laura Kreidberg, Michael Radica, Michael Zhang

Abstract: Eclipse mapping uses the shape of the eclipse of an exoplanet to measure its two-dimensional structure. Light curves are mostly composed of longitudinal information, with the latitudinal information only contained in the brief ingress and egress of the eclipse. This imbalance can lead to a spuriously confident map, where the longitudinal structure is constrained by out-of-eclipse data and the lati… ▽ More Eclipse mapping uses the shape of the eclipse of an exoplanet to measure its two-dimensional structure. Light curves are mostly composed of longitudinal information, with the latitudinal information only contained in the brief ingress and egress of the eclipse. This imbalance can lead to a spuriously confident map, where the longitudinal structure is constrained by out-of-eclipse data and the latitudinal structure is wrongly determined by the priors on the map. We present a new method to address this issue. The method tests for the presence of an eclipse mapping signal by using k-fold cross-validation to compare the performance of a simple mapping model to the null hypothesis of a uniform disk. If a signal is found, the method fits a map with more degrees of freedom, optimising its information content. The information content is varied by penalising the model likelihood by a factor proportional to the spatial entropy of the map, optimised by cross-validation. We demonstrate this method for simulated datasets then apply it to three observational datasets. The method identifies an eclipse mapping signal for JWST MIRI/LRS observations of WASP-43b but does not identify a signal for JWST NIRISS/SOSS observations of WASP-18b or Spitzer Space Telescope observations of HD 189733b. It is possible to fit eclipse maps to these datasets, but we suggest that these maps are overfitting the eclipse shape. We fit a new map with more spatial freedom to the WASP-43b dataset and show a flatter east-west structure than previously derived. △ Less

Submitted 31 May, 2024; originally announced May 2024.

Comments: Submitted to MNRAS

Two-Dimensional Eclipse Mapping of the Hot Jupiter WASP-43b with JWST MIRI/LRS

Authors: Mark Hammond, Taylor J. Bell, Ryan C. Challener, Neil T. Lewis, Megan Weiner Mansfield, Isaac Malsky, Emily Rauscher, Jacob L. Bean, Ludmila Carone, João M. Mendonça, Lucas Teinturier, Xianyu Tan, Nicolas Crouzet, Laura Kreidberg, Giuseppe Morello, Vivien Parmentier, Jasmina Blecic, Jean-Michel Désert, Christiane Helling, Pierre-Olivier Lagage, Karan Molaverdikhani, Matthew C. Nixon, Benjamin V. Rackham, Jingxuan Yang

Abstract: We present eclipse maps of the two-dimensional thermal emission from the dayside of the hot Jupiter WASP-43b, derived from an observation of a phase curve with the JWST MIRI/LRS instrument. The observed eclipse shapes deviate significantly from those expected for a planet emitting uniformly over its surface. We fit a map to this deviation, constructed from spherical harmonics up to order… ▽ More We present eclipse maps of the two-dimensional thermal emission from the dayside of the hot Jupiter WASP-43b, derived from an observation of a phase curve with the JWST MIRI/LRS instrument. The observed eclipse shapes deviate significantly from those expected for a planet emitting uniformly over its surface. We fit a map to this deviation, constructed from spherical harmonics up to order $\ell_{\rm max}=2$, alongside the planetary, orbital, stellar, and systematic parameters. This yields a map with a meridionally-averaged eastward hot-spot shift of $(7.75 \pm 0.36)^{\circ}$, with no significant degeneracy between the map and the additional parameters. We show the latitudinal and longitudinal contributions of the day-side emission structure to the eclipse shape, finding a latitudinal signal of $\sim$200 ppm and a longitudinal signal of $\sim$250 ppm. To investigate the sensitivity of the map to the method, we fix the non-mapping parameters and derive an "eigenmap" fitted with an optimised number of orthogonal phase curves, which yields a similar map to the $\ell_{\rm max}=2$ map. We also fit a map up to $\ell_{\rm max}=3$, which shows a smaller hot-spot shift, with a larger uncertainty. These maps are similar to those produced by atmospheric simulations. We conclude that there is a significant mapping signal which constrains the spherical harmonic components of our model up to $\ell_{\rm max}=2$. Alternative mapping models may derive different structures with smaller-scale features; we suggest that further observations of WASP-43b and other planets will drive the development of more robust methods and more accurate maps. △ Less

Submitted 25 April, 2024; originally announced April 2024.

Comments: Accepted for publication in The Astronomical Journal

Nightside clouds and disequilibrium chemistry on the hot Jupiter WASP-43b

Authors: Taylor J. Bell, Nicolas Crouzet, Patricio E. Cubillos, Laura Kreidberg, Anjali A. A. Piette, Michael T. Roman, Joanna K. Barstow, Jasmina Blecic, Ludmila Carone, Louis-Philippe Coulombe, Elsa Ducrot, Mark Hammond, João M. Mendonça, Julianne I. Moses, Vivien Parmentier, Kevin B. Stevenson, Lucas Teinturier, Michael Zhang, Natalie M. Batalha, Jacob L. Bean, Björn Benneke, Benjamin Charnay, Katy L. Chubb, Brice-Olivier Demory, Peter Gao , et al. (58 additional authors not shown)

Abstract: Hot Jupiters are among the best-studied exoplanets, but it is still poorly understood how their chemical composition and cloud properties vary with longitude. Theoretical models predict that clouds may condense on the nightside and that molecular abundances can be driven out of equilibrium by zonal winds. Here we report a phase-resolved emission spectrum of the hot Jupiter WASP-43b measured from 5… ▽ More Hot Jupiters are among the best-studied exoplanets, but it is still poorly understood how their chemical composition and cloud properties vary with longitude. Theoretical models predict that clouds may condense on the nightside and that molecular abundances can be driven out of equilibrium by zonal winds. Here we report a phase-resolved emission spectrum of the hot Jupiter WASP-43b measured from 5-12 $μ$m with JWST's Mid-Infrared Instrument (MIRI). The spectra reveal a large day-night temperature contrast (with average brightness temperatures of 1524$\pm$35 and 863$\pm$23 Kelvin, respectively) and evidence for water absorption at all orbital phases. Comparisons with three-dimensional atmospheric models show that both the phase curve shape and emission spectra strongly suggest the presence of nightside clouds which become optically thick to thermal emission at pressures greater than ~100 mbar. The dayside is consistent with a cloudless atmosphere above the mid-infrared photosphere. Contrary to expectations from equilibrium chemistry but consistent with disequilibrium kinetics models, methane is not detected on the nightside (2$σ$ upper limit of 1-6 parts per million, depending on model assumptions). △ Less

Submitted 23 January, 2024; originally announced January 2024.

Comments: 61 pages, 13 figures, 4 tables. This preprint has been submitted to and accepted in principle for publication in Nature Astronomy without significant changes

Global Chemical Transport on Hot Jupiters: Insights from 2D VULCAN photochemical model

Authors: Shang-Min Tsai, Vivien Parmentier, João M. Mendonça, Xianyu Tan, Russell Deitrick, Mark Hammond, Arjun B. Savel, Xi Zhang, Raymond T. Pierrehumbert, Edward W. Schwieterman

Abstract: The atmospheric dynamics of tidally-locked hot Jupiters is characterized by strong equatorial winds. Understanding the interaction between global circulation and chemistry is crucial in atmospheric studies and interpreting observations. Two-dimensional (2D) photochemical transport models shed light on how the atmospheric composition depends on circulation. In this paper, we introduce the 2D photoc… ▽ More The atmospheric dynamics of tidally-locked hot Jupiters is characterized by strong equatorial winds. Understanding the interaction between global circulation and chemistry is crucial in atmospheric studies and interpreting observations. Two-dimensional (2D) photochemical transport models shed light on how the atmospheric composition depends on circulation. In this paper, we introduce the 2D photochemical (horizontal and vertical) transport model, VULCAN 2D, which improves on the pseudo-2D approaches by allowing for non-uniform zonal winds. We extensively validate our VULCAN 2D with analytical solutions and benchmark comparisons. Applications to HD 189733 b and HD 209458 b reveal a transition in mixing regimes: horizontal transport predominates below $\sim$0.1 mbar while vertical mixing is more important at higher altitudes above 0.1 mbar. Motivated by the previously inferred carbon-rich atmosphere, we find that HD 209458 b with super-solar carbon-to-oxygen ratio (C/O) exhibits pronounced C$_2$H$_4$ absorption on the morning limb but not on the evening limb, owing to horizontal transport from the nightside. We discuss when a pseudo-2D approach is a valid assumption and its inherent limitations. Finally, we demonstrate the effect of horizontal transport in transmission observations and its impact on the morning-evening limb asymmetry with synthetic spectra, highlighting the need to consider global transport when interpreting exoplanet atmospheres. △ Less

Submitted 25 March, 2024; v1 submitted 26 October, 2023; originally announced October 2023.

Comments: 19 pages, 20 figures, published in ApJ

Journal ref: Tsai et al. (2024), ApJ, 963, 41

Transmission spectroscopy of the lowest-density gas giant: metals and a potential extended outflow in HAT-P-67b

Authors: Aaron Bello-Arufe, Heather A. Knutson, João M. Mendonça, Michael M. Zhang, Samuel H. C. Cabot, Alexander D. Rathcke, Ana Ulla, Shreyas Vissapragada, Lars A. Buchhave

Abstract: Extremely low-density exoplanets are tantalizing targets for atmospheric characterization because of their promisingly large signals in transmission spectroscopy. We present the first analysis of the atmosphere of the lowest-density gas giant currently known, HAT-P-67 b. This inflated Saturn-mass exoplanet sits at the boundary between hot and ultrahot gas giants, where thermal dissociation of mole… ▽ More Extremely low-density exoplanets are tantalizing targets for atmospheric characterization because of their promisingly large signals in transmission spectroscopy. We present the first analysis of the atmosphere of the lowest-density gas giant currently known, HAT-P-67 b. This inflated Saturn-mass exoplanet sits at the boundary between hot and ultrahot gas giants, where thermal dissociation of molecules begins to dominate atmospheric composition. We observed a transit of HAT-P-67 b at high spectral resolution with CARMENES and searched for atomic and molecular species using cross-correlation and likelihood mapping. Furthermore, we explored potential atmospheric escape by targeting H$α$ and the metastable helium line. We detect Ca II and Na I with significances of 13.2$σ$ and 4.6$σ$, respectively. Unlike in several ultrahot Jupiters, we do not measure a day-to-night wind. The large line depths of Ca II suggest that the upper atmosphere may be more ionized than models predict. We detect strong variability in H$α$ and the helium triplet during the observations. These signals suggest the possible presence of an extended planetary outflow that causes an early ingress and late egress. In the averaged transmission spectrum, we measure redshifted absorption at the $\sim 3.8\%$ and $\sim 4.5\%$ level in the H$α$ and He I triplet lines, respectively. From an isothermal Parker wind model, we derive a mass loss rate of $\dot{M} \sim 10^{13}~\rm{g/s}$ and an outflow temperature of $T \sim 9900~\rm{K}$. However, due to the lack of a longer out-of-transit baseline in our data, additional observations are needed to rule out stellar variability as the source of the H$α$ and He signals. △ Less

Submitted 12 July, 2023; originally announced July 2023.

Comments: The Astronomical Journal, in press. 17 pages, 9 figures

Multi-epoch detections of the extended atmosphere and transmission spectra of KELT-9b with a 1.5 m telescope

Authors: Nataliea Lowson, George Zhou, Duncan J. Wright, Chelsea X. Huang, Joao M. Mendonca, Samuel H. C. Cabot, Christa Pudmenzky, Robert A. Wittenmyer, David W. Latham, Allyson Bieryla, Gilbert A. Esquerdo, Perry Berlind, Michael L. Calkins

Abstract: Irradiated Jovian atmospheres are complex, dynamic, and can undergo temporal variations due to the close proximity of their parent stars. Of the Jovian planets that have been catalogued to date, KELT-9b is the hottest Gas Giant known, with an equilibrium temperature of 4050 K. We probe the temporal variability of transmission spectroscopic signatures from KELT-9b via a set of archival multi-year g… ▽ More Irradiated Jovian atmospheres are complex, dynamic, and can undergo temporal variations due to the close proximity of their parent stars. Of the Jovian planets that have been catalogued to date, KELT-9b is the hottest Gas Giant known, with an equilibrium temperature of 4050 K. We probe the temporal variability of transmission spectroscopic signatures from KELT-9b via a set of archival multi-year ground-based transit observations, performed with the TRES facility on the 1.5 m reflector at the Fred Lawrence Whipple Observatory. Our observations confirm past detections of Fe I, Fe II and Mg I over multiple epochs, in addition to excess absorption at H-alpha, which is an indicator for ongoing mass-loss. From our multi-year dataset, the H-alpha light curve consistently deviates from a standard transit, and follows a 'W' shape that is deeper near ingress and egress, and shallower mid-transit. To search for and quantify any seasonal variations that may be present, we parameterise a 'cometary tail' model to fit for the H-alpha transit. We find no detectable variations between the different observed epochs. Though a 'cometary tail' describes the H-alpha flux variations well, we note that such a scenario requires a high density of neutral hydrogen in the n = 2 excited state far beyond the planetary atmosphere. Other scenarios, such as centre-to-limb variations larger than that expected from 1-D atmosphere models, may also contribute to the observed H-alpha transit shape. These multi-epoch observations highlight the capabilities of small telescopes to provide temporal monitoring of the dynamics of exoplanet atmospheres. △ Less

Submitted 24 January, 2023; originally announced January 2023.

Comments: 19 pages, 11 figures, accepted for publication in AJ

Ground-based Optical Transmission Spectroscopy of the Nearby Terrestrial Exoplanet LTT 1445Ab

Authors: Hannah Diamond-Lowe, João M. Mendonça, David Charbonneau, Lars A. Buchhave

Abstract: Nearby M dwarf systems currently offer the most favorable opportunities for spectroscopic investigations of terrestrial exoplanet atmospheres. The LTT~1445 system is a hierarchical triple of M dwarfs with two known planets orbiting the primary star, LTT~1445A. We observe four transits of the terrestrial world LTT~1445Ab ($R=1.3$ R$_\oplus$, $M=2.9$ M$_\oplus$) at low resolution with Magellan II/LD… ▽ More Nearby M dwarf systems currently offer the most favorable opportunities for spectroscopic investigations of terrestrial exoplanet atmospheres. The LTT~1445 system is a hierarchical triple of M dwarfs with two known planets orbiting the primary star, LTT~1445A. We observe four transits of the terrestrial world LTT~1445Ab ($R=1.3$ R$_\oplus$, $M=2.9$ M$_\oplus$) at low resolution with Magellan II/LDSS3C. We use the combined flux of the LTT~1445BC pair as a comparison star, marking the first time that an M dwarf is used to remove telluric variability from time-series observations of another M dwarf. We find H$α$ in emission from both LTT~1445B and C, as well as a flare in one of the data sets from LTT~1445C. These contaminated data are removed from the analysis. We construct a broadband transit light curve of LTT~1445Ab from 620--1020 nm. Binned to 3-minute time bins we achieve an rms of 49 ppm for the combined broadband light curve. We construct a transmission spectrum with 20 spectrophotometric bins each spanning 20 nm and compare it to models of clear, 1$\times$ solar composition atmospheres. We rule out this atmospheric case with a surface pressure of 10 bars to $3.2 σ$ confidence, and with a surface pressure of 1 bar to $3.1 σ$ confidence. Upcoming secondary eclipse observations of LTT~1445Ab with JWST will further probe the cases of a high mean molecular weight atmosphere, a hazy or cloudy atmosphere, or no atmosphere at all on this terrestrial world. △ Less

Submitted 25 February, 2023; v1 submitted 21 October, 2022; originally announced October 2022.

Comments: Accepted for publication in AJ

Is ozone a reliable proxy for molecular oxygen? I. The O2-O3 relationship for Earth-like atmospheres

Authors: Thea Kozakis, João M. Mendonça, Lars A. Buchhave

Abstract: Molecular oxygen (O2) paired with a reducing gas is regarded as a promising biosignature pair for the atmospheric characterization of terrestrial exoplanets. In circumstances when O2 may not be detectable in a planetary atmosphere (e.g., at mid-IR wavelengths) it has been suggested that ozone (O3), the photochemical product of O2, could be used as a proxy to infer the presence of O2. However, O3 p… ▽ More Molecular oxygen (O2) paired with a reducing gas is regarded as a promising biosignature pair for the atmospheric characterization of terrestrial exoplanets. In circumstances when O2 may not be detectable in a planetary atmosphere (e.g., at mid-IR wavelengths) it has been suggested that ozone (O3), the photochemical product of O2, could be used as a proxy to infer the presence of O2. However, O3 production has a nonlinear dependence on O2 and is strongly influenced by the UV spectrum of the host star. To evaluate the reliability of O3 as a proxy for O2, we used Atmos, a 1D coupled climate/photochemistry code, to study the O2-O3 relationship for "Earth-like" habitable zone planets around a variety of stellar hosts (G0V-M5V) and O2 abundances. Overall, we found that the O2-O3 relationship differed significantly with stellar hosts and resulted in different trends for hotter stars (G0V-K2V) vs cooler stars (K5V-M5V). Planets orbiting hotter host stars counter-intuitively experience an increase in O3 when O2 levels are initially decreased from 100% Earth's present atmospheric level (PAL), with a maximum O3 abundance occurring at 25-55% PAL O2. As O2 abundance initially decreases, larger amounts of UV photons capable of O2 photolysis reach the lower (denser) regions of the atmosphere where O3 production is more efficient, resulting in these increased O3 levels. This effect does not occur for cooler host stars (K5V-M5V), as the weaker incident UV flux does not allow O3 formation to occur at dense enough regions of the atmosphere where the faster O3 production can outweigh a smaller source of O2 from which to create O3. Overall it will be extremely difficult (or impossible) to infer precise O2 levels from an O3 measurement, however, with information about the UV spectrum of the host star and context clues, O3 will provide valuable information about potential surface habitability of an exoplanet. △ Less

Submitted 19 August, 2022; originally announced August 2022.

Comments: 19 pages, 7 figures, 4 tables; in press with A & A

Exoplanet atmospheres at high resolution through a modest-size telescope. Fe II in MASCARA-2b and KELT-9b with FIES on the Nordic Optical Telescope

Authors: Aaron Bello-Arufe, Lars A. Buchhave, João M. Mendonça, René Tronsgaard, Kevin Heng, H. Jens Hoeijmakers, Andrew W. Mayo

Abstract: Ground-based, high-resolution spectrographs are providing us with an unprecedented view of the dynamics and chemistry of the atmospheres of planets outside the Solar System. While there is a large number of stable and precise high-resolution spectrographs on modest-size telescopes, it is the spectrographs at observatories with apertures larger than 3.5 metres that dominate the atmospheric follow-u… ▽ More Ground-based, high-resolution spectrographs are providing us with an unprecedented view of the dynamics and chemistry of the atmospheres of planets outside the Solar System. While there is a large number of stable and precise high-resolution spectrographs on modest-size telescopes, it is the spectrographs at observatories with apertures larger than 3.5 metres that dominate the atmospheric follow-up of exoplanets. In this work, we explore the potential of characterising exoplanetary atmospheres with FIES, a high-resolution spectrograph at the 2.56 metre Nordic Optical Telescope. We observed two transits of MASCARA-2 b (also known as KELT-20 b) and one transit of KELT-9 b to search for atomic iron, a species that has been recently discovered in both neutral and ionised forms in the atmospheres of these ultra-hot Jupiters using large telescopes. Using a cross-correlation method, we detect a signal of Fe II at the $4.5σ$ and $4.0σ$ level in the transits of MASCARA-2 b. We also detect Fe II in the transit of KELT-9 b at the $8.5σ$ level. Although we do not find any significant Doppler shift in the signal of MASCARA-2 b, we do measure a moderate blueshift (3-6 km/s) of the feature in KELT-9 b, which might be a manifestation of high-velocity winds transporting Fe II from the planetary dayside to the nightside. Our work demonstrates the feasibility of investigating exoplanet atmospheres with FIES, potentially unlocking a wealth of additional atmosphere detections with this and other high-resolution spectrographs mounted on similar-size telescopes. △ Less

Submitted 14 March, 2022; v1 submitted 9 March, 2022; originally announced March 2022.

Comments: 9 pages, 6 figures, accepted to A&A

Journal ref: A&A 662, A51 (2022)

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

Authors: Russell Deitrick, Kevin Heng, Urs Schroffenegger, Daniel Kitzman, Simon L. Grimm, Matej Malik, João M. Mendonça, Brett M. Morris

Abstract: 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 a… ▽ More 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. △ Less

Submitted 7 March, 2022; v1 submitted 4 March, 2022; originally announced March 2022.

Comments: 31 pages, 24 figures, accepted for publication at MNRAS

Mining the Ultra-Hot Skies of HAT-P-70b: Detection of a Profusion of Neutral and Ionized Species

Authors: Aaron Bello-Arufe, Samuel H. C. Cabot, João M. Mendonça, Lars A. Buchhave, Alexander D. Rathcke

Abstract: With an equilibrium temperature above 2500 K, the recently discovered HAT-P-70 b belongs to a new class of exoplanets known as ultra-hot Jupiters: extremely irradiated gas giants with day-side temperatures that resemble those found in stars. These ultra-hot Jupiters are among the most amenable targets for follow-up atmospheric characterization through transmission spectroscopy. Here, we present th… ▽ More With an equilibrium temperature above 2500 K, the recently discovered HAT-P-70 b belongs to a new class of exoplanets known as ultra-hot Jupiters: extremely irradiated gas giants with day-side temperatures that resemble those found in stars. These ultra-hot Jupiters are among the most amenable targets for follow-up atmospheric characterization through transmission spectroscopy. Here, we present the first analysis of the transmission spectrum of HAT-P-70 b using high-resolution data from the HARPS-N spectrograph of a single transit event. We use a cross-correlation analysis and transmission spectroscopy to look for atomic and molecular species in the planetary atmosphere. We detect absorption by Ca II, Cr I, Cr II, Fe I, Fe II, H I, Mg I, Na I and V I, and we find tentative evidence of Ca I and Ti II. Overall, these signals appear blue-shifted by a few km s$^{-1}$, suggestive of winds flowing at high velocity from the day-side to the night-side. We individually resolve the Ca II H & K lines, the Na I doublet, and the H$α$, H$β$ and H$γ$ Balmer lines. The cores of the Ca II and H I lines form well above the continuum, indicating the existence of an extended envelope. We refine the obliquity of this highly misaligned planet to $107.9^{+2.0}_{-1.7}$ degrees by examining the Doppler shadow that the planet casts on its A-type host star. These results place HAT-P-70 b as one of the exoplanets with the highest number of species detected in its atmosphere. △ Less

Submitted 6 December, 2021; originally announced December 2021.

Comments: 17 pages, 10 figures, accepted to AJ

General Circulation Model Errors are Variable across Exoclimate Parameter Spaces

Authors: Pushkar Kopparla, Russell Deitrick, Kevin Heng, João M. Mendonça, Mark Hammond

Abstract: General circulation models are often used to explore exoclimate parameter spaces and classify atmospheric circulation regimes. Models are tuned to give reasonable climate states for standard test cases, such as the Held-Suarez test, and then used to simulate diverse exoclimates by varying input parameters such as rotation rates, instellation, atmospheric optical properties, frictional timescales a… ▽ More General circulation models are often used to explore exoclimate parameter spaces and classify atmospheric circulation regimes. Models are tuned to give reasonable climate states for standard test cases, such as the Held-Suarez test, and then used to simulate diverse exoclimates by varying input parameters such as rotation rates, instellation, atmospheric optical properties, frictional timescales and so on. In such studies, there is an implicit assumption that the model which works reasonably well for the standard test case will be credible at all points in an arbitrarily wide parameter space. Here, we test this assumption using the open-source general circulation model THOR to simulate atmospheric circulation on tidally locked Earth-like planets with rotation periods of 0.1 to 100 days. We find that the model error, as quantified by the ratio between physical and spurious numerical contributions to the angular momentum balance, is extremely variable across this range of rotation periods with some cases where numerical errors are the dominant component. Increasing model grid resolution does improve errors but using a higher-order numerical diffusion scheme can sometimes magnify errors for finite-volume dynamical solvers. We further show that to minimize error and make the angular momentum balance more physical within our model, the surface friction timescale must be smaller than the rotational timescale. △ Less

Submitted 21 October, 2021; originally announced October 2021.

Comments: 8 pages, 5 figures, accepted by ApJ

Mass Transport in a Moist Planetary Climate Model

Authors: João M. Mendonça

Abstract: Planetary Climate Models (PCMs) are developed to explore planetary climates other than the Earth. Therefore, the methods implemented need to be suitable for a large diversity of conditions. Every planet with a significant atmosphere has condensible cycles (e.g., hydrological cycle), which can play an essential role in the planet's appearance and environment. We must accurately represent a condensi… ▽ More Planetary Climate Models (PCMs) are developed to explore planetary climates other than the Earth. Therefore, the methods implemented need to be suitable for a large diversity of conditions. Every planet with a significant atmosphere has condensible cycles (e.g., hydrological cycle), which can play an essential role in the planet's appearance and environment. We must accurately represent a condensible cycle in our planet simulations to build a powerful planetary climate predictor. OASIS is a 3D PCM capable of self-consistently representing the main physical processes that drive a planet's environment. In this work, we improve the representation of mass transport in OASIS, which is the first step towards a complete and flexible implementation of a condensible cycle. We implement an upwind-biased scheme on a piece-wise linear approximation with a flux-limiter to solve the mass transport equation. We first benchmark the new scheme on a 2D problem that confirms the superior properties of the new method over the central finite-volume method in terms of performance, accuracy and shape-preserving mass distribution. Due to the new scheme's less dispersive nature, we do not have to apply any unphysical diffusion to maintain the model stable. OASIS includes the new improved solver in the total mass and the tracers (e.g., clouds/hazes and individual gas chemical species) transport. We couple the new formulation with physical schemes and validate the new code on two 3D simulations of an ocean Earth-like planet and an ocean tidally-locked planet. The new OASIS simulations are robust and do not show any known problems from the dynamics-physics coupling. We show that the two simulations capture the main characteristics of ocean planet atmospheres and are easy to set up. We propose these two simulations to be the first standard benchmark tests for models built to explore moist planetary environments. △ Less

Submitted 1 April, 2022; v1 submitted 7 October, 2021; originally announced October 2021.

Comments: Published version in A&A

Journal ref: A&A 659, A43 (2022)

TOI-1518b: A Misaligned Ultra-hot Jupiter with Iron in its Atmosphere

Authors: Samuel H. C. Cabot, Aaron Bello-Arufe, João M. Mendonça, René Tronsgaard, Ian Wong, George Zhou, Lars A. Buchhave, Debra A. Fischer, Keivan G. Stassun, Victoria Antoci, David Baker, Alexander A. Belinski, Björn Benneke, Luke G. Bouma, Jessie L. Christiansen, Karen A. Collins, Maria V. Goliguzova, Simone Hagey, Jon M. Jenkins, Eric L. N. Jensen, Richard C. Kidwell Jr, Didier Laloum, Bob Massey, Kim K. McLeod, David W. Latham , et al. (14 additional authors not shown)

Abstract: We present the discovery of TOI-1518b -- an ultra-hot Jupiter orbiting a bright star $V = 8.95$. The transiting planet is confirmed using high-resolution optical transmission spectra from EXPRES. It is inflated, with $R_p = 1.875\pm0.053\,R_{\rm J}$, and exhibits several interesting properties, including a misaligned orbit (${240.34^{+0.93}_{-0.98}}$ degrees) and nearly grazing transit (… ▽ More We present the discovery of TOI-1518b -- an ultra-hot Jupiter orbiting a bright star $V = 8.95$. The transiting planet is confirmed using high-resolution optical transmission spectra from EXPRES. It is inflated, with $R_p = 1.875\pm0.053\,R_{\rm J}$, and exhibits several interesting properties, including a misaligned orbit (${240.34^{+0.93}_{-0.98}}$ degrees) and nearly grazing transit ($b =0.9036^{+0.0061}_{-0.0053}$). The planet orbits a fast-rotating F0 host star ($T_{\mathrm{eff}} \simeq 7300$ K) in 1.9 days and experiences intense irradiation. Notably, the TESS data show a clear secondary eclipse with a depth of $364\pm28$ ppm and a significant phase curve signal, from which we obtain a relative day-night planetary flux difference of roughly 320 ppm and a 5.2$σ$ detection of ellipsoidal distortion on the host star. Prompted by recent detections of atomic and ionized species in ultra-hot Jupiter atmospheres, we conduct an atmospheric cross-correlation analysis. We detect neutral iron (${5.2σ}$), at $K_p = 157^{+68}_{-44}$ km s$^{-1}$ and $V_{\rm sys} = -16^{+2}_{-4}$ km s$^{-1}$, adding another object to the small sample of highly irradiated gas-giant planets with Fe detections in transmission. Detections so far favor particularly inflated gas giants with radii $\gtrsim 1.78\,R_{\rm J}$; although this may be due to observational bias. With an equilibrium temperature of $T_{\rm eq}=2492\pm38$ K and a measured dayside brightness temperature of $3237\pm59$ K (assuming zero geometric albedo), TOI-1518b is a promising candidate for future emission spectroscopy to probe for a thermal inversion. △ Less

Submitted 25 August, 2021; originally announced August 2021.

Comments: 25 pages, 11 figures, accepted to AJ

The obliquity and atmosphere of the ultra-hot Jupiter TOI-1431b (MASCARA-5b): A misaligned orbit and no signs of atomic ormolecular absorptions

Authors: M. Stangret, E. Pallé, N. Casasayas-Barris, M. Oshagh, A. Bello-Arufe, R. Luque, V. Nascimbeni, F. Yan, J. Orell-Miquel, D. Sicilia, L. Malavolta, B. C. Addison, L. A. Buchhave, A. S. Bonomo, F. Borsa, S. H. C. Cabot, M. Cecconi, D. A. Fischer, A. Harutyunyan, J. M. Mendonça, G. Nowak, H. Parviainen, A. Sozzetti, R. Tronsgaard

Abstract: Ultra-hot Jupiters are defined as giant planets with equilibrium temperatures larger than 2000 K. Most of them are found orbiting bright A-F type stars, making them extremely suitable objects to study their atmospheres using high-resolution spectroscopy. Recent studies show a variety of atoms and molecules detected in the atmospheres of this type of planets. Here we present our analysis of the new… ▽ More Ultra-hot Jupiters are defined as giant planets with equilibrium temperatures larger than 2000 K. Most of them are found orbiting bright A-F type stars, making them extremely suitable objects to study their atmospheres using high-resolution spectroscopy. Recent studies show a variety of atoms and molecules detected in the atmospheres of this type of planets. Here we present our analysis of the newly discovered ultra-hot Jupiter TOI-1431b/MASCARA-5b, using two transit observations with the HARPS-N spectrograph and one transit observation with the EXPRES spectrograph. Analysis of the Rossiter-McLaughlin effect shows that the planet is in a polar orbit, with a projected obliquity $ λ= -155^{+20}_{-10}$ degrees. Combining the nights and applying both cross-correlation methods and transmission spectroscopy, we find no evidences of CaI, FeI, FeII, MgI, NaI, VI, TiO, VO or H$α$ in the atmosphere of the planet. Our most likely explanation for the lack of atmospheric features is the large surface gravity of the planet. △ Less

Submitted 26 April, 2021; originally announced April 2021.

Comments: Accepted for publication in A&A (14 pages, 13 figures)

Journal ref: A&A 654, A73 (2021)

HST PanCET Program: A Complete Near-UV to Infrared Transmission Spectrum for the Hot Jupiter WASP-79b

Authors: Alexander D. Rathcke, Ryan J. MacDonald, Joanna K. Barstow, Jayesh M. Goyal, Mercedes Lopez-Morales, João M. Mendonça, Jorge Sanz-Forcada, Gregory W. Henry, David K. Sing, Munazza K. Alam, Nikole K. Lewis, Katy L. Chubb, Jake Taylor, Nikolay Nikolov, Lars A. Buchhave

Abstract: We present a new optical transmission spectrum of the hot Jupiter WASP-79b. We observed three transits with the STIS instrument mounted on HST, spanning 0.3 - 1.0 um. Combining these transits with previous observations, we construct a complete 0.3 - 5.0 um transmission spectrum of WASP-79b. Both HST and ground-based observations show decreasing transit depths towards blue wavelengths, contrary to… ▽ More We present a new optical transmission spectrum of the hot Jupiter WASP-79b. We observed three transits with the STIS instrument mounted on HST, spanning 0.3 - 1.0 um. Combining these transits with previous observations, we construct a complete 0.3 - 5.0 um transmission spectrum of WASP-79b. Both HST and ground-based observations show decreasing transit depths towards blue wavelengths, contrary to expectations from Rayleigh scattering or hazes. We infer atmospheric and stellar properties from the full near-UV to infrared transmission spectrum of WASP-79b using three independent retrieval codes, all of which yield consistent results. Our retrievals confirm previous detections of H$_{2}$O (at 4.0$σ$ confidence), while providing moderate evidence of H$^{-}$ bound-free opacity (3.3$σ$) and strong evidence of stellar contamination from unocculted faculae (4.7$σ$). The retrieved H$_{2}$O abundance ($\sim$ 1$\%$) suggests a super-stellar atmospheric metallicity, though stellar or sub-stellar abundances remain consistent with present observations (O/H = 0.3 - 34$\times$ stellar). All three retrieval codes obtain a precise H$^{-}$ abundance constraint: log(X$_{\rm{H^{-}}}$) $\approx$ -8.0 $\pm$ 0.7. The potential presence of H$^{-}$ suggests that JWST observations may be sensitive to ionic chemistry in the atmosphere of WASP-79b. The inferred faculae are $\sim$ 500 K hotter than the stellar photosphere, covering $\sim$ 15$\%$ of the stellar surface. Our analysis underscores the importance of observing UV - optical transmission spectra in order to disentangle the influence of unocculted stellar heterogeneities from planetary transmission spectra. △ Less

Submitted 21 April, 2021; originally announced April 2021.

Comments: 27 pages, 10 figures. Resubmitted to AJ after referee report

A survey of exoplanet phase curves with Ariel

Authors: Benjamin Charnay, Joao M. Mendonça, Laura Kreidberg, Nicolas B. Cowan, Jake Taylor, Taylor J. Bell, Olivier Demangeon, Billy Edwards, Carole A. Haswell, Giuseppe Morello, Lorenzo V. Mugnai, Enzo Pascale, Giovanna Tinetti, Pascal Tremblin, Robert T. Zellem

Abstract: The ESA-Ariel mission will include a tier dedicated to exoplanet phase curves corresponding to ~10% of the science time. We present here the current observing strategy for studying exoplanet phase curves with Ariel. We define science questions, requirements and a list of potential targets. We also estimate the precision of phase curve reconstruction and atmospheric retrieval using simulated phase… ▽ More The ESA-Ariel mission will include a tier dedicated to exoplanet phase curves corresponding to ~10% of the science time. We present here the current observing strategy for studying exoplanet phase curves with Ariel. We define science questions, requirements and a list of potential targets. We also estimate the precision of phase curve reconstruction and atmospheric retrieval using simulated phase curves. Based on this work, we found that full-orbit phase variations for 35-40 exoplanets could be observed during the 3.5-yr mission. This statistical sample would provide key constraints on atmospheric dynamics, composition, thermal structure and clouds of warm exoplanets, complementary to the scientific yield from spectroscopic transits/eclipses measurements. △ Less

Submitted 15 February, 2021; v1 submitted 10 February, 2021; originally announced February 2021.

Comments: 27 pages, 10 figures, accepted for publication in Experimental Astronomy, Ariel Special Issue

HELIOS-K 2.0 Opacity Calculator and Open-source Opacity Database for Exoplanetary Atmospheres

Authors: Simon L. Grimm, Matej Malik, Daniel Kitzmann, Andrea Guzmán-Mesa, H. Jens Hoeijmakers, Chloe Fisher, João M. Mendonça, Sergey N. Yurchenko, Jonathan Tennyson, Fabien Alesina, Nicolas Buchschacher, Julien Burnier, Damien Segransan, Robert L. Kurucz, Kevin Heng

Abstract: Computing and using opacities is a key part of modeling and interpreting data of exoplanetary atmospheres. Since the underlying spectroscopic line lists are constantly expanding and currently include up to ~ 10^10 - 10^11 transition lines, the opacity calculator codes need to become more powerful. Here we present major upgrades to the HELIOS-K GPU-accelerated opacity calculator and describe the ne… ▽ More Computing and using opacities is a key part of modeling and interpreting data of exoplanetary atmospheres. Since the underlying spectroscopic line lists are constantly expanding and currently include up to ~ 10^10 - 10^11 transition lines, the opacity calculator codes need to become more powerful. Here we present major upgrades to the HELIOS-K GPU-accelerated opacity calculator and describe the necessary steps to process large line lists within a reasonable amount of time. Besides performance improvements, we include more capabilities and present a toolbox for handling different atomic and molecular data sets: from downloading and pre-processing the data to performing the opacity calculations in a user-friendly way. HELIOS-K supports line lists from ExoMol, HITRAN, HITEMP, NIST, Kurucz and VALD3. By matching the resolution of 0.1 cm^-1 and cutting length of 25 cm^-1 used by the ExoCross code for timing performance (251 seconds excluding data read-in time), HELIOS-K can process the ExoMol BT2 water line list in 12.5 seconds. Using a resolution of 0.01 cm^-1, it takes 45 seconds - equivalent to about 10^7 lines per second. As a wavenumber resolution of 0.01 cm^-1 suffices for most exoplanetary atmosphere spectroscopic calculations, we adopt this resolution in calculating opacity functions for several hundred atomic and molecular species, and make them freely available on the open-access DACE database. For the opacity calculations of the database, we use a cutting length of 100 cm^-1 for molecules and no cutting length for atoms. Our opacities are available for downloading from https://dace.unige.ch/opacityDatabase and may be visualized using https://dace.unige.ch/opacity. △ Less

Submitted 22 March, 2021; v1 submitted 6 January, 2021; originally announced January 2021.

Comments: Published in The Astrophysical Journal Supplement Series

Modelling the 3D Climate of Venus with OASIS

Authors: João M. Mendonça, Lars A. Buchhave

Abstract: Flexible 3D models to explore the vast diversity of terrestrial planets and interpret observational data are still in their early stages. In this work, we present OASIS: a novel and flexible 3D virtual planet laboratory. With OASIS we envision a platform that couples self-consistently seven individual modules representing the main physical and chemical processes that shape planetary environments.… ▽ More Flexible 3D models to explore the vast diversity of terrestrial planets and interpret observational data are still in their early stages. In this work, we present OASIS: a novel and flexible 3D virtual planet laboratory. With OASIS we envision a platform that couples self-consistently seven individual modules representing the main physical and chemical processes that shape planetary environments. Additionally, OASIS is capable of producing simulated spectra from different instruments and observational techniques. In this work we focus on the benchmark test of coupling four of the physical modules: fluid dynamics, radiation, turbulence and surface/soil. To test the OASIS platform, we produced 3D simulations of the Venus climate and its atmospheric circulation and study how the modeled atmosphere changes with various cloud covers, atmospheric heat capacity, and surface friction. 3D simulations of Venus are challenging because they require long integration times with a computationally expensive radiative transfer code. By comparing OASIS results with observational data, we verify that the new model is able to successfully simulate Venus. With simulated spectra produced directly from the 3D simulations, we explore the capabilities of future missions, like LUVOIR, to observe Venus analogs located at a distance of 10 pc. With OASIS, we have taken the first steps to build a sophisticated and very flexible platform capable of studying the environment of terrestrial planets, which will be an essential tool to characterize observed terrestrial planets and plan future observations. △ Less

Submitted 17 July, 2020; v1 submitted 21 February, 2020; originally announced February 2020.

Comments: MNRAS published version

THOR 2.0: Major Improvements to the Open-Source General Circulation Model

Authors: Russell Deitrick, João M. Mendonça, Urs Schroffenegger, Simon L. Grimm, Shang-Min Tsai, Kevin Heng

Abstract: THOR is the first open-source general circulation model (GCM) developed from scratch to study the atmospheres and climates of exoplanets, free from Earth- or Solar System-centric tunings. It solves the general non-hydrostatic Euler equations (instead of the primitive equations) on a sphere using the icosahedral grid. In the current study, we report major upgrades to THOR, building upon the work of… ▽ More THOR is the first open-source general circulation model (GCM) developed from scratch to study the atmospheres and climates of exoplanets, free from Earth- or Solar System-centric tunings. It solves the general non-hydrostatic Euler equations (instead of the primitive equations) on a sphere using the icosahedral grid. In the current study, we report major upgrades to THOR, building upon the work of Mendonça et al. (2016). First, while the Horizontally Explicit Vertically Implicit (HEVI) integration scheme is the same as that described in Mendonça et al. (2016), we provide a clearer description of the scheme and improved its implementation in the code. The differences in implementation between the hydrostatic shallow (HSS), quasi-hydrostatic deep (QHD) and non-hydrostatic deep (NHD) treatments are fully detailed. Second, standard physics modules are added: two-stream, double-gray radiative transfer and dry convective adjustment. Third, THOR is tested on additional benchmarks: tidally-locked Earth, deep hot Jupiter, acoustic wave, and gravity wave. Fourth, we report that differences between the hydrostatic and non-hydrostatic simulations are negligible in the Earth case, but pronounced in the hot Jupiter case. Finally, the effects of the so-called "sponge layer", a form of drag implemented in most GCMs to provide numerical stability, are examined. Overall, these upgrades have improved the flexibility, user-friendliness, and stability of THOR. △ Less

Submitted 25 May, 2020; v1 submitted 29 November, 2019; originally announced November 2019.

Comments: 57 pages, 31 figures, revised, accepted for publication in ApJS

Angular Momentum and Heat Transport on Tidally Locked Hot Jupiter Planets

Authors: João M. Mendonça

Abstract: The atmospheric circulation in the upper atmosphere of hot Jupiter planets is strongly influenced by the incoming stellar radiation. In this work we explore the results from a 3D atmospheric model and revisit the main processes driving the circulation in hot Jupiter planets. We use the angular momentum transport as a diagnostic and carry out a Fourier analysis to identify the atmospheric waves inv… ▽ More The atmospheric circulation in the upper atmosphere of hot Jupiter planets is strongly influenced by the incoming stellar radiation. In this work we explore the results from a 3D atmospheric model and revisit the main processes driving the circulation in hot Jupiter planets. We use the angular momentum transport as a diagnostic and carry out a Fourier analysis to identify the atmospheric waves involved. We find that the coupling between the angular momentum transported horizontally by the semi-diurnal tide and the mean circulation is the mechanism responsible for producing the strong jet at low latitudes. Our simulations indicate the possible formation of atmospheric indirect cells at low latitudes. The formation of these cells is induced by the presence of the semi-diurnal tide that is driven by the stellar irradiation. The tropical circulation has an important impact transporting heat and momentum from the upper towards the lower atmosphere. One of the consequences of this heat and momentum transport is a global increase of the temperature. We show that the initial conditions do not affect the output of the reference simulation. However, when the period of rotation of the planet was increased ($P_{rot} > 5$ Earth days), vertical transport by stationary waves became stronger, transient waves became non-negligible, and Coriolis influence less dominant, which allowed a steady state with a strong retrograde jet to be stable. We found that at least two statically steady state solutions exist for the same planet parameters. △ Less

Submitted 23 October, 2019; originally announced October 2019.

Comments: Accepted for publication in MNRAS

Self-luminous and irradiated exoplanetary atmospheres explored with HELIOS

Authors: Matej Malik, Daniel Kitzmann, João M. Mendonça, Simon L. Grimm, Gabriel-Dominique Marleau, Esther F. Linder, Shang-Min Tsai, Kevin Heng

Abstract: We present new methodological features and physical ingredients included in the 1D radiative transfer code HELIOS, improving the hemispheric two-stream formalism. We conduct a thorough intercomparison survey with several established forward models, including COOLTLUSTY, PHOENIX, and find satisfactory consistency with their results. Then, we explore the impact of (i) different groups of opacity sou… ▽ More We present new methodological features and physical ingredients included in the 1D radiative transfer code HELIOS, improving the hemispheric two-stream formalism. We conduct a thorough intercomparison survey with several established forward models, including COOLTLUSTY, 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, the alkali atoms (and ionized hydrogen) playing an important role for the absorption of shortwave radiation (in very hot surroundings). Moreover, if these species are sufficiently abundant, they are likely to induce non-monotonic 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. △ Less

Submitted 15 March, 2019; originally announced March 2019.

Comments: 33 pages, 18 figures, accepted for publication in AJ

Three-Dimensional Circulation Driving Chemical Disequilibrium in WASP-43b

Authors: João M. Mendonça, Shang-Min Tsai, Matej Malik, Simon L. Grimm, Kevin Heng

Abstract: Spectral features in the observed spectra of exoplanets depend on the composition of their atmospheres. A good knowledge of the main atmospheric processes that drive the chemical distribution is therefore essential to interpret exoplanetary spectra. An atmosphere reaches chemical equilibrium if the rates of the forward and backward chemical reactions converge to the same value. However, there are… ▽ More Spectral features in the observed spectra of exoplanets depend on the composition of their atmospheres. A good knowledge of the main atmospheric processes that drive the chemical distribution is therefore essential to interpret exoplanetary spectra. An atmosphere reaches chemical equilibrium if the rates of the forward and backward chemical reactions converge to the same value. However, there are atmospheric processes, such as atmospheric transport, that destabilize this equilibrium. In this work we study the changes in composition driven by a 3D wind field in WASP-43b using our Global Circulation Model, THOR. Our model uses validated temperature- and pressure-dependent chemical timescales that allow us to explore the disequilibrium chemistry of CO, CO$_2$, H$_2$O and CH$_4$. In WASP-43b the formation of the equatorial jet has an important impact in the chemical distribution of the different species across the atmosphere. At low latitudes the chemistry is longitudinally quenched, except for CO$_2$ at solar abundances. The polar vortexes have a distinct chemical distribution since these are regions with lower temperature and atmospheric mixing. Vertical and latitudinal mixing have a secondary impact in the chemical transport. We determine graphically the effect of disequilibrium on observed emission spectra. Our results do not show any significant differences in the emission spectra between the equilibrium and disequilibrium solutions for C/O = 0.5. However, if C/O is increased to 2.0, differences in the spectra due to the disequilibrium chemistry of CH$_4$ become non-negligible. In some spectral ranges the emission spectra can have more than 15$\%$ departures from the equilibrium solution. △ Less

Submitted 30 October, 2018; v1 submitted 1 August, 2018; originally announced August 2018.

Comments: Accepted for publication in ApJ

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

Authors: João M. Mendonça, Matej Malik, Brice-Olivier Demory, Kevin Heng

Abstract: 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… ▽ More 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 cloudfree 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 cloudfree 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. △ Less

Submitted 8 February, 2018; originally announced February 2018.

Comments: Accepted for publication in ApJ

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

Authors: Baptiste Lavie, João M. Mendonça, Christoph Mordasini, Matej Malik, Mickaël Bonnefoy, Brice-Olivier Demory, Maria Oreshenko, Simon L. Grimm, David Ehrenreich, Kevin Heng

Abstract: We present an open-source retrieval code named HELIOS-Retrieval (hereafter HELIOS-R), designed to obtain chemical abundances and temperature-pressure profiles from inverting the measured spectra of exoplanetary atmospheres. In the current implementation, we use an exact solution of the radiative transfer equation, in the pure absorption limit, in our forward model, which allows us to analytically… ▽ More We present an open-source retrieval code named HELIOS-Retrieval (hereafter HELIOS-R), designed to obtain chemical abundances and temperature-pressure profiles from inverting the measured spectra of exoplanetary atmospheres. In the current implementation, we use an exact solution of the radiative transfer equation, in the pure absorption limit, in our forward model, which allows us to analytically integrate over all of the outgoing rays (instead of performing Gaussian quadrature). Two chemistry models are considered: unconstrained chemistry (where the mixing ratios are treated as free parameters) and equilibrium chemistry (enforced via analytical formulae, where only the elemental abundances are free parameters). 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 HR 8799b, c, d and e directly imaged exoplanets. Chemical equilibrium is disfavored by the Bayesian evidence for HR 8799b, c and d. We find supersolar C/O, C/H and O/H values for the outer HR 8799b and c exoplanets, while the inner HR 8799d and e exoplanets have substellar C/O, substellar C/H and superstellar O/H values. 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-R is publicly available as part of the Exoclimes Simulation Platform (ESP; www.exoclime.org). △ Less

Submitted 21 August, 2017; v1 submitted 11 October, 2016; originally announced October 2016.

Comments: 27 pages, 21 figures, 3 tables, published in AJ

Journal ref: AJ 154 91 (2017)

Exploring the Venus global super-rotation using a comprehensive General Circulation Model

Authors: João M. Mendonça, Peter L. Read

Abstract: The atmospheric circulation in Venus is well known to exhibit strong super-rotation. However, the atmospheric mechanisms responsible for the formation of this super-rotation are still not fully understood. In this work, we developed a new Venus general circulation model to study the most likely mechanisms driving the atmosphere to the current observed circulation. Our model includes a new radiativ… ▽ More The atmospheric circulation in Venus is well known to exhibit strong super-rotation. However, the atmospheric mechanisms responsible for the formation of this super-rotation are still not fully understood. In this work, we developed a new Venus general circulation model to study the most likely mechanisms driving the atmosphere to the current observed circulation. Our model includes a new radiative transfer, convection and suitably adapted boundary layer schemes and a dynamical core that takes into account the dependence of the heat capacity at constant pressure with temperature. The new Venus model is able to simulate a super-rotation phenomenon in the cloud region quantitatively similar to the one observed. The mechanisms maintaining the strong winds in the cloud region were found in the model results to be a combination of zonal mean circulation, thermal tides and transient waves. In this process, the semi-diurnal tide excited in the upper clouds has a key contribution in transporting axial angular momentum mainly from the upper atmosphere towards the cloud region. The magnitude of the super-rotation in the cloud region is sensitive to various radiative parameters such as the amount of solar radiative energy absorbed by the surface, which controls the static stability near the surface. In this work, we also discuss the main difficulties in representing the flow below the cloud base in Venus atmospheric models. Our new radiative scheme is more suitable for 3D Venus climate models than those used in previous work due to its easy adaptability to different atmospheric conditions. This flexibility of the model was crucial to explore the uncertainties in the lower atmospheric conditions and may also be used in the future to explore, for example, dynamical-radiative-microphysical feedbacks. △ Less

Submitted 21 September, 2016; originally announced September 2016.

Comments: Accepted for publication in Planet. Space Sci

THOR: A New and Flexible Global Circulation Model to Explore Planetary Atmospheres

Authors: João M. Mendonça, Simon L. Grimm, Luc Grosheintz, Kevin Heng

Abstract: We have designed and developed, from scratch, a global circulation model named THOR that solves the three-dimensional non-hydrostatic Euler equations. Our general approach lifts the commonly used assumptions of a shallow atmosphere and hydrostatic equilibrium. We solve the "pole problem" (where converging meridians on a sphere lead to increasingly smaller time steps near the poles) by implementing… ▽ More We have designed and developed, from scratch, a global circulation model named THOR that solves the three-dimensional non-hydrostatic Euler equations. Our general approach lifts the commonly used assumptions of a shallow atmosphere and hydrostatic equilibrium. We solve the "pole problem" (where converging meridians on a sphere lead to increasingly smaller time steps near the poles) by implementing an icosahedral grid. Irregularities in the grid, which lead to grid imprinting, are smoothed using the "spring dynamics" technique. We validate our implementation of spring dynamics by examining calculations of the divergence and gradient of test functions. To prevent the computational time step from being bottlenecked by having to resolve sound waves, we implement a split-explicit method together with a horizontally explicit and vertically implicit integration. We validate our global circulation model by reproducing the Earth and also the hot Jupiter-like benchmark tests. THOR was designed to run on Graphics Processing Units (GPUs), which allows for physics modules (radiative transfer, clouds, chemistry) to be added in the future, and is part of the open-source Exoclimes Simulation Platform (ESP; www.exoclime.org). △ Less

Submitted 19 July, 2016; originally announced July 2016.

Comments: Accepted for publication in ApJ

HELIOS: An Open-source, GPU-accelerated Radiative Transfer Code For Self-consistent Exoplanetary Atmospheres

Authors: Matej Malik, Luc Grosheintz, João M. Mendonça, Simon L. Grimm, Baptiste Lavie, Daniel Kitzmann, Shang-Min Tsai, Adam Burrows, Laura Kreidberg, Megan Bedell, Jacob L. Bean, Kevin B. Stevenson, Kevin Heng

Abstract: 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 non-isotropic scattering. A small set of the main infrared absorbers is employed, computed wi… ▽ More 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 non-isotropic 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 & 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 $\gtrsim 0.01$ cm$^{-1}$ resolution in the opacity function ($ \lesssim 10^3$ points per wavenumber bin) may result in errors $\gtrsim 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 element abundances) for 6 hot Jupiters. We find that the dayside emission spectra of HD 189733b and WASP-43b are consistent with the null hypothesis, while it consistently under-predicts 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 (ESP; exoclime.net). △ Less

Submitted 19 November, 2016; v1 submitted 17 June, 2016; originally announced June 2016.

Comments: 24 pages, 14 figures, accepted for publication in ApJ

Analytical Models of Exoplanetary Atmospheres. II. Radiative Transfer via the Two-stream Approximation

Authors: Kevin Heng, João M. Mendonça, Jae-Min Lee

Abstract: We present a comprehensive analytical study of radiative transfer using the method of moments and include the effects of non-isotropic scattering in the coherent limit. Within this unified formalism, we derive the governing equations and solutions describing two-stream radiative transfer (which approximates the passage of radiation as a pair of outgoing and incoming fluxes), flux-limited diffusion… ▽ More We present a comprehensive analytical study of radiative transfer using the method of moments and include the effects of non-isotropic scattering in the coherent limit. Within this unified formalism, we derive the governing equations and solutions describing two-stream radiative transfer (which approximates the passage of radiation as a pair of outgoing and incoming fluxes), flux-limited diffusion (which describes radiative transfer in the deep interior) and solutions for the temperature-pressure profiles. Generally, the problem is mathematically under-determined unless a set of closures (Eddington coefficients) is specified. We demonstrate that the hemispheric (or hemi-isotropic) closure naturally derives from the radiative transfer equation if energy conservation is obeyed, while the Eddington closure produces spurious enhancements of both reflected light and thermal emission. We concoct recipes for implementing two-stream radiative transfer in stand-alone numerical calculations and general circulation models. We use our two-stream solutions to construct toy models of the runaway greenhouse effect. We present a new solution for temperature-pressure profiles with a non-constant optical opacity and elucidate the effects of non-isotropic scattering in the optical and infrared. We derive generalized expressions for the spherical and Bond albedos and the photon deposition depth. We demonstrate that the value of the optical depth corresponding to the photosphere is not always 2/3 (Milne's solution) and depends on a combination of stellar irradiation, internal heat and the properties of scattering both in optical and infrared. Finally, we derive generalized expressions for the total, net, outgoing and incoming fluxes in the convective regime. △ Less

Submitted 27 October, 2014; v1 submitted 30 April, 2014; originally announced May 2014.

Comments: Accepted by ApJS. 23 pages, 11 figures, 3 tables, 158 equations. No change from previous version except for title (to match ApJS convention)