Showing 1–50 of 175 results for author: Lee, K H

Effects of the internal temperature on vertical mixing and on cloud structures in Ultra Hot Jupiters

Authors: Pascal A. Noti, Elspeth K. H. Lee

Abstract: The vertical mixing in hot Jupiter atmospheres plays a critical role in the formation and spacial distribution of cloud particles in their atmospheres. This affects the observed spectra of a planet through cloud opacity, which can be influenced by the degree of cold trapping of refractory species in the deep atmosphere. We aim to isolate the effects of the internal temperature on the mixing effici… ▽ More The vertical mixing in hot Jupiter atmospheres plays a critical role in the formation and spacial distribution of cloud particles in their atmospheres. This affects the observed spectra of a planet through cloud opacity, which can be influenced by the degree of cold trapping of refractory species in the deep atmosphere. We aim to isolate the effects of the internal temperature on the mixing efficiency in the atmospheres of Ultra Hot Jupiters (UHJ) and the spacial distribution of cloud particles across the globe. We couple a simplified tracer based cloud model, picket fence radiative-transfer scheme and mixing length theory to the Exo-FMS general circulation model. We run the model for five different internal temperatures at typical UHJ atmosphere system parameters. Our results show the convective eddy diffusion coefficient remains low throughout the vast majority of the atmosphere, with mixing dominated by advective flows. However, some regions can show convective mixing in the upper atmosphere for colder interior temperatures. The vertical extent of the clouds is reduced as the internal temperature is increased. Additionally, a global cloud layer gets formed below the radiative-convective boundary (RCB) in the cooler cases. Convection is generally strongly inhibited in UHJ atmospheres above the RCB due to their strong irradiation. Convective mixing plays a minor role compared to advective mixing in keeping cloud particles aloft in ultra hot Jupiters with warm interiors. Higher vertical turbulent heat fluxes and the advection of potential temperature inhibit convection in warmer interiors. Our results suggest isolated upper atmosphere regions above cold interiors may exhibit strong convective mixing in isolated regions around Rossby gyres, allowing aerosols to be better retained in these areas. △ Less

Submitted 25 September, 2024; originally announced September 2024.

Comments: accepted A&A (23 September 2024)

A photochemical PHO network for hydrogen-dominated exoplanet atmospheres

Authors: Elspeth K. H. Lee, Shang-Min Tsai, Julianne I. Moses, John M. C. Plane, Channon Visscher, Stephen J. Klippenstein

Abstract: Due to the detection of phosphine PH3 in the Solar System gas giants Jupiter and Saturn, PH3 has long been suggested to be detectable in exosolar substellar atmospheres too. However, to date, a direct detection of phosphine has proven to be elusive in exoplanet atmosphere surveys. We construct an updated phosphorus-hydrogen-oxygen (PHO) photochemical network suitable for simulation of gas giant hy… ▽ More Due to the detection of phosphine PH3 in the Solar System gas giants Jupiter and Saturn, PH3 has long been suggested to be detectable in exosolar substellar atmospheres too. However, to date, a direct detection of phosphine has proven to be elusive in exoplanet atmosphere surveys. We construct an updated phosphorus-hydrogen-oxygen (PHO) photochemical network suitable for simulation of gas giant hydrogen-dominated atmospheres. Using this network, we examine PHO photochemistry in hot Jupiter and warm Neptune exoplanet atmospheres at Solar and enriched metallicities. Our results show for HD 189733b-like hot Jupiters that HOPO, PO and P2 are typically the dominant P carriers at pressures important for transit and emission spectra, rather than PH3. For GJ1214b-like warm Neptune atmospheres our results suggest that at Solar metallicity PH3 is dominant in the absence of photochemistry, but is generally not in high abundance for all other chemical environments. At 10 and 100 times Solar, small oxygenated phosphorus molecules such as HOPO and PO dominate for both thermochemical and photochemical simulations. The network is able to reproduce well the observed PH3 abundances on Jupiter and Saturn. Despite progress in improving the accuracy of the PHO network, large portions of the reaction rate data remain with approximate, uncertain or missing values, which could change the conclusions of the current study significantly. Improving understanding of the kinetics of phosphorus-bearing chemical reactions will be a key undertaking for astronomers aiming to detect phosphine and other phosphorus species in both rocky and gaseous exoplanetary atmospheres in the near future. △ Less

Submitted 20 October, 2024; v1 submitted 10 September, 2024; originally announced September 2024.

Comments: Submitted to ApJ (12 July 2024) - Accepted ApJ (20 Oct 2024)

Multi-stream deep learning framework to predict mild cognitive impairment with Rey Complex Figure Test

Authors: Junyoung Park, Eun Hyun Seo, Sunjun Kim, SangHak Yi, Kun Ho Lee, Sungho Won

Abstract: Drawing tests like the Rey Complex Figure Test (RCFT) are widely used to assess cognitive functions such as visuospatial skills and memory, making them valuable tools for detecting mild cognitive impairment (MCI). Despite their utility, existing predictive models based on these tests often suffer from limitations like small sample sizes and lack of external validation, which undermine their reliab… ▽ More Drawing tests like the Rey Complex Figure Test (RCFT) are widely used to assess cognitive functions such as visuospatial skills and memory, making them valuable tools for detecting mild cognitive impairment (MCI). Despite their utility, existing predictive models based on these tests often suffer from limitations like small sample sizes and lack of external validation, which undermine their reliability. We developed a multi-stream deep learning framework that integrates two distinct processing streams: a multi-head self-attention based spatial stream using raw RCFT images and a scoring stream employing a previously developed automated scoring system. Our model was trained on data from 1,740 subjects in the Korean cohort and validated on an external hospital dataset of 222 subjects from Korea. The proposed multi-stream model demonstrated superior performance over baseline models (AUC = 0.872, Accuracy = 0.781) in external validation. The integration of both spatial and scoring streams enables the model to capture intricate visual details from the raw images while also incorporating structured scoring data, which together enhance its ability to detect subtle cognitive impairments. This dual approach not only improves predictive accuracy but also increases the robustness of the model, making it more reliable in diverse clinical settings. Our model has practical implications for clinical settings, where it could serve as a cost-effective tool for early MCI screening. △ Less

Submitted 4 September, 2024; originally announced September 2024.

Comments: 20 pages, 3 figures, 2 tables

Fast Private Location-based Information Retrieval Over the Torus

Authors: Joon Soo Yoo, Mi Yeon Hong, Ji Won Heo, Kang Hoon Lee, Ji Won Yoon

Abstract: Location-based services offer immense utility, but also pose significant privacy risks. In response, we propose LocPIR, a novel framework using homomorphic encryption (HE), specifically the TFHE scheme, to preserve user location privacy when retrieving data from public clouds. Our system employs TFHE's expertise in non-polynomial evaluations, crucial for comparison operations. LocPIR showcases min… ▽ More Location-based services offer immense utility, but also pose significant privacy risks. In response, we propose LocPIR, a novel framework using homomorphic encryption (HE), specifically the TFHE scheme, to preserve user location privacy when retrieving data from public clouds. Our system employs TFHE's expertise in non-polynomial evaluations, crucial for comparison operations. LocPIR showcases minimal client-server interaction, reduced memory overhead, and efficient throughput. Performance tests confirm its computational speed, making it a viable solution for practical scenarios, demonstrated via application to a COVID-19 alert model. Thus, LocPIR effectively addresses privacy concerns in location-based services, enabling secure data sharing from the public cloud. △ Less

Submitted 29 July, 2024; originally announced July 2024.

Comments: Accepted at the IEEE International Conference on Advanced Video and Signal-Based Surveillance (AVSS) 2024

On selection of centroids of fuzzy clusters for color classification

Authors: Dae-Won Kim, Kwang H. Lee

Abstract: A novel initialization method in the fuzzy c-means (FCM) algorithm is proposed for the color clustering problem. Given a set of color points, the proposed initialization extracts dominant colors that are the most vivid and distinguishable colors. Color points closest to the dominant colors are selected as initial centroids in the FCM. To obtain the dominant colors and their closest color points, w… ▽ More A novel initialization method in the fuzzy c-means (FCM) algorithm is proposed for the color clustering problem. Given a set of color points, the proposed initialization extracts dominant colors that are the most vivid and distinguishable colors. Color points closest to the dominant colors are selected as initial centroids in the FCM. To obtain the dominant colors and their closest color points, we introduce reference colors and define a fuzzy membership model between a color point and a reference color. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Sulphur dioxide in the mid-infrared transmission spectrum of WASP-39b

Authors: Diana Powell, Adina D. Feinstein, Elspeth K. H. Lee, Michael Zhang, Shang-Min Tsai, Jake Taylor, James Kirk, Taylor Bell, Joanna K. Barstow, Peter Gao, Jacob L. Bean, Jasmina Blecic, Katy L. Chubb, Ian J. M. Crossfield, Sean Jordan, Daniel Kitzmann, Sarah E. Moran, Giuseppe Morello, Julianne I. Moses, Luis Welbanks, Jeehyun Yang, Xi Zhang, Eva-Maria Ahrer, Aaron Bello-Arufe, Jonathan Brande , et al. (48 additional authors not shown)

Abstract: The recent inference of sulphur dioxide (SO$_2$) in the atmosphere of the hot ($\sim$1100 K), Saturn-mass exoplanet WASP-39b from near-infrared JWST observations suggests that photochemistry is a key process in high temperature exoplanet atmospheres. This is due to the low ($<$1 ppb) abundance of SO$_2$ under thermochemical equilibrium, compared to that produced from the photochemistry of H$_2$O a… ▽ More The recent inference of sulphur dioxide (SO$_2$) in the atmosphere of the hot ($\sim$1100 K), Saturn-mass exoplanet WASP-39b from near-infrared JWST observations suggests that photochemistry is a key process in high temperature exoplanet atmospheres. This is due to the low ($<$1 ppb) abundance of SO$_2$ under thermochemical equilibrium, compared to that produced from the photochemistry of H$_2$O and H$_2$S (1-10 ppm). However, the SO$_2$ inference was made from a single, small molecular feature in the transmission spectrum of WASP-39b at 4.05 $μ$m, and therefore the detection of other SO$_2$ absorption bands at different wavelengths is needed to better constrain the SO$_2$ abundance. Here we report the detection of SO$_2$ spectral features at 7.7 and 8.5 $μ$m in the 5-12 $μ$m transmission spectrum of WASP-39b measured by the JWST Mid-Infrared Instrument (MIRI) Low Resolution Spectrometer (LRS). Our observations suggest an abundance of SO$_2$ of 0.5-25 ppm (1$σ$ range), consistent with previous findings. In addition to SO$_2$, we find broad water vapour absorption features, as well as an unexplained decrease in the transit depth at wavelengths longer than 10 $μ$m. Fitting the spectrum with a grid of atmospheric forward models, we derive an atmospheric heavy element content (metallicity) for WASP-39b of $\sim$7.1-8.0 $\times$ solar and demonstrate that photochemistry shapes the spectra of WASP-39b across a broad wavelength range. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Comments: Published in Nature

Journal ref: Nature 626, 979-983 (2024)

Fuzzy color model and clustering algorithm for color clustering problem

Authors: Dae-Won Kim, Kwang H. Lee

Abstract: The research interest of this paper is focused on the efficient clustering task for an arbitrary color data. In order to tackle this problem, we have tried to model the inherent uncertainty and vagueness of color data using fuzzy color model. By taking fuzzy approach to color modeling, we could make a soft decision for the vague regions between neighboring colors. The proposed fuzzy color model de… ▽ More The research interest of this paper is focused on the efficient clustering task for an arbitrary color data. In order to tackle this problem, we have tried to model the inherent uncertainty and vagueness of color data using fuzzy color model. By taking fuzzy approach to color modeling, we could make a soft decision for the vague regions between neighboring colors. The proposed fuzzy color model defined a three dimensional fuzzy color ball and color membership computation method with two inter-color distances. With the fuzzy color model, we developed a new fuzzy clustering algorithm for an efficient partition of color data. Each fuzzy cluster set has a cluster prototype which is represented by fuzzy color centroid. △ Less

