Showing 1–50 of 185 results for author: Lee, E

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)

First laboratory and on-sky results of an adaptive secondary mirror with TNO-style actuators on the NASA Infrared Telescope Facility

Authors: Ellen Lee, Mark Chun, Olivier Lai, Ruihan Zhang, Max Baeten, Arjo Bos, Matias Kidron, Fred Kamphues, Stefan Kuiper, Wouter Jonker, Michael Connelley, John Rayner, Alan Ryan, Philip Hinz, Rachel Bowens-Rubin, Charles Lockhart, Michael Kelii

Abstract: We are developing an adaptive secondary mirror (ASM) that uses a new actuator technology created by the Netherlands Organization for Applied Scientific Research (TNO). The TNO hybrid variable reluctance actuators have more than an order of magnitude better efficiency over the traditional voice coil actuators that have been used on existing ASMs and show potential for improving the long-term robust… ▽ More We are developing an adaptive secondary mirror (ASM) that uses a new actuator technology created by the Netherlands Organization for Applied Scientific Research (TNO). The TNO hybrid variable reluctance actuators have more than an order of magnitude better efficiency over the traditional voice coil actuators that have been used on existing ASMs and show potential for improving the long-term robustness and reliability of ASMs. To demonstrate the performance, operations, and serviceability of TNO's actuators in an observatory, we have developed a 36-actuator prototype ASM for the NASA Infrared Telescope Facility (IRTF) called IRTF-ASM-1. IRTF-ASM-1 provides the first on-sky demonstration of this approach and will help us evaluate the long-term performance and use of this technology in an astronomical facility environment. We present calibration and performance results with the ASM in a Meniscus Hindle Sphere lens setup as well as preliminary on-sky results on IRTF. IRTF-ASM-1 achieved stable closed-loop performance on-sky with H-band Strehl ratios of 35-40% in long-exposure images under a variety of seeing conditions. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Comments: 15 pages, 12 figures, SPIE Astronomical Telescopes + Instrumentation: Adaptive Optics Systems IX

Probing Dust and Gas Properties Using Ringed Disks

Authors: Eve J. Lee

Abstract: How rapidly a planet grows in mass and how far they may park from the host star depend sensitively on two non-dimensional parameters: Stokes number St and turbulent $α$. Yet, these parameters remain highly uncertain being difficult or impossible to measure directly. Here, we demonstrate how the ringed disks can be leveraged to obtain St and $α$ separately by constructing a simple toy model that co… ▽ More How rapidly a planet grows in mass and how far they may park from the host star depend sensitively on two non-dimensional parameters: Stokes number St and turbulent $α$. Yet, these parameters remain highly uncertain being difficult or impossible to measure directly. Here, we demonstrate how the ringed disks can be leveraged to obtain St and $α$ separately by constructing a simple toy model that combines dust radial equation of motion under aerodynamic drag and coupling to gas motion with the measured distribution of dust masses in Class 0/I disks. Focusing on known systems with well-resolved dust rings, we find that the range of St and $α$ that are consistent with the measured properties of the rings are small: $10^{-4} \lesssim {\rm St} \lesssim 10^{-2}$ and $10^{-5} \lesssim α\lesssim 10^{-3}$. These low St and $α$ ensure the observed rings are stable against clumping. Even in one marginal case where the formation of bound clumps is possible, further mass growth by pebble accretion is inhibited. Furthermore, the derived low $α$ is consistent with the nearly inviscid regime where Type I migration can be prematurely halted. Our analysis predicts minimal planet population beyond $\sim$10s of au where we observe dust rings and significantly more vigorous planet formation inside $\sim$10 AU, consistent with current exo-giant statistics. We close with discussions on the implications of our results on small planet statistics at large orbital distances. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: Accepted to ApJL

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 19 June, 2024; v1 submitted 13 June, 2024; originally announced June 2024.

Comments: 24 pages, 16 figures, resubmitted to PASP (made a few minor changes to the text w.r.t. v1)

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

Giant Outer Transiting Exoplanet Mass (GOT `EM) Survey. V. Two Giant Planets in Kepler-511 but Only One Ran Away

Authors: Yayaati Chachan, Paul A. Dalba, Daniel P. Thorngren, Stephen R. Kane, Eve J. Lee, Edward W. Schwieterman, Howard Isaacson, Andrew W. Howard, Matthew J. Payne

Abstract: Systems hosting multiple giant planets are important laboratories for understanding planetary formation and migration processes. We present a nearly decade-long Doppler spectroscopy campaign from the Keck-I telescope to characterize the two transiting giant planets orbiting Kepler-511 on orbits of 27 days and 297 days. The radial velocity measurements yield precise masses for both planets, which w… ▽ More Systems hosting multiple giant planets are important laboratories for understanding planetary formation and migration processes. We present a nearly decade-long Doppler spectroscopy campaign from the Keck-I telescope to characterize the two transiting giant planets orbiting Kepler-511 on orbits of 27 days and 297 days. The radial velocity measurements yield precise masses for both planets, which we use to infer their bulk heavy element content. Both planets contain approximately 30 Earth masses of heavy elements, but their bulk metallicities (i.e., the ratio between metal mass and total mass) are drastically different ($0.86 \pm 0.04$ and $0.22 \pm 0.04$ respectively). Envelope mass loss cannot account for this difference due to the relatively large orbital distance and mass of the inner planet. We conclude that the outer planet underwent runaway gas accretion while the inner planet did not. This bifurcation in accretion histories is likely a result of the accretion of gas with very different metallicities by the two planets or the late formation of the inner planet from a merger of sub-Neptunes. Kepler-511 uniquely demonstrates how giant planet formation can produce strikingly different outcomes even for planets in the same system. △ Less

Submitted 17 May, 2024; originally announced May 2024.

Comments: submitted to AAS journals, comments welcome

Figuring Out Gas & Galaxies In Enzo (FOGGIE) VIII: Complex and Stochastic Metallicity Gradients at z > 2

Authors: Ayan Acharyya, Molly S. Peeples, Jason Tumlinson, Brian W. O Shea, Cassandra Lochhaas, Anna C. Wright, Raymond C. Simons, Ramona Augustin, Britton D. Smith, Eugene Hyeonmin Lee

Abstract: Gas-phase metallicity gradients are a crucial element in understanding the chemical evolution of galaxies. We use the FOGGIE simulations to study the metallicity gradients ($\nabla Z$) of six Milky Way-like galaxies throughout their evolution. FOGGIE galaxies generally exhibit steep negative gradients for most of their history, with only a few short-lived instances reaching positive slopes that ap… ▽ More Gas-phase metallicity gradients are a crucial element in understanding the chemical evolution of galaxies. We use the FOGGIE simulations to study the metallicity gradients ($\nabla Z$) of six Milky Way-like galaxies throughout their evolution. FOGGIE galaxies generally exhibit steep negative gradients for most of their history, with only a few short-lived instances reaching positive slopes that appear to arise mainly from interactions with other galaxies. FOGGIE concurs with other simulation results but disagrees with the robust observational finding that flat and positive gradients are common at $z>1$. By tracking the metallicity gradient at a rapid cadence of simulation outputs ($\sim 5$--10 Myr), we find that theoretical gradients are highly stochastic: the FOGGIE galaxies spend $\sim 30-50$\% of their time far away from a smoothed trajectory inferred from analytic models or other, less high-cadence simulations. This rapid variation makes instantaneous gradients from observations more difficult to interpret in terms of physical processes. Because of these geometric and stochastic complications, we explore non-parametric methods of quantifying the evolving metallicity distribution at $z > 1$. We investigate how efficiently non-parametric measures of the 2-D metallicity distribution respond to metal production and mixing. Our results suggest that new methods of quantifying and interpreting gas-phase metallicity will be needed to relate trends in upcoming high-$z$ {\it JWST} observations with the underlying physics of gas accretion, expulsion, and recycling in early galaxies. △ Less

Submitted 9 April, 2024; originally announced April 2024.

Comments: Submitted to ApJ

Cosmological inference from combining Planck and ACT cluster counts

Authors: Eunseong Lee, Richard Battye, Boris Bolliet

Abstract: We have adapted the Planck cluster likelihood in such a way that it can be applied to the sample of clusters detected by the Atacama Cosmology Telescope (ACT). Applying it to the 2016 sample from Planck and the 2018 sample from ACT we find, by fixing the cosmology using CMB observations and the cluster model adopted by Planck, that the mass bias required by the two are… ▽ More We have adapted the Planck cluster likelihood in such a way that it can be applied to the sample of clusters detected by the Atacama Cosmology Telescope (ACT). Applying it to the 2016 sample from Planck and the 2018 sample from ACT we find, by fixing the cosmology using CMB observations and the cluster model adopted by Planck, that the mass bias required by the two are $1-b_{\rm Planck}=0.61\pm 0.03$ and $1-b_{\rm ACT}=0.75\pm 0.06$. These are broadly in agreement but hint that the model could be adapted to reach a better agreement. By normalizing the cluster model using weak lensing observations, we find evidence for either evolution in the cluster model, quantified by the cluster modeling parameter describing redshift dependence $β=0.86 \pm 0.07$ using an updated CCCP-based normalization, or evolution in the cosmological model quantified by the dark energy equation of state parameter $w=-0.82 \pm 0.07$. △ Less

Submitted 28 March, 2024; originally announced March 2024.

Comments: 17 pages, 14 figures

The SZ effect with anisotropic distributions and high energy electrons

Authors: Elizabeth Lee, Jens Chluba

Abstract: Future observations of the Sunyaev-Zeldovich (SZ) effect promise ever improving measurements in terms of both sensitivity and angular resolution. As such, it is increasingly relevant to model `higher-order' contributions to the SZ effect. This work examines the effects of high-energy non-thermal electron distributions and those of anisotropic electron and photon distributions on the SZ signals. An… ▽ More Future observations of the Sunyaev-Zeldovich (SZ) effect promise ever improving measurements in terms of both sensitivity and angular resolution. As such, it is increasingly relevant to model `higher-order' contributions to the SZ effect. This work examines the effects of high-energy non-thermal electron distributions and those of anisotropic electron and photon distributions on the SZ signals. Analytic forms of the anisotropic scattering kernels for photons and electrons have been derived and investigated. We present a method for determining the anisotropic contributions through a spherical harmonic decomposition to arbitrary angular multipoles, and discuss the behaviour of these scattering kernels. We then carry out an exploration of various simplistic models of high energy non-thermal electron distributions, and examine their anisotropic behaviour. The kinematic SZ in the relativistic regime is studied using the kernel formulation allowing us to clarifying the role of kinematic corrections to the scattering optical depth. We finally present a release of an updated and refined version of SZpack including a new integrated Python interface and new modules for the calculation of various SZ signals, including those described in this paper. △ Less

Submitted 27 March, 2024; originally announced March 2024.

Comments: 43 pages, 16 figure, comments welcome and to be submitted to JCAP

Zooming by in the CARPoolGP lane: new CAMELS-TNG simulations of zoomed-in massive halos

Authors: Max E. Lee, Shy Genel, Benjamin D. Wandelt, Benjamin Zhang, Ana Maria Delgado, Shivam Pandey, Erwin T. Lau, Christopher Carr, Harrison Cook, Daisuke Nagai, Daniel Angles-Alcazar, Francisco Villaescusa-Navarro, Greg L. Bryan

Abstract: Galaxy formation models within cosmological hydrodynamical simulations contain numerous parameters with non-trivial influences over the resulting properties of simulated cosmic structures and galaxy populations. It is computationally challenging to sample these high dimensional parameter spaces with simulations, particularly for halos in the high-mass end of the mass function. In this work, we dev… ▽ More Galaxy formation models within cosmological hydrodynamical simulations contain numerous parameters with non-trivial influences over the resulting properties of simulated cosmic structures and galaxy populations. It is computationally challenging to sample these high dimensional parameter spaces with simulations, particularly for halos in the high-mass end of the mass function. In this work, we develop a novel sampling and reduced variance regression method, CARPoolGP, which leverages built-in correlations between samples in different locations of high dimensional parameter spaces to provide an efficient way to explore parameter space and generate low variance emulations of summary statistics. We use this method to extend the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) to include a set of 768 zoom-in simulations of halos in the mass range of $10^{13} - 10^{14.5} M_\odot\,h^{-1}$ that span a 28-dimensional parameter space in the IllustrisTNG model. With these simulations and the CARPoolGP emulation method, we explore parameter trends in the Compton $Y-M$, black hole mass-halo mass, and metallicity-mass relations, as well as thermodynamic profiles and quenched fractions of satellite galaxies. We use these emulations to provide a physical picture of the complex interplay between supernova and active galactic nuclei feedback. We then use emulations of the $Y-M$ relation of massive halos to perform Fisher forecasts on astrophysical parameters for future Sunyaev-Zeldovich observations and find a significant improvement in forecasted constraints. We publicly release both the simulation suite and CARPoolGP software package. △ Less

Submitted 15 March, 2024; originally announced March 2024.

Comments: The manuscript was submitted to arxiv after receiving and responding to comments from the first referee report

Friends not Foes: Strong Correlation between Inner Super-Earths and Outer Gas Giants

Authors: Marta L. Bryan, Eve J. Lee

Abstract: The connection between outer gas giants and inner super-Earths reflects their formation and evolutionary histories. Past work exploring this link has suggested a tentative positive correlation between these two populations, but these studies have been limited by small sample sizes and in some cases sample biases. Here we take a new look at this connection with a sample of 184 super-Earth systems w… ▽ More The connection between outer gas giants and inner super-Earths reflects their formation and evolutionary histories. Past work exploring this link has suggested a tentative positive correlation between these two populations, but these studies have been limited by small sample sizes and in some cases sample biases. Here we take a new look at this connection with a sample of 184 super-Earth systems with publicly available radial velocity data and resolved outer gas giants. We calculate the frequency of outer gas giants (GG) in super-Earth (SE) systems, dividing our sample into metal-rich ([Fe/H] $>$ 0) and metal-poor ([Fe/H]$\leq$0) hosts. We find P(GG$|$SE, [Fe/H]$>$0) = 28.0$^{+4.9}_{-4.6}\%$ and P(GG$|$SE, [Fe/H]$\leq$0) = 4.5$^{+2.6}_{-1.9}\%$. Comparing these conditional occurrence rates to field giant occurrence rates from Rosenthal et al 2021, we show that there is a distinct positive correlation between inner super-Earths and outer gas giants for metal-rich host stars at the 2.7$σ$ level, but this correlation disappears for metal-poor systems. We further find that, around metal-rich stars, the GG/SE correlation enhances slightly for systems with giants that are more distant (beyond 3 AU), more eccentric ($e > 0.2$), and/or in multi-gas giant systems. Such trends disappear around metal-poor stars with the exception of systems of multiple giants in which we observe a tentative anti-correlation. Our findings highlight the critical role metallicity (disk solid budget) plays in shaping the overall planetary architecture. △ Less

Submitted 21 May, 2024; v1 submitted 13 March, 2024; originally announced March 2024.

Comments: accepted to ApJL

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

Second gadolinium loading to Super-Kamiokande

Authors: K. Abe, C. Bronner, Y. Hayato, K. Hiraide, K. Hosokawa, K. Ieki, M. Ikeda, J. Kameda, Y. Kanemura, R. Kaneshima, Y. Kashiwagi, Y. Kataoka, S. Miki, S. Mine, M. Miura, S. Moriyama, Y. Nakano, M. Nakahata, S. Nakayama, Y. Noguchi, K. Sato, H. Sekiya, H. Shiba, K. Shimizu, M. Shiozawa , et al. (225 additional authors not shown)

Abstract: The first loading of gadolinium (Gd) into Super-Kamiokande in 2020 was successful, and the neutron capture efficiency on Gd reached 50\%. To further increase the Gd neutron capture efficiency to 75\%, 26.1 tons of $\rm Gd_2(\rm SO_4)_3\cdot \rm 8H_2O$ was additionally loaded into Super-Kamiokande (SK) from May 31 to July 4, 2022. As the amount of loaded $\rm Gd_2(\rm SO_4)_3\cdot \rm 8H_2O$ was do… ▽ More The first loading of gadolinium (Gd) into Super-Kamiokande in 2020 was successful, and the neutron capture efficiency on Gd reached 50\%. To further increase the Gd neutron capture efficiency to 75\%, 26.1 tons of $\rm Gd_2(\rm SO_4)_3\cdot \rm 8H_2O$ was additionally loaded into Super-Kamiokande (SK) from May 31 to July 4, 2022. As the amount of loaded $\rm Gd_2(\rm SO_4)_3\cdot \rm 8H_2O$ was doubled compared to the first loading, the capacity of the powder dissolving system was doubled. We also developed new batches of gadolinium sulfate with even further reduced radioactive impurities. In addition, a more efficient screening method was devised and implemented to evaluate these new batches of $\rm Gd_2(\rm SO_4)_3\cdot \rm 8H_2O$. Following the second loading, the Gd concentration in SK was measured to be $333.5\pm2.5$ ppm via an Atomic Absorption Spectrometer (AAS). From the mean neutron capture time constant of neutrons from an Am/Be calibration source, the Gd concentration was independently measured to be 332.7 $\pm$ 6.8(sys.) $\pm$ 1.1(stat.) ppm, consistent with the AAS result. Furthermore, during the loading the Gd concentration was monitored continually using the capture time constant of each spallation neutron produced by cosmic-ray muons,and the final neutron capture efficiency was shown to become 1.5 times higher than that of the first loaded phase, as expected. △ Less

Submitted 18 June, 2024; v1 submitted 12 March, 2024; originally announced March 2024.

Comments: 34 pages, 13 figures, submitted to Nuclear Inst. and Methods in Physics Research, A

Journal ref: Nuclear Inst. and Methods in Physics Research, A 1065 (2024) 169480

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)

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

The 2022-2023 accretion outburst of the young star V1741 Sgr

Authors: Michael A. Kuhn, Lynne A. Hillenbrand, Michael S. Connelley, R. Michael Rich, Bart Staels, Adolfo S. Carvalho, Philip W. Lucas, Christoffer Fremling, Viraj R. Karambelkar, Ellen Lee, Tomás Ahumada, Emille E. O. Ishida, Kishalay De, Rafael S. de Souza, Mansi Kasliwal

Abstract: V1741 Sgr (= SPICY 71482/Gaia22dtk) is a Classical T Tauri star on the outskirts of the Lagoon Nebula. After at least a decade of stability, in mid-2022, the optical source brightened by ~3 mag over two months, remained bright until early 2023, then dimmed erratically over the next four months. This event was monitored with optical and infrared spectroscopy and photometry. Spectra from the peak (O… ▽ More V1741 Sgr (= SPICY 71482/Gaia22dtk) is a Classical T Tauri star on the outskirts of the Lagoon Nebula. After at least a decade of stability, in mid-2022, the optical source brightened by ~3 mag over two months, remained bright until early 2023, then dimmed erratically over the next four months. This event was monitored with optical and infrared spectroscopy and photometry. Spectra from the peak (October 2022) indicate an EX Lup-type (EXor) accretion outburst, with strong emission from H I, He I, and Ca II lines and CO bands. At this stage, spectroscopic absorption features indicated a temperature of T ~ 4750 K with low-gravity lines (e.g., Ba II and Sr II). By April 2023, with the outburst beginning to dim, strong TiO absorption appeared, indicating a cooler T ~ 3600 K temperature. However, once the source had returned to its pre-outburst flux in August 2023, the TiO absorption and the CO emission disappeared. When the star went into outburst, the source's spectral energy distribution became flatter, leading to bluer colours at wavelengths shorter than ~1.6 microns and redder colours at longer wavelengths. The brightening requires a continuum emitting area larger than the stellar surface, likely from optically thick circumstellar gas with cooler surface layers producing the absorption features. Additional contributions to the outburst spectrum may include blue excess from hotspots on the stellar surface, emission lines from diffuse gas, and reprocessed emission from the dust disc. Cooling of the circumstellar gas would explain the appearance of TiO, which subsequently disappeared once this gas had faded and the stellar spectrum reemerged. △ Less

Submitted 17 January, 2024; originally announced January 2024.

Comments: Accepted for publication in MNRAS; 17 pages, 16 figures, and 2 tables

New Eruptive YSOs from SPICY and WISE

Authors: C. Contreras Peña, M. Ashraf, J. E. Lee, G. Herczeg, P. W. Lucas, Z. Guo, D. Johnstone, H. G. Lee, J. Jose

Abstract: This work presents four high-amplitude variable YSOs ($\simeq$ 3 mag at near- or mid-IR wavelengths) arising from the SPICY catalog. Three outbursts show a duration that is longer than 1 year, and are still ongoing. And additional YSO brightened over the last two epochs of NEOWISE observations and the duration of the outburst is thus unclear. Analysis of the spectra of the four sources confirms th… ▽ More This work presents four high-amplitude variable YSOs ($\simeq$ 3 mag at near- or mid-IR wavelengths) arising from the SPICY catalog. Three outbursts show a duration that is longer than 1 year, and are still ongoing. And additional YSO brightened over the last two epochs of NEOWISE observations and the duration of the outburst is thus unclear. Analysis of the spectra of the four sources confirms them as new members of the eruptive variable class. We find two YSOs that can be firmly classified as bona fide FUors and one object that falls in the V1647 Ori-like class. Given the uncertainty in the duration of its outburst, an additional YSO can only be classified as a candidate FUor. Continued monitoring and follow-up of these particular sources is important to better understand the accretion process of YSOs. △ Less

Submitted 29 November, 2023; originally announced November 2023.

Comments: 10 pages, 10 figures, Accepted for publication at the Journal of the Korean Astronomical Society

Journal ref: JKAS 2023, 56, 2, 253

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 Uniform Analysis of Debris Disks with the Gemini Planet Imager I: An Empirical Search for Perturbations from Planetary Companions in Polarized Light Images

Authors: Katie A. Crotts, Brenda C. Matthews, Gaspard Duchêne, Thomas M. Esposito, Ruobing Dong, Justin Hom, Rebecca Oppenheimer, Malena Rice, Schuyler G. Wolff, Christine H. Chen, Clarissa R. Do Ó, Paul Kalas, Briley L. Lewis, Alycia J. Weinberger, David J. Wilner, Mark Ammons, Pauline Arriaga, Robert J. De Rosa, John H. Debes, Michael P. Fitzgerald, Eileen C. Gonzales, Dean C. Hines, Sasha Hinkley, A. Meredith Hughes, Ludmilla Kolokolova , et al. (15 additional authors not shown)

Abstract: The Gemini Planet Imager (GPI) has excelled in imaging debris disks in the near-infrared. The GPI Exoplanet Survey (GPIES) imaged twenty-four debris disks in polarized $H$-band light, while other programs observed half of these disks in polarized $J$- and/or $K1$-bands. Using these data, we present a uniform analysis of the morphology of each disk to find asymmetries suggestive of perturbations, p… ▽ More The Gemini Planet Imager (GPI) has excelled in imaging debris disks in the near-infrared. The GPI Exoplanet Survey (GPIES) imaged twenty-four debris disks in polarized $H$-band light, while other programs observed half of these disks in polarized $J$- and/or $K1$-bands. Using these data, we present a uniform analysis of the morphology of each disk to find asymmetries suggestive of perturbations, particularly those due to planet-disk interactions. The multi-wavelength surface brightness, the disk color and geometry permit identification of any asymmetries such as warps or disk offsets from the central star. We find that nineteen of the disks in this sample exhibit asymmetries in surface brightness, disk color, disk geometry, or a combination of the three, suggesting that for this sample, perturbations, as seen in scattered light, are common. The relationship between these perturbations and potential planets in the system are discussed. We also explore correlations among stellar temperatures, ages, disk properties, and observed perturbations. We find significant trends between the vertical aspect ratio and the stellar temperature, disk radial extent, and the dust grain size distribution power-law, $q$. We also confirm a trend between the disk color and stellar effective temperature, where the disk becomes increasingly red/neutral with increasing temperature. Such results have important implications on the evolution of debris disk systems around stars of various spectral types. △ Less

Submitted 24 November, 2023; originally announced November 2023.

Comments: 46 pages, 20 figures, 6 tables, accepted for publication in ApJ

Alpha backgrounds in NaI(Tl) crystals of COSINE-100

Authors: G. Adhikari, N. Carlin, D. F. F. S. Cavalcante, J. Y. Cho, J. J. Choi, S. Choi, A. C. Ezeribe, L. E. Franca, C. Ha, I. S. Hahn, S. J. Hollick, E. J. Jeon, H. W. Joo, W. G. Kang, M. Kauer, B. H. Kim, H. J. Kim, J. Kim, K. W. Kim, S. H. Kim, S. K. Kim, S. W. Kim, W. K. Kim, Y. D. Kim, Y. H. Kim , et al. (38 additional authors not shown)

Abstract: COSINE-100 is a dark matter direct detection experiment with 106 kg NaI(Tl) as the target material. 210Pb and daughter isotopes are a dominant background in the WIMP region of interest and are detected via beta decay and alpha decay. Analysis of the alpha channel complements the background model as observed in the beta/gamma channel. We present the measurement of the quenching factors and Monte Ca… ▽ More COSINE-100 is a dark matter direct detection experiment with 106 kg NaI(Tl) as the target material. 210Pb and daughter isotopes are a dominant background in the WIMP region of interest and are detected via beta decay and alpha decay. Analysis of the alpha channel complements the background model as observed in the beta/gamma channel. We present the measurement of the quenching factors and Monte Carlo simulation results and activity quantification of the alpha decay components of the COSINE-100 NaI(Tl) crystals. The data strongly indicate that the alpha decays probabilistically undergo two possible quenching factors but require further investigation. The fitted results are consistent with independent measurements and improve the overall understanding of the COSINE-100 backgrounds. Furthermore, the half-life of 216Po has been measured to be 143.4 +/- 1.2 ms, which is consistent with and more precise than recent measurements. △ Less

Submitted 30 January, 2024; v1 submitted 8 November, 2023; originally announced November 2023.

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

Calibrating approximate Bayesian credible intervals of gravitational-wave parameters

Authors: Ruiting Mao, Jeong Eun Lee, Ollie Burke, Alvin J. K. Chua, Matthew C. Edwards, Renate Meyer

Abstract: Approximations are commonly employed in realistic applications of scientific Bayesian inference, often due to convenience if not necessity. In the field of gravitational-wave (GW) data analysis, fast-to-evaluate but approximate waveform models of astrophysical GW signals are sometimes used in lieu of more accurate models to infer properties of a true GW signal buried within detector noise. In addi… ▽ More Approximations are commonly employed in realistic applications of scientific Bayesian inference, often due to convenience if not necessity. In the field of gravitational-wave (GW) data analysis, fast-to-evaluate but approximate waveform models of astrophysical GW signals are sometimes used in lieu of more accurate models to infer properties of a true GW signal buried within detector noise. In addition, a Fisher-information-based normal approximation to the posterior distribution can also be used to conduct inference in bulk, without the need for extensive numerical calculations such as Markov chain Monte Carlo (MCMC) simulations. Such approximations can generally lead to an inaccurate posterior distribution with poor statistical coverage of the true posterior. In this article, we present a novel calibration procedure that calibrates the credible sets for a family of approximate posterior distributions, to ensure coverage of the true posterior at a level specified by the analyst. Tools such as autoencoders and artificial neural networks are used within our calibration model to compress the data (for efficiency) and to perform tasks such as logistic regression. As a proof of principle, we demonstrate our formalism on the GW signal from a high-mass binary black hole merger, a promising source for the near-future space-based GW observatory LISA. △ Less

Submitted 30 January, 2024; v1 submitted 10 October, 2023; originally announced October 2023.

Comments: 24 pages, 12 figures

Journal ref: Phys. Rev. D 109, 083002 (2024)

The Not-So Dramatic Effect of Advective Flows on Gas Accretion

Authors: Vincent Savignac, Eve J. Lee

Abstract: Super-Earths and mini-Neptunes are the most common types of exoplanets discovered, yet the physics of their formation are still debated. Standard core accretion models in gas-rich environment find that typical mini-Neptune mass planets would blow up into Jupiters before the underlying disk gas dissipates away. The injection of entropy from the protoplanetary disk into forming gaseous envelopes has… ▽ More Super-Earths and mini-Neptunes are the most common types of exoplanets discovered, yet the physics of their formation are still debated. Standard core accretion models in gas-rich environment find that typical mini-Neptune mass planets would blow up into Jupiters before the underlying disk gas dissipates away. The injection of entropy from the protoplanetary disk into forming gaseous envelopes has recently been put forward as a mechanism to delay this runaway accretion, specifically at short orbital distances. Here, we reevaluate this line of reasoning by incorporating recycling flows of gas into a numerical one-dimensional thermodynamic model with more realistic equation of state and opacities and the thermal state of the advective flow. At 0.1 AU, we find that advective flows are only able to produce mini-Neptunes if they can penetrate below ~0.25 of the planet's gravitational sphere of influence. Otherwise, the gas-to-core mass ratio (GCR) reaches above ~10% which is too large to explain the measured properties of mini-Neptunes, necessitating other gas-limiting processes such as late-time core assembly. The effect of entropy advection on gas accretion weakens even further beyond 0.1 AU. We present an updated scaling relation between GCR and the penetration depth of the advective flows which varies non-trivially with orbital distances, core masses and dusty vs. dust-free opacity. We further demonstrate how measurements of planet mass distribution beyond ~1 AU using future instruments such as the Nancy Grace Roman Space Telescope could be used to disambiguate between different formation conditions of gas-poor planets. △ Less

Submitted 9 July, 2024; v1 submitted 9 October, 2023; originally announced October 2023.

Comments: 18 pages, 9 figures. Accepted to ApJ

A Joint Fermi-GBM and Swift-BAT Analysis of Gravitational-Wave Candidates from the Third Gravitational-wave Observing Run

Authors: C. Fletcher, J. Wood, R. Hamburg, P. Veres, C. M. Hui, E. Bissaldi, M. S. Briggs, E. Burns, W. H. Cleveland, M. M. Giles, A. Goldstein, B. A. Hristov, D. Kocevski, S. Lesage, B. Mailyan, C. Malacaria, S. Poolakkil, A. von Kienlin, C. A. Wilson-Hodge, The Fermi Gamma-ray Burst Monitor Team, M. Crnogorčević, J. DeLaunay, A. Tohuvavohu, R. Caputo, S. B. Cenko , et al. (1674 additional authors not shown)

Abstract: We present Fermi Gamma-ray Burst Monitor (Fermi-GBM) and Swift Burst Alert Telescope (Swift-BAT) searches for gamma-ray/X-ray counterparts to gravitational wave (GW) candidate events identified during the third observing run of the Advanced LIGO and Advanced Virgo detectors. Using Fermi-GBM on-board triggers and sub-threshold gamma-ray burst (GRB) candidates found in the Fermi-GBM ground analyses,… ▽ More We present Fermi Gamma-ray Burst Monitor (Fermi-GBM) and Swift Burst Alert Telescope (Swift-BAT) searches for gamma-ray/X-ray counterparts to gravitational wave (GW) candidate events identified during the third observing run of the Advanced LIGO and Advanced Virgo detectors. Using Fermi-GBM on-board triggers and sub-threshold gamma-ray burst (GRB) candidates found in the Fermi-GBM ground analyses, the Targeted Search and the Untargeted Search, we investigate whether there are any coincident GRBs associated with the GWs. We also search the Swift-BAT rate data around the GW times to determine whether a GRB counterpart is present. No counterparts are found. Using both the Fermi-GBM Targeted Search and the Swift-BAT search, we calculate flux upper limits and present joint upper limits on the gamma-ray luminosity of each GW. Given these limits, we constrain theoretical models for the emission of gamma-rays from binary black hole mergers. △ Less

Submitted 25 August, 2023; originally announced August 2023.

Identification of the Top TESS Objects of Interest for Atmospheric Characterization of Transiting Exoplanets with JWST

Authors: Benjamin J. Hord, Eliza M. -R. Kempton, Thomas Mikal-Evans, David W. Latham, David R. Ciardi, Diana Dragomir, Knicole D. Colón, Gabrielle Ross, Andrew Vanderburg, Zoe L. de Beurs, Karen A. Collins, Cristilyn N. Watkins, Jacob Bean, Nicolas B. Cowan, Tansu Daylan, Caroline V. Morley, Jegug Ih, David Baker, Khalid Barkaoui, Natalie M. Batalha, Aida Behmard, Alexander Belinski, Zouhair Benkhaldoun, Paul Benni, Krzysztof Bernacki , et al. (120 additional authors not shown)

Abstract: JWST has ushered in an era of unprecedented ability to characterize exoplanetary atmospheres. While there are over 5,000 confirmed planets, more than 4,000 TESS planet candidates are still unconfirmed and many of the best planets for atmospheric characterization may remain to be identified. We present a sample of TESS planets and planet candidates that we identify as "best-in-class" for transmissi… ▽ More JWST has ushered in an era of unprecedented ability to characterize exoplanetary atmospheres. While there are over 5,000 confirmed planets, more than 4,000 TESS planet candidates are still unconfirmed and many of the best planets for atmospheric characterization may remain to be identified. We present a sample of TESS planets and planet candidates that we identify as "best-in-class" for transmission and emission spectroscopy with JWST. These targets are sorted into bins across equilibrium temperature $T_{\mathrm{eq}}$ and planetary radius $R{_\mathrm{p}}$ and are ranked by transmission and emission spectroscopy metric (TSM and ESM, respectively) within each bin. In forming our target sample, we perform cuts for expected signal size and stellar brightness, to remove sub-optimal targets for JWST. Of the 194 targets in the resulting sample, 103 are unconfirmed TESS planet candidates, also known as TESS Objects of Interest (TOIs). We perform vetting and statistical validation analyses on these 103 targets to determine which are likely planets and which are likely false positives, incorporating ground-based follow-up from the TESS Follow-up Observation Program (TFOP) to aid the vetting and validation process. We statistically validate 23 TOIs, marginally validate 33 TOIs to varying levels of confidence, deem 29 TOIs likely false positives, and leave the dispositions for 4 TOIs as inconclusive. 14 of the 103 TOIs were confirmed independently over the course of our analysis. We provide our final best-in-class sample as a community resource for future JWST proposals and observations. We intend for this work to motivate formal confirmation and mass measurements of each validated planet and encourage more detailed analysis of individual targets by the community. △ Less

Submitted 18 August, 2023; originally announced August 2023.

Comments: Submitted to AJ. Machine-readable versions of Tables 2 and 3 are included. 40 pages, 7 figures, 3 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

The Atacama Cosmology Telescope: High-resolution component-separated maps across one-third of the sky

Authors: William R. Coulton, Mathew S. Madhavacheril, Adriaan J. Duivenvoorden, J. Colin Hill, Irene Abril-Cabezas, Peter A. R. Ade, Simone Aiola, Tommy Alford, Mandana Amiri, Stefania Amodeo, Rui An, Zachary Atkins, Jason E. Austermann, Nicholas Battaglia, Elia Stefano Battistelli, James A. Beall, Rachel Bean, Benjamin Beringue, Tanay Bhandarkar, Emily Biermann, Boris Bolliet, J Richard Bond, Hongbo Cai, Erminia Calabrese, Victoria Calafut , et al. (129 additional authors not shown)

Abstract: Observations of the millimeter sky contain valuable information on a number of signals, including the blackbody cosmic microwave background (CMB), Galactic emissions, and the Compton-$y$ distortion due to the thermal Sunyaev-Zel'dovich (tSZ) effect. Extracting new insight into cosmological and astrophysical questions often requires combining multi-wavelength observations to spectrally isolate one… ▽ More Observations of the millimeter sky contain valuable information on a number of signals, including the blackbody cosmic microwave background (CMB), Galactic emissions, and the Compton-$y$ distortion due to the thermal Sunyaev-Zel'dovich (tSZ) effect. Extracting new insight into cosmological and astrophysical questions often requires combining multi-wavelength observations to spectrally isolate one component. In this work, we present a new arcminute-resolution Compton-$y$ map, which traces out the line-of-sight-integrated electron pressure, as well as maps of the CMB in intensity and E-mode polarization, across a third of the sky (around 13,000 sq.~deg.). We produce these through a joint analysis of data from the Atacama Cosmology Telescope (ACT) Data Release 4 and 6 at frequencies of roughly 93, 148, and 225 GHz, together with data from the \textit{Planck} satellite at frequencies between 30 GHz and 545 GHz. We present detailed verification of an internal linear combination pipeline implemented in a needlet frame that allows us to efficiently suppress Galactic contamination and account for spatial variations in the ACT instrument noise. These maps provide a significant advance, in noise levels and resolution, over the existing \textit{Planck} component-separated maps and will enable a host of science goals including studies of cluster and galaxy astrophysics, inferences of the cosmic velocity field, primordial non-Gaussianity searches, and gravitational lensing reconstruction of the CMB. △ Less

Submitted 3 July, 2023; originally announced July 2023.

Comments: The Compton-y map and associated products will be made publicly available upon publication of the paper. The CMB T and E mode maps will be made available when the DR6 maps are made public

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)

Coronal Heating as Determined by the Solar Flare Frequency Distribution Obtained by Aggregating Case Studies

Authors: James Paul Mason, Alexandra Werth, Colin G. West, Allison A. Youngblood, Donald L. Woodraska, Courtney Peck, Kevin Lacjak, Florian G. Frick, Moutamen Gabir, Reema A. Alsinan, Thomas Jacobsen, Mohammad Alrubaie, Kayla M. Chizmar, Benjamin P. Lau, Lizbeth Montoya Dominguez, David Price, Dylan R. Butler, Connor J. Biron, Nikita Feoktistov, Kai Dewey, N. E. Loomis, Michal Bodzianowski, Connor Kuybus, Henry Dietrick, Aubrey M. Wolfe , et al. (977 additional authors not shown)

Abstract: Flare frequency distributions represent a key approach to addressing one of the largest problems in solar and stellar physics: determining the mechanism that counter-intuitively heats coronae to temperatures that are orders of magnitude hotter than the corresponding photospheres. It is widely accepted that the magnetic field is responsible for the heating, but there are two competing mechanisms th… ▽ More Flare frequency distributions represent a key approach to addressing one of the largest problems in solar and stellar physics: determining the mechanism that counter-intuitively heats coronae to temperatures that are orders of magnitude hotter than the corresponding photospheres. It is widely accepted that the magnetic field is responsible for the heating, but there are two competing mechanisms that could explain it: nanoflares or Alfvén waves. To date, neither can be directly observed. Nanoflares are, by definition, extremely small, but their aggregate energy release could represent a substantial heating mechanism, presuming they are sufficiently abundant. One way to test this presumption is via the flare frequency distribution, which describes how often flares of various energies occur. If the slope of the power law fitting the flare frequency distribution is above a critical threshold, $α=2$ as established in prior literature, then there should be a sufficient abundance of nanoflares to explain coronal heating. We performed $>$600 case studies of solar flares, made possible by an unprecedented number of data analysts via three semesters of an undergraduate physics laboratory course. This allowed us to include two crucial, but nontrivial, analysis methods: pre-flare baseline subtraction and computation of the flare energy, which requires determining flare start and stop times. We aggregated the results of these analyses into a statistical study to determine that $α= 1.63 \pm 0.03$. This is below the critical threshold, suggesting that Alfvén waves are an important driver of coronal heating. △ Less

Submitted 9 May, 2023; originally announced May 2023.

Comments: 1,002 authors, 14 pages, 4 figures, 3 tables, published by The Astrophysical Journal on 2023-05-09, volume 948, page 71

Small Planets Around Cool Dwarfs: Enhanced Formation Efficiency of Super-Earths around M dwarfs

Authors: Yayaati Chachan, Eve J. Lee

Abstract: Current measurements of planet population as a function of stellar mass show three seemingly contradictory signatures: close-in super-Earths are more prevalent around M dwarfs than FGK dwarfs; inner super-Earths are correlated with outer giants; and outer giants are less common around M dwarfs than FGK dwarfs. Here, we build a simple framework that combines the theory of pebble accretion with the… ▽ More Current measurements of planet population as a function of stellar mass show three seemingly contradictory signatures: close-in super-Earths are more prevalent around M dwarfs than FGK dwarfs; inner super-Earths are correlated with outer giants; and outer giants are less common around M dwarfs than FGK dwarfs. Here, we build a simple framework that combines the theory of pebble accretion with the measurements of dust masses in protoplanetary disks to reconcile all three observations. First, we show that cooler stars are more efficient at converting pebbles into planetary cores at short orbital periods. Second, when disks are massive enough to nucleate a heavy core at 5 AU, more than enough dust can drift in to assemble inner planets, establishing the correlation between inner planets and outer giants. Finally, while stars of varying masses are similarly capable of converting pebbles into cores at long orbital periods, hotter stars are much more likely to harbor more massive dust disks so that the giant planet occurrence rate rises around hotter stars. Our results are valid over a wide range of parameter space for a disk accretion rate that follows $\dot{M}_\star \sim 10^{-8}\,M_\odot\,{\rm yr}^{-1}(M_\star/M_\odot)^2$. We predict a decline in mini-Neptune population (but not necessarily terrestrial planets) around stars lighter than $\sim 0.3-0.5 M_\odot$. Cold giants ($\gtrsim$5 AU), if they exist, should remain correlated with inner planets even around lower mass stars. △ Less

Submitted 28 June, 2023; v1 submitted 1 May, 2023; originally announced May 2023.

Comments: accepted to ApJL

Search for gravitational-lensing signatures in the full third observing run of the LIGO-Virgo network

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, C. Alléné, A. Allocca, P. A. Altin , et al. (1670 additional authors not shown)

Abstract: Gravitational lensing by massive objects along the line of sight to the source causes distortions of gravitational wave-signals; such distortions may reveal information about fundamental physics, cosmology and astrophysics. In this work, we have extended the search for lensing signatures to all binary black hole events from the third observing run of the LIGO--Virgo network. We search for repeated… ▽ More Gravitational lensing by massive objects along the line of sight to the source causes distortions of gravitational wave-signals; such distortions may reveal information about fundamental physics, cosmology and astrophysics. In this work, we have extended the search for lensing signatures to all binary black hole events from the third observing run of the LIGO--Virgo network. We search for repeated signals from strong lensing by 1) performing targeted searches for subthreshold signals, 2) calculating the degree of overlap amongst the intrinsic parameters and sky location of pairs of signals, 3) comparing the similarities of the spectrograms amongst pairs of signals, and 4) performing dual-signal Bayesian analysis that takes into account selection effects and astrophysical knowledge. We also search for distortions to the gravitational waveform caused by 1) frequency-independent phase shifts in strongly lensed images, and 2) frequency-dependent modulation of the amplitude and phase due to point masses. None of these searches yields significant evidence for lensing. Finally, we use the non-detection of gravitational-wave lensing to constrain the lensing rate based on the latest merger-rate estimates and the fraction of dark matter composed of compact objects. △ Less

Submitted 17 April, 2023; originally announced April 2023.

Comments: 28 pages, 11 figures

Report number: LIGO-P2200031

The Atacama Cosmology Telescope: DR6 Gravitational Lensing Map and Cosmological Parameters

Authors: Mathew S. Madhavacheril, Frank J. Qu, Blake D. Sherwin, Niall MacCrann, Yaqiong Li, Irene Abril-Cabezas, Peter A. R. Ade, Simone Aiola, Tommy Alford, Mandana Amiri, Stefania Amodeo, Rui An, Zachary Atkins, Jason E. Austermann, Nicholas Battaglia, Elia Stefano Battistelli, James A. Beall, Rachel Bean, Benjamin Beringue, Tanay Bhandarkar, Emily Biermann, Boris Bolliet, J Richard Bond, Hongbo Cai, Erminia Calabrese , et al. (134 additional authors not shown)

Abstract: We present cosmological constraints from a gravitational lensing mass map covering 9400 sq. deg. reconstructed from CMB measurements made by the Atacama Cosmology Telescope (ACT) from 2017 to 2021. In combination with BAO measurements (from SDSS and 6dF), we obtain the amplitude of matter fluctuations $σ_8 = 0.819 \pm 0.015$ at 1.8% precision, $S_8\equivσ_8({Ω_{\rm m}}/0.3)^{0.5}=0.840\pm0.028$ an… ▽ More We present cosmological constraints from a gravitational lensing mass map covering 9400 sq. deg. reconstructed from CMB measurements made by the Atacama Cosmology Telescope (ACT) from 2017 to 2021. In combination with BAO measurements (from SDSS and 6dF), we obtain the amplitude of matter fluctuations $σ_8 = 0.819 \pm 0.015$ at 1.8% precision, $S_8\equivσ_8({Ω_{\rm m}}/0.3)^{0.5}=0.840\pm0.028$ and the Hubble constant $H_0= (68.3 \pm 1.1)\, \text{km}\,\text{s}^{-1}\,\text{Mpc}^{-1}$ at 1.6% precision. A joint constraint with CMB lensing measured by the Planck satellite yields even more precise values: $σ_8 = 0.812 \pm 0.013$, $S_8\equivσ_8({Ω_{\rm m}}/0.3)^{0.5}=0.831\pm0.023$ and $H_0= (68.1 \pm 1.0)\, \text{km}\,\text{s}^{-1}\,\text{Mpc}^{-1}$. These measurements agree well with $Λ$CDM-model extrapolations from the CMB anisotropies measured by Planck. To compare these constraints to those from the KiDS, DES, and HSC galaxy surveys, we revisit those data sets with a uniform set of assumptions, and find $S_8$ from all three surveys are lower than that from ACT+Planck lensing by varying levels ranging from 1.7-2.1$σ$. These results motivate further measurements and comparison, not just between the CMB anisotropies and galaxy lensing, but also between CMB lensing probing $z\sim 0.5-5$ on mostly-linear scales and galaxy lensing at $z\sim 0.5$ on smaller scales. We combine our CMB lensing measurements with CMB anisotropies to constrain extensions of $Λ$CDM, limiting the sum of the neutrino masses to $\sum m_ν < 0.12$ eV (95% c.l.), for example. Our results provide independent confirmation that the universe is spatially flat, conforms with general relativity, and is described remarkably well by the $Λ$CDM model, while paving a promising path for neutrino physics with gravitational lensing from upcoming ground-based CMB surveys. △ Less

Submitted 11 April, 2023; originally announced April 2023.

Comments: 30 pages, 16 figures, prepared for submission to ApJ. Cosmological likelihood data is here: https://lambda.gsfc.nasa.gov/product/act/actadv_prod_table.html ; likelihood software is here: https://github.com/ACTCollaboration/act_dr6_lenslike . Also see companion papers Qu et al and MacCrann et al. Mass maps will be released when papers are published

The Atacama Cosmology Telescope: A Measurement of the DR6 CMB Lensing Power Spectrum and its Implications for Structure Growth

Authors: Frank J. Qu, Blake D. Sherwin, Mathew S. Madhavacheril, Dongwon Han, Kevin T. Crowley, Irene Abril-Cabezas, Peter A. R. Ade, Simone Aiola, Tommy Alford, Mandana Amiri, Stefania Amodeo, Rui An, Zachary Atkins, Jason E. Austermann, Nicholas Battaglia, Elia Stefano Battistelli, James A. Beall, Rachel Bean, Benjamin Beringue, Tanay Bhandarkar, Emily Biermann, Boris Bolliet, J Richard Bond, Hongbo Cai, Erminia Calabrese , et al. (133 additional authors not shown)

Abstract: We present new measurements of cosmic microwave background (CMB) lensing over $9400$ sq. deg. of the sky. These lensing measurements are derived from the Atacama Cosmology Telescope (ACT) Data Release 6 (DR6) CMB dataset, which consists of five seasons of ACT CMB temperature and polarization observations. We determine the amplitude of the CMB lensing power spectrum at $2.3\%$ precision ($43σ$ sign… ▽ More We present new measurements of cosmic microwave background (CMB) lensing over $9400$ sq. deg. of the sky. These lensing measurements are derived from the Atacama Cosmology Telescope (ACT) Data Release 6 (DR6) CMB dataset, which consists of five seasons of ACT CMB temperature and polarization observations. We determine the amplitude of the CMB lensing power spectrum at $2.3\%$ precision ($43σ$ significance) using a novel pipeline that minimizes sensitivity to foregrounds and to noise properties. To ensure our results are robust, we analyze an extensive set of null tests, consistency tests, and systematic error estimates and employ a blinded analysis framework. The baseline spectrum is well fit by a lensing amplitude of $A_{\mathrm{lens}}=1.013\pm0.023$ relative to the Planck 2018 CMB power spectra best-fit $Λ$CDM model and $A_{\mathrm{lens}}=1.005\pm0.023$ relative to the $\text{ACT DR4} + \text{WMAP}$ best-fit model. From our lensing power spectrum measurement, we derive constraints on the parameter combination $S^{\mathrm{CMBL}}_8 \equiv σ_8 \left({Ω_m}/{0.3}\right)^{0.25}$ of $S^{\mathrm{CMBL}}_8= 0.818\pm0.022$ from ACT DR6 CMB lensing alone and $S^{\mathrm{CMBL}}_8= 0.813\pm0.018$ when combining ACT DR6 and Planck NPIPE CMB lensing power spectra. These results are in excellent agreement with $Λ$CDM model constraints from Planck or $\text{ACT DR4} + \text{WMAP}$ CMB power spectrum measurements. Our lensing measurements from redshifts $z\sim0.5$--$5$ are thus fully consistent with $Λ$CDM structure growth predictions based on CMB anisotropies probing primarily $z\sim1100$. We find no evidence for a suppression of the amplitude of cosmic structure at low redshifts △ Less

Submitted 28 May, 2024; v1 submitted 11 April, 2023; originally announced April 2023.

Comments: 45+22 pages, 50 figures. v2 matches with published version in ApJ. Cosmological likelihood data and lensing maps are here: https://lambda.gsfc.nasa.gov/product/act/actadv_prod_table.html ; likelihood software is here: https://github.com/ACTCollaboration/act_dr6_lenslike . Also see companion papers Madhavacheril et al and MacCrann et al

Report number: FERMILAB-PUB-23-237-PPD

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

Mass production of ultra-pure NaI powder for COSINE-200

Authors: KeonAh Shin, JunSeok Choe, Olga Gileva, Alain Iltis, Yena Kim, Yeongduk Kim, Cheolho Lee, Eunkyung Lee, HyunSu Lee, Moo Hyun Lee

Abstract: COSINE-200 is the next phase experiment of the ongoing COSINE-100 that aims to unambiguously verify the annual modulation signals observed by the DAMA experiment and to reach the world competitive sensitivity on the low-mass dark matter search. To achieve the physics goal of the COSINE-200, the successful production of the low-background NaI(Tl) detectors is crucial and it must begin from the mass… ▽ More COSINE-200 is the next phase experiment of the ongoing COSINE-100 that aims to unambiguously verify the annual modulation signals observed by the DAMA experiment and to reach the world competitive sensitivity on the low-mass dark matter search. To achieve the physics goal of the COSINE-200, the successful production of the low-background NaI(Tl) detectors is crucial and it must begin from the mass production of the ultra-low background NaI powder. A clean facility for mass-producing the pure-NaI powder has been constructed at the Center for Underground Physics (CUP) in Korea. Two years of operation determined efficient parameters of the mass purification and provided a total of 480 kg of the ultra-pure NaI powder in hand. The potassium concentration in the produced powders varied from 5.4 to 11 ppb, and the maximum production capacity of 35 kg per two weeks was achieved. Here, we report our operational practice with the mass purification of the NaI powder, which includes raw powder purification, recycling of the mother solution, and recovery of NaI from the residual melt that remained after crystal growth. △ Less

Submitted 13 January, 2023; originally announced January 2023.

Search for subsolar-mass black hole binaries in the second part of Advanced LIGO's and Advanced Virgo's third observing run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, C. Alléné, A. Allocca, P. A. Altin , et al. (1680 additional authors not shown)

Abstract: We describe a search for gravitational waves from compact binaries with at least one component with mass 0.2 $M_\odot$ -- $1.0 M_\odot$ and mass ratio $q \geq 0.1$ in Advanced LIGO and Advanced Virgo data collected between 1 November 2019, 15:00 UTC and 27 March 2020, 17:00 UTC. No signals were detected. The most significant candidate has a false alarm rate of 0.2 $\mathrm{yr}^{-1}$. We estimate t… ▽ More We describe a search for gravitational waves from compact binaries with at least one component with mass 0.2 $M_\odot$ -- $1.0 M_\odot$ and mass ratio $q \geq 0.1$ in Advanced LIGO and Advanced Virgo data collected between 1 November 2019, 15:00 UTC and 27 March 2020, 17:00 UTC. No signals were detected. The most significant candidate has a false alarm rate of 0.2 $\mathrm{yr}^{-1}$. We estimate the sensitivity of our search over the entirety of Advanced LIGO's and Advanced Virgo's third observing run, and present the most stringent limits to date on the merger rate of binary black holes with at least one subsolar-mass component. We use the upper limits to constrain two fiducial scenarios that could produce subsolar-mass black holes: primordial black holes (PBH) and a model of dissipative dark matter. The PBH model uses recent prescriptions for the merger rate of PBH binaries that include a rate suppression factor to effectively account for PBH early binary disruptions. If the PBHs are monochromatically distributed, we can exclude a dark matter fraction in PBHs $f_\mathrm{PBH} \gtrsim 0.6$ (at 90% confidence) in the probed subsolar-mass range. However, if we allow for broad PBH mass distributions we are unable to rule out $f_\mathrm{PBH} = 1$. For the dissipative model, where the dark matter has chemistry that allows a small fraction to cool and collapse into black holes, we find an upper bound $f_{\mathrm{DBH}} < 10^{-5}$ on the fraction of atomic dark matter collapsed into black holes. △ Less

Submitted 26 January, 2024; v1 submitted 2 December, 2022; originally announced December 2022.

Comments: https://dcc.ligo.org/P2200139

Removal of Hot Saturns in Mass-Radius Plane by Runaway Mass Loss

Authors: Daniel P. Thorngren, Eve J. Lee, Eric D. Lopez

Abstract: The hot Saturn population exhibits a boundary in mass-radius space, such that no planets are observed at a density less than $\sim$0.1 g cm$^{-3}$. Yet, planet interior structure models can readily construct such objects as the natural result of radius inflation. Here, we investigate the role XUV-driven mass-loss plays in sculpting the density boundary by constructing interior structure models tha… ▽ More The hot Saturn population exhibits a boundary in mass-radius space, such that no planets are observed at a density less than $\sim$0.1 g cm$^{-3}$. Yet, planet interior structure models can readily construct such objects as the natural result of radius inflation. Here, we investigate the role XUV-driven mass-loss plays in sculpting the density boundary by constructing interior structure models that include radius inflation, photoevaporative mass loss and a simple prescription of Roche lobe overflow. We demonstrate that planets puffier than $\sim$0.1 g cm$^{-3}$ experience a runaway mass loss caused by adiabatic radius expansion as the gas layer is stripped away, providing a good explanation of the observed edge in mass-radius space. The process is also visible in the radius-period and mass-period spaces, though smaller, high-bulk-metallicity planets can still survive at short periods, preserving a partial record of the population distribution at formation. △ Less

Submitted 8 March, 2023; v1 submitted 21 November, 2022; originally announced November 2022.

Comments: 11 pages, 5 figures, Accepted to ApJ Letters. This revision includes the effect of orbital evolution resulting from Roche lobe overflow; the overall conclusions are the same

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

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

Search for gravitational-wave transients associated with magnetar bursts in Advanced LIGO and Advanced Virgo data from the third observing run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, A. Allocca, P. A. Altin , et al. (1645 additional authors not shown)

Abstract: Gravitational waves are expected to be produced from neutron star oscillations associated with magnetar giant flares and short bursts. We present the results of a search for short-duration (milliseconds to seconds) and long-duration ($\sim$ 100 s) transient gravitational waves from 13 magnetar short bursts observed during Advanced LIGO, Advanced Virgo and KAGRA's third observation run. These 13 bu… ▽ More Gravitational waves are expected to be produced from neutron star oscillations associated with magnetar giant flares and short bursts. We present the results of a search for short-duration (milliseconds to seconds) and long-duration ($\sim$ 100 s) transient gravitational waves from 13 magnetar short bursts observed during Advanced LIGO, Advanced Virgo and KAGRA's third observation run. These 13 bursts come from two magnetars, SGR 1935$+$2154 and Swift J1818.0$-$1607. We also include three other electromagnetic burst events detected by Fermi GBM which were identified as likely coming from one or more magnetars, but they have no association with a known magnetar. No magnetar giant flares were detected during the analysis period. We find no evidence of gravitational waves associated with any of these 16 bursts. We place upper bounds on the root-sum-square of the integrated gravitational-wave strain that reach $2.2 \times 10^{-23}$ $/\sqrt{\text{Hz}}$ at 100 Hz for the short-duration search and $8.7 \times 10^{-23}$ $/\sqrt{\text{Hz}}$ at $450$ Hz for the long-duration search, given a detection efficiency of 50%. For a ringdown signal at 1590 Hz targeted by the short-duration search the limit is set to $1.8 \times 10^{-22}$ $/\sqrt{\text{Hz}}$. Using the estimated distance to each magnetar, we derive upper bounds on the emitted gravitational-wave energy of $3.2 \times 10^{43}$ erg ($7.3 \times 10^{43}$ erg) for SGR 1935$+$2154 and $8.2 \times 10^{42}$ erg ($2.8 \times 10^{43}$ erg) for Swift J1818.0$-$1607, for the short-duration (long-duration) search. Assuming isotropic emission of electromagnetic radiation of the burst fluences, we constrain the ratio of gravitational-wave energy to electromagnetic energy for bursts from SGR 1935$+$2154 with available fluence information. The lowest of these ratios is $3 \times 10^3$. △ Less

Submitted 19 October, 2022; originally announced October 2022.

Comments: 30 pages with appendices, 5 figures, 10 tables

Report number: LIGO-P2100387

Model-based cross-correlation search for gravitational waves from the low-mass X-ray binary Scorpius X-1 in LIGO O3 data

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, C. Alléné, A. Allocca, P. A. Altin , et al. (1670 additional authors not shown)

Abstract: We present the results of a model-based search for continuous gravitational waves from the low-mass X-ray binary Scorpius X-1 using LIGO detector data from the third observing run of Advanced LIGO, Advanced Virgo and KAGRA. This is a semicoherent search which uses details of the signal model to coherently combine data separated by less than a specified coherence time, which can be adjusted to bala… ▽ More We present the results of a model-based search for continuous gravitational waves from the low-mass X-ray binary Scorpius X-1 using LIGO detector data from the third observing run of Advanced LIGO, Advanced Virgo and KAGRA. This is a semicoherent search which uses details of the signal model to coherently combine data separated by less than a specified coherence time, which can be adjusted to balance sensitivity with computing cost. The search covered a range of gravitational-wave frequencies from 25Hz to 1600Hz, as well as ranges in orbital speed, frequency and phase determined from observational constraints. No significant detection candidates were found, and upper limits were set as a function of frequency. The most stringent limits, between 100Hz and 200Hz, correspond to an amplitude h0 of about 1e-25 when marginalized isotropically over the unknown inclination angle of the neutron star's rotation axis, or less than 4e-26 assuming the optimal orientation. The sensitivity of this search is now probing amplitudes predicted by models of torque balance equilibrium. For the usual conservative model assuming accretion at the surface of the neutron star, our isotropically-marginalized upper limits are close to the predicted amplitude from about 70Hz to 100Hz; the limits assuming the neutron star spin is aligned with the most likely orbital angular momentum are below the conservative torque balance predictions from 40Hz to 200Hz. Assuming a broader range of accretion models, our direct limits on gravitational-wave amplitude delve into the relevant parameter space over a wide range of frequencies, to 500Hz or more. △ Less

Submitted 2 January, 2023; v1 submitted 6 September, 2022; originally announced September 2022.

Comments: 19 pages, Open Access Journal PDF

Report number: LIGO-P2100110-v13

Journal ref: The Astrophysical Journal Letters, 941, L30 (2022)