Showing 1–50 of 123 results for author: Tsai, S

Probing Cold-to-Temperate Exoplanetary Atmospheres: The Role of Water Condensation on Surface Identification with JWST

Authors: Ziyu Huang, Xinting Yu, Shang-Min Tsai, Julianne I. Moses, Kazumasa Ohno, Joshua Krissansen-Totton, Xi Zhang, Jonathan Fortney

Abstract: Understanding the surface temperature and interior structure of cold-to-temperate sub-Neptunes is critical for assessing their habitability, yet direct observations are challenging. In this study, we investigate the impact of water condensation on the atmospheric compositions of sub-Neptunes, focusing on the implications for JWST spectroscopic observations. By modeling the atmospheric photochemist… ▽ More Understanding the surface temperature and interior structure of cold-to-temperate sub-Neptunes is critical for assessing their habitability, yet direct observations are challenging. In this study, we investigate the impact of water condensation on the atmospheric compositions of sub-Neptunes, focusing on the implications for JWST spectroscopic observations. By modeling the atmospheric photochemistry of two canonical sub-Neptunes, K2-18 b and LHS 1140 b, both with and without water condensation and with and without thick atmospheres, we demonstrate that water condensation can significantly affect the predicted atmospheric compositions. This effect is driven by oxygen depletion from the condensation of water vapor and primarily manifests as an increase in the C/O ratio within the photochemically active regions of the atmosphere. This change in composition particularly affects planets with thin H$_2$-dominated atmospheres, leading to a transition in dominant nitrogen and carbon carriers from N$_2$ and oxygen-rich species like CO/CO$_2$ towards heavier hydrocarbons and nitriles. While our models do not fully account for the loss mechanisms of these higher-order species, such molecules can go on to form more refractory molecules or hazes. Planets with thin H2-rich atmospheres undergoing significant water condensation are thus likely to exhibit very hazy atmospheres. The relatively flat JWST spectra observed for LHS 1140 b could be consistent with such a scenario, suggesting a shallow surface with extensive water condensation or a high atmospheric C/O ratio. Conversely, the JWST observations of K2-18 b are better aligned with a volatile-rich mini-Neptune with a thick atmosphere. △ Less

Submitted 12 July, 2024; originally announced July 2024.

Comments: 21 pages, 7 figures

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)

Searching Accretion-Enhanced Dark Matter Annihilation Signals in the Galactic Centre

Authors: Mei-Wen Yang, Zhi-Qi Guo, Xiao-Yi Luo, Zhao-Qiang Shen, Zi-Qing Xia, Chih-Ting Lu, Yue-Lin Sming Tsai, Yi-Zhong Fan

Abstract: This study reanalyzes the detection prospects of dark matter (DM) annihilation signals in the Galactic Center, focusing on velocity-dependent dynamics within a spike density near the supermassive black hole (Sgr~A$^{\star}$). We investigate three annihilation processes -- $p$-wave, resonance, and forbidden annihilation -- under semi-relativistic velocities, leveraging gamma-ray data from Fermi and… ▽ More This study reanalyzes the detection prospects of dark matter (DM) annihilation signals in the Galactic Center, focusing on velocity-dependent dynamics within a spike density near the supermassive black hole (Sgr~A$^{\star}$). We investigate three annihilation processes -- $p$-wave, resonance, and forbidden annihilation -- under semi-relativistic velocities, leveraging gamma-ray data from Fermi and DAMPE telescopes. Our analysis integrates a fermionic DM model with an electroweak axion-like particle (ALP) portal, exploring annihilation into two or four photons. Employing a comprehensive six-dimensional integration, we precisely calculate DM-induced gamma-ray fluxes near Sgr~A$^{\star}$, incorporating velocity and positional dependencies in the annihilation cross-section and photon yield spectra. Our findings highlight scenarios of resonance and forbidden annihilation, where the larger ALP-DM-DM coupling constant $C_{aχχ}$ can affect spike density, potentially yielding detectable gamma-ray line spectra within Fermi and DAMPE energy resolution. We set upper limits for $C_{aχχ}$ across these scenarios, offering insights into the detectability and spectral characteristics of DM annihilation signals from the Galactic Center. △ Less

Submitted 9 July, 2024; originally announced July 2024.

A new lever on exoplanetary B fields: measuring heavy ion velocities

Authors: Arjun B. Savel, Hayley Beltz, Thaddeus D. Komacek, Shang-Min Tsai, Eliza M. -R. Kempton

Abstract: Magnetic fields connect an array of planetary processes, from atmospheric escape to interior convection. Despite their importance, exoplanet magnetic fields are largely unconstrained by both theory and observation. In this Letter, we propose a novel method for constraining the B field strength of hot gas giants: comparing the velocities of heavy ions and neutral gas with high-resolution spectrosco… ▽ More Magnetic fields connect an array of planetary processes, from atmospheric escape to interior convection. Despite their importance, exoplanet magnetic fields are largely unconstrained by both theory and observation. In this Letter, we propose a novel method for constraining the B field strength of hot gas giants: comparing the velocities of heavy ions and neutral gas with high-resolution spectroscopy. The core concept of this method is that ions are directly deflected by magnetic fields. While neutrals are also affected by B fields via friction with field-accelerated ions, ionic gas should be more strongly coupled to the underlying magnetic field than bulk neutral flow. Hence, measuring the difference between the two velocities yields rough constraints on the B field, provided an estimate of the stellar UV flux is known. We demonstrate that heavy ions are particularly well suited for this technique, because they are less likely to be entrained in complex hydrodynamic outflows than their lighter counterparts. We perform a proof-of-concept calculation with Ba II, an ion whose velocity has been repeatedly measured at high confidence with high-resolution spectroscopy. Our work shows that a 10G magnetic field would produce ~ km/s ion--neutral velocity differences at a microbar, whereas a 50G magnetic field would produce ~20km/s velocity difference. With new leverage on magnetic fields, we will be able to investigate magnetic field generation in the extreme edge cases of hot gas giants, with wide-ranging consequences for planetary interior structure, dynamo theory, and habitability. △ Less

Submitted 24 June, 2024; v1 submitted 18 June, 2024; originally announced June 2024.

Comments: Accepted at ApJ Letters. 12 pages; 4 figures; 2 tables

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)

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

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

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

Submitted 5 June, 2024; originally announced June 2024.

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

Debris Disks can Contaminate Mid-Infrared Exoplanet Spectra: Evidence for a Circumstellar Debris Disk around Exoplanet Host WASP-39

Authors: Laura Flagg, Alycia J. Weinberger, Taylor J. Bell, Luis Welbanks, Giuseppe Morello, Diana Powell, Jacob L. Bean, Jasmina Blecic, Nicolas Crouzet, Peter Gao, Julie Inglis, James Kirk, Mercedes Lopez-Morales, Karan Molaverdikhani, Nikolay Nikolov, Apurva V. Oza, Benjamin V. Rackham, Seth Redfield, Shang-Min Tsai, Ray Jayawardhana, Laura Kreidberg, Matthew C. Nixon, Kevin B. Stevenson, Jake D. Turner

Abstract: The signal from a transiting planet can be diluted by astrophysical contamination. In the case of circumstellar debris disks, this contamination could start in the mid-infrared and vary as a function of wavelength, which would then change the observed transmission spectrum for any planet in the system. The MIRI/LRS WASP-39b transmission spectrum shows an unexplained dip starting at $\sim$10 $μ$m t… ▽ More The signal from a transiting planet can be diluted by astrophysical contamination. In the case of circumstellar debris disks, this contamination could start in the mid-infrared and vary as a function of wavelength, which would then change the observed transmission spectrum for any planet in the system. The MIRI/LRS WASP-39b transmission spectrum shows an unexplained dip starting at $\sim$10 $μ$m that could be caused by astrophysical contamination. The spectral energy distribution displays excess flux at similar levels to that which are needed to create the dip in the transmission spectrum. In this article, we show that this dip is consistent with the presence of a bright circumstellar debris disk, at a distance of $>$2 au. We discuss how a circumstellar debris disk like that could affect the atmosphere of WASP-39b. We also show that even faint debris disks can be a source of contamination in MIRI exoplanet spectra. △ Less

Submitted 4 June, 2024; originally announced June 2024.

Comments: accepted to ApJL

The relic density and temperature evolution of light dark sector

Authors: Xin-Chen Duan, Raymundo Ramos, Yue-Lin Sming Tsai

Abstract: We have developed a set of four fully coupled Boltzmann equations to precisely determine the relic density and temperature of dark matter by including three distinct sectors: dark matter, light scalar, and standard model sectors. The intricacies of heat transfer between DM and the SM sector through a light scalar particle are explored, inspired by stringent experimental constraints on the scalar-H… ▽ More We have developed a set of four fully coupled Boltzmann equations to precisely determine the relic density and temperature of dark matter by including three distinct sectors: dark matter, light scalar, and standard model sectors. The intricacies of heat transfer between DM and the SM sector through a light scalar particle are explored, inspired by stringent experimental constraints on the scalar-Higgs mixing angle and the DM-scalar coupling. Three distinct sectors emerge prior to DM freeze-out, requiring fully coupled Boltzmann equations to accurately compute relic density. Investigation of forbidden, resonance, and secluded DM scenarios demonstrates significant deviations between established methods and the novel approach with fully coupled Boltzmann equations. Despite increased computational demands, this emphasizes the need for improved precision in relic density calculations, underlining the importance of incorporating these equations in comprehensive analyses. △ Less

Submitted 18 April, 2024; originally announced April 2024.

Comments: 38 pages, 10 figures

Biogenic sulfur gases as biosignatures on temperate sub-Neptune waterworlds

Authors: Shang-Min Tsai, Hamish Innes, Nicholas F. Wogan, Edward W. Schwieterman

Abstract: Theoretical predictions and observational data indicate a class of sub-Neptune exoplanets may have water-rich interiors covered by hydrogen-dominated atmospheres. Provided suitable climate conditions, such planets could host surface liquid oceans. Motivated by recent JWST observations of K2-18 b, we self-consistently model the photochemistry and potential detectability of biogenic sulfur gases in… ▽ More Theoretical predictions and observational data indicate a class of sub-Neptune exoplanets may have water-rich interiors covered by hydrogen-dominated atmospheres. Provided suitable climate conditions, such planets could host surface liquid oceans. Motivated by recent JWST observations of K2-18 b, we self-consistently model the photochemistry and potential detectability of biogenic sulfur gases in the atmospheres of temperate sub-Neptune waterworlds for the first time. On Earth today, organic sulfur compounds produced by marine biota are rapidly destroyed by photochemical processes before they can accumulate to significant levels. Domagal-Goldman et al. (2011) suggest that detectable biogenic sulfur signatures could emerge in Archean-like atmospheres with higher biological production or low UV flux. In this study, we explore biogenic sulfur across a wide range of biological fluxes and stellar UV environments. Critically, the main photochemical sinks are absent on the nightside of tidally locked planets. To address this, we further perform experiments with a 3D GCM and a 2D photochemical model (VULCAN 2D (Tsai et al. 2024)) to simulate the global distribution of biogenic gases to investigate their terminator concentrations as seen via transmission spectroscopy. Our models indicate that biogenic sulfur gases can rise to potentially detectable levels on hydrogen-rich waterworlds, but only for enhanced global biosulfur flux ($\gtrsim$20 times modern Earth's flux). We find that it is challenging to identify DMS at 3.4 $μm$ where it strongly overlaps with CH$_4$, whereas it is more plausible to detect DMS and companion byproducts, ethylene (C$_2$H$_4$) and ethane (C$_2$H$_6$), in the mid-infrared between 9 and 13 $μm$. △ Less

Submitted 26 March, 2024; v1 submitted 21 March, 2024; originally announced March 2024.

Comments: 9 pages, 4 figures, accepted for publication in ApJL

Light Thermal Dark Matter Beyond $p$-Wave Annihilation in Minimal Higgs Portal Model

Authors: Yu-Tong Chen, Shigeki Matsumoto, Tian-Peng Tang, Yue-Lin Sming Tsai, Lei Wu

Abstract: This study explores a minimal renormalizable dark matter (DM) model, incorporating a sub-GeV Majorana DM and a singlet scalar particle $φ$. Using scalar and pseudo-scalar interactions (couplings $c_s$ and $c_p$), we investigate implications for DM detection, considering $s$-wave, $p$-wave, and combined ($s$+$p$ wave) contributions in DM annihilation cross-section, as well as loop-correction contri… ▽ More This study explores a minimal renormalizable dark matter (DM) model, incorporating a sub-GeV Majorana DM and a singlet scalar particle $φ$. Using scalar and pseudo-scalar interactions (couplings $c_s$ and $c_p$), we investigate implications for DM detection, considering $s$-wave, $p$-wave, and combined ($s$+$p$ wave) contributions in DM annihilation cross-section, as well as loop-correction contributions to DM-nucleon elastic scattering. Identifying a broad parameter space ($10 \,\rm{MeV} < m_χ\lesssim m_φ$) within the $2σ$ allowed region, we explore scenarios ($\left|c_s\right|\gg \left|c_p\right|$, $\left|c_s\right|\ll \left|c_p\right|$, and $\left|c_s\right|\approx \left|c_p\right|$). We find that (i) a non-zero pseudo-scalar coupling alleviates direct detection constraints as a comparison with the previous pure scalar coupling case; (ii) CMB observations set stringent limits on pseudo-scalar interaction dominant cases, making $s$-wave annihilation viable only for $m_χ>1\,\rm{GeV}$; (iii) the preferred $φ$-resonance region can be tested in the future indirect detection experiments, such as e-ASTROGAM. △ Less

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

Comments: 35 pages, 4 figures

Modeling the Progenitor Stars of Observed IIP Supernovae

Authors: Kai-An You, Ke-Jung Chen, Yen-Chen Pan, Sung-Han Tsai, Po-Sheng Ou

Abstract: Type IIP supernovae (SNe IIP) are thought to originate from the explosion of massive stars > 10 Msun. Their luminosity is primarily powered by the explosion energy and the radioactive decay energy of 56Co, with the photosphere location regulated by hydrogen recombination. However, the physical connections between SNe IIP and their progenitor stars remain unclear. This paper presents a comprehensiv… ▽ More Type IIP supernovae (SNe IIP) are thought to originate from the explosion of massive stars > 10 Msun. Their luminosity is primarily powered by the explosion energy and the radioactive decay energy of 56Co, with the photosphere location regulated by hydrogen recombination. However, the physical connections between SNe IIP and their progenitor stars remain unclear. This paper presents a comprehensive study of SNe IIP and their progenitor stars by using the one-dimensional stellar evolution code, MESA. Our model grids consider the effects of stellar metallicity, mass, and rotation in the evolution of massive stars, as well as explosion energy and 56Ni production in modeling supernovae. To elucidate the observed SNe IIP and their origins, we compare their light curves (LCs) with our models. Furthermore, we investigate the impact of stellar parameters on LCs by considering stellar mass metallicity, rotation, explosion energy, and 56Ni production. We find that more massive stars exhibit longer plateaus due to increased photon diffusion time caused by massive ejecta. Higher metallicity leads to increased opacity and mass loss of progenitor stars. Rapid rotation affects internal stellar structures, enhancing convective mixing and mass loss, potentially affecting the plateau's brightness and duration. Higher explosion energy results in brighter but shorter plateaus due to faster-moving ejecta. 56Ni mass affects late-time luminosity and plateau duration, with larger masses leading to slower declines. △ Less

Submitted 27 May, 2024; v1 submitted 29 February, 2024; originally announced February 2024.

Comments: 22 pages, 7 figures. Accepted for publication in ApJ

Weak Lensing Constraints on Dark Matter-Baryon Interactions with $N$-Body Simulations and Machine Learning

Authors: Chi Zhang, Lei Zu, Hou-Zun Chen, Yue-Lin Sming Tsai, Yi-Zhong Fan

Abstract: We investigate the elastic scattering cross section between dark matter and protons using the DES Year 3 weak lensing data. This scattering induces a dark acoustic oscillation structure in the matter power spectra. To address non-linear effects at low redshift, we utilize principal component analysis alongside a limited set of $N$-body simulations, improving the reliability of our matter power spe… ▽ More We investigate the elastic scattering cross section between dark matter and protons using the DES Year 3 weak lensing data. This scattering induces a dark acoustic oscillation structure in the matter power spectra. To address non-linear effects at low redshift, we utilize principal component analysis alongside a limited set of $N$-body simulations, improving the reliability of our matter power spectrum prediction. We further perform a robust Markov Chain Monte Carlo analysis to derive the upper bounds on the DM-proton elastic scattering cross-section, assuming different velocity dependencies. Our results, presented as the first Frequentist upper limits, are compared with the ones obtained by Bayesian approach. Compared with the upper limits derived from the Planck cosmic microwave background data, our findings from DES Year 3 data exhibit improvements of up to a factor of five. In addition, we forecast the future sensitivities of the China Space Station Telescope, the upcoming capabilities of this telescope could improve the current limits by approximately one order of magnitude. △ Less

Submitted 29 February, 2024; originally announced February 2024.

Comments: 10 pages, 5 figures

JWST observations of K2-18b can be explained by a gas-rich mini-Neptune with no habitable surface

Authors: Nicholas F. Wogan, Natasha E. Batalha, Kevin Zahnle, Joshua Krissansen-Totton, Shang-Min Tsai, Renyu Hu

Abstract: JWST recently measured the transmission spectrum of K2-18b, a habitable-zone sub-Neptune exoplanet, detecting CH$_4$ and CO$_2$ in its atmosphere. The discovery paper argued the data are best explained by a habitable "Hycean" world, consisting of a relatively thin H$_2$-dominated atmosphere overlying a liquid water ocean. Here, we use photochemical and climate models to simulate K2-18b as both a H… ▽ More JWST recently measured the transmission spectrum of K2-18b, a habitable-zone sub-Neptune exoplanet, detecting CH$_4$ and CO$_2$ in its atmosphere. The discovery paper argued the data are best explained by a habitable "Hycean" world, consisting of a relatively thin H$_2$-dominated atmosphere overlying a liquid water ocean. Here, we use photochemical and climate models to simulate K2-18b as both a Hycean planet and a gas-rich mini-Neptune with no defined surface. We find that a lifeless Hycean world is hard to reconcile with the JWST observations because photochemistry only supports $< 1$ part-per-million CH$_4$ in such an atmosphere while the data suggest about $\sim 1\%$ of the gas is present. Sustaining %-level CH$_4$ on a Hycean K2-18b may require the presence of a methane-producing biosphere, similar to microbial life on Earth $\sim 3$ billion years ago. On the other hand, we predict that a gas-rich mini-Neptune with $100 \times$ solar metallicity should have 4% CH$_4$ and nearly 0.1% CO$_2$, which are compatible with the JWST data. The CH$_4$ and CO$_2$ are produced thermochemically in the deep atmosphere and mixed upward to the low pressures sensitive to transmission spectroscopy. The model predicts H$_2$O, NH$_3$ and CO abundances broadly consistent with the non-detections. Given the additional obstacles to maintaining a stable temperate climate on Hycean worlds due to H$_2$ escape and potential supercriticality at depth, we favor the mini-Neptune interpretation because of its relative simplicity and because it does not need a biosphere or other unknown source of methane to explain the data. △ Less

Submitted 3 February, 2024; v1 submitted 19 January, 2024; originally announced January 2024.

Comments: Accepted for publication at ApJL

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

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

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

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

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

Comments: 19 pages, 20 figures, published in ApJ

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

Inelastic Scattering of Dark Matter with Heavy Cosmic Rays

Authors: Keyu Lu, Yue-Lin Sming Tsai, Qiang Yuan, Le Zhang

Abstract: We investigate the impact of inelastic collisions between dark matter (DM) and heavy cosmic ray (CR) nuclei on CR propagation. We approximate the fragmentation cross-sections for DM-CR collisions using collider-measured proton-nuclei scattering cross-sections, allowing us to assess how these collisions affect the spectra of CR Boron and Carbon. We derive new CR spectra from DM-CR collisions by inc… ▽ More We investigate the impact of inelastic collisions between dark matter (DM) and heavy cosmic ray (CR) nuclei on CR propagation. We approximate the fragmentation cross-sections for DM-CR collisions using collider-measured proton-nuclei scattering cross-sections, allowing us to assess how these collisions affect the spectra of CR Boron and Carbon. We derive new CR spectra from DM-CR collisions by incorporating their cross-sections into the source terms and solving the diffusion equation for the complete network of reactions involved in generating secondary species. In a specific example with a coupling strength of $b_χ=0.1$ and a DM mass of $m_χ=0.1$ GeV, considering a simplified scenario where DM interacts exclusively with Oxygen, a notable modification in the Boron-to-Carbon spectrum due to the DM-CR interaction is observed. Particularly, the peak within the spectrum, spanning from $0.1$ GeV to $10$ GeV, experiences an enhancement of approximately 1.5 times. However, in a more realistic scenario where DM particles interact with all CRs, this peak can be amplified to twice its original value.Utilizing the latest data from AMS-02 and DAMPE on the Boron-to-Carbon ratio, we estimate a 95\% upper limit for the effective inelastic cross-section of DM-proton as a function of DM mass. Our findings reveal that at $m_χ\simeq 2$ MeV, the effective inelastic cross-section between DM and protons must be less than $\mathcal{O}(10^{-32})~{\rm cm}^2$. △ Less

Submitted 7 June, 2024; v1 submitted 19 October, 2023; originally announced October 2023.

Comments: 19 pages, 8 figures

Journal ref: Research in Astronomy and Astrophysics 24, 065007 (2024)

Primordial magnetic field as a common solution of nanohertz gravitational waves and the Hubble tension

Authors: Yao-Yu Li, Chi Zhang, Ziwei Wang, Ming-Yang Cui, Yue-Lin Sming Tsai, Qiang Yuan, Yi-Zhong Fan

Abstract: The origin of interstellar and intergalactic magnetic fields remains largely unknown. One possibility is that they are related to the primordial magnetic fields (PMFs) produced by, for instance, the phase transitions of the early Universe. In this paper, we show that the PMF-induced turbulence generated at around the QCD phase transition epoch--the characteristic magnetic field strength… ▽ More The origin of interstellar and intergalactic magnetic fields remains largely unknown. One possibility is that they are related to the primordial magnetic fields (PMFs) produced by, for instance, the phase transitions of the early Universe. In this paper, we show that the PMF-induced turbulence generated at around the QCD phase transition epoch--the characteristic magnetic field strength $B_{\rm ch}^* \sim \mathcal{O}(1)~\rm{μG}$ and coherent length scale $\ell_{\rm ch}^* \sim \mathcal{O}(1)~\rm{pc}$--can naturally accommodate nanohertz gravitational waves reported by pulsar timing array (PTA) collaborations. Moreover, the evolution of the PMFs to the recombination era with the form of $B_{\rm ch}\sim \ell_{\rm ch}^{-α}$ can induce baryon density inhomogeneities, alter the recombination history, and alleviate the tension of the Hubble parameter $H_0$ and the matter clumpiness parameter $S_8$ between early- and late-time measurements for $0.88\leq α\leq 1.17$ (approximate 95\% credible region based on three PTA likelihoods). The further evolved PMFs may account for the $\sim {\cal O}(10^{-16})$ Gauss extragalactic magnetic field inferred with GRB 221009A. △ Less

Submitted 4 March, 2024; v1 submitted 29 June, 2023; originally announced June 2023.

Comments: 7 pages, 3 figures, 1 table

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

Mirror QCD phase transition as the origin of the nanohertz Stochastic Gravitational-Wave Background

Authors: Lei Zu, Chi Zhang, Yao-Yu Li, Yu-Chao Gu, Yue-Lin Sming Tsai, Yi-Zhong Fan

Abstract: Several Pulsar Timing Array (PTA) collaborations have recently provided strong evidence for a nHz Stochastic Gravitational-Wave Background (SGWB). Here we investigate the implications of a first-order phase transition occurring within the early universe's dark quantum chromodynamics (dQCD) epoch, specifically within the framework of the mirror twin Higgs dark sector model. Our analysis indicates a… ▽ More Several Pulsar Timing Array (PTA) collaborations have recently provided strong evidence for a nHz Stochastic Gravitational-Wave Background (SGWB). Here we investigate the implications of a first-order phase transition occurring within the early universe's dark quantum chromodynamics (dQCD) epoch, specifically within the framework of the mirror twin Higgs dark sector model. Our analysis indicates a distinguishable SGWB signal originating from this phase transition, which can explain the measurements obtained by PTAs. Remarkably, a significant portion of the parameter space for the SGWB signal also effectively resolves the existing tensions in both the $H_0$ and $S_8$ measurements in Cosmology. This intriguing correlation suggests a possible common origin of these three phenomena for $0.2 < ΔN_{\rm eff} < 0.5$, where the mirror dark matter component constitutes about $30\%$ of the total dark matter abundance. Next generation CMB experiment such as CMB-S4 is able to test this parameter region. △ Less

Submitted 2 July, 2024; v1 submitted 29 June, 2023; originally announced June 2023.

Comments: 7 pages, 4 figures

Journal ref: Sci.Bull. 69 (2024) 741-746

R-process beta-decay neutrino flux from binary neutron star mergers and collapsars

Authors: Yu An, Meng-Ru Wu, Gang Guo, Yue-Lin Sming Tsai, Shih-Jie Huang, Yi-Zhong Fan

Abstract: This study investigates the antineutrinos production by $β$-decay of $r$-process nuclei in two astrophysical sites that are capable of producing gamma-ray bursts (GRBs): binary neutron star mergers (BNSMs) and collapsars, which are promising sites for heavy element nucleosynthesis. We employ a simplified method to compute the $β$-decay $\barν_e$ energy spectrum and consider a number of different r… ▽ More This study investigates the antineutrinos production by $β$-decay of $r$-process nuclei in two astrophysical sites that are capable of producing gamma-ray bursts (GRBs): binary neutron star mergers (BNSMs) and collapsars, which are promising sites for heavy element nucleosynthesis. We employ a simplified method to compute the $β$-decay $\barν_e$ energy spectrum and consider a number of different representative thermodynamic trajectories for $r$-process simulations, each with four sets of $Y_e$ distribution. The time evolution of the $\barν_e$ spectrum is derived for both the dynamical ejecta and the disk wind for BNSMs and collapsar outflow, based on approximated mass outflow rates. Our results show that the $\barν_e$ has an average energy of approximately 3 to 9~MeV, with a high energy tail of up to 20 MeV. The $\barν_e$ flux evolution is primarily determined by the outflow duration, and can thus remain large for $\mathcal{O}(10)$~s and $\mathcal{O}(100)$~s for BNSMs and collapsars, respectively. For a single merger or collapsar at 40~Mpc, the $\barν_e$ flux is $\mathcal{O}(10-100)$~cm$^{-2}$~s$^{-1}$, indicating a possible detection horizon up to $0.1-1$~Mpc for Hyper-Kamiokande. We also estimate their contributions to the diffuse $\barν_e$ background, and find that both sources should only contribute subdominantly to the diffuse background when compared to that expected from core-collapse supernovae. △ Less

Submitted 15 December, 2023; v1 submitted 13 June, 2023; originally announced June 2023.

Comments: 14 pages, 8 figures

Journal ref: Phys. Rev. D 108, 123038 (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)

The Evolution of Population III and Extremely Metal-Poor Binary Stars

Authors: Sung-Han Tsai, Ke-Jung Chen, Daniel Whalen, Po-Sheng Ou, Tyrone E. Woods

Abstract: Numerical simulations have now shown that Population III (Pop III) stars can form in binaries and small clusters and that these stars can be in close proximity to each other. If so, they could be subject to binary interactions such as mass exchange that could profoundly alter their evolution, ionizing UV and Lyman-Werner (LW) photon emission and explosion yields, with important consequences for ea… ▽ More Numerical simulations have now shown that Population III (Pop III) stars can form in binaries and small clusters and that these stars can be in close proximity to each other. If so, they could be subject to binary interactions such as mass exchange that could profoundly alter their evolution, ionizing UV and Lyman-Werner (LW) photon emission and explosion yields, with important consequences for early cosmological reionization and chemical enrichment. Here we investigate the evolution of Pop III and extremely metal-poor binary stars with the MESA code. We find that interactions ranging from stable mass transfer to common envelope evolution can occur in these binaries for a wide range of mass ratios and initial separations. Mass transfer can nearly double UV photon yields in some of these binaries with respect to their individual stars by extending the life of the companion star, which in turn can enhance early cosmological reionization but also suppress the formation of later generations of primordial stars. Binary interactions can also have large effects on the nucleosynthetic yields of the stars by promoting or removing them into or out of mass ranges for specific SN types. We provide fits to total photon yields for the binaries in our study for use in cosmological simulations. △ Less

Submitted 29 May, 2023; v1 submitted 20 April, 2023; originally announced April 2023.

Comments: 53 pages, 17 tables. Accepted for publication in ApJ

Exploring Mirror Twin Higgs Cosmology with Present and Future Weak Lensing Surveys

Authors: Lei Zu, Chi Zhang, Hou-Zun Chen, Wei Wang, Yue-Lin Sming Tsai, Yuhsin Tsai, Wentao Luo, Yi-Zhong Fan

Abstract: We explore the potential of precision cosmological data to study non-minimal dark sectors by updating the cosmological constraint on the mirror twin Higgs model (MTH). The MTH model addresses the Higgs little hierarchy problem by introducing dark sector particles. In this work, we perform a Bayesian global analysis that includes the latest cosmic shear measurement from the DES three-year survey an… ▽ More We explore the potential of precision cosmological data to study non-minimal dark sectors by updating the cosmological constraint on the mirror twin Higgs model (MTH). The MTH model addresses the Higgs little hierarchy problem by introducing dark sector particles. In this work, we perform a Bayesian global analysis that includes the latest cosmic shear measurement from the DES three-year survey and the Planck CMB and BAO data. In the early Universe, the mirror baryon and mirror radiation behave as dark matter and dark radiation, and their presence modifies the Universe's expansion history. Additionally, the scattering between mirror baryon and photon generates the dark acoustic oscillation process, suppressing the matter power spectrum from the cosmic shear measurement. We demonstrate how current data constrain these corrections to the $Λ$CDM cosmology and find that for a viable solution to the little hierarchy problem, the proportion of MTH dark matter cannot exceed about $30\%$ of the total dark matter density, unless the temperature of twin photon is less than $30\%$ of that of the standard model photon. While the MTH model is presently not a superior solution to the observed $H_0$ tension compared to the $Λ$CDM+$ΔN_{\rm eff}$ model, we demonstrate that it has the potential to alleviate both the $H_0$ and $S_8$ tensions, especially if the $S_8$ tension persists in the future and approaches the result reported by the Planck SZ (2013) analysis. In this case, the MTH model can relax the tensions while satisfying the DES power spectrum constraint up to $k \lesssim 10~h\rm {Mpc}^{-1}$. If the MTH model is indeed accountable for the $S_8$ and $H_0$ tensions, we show that the future China Space Station Telescope (CSST) can determine the twin baryon abundance with a $10\%$ level precision. △ Less

Submitted 5 October, 2023; v1 submitted 13 April, 2023; originally announced April 2023.

Comments: 32 pages, 12 figures, 4 tables. Published in JCAP

Journal ref: JCAP08(2023)023

The Runaway Greenhouse Effect on Hycean Worlds

Authors: Hamish Innes, Shang-Min Tsai, Raymond T. Pierrehumbert

Abstract: Hycean worlds are a proposed subset of sub-Neptune exoplanets with substantial water inventories, liquid surface oceans and extended hydrogen-dominated atmospheres that could be favourable for habitability. In this work, we aim to quantitatively define the inner edge of the Hycean habitable zone using a 1D radiative-convective model. As a limiting case, we model a dry hydrogen-helium envelope abov… ▽ More Hycean worlds are a proposed subset of sub-Neptune exoplanets with substantial water inventories, liquid surface oceans and extended hydrogen-dominated atmospheres that could be favourable for habitability. In this work, we aim to quantitatively define the inner edge of the Hycean habitable zone using a 1D radiative-convective model. As a limiting case, we model a dry hydrogen-helium envelope above a surface ocean. We find that 10 to 20 bars of atmosphere produces enough greenhouse effect to drive a liquid surface ocean supercritical when forced with current Earth-like instellation. Introducing water vapour into the atmosphere, we show the runaway greenhouse instellation limit is greatly reduced due to the presence of superadiabatic layers where convection is inhibited. This moves the inner edge of the habitable zone from $\approx$ 1 AU for a G-star to 1.6 AU (3.85 AU) for a Hycean world with a H$_2$-He inventory of 1 bar (10 bar). For an M-star, the inner edge is equivalently moved from 0.17 AU to 0.28 AU (0.54 AU). Our results suggest that most of the current Hycean world observational targets are not likely to sustain a liquid water ocean. We present an analytical framework for interpreting our results, finding that the maximum possible OLR scales approximately inversely with the dry mass inventory of the atmosphere. We discuss the possible limitations of our 1D modelling and recommend the use of 3D convection-resolving models to explore the robustness of superadiabatic layers. △ Less

Submitted 5 April, 2023; originally announced April 2023.

Comments: Submitted to ApJ, comments welcome

High atmospheric metal enrichment for a Saturn-mass planet

Authors: Jacob L. Bean, Qiao Xue, Prune C. August, Jonathan Lunine, Michael Zhang, Daniel Thorngren, Shang-Min Tsai, Keivan G. Stassun, Everett Schlawin, Eva-Maria Ahrer, Jegug Ih, Megan Mansfield

Abstract: Atmospheric metal enrichment (i.e., elements heavier than helium, also called "metallicity") is a key diagnostic of the formation of giant planets. The giant planets of the solar system exhibit an inverse relationship between mass and both their bulk metallicities and atmospheric metallicities. Extrasolar giant planets also display an inverse relationship between mass and bulk metallicity. However… ▽ More Atmospheric metal enrichment (i.e., elements heavier than helium, also called "metallicity") is a key diagnostic of the formation of giant planets. The giant planets of the solar system exhibit an inverse relationship between mass and both their bulk metallicities and atmospheric metallicities. Extrasolar giant planets also display an inverse relationship between mass and bulk metallicity. However, there is significant scatter in the relationship and it is not known how atmospheric metallicity correlates with either planet mass or bulk metallicity. Here we show that the Saturn-mass exoplanet HD 149026b has an atmospheric metallicity 59 - 276 times solar (at 1 $σ$), which is greater than Saturn's atmospheric metallicity of ~7.5 times solar at >4 $σ$ confidence. This result is based on modeling CO$_2$ and H$_2$O absorption features in the thermal emission spectrum of the planet measured by JWST. HD 149026b is the most metal-rich giant planet known, with an estimated bulk heavy element abundance of 66 $\pm$ 2% by mass. We find that the atmospheric metallicities of both HD 149026b and the solar system giant planets are more correlated with bulk metallicity than planet mass. △ Less

Submitted 24 March, 2023; originally announced March 2023.

Comments: Published online in Nature on March 27, 2023; a JWST thermal emission spectrum of a planet that definitely has an atmosphere ;)

Exploring dark matter spike distribution around the Galactic centre with stellar orbits

Authors: Zhao-Qiang Shen, Guan-Wen Yuan, Cheng-Zi Jiang, Yue-Lin Sming Tsai, Qiang Yuan, Yi-Zhong Fan

Abstract: Precise measurements of the stellar orbits around Sagittarius A* have established the existence of a supermassive black hole (SMBH) at the Galactic centre (GC). Due to the interplay between the SMBH and dark matter (DM), the DM density profile in the innermost region of the Galaxy, which is crucial for the DM indirect detection, is still an open question. Among the most popular models in the liter… ▽ More Precise measurements of the stellar orbits around Sagittarius A* have established the existence of a supermassive black hole (SMBH) at the Galactic centre (GC). Due to the interplay between the SMBH and dark matter (DM), the DM density profile in the innermost region of the Galaxy, which is crucial for the DM indirect detection, is still an open question. Among the most popular models in the literature, the theoretical spike profile proposed by Gondolo and Silk (1999; GS hereafter) is well adopted. In this work, we investigate the DM spike profile using updated data from the Keck and VLT telescopes considering that the presence of such an extended mass component may affect the orbits of the S-stars in the Galactic center. We examine the radius and slope of the generalized NFW spike profile, analyze the Einasto spike, and discuss the influence of DM annihilation on the results. Our findings indicate that an initial slope of $γ\gtrsim 0.92$ for the generalized NFW spike profile is ruled out at a 95% confidence level. Additionally, the spike radius $R_{\rm sp}$ larger than 21.5 pc is rejected at 95% probability for the Einasto spike with $α=0.17$, which also contradicts the GS spike model. The constraints with the VLT/GRAVITY upper limits are also projected. Although the GS NFW spike is well constrained by the Keck and VLT observation of S2, an NFW spike with a weak annihilation cusp may still be viable, as long as the DM annihilation cross section satisfies $\left< σv \right> \gtrsim 7.7\times 10^{-27}~{\rm cm^3\,s^{-1}} (m_{\rm DM}/100~{\rm GeV})$ at 95% level. △ Less

Submitted 24 October, 2023; v1 submitted 16 March, 2023; originally announced March 2023.

Comments: 12 pages, 10 figures and 1 table. Accepted for publication in MNRAS

Journal ref: MNRAS 527, 3196 (2024)

Multi-frequency test of dark matter annihilation into long-lived particles in Sirius

Authors: Yu-Xuan Chen, Lei Zu, Zi-Qing Xia, Yue-Lin Sming Tsai, Yi-Zhong Fan

Abstract: New long-lived particles produced at the colliders may escape from conventional particle detectors. Using satellites or ground telescopes, we can detect the photons generated from the annihilation of the star-captured dark matter into a pair of long-lived particles. When the propagation length of these long-lived particles surpasses the interplanetary distance between the Sun and Jupiter, it becom… ▽ More New long-lived particles produced at the colliders may escape from conventional particle detectors. Using satellites or ground telescopes, we can detect the photons generated from the annihilation of the star-captured dark matter into a pair of long-lived particles. When the propagation length of these long-lived particles surpasses the interplanetary distance between the Sun and Jupiter, it becomes unfeasible to detect such dark matter signals originating from the Sun or Jupiter on Earth. Our analysis of the dark matter-induced photons produced by prompt radiation, inverse Compton scattering, and synchrotron radiation mechanisms reveals that a decay length of about $10^{-3}$ pc for long-lived particles is required for maximum detectability. We investigate the parameters that allow the long-lived particle's lifetime to be consistent with Big Bang nucleosynthesis while also allowing it to escape the confines of our solar system. The Sirius system is proposed as a promising target for the indirect detection of such long-lived particles. Utilizing the prompt, inverse Compton scattering, and synchrotron radiation, upper limits on the dark matter-proton spin-independent and spin-dependent cross section are estimated with the Fermi-LAT null-signal observation and the capabilities of the upcoming Square Kilometre Array radio telescope. △ Less

Submitted 28 July, 2023; v1 submitted 20 February, 2023; originally announced February 2023.

Comments: 30 pages, 6 figures, JCAP accepted

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

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 NIRCam

Authors: Eva-Maria Ahrer, Kevin B. Stevenson, Megan Mansfield, Sarah E. Moran, Jonathan Brande, Giuseppe Morello, Catriona A. Murray, Nikolay K. Nikolov, Dominique J. M. Petit dit de la Roche, Everett Schlawin, Peter J. Wheatley, Sebastian Zieba, Natasha E. Batalha, Mario Damiano, Jayesh M Goyal, Monika Lendl, Joshua D. Lothringer, Sagnick Mukherjee, Kazumasa Ohno, Natalie M. Batalha, Matthew P. Battley, Jacob L. Bean, Thomas G. Beatty, Björn Benneke, Zachory K. Berta-Thompson , et al. (74 additional authors not shown)

Abstract: Measuring the metallicity and carbon-to-oxygen (C/O) ratio in exoplanet atmospheres is a fundamental step towards constraining the dominant chemical processes at work and, if in equilibrium, revealing planet formation histories. Transmission spectroscopy provides the necessary means by constraining the abundances of oxygen- and carbon-bearing species; however, this requires broad wavelength covera… ▽ More Measuring the metallicity and carbon-to-oxygen (C/O) ratio in exoplanet atmospheres is a fundamental step towards constraining the dominant chemical processes at work and, if in equilibrium, revealing planet formation histories. Transmission spectroscopy provides the necessary means by constraining the abundances of oxygen- and carbon-bearing species; however, this requires broad wavelength coverage, moderate spectral resolution, and high precision that, together, are not achievable with previous observatories. Now that JWST has commenced science operations, we are able to observe exoplanets at previously uncharted wavelengths and spectral resolutions. Here we report time-series observations of the transiting exoplanet WASP-39b using JWST's Near InfraRed Camera (NIRCam). The long-wavelength spectroscopic and short-wavelength photometric light curves span 2.0 - 4.0 $μ$m, exhibit minimal systematics, and reveal well-defined molecular absorption features in the planet's spectrum. Specifically, we detect gaseous H$_2$O in the atmosphere and place an upper limit on the abundance of CH$_4$. The otherwise prominent CO$_2$ feature at 2.8 $μ$m is largely masked by H$_2$O. The best-fit chemical equilibrium models favour an atmospheric metallicity of 1-100$\times$ solar (i.e., an enrichment of elements heavier than helium relative to the Sun) and a sub-stellar carbon-to-oxygen (C/O) ratio. The inferred high metallicity and low C/O ratio may indicate significant accretion of solid materials during planet formation or disequilibrium processes in the upper atmosphere. △ Less

Submitted 18 November, 2022; originally announced November 2022.

Comments: 35 pages, 13 figures, 3 tables, Nature, accepted

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

How Population III Supernovae Determined the Properties of the First Galaxies

Authors: Ke-Jung Chen, Ching-Yao Tang, Daniel J. Whalen, Meng-Yuan Ho, Sung-Han Tsai, Po-Sheng Ou, Masaomi Ono

Abstract: Massive Pop III stars can die as energetic supernovae that enrich the early universe with metals and determine the properties of the first galaxies. With masses of about $10^9$ Msun at $z \gtrsim 10$, these galaxies are believed to be the ancestors of the Milky Way. This paper investigates the impact of Pop III supernova remnants (SNRs) from both Salpeter-like and top-heavy initial mass functions… ▽ More Massive Pop III stars can die as energetic supernovae that enrich the early universe with metals and determine the properties of the first galaxies. With masses of about $10^9$ Msun at $z \gtrsim 10$, these galaxies are believed to be the ancestors of the Milky Way. This paper investigates the impact of Pop III supernova remnants (SNRs) from both Salpeter-like and top-heavy initial mass functions (IMFs) on the formation of first galaxies with high-resolution radiation-hydrodynamical simulations with the ENZO code. Our findings indicate that SNRs from a top-heavy Pop III IMF produce more metals, leading to more efficient gas cooling and earlier Pop II star formation in the first galaxies. From a few hundred to a few thousand Pop II stars can form in the central regions of these galaxies. These stars have metallicities of $10^{-3}$ to $10^{-2}$, Zsun, greater than those of extremely metal-poor (EMP) stars. Their mass function follows a power-law distribution with $dN(M_*)/dM_* \propto M_*^α$, where $M_*$ is stellar mass and $α= 2.66 - 5.83$ and is steeper for a top-heavy IMF. We thus find that EMP stars were not typical of most primitive galaxies. △ Less

Submitted 12 February, 2024; v1 submitted 11 November, 2022; originally announced November 2022.

Comments: 14 pages, 9 figures. Accepted for publication in ApJ. arXiv admin note: substantial text overlap with arXiv:2010.02212

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

Critical Metallicity of Cool Supergiant Formation. I. Effects on Stellar Mass Loss and Feedback

Authors: Po-Sheng Ou, Ke-Jung Chen, You-Hua Chu, Sung-Han Tsai

Abstract: This paper systematically studies the relation between metallicity and mass loss of massive stars. We perform one-dimensional stellar evolution simulations and build a grid of $\sim$2000 models with initial masses ranging between 11 and 60 $M_{\odot}$ and absolute metallicities $Z$ between 0.00001 and 0.02. Steady-state winds, comprising hot main-sequence winds and cool supergiant winds, are the m… ▽ More This paper systematically studies the relation between metallicity and mass loss of massive stars. We perform one-dimensional stellar evolution simulations and build a grid of $\sim$2000 models with initial masses ranging between 11 and 60 $M_{\odot}$ and absolute metallicities $Z$ between 0.00001 and 0.02. Steady-state winds, comprising hot main-sequence winds and cool supergiant winds, are the main drivers of the mass loss of massive stars in our models. We calculate the total mass loss over the stellar lifetime for each model. Our results reveal the existence of a critical metallicity $Z_{\rm{c}}$ at $Z \sim 10^{-3}$, where the mass loss exhibits a dramatic jump. If $Z>Z_{\rm{c}}$, massive stars tend to evolve into cool supergiants, and a robust cool wind is operational. In contrast, if $Z<Z_{\rm{c}}$, massive stars usually remain as blue supergiants, wherein the cool wind is not activated and the mass loss is generally weak. Moreover, we calculate the wind feedback in a $10^5$ $M_{\odot}$ star cluster with the Salpeter initial mass function. The kinetic energy released by winds does not exhibit any significant transition at $Z_{\rm{c}}$ because the wind velocity of a cool supergiant wind is low and contributes little to the kinetic energy. The effects of critical metallicity provide implications for the fates of metal-poor stars in the early universe. △ Less

Submitted 19 December, 2022; v1 submitted 18 October, 2022; originally announced October 2022.

Comments: 22 pages, 12 figures, accepted for publication in the Astrophysical Journal

The Climate and Compositional Variation of the Highly Eccentric Planet HD 80606 b -- the rise and fall of carbon monoxide and elemental sulfur

Authors: Shang-Min Tsai, Maria Steinrueck, Vivien Parmentier, Nikole Lewis, Raymond Pierrehumbert

Abstract: The gas giant HD 80606 b has a highly eccentric orbit (e $\sim$ 0.93). The variation due to the rapid shift of stellar irradiation provides a unique opportunity to probe the physical and chemical timescales and to study the interplay between climate dynamics and atmospheric chemistry. In this work, we present integrated models to study the atmospheric responses and the underlying physical and chem… ▽ More The gas giant HD 80606 b has a highly eccentric orbit (e $\sim$ 0.93). The variation due to the rapid shift of stellar irradiation provides a unique opportunity to probe the physical and chemical timescales and to study the interplay between climate dynamics and atmospheric chemistry. In this work, we present integrated models to study the atmospheric responses and the underlying physical and chemical mechanisms of HD 80606 b. We first run three-dimensional general circulation models (GCMs) to establish the atmospheric thermal and dynamical structures for different atmospheric metallicities and internal heat. Based on the GCM output, we then adopted a 1D time-dependent photochemical model to investigate the compositional variation along the eccentric orbit. The transition of the circulation patterns of HD 80606 b matched the dynamics regimes in previous works. Our photochemical models show that efficient vertical mixing leads to deep quench levels of the major carbon and nitrogen species and the quenching behavior does not change throughout the eccentric orbit. Instead, photolysis is the main driver of the time-dependent chemistry. While CH$_4$ dominates over CO through most of the orbits, a transient state of [CO]/[CH$_4$}] $>$ 1 after periastron is confirmed for all metallicity and internal heat cases. The upcoming JWST Cycle 1 GO program will be able to track this real-time CH$_4$--CO conversion and infer the chemical timescale. Furthermore, sulfur species initiated by sudden heating and photochemical forcing exhibit both short-term and long-term cycles, opening an interesting avenue for detecting sulfur on exoplanets. △ Less

Submitted 16 March, 2023; v1 submitted 13 September, 2022; originally announced September 2022.

Comments: 18 pages, 18 figures, published in MNRAS

Identification of carbon dioxide in an exoplanet atmosphere

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

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

Submitted 24 August, 2022; originally announced August 2022.

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

Constraining ultralight bosonic dark matter with Keck observations of S2's orbit and kinematics

Authors: Guan-Wen Yuan, Zhao-Qiang Shen, Yue-Lin Sming Tsai, Qiang Yuan, Yi-Zhong Fan

Abstract: Ultralight bosonic dark matter is expected to be able to form a cloud surrounding the supermassive black hole (SMBH) in the Galactic center. With increasing precision of the observations of the stellar kinematics around the SMBH, tiny effects from such a dark matter cloud, including its gravitational perturbation and the direct coupling with the ordinary matter may be detectable. In this work, we… ▽ More Ultralight bosonic dark matter is expected to be able to form a cloud surrounding the supermassive black hole (SMBH) in the Galactic center. With increasing precision of the observations of the stellar kinematics around the SMBH, tiny effects from such a dark matter cloud, including its gravitational perturbation and the direct coupling with the ordinary matter may be detectable. In this work, we search for possible evidence of the scalar cloud using accurate orbital measurements of the S2 star around Sgr~A*. We solve the first order Post-Newtonian equation, considering simultaneously the extended mass distribution of the scalar cloud and the frequency shift induced by the additional coupling via Higgs portal or photon portal interaction. Furthermore, we also investigate the impact of an astrophysical power-law component from the gas and stellar remnants. We find that the astrometric and spectroscopic data of the S2 star are well consistent with the scenario of a point-like mass of Sgr~A*. We thus derive upper limits of the coupling of the new interaction and the extended mass, with and without the contribution from the astrophysical component. The limits of the Higgs/photon coupling and the extended mass of the scalar cloud are the most stringent ones for the scalar mass window between $3.2\times 10^{-19}$~eV and $1.6\times 10^{-18}$~eV. △ Less

Submitted 21 November, 2022; v1 submitted 10 May, 2022; originally announced May 2022.

Comments: 12 pages, 6 figures, published at PRD

Journal ref: Phys.Rev.D 106 (2022) 10, 103024

The impact of time-dependent stellar activity on exoplanet atmospheres

Authors: Amy J. Louca, Yamila Miguel, Shang-Min Tsai, Cynthia S. Froning, R. O. Parke Loyd, Kevin France

Abstract: M-dwarfs are thought to be hostile environments for exoplanets. Stellar events are very common on such stars. These events might cause the atmospheres of exoplanets to change significantly over time. It is not only the major stellar flare events that contribute to this disequilibrium, but the smaller flares might also affect the atmospheres in an accumulating manner. In this study, we aim to inves… ▽ More M-dwarfs are thought to be hostile environments for exoplanets. Stellar events are very common on such stars. These events might cause the atmospheres of exoplanets to change significantly over time. It is not only the major stellar flare events that contribute to this disequilibrium, but the smaller flares might also affect the atmospheres in an accumulating manner. In this study, we aim to investigate the effects of time-dependent stellar activity on the atmospheres of known exoplanets. We simulate the chemistry of GJ876c, GJ581c, and GJ832c that go from H$_2$-dominated to N$_2$-dominated atmospheres using observed stellar spectra from the MUSCLES-collaboration. We make use of the chemical kinetics code VULCAN and implement a flaring routine that stochastically generates synthetic flares based on observed flare statistics. Using the radiative transfer code petitRADTrans we also simulate the evolution of emission and transmission spectra. We investigate the effect of recurring flares for a total of 11 days covering 515 flares. Results show a significant change in abundance for some relevant species such as H, OH, and CH$_4$, with factors going up to 3 orders of magnitude difference with respect to the preflare abundances. We find a maximum change of $\sim$12 ppm for CH$_4$ in transmission spectra on GJ876c. These changes in the spectra remain too small to observe. We also find that the change in abundance and spectra of the planets accumulate throughout time, causing permanent changes in the chemistry. We conclude this small but gradual change in chemistry arises due to the recurring flares. △ Less

Submitted 22 April, 2022; originally announced April 2022.

Comments: Accepted for publication in MNRAS

NMSSM neutralino dark matter for CDF II $W$-boson mass and muon $g-2$ and the promising prospect of direct detection

Authors: Tian-Peng Tang, Murat Abdughani, Lei Feng, Yue-Lin Sming Tsai, Jian Wu, Yi-Zhong Fan

Abstract: Two experiments from the Fermilab, E989 and CDF II, have reported two anomalies for muon $g-2$ and $W$-boson mass that may indicate the new physics at the low energy scale. Here we examine the possibility of a common origin of these two anomalies in the Next-to-Minimal Supersymmetric Standard Model. Considering various experimental and astrophysical constraints such as the Higgs mass, collider dat… ▽ More Two experiments from the Fermilab, E989 and CDF II, have reported two anomalies for muon $g-2$ and $W$-boson mass that may indicate the new physics at the low energy scale. Here we examine the possibility of a common origin of these two anomalies in the Next-to-Minimal Supersymmetric Standard Model. Considering various experimental and astrophysical constraints such as the Higgs mass, collider data, flavor physics, dark matter relic density, and direct detection experiments, we find that lighter electroweakinos and sleptons can generate sufficient contributions to muon $g-2$ and $m_W$. Moreover, the corresponding bino-like neutralino dark matter mass is in the $\sim 180-280$ GeV range. Interestingly, the favored DM mass region can soon be entirely probed by ongoing direct detection experiments like PandaX-4T, XENONnT, LUX-ZEPLIN, and DARWIN. △ Less

Submitted 21 February, 2023; v1 submitted 8 April, 2022; originally announced April 2022.

Comments: 21 pages, 7 figures

Journal ref: Sci.China Phys.Mech.Astron. 66 (2023) 3, 239512

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

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

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

Submitted 8 April, 2022; originally announced April 2022.

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

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

Linking atmospheric chemistry of the hot Jupiter HD 209458b to its formation location through infrared transmission and emission spectra

Authors: Spandan Dash, Liton Majumdar, Karen Willacy, Shang-Min Tsai, Neal Turner, P. B. Rimmer, Murthy S. Gudipati, Wladimir Lyra, Anil Bhardwaj

Abstract: The elemental ratios of carbon, nitrogen, and oxygen in the atmospheres of hot Jupiters may hold clues to their formation locations in the protostellar disc. In this work, we adopt gas phase chemical abundances of C, N and O from several locations in a disc chemical kinetics model as sources for the envelope of the hot Jupiter HD 209458b and evolve the planet's atmospheric composition using a 1D c… ▽ More The elemental ratios of carbon, nitrogen, and oxygen in the atmospheres of hot Jupiters may hold clues to their formation locations in the protostellar disc. In this work, we adopt gas phase chemical abundances of C, N and O from several locations in a disc chemical kinetics model as sources for the envelope of the hot Jupiter HD 209458b and evolve the planet's atmospheric composition using a 1D chemical kinetics model, treating both vertical mixing and photochemistry. We consider two atmospheric pressure-temperature profiles, one with and one without a thermal inversion. From each of the resulting 32 atmospheric composition profiles, we find that the molecules CH4, NH3, HCN, and C2H2 are more prominent in the atmospheres computed using a realistic non-inverted P-T profile in comparison to a prior equilibrium chemistry based work which used an analytical P-T profile. We also compute the synthetic transmission and emission spectra for these atmospheres and find that many spectral features vary with the location in the disc where the planet's envelope was accreted. By comparing with the species detected using the latest high-resolution ground-based observations, our model suggests HD 209458b could have accreted most of its gas between the CO2 and CH4 icelines with a super solar C/O ratio from its protostellar disc, which in turn directly inherited its chemical abundances from the protostellar cloud. Finally, we simulate observing the planet with the James Webb Space Telescope (JWST) and show that differences in spectral signatures of key species can be recognized. Our study demonstrates the enormous importance of JWST in providing new insights into hot Jupiter's formation environments. △ Less

Submitted 11 May, 2022; v1 submitted 8 April, 2022; originally announced April 2022.

Comments: Accepted for publication in ApJ, 28 Pages, 11 figures, 2 Tables (Appendix: 5 Figures, 1 Table)

Inert Higgs Dark Matter for CDF-II W-boson Mass and Detection Prospects

Authors: Yi-Zhong Fan, Tian-Peng Tang, Yue-Lin Sming Tsai, Lei Wu

Abstract: The $W$-boson mass, which was recently measured at FermiLab with an unprecedented precision, suggests the presence of new multiplets beyond the Standard Model (SM). One of the minimal extensions of the SM is to introduce an additional scalar doublet, in which the non-SM scalars can enhance $W$-boson mass via the loop corrections. On the other hand, with a proper discrete symmetry, the lightest new… ▽ More The $W$-boson mass, which was recently measured at FermiLab with an unprecedented precision, suggests the presence of new multiplets beyond the Standard Model (SM). One of the minimal extensions of the SM is to introduce an additional scalar doublet, in which the non-SM scalars can enhance $W$-boson mass via the loop corrections. On the other hand, with a proper discrete symmetry, the lightest new scalar in the doublet can be stable and play the role of dark matter particle. We show that the inert two Higgs doublet model can naturally handle the new $W$-boson mass without violating other constraints, and the preferred dark matter mass is between $54$ and $74$ GeV. We identify three feasible parameter regions for the thermal relic density: the $SA$ co-annihilation, the Higgs resonance, and the $SS \to WW^*$ annihilation. We find that the first region can be fully tested by the HL-LHC, the second region will be tightly constrained by direct detection experiments, and the third region could yield detectable GeV gamma-ray and antiproton signals in the Galaxy that may have been observed by Fermi-LAT and AMS-02. △ Less

Submitted 25 August, 2022; v1 submitted 7 April, 2022; originally announced April 2022.

Comments: 9 pages, 7 figures, Replaced to match published version in PRL

Journal ref: Phys.Rev.Lett. 129 (2022) 9, 091802

Search for Gravitational Waves Associated with Fast Radio Bursts Detected by CHIME/FRB During the LIGO--Virgo Observing Run O3a

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, the CHIME/FRB Collaboration, :, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, N. Adhikari, R. X. Adhikari, 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, A. Allocca , et al. (1633 additional authors not shown)

Abstract: We search for gravitational-wave transients associated with fast radio bursts (FRBs) detected by the Canadian Hydrogen Intensity Mapping Experiment Fast Radio Burst Project (CHIME/FRB), during the first part of the third observing run of Advanced LIGO and Advanced Virgo (1 April 2019 15:00 UTC-1 Oct 2019 15:00 UTC). Triggers from 22 FRBs were analyzed with a search that targets compact binary coal… ▽ More We search for gravitational-wave transients associated with fast radio bursts (FRBs) detected by the Canadian Hydrogen Intensity Mapping Experiment Fast Radio Burst Project (CHIME/FRB), during the first part of the third observing run of Advanced LIGO and Advanced Virgo (1 April 2019 15:00 UTC-1 Oct 2019 15:00 UTC). Triggers from 22 FRBs were analyzed with a search that targets compact binary coalescences with at least one neutron star component. A targeted search for generic gravitational-wave transients was conducted on 40 FRBs. We find no significant evidence for a gravitational-wave association in either search. Given the large uncertainties in the distances of the FRBs inferred from the dispersion measures in our sample, however, this does not conclusively exclude any progenitor models that include emission of a gravitational wave of the types searched for from any of these FRB events. We report $90\%$ confidence lower bounds on the distance to each FRB for a range of gravitational-wave progenitor models. By combining the inferred maximum distance information for each FRB with the sensitivity of the gravitational-wave searches, we set upper limits on the energy emitted through gravitational waves for a range of emission scenarios. We find values of order $10^{51}$-$10^{57}$ erg for a range of different emission models with central gravitational wave frequencies in the range 70-3560 Hz. Finally, we also found no significant coincident detection of gravitational waves with the repeater, FRB 20200120E, which is the closest known extragalactic FRB. △ Less

Submitted 22 March, 2022; originally announced March 2022.

Comments: 35 pages, 6 figures, 8 tables

Report number: P2100124

Narrowband searches for continuous and long-duration transient gravitational waves from known pulsars in the LIGO-Virgo third observing run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, N. Adhikari, R. X. Adhikari, 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, A. Allocca, P. A. Altin, A. Amato , et al. (1636 additional authors not shown)

Abstract: Isolated neutron stars that are asymmetric with respect to their spin axis are possible sources of detectable continuous gravitational waves. This paper presents a fully-coherent search for such signals from eighteen pulsars in data from LIGO and Virgo's third observing run (O3). For known pulsars, efficient and sensitive matched-filter searches can be carried out if one assumes the gravitational… ▽ More Isolated neutron stars that are asymmetric with respect to their spin axis are possible sources of detectable continuous gravitational waves. This paper presents a fully-coherent search for such signals from eighteen pulsars in data from LIGO and Virgo's third observing run (O3). For known pulsars, efficient and sensitive matched-filter searches can be carried out if one assumes the gravitational radiation is phase-locked to the electromagnetic emission. In the search presented here, we relax this assumption and allow the frequency and frequency time-derivative of the gravitational waves to vary in a small range around those inferred from electromagnetic observations. We find no evidence for continuous gravitational waves, and set upper limits on the strain amplitude for each target. These limits are more constraining for seven of the targets than the spin-down limit defined by ascribing all rotational energy loss to gravitational radiation. In an additional search we look in O3 data for long-duration (hours-months) transient gravitational waves in the aftermath of pulsar glitches for six targets with a total of nine glitches. We report two marginal outliers from this search, but find no clear evidence for such emission either. The resulting duration-dependent strain upper limits do not surpass indirect energy constraints for any of these targets. △ Less

Submitted 27 June, 2022; v1 submitted 21 December, 2021; originally announced December 2021.

Comments: 37 pages, 9 figures, submitted to ApJ

Report number: LIGO-P2100267

Journal ref: ApJ, 932, 133 (2022)

Inferring Shallow Surfaces on sub-Neptune Exoplanets with JWST

Authors: Shang-Min Tsai, Hamish Innes, Tim Lichtenberg, Jake Taylor, Matej Malik, Katy Chubb, Raymond Pierrehumbert

Abstract: Planets smaller than Neptune and larger than Earth make up the majority of the discovered exoplanets. Those with H$_2$-rich atmospheres are prime targets for atmospheric characterization. The transition between the two main classes, super-Earths and sub-Neptunes, is not clearly understood as the rocky surface is likely not accessible to observations. Tracking several trace gases (specifically the… ▽ More Planets smaller than Neptune and larger than Earth make up the majority of the discovered exoplanets. Those with H$_2$-rich atmospheres are prime targets for atmospheric characterization. The transition between the two main classes, super-Earths and sub-Neptunes, is not clearly understood as the rocky surface is likely not accessible to observations. Tracking several trace gases (specifically the loss of ammonia (NH$_3$) and hydrogen cyanide (HCN)) has been proposed as a proxy for the presence of a shallow surface. In this work, we revisit the proposed mechanism of nitrogen conversion in detail and find its timescale on the order of a million years. NH$_3$ exhibits dual paths converting to N$_2$ or HCN, depending on the UV radiation of the star and the stage of the system. In addition, methanol (CH$_3$OH) is identified as a robust and complementary proxy for a shallow surface. We follow the fiducial example of K2-18b with a 2D photochemical model (VULCAN) on an equatorial plane. We find a fairly uniform composition distribution below 0.1 mbar controlled by the dayside, as a result of slow chemical evolution. NH$_3$ and CH$_3$OH are concluded to be the most unambiguous proxies to infer surfaces on sub-Neptunes in the era of the James Webb Space Telescope (JWST). △ Less

Submitted 11 November, 2021; originally announced November 2021.

Comments: Accepted for publication in ApJL

The population of merging compact binaries inferred using gravitational waves through GWTC-3

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, N. Adhikari, R. X. Adhikari, 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, A. Allocca, P. A. Altin, A. Amato , et al. (1612 additional authors not shown)

Abstract: We report on the population properties of 76 compact binary mergers detected with gravitational waves below a false alarm rate of 1 per year through GWTC-3. The catalog contains three classes of binary mergers: BBH, BNS, and NSBH mergers. We infer the BNS merger rate to be between 10 $\rm{Gpc^{-3} yr^{-1}}$ and 1700 $\rm{Gpc^{-3} yr^{-1}}$ and the NSBH merger rate to be between 7.8… ▽ More We report on the population properties of 76 compact binary mergers detected with gravitational waves below a false alarm rate of 1 per year through GWTC-3. The catalog contains three classes of binary mergers: BBH, BNS, and NSBH mergers. We infer the BNS merger rate to be between 10 $\rm{Gpc^{-3} yr^{-1}}$ and 1700 $\rm{Gpc^{-3} yr^{-1}}$ and the NSBH merger rate to be between 7.8 $\rm{Gpc^{-3}\, yr^{-1}}$ and 140 $\rm{Gpc^{-3} yr^{-1}}$ , assuming a constant rate density versus comoving volume and taking the union of 90% credible intervals for methods used in this work. Accounting for the BBH merger rate to evolve with redshift, we find the BBH merger rate to be between 17.9 $\rm{Gpc^{-3}\, yr^{-1}}$ and 44 $\rm{Gpc^{-3}\, yr^{-1}}$ at a fiducial redshift (z=0.2). We obtain a broad neutron star mass distribution extending from $1.2^{+0.1}_{-0.2} M_\odot$ to $2.0^{+0.3}_{-0.3} M_\odot$. We can confidently identify a rapid decrease in merger rate versus component mass between neutron star-like masses and black-hole-like masses, but there is no evidence that the merger rate increases again before 10 $M_\odot$. We also find the BBH mass distribution has localized over- and under-densities relative to a power law distribution. While we continue to find the mass distribution of a binary's more massive component strongly decreases as a function of primary mass, we observe no evidence of a strongly suppressed merger rate above $\sim 60 M_\odot$. The rate of BBH mergers is observed to increase with redshift at a rate proportional to $(1+z)^κ$ with $κ= 2.9^{+1.7}_{-1.8}$ for $z\lesssim 1$. Observed black hole spins are small, with half of spin magnitudes below $χ_i \simeq 0.25$. We observe evidence of negative aligned spins in the population, and an increase in spin magnitude for systems with more unequal mass ratio. △ Less

Submitted 23 February, 2022; v1 submitted 5 November, 2021; originally announced November 2021.

Comments: v2: minor edits, most to Table 1 and caption; v3: rerun with public data; Data release: https://zenodo.org/record/5655785; v4: update Fig 14

Report number: LIGO-P2100239

Search for Gravitational Waves Associated with Gamma-Ray Bursts Detected by Fermi and Swift During the LIGO-Virgo Run O3b

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, N. Adhikari, R. X. Adhikari, 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, A. Allocca, P. A. Altin, A. Amato , et al. (1610 additional authors not shown)

Abstract: We search for gravitational-wave signals associated with gamma-ray bursts detected by the Fermi and Swift satellites during the second half of the third observing run of Advanced LIGO and Advanced Virgo (1 November 2019 15:00 UTC-27 March 2020 17:00 UTC).We conduct two independent searches: a generic gravitational-wave transients search to analyze 86 gamma-ray bursts and an analysis to target bina… ▽ More We search for gravitational-wave signals associated with gamma-ray bursts detected by the Fermi and Swift satellites during the second half of the third observing run of Advanced LIGO and Advanced Virgo (1 November 2019 15:00 UTC-27 March 2020 17:00 UTC).We conduct two independent searches: a generic gravitational-wave transients search to analyze 86 gamma-ray bursts and an analysis to target binary mergers with at least one neutron star as short gamma-ray burst progenitors for 17 events. We find no significant evidence for gravitational-wave signals associated with any of these gamma-ray bursts. A weighted binomial test of the combined results finds no evidence for sub-threshold gravitational wave signals associated with this GRB ensemble either. We use several source types and signal morphologies during the searches, resulting in lower bounds on the estimated distance to each gamma-ray burst. Finally, we constrain the population of low luminosity short gamma-ray bursts using results from the first to the third observing runs of Advanced LIGO and Advanced Virgo. The resulting population is in accordance with the local binary neutron star merger rate. △ Less

Submitted 5 November, 2021; originally announced November 2021.

Comments: 26 pages, 6 figures, 2 tables

Report number: P2100091

GWTC-3: Compact Binary Coalescences Observed by LIGO and Virgo During the Second Part of the Third Observing Run

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

Abstract: The third Gravitational-Wave Transient Catalog (GWTC-3) describes signals detected with Advanced LIGO and Advanced Virgo up to the end of their third observing run. Updating the previous GWTC-2.1, we present candidate gravitational waves from compact binary coalescences during the second half of the third observing run (O3b) between 1 November 2019, 15:00 UTC and 27 March 2020, 17:00 UTC. There ar… ▽ More The third Gravitational-Wave Transient Catalog (GWTC-3) describes signals detected with Advanced LIGO and Advanced Virgo up to the end of their third observing run. Updating the previous GWTC-2.1, we present candidate gravitational waves from compact binary coalescences during the second half of the third observing run (O3b) between 1 November 2019, 15:00 UTC and 27 March 2020, 17:00 UTC. There are 35 compact binary coalescence candidates identified by at least one of our search algorithms with a probability of astrophysical origin $p_\mathrm{astro} > 0.5$. Of these, 18 were previously reported as low-latency public alerts, and 17 are reported here for the first time. Based upon estimates for the component masses, our O3b candidates with $p_\mathrm{astro} > 0.5$ are consistent with gravitational-wave signals from binary black holes or neutron star-black hole binaries, and we identify none from binary neutron stars. However, from the gravitational-wave data alone, we are not able to measure matter effects that distinguish whether the binary components are neutron stars or black holes. The range of inferred component masses is similar to that found with previous catalogs, but the O3b candidates include the first confident observations of neutron star-black hole binaries. Including the 35 candidates from O3b in addition to those from GWTC-2.1, GWTC-3 contains 90 candidates found by our analysis with $p_\mathrm{astro} > 0.5$ across the first three observing runs. These observations of compact binary coalescences present an unprecedented view of the properties of black holes and neutron stars. △ Less

Submitted 23 October, 2023; v1 submitted 5 November, 2021; originally announced November 2021.

Comments: 88 pages (10 pages author list, 31 pages main text, 1 page acknowledgements, 24 pages appendices, 22 pages bibliography), 17 figures, 16 tables. Update to match version to be published in Physical Review X. Data products available from https://gwosc.org/GWTC-3/

Report number: LIGO-P2000318

Journal ref: Phys. Rev. X; 13(4):041039; 2023

Muonphilic Dark Matter explanation of gamma-ray galactic center excess: a comprehensive analysis

Authors: Murat Abdughani, Yi-Zhong Fan, Chih-Ting Lu, Tian-Peng Tang, Yue-Lin Sming Tsai

Abstract: The Galactic center gamma-ray excess (GCE) is a long-standing unsolved problem. One of candidate solutions, the dark matter (DM) annihilation, has been recently tested with other astrophysical observations, such as AMS-02 electron-positron spectra, Fermi Dwarf spheroidal galaxies gamma-ray data, and so on. By assuming that the DM particles annihilate purely into a normal charged fermion pair, Di M… ▽ More The Galactic center gamma-ray excess (GCE) is a long-standing unsolved problem. One of candidate solutions, the dark matter (DM) annihilation, has been recently tested with other astrophysical observations, such as AMS-02 electron-positron spectra, Fermi Dwarf spheroidal galaxies gamma-ray data, and so on. By assuming that the DM particles annihilate purely into a normal charged fermion pair, Di Mauro and Winkle (2021) claimed that only a muon-pair is compatible with the null detection of all the corresponding astrophysical measurements and can explain GCE simultaneously. On the other hand, a muonphilic DM model may also lead to a signal in the recent Fermilab muon $g-2$ measurement or be constrained by the latest PandaX-4T limit. In this work, we comprehensively study interactions between DM and muon, including various combinations of DM and mediator spins. In agreement with GCE (not only $2μ$ but also $4μ$ final states), we test these interactions against all the thermal DM constraints. Our results show that only the parameter space near the resonance region of mediator can explain GCE and relic density simultaneously, and larger parameter spaces are still allowed if other poorly-known systematic uncertainties are included. Regardless of the DM spin, only the interactions with the spin-0 mediator can explain the recent muon $g-2$ excess on top of GCE, relic density, and other DM and mediator constraints. △ Less

Submitted 15 July, 2022; v1 submitted 4 November, 2021; originally announced November 2021.

Comments: 34 pages, 8 figures, Matches JHEP accepted version. Major revisions to address referee's comments