Submitted 9 July, 2024; originally announced July 2024.

A new validity measure for fuzzy c-means clustering

Authors: Dae-Won Kim, Kwang H. Lee

Abstract: A new cluster validity index is proposed for fuzzy clusters obtained from fuzzy c-means algorithm. The proposed validity index exploits inter-cluster proximity between fuzzy clusters. Inter-cluster proximity is used to measure the degree of overlap between clusters. A low proximity value refers to well-partitioned clusters. The best fuzzy c-partition is obtained by minimizing inter-cluster proximi… ▽ More A new cluster validity index is proposed for fuzzy clusters obtained from fuzzy c-means algorithm. The proposed validity index exploits inter-cluster proximity between fuzzy clusters. Inter-cluster proximity is used to measure the degree of overlap between clusters. A low proximity value refers to well-partitioned clusters. The best fuzzy c-partition is obtained by minimizing inter-cluster proximity with respect to c. Well-known data sets are tested to show the effectiveness and reliability of the proposed index. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: Accepted at FIP-2002

Isotropy of cosmic rays beyond $10^{20}$ eV favors their heavy mass composition

Authors: Telescope Array Collaboration, R. U. Abbasi, Y. Abe, T. Abu-Zayyad, M. Allen, Y. Arai, R. Arimura, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, I. Buckland, B. G. Cheon, M. Chikawa, T. Fujii, K. Fujisue, K. Fujita, R. Fujiwara, M. Fukushima, G. Furlich, N. Globus, R. Gonzalez, W. Hanlon, N. Hayashida, H. He , et al. (118 additional authors not shown)

Abstract: We report an estimation of the injected mass composition of ultra-high energy cosmic rays (UHECRs) at energies higher than 10 EeV. The composition is inferred from an energy-dependent sky distribution of UHECR events observed by the Telescope Array surface detector by comparing it to the Large Scale Structure of the local Universe. In the case of negligible extra-galactic magnetic fields the resul… ▽ More We report an estimation of the injected mass composition of ultra-high energy cosmic rays (UHECRs) at energies higher than 10 EeV. The composition is inferred from an energy-dependent sky distribution of UHECR events observed by the Telescope Array surface detector by comparing it to the Large Scale Structure of the local Universe. In the case of negligible extra-galactic magnetic fields the results are consistent with a relatively heavy injected composition at E ~ 10 EeV that becomes lighter up to E ~ 100 EeV, while the composition at E > 100 EeV is very heavy. The latter is true even in the presence of highest experimentally allowed extra-galactic magnetic fields, while the composition at lower energies can be light if a strong EGMF is present. The effect of the uncertainty in the galactic magnetic field on these results is subdominant. △ Less

Submitted 3 July, 2024; v1 submitted 27 June, 2024; originally announced June 2024.

Comments: 8 pages, 3 figures, accepted for publication in PRL

Mass composition of ultra-high energy cosmic rays from distribution of their arrival directions with the Telescope Array

Authors: Telescope Array Collaboration, R. U. Abbasi, Y. Abe, T. Abu-Zayyad, M. Allen, Y. Arai, R. Arimura, E. Barcikowski, J. W. Belz, D. R. Bergman, S. A. Blake, I. Buckland, B. G. Cheon, M. Chikawa, T. Fujii, K. Fujisue, K. Fujita, R. Fujiwara, M. Fukushima, G. Furlich, N. Globus, R. Gonzalez, W. Hanlon, N. Hayashida, H. He , et al. (118 additional authors not shown)

Abstract: We use a new method to estimate the injected mass composition of ultrahigh cosmic rays (UHECRs) at energies higher than 10 EeV. The method is based on comparison of the energy-dependent distribution of cosmic ray arrival directions as measured by the Telescope Array experiment (TA) with that calculated in a given putative model of UHECR under the assumption that sources trace the large-scale struc… ▽ More We use a new method to estimate the injected mass composition of ultrahigh cosmic rays (UHECRs) at energies higher than 10 EeV. The method is based on comparison of the energy-dependent distribution of cosmic ray arrival directions as measured by the Telescope Array experiment (TA) with that calculated in a given putative model of UHECR under the assumption that sources trace the large-scale structure (LSS) of the Universe. As we report in the companion letter, the TA data show large deflections with respect to the LSS which can be explained, assuming small extra-galactic magnetic fields (EGMF), by an intermediate composition changing to a heavy one (iron) in the highest energy bin. Here we show that these results are robust to uncertainties in UHECR injection spectra, the energy scale of the experiment and galactic magnetic fields (GMF). The assumption of weak EGMF, however, strongly affects this interpretation at all but the highest energies E > 100 EeV, where the remarkable isotropy of the data implies a heavy injected composition even in the case of strong EGMF. This result also holds if UHECR sources are as rare as $2 \times 10^{-5}$ Mpc$^{-3}$, that is the conservative lower limit for the source number density. △ Less

Submitted 3 July, 2024; v1 submitted 27 June, 2024; originally announced June 2024.

Comments: 18 pages, 11 figures, accepted for publication in PRD

Latitudinal Asymmetry in the Dayside Atmosphere of WASP-43b

Authors: Ryan C. Challener, Zafar Rustamkulov, Elspeth K. H. Lee, Nikole Lewis, David K. Sing, Stephan M. Birkmann, Nicolas Crouzet, Néstor Espinoza, Elena Manjavacas, Natalia Oliveros-Gomez, Jeff A. Valenti, Jingxuan Yang

Abstract: We present two-dimensional near-infrared temperature maps of the canonical hot Jupiter WASP-43b using a phase-curve observation with JWST NIRSpec/G395H. From the white-light planetary transit, we improve constraints on the planet's orbital parameters and measure a planet-to-star radius ratio of $0.15883^{+0.00056}_{-0.00053}$. Using the white-light phase curve, we measure a longitude of maximum br… ▽ More We present two-dimensional near-infrared temperature maps of the canonical hot Jupiter WASP-43b using a phase-curve observation with JWST NIRSpec/G395H. From the white-light planetary transit, we improve constraints on the planet's orbital parameters and measure a planet-to-star radius ratio of $0.15883^{+0.00056}_{-0.00053}$. Using the white-light phase curve, we measure a longitude of maximum brightness of $6.9^{+0^\circ.5}_{-0^\circ.5}$ east of the substellar point and a phase-curve offset of $10.0^{+0^\circ.8}_{-0^\circ.8}$. We also find an $\approx4σ$ detection of a latitudinal hotspot offset of $-13.4^{+3^\circ.2}_{-1^\circ.7}$, the first significant detection of a non-equatorial hotspot in an exoplanet atmosphere. We show that this detection is robust to variations within planetary parameter uncertainties, but only if the transit is used to improve constraints, showing the importance of transit observations to eclipse mapping. Maps retrieved from the NRS1 and NRS2 detectors are similar, with hotspot locations consistent between the two detectors at the $1σ$ level. Our JWST data show brighter (hotter) nightsides and a dimmer (colder) dayside at the shorter wavelengths relative to fits to \textit{Spitzer} 3.6 and 4.5 \microns\ phase curves. Through comparison between our phase curves and a set of general circulation models, we find evidence for clouds on the nightside and atmospheric drag or high metallicity reducing the eastward hotspot offset. △ Less

Submitted 14 June, 2024; originally announced June 2024.

Comments: 15 pages, 8 figures. In review at ApJ Letters

Phase-resolving the absorption signatures of water and carbon monoxide in the atmosphere of the ultra-hot Jupiter WASP-121b with GEMINI-S/IGRINS

Authors: Joost P. Wardenier, Vivien Parmentier, Michael R. Line, Megan Weiner Mansfield, Xianyu Tan, Shang-Min Tsai, Jacob L. Bean, Jayne L. Birkby, Matteo Brogi, Jean-Michel Désert, Siddharth Gandhi, Elspeth K. H. Lee, Colette I. Levens, Lorenzo Pino, Peter C. B. Smith

Abstract: Ultra-hot Jupiters are among the best targets for atmospheric characterization at high spectral resolution. Resolving their transmission spectra as a function of orbital phase offers a unique window into the 3D nature of these objects. In this work, we present three transits of the ultra-hot Jupiter WASP-121b observed with Gemini-S/IGRINS. For the first time, we measure the phase-dependent absorpt… ▽ More Ultra-hot Jupiters are among the best targets for atmospheric characterization at high spectral resolution. Resolving their transmission spectra as a function of orbital phase offers a unique window into the 3D nature of these objects. In this work, we present three transits of the ultra-hot Jupiter WASP-121b observed with Gemini-S/IGRINS. For the first time, we measure the phase-dependent absorption signals of CO and H$_{\text{2}}$O in the atmosphere of an exoplanet, and we find that they are different. While the blueshift of CO increases during the transit, the absorption lines of H$_{\text{2}}$O become less blueshifted with phase, and even show a redshift in the second half of the transit. These measurements reveal the distinct spatial distributions of both molecules across the atmospheres of ultra-hot Jupiters. Also, we find that the H$_{\text{2}}$O signal is absent in the first quarter of the transit, potentially hinting at cloud formation on the evening terminator of WASP-121b. To further interpret the absorption trails of CO and H$_{\text{2}}$O, as well as the Doppler shifts of Fe previously measured with VLT/ESPRESSO, we compare the data to simulated transits of WASP-121b. To this end, we post-processes the outputs of global circulation models with a 3D Monte-Carlo radiative transfer code. Our analysis shows that the atmosphere of WASP-121b is subject to atmospheric drag, as previously suggested by small hotspot offsets inferred from phase-curve observations. Our study highlights the importance of phase-resolved spectroscopy in unravelling the complex atmospheric structure of ultra-hot Jupiters and sets the stage for further investigations into their chemistry and dynamics. △ Less

Submitted 18 July, 2024; v1 submitted 13 June, 2024; originally announced June 2024.

Comments: 24 pages, 16 figures, accepted for publication in PASP

Observation of Declination Dependence in the Cosmic Ray Energy Spectrum

Authors: The Telescope Array Collaboration, R. U. Abbasi, T. Abu-Zayyad, M. Allen, J. W. Belz, D. R. Bergman, I. Buckland, W. Campbell, B. G. Cheon, K. Endo, A. Fedynitch, T. Fujii, K. Fujisue, K. Fujita, M. Fukushima, G. Furlich, Z. Gerber, N. Globus, W. Hanlon, N. Hayashida, H. He, K. Hibino, R. Higuchi, D. Ikeda, T. Ishii , et al. (101 additional authors not shown)

Abstract: We report on an observation of the difference between northern and southern skies of the ultrahigh energy cosmic ray energy spectrum with a significance of ${\sim}8σ$. We use measurements from the two largest experiments$\unicode{x2014}$the Telescope Array observing the northern hemisphere and the Pierre Auger Observatory viewing the southern hemisphere. Since the comparison of two measurements fr… ▽ More We report on an observation of the difference between northern and southern skies of the ultrahigh energy cosmic ray energy spectrum with a significance of ${\sim}8σ$. We use measurements from the two largest experiments$\unicode{x2014}$the Telescope Array observing the northern hemisphere and the Pierre Auger Observatory viewing the southern hemisphere. Since the comparison of two measurements from different observatories introduces the issue of possible systematic differences between detectors and analyses, we validate the methodology of the comparison by examining the region of the sky where the apertures of the two observatories overlap. Although the spectra differ in this region, we find that there is only a $1.8σ$ difference between the spectrum measurements when anisotropic regions are removed and a fiducial cut in the aperture is applied. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: 8 pages, 6 figures

A warm Neptune's methane reveals core mass and vigorous atmospheric mixing

Authors: David K. Sing, Zafar Rustamkulov, Daniel P. Thorngren, Joanna K. Barstow, Pascal Tremblin, Catarina Alves de Oliveira, Tracy L. Beck, Stephan M. Birkmann, Ryan C. Challener, Nicolas Crouzet, Néstor Espinoza, Pierre Ferruit, Giovanna Giardino, Amélie Gressier, Elspeth K. H. Lee, Nikole K. Lewis, Roberto Maiolino, Elena Manjavacas, Bernard J. Rauscher, Marco Sirianni, Jeff A. Valenti

Abstract: Observations of transiting gas giant exoplanets have revealed a pervasive depletion of methane, which has only recently been identified atmospherically. The depletion is thought to be maintained by disequilibrium processes such as photochemistry or mixing from a hotter interior. However, the interiors are largely unconstrained along with the vertical mixing strength and only upper limits on the CH… ▽ More Observations of transiting gas giant exoplanets have revealed a pervasive depletion of methane, which has only recently been identified atmospherically. The depletion is thought to be maintained by disequilibrium processes such as photochemistry or mixing from a hotter interior. However, the interiors are largely unconstrained along with the vertical mixing strength and only upper limits on the CH$_4$ depletion have been available. The warm Neptune WASP-107 b stands out among exoplanets with an unusually low density, reported low core mass, and temperatures amenable to CH$_4$ though previous observations have yet to find the molecule. Here we present a JWST NIRSpec transmission spectrum of WASP-107 b which shows features from both SO$_2$ and CH$_4$ along with H$_2$O, CO$_2$, and CO. We detect methane with 4.2$σ$ significance at an abundance of 1.0$\pm$0.5 ppm, which is depleted by 3 orders of magnitude relative to equilibrium expectations. Our results are highly constraining for the atmosphere and interior, which indicate the envelope has a super-solar metallicity of 43$\pm$8$\times$ solar, a hot interior with an intrinsic temperature of T$_{\rm int}$=460$\pm$40 K, and vigorous vertical mixing which depletes CH4 with a diffusion coefficient of Kzz = 10$^{11.6\pm0.1}$ cm$^2$/s. Photochemistry has a negligible effect on the CH$_4$ abundance, but is needed to account for the SO$_2$. We infer a core mass of 11.5$_{-3.6}^{+3.0}$ M$_{\odot}$, which is much higher than previous upper limits, releasing a tension with core-accretion models. △ Less

Submitted 17 May, 2024; originally announced May 2024.

Comments: Published in Nature at this URL, https://www.nature.com/articles/s41586-024-07395-z . This is the authors version of the manuscript, 20 pages including Methods. Data from Figs. 1, 2, and 3 are available at this URL, https://doi.org/10.5281/zenodo.10891400

An atlas of resolved spectral features in the transmission spectrum of WASP-189 b with MAROON-X

Authors: B. Prinoth, H. J. Hoeijmakers, B. M. Morris, M. Lam, D. Kitzmann, E. Sedaghati, J. V. Seidel, E. K. H. Lee, B. Thorsbro, N. W. Borsato, Y. C. Damasceno, S. Pelletier, A. Seifahrt

Abstract: Exoplanets in the ultra-hot Jupiter regime provide an excellent laboratory for testing the impact of stellar irradiation on the dynamics and chemical composition of gas giant atmospheres. In this study, we observed two transits of the ultra-hot Jupiter WASP-189 b with MAROON-X/Gemini-North to probe its high-altitude atmospheric layers, using strong absorption lines. We derived posterior probabilit… ▽ More Exoplanets in the ultra-hot Jupiter regime provide an excellent laboratory for testing the impact of stellar irradiation on the dynamics and chemical composition of gas giant atmospheres. In this study, we observed two transits of the ultra-hot Jupiter WASP-189 b with MAROON-X/Gemini-North to probe its high-altitude atmospheric layers, using strong absorption lines. We derived posterior probability distributions for the planetary and stellar parameters by calculating the stellar spectrum behind the planet at every orbital phase during the transit. This was used to correct the Rossiter-McLaughlin imprint on the transmission spectra. Using differential transmission spectroscopy, we detect strong absorption lines of Ca+, Ba+, Na, H$α$, Mg, Fe, and Fe+, providing an unprecedented and detailed view of the atmospheric chemical composition. Ca+ absorption is particularly well suited for analysis through time-resolved narrow-band spectroscopy, owing to its transition lines formed in high-altitude layers. The spectral absorption lines show no significant blueshifts that would indicate high-altitude day-to-night winds, and further analysis is needed to investigate the implications for atmospheric dynamics. These high signal-to-noise observations provide a benchmark data set for testing high-resolution retrievals and the assumptions of atmospheric models. We also simulate observations of WASP-189 b with ANDES/ELT, and show that ANDES will be highly sensitive to the individual absorption lines of a myriad of elements and molecules, including TiO and CO. △ Less

Submitted 13 March, 2024; originally announced March 2024.

Comments: 34 pages, 31 figures, accepted for publication in A&A on 16 February 2024

Testing approximate infrared scattering radiative-transfer methods for hot Jupiter atmospheres

Authors: Elspeth K. H. Lee

Abstract: The calculation of internal atmospheric (longwave) fluxes is a key component of any model of exoplanet atmospheres that requires radiative-transfer (RT) calculations. For atmospheres containing a strong scattering component such as cloud particles, most 1D multiple-scattering RT methods typically involve numerically expensive matrix inversions. This computational bottleneck is exacerbated when mul… ▽ More The calculation of internal atmospheric (longwave) fluxes is a key component of any model of exoplanet atmospheres that requires radiative-transfer (RT) calculations. For atmospheres containing a strong scattering component such as cloud particles, most 1D multiple-scattering RT methods typically involve numerically expensive matrix inversions. This computational bottleneck is exacerbated when multitudes of RT calculations are required, such as in general circulation models (GCMs) and retrieval methods. In an effort to increase the speed of RT calculations without sacrificing too much accuracy, we investigate the applicability of approximate longwave scattering methods developed for the Earth science community to hot Jupiter atmospheres. We test the absorption approximation (AA) and variational iteration method (VIM) applied to typical cloudy hot Jupiter scenarios, using 64 stream DISORT calculations as reference solutions. We find the four-stream VIM variant is a highly promising method to explore using for hot Jupiter GCM and retrieval modelling, showing excellent speed characteristics, with typical errors $\sim$1\% for outgoing fluxes and within $\sim$50\%, but with larger errors in the deep cloud layer test case, for vertical heating rates. Other methods explored in this study were found to typically produce similar error characteristics in vertical heating rates. △ Less

Submitted 27 February, 2024; originally announced February 2024.

Comments: Accepted ApJ (27th Feb 2024)

Threads of Subtlety: Detecting Machine-Generated Texts Through Discourse Motifs

Authors: Zae Myung Kim, Kwang Hee Lee, Preston Zhu, Vipul Raheja, Dongyeop Kang

Abstract: With the advent of large language models (LLM), the line between human-crafted and machine-generated texts has become increasingly blurred. This paper delves into the inquiry of identifying discernible and unique linguistic properties in texts that were written by humans, particularly uncovering the underlying discourse structures of texts beyond their surface structures. Introducing a novel metho… ▽ More With the advent of large language models (LLM), the line between human-crafted and machine-generated texts has become increasingly blurred. This paper delves into the inquiry of identifying discernible and unique linguistic properties in texts that were written by humans, particularly uncovering the underlying discourse structures of texts beyond their surface structures. Introducing a novel methodology, we leverage hierarchical parse trees and recursive hypergraphs to unveil distinctive discourse patterns in texts produced by both LLMs and humans. Empirical findings demonstrate that, although both LLMs and humans generate distinct discourse patterns influenced by specific domains, human-written texts exhibit more structural variability, reflecting the nuanced nature of human writing in different domains. Notably, incorporating hierarchical discourse features enhances binary classifiers' overall performance in distinguishing between human-written and machine-generated texts, even on out-of-distribution and paraphrased samples. This underscores the significance of incorporating hierarchical discourse features in the analysis of text patterns. The code and dataset are available at https://github.com/minnesotanlp/threads-of-subtlety. △ Less

Submitted 6 June, 2024; v1 submitted 16 February, 2024; originally announced February 2024.

Comments: 26 pages, accepted at ACL 2024 (Main)

Modeling Atmospheric Lines By the Exoplanet Community (MALBEC) version 1.0: A CUISINES radiative transfer intercomparison project

Authors: Geronimo L. Villanueva, Thomas J. Fauchez, Vincent Kofman, Eleonora Alei, Elspeth K. H. Lee, Estelle Janin, Michael D. Himes, Jeremy Leconte, Michaela Leung, Sara Faggi, Mei Ting Mak, Denis E. Sergeev, Thea Kozakis, James Manners, Nathan Mayne, Edward W. Schwieterman, Alex R. Howe, Natasha Batalha

Abstract: Radiative transfer (RT) models are critical in the interpretation of exoplanetary spectra, in simulating exoplanet climates and when designing the specifications of future flagship observatories. However, most models differ in methodologies and input data, which can lead to significantly different spectra. In this paper, we present the experimental protocol of the MALBEC (Modeling Atmospheric Line… ▽ More Radiative transfer (RT) models are critical in the interpretation of exoplanetary spectra, in simulating exoplanet climates and when designing the specifications of future flagship observatories. However, most models differ in methodologies and input data, which can lead to significantly different spectra. In this paper, we present the experimental protocol of the MALBEC (Modeling Atmospheric Lines By the Exoplanet Community) project. MALBEC is an exoplanet model intercomparison project (exoMIP) that belongs to the CUISINES (Climates Using Interactive Suites of Intercomparisons Nested for Exoplanet Studies) framework which aims to provide the exoplanet community with a large and diverse set of comparison and validation of models. The proposed protocol tests include a large set of initial participating RT models, a broad range of atmospheres (from Hot Jupiters to temperate terrestrials) and several observation geometries, which would allow us to quantify and compare the differences between different RT models used by the exoplanetary community. Two types of tests are proposed: transit spectroscopy and direct imaging modeling, with results from the proposed tests to be published in dedicated follow-up papers. To encourage the community to join this comparison effort and as an example, we present simulation results for one specific transit case (GJ-1214 b), in which we find notable differences in how the various codes handle the discretization of the atmospheres (e.g., sub-layering), the treatment of molecular opacities (e.g., correlated-k, line-by-line) and the default spectroscopic repositories generally used by each model (e.g., HITRAN, HITEMP, ExoMol). △ Less

Submitted 6 February, 2024; originally announced February 2024.

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

Group lasso priors for Bayesian accelerated failure time models with left-truncated and interval-censored data

Authors: Harrison T. Reeder, Sebastien Haneuse, Kyu Ha Lee

Abstract: An important task in health research is to characterize time-to-event outcomes such as disease onset or mortality in terms of a potentially high-dimensional set of risk factors. For example, prospective cohort studies of Alzheimer's disease typically enroll older adults for observation over several decades to assess the long-term impact of genetic and other factors on cognitive decline and mortali… ▽ More An important task in health research is to characterize time-to-event outcomes such as disease onset or mortality in terms of a potentially high-dimensional set of risk factors. For example, prospective cohort studies of Alzheimer's disease typically enroll older adults for observation over several decades to assess the long-term impact of genetic and other factors on cognitive decline and mortality. The accelerated failure time model is particularly well-suited to such studies, structuring covariate effects as `horizontal' changes to the survival quantiles that conceptually reflect shifts in the outcome distribution due to lifelong exposures. However, this modeling task is complicated by the enrollment of adults at differing ages, and intermittent followup visits leading to interval censored outcome information. Moreover, genetic and clinical risk factors are not only high-dimensional, but characterized by underlying grouping structure, such as by function or gene location. Such grouped high-dimensional covariates require shrinkage methods that directly acknowledge this structure to facilitate variable selection and estimation. In this paper, we address these considerations directly by proposing a Bayesian accelerated failure time model with a group-structured lasso penalty, designed for left-truncated and interval-censored time-to-event data. We develop a custom Markov chain Monte Carlo sampler for efficient estimation, and investigate the impact of various methods of penalty tuning and thresholding for variable selection. We present a simulation study examining the performance of this method relative to models with an ordinary lasso penalty, and apply the proposed method to identify groups of predictive genetic and clinical risk factors for Alzheimer's disease in the Religious Orders Study and Memory and Aging Project (ROSMAP) prospective cohort studies of AD and dementia. △ Less

Submitted 11 January, 2024; v1 submitted 9 January, 2024; originally announced January 2024.

Dynamically coupled kinetic chemistry in brown dwarf atmospheres -- II. Cloud and chemistry connections in directly imaged sub-Jupiter exoplanets

Authors: Elspeth K. H. Lee, Xianyu Tan, Shang-Min Tsai

Abstract: With JWST slated to gain high fidelity time dependent data on brown dwarf atmospheres, it is highly anticipated to do the same for directly imaged, sub-Jupiter exoplanets. With this new capability, the need for a full 3D understanding to explain spectral features and their time dependence is becoming a vital aspect for consideration. To examine the atmospheric properties of directly imaged sub-Jup… ▽ More With JWST slated to gain high fidelity time dependent data on brown dwarf atmospheres, it is highly anticipated to do the same for directly imaged, sub-Jupiter exoplanets. With this new capability, the need for a full 3D understanding to explain spectral features and their time dependence is becoming a vital aspect for consideration. To examine the atmospheric properties of directly imaged sub-Jupiter exoplanets, we use the three dimensional Exo-FMS general circulation model (GCM) to simulate a metal enhanced generic young sub-Jupiter object. We couple Exo-FMS to a kinetic chemistry scheme, a tracer based cloud formation scheme and a spectral radiative-transfer model to take into account the chemical and cloud feedback on the atmospheric thermochemical and dynamical properties. Our results show a highly complex feedback between clouds and chemistry onto the 3D temperature structure of the atmosphere, bringing about latitudinal differences and inducing time-dependent stormy features at photospheric pressures. This suggests a strong connection and feedback between the spatial cloud coverage and chemical composition of the atmosphere, with the temperature changes and dynamical motions induced by cloud opacity and triggered convection feedback driving chemical species behaviour. In addition, we also produce synthetic latitude dependent and time dependent spectra of our model to investigate atmospheric variability and periodicity in commonly used photometric bands. Overall, our efforts put the included physics in 3D simulations of exoplanets on par with contemporary 1D radiative-convective equilibrium modelling. △ Less

Submitted 15 February, 2024; v1 submitted 28 November, 2023; originally announced November 2023.

Comments: Accepted MNRAS (Feb 2024), 17 pages

A Direct Comparison between the use of Double Gray and Multiwavelength Radiative Transfer in a General Circulation Model with and without Radiatively Active Clouds

Authors: Isaac Malsky, Emily Rauscher, Michael T. Roman, Elspeth K. H. Lee, Hayley Beltz, Arjun Savel, Eliza M. R. Kempton, L. Cinque

Abstract: Inhomogeneous cloud formation and wavelength-dependent phenomena are expected to shape hot Jupiter atmospheres. We present a General Circulation Model (GCM) with multiwavelength "picket fence" radiative transfer and radiatively active, temperature dependent clouds, and compare the results to a double gray routine. The double gray method inherently fails to model polychromatic effects in hot Jupite… ▽ More Inhomogeneous cloud formation and wavelength-dependent phenomena are expected to shape hot Jupiter atmospheres. We present a General Circulation Model (GCM) with multiwavelength "picket fence" radiative transfer and radiatively active, temperature dependent clouds, and compare the results to a double gray routine. The double gray method inherently fails to model polychromatic effects in hot Jupiter atmospheres, while picket fence captures these non-gray aspects and performs well compared to fully wavelength-dependent methods. We compare both methods with radiatively active clouds and cloud-free models, assessing the limitations of the double gray method. Although there are broad similarities, the picket fence models have larger day-night side temperature differences, non-isothermal upper atmospheres, and multiwavelength effects in the presence of radiatively active clouds. We model the well-known hot Jupiters HD 189733 b and HD 209458 b. For the hotter HD 209458 b, the picket fence method prevents clouds from thermostating dayside temperatures, resulting in hotter upper atmospheres and the dissipation of dayside clouds. Differences in the temperature structures are then associated with nuanced differences in the circulation patterns and clouds. Models of the cooler HD 189733 b have global cloud coverage, regardless of radiative transfer scheme, whereas there are larger differences in the models of HD 209458 b, particularly in the extent of the partial cloud coverage on its dayside. This results in minor changes to the thermal and reflected light phase curves of HD 189733 b, but more significant differences for the picket fence and double gray versions of HD 209458 b. △ Less

Submitted 26 February, 2024; v1 submitted 2 November, 2023; originally announced November 2023.

Comments: Accepted in ApJ, 31 pages

Understanding an Acquisition Function Family for Bayesian Optimization

Authors: Jiajie Kong, Tony Pourmohamad, Herbert K. H. Lee

Abstract: Bayesian optimization (BO) developed as an approach for the efficient optimization of expensive black-box functions without gradient information. A typical BO paper introduces a new approach and compares it to some alternatives on simulated and possibly real examples to show its efficacy. Yet on a different example, this new algorithm might not be as effective as the alternatives. This paper looks… ▽ More Bayesian optimization (BO) developed as an approach for the efficient optimization of expensive black-box functions without gradient information. A typical BO paper introduces a new approach and compares it to some alternatives on simulated and possibly real examples to show its efficacy. Yet on a different example, this new algorithm might not be as effective as the alternatives. This paper looks at a broader family of approaches to explain the strengths and weaknesses of algorithms in the family, with guidance on what choices might work best on different classes of problems. △ Less

Submitted 16 October, 2023; originally announced October 2023.

A Quantitatively Interpretable Model for Alzheimer's Disease Prediction Using Deep Counterfactuals

Authors: Kwanseok Oh, Da-Woon Heo, Ahmad Wisnu Mulyadi, Wonsik Jung, Eunsong Kang, Kun Ho Lee, Heung-Il Suk

Abstract: Deep learning (DL) for predicting Alzheimer's disease (AD) has provided timely intervention in disease progression yet still demands attentive interpretability to explain how their DL models make definitive decisions. Recently, counterfactual reasoning has gained increasing attention in medical research because of its ability to provide a refined visual explanatory map. However, such visual explan… ▽ More Deep learning (DL) for predicting Alzheimer's disease (AD) has provided timely intervention in disease progression yet still demands attentive interpretability to explain how their DL models make definitive decisions. Recently, counterfactual reasoning has gained increasing attention in medical research because of its ability to provide a refined visual explanatory map. However, such visual explanatory maps based on visual inspection alone are insufficient unless we intuitively demonstrate their medical or neuroscientific validity via quantitative features. In this study, we synthesize the counterfactual-labeled structural MRIs using our proposed framework and transform it into a gray matter density map to measure its volumetric changes over the parcellated region of interest (ROI). We also devised a lightweight linear classifier to boost the effectiveness of constructed ROIs, promoted quantitative interpretation, and achieved comparable predictive performance to DL methods. Throughout this, our framework produces an ``AD-relatedness index'' for each ROI and offers an intuitive understanding of brain status for an individual patient and across patient groups with respect to AD progression. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: 15 pages, 5 figures, 4 tables

Modelling the effect of 3D temperature and chemistry on the cross-correlation signal of transiting ultra-hot Jupiters: A study of 5 chemical species on WASP-76b

Authors: Joost P. Wardenier, Vivien Parmentier, Michael R. Line, Elspeth K. H. Lee

Abstract: Ultra-hot Jupiters are perfect targets for transmission spectroscopy. However, their atmospheres feature strong spatial variations in temperature, chemistry, dynamics, cloud coverage, and scale height. This makes transit observations at high spectral resolution challenging to interpret. In this work, we model the cross-correlation signal of five chemical species (Fe, CO, H$_\text{2}$O, OH, and TiO… ▽ More Ultra-hot Jupiters are perfect targets for transmission spectroscopy. However, their atmospheres feature strong spatial variations in temperature, chemistry, dynamics, cloud coverage, and scale height. This makes transit observations at high spectral resolution challenging to interpret. In this work, we model the cross-correlation signal of five chemical species (Fe, CO, H$_\text{2}$O, OH, and TiO) on WASP-76b, a benchmark ultra-hot Jupiter. We compute phase-dependent high-resolution transmission spectra of 3D SPARC/MITgcm models. The spectra are obtained with gCMCRT, a 3D Monte-Carlo radiative-transfer code. We find that, on top of atmospheric dynamics, the phase-dependent Doppler shift of the absorption lines in the planetary rest frame is shaped by the combined effect of planetary rotation and the unique 3D spatial distribution of chemical species. For species probing the dayside (e.g., refractories or molecules like CO and OH), the two effects act in tandem, leading to increasing blueshifts with orbital phase. For species that are depleted on the dayside (e.g., H$_\text{2}$O and TiO), the two effects act in an opposite manner, and could lead to increasing redshifts during the transit. This behaviour yields species-dependent offsets from a planet's expected $K_\text{p}$ value that can be much larger than planetary wind speeds. The offsets are usually negative for refractory species. We provide an analytical formula to estimate the size of a planet's $K_\text{p}$ offsets, which can serve as a prior for atmospheric retrievals. We conclude that observing the phase-resolved absorption signal of multiple species is key to constraining the 3D thermochemical structure and dynamics of ultra-hot Jupiters. △ Less

Submitted 4 September, 2023; v1 submitted 10 July, 2023; originally announced July 2023.

Comments: 21 pages, 14 figures, accepted for publication in MNRAS

Modelling dynamically driven global cloud formation microphysics in the HAT-P-1b atmosphere

Authors: Elspeth K. H. Lee

Abstract: Insight into the formation and global distribution of cloud particles in exoplanet atmospheres continues to be a key problem to tackle going into the JWST era. Understanding microphysical cloud processes and atmospheric feedback mechanisms in 3D has proven to be a challenging prospect for exoplaneteers. In an effort to address the large computational burden of coupling these models in 3D simulatio… ▽ More Insight into the formation and global distribution of cloud particles in exoplanet atmospheres continues to be a key problem to tackle going into the JWST era. Understanding microphysical cloud processes and atmospheric feedback mechanisms in 3D has proven to be a challenging prospect for exoplaneteers. In an effort to address the large computational burden of coupling these models in 3D simulations, we develop an open source, lightweight and efficient microphysical cloud model for exoplanet atmospheres. `Mini-cloud' is a microphysical based cloud model for exoplanet condensate clouds that can be coupled to contemporary general circulation models (GCMs) and other time dependent simulations. We couple mini-cloud to the Exo-FMS GCM and use a prime JWST target, the hot Jupiter HAT-P-1b, as a test case for the cloud formation module. After 1000+ of days of integration with mini-cloud, our results show a complex 3D cloud structure with cloud properties relating closely the dynamical and temperature properties of the atmosphere. Current transit and emission spectra data are best fit with a reduced cloud particle number density compared to the nominal simulation, with our simulated JWST NIRISS SOSS spectra showing promising prospects for characterising the atmosphere in detail. Overall, our study is another small step in first principles 3D exoplanet cloud formation microphysical modelling. We suggest that additional physics not included in the present model, such as coagulation, are required to reduce the number density of particles to appropriately observed levels. △ Less

Submitted 5 July, 2023; originally announced July 2023.

Comments: MNRAS Accepted

Examining NHD vs QHD in the GCM THOR with non-grey radiative transfer for the hot Jupiter regime

Authors: Pascal A. Noti, Elspeth K. H. Lee, Russell Deitrick, Mark Hammond

Abstract: Global circulation models (GCMs) play an important role in contemporary investigations of exoplanet atmospheres. Different GCMs evolve various sets of dynamical equations which can result in obtaining different atmospheric properties between models. In this study, we investigate the effect of different dynamical equation sets on the atmospheres of hot Jupiter exoplanets. We compare GCM simulations… ▽ More Global circulation models (GCMs) play an important role in contemporary investigations of exoplanet atmospheres. Different GCMs evolve various sets of dynamical equations which can result in obtaining different atmospheric properties between models. In this study, we investigate the effect of different dynamical equation sets on the atmospheres of hot Jupiter exoplanets. We compare GCM simulations using the quasi-primitive dynamical equations (QHD) and the deep Navier-Stokes equations (NHD) in the GCM THOR. We utilise a two-stream non-grey "picket-fence" scheme to increase the realism of the radiative transfer calculations. We perform GCM simulations covering a wide parameter range grid of system parameters in the population of exoplanets. Our results show significant differences between simulations with the NHD and QHD equation sets at lower gravity, higher rotation rates or at higher irradiation temperatures. The chosen parameter range shows the relevance of choosing dynamical equation sets dependent on system and planetary properties. Our results show the climate states of hot Jupiters seem to be very diverse, where exceptions to prograde superrotation can often occur. Overall, our study shows the evolution of different climate states which arise just due to different selections of Navier-Stokes equations and approximations. We show the divergent behaviour of approximations used in GCMs for Earth, but applied for non Earth-like planets. △ Less

Submitted 3 July, 2023; originally announced July 2023.

Comments: accepted MNRAS (21 June 2023)

Dynamically coupled kinetic chemistry in brown dwarf atmospheres I. Performing global scale kinetic modelling

Authors: Elspeth K. H. Lee, Xianyu Tan, Shang-Min Tsai

Abstract: The atmospheres of brown dwarfs have been long observed to exhibit a multitude of non-equilibrium chemical signatures and spectral variability across the L, T and Y spectral types. We aim to investigate the link between the large-scale 3D atmospheric dynamics and time-dependent chemistry in the brown dwarf regime, and to assess its impact on spectral variability. We couple the miniature kinetic ch… ▽ More The atmospheres of brown dwarfs have been long observed to exhibit a multitude of non-equilibrium chemical signatures and spectral variability across the L, T and Y spectral types. We aim to investigate the link between the large-scale 3D atmospheric dynamics and time-dependent chemistry in the brown dwarf regime, and to assess its impact on spectral variability. We couple the miniature kinetic chemistry module `mini-chem' to the Exo-FMS general circulation model (GCM). We then perform a series of idealised brown dwarf regime atmospheric models to investigate the dynamical 3D chemical structures produced by our simulations. The GCM output is post-processed using a 3D radiative-transfer model to investigate hemisphere-dependent spectral signatures and rotational variability. Our results show the expected strong non-equilibrium chemical behaviour brought on by vertical mixing as well as global spacial variations due to zonal flows. Chemical species are generally globally homogenised, showing variations of $\pm$10\% or less, dependent on pressure level, and follow the dynamical structures present in the atmosphere. However, we find localised storm regions and eddies can show higher contrasts, up to $\pm$100\%, in mixing ratio compared to the background global mean. This initial study represents another step in understanding the connection between three-dimensional atmospheric flows in brown dwarfs and their rich chemical inventories. △ Less

Submitted 6 June, 2023; originally announced June 2023.

Comments: MNRAS Accepted: 5 June 2023

Day-night transport induced chemistry and clouds on WASP-39b I: Gas-phase composition

Authors: Shang-Min Tsai, Julianne I. Moses, Diana Powell, Elspeth K. H. Lee

Abstract: JWST has recently detected the first robust photochemical product on an exoplanet: sulfur dioxide (SO$_2$) on WASP-39b (Rustamkulov et al. 2023; Alderson et al. 2023; Tsai et al. 2023b). The data from the NIRISS instrument also reveal signs of partial coverage of clouds (Feinstein et al. 2023). Most of the previous studies have focused on interpreting spectral data with 1D models. To explore how t… ▽ More JWST has recently detected the first robust photochemical product on an exoplanet: sulfur dioxide (SO$_2$) on WASP-39b (Rustamkulov et al. 2023; Alderson et al. 2023; Tsai et al. 2023b). The data from the NIRISS instrument also reveal signs of partial coverage of clouds (Feinstein et al. 2023). Most of the previous studies have focused on interpreting spectral data with 1D models. To explore how the chemical species and cloud particles are altered by global circulation, we applied a 2D photochemical model and a 2D microphysical cloud model separately to post-process the thermal and dynamical structures simulated by a 3D general circulation model (GCM) of WASP-39b. We found that SO$_2$ produced by photochemistry on the dayside can be transported to the nightside owing to the efficient replenishment of horizontal transport. The morning-evening limb differences in methane (CH$_4$) abundances predicted by the 1D models disappeared after horizontal transport is included. Similarly, the inclusion of horizontal transport also reduced the limb differences in SO$_2$. Our modeling results suggest that the fast zonal wind results in minimal or negligible limb asymmetry in composition. Based on the synthetic spectra generated by our 2D atmosphere simulations, we propose that observing SO$_2$ absorption in the emission spectra of WASP-39b at different phases may offer opportunities to probe the horizontal quenching process of photochemical products. We will focus on the gas-phase chemistry in this paper and leave the results regarding clouds in the subsequent paper as part of the series. △ Less

Submitted 20 March, 2024; v1 submitted 30 May, 2023; originally announced May 2023.

Comments: 10 pages, 6 figures, published in ApJL

Journal ref: ApJL 959 L30 (2023)

Deposition and alignment of fiber suspensions by dip coating

Authors: Deok-Hoon Jeong, Langqi Xing, Michael Ka Ho Lee, Nathan Vani, Alban Sauret

Abstract: The dip coating of suspensions made of monodisperse non-Brownian spherical particles dispersed in a Newtonian fluid leads to different coating regimes depending on the ratio of the particle diameter to the thickness of the film entrained on the substrate. In particular, dilute particles dispersed in the liquid are entrained only above a threshold value of film thickness. In the case of anisotropic… ▽ More The dip coating of suspensions made of monodisperse non-Brownian spherical particles dispersed in a Newtonian fluid leads to different coating regimes depending on the ratio of the particle diameter to the thickness of the film entrained on the substrate. In particular, dilute particles dispersed in the liquid are entrained only above a threshold value of film thickness. In the case of anisotropic particles, in particular fibers, the smallest characteristic dimension will control the entrainment of the particle. Furthermore, it is possible to control the orientation of the anisotropic particles depending on the substrate geometry. To test the hypotheses, we performed dip-coating experiments with dilute suspensions of non-Brownian fibers with different length-to-diameter aspect ratios. We characterize the number of fibers entrained on the surface of the substrate as a function of the withdrawal velocity, allowing us to estimate a threshold capillary number below which all the particles remain in the liquid bath. Besides, we measure the angular distribution of the entrained fibers for two different substrate geometries: flat plates and cylindrical rods. We then measure the film thickness for more concentrated fiber suspensions. The entrainment of the fibers on a flat plate and a cylindrical rod is primarily controlled by the smaller characteristic length of the fibers: their diameter. At first order, the entrainment threshold scales similarly to that of spherical particles. The length of the fibers only appears to have a minor influence on the entrainment threshold. No preferential alignment is observed for non-Brownian fibers on a flat plate, except for very thin films, whereas the fibers tend to align themselves along the axis of a cylindrical rod for a large enough ratio of the fiber length to the radius of the cylindrical rod. △ Less

Submitted 1 May, 2023; originally announced May 2023.

Tetrahedral triple-Q magnetic ordering and large spontaneous Hall conductivity in the metallic triangular antiferromagnet Co1/3TaS2

Authors: Pyeongjae Park, Woonghee Cho, Chaebin Kim, Yeochan An, Yoon-Gu Kang, Maxim Avdeev, Romain Sibille, Kazuki Iida, Ryoichi Kajimoto, Ki Hoon Lee, Woori Ju, En-Jin Cho, Han-Jin Noh, Myung Joon Han, Shang-Shun Zhang, Cristian D. Batista, Je-Geun Park

Abstract: The triangular lattice antiferromagnet (TLAF) has been the standard paradigm of frustrated magnetism for several decades. The most common magnetic ordering in insulating TLAFs is the 120 structure. However, a new triple-Q chiral ordering can emerge in metallic TLAFs, representing the short wavelength limit of magnetic skyrmion crystals. We report the metallic TLAF Co1/3TaS2 as the first example of… ▽ More The triangular lattice antiferromagnet (TLAF) has been the standard paradigm of frustrated magnetism for several decades. The most common magnetic ordering in insulating TLAFs is the 120 structure. However, a new triple-Q chiral ordering can emerge in metallic TLAFs, representing the short wavelength limit of magnetic skyrmion crystals. We report the metallic TLAF Co1/3TaS2 as the first example of tetrahedral triple-Q magnetic ordering with the associated topological Hall effect (non-zero σ_{xy}(H=0)). We also present a theoretical framework that describes the emergence of this magnetic ground state, which is further supported by the electronic structure measured by angle-resolved photoemission spectroscopy. Additionally, our measurements of the inelastic neutron scattering cross section are consistent with the calculated dynamical structure factor of the tetrahedral triple-Q state. △ Less

Submitted 12 November, 2023; v1 submitted 7 March, 2023; originally announced March 2023.

Comments: 23 pages, 4 figures, SI not included. Accepted for publication in Nature Communications

Journal ref: Nature Communications 14, 8346 (2023)

A Mini-Chemical Scheme with Net Reactions for 3D GCMs II. 3D thermochemical modelling of WASP-39b and HD 189733b

Authors: Elspeth K. H. Lee, Shang-Min Tsai, Mark Hammond, Xianyu Tan

Abstract: The chemical inventory of hot Jupiter (HJ) exoplanets atmospheres continue to be observed by various ground and space based instruments in increasing detail and precision. It is expected that some HJs will exhibit strong non-equilibrium chemistry characteristics in their atmospheres, which might be inferred from spectral observations. We aim to model the three dimensional thermochemical non-equili… ▽ More The chemical inventory of hot Jupiter (HJ) exoplanets atmospheres continue to be observed by various ground and space based instruments in increasing detail and precision. It is expected that some HJs will exhibit strong non-equilibrium chemistry characteristics in their atmospheres, which might be inferred from spectral observations. We aim to model the three dimensional thermochemical non-equilibrium chemistry in the atmospheres of the HJs WASP-39b and HD 189733b. We couple a lightweight, reduced chemical network `mini-chem' that utilises net reaction rate tables to the Exo-FMS General Circulation Model (GCM). We perform GCM models of the exoplanets WASP-39b and HD 189733b as case studies of the coupled mini-chem scheme. The GCM results are then post-processed using the 3D radiative-transfer model gCMCRT to produce transmission and emission spectra to assess the impact of non-equilibrium chemistry on their observable properties. Both simulations show significant departures from chemical equilibrium (CE) due to the dynamical motions of the atmosphere. The spacial distribution of species generally follows closely the dynamical features of the atmosphere rather than the temperature field. Each molecular species exhibits a different quench level in the simulations, also dependent on the latitude of the planet. Major differences are seen in the transmission and emission spectral features between the CE and kinetic models. Our simulations indicate that considering the 3D kinetic chemical structures of HJ atmospheres has an important impact on physical interpretation of observational data. Drawing bulk atmospheric parameters from fitting feature strengths may lead to inaccurate interpretation of chemical conditions in the atmosphere of HJs. Our open source mini-chem module is simple to couple with contemporary HJ GCM models without substantially increasing required computational resources. △ Less

Submitted 19 February, 2023; originally announced February 2023.

Comments: A&A Accepted 17 Feb 2023

Journal ref: A&A 672, A110 (2023)

A broadband thermal emission spectrum of the ultra-hot Jupiter WASP-18b

Authors: Louis-Philippe Coulombe, Björn Benneke, Ryan Challener, Anjali A. A. Piette, Lindsey S. Wiser, Megan Mansfield, Ryan J. MacDonald, Hayley Beltz, Adina D. Feinstein, Michael Radica, Arjun B. Savel, Leonardo A. Dos Santos, Jacob L. Bean, Vivien Parmentier, Ian Wong, Emily Rauscher, Thaddeus D. Komacek, Eliza M. -R. Kempton, Xianyu Tan, Mark Hammond, Neil T. Lewis, Michael R. Line, Elspeth K. H. Lee, Hinna Shivkumar, Ian J. M. Crossfield , et al. (51 additional authors not shown)

Abstract: Close-in giant exoplanets with temperatures greater than 2,000 K (''ultra-hot Jupiters'') have been the subject of extensive efforts to determine their atmospheric properties using thermal emission measurements from the Hubble and Spitzer Space Telescopes. However, previous studies have yielded inconsistent results because the small sizes of the spectral features and the limited information conten… ▽ More Close-in giant exoplanets with temperatures greater than 2,000 K (''ultra-hot Jupiters'') have been the subject of extensive efforts to determine their atmospheric properties using thermal emission measurements from the Hubble and Spitzer Space Telescopes. However, previous studies have yielded inconsistent results because the small sizes of the spectral features and the limited information content of the data resulted in high sensitivity to the varying assumptions made in the treatment of instrument systematics and the atmospheric retrieval analysis. Here we present a dayside thermal emission spectrum of the ultra-hot Jupiter WASP-18b obtained with the NIRISS instrument on JWST. The data span 0.85 to 2.85 $μ$m in wavelength at an average resolving power of 400 and exhibit minimal systematics. The spectrum shows three water emission features (at $>$6$σ$ confidence) and evidence for optical opacity, possibly due to H$^-$, TiO, and VO (combined significance of 3.8$σ$). Models that fit the data require a thermal inversion, molecular dissociation as predicted by chemical equilibrium, a solar heavy element abundance (''metallicity'', M/H = 1.03$_{-0.51}^{+1.11}$ $\times$ solar), and a carbon-to-oxygen (C/O) ratio less than unity. The data also yield a dayside brightness temperature map, which shows a peak in temperature near the sub-stellar point that decreases steeply and symmetrically with longitude toward the terminators. △ Less

Submitted 20 January, 2023; v1 submitted 19 January, 2023; originally announced January 2023.

Comments: JWST ERS bright star observations. Uploaded to inform JWST Cycle 2 proposals. Manuscript under review. 50 pages, 14 figures, 2 tables

Characterizing quantile-varying covariate effects under the accelerated failure time model

Authors: Harrison T. Reeder, Kyu Ha Lee, Sebastien Haneuse

Abstract: An important task in survival analysis is choosing a structure for the relationship between covariates of interest and the time-to-event outcome. For example, the accelerated failure time (AFT) model structures each covariate effect as a constant multiplicative shift in the outcome distribution across all survival quantiles. Though parsimonious, this structure cannot detect or capture effects that… ▽ More An important task in survival analysis is choosing a structure for the relationship between covariates of interest and the time-to-event outcome. For example, the accelerated failure time (AFT) model structures each covariate effect as a constant multiplicative shift in the outcome distribution across all survival quantiles. Though parsimonious, this structure cannot detect or capture effects that differ across quantiles of the distribution, a limitation that is analogous to only permitting proportional hazards in the Cox model. To address this, we propose a general framework for quantile-varying multiplicative effects under the AFT model. Specifically, we embed flexible regression structures within the AFT model, and derive a novel formula for interpretable effects on the quantile scale. A regression standardization scheme based on the g-formula is proposed to enable estimation of both covariate-conditional and marginal effects for an exposure of interest. We implement a user-friendly Bayesian approach for estimation and quantification of uncertainty, while accounting for left truncation and complex censoring. We emphasize the intuitive interpretation of this model through numerical and graphical tools, and illustrate its performance by application to a study of Alzheimer's disease and dementia. △ Less

Submitted 8 January, 2023; originally announced January 2023.

Comments: This is the pre-peer reviewed, "submitted" version of the manuscript published in final form in Biostatistics by Oxford University Press at the below citation/doi. This upload will be updated with the final peer-reviewed "accepted" version of the manuscript following a 24 month embargo period

Journal ref: Biostatistics (Oxford, England), kxac052 (2023)

Twisted carotenoids do not support efficient intramolecular singlet fission in the orange carotenoid protein

Authors: George A. Sutherland, James P. Pidgeon, Harrison Ka Hin Lee, Matthew S. Proctor, Andrew Hitchcock, Shuangqing Wang, Dimitri Chekulaev, Wing Chung Tsoi, Matthew P. Johnson, C. Neil Hunter, Jenny Clark

Abstract: Singlet exciton fission is the spin-allowed generation of two triplet electronic excited states from a singlet state. Intramolecular singlet fission has been suggested to occur on individual carotenoid molecules within protein complexes, provided the conjugated backbone is twisted out-of-plane. However, this hypothesis has only been forwarded in protein complexes containing multiple carotenoids an… ▽ More Singlet exciton fission is the spin-allowed generation of two triplet electronic excited states from a singlet state. Intramolecular singlet fission has been suggested to occur on individual carotenoid molecules within protein complexes, provided the conjugated backbone is twisted out-of-plane. However, this hypothesis has only been forwarded in protein complexes containing multiple carotenoids and bacteriochlorophylls in close contact. To test the hypothesis on twisted carotenoids in a 'minimal' one-carotenoid system, we study the orange carotenoid protein (OCP). OCP exists in two forms: in its orange form (OCPo), the single bound carotenoid is twisted, whereas in its red form (OCPr), the carotenoid is planar. To enable room-temperature spectroscopy on canthaxanthin-binding OCPo and OCPr without laser-induced photoconversion, we trap them in trehalose glass. Using transient absorption spectroscopy, we show that there is no evidence of long-lived triplet generation through intramolecular singlet fission, despite the canthaxanthin twist in OCPo. △ Less

Submitted 24 November, 2022; originally announced November 2022.

Photochemically-produced SO$_2$ in the atmosphere of WASP-39b

Authors: Shang-Min Tsai, Elspeth K. H. Lee, Diana Powell, Peter Gao, Xi Zhang, Julianne Moses, Eric Hébrard, Olivia Venot, Vivien Parmentier, Sean Jordan, Renyu Hu, Munazza K. Alam, Lili Alderson, Natalie M. Batalha, Jacob L. Bean, Björn Benneke, Carver J. Bierson, Ryan P. Brady, Ludmila Carone, Aarynn L. Carter, Katy L. Chubb, Julie Inglis, Jérémy Leconte, Mercedes Lopez-Morales, Yamila Miguel , et al. (60 additional authors not shown)

Abstract: Photochemistry is a fundamental process of planetary atmospheres that regulates the atmospheric composition and stability. However, no unambiguous photochemical products have been detected in exoplanet atmospheres to date. Recent observations from the JWST Transiting Exoplanet Early Release Science Program found a spectral absorption feature at 4.05 $μ$m arising from SO$_2$ in the atmosphere of WA… ▽ More Photochemistry is a fundamental process of planetary atmospheres that regulates the atmospheric composition and stability. However, no unambiguous photochemical products have been detected in exoplanet atmospheres to date. Recent observations from the JWST Transiting Exoplanet Early Release Science Program found a spectral absorption feature at 4.05 $μ$m arising from SO$_2$ in the atmosphere of WASP-39b. WASP-39b is a 1.27-Jupiter-radii, Saturn-mass (0.28 M$_J$) gas giant exoplanet orbiting a Sun-like star with an equilibrium temperature of $\sim$1100 K. The most plausible way of generating SO$_2$ in such an atmosphere is through photochemical processes. Here we show that the SO$_2$ distribution computed by a suite of photochemical models robustly explains the 4.05 $μ$m spectral feature identified by JWST transmission observations with NIRSpec PRISM (2.7$σ$) and G395H (4.5$σ$). SO$_2$ is produced by successive oxidation of sulphur radicals freed when hydrogen sulphide (H$_2$S) is destroyed. The sensitivity of the SO$_2$ feature to the enrichment of the atmosphere by heavy elements (metallicity) suggests that it can be used as a tracer of atmospheric properties, with WASP-39b exhibiting an inferred metallicity of $\sim$10$\times$ solar. We further point out that SO$_2$ also shows observable features at ultraviolet and thermal infrared wavelengths not available from the existing observations. △ Less

Submitted 24 March, 2023; v1 submitted 18 November, 2022; originally announced November 2022.

Comments: 39 pages, 14 figures, accepted to be published in Nature

Early Release Science of the Exoplanet WASP-39b with JWST NIRSpec G395H

Authors: Lili Alderson, Hannah R. Wakeford, Munazza K. Alam, Natasha E. Batalha, Joshua D. Lothringer, Jea Adams Redai, Saugata Barat, Jonathan Brande, Mario Damiano, Tansu Daylan, Néstor Espinoza, Laura Flagg, Jayesh M. Goyal, David Grant, Renyu Hu, Julie Inglis, Elspeth K. H. Lee, Thomas Mikal-Evans, Lakeisha Ramos-Rosado, Pierre-Alexis Roy, Nicole L. Wallack, Natalie M. Batalha, Jacob L. Bean, Björn Benneke, Zachory K. Berta-Thompson , et al. (67 additional authors not shown)

Abstract: Measuring the abundances of carbon and oxygen in exoplanet atmospheres is considered a crucial avenue for unlocking the formation and evolution of exoplanetary systems. Access to an exoplanet's chemical inventory requires high-precision observations, often inferred from individual molecular detections with low-resolution space-based and high-resolution ground-based facilities. Here we report the m… ▽ More Measuring the abundances of carbon and oxygen in exoplanet atmospheres is considered a crucial avenue for unlocking the formation and evolution of exoplanetary systems. Access to an exoplanet's chemical inventory requires high-precision observations, often inferred from individual molecular detections with low-resolution space-based and high-resolution ground-based facilities. Here we report the medium-resolution (R$\sim$600) transmission spectrum of an exoplanet atmosphere between 3-5 $μ$m covering multiple absorption features for the Saturn-mass exoplanet WASP-39b, obtained with JWST NIRSpec G395H. Our observations achieve 1.46x photon precision, providing an average transit depth uncertainty of 221 ppm per spectroscopic bin, and present minimal impacts from systematic effects. We detect significant absorption from CO$_2$ (28.5$σ$) and H$_2$O (21.5$σ$), and identify SO$_2$ as the source of absorption at 4.1 $μ$m (4.8$σ$). Best-fit atmospheric models range between 3 and 10x solar metallicity, with sub-solar to solar C/O ratios. These results, including the detection of SO$_2$, underscore the importance of characterising the chemistry in exoplanet atmospheres, and showcase NIRSpec G395H as an excellent mode for time series observations over this critical wavelength range. △ Less

Submitted 18 November, 2022; originally announced November 2022.

Comments: 44 pages, 11 figures, 3 tables. Resubmitted after revision to Nature

CAMEMBERT: A Mini-Neptunes GCM Intercomparison, Protocol Version 1.0. A CUISINES Model Intercomparison Project

Authors: Duncan A. Christie, Elspeth K. H. Lee, Hamish Innes, Pascal A. Noti, Benjamin Charnay, Thomas J. Fauchez, Nathan J. Mayne, Russell Deitrick, Feng Ding, Jennifer J. Greco, Mark Hammond, Isaac Malsky, Avi Mandell, Emily Rauscher, Michael T. Roman, Denis E. Sergeev, Linda Sohl, Maria E. Steinrueck, Martin Turbet, Eric T. Wolf, Maria Zamyatina, Ludmila Carone

Abstract: With an increased focus on the observing and modelling of mini-Neptunes, there comes a need to better understand the tools we use to model their atmospheres. In this paper, we present the protocol for the CAMEMBERT (Comparing Atmospheric Models of Extrasolar Mini-neptunes Building and Envisioning Retrievals and Transits) project, an intercomparison of general circulation models (GCMs) used by the… ▽ More With an increased focus on the observing and modelling of mini-Neptunes, there comes a need to better understand the tools we use to model their atmospheres. In this paper, we present the protocol for the CAMEMBERT (Comparing Atmospheric Models of Extrasolar Mini-neptunes Building and Envisioning Retrievals and Transits) project, an intercomparison of general circulation models (GCMs) used by the exoplanetary science community to simulate the atmospheres of mini-Neptunes. We focus on two targets well studied both observationally and theoretically with planned JWST Cycle 1 observations: the warm GJ~1214b and the cooler K2-18b. For each target, we consider a temperature-forced case, a clear sky dual-grey radiative transfer case, and a clear sky multi band radiative transfer case, covering a range of complexities and configurations where we know differences exist between GCMs in the literature. This paper presents all the details necessary to participate in the intercomparison, with the intention of presenting the results in future papers. Currently, there are eight GCMs participating (ExoCAM, Exo-FMS, FMS PCM, Generic PCM, MITgcm, RM-GCM, THOR, and the UM), and membership in the project remains open. Those interested in participating are invited to contact the authors. △ Less

Submitted 8 November, 2022; originally announced November 2022.

Comments: Accepted to PSJ

Fractals Generated by Modifying Aperiodic Substitution Tilings

Authors: Kah Heng Lee

Abstract: This study proposes a method for producing an infinite number of fractals using aperiodic substitution tilings, exemplified by the Ammann Chair tiling. Higher-order substitutions of aperiodic tilings are utilized in relation to the Sierpinski carpet concept. The similarity dimensions of the fractals generated by the Ammann Chair tiling are calculated and shown to be dense. A fractal image generato… ▽ More This study proposes a method for producing an infinite number of fractals using aperiodic substitution tilings, exemplified by the Ammann Chair tiling. Higher-order substitutions of aperiodic tilings are utilized in relation to the Sierpinski carpet concept. The similarity dimensions of the fractals generated by the Ammann Chair tiling are calculated and shown to be dense. A fractal image generator was implemented in the Java programming language and is freely available for public use at https://github.com/KahHengLee/Ammann-Chair-Fractal.git △ Less

Submitted 18 May, 2023; v1 submitted 18 October, 2022; originally announced October 2022.

Comments: 10 pages, 12 figures

The Mantis Network IV: A titanium cold-trap on the ultra-hot Jupiter WASP-121 b

Authors: H. J. Hoeijmakers, D. Kitzmann, B. M. Morris, B. Prinoth, N. Borsato, B. Thorsbro, L. Pino, E. K. H. Lee, C. Akın, J. V. Seidel, J. L. Birkby, R. Allart, K. Heng

Abstract: Observations of WASP-121 b have suggested an under-abundance of titanium and titanium-oxide from its terminator region. In this study, we aim to determine whether this depletion is global by investigating the day-side emission spectrum. We analyse 8 epochs of high-resolution spectra obtained with ESPRESSO, targeting orbital phases when the day-side is in view. We use a cross-correlation method to… ▽ More Observations of WASP-121 b have suggested an under-abundance of titanium and titanium-oxide from its terminator region. In this study, we aim to determine whether this depletion is global by investigating the day-side emission spectrum. We analyse 8 epochs of high-resolution spectra obtained with ESPRESSO, targeting orbital phases when the day-side is in view. We use a cross-correlation method to search for various atoms, TiO and VO and compare to models. We constrain the velocities and phase-function of the emission signal using a Bayesian framework. We report significant detections of Ca I, V I, Cr I, Mn I, Fe I, Co I and Ni I, but not T i or TiO. Models containing Ti are unable to reproduce the data. The detected signals are consistent with the known orbital and systemic velocities and with peak emission originating from the sub-stellar point. We find that Ti is depleted from regions of the atmosphere where transmission and emission spectroscopy are sensitive. We interpret this as evidence for the night-side condensation of titanium, preventing it from being mixed back into the upper layers of the atmosphere elsewhere on the planet. Species with lower condensation temperatures are unaffected, implying sharp chemical transitions exist between ultra-hot Jupiters that have slight differences in temperature or dynamical properties. As TiO can act as a strong source of stratospheric heating, cold-trapping creates a coupling between the thermal structures on the day-side and night-side, and thus condensation chemistry needs to be included in global circulation models. Observed elemental abundances in hot Jupiters are not reliably representative of bulk abundances unless night-side condensation is accounted for or the planet is hot enough to avoid night-side cold-traps entirely. Planetary rotation may significantly lower the apparent orbital velocity of emission signals. △ Less

Submitted 16 January, 2024; v1 submitted 23 October, 2022; originally announced October 2022.

Comments: 23 pages, submitted to A\&A on September 13, 2022, accepted on October 6, 2023

The Mantis Network II: Examining the 3D high-resolution observable properties of the UHJs WASP-121b and WASP-189b through GCM modelling

Authors: Elspeth K. H. Lee, Bibiana Prinoth, Daniel Kitzmann, Shang-Min Tsai, Jens Hoeijmakers, Nicholas W. Borsato, Kevin Heng

Abstract: The atmospheres of ultra hot Jupiters (UHJs) are prime targets for the detection of molecules and atoms at both low and high spectral resolution. We study the atmospheres of the UHJs WASP-121b and WASP-189b by performing 3D general circulation models (GCMs) of these planets using high temperature correlated-k opacity schemes with ultra-violet (UV) absorbing species included. The GCM results are th… ▽ More The atmospheres of ultra hot Jupiters (UHJs) are prime targets for the detection of molecules and atoms at both low and high spectral resolution. We study the atmospheres of the UHJs WASP-121b and WASP-189b by performing 3D general circulation models (GCMs) of these planets using high temperature correlated-k opacity schemes with ultra-violet (UV) absorbing species included. The GCM results are then post-processed at low and high spectral resolutions and compared to available data. The high resolution results are cross-correlated with molecular and atomic templates to produce mock molecular detections. Our GCM models produce similar temperature-pressure (T-p) structure trends to previous 1D radiative-convective equilibrium models of UHJs. Furthermore, the inclusion of UV opacities greatly shapes the thermal and dynamical properties of the high-altitude, low-pressure regions of the UHJ atmospheres, with sharp T-p inversions due to the absorption of UV light. This suggests that optical wavelength, high-resolution observations probe a dynamically distinct upper atmospheric region, rather than the deeper jet forming layers. △ Less

Submitted 21 October, 2022; originally announced October 2022.

Comments: MNRAS accepted: 8 August 2022, published 11 August 2022

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 517, Issue 1, November 2022, Pages 240 256

Designing topological acoustic lattices via electroacoustic analogies

Authors: Hasan Al Ba'ba'a, Kyung Hoon Lee, Qiming Wang

Abstract: Topological acoustics has recently witnessed a spurt in research activity, owing to their unprecedented properties transcending typical wave phenomena. In recent years, the use of coupled arrays of acoustic chambers has gained popularity in designing topological acoustic systems. In their common form, an array of acoustic chambers with relatively large volume is coupled via narrow channels. Such c… ▽ More Topological acoustics has recently witnessed a spurt in research activity, owing to their unprecedented properties transcending typical wave phenomena. In recent years, the use of coupled arrays of acoustic chambers has gained popularity in designing topological acoustic systems. In their common form, an array of acoustic chambers with relatively large volume is coupled via narrow channels. Such configuration is generally modeled as a full three-dimensional system, requiring extended computational time for simulating its harmonic response. To this end, this paper establishes a comprehensive mathematical treatment of the use of electroacoustic analogies for designing topological acoustic lattices. The potential of such analytical approach is demonstrated via two types of topological systems: (i) edge states with quantized winding numbers in an acoustic diatomic lattice and (ii) valley Hall transition in an acoustic honeycomb lattice that leads to robust waveguiding. In both cases, the established analytical approach exhibits an excellent agreement with the full three-dimensional model, whether in dispersion analyses or the response of an acoustic system with a finite number of cells. The established analytical framework is invaluable for designing a variety of acoustic topological insulators with minimal computational cost. △ Less

Submitted 3 May, 2023; v1 submitted 14 September, 2022; originally announced September 2022.

Modelling of physical systems with a Hopf bifurcation using mechanistic models and machine learning

Authors: K. H. Lee, D. A. W. Barton, L. Renson

Abstract: We propose a new hybrid modelling approach that combines a mechanistic model with a machine-learnt model to predict the limit cycle oscillations of physical systems with a Hopf bifurcation. The mechanistic model is an ordinary differential equation normal-form model capturing the bifurcation structure of the system. A data-driven mapping from this model to the experimental observations is then ide… ▽ More We propose a new hybrid modelling approach that combines a mechanistic model with a machine-learnt model to predict the limit cycle oscillations of physical systems with a Hopf bifurcation. The mechanistic model is an ordinary differential equation normal-form model capturing the bifurcation structure of the system. A data-driven mapping from this model to the experimental observations is then identified based on experimental data using machine learning techniques. The proposed method is first demonstrated numerically on a Van der Pol oscillator and a three-degree-of-freedom aeroelastic model. It is then applied to model the behaviour of a physical aeroelastic structure exhibiting limit cycle oscillations during wind tunnel tests. The method is shown to be general, data-efficient and to offer good accuracy without any prior knowledge about the system other than its bifurcation structure. △ Less

Submitted 7 September, 2022; originally announced September 2022.

Identification of carbon dioxide in an exoplanet atmosphere

Authors: The JWST Transiting Exoplanet Community Early Release Science Team, Eva-Maria Ahrer, Lili Alderson, Natalie M. Batalha, Natasha E. Batalha, Jacob L. Bean, Thomas G. Beatty, Taylor J. Bell, Björn Benneke, Zachory K. Berta-Thompson, Aarynn L. Carter, Ian J. M. Crossfield, Néstor Espinoza, Adina D. Feinstein, Jonathan J. Fortney, Neale P. Gibson, Jayesh M. Goyal, Eliza M. -R. Kempton, James Kirk, Laura Kreidberg, Mercedes López-Morales, Michael R. Line, Joshua D. Lothringer, Sarah E. Moran, Sagnick Mukherjee , et al. (107 additional authors not shown)

Abstract: Carbon dioxide (CO2) is a key chemical species that is found in a wide range of planetary atmospheres. In the context of exoplanets, CO2 is an indicator of the metal enrichment (i.e., elements heavier than helium, also called "metallicity"), and thus formation processes of the primary atmospheres of hot gas giants. It is also one of the most promising species to detect in the secondary atmospheres… ▽ More Carbon dioxide (CO2) is a key chemical species that is found in a wide range of planetary atmospheres. In the context of exoplanets, CO2 is an indicator of the metal enrichment (i.e., elements heavier than helium, also called "metallicity"), and thus formation processes of the primary atmospheres of hot gas giants. It is also one of the most promising species to detect in the secondary atmospheres of terrestrial exoplanets. Previous photometric measurements of transiting planets with the Spitzer Space Telescope have given hints of the presence of CO2 but have not yielded definitive detections due to the lack of unambiguous spectroscopic identification. Here we present the detection of CO2 in the atmosphere of the gas giant exoplanet WASP-39b from transmission spectroscopy observations obtained with JWST as part of the Early Release Science Program (ERS). The data used in this study span 3.0 to 5.5 μm in wavelength and show a prominent CO2 absorption feature at 4.3 μm (26σ significance). The overall spectrum is well matched by one-dimensional, 10x solar metallicity models that assume radiative-convective-thermochemical equilibrium and have moderate cloud opacity. These models predict that the atmosphere should have water, carbon monoxide, and hydrogen sulfide in addition to CO2, but little methane. Furthermore, we also tentatively detect a small absorption feature near 4.0 μm that is not reproduced by these models. △ Less

Submitted 24 August, 2022; originally announced August 2022.

Comments: 27 pages, 6 figures, Accepted for publication in Nature, data and models available at https://doi.10.5281/zenodo.6959427

Significant thermal Hall effect in the 3$d$ cobalt Kitaev system $Na_2Co_2TeO_6$

Authors: Heejun Yang, Chaebin Kim, Ysun Choi, Jun Han Lee, Gaoting Lin, Jie Ma, Marie Kratochvílová, Petr Proschek, Eun-Gook Moon, Ki Hoon Lee, Yoon Seok Oh, Je-Geun Park

Abstract: Kitaev physics has recently attracted attention in condensed matter for its anticipated novel quantum spin liquid state. The thermal transport measurement is crucial for probing the novel features of charge-neutral quasiparticles. In this letter, we report a significant thermal Hall effect in $Na_2Co_2TeO_6$ (NCTO), a Kitaev quantum spin liquid candidate, when the magnetic field is applied along t… ▽ More Kitaev physics has recently attracted attention in condensed matter for its anticipated novel quantum spin liquid state. The thermal transport measurement is crucial for probing the novel features of charge-neutral quasiparticles. In this letter, we report a significant thermal Hall effect in $Na_2Co_2TeO_6$ (NCTO), a Kitaev quantum spin liquid candidate, when the magnetic field is applied along the out-of-plane direction of the honeycomb plane. The thermal conductivity ($κ_{xx}$) and thermal Hall conductivity ($κ_{xy}$) in NCTO reveal distinct magnetic field dependences below and above the Neel temperature ($T_N$) of 27 K. For $T>T_N$, $κ_{xx}$ has a monotonic decrease in the field dependence, while $κ_{xy}$ persists up to $T^*$ = 150 K. On the other hand, both $κ_{xx}$ and $κ_{xy}$ exhibit complex field dependence for $T<T_N$. △ Less

Submitted 25 August, 2022; v1 submitted 3 August, 2022; originally announced August 2022.

Comments: 13 pages, 3 figures, accepted for publication in Physical Review B

Journal ref: Phys. Rev. B 106, L081116 (2022)

Patchy nightside clouds on ultra-hot Jupiters: General Circulation Model simulations with radiatively active cloud tracers

Authors: Thaddeus D. Komacek, Xianyu Tan, Peter Gao, Elspeth K. H. Lee

Abstract: The atmospheres of ultra-hot Jupiters have been characterized in detail through recent phase curve and low- and high-resolution emission and transmission spectroscopic observations. Previous numerical studies have analyzed the effect of the localized recombination of hydrogen on the atmospheric dynamics and heat transport of ultra-hot Jupiters, finding that hydrogen dissociation and recombination… ▽ More The atmospheres of ultra-hot Jupiters have been characterized in detail through recent phase curve and low- and high-resolution emission and transmission spectroscopic observations. Previous numerical studies have analyzed the effect of the localized recombination of hydrogen on the atmospheric dynamics and heat transport of ultra-hot Jupiters, finding that hydrogen dissociation and recombination lead to a reduction in the day-to-night contrasts of ultra-hot Jupiters relative to previous expectations. In this work, we add to previous efforts by also considering the localized condensation of clouds in the atmospheres of ultra-hot Jupiters, their resulting transport by the atmospheric circulation, and the radiative feedback of clouds on the atmospheric dynamics. To do so, we include radiatively active cloud tracers into the existing MITgcm framework for simulating the atmospheric dynamics of ultra-hot Jupiters. We take cloud condensate properties appropriate for the high-temperature condensate corundum from CARMA cloud microphysics models. We conduct a suite of GCM simulations with varying cloud microphysical and radiative properties, and we find that partial cloud coverage is a ubiquitous outcome of our simulations. This patchy cloud distribution is inherently set by atmospheric dynamics in addition to equilibrium cloud condensation, and causes a cloud greenhouse effect that warms the atmosphere below the cloud deck. Nightside clouds are further sequestered at depth due to a dynamically induced high-altitude thermal inversion. We post-process our GCMs with the Monte Carlo radiative transfer code gCMCRT and find that the patchy clouds on ultra-hot Jupiters do not significantly impact transmission spectra but can affect their phase-dependent emission spectra. △ Less

Submitted 16 May, 2022; originally announced May 2022.

Comments: 41 pages, 23 figures, 2 tables. Re-submitted to ApJ. Co-first authors

Journal ref: The Astrophysical Journal, 934:79, 2022 July 20

First High-speed Video Camera Observations of a Lightning Flash Associated with a Downward Terrestrial Gamma-ray Flash

Authors: R. U. Abbasi, M. M. F. Saba, J. W. Belz, P. R. Krehbiel, W. Rison, N. Kieu, D. R. da Silva, Dan Rodeheffer, M. A. Stanley, J. Remington, J. Mazich, R. LeVon, K. Smout, A. Petrizze, T. Abu-Zayyad, M. Allen, Y. Arai, R. Arimura, E. Barcikowski, D. R. Bergman, S. A. Blake, I. Buckland, B. G. Cheon, M. Chikawa, T. Fujii , et al. (127 additional authors not shown)

Abstract: In this paper, we present the first high-speed video observation of a cloud-to-ground lightning flash and its associated downward-directed Terrestrial Gamma-ray Flash (TGF). The optical emission of the event was observed by a high-speed video camera running at 40,000 frames per second in conjunction with the Telescope Array Surface Detector, Lightning Mapping Array, interferometer, electric-field… ▽ More In this paper, we present the first high-speed video observation of a cloud-to-ground lightning flash and its associated downward-directed Terrestrial Gamma-ray Flash (TGF). The optical emission of the event was observed by a high-speed video camera running at 40,000 frames per second in conjunction with the Telescope Array Surface Detector, Lightning Mapping Array, interferometer, electric-field fast antenna, and the National Lightning Detection Network. The cloud-to-ground flash associated with the observed TGF was formed by a fast downward leader followed by a very intense return stroke peak current of -154 kA. The TGF occurred while the downward leader was below cloud base, and even when it was halfway in its propagation to ground. The suite of gamma-ray and lightning instruments, timing resolution, and source proximity offer us detailed information and therefore a unique look at the TGF phenomena. △ Less

Submitted 9 August, 2023; v1 submitted 10 May, 2022; originally announced May 2022.

Journal ref: Geophysical Research Letters, 50, e2023GL102958 (2023)

A Mini-Chemical Scheme with Net Reactions for 3D GCMs I.: Thermochemical Kinetics

Authors: Shang-Min Tsai, Elspeth K. H. Lee, Raymond Pierrehumbert

Abstract: Growing evidence has indicated that the global composition distribution plays an indisputable role in interpreting observational data. 3D general circulation models (GCMs) with a reliable treatment of chemistry and clouds are particularly crucial in preparing for the upcoming observations. In the effort of achieving 3D chemistry-climate modeling, the challenge mainly lies in the expensive computin… ▽ More Growing evidence has indicated that the global composition distribution plays an indisputable role in interpreting observational data. 3D general circulation models (GCMs) with a reliable treatment of chemistry and clouds are particularly crucial in preparing for the upcoming observations. In the effort of achieving 3D chemistry-climate modeling, the challenge mainly lies in the expensive computing power required for treating a large number of chemical species and reactions. Motivated by the need for a robust and computationally efficient chemical scheme, we devise a mini-chemical network with a minimal number of species and reactions for H$_2$-dominated atmospheres. We apply a novel technique to simplify the chemical network from a full kinetics model -- VULCAN by replacing a large number of intermediate reactions with net reactions. The number of chemical species is cut down from 67 to 12, with the major species of thermal and observational importance retained, including H$_2$O, CH$_4$, CO, CO$_2$, C$_2$H$_2$, NH$_3$, and HCN. The size of the total reactions is greatly reduced from $\sim$ 800 to 20. The mini-chemical scheme is validated by verifying the temporal evolution and benchmarking the predicted compositions in four exoplanet atmospheres (GJ 1214b, GJ 436b, HD 189733b, HD 209458b) against the full kinetics of VULCAN. It reproduces the chemical timescales and composition distributions of the full kinetics well within an order of magnitude for the major species in the pressure range of 1 bar -- 0.1 mbar across various metallicities and carbon-to-oxygen (C/O) ratios. The small scale of the mini-chemical scheme permits simple use and fast computation, which is optimal for implementation in a 3D GCM or a retrieval framework. We focus on the thermochemical kinetics of net reactions in this paper and address photochemistry in a follow-up paper. △ Less

Submitted 8 April, 2022; originally announced April 2022.

Comments: 9 pages, 5 figures, accepted for publication in A&A

Journal ref: A&A 664, A82 (2022)

Carbon monoxide emission lines reveal an inverted atmosphere in the ultra hot Jupiter WASP-33 b consistent with an eastward hot spot

Authors: Lennart van Sluijs, Jayne L. Birkby, Joshua Lothringer, Elspeth K. H. Lee, Ian J. M. Crossfield, Vivien Parmentier, Matteo Brogi, Craig Kulesa, Don McCarthy, David Charbonneau

Abstract: We report the first detection of CO emission at high spectral resolution in the day-side infrared thermal spectrum of an exoplanet. These emission lines, found in the atmosphere of the transiting ultra hot Jupiter (UHJ) WASP-33 b, provide unambiguous evidence of its thermal inversion. Using spectra from the MMT Exoplanet Atmosphere Survey (MEASURE, $R\sim15,000$), covering pre- and post-eclipse ph… ▽ More We report the first detection of CO emission at high spectral resolution in the day-side infrared thermal spectrum of an exoplanet. These emission lines, found in the atmosphere of the transiting ultra hot Jupiter (UHJ) WASP-33 b, provide unambiguous evidence of its thermal inversion. Using spectra from the MMT Exoplanet Atmosphere Survey (MEASURE, $R\sim15,000$), covering pre- and post-eclipse phases, we cross-correlate with 1D PHOENIX spectral templates to detect CO at S/N = 7.9 ($v_{sys}=0.15^{+0.64}_{-0.65}$ km/s, $K_{p}=229.5^{+1.1}_{-1.0}$ km/s). Moreover, using cross-correlation-to-log-likelihood mapping, we find that the scaling parameter which controls the spectral line contrast changes with phase. We thus use the general circulation model SPARC/MITgcm post-processed by the 3D gCMCRT radiative transfer code to interpret this variation, finding it consistent with an eastward-shifted hot spot. Pre-eclipse, when the hot spot faces Earth, the thermal profiles are shallower leading to smaller line contrast despite greater overall flux. Post-eclipse, the western part of the day-side faces Earth and has much steeper thermal profiles, leading to larger line contrast despite less overall flux. This demonstrates that within the log-likelihood framework, even relatively moderate resolution spectra can be used to understand the 3D nature of close-in exoplanets, and that resolution can be traded for photon-collecting power when the induced Doppler-shift is sufficiently large. We highlight CO as a good probe of UHJ thermal structure and dynamics that does not suffer from stellar activity, unlike species that are also present in the host star e.g. iron lines. △ Less

Submitted 26 April, 2023; v1 submitted 24 March, 2022; originally announced March 2022.

Comments: 24 pages, 21 figures, accepted to MNRAS

Journal ref: MNRAS, Volume 522, Issue 2, June 2023, Pages, 2145-2170

Sunbathing under white light -- 3D modelling of brown dwarf - white dwarf atmospheres with strong UV irradiation

Authors: Elspeth K. H. Lee, Joshua D. Lothringer, Sarah L. Casewell, Daniel Kitzmann, Ben W. P. Lew, Yifan Zhou

Abstract: The atmospheres of brown dwarfs orbiting in close proximity to their parent white dwarf represent some of the most extreme irradiated environments known. Understanding their complex dynamical mechanisms pushes the limits of theoretical and modelling efforts, making them valuable objets to study to test contemporary understanding of irradiated atmospheres. We use the Exo-FMS GCM to simulate the bro… ▽ More The atmospheres of brown dwarfs orbiting in close proximity to their parent white dwarf represent some of the most extreme irradiated environments known. Understanding their complex dynamical mechanisms pushes the limits of theoretical and modelling efforts, making them valuable objets to study to test contemporary understanding of irradiated atmospheres. We use the Exo-FMS GCM to simulate the brown dwarfs WD0137-349B, SDSS J141126.20+200911.1B and EPIC212235321B, first coupled to a multi-banded grey radiative-transfer scheme then a spectral correlated-k scheme with high temperature opacity tables. We then post-process the GCM results using gCMCRT to compare to available observational data. Our GCM models predict strongly temperature inverted atmospheres, spanning many decades in pressure due to impact of UV band heating. Post-processing of our models suggest that the day-night contrast is too small in the GCM results. We therefore suggest that the formation of cloud particles as well as atmospheric drag effects such as magnetic drag are important considerations in setting the day-night temperature contrast for these objects. △ Less

Submitted 18 March, 2022; originally announced March 2022.

Comments: Submitted to MNRAS: 18 March 2020