-
Revisiting the Mysterious Origin of FRB 20121102A with Machine-learning Classification
Authors:
Leah Ya-Ling Lin,
Tetsuya Hashimoto,
Tomotsugu Goto,
Bjorn Jasper Raquel,
Simon C. -C. Ho,
Bo-Han Chen,
Seong Jin Kim,
Chih-Teng Ling
Abstract:
Fast radio bursts (FRBs) are millisecond-duration radio waves from the Universe. Even though more than 50 physical models have been proposed, the origin and physical mechanism of FRB emissions are still unknown. The classification of FRBs is one of the primary approaches to understanding their mechanisms, but previous studies classified conventionally using only a few observational parameters, suc…
▽ More
Fast radio bursts (FRBs) are millisecond-duration radio waves from the Universe. Even though more than 50 physical models have been proposed, the origin and physical mechanism of FRB emissions are still unknown. The classification of FRBs is one of the primary approaches to understanding their mechanisms, but previous studies classified conventionally using only a few observational parameters, such as fluence and duration, which might be incomplete. To overcome this problem, we use an unsupervised machine-learning model, the Uniform Manifold Approximation and Projection (UMAP) to handle seven parameters simultaneously, including amplitude, linear temporal drift, time duration, central frequency, bandwidth, scaled energy, and fluence.
We test the method for homogeneous 977 sub-bursts of FRB 20121102A detected by the Arecibo telescope. Our machine-learning analysis identified five distinct clusters, suggesting the possible existence of multiple different physical mechanisms responsible for the observed FRBs from the FRB 20121102A source. The geometry of the emission region and the propagation effect of FRB signals could also make such distinct clusters.
This research will be a benchmark for future FRB classifications when dedicated radio telescopes such as the Square Kilometer Array (SKA) or Bustling Universe Radio Survey Telescope in Taiwan (BURSTT) discover more FRBs than before.
△ Less
Submitted 1 October, 2024;
originally announced October 2024.
-
Accelerating Giant Impact Simulations with Machine Learning
Authors:
Caleb Lammers,
Miles Cranmer,
Sam Hadden,
Shirley Ho,
Norman Murray,
Daniel Tamayo
Abstract:
Constraining planet formation models based on the observed exoplanet population requires generating large samples of synthetic planetary systems, which can be computationally prohibitive. A significant bottleneck is simulating the giant impact phase, during which planetary embryos evolve gravitationally and combine to form planets, which may themselves experience later collisions. To accelerate gi…
▽ More
Constraining planet formation models based on the observed exoplanet population requires generating large samples of synthetic planetary systems, which can be computationally prohibitive. A significant bottleneck is simulating the giant impact phase, during which planetary embryos evolve gravitationally and combine to form planets, which may themselves experience later collisions. To accelerate giant impact simulations, we present a machine learning (ML) approach to predicting collisional outcomes in multiplanet systems. Trained on more than 500,000 $N$-body simulations of three-planet systems, we develop an ML model that can accurately predict which two planets will experience a collision, along with the state of the post-collision planets, from a short integration of the system's initial conditions. Our model greatly improves on non-ML baselines that rely on metrics from dynamics theory, which struggle to accurately predict which pair of planets will experience a collision. By combining with a model for predicting long-term stability, we create an ML-based giant impact emulator, which can predict the outcomes of giant impact simulations with reasonable accuracy and a speedup of up to four orders of magnitude. We expect our model to enable analyses that would not otherwise be computationally feasible. As such, we release our training code, along with an easy-to-use API for our collision outcome model and giant impact emulator.
△ Less
Submitted 25 September, 2024; v1 submitted 16 August, 2024;
originally announced August 2024.
-
Field-level Emulation of Cosmic Structure Formation with Cosmology and Redshift Dependence
Authors:
Drew Jamieson,
Yin Li,
Francisco Villaescusa-Navarro,
Shirley Ho,
David N. Spergel
Abstract:
We present a field-level emulator for large-scale structure, capturing the cosmology dependence and the time evolution of cosmic structure formation. The emulator maps linear displacement fields to their corresponding nonlinear displacements from N-body simulations at specific redshifts. Designed as a neural network, the emulator incorporates style parameters that encode dependencies on…
▽ More
We present a field-level emulator for large-scale structure, capturing the cosmology dependence and the time evolution of cosmic structure formation. The emulator maps linear displacement fields to their corresponding nonlinear displacements from N-body simulations at specific redshifts. Designed as a neural network, the emulator incorporates style parameters that encode dependencies on $Ω_{\rm m}$ and the linear growth factor $D(z)$ at redshift $z$. We train our model on the six-dimensional N-body phase space, predicting particle velocities as the time derivative of the model's displacement outputs. This innovation results in significant improvements in training efficiency and model accuracy. Tested on diverse cosmologies and redshifts not seen during training, the emulator achieves percent-level accuracy on scales of $k\sim~1~{\rm Mpc}^{-1}~h$ at $z=0$, with improved performance at higher redshifts. We compare predicted structure formation histories with N-body simulations via merger trees, finding consistent merger event sequences and statistical properties.
△ Less
Submitted 14 August, 2024;
originally announced August 2024.
-
Finding dusty AGNs from the JWST CEERS survey with mid-infrared photometry
Authors:
Tom C. -C. Chien,
Chih-Teng Ling,
Tomotsugu Goto,
Cossas K. -W. Wu,
Seong Jin Kim,
Tetsuya Hashimoto,
Yu-Wei Lin,
Ece Kilerci,
Simon C. -C. Ho,
Po-Ya Wang,
Bjorn Jasper R. Raquel
Abstract:
The nature of the interaction between active galactic nuclei (AGNs) and their host galaxies remains an unsolved question. Therefore, conducting an AGN census is valuable to AGN research. Nevertheless, a significant fraction of AGNs are obscured by their environment, which blocks UV and optical emissions due to the dusty torus surrounding the central supermassive black hole (SMBH). To overcome this…
▽ More
The nature of the interaction between active galactic nuclei (AGNs) and their host galaxies remains an unsolved question. Therefore, conducting an AGN census is valuable to AGN research. Nevertheless, a significant fraction of AGNs are obscured by their environment, which blocks UV and optical emissions due to the dusty torus surrounding the central supermassive black hole (SMBH). To overcome this challenge, mid-infrared (IR) surveys have emerged as a valuable tool for identifying obscured AGNs, as the obscured light is re-emitted in this range. With its high sensitivity, the James Webb Space Telescope (JWST) uncovered more fainter objects than previous telescopes. By applying the SED fitting, this work investigates AGN candidates in JWST Cosmic Evolution Early Release Science (CEERS) fields. We identified 42 candidates, 30 of them are classified as composites ($0.2\leq f_{\rm AGN, IR}< 0.5$), and 12 of them are AGNs ($f_{\rm AGN, IR}\geq 0.5$). We report the AGN luminosity contributions and AGN number fractions as a function of redshift and total infrared luminosity, showing that previously reported increasing relations are not apparent in our sample due to the sample size. We also extend the previous results on ultra-luminous infrared galaxies (ULIRGs, $L_{\rm TIR}\geq 10^{12} L_{\odot}$) to less luminous AGNs, highlighting the power of JWST.
△ Less
Submitted 21 June, 2024;
originally announced June 2024.
-
{\sc SimBIG}: Cosmological Constraints using Simulation-Based Inference of Galaxy Clustering with Marked Power Spectra
Authors:
Elena Massara,
ChangHoon Hahn,
Michael Eickenberg,
Shirley Ho,
Jiamin Hou,
Pablo Lemos,
Chirag Modi,
Azadeh Moradinezhad Dizgah,
Liam Parker,
Bruno Régaldo-Saint Blancard
Abstract:
We present the first $Λ$CDM cosmological analysis performed on a galaxy survey using marked power spectra. The marked power spectrum is the two-point function of a marked field, where galaxies are weighted by a function that depends on their local density. The presence of the mark leads these statistics to contain higher-order information of the original galaxy field, making them a good candidate…
▽ More
We present the first $Λ$CDM cosmological analysis performed on a galaxy survey using marked power spectra. The marked power spectrum is the two-point function of a marked field, where galaxies are weighted by a function that depends on their local density. The presence of the mark leads these statistics to contain higher-order information of the original galaxy field, making them a good candidate to exploit the non-Gaussian information of a galaxy catalog. In this work we make use of \simbig, a forward modeling framework for galaxy clustering analyses, and perform simulation-based inference using normalizing flows to infer the posterior distribution of the $Λ$CDM cosmological parameters. We consider different mark configurations (ways to weight the galaxy field) and deploy them in the \simbig~pipeline to analyze the corresponding marked power spectra measured from a subset of the BOSS galaxy sample. We analyze the redshift-space mark power spectra decomposed in $\ell = 0, 2, 4$ multipoles and include scales up to the non-linear regime. Among the various mark configurations considered, the ones that give the most stringent cosmological constraints produce posterior median and $68\%$ confidence limits on the growth of structure parameters equal to $Ω_m=0.273^{+0.040}_{-0.030}$ and $σ_8=0.777^{+0.077}_{-0.071}$. Compared to a perturbation theory analysis using the power spectrum of the same dataset, the \simbig~marked power spectra constraints on $σ_8$ are up to $1.2\times$ tighter, while no improvement is seen for the other cosmological parameters.
△ Less
Submitted 5 April, 2024;
originally announced April 2024.
-
Effects of galaxy environment on merger fraction
Authors:
W. J. Pearson,
D. J. D. Santos,
T. Goto,
T. -C. Huang,
S. J. Kim,
H. Matsuhara,
A. Pollo,
S. C. -C. Ho,
H. S. Hwang,
K. Małek,
T. Nakagawa,
M. Romano,
S. Serjeant,
L. Suelves,
H. Shim,
G. J. White
Abstract:
Aims. In this work, we intend to examine how environment influences the merger fraction, from the low density field environment to higher density groups and clusters. We also aim to study how the properties of a group or cluster, as well as the position of a galaxy in the group or cluster, influences the merger fraction.
Methods. We identified galaxy groups and clusters in the North Ecliptic Pol…
▽ More
Aims. In this work, we intend to examine how environment influences the merger fraction, from the low density field environment to higher density groups and clusters. We also aim to study how the properties of a group or cluster, as well as the position of a galaxy in the group or cluster, influences the merger fraction.
Methods. We identified galaxy groups and clusters in the North Ecliptic Pole using a friends-of-friends algorithm and the local density. Once identified, we determined the central galaxies, group radii, velocity dispersions, and group masses of these groups and clusters. Merging systems were identified with a neural network as well as visually. With these, we examined how the merger fraction changes as the local density changes for all galaxies as well as how the merger fraction changes as the properties of the groups or clusters change.
Results. We find that the merger fraction increases as local density increases and decreases as the velocity dispersion increases, as is often found in literature. A decrease in merger fraction as the group mass increases is also found. We also find groups with larger radii have higher merger fractions. The number of galaxies in a group does not influence the merger fraction.
Conclusions. The decrease in merger fraction as group mass increases is a result of the link between group mass and velocity dispersion. Hence, this decrease of merger fraction with increasing mass is a result of the decrease of merger fraction with velocity dispersion. The increasing relation between group radii and merger fraction may be a result of larger groups having smaller velocity dispersion at a larger distance from the centre or larger groups hosting smaller, infalling groups with more mergers. However, we do not find evidence of smaller groups having higher merger fractions.
△ Less
Submitted 18 March, 2024;
originally announced March 2024.
-
Cloud-by-cloud Multiphase Investigation of the Circumgalactic Medium of Low-redshift Galaxies
Authors:
Sameer,
Jane C. Charlton,
Bart P. Wakker,
Glenn G. Kacprzak,
Nikole M. Nielsen,
Christopher W. Churchill,
Philipp Richter,
Sowgat Muzahid,
Stephanie H. Ho,
Hasti Nateghi,
Benjamin Rosenwasser,
Anand Narayanan,
Rajib Ganguly
Abstract:
The pervasive presence of warm gas in galaxy halos suggests that the circumgalactic medium (CGM) is multiphase in its ionization structure and complex in its kinematics. Some recent state-of-the-art cosmological galaxy simulations predict an azimuthal dependence of CGM metallicities. We investigate the presence of such a trend by analyzing the distribution of gas properties in the CGM around 47…
▽ More
The pervasive presence of warm gas in galaxy halos suggests that the circumgalactic medium (CGM) is multiphase in its ionization structure and complex in its kinematics. Some recent state-of-the-art cosmological galaxy simulations predict an azimuthal dependence of CGM metallicities. We investigate the presence of such a trend by analyzing the distribution of gas properties in the CGM around 47 $z <$ 0.7 galaxies from the Multiphase Galaxy Halos Survey determined using a cloud-by-cloud, multiphase, ionization modelling approach. We identify three distinct populations of absorbers: cool clouds ($T \sim$ 10$^{4.1}$ K) in photoionization equilibrium, warm-hot collisionally ionized clouds ($T \sim$ 10$^{4.5-5}$ K) affected by time-dependent photoionization, and hotter clouds ($T \sim$ 10$^{5.4-6}$ K) with broad OVI and Lya absorption consistent with collisional ionization. We find that fragmentation can play a role in the origin of cool clouds, that warm-hot clouds are out of equilibrium due to rapid cooling, and that hotter clouds are representative of virialized halo gas in all but the lowest mass galaxies. The metallicities of clouds do not depend on the azimuthal angle or other galaxy properties for any of these populations. At face value, this disagrees with the simplistic model of the CGM with bipolar outflows and cold-mode planar accretion. However, the number of clouds per sightline is significantly larger close to the minor and major axes. This implies that the processes of outflows and accretion are contributing to these CGM cloud populations, and our sightlines are probing gas of mixed origins at all azimuthal angles in these low redshift galaxies.
△ Less
Submitted 5 April, 2024; v1 submitted 8 March, 2024;
originally announced March 2024.
-
Exploring the faintest end of mid-infrared luminosity functions up to $z\simeq 5$ with the JWST CEERS survey
Authors:
Chih-Teng Ling,
Tomotsugu Goto,
Seong Jin Kim,
Cossas K. -W. Wu,
Tetsuya Hashimoto,
Tom C. -C. Chien,
Yu-Wei Lin,
Simon C. -C. Ho,
Ece Kilerci
Abstract:
Mid-infrared (MIR) light from galaxies is sensitive to dust-obscured star-formation activities because it traces the characteristic emission of dust heated by young, massive stars. By constructing the MIR luminosity functions (LFs), we are able to quantify the overall dusty star formation history and the evolution of galaxies over cosmic time. In this work, we report the first rest-frame MIR LFs a…
▽ More
Mid-infrared (MIR) light from galaxies is sensitive to dust-obscured star-formation activities because it traces the characteristic emission of dust heated by young, massive stars. By constructing the MIR luminosity functions (LFs), we are able to quantify the overall dusty star formation history and the evolution of galaxies over cosmic time. In this work, we report the first rest-frame MIR LFs at 7.7, 10, 12.8, 15, 18, and 21 $μ$m as well as the total IR LF from the James Webb Space Telescope (JWST) Cosmic Evolution Early Release Science (CEERS) survey. We identify 506 galaxies at $z=0-5.1$ in the CEERS survey that also have optical photometry from the Hubble Space Telescope. With the unprecedented sensitivity of the JWST, we probe the faintest end of the LFs at $z=0-1$ down to $L^* \sim 10^7 L_\odot$, $\sim 2$ orders of magnitude fainter than those from the previous generation of IR space telescopes. Our findings connect well with and continue the faint end of the MIR LFs from the deepest observations in past works. As a proxy of star formation history, we present the MIR-based luminosity density up to $z\simeq4.0$, marking the first probe of the early Universe by JWST MIRI.
△ Less
Submitted 7 February, 2024;
originally announced February 2024.
-
${\rm S{\scriptsize IM}BIG}$: Cosmological Constraints from the Redshift-Space Galaxy Skew Spectra
Authors:
Jiamin Hou,
Azadeh Moradinezhad Dizgah,
ChangHoon Hahn,
Michael Eickenberg,
Shirley Ho,
Pablo Lemos,
Elena Massara,
Chirag Modi,
Liam Parker,
Bruno Régaldo-Saint Blancard
Abstract:
Extracting the non-Gaussian information of the cosmic large-scale structure (LSS) is vital in unlocking the full potential of the rich datasets from the upcoming stage-IV galaxy surveys. Galaxy skew spectra serve as efficient beyond-two-point statistics, encapsulating essential bispectrum information with computational efficiency akin to power spectrum analysis. This paper presents the first cosmo…
▽ More
Extracting the non-Gaussian information of the cosmic large-scale structure (LSS) is vital in unlocking the full potential of the rich datasets from the upcoming stage-IV galaxy surveys. Galaxy skew spectra serve as efficient beyond-two-point statistics, encapsulating essential bispectrum information with computational efficiency akin to power spectrum analysis. This paper presents the first cosmological constraints from analyzing the full set of redshift-space galaxy skew spectra of the data from the SDSS-III BOSS, accessing cosmological information down to nonlinear scales. Employing the ${\rm S{\scriptsize IM}BIG}$ forward modeling framework and simulation-based inference via normalizing flows, we analyze the CMASS-SGC sub-sample, which constitute approximately 10\% of the full BOSS data. Analyzing the scales up to $k_{\rm max}=0.5 \, {\rm Mpc}^{-1}h$, we find that the skew spectra improve the constraints on $Ω_{\rm m}, Ω_{\rm b}, h$, and $n_s$ by 34\%, 35\%, 18\%, 10\%, respectively, compared to constraints from previous ${\rm S{\scriptsize IM}BIG}$ power spectrum multipoles analysis, yielding $Ω_{\rm m}=0.288^{+0.024}_{-0.034}$, $Ω_{\rm b}= 0.043^{+0.005}_{-0.007}$, $h=0.759^{+0.104}_{-0.050}$, $n_{\rm s} = 0.918^{+0.041}_{-0.090}$ (at 68\% confidence limit). On the other hand, the constraints on $σ_8$ are weaker than from the power spectrum. Including the Big Bang Nucleosynthesis (BBN) prior on baryon density reduces the uncertainty on the Hubble parameter further, achieving $h=0.750^{+0.034}_{-0.032}$, which is a 38\% improvement over the constraint from the power spectrum with the same prior. Compared to the ${\rm S{\scriptsize IM}BIG}$ bispectrum (monopole) analysis, skew spectra offer comparable constraints on larger scales ($k_{\rm max}<0.3\, {\rm Mpc}^{-1}h$) for most parameters except for $σ_8$.
△ Less
Submitted 26 January, 2024;
originally announced January 2024.
-
Shallower radius valley around low-mass hosts: Evidence for icy planets, collisions or high-energy radiation scatter
Authors:
Cynthia S. K. Ho,
James G. Rogers,
Vincent Van Eylen,
James E. Owen,
Hilke E. Schlichting
Abstract:
The radius valley, i.e., a dearth of planets with radii between 1.5 and 2 Earth radii, provides insights into planetary formation and evolution. Using homogenously revised planetary parameters from Kepler 1-minute short cadence light curves, we remodel transits of 72 small planets mostly orbiting low-mass stars, improving the precision and accuracy of planet parameters. By combining this sample wi…
▽ More
The radius valley, i.e., a dearth of planets with radii between 1.5 and 2 Earth radii, provides insights into planetary formation and evolution. Using homogenously revised planetary parameters from Kepler 1-minute short cadence light curves, we remodel transits of 72 small planets mostly orbiting low-mass stars, improving the precision and accuracy of planet parameters. By combining this sample with a similar sample of planets around higher-mass stars, we determine the depth of the radius valley as a function of stellar mass. We find that the radius valley is shallower for low-mass stars compared to their higher mass counterparts. Upon comparison, we find that theoretical models of photoevaporation under-predict the number of planets observed inside the radius valley for low-mass stars: with decreasing stellar mass, the predicted fraction of planets inside the valley remains approximately constant whereas the observed fraction increases. We argue that this provides evidence for the presence of icy planets around low-mass stars. Alternatively, planets orbiting low-mass stars undergo more frequent collisions and scatter in the stars' high-energy output may also cause planets to fill the valley. We predict that more precise mass measurements for planets orbiting low mass stars would be able to distinguish between these scenarios.
△ Less
Submitted 11 June, 2024; v1 submitted 22 January, 2024;
originally announced January 2024.
-
Polycyclic aromatic hydrocarbon (PAH) luminous galaxies in JWST CEERS data
Authors:
Yu-Wei Lin,
Cossas K. -W. Wu,
Chih-Teng Ling,
Tomotsugu Goto,
Seong Jin Kim,
Ece Kilerci,
Tetsuya Hashimoto,
Po-Ya Wang,
Simon C. -C. Ho,
Tiger Yu-Yang Hsiao,
Bjorn Jasper R. Raquel,
Yuri Uno
Abstract:
It has been an unanswered question how many dusty galaxies have been undetected from the state-of-the-art observational surveys. JWST enables us to detect faint IR galaxies that have prominent polycyclic aromatic hydrocarbon (PAH) features in the mid-IR wavelengths. PAH is a valuable tracer of star formation and dust properties in the mid-infrared wavelength. The JWST Cosmic Evolution Early Releas…
▽ More
It has been an unanswered question how many dusty galaxies have been undetected from the state-of-the-art observational surveys. JWST enables us to detect faint IR galaxies that have prominent polycyclic aromatic hydrocarbon (PAH) features in the mid-IR wavelengths. PAH is a valuable tracer of star formation and dust properties in the mid-infrared wavelength. The JWST Cosmic Evolution Early Release Science (CEERS) fields provide us with wavelength coverage from 7.7 to 21 $μ$m using six photometric bands of the mid-infrared instrument (MIRI). We have identified galaxies dominated by mid-IR emission from PAHs, termed PAH galaxies. From our multi-band photometry catalogue, we selected ten PAH galaxies displaying high flux ratios of $\log(S_{15}/S_{10}) > 0.8$. The SED fitting analysis indicates that these galaxies are star-forming galaxies with total IR luminosities of $10^{10}$ $\sim$ $10^{11.5}$ $L_{\odot}$ at z $\sim 1$. The morphology of PAH galaxies does not show any clear signatures of major merging or interaction within the MIRI resolution. The majority of them are on the star-formation main sequence at $z \sim 1$. Our result demonstrates that JWST can detect PAH emissions from normal star-forming galaxies at $z \sim 1$, in addition to ultra-luminous infrared galaxies (ULIRGs) or luminous infrared galaxies (LIRGs).
△ Less
Submitted 2 January, 2024;
originally announced January 2024.
-
Cosmic star-formation history and black hole accretion history inferred from the JWST mid-infrared source counts
Authors:
Seong Jin Kim,
Tomotsugu Goto,
Chih-Teng Ling,
Cossas K. -W. Wu,
Tetsuya Hashimoto,
Ece Kilerci,
Simon C. -C. Ho,
Yuri Uno,
Po-Ya Wang,
Yu-Wei Lin
Abstract:
With the advent of the James Webb Space Telescope (JWST), extra-galactic source count studies were conducted down to sub-microJy in the mid-infrared (MIR), which is several tens of times fainter than what the previous-generation infrared (IR) telescopes achieved in the MIR. In this work, we aim to interpret the JWST source counts and constrain cosmic star-formation history (CSFH) and black hole ac…
▽ More
With the advent of the James Webb Space Telescope (JWST), extra-galactic source count studies were conducted down to sub-microJy in the mid-infrared (MIR), which is several tens of times fainter than what the previous-generation infrared (IR) telescopes achieved in the MIR. In this work, we aim to interpret the JWST source counts and constrain cosmic star-formation history (CSFH) and black hole accretion history (BHAH). We employ the backward evolution of local luminosity functions (LLFs) of galaxies to reproduce the observed source counts from sub-microJy to a few tens of mJy in the MIR bands of the JWST. The shapes of the LLFs at the MIR bands are determined using the model templates of the spectral energy distributions (SEDs) for five representative galaxy types (star-forming galaxies, starbursts, composite, AGN type 2 and 1). By simultaneously fitting our model to all the source counts in the six MIR bands, along with the previous results, we determine the best-fit evolutions of MIR LFs for each of the five galaxy types, and subsequently estimate the CSFH and BHAH. Thanks to the JWST, our estimates are based on several tens of times fainter MIR sources, the existence of which was merely an extrapolation in previous studies.
△ Less
Submitted 4 December, 2023;
originally announced December 2023.
-
Leaving No Branches Behind: Predicting Baryonic Properties of Galaxies from Merger Trees
Authors:
Chen-Yu Chuang,
Christian Kragh Jespersen,
Yen-Ting Lin,
Shirley Ho,
Shy Genel
Abstract:
Galaxies play a key role in our endeavor to understand how structure formation proceeds in the Universe. For any precision study of cosmology or galaxy formation, there is a strong demand for huge sets of realistic mock galaxy catalogs, spanning cosmologically significant volumes. For such a daunting task, methods that can produce a direct mapping between dark matter halos from dark matter-only si…
▽ More
Galaxies play a key role in our endeavor to understand how structure formation proceeds in the Universe. For any precision study of cosmology or galaxy formation, there is a strong demand for huge sets of realistic mock galaxy catalogs, spanning cosmologically significant volumes. For such a daunting task, methods that can produce a direct mapping between dark matter halos from dark matter-only simulations and galaxies are strongly preferred, as producing mocks from full-fledged hydrodynamical simulations or semi-analytical models is too expensive. Here we present a Graph Neural Network-based model that is able to accurately predict key properties of galaxies such as stellar mass, $g-r$ color, star formation rate, gas mass, stellar metallicity, and gas metallicity, purely from dark matter properties extracted from halos along the full assembly history of the galaxies. Tests based on the TNG300 simulation of the IllustrisTNG project show that our model can recover the baryonic properties of galaxies to high accuracy, over a wide redshift range ($z = 0-5$), for all galaxies with stellar masses more massive than $10^9\,M_\odot$ and their progenitors, with strong improvements over the state-of-the-art methods. We further show that our method makes substantial strides toward providing an understanding of the implications of the IllustrisTNG galaxy formation model.
△ Less
Submitted 15 November, 2023;
originally announced November 2023.
-
Surrogate Modeling for Computationally Expensive Simulations of Supernovae in High-Resolution Galaxy Simulations
Authors:
Keiya Hirashima,
Kana Moriwaki,
Michiko S. Fujii,
Yutaka Hirai,
Takayuki R. Saitoh,
Junichiro Makino,
Shirley Ho
Abstract:
Some stars are known to explode at the end of their lives, called supernovae (SNe). The substantial amount of matter and energy that SNe release provides significant feedback to star formation and gas dynamics in a galaxy. SNe release a substantial amount of matter and energy to the interstellar medium, resulting in significant feedback to star formation and gas dynamics in a galaxy. While such fe…
▽ More
Some stars are known to explode at the end of their lives, called supernovae (SNe). The substantial amount of matter and energy that SNe release provides significant feedback to star formation and gas dynamics in a galaxy. SNe release a substantial amount of matter and energy to the interstellar medium, resulting in significant feedback to star formation and gas dynamics in a galaxy. While such feedback has a crucial role in galaxy formation and evolution, in simulations of galaxy formation, it has only been implemented using simple {\it sub-grid models} instead of numerically solving the evolution of gas elements around SNe in detail due to a lack of resolution. We develop a method combining machine learning and Gibbs sampling to predict how a supernova (SN) affects the surrounding gas. The fidelity of our model in the thermal energy and momentum distribution outperforms the low-resolution SN simulations. Our method can replace the SN sub-grid models and help properly simulate un-resolved SN feedback in galaxy formation simulations. We find that employing our new approach reduces the necessary computational cost to $\sim$ 1 percent compared to directly resolving SN feedback.
△ Less
Submitted 14 November, 2023;
originally announced November 2023.
-
SimBIG: Field-level Simulation-Based Inference of Galaxy Clustering
Authors:
Pablo Lemos,
Liam Parker,
ChangHoon Hahn,
Shirley Ho,
Michael Eickenberg,
Jiamin Hou,
Elena Massara,
Chirag Modi,
Azadeh Moradinezhad Dizgah,
Bruno Regaldo-Saint Blancard,
David Spergel
Abstract:
We present the first simulation-based inference (SBI) of cosmological parameters from field-level analysis of galaxy clustering. Standard galaxy clustering analyses rely on analyzing summary statistics, such as the power spectrum, $P_\ell$, with analytic models based on perturbation theory. Consequently, they do not fully exploit the non-linear and non-Gaussian features of the galaxy distribution.…
▽ More
We present the first simulation-based inference (SBI) of cosmological parameters from field-level analysis of galaxy clustering. Standard galaxy clustering analyses rely on analyzing summary statistics, such as the power spectrum, $P_\ell$, with analytic models based on perturbation theory. Consequently, they do not fully exploit the non-linear and non-Gaussian features of the galaxy distribution. To address these limitations, we use the {\sc SimBIG} forward modelling framework to perform SBI using normalizing flows. We apply SimBIG to a subset of the BOSS CMASS galaxy sample using a convolutional neural network with stochastic weight averaging to perform massive data compression of the galaxy field. We infer constraints on $Ω_m = 0.267^{+0.033}_{-0.029}$ and $σ_8=0.762^{+0.036}_{-0.035}$. While our constraints on $Ω_m$ are in-line with standard $P_\ell$ analyses, those on $σ_8$ are $2.65\times$ tighter. Our analysis also provides constraints on the Hubble constant $H_0=64.5 \pm 3.8 \ {\rm km / s / Mpc}$ from galaxy clustering alone. This higher constraining power comes from additional non-Gaussian cosmological information, inaccessible with $P_\ell$. We demonstrate the robustness of our analysis by showcasing our ability to infer unbiased cosmological constraints from a series of test simulations that are constructed using different forward models than the one used in our training dataset. This work not only presents competitive cosmological constraints but also introduces novel methods for leveraging additional cosmological information in upcoming galaxy surveys like DESI, PFS, and Euclid.
△ Less
Submitted 23 October, 2023;
originally announced October 2023.
-
Galaxy Clustering Analysis with SimBIG and the Wavelet Scattering Transform
Authors:
Bruno Régaldo-Saint Blancard,
ChangHoon Hahn,
Shirley Ho,
Jiamin Hou,
Pablo Lemos,
Elena Massara,
Chirag Modi,
Azadeh Moradinezhad Dizgah,
Liam Parker,
Yuling Yao,
Michael Eickenberg
Abstract:
The non-Gaussisan spatial distribution of galaxies traces the large-scale structure of the Universe and therefore constitutes a prime observable to constrain cosmological parameters. We conduct Bayesian inference of the $Λ$CDM parameters $Ω_m$, $Ω_b$, $h$, $n_s$, and $σ_8$ from the BOSS CMASS galaxy sample by combining the wavelet scattering transform (WST) with a simulation-based inference approa…
▽ More
The non-Gaussisan spatial distribution of galaxies traces the large-scale structure of the Universe and therefore constitutes a prime observable to constrain cosmological parameters. We conduct Bayesian inference of the $Λ$CDM parameters $Ω_m$, $Ω_b$, $h$, $n_s$, and $σ_8$ from the BOSS CMASS galaxy sample by combining the wavelet scattering transform (WST) with a simulation-based inference approach enabled by the ${\rm S{\scriptsize IM}BIG}$ forward model. We design a set of reduced WST statistics that leverage symmetries of redshift-space data. Posterior distributions are estimated with a conditional normalizing flow trained on 20,000 simulated ${\rm S{\scriptsize IM}BIG}$ galaxy catalogs with survey realism. We assess the accuracy of the posterior estimates using simulation-based calibration and quantify generalization and robustness to the change of forward model using a suite of 2,000 test simulations. When probing scales down to $k_{\rm max}=0.5~h/\text{Mpc}$, we are able to derive accurate posterior estimates that are robust to the change of forward model for all parameters, except $σ_8$. We mitigate the robustness issues with $σ_8$ by removing the WST coefficients that probe scales smaller than $k \sim 0.3~h/\text{Mpc}$. Applied to the BOSS CMASS sample, our WST analysis yields seemingly improved constraints obtained from a standard PT-based power spectrum analysis with $k_{\rm max}=0.25~h/\text{Mpc}$ for all parameters except $h$. However, we still raise concerns on these results. The observational predictions significantly vary across different normalizing flow architectures, which we interpret as a form of model misspecification. This highlights a key challenge for forward modeling approaches when using summary statistics that are sensitive to detailed model-specific or observational imprints on galaxy clustering.
△ Less
Submitted 18 July, 2024; v1 submitted 23 October, 2023;
originally announced October 2023.
-
${\rm S{\scriptsize IM}BIG}$: The First Cosmological Constraints from Non-Gaussian and Non-Linear Galaxy Clustering
Authors:
ChangHoon Hahn,
Pablo Lemos,
Liam Parker,
Bruno Régaldo-Saint Blancard,
Michael Eickenberg,
Shirley Ho,
Jiamin Hou,
Elena Massara,
Chirag Modi,
Azadeh Moradinezhad Dizgah,
David Spergel
Abstract:
The 3D distribution of galaxies encodes detailed cosmological information on the expansion and growth history of the Universe. We present the first cosmological constraints that exploit non-Gaussian cosmological information on non-linear scales from galaxy clustering, inaccessible with current standard analyses. We analyze a subset of the BOSS galaxy survey using ${\rm S{\scriptsize IM}BIG}$, a ne…
▽ More
The 3D distribution of galaxies encodes detailed cosmological information on the expansion and growth history of the Universe. We present the first cosmological constraints that exploit non-Gaussian cosmological information on non-linear scales from galaxy clustering, inaccessible with current standard analyses. We analyze a subset of the BOSS galaxy survey using ${\rm S{\scriptsize IM}BIG}$, a new framework for cosmological inference that leverages high-fidelity simulations and deep generative models. We use two clustering statistics beyond the standard power spectrum: the bispectrum and a convolutional neural network based summary of the galaxy field. We infer constraints on $Λ$CDM parameters, $Ω_b$, $h$, $n_s$, $Ω_m$, and $σ_8$, that are 1.6, 1.5, 1.7, 1.2, and 2.3$\times$ tighter than power spectrum analyses. With this increased precision, we derive constraints on the Hubble constant, $H_0$, and $S_8 = σ_8 \sqrt{Ω_m/0.3}$ that are competitive with other cosmological probes, even with a sample that only spans 10% of the full BOSS volume. Our $H_0$ constraints, imposing the Big Bang Nucleosynthesis prior on the baryon density, are consistent with the early time constraints from the cosmic microwave background (CMB). Meanwhile, our $S_8$ constraints are consistent with weak lensing experiments and similarly lie below CMB constraints. Lastly, we present forecasts to show that future work extending ${\rm S{\scriptsize IM}BIG}$ to upcoming spectroscopic galaxy surveys (DESI, PFS, Euclid) will produce leading $H_0$ and $S_8$ constraints that bridge the gap between early and late time measurements and shed light on current cosmic tensions.
△ Less
Submitted 23 October, 2023;
originally announced October 2023.
-
${\rm S{\scriptsize IM}BIG}$: The First Cosmological Constraints from the Non-Linear Galaxy Bispectrum
Authors:
ChangHoon Hahn,
Michael Eickenberg,
Shirley Ho,
Jiamin Hou,
Pablo Lemos,
Elena Massara,
Chirag Modi,
Azadeh Moradinezhad Dizgah,
Liam Parker,
Bruno Régaldo-Saint Blancard
Abstract:
We present the first cosmological constraints from analyzing higher-order galaxy clustering on non-linear scales. We use ${\rm S{\scriptsize IM}BIG}$, a forward modeling framework for galaxy clustering analyses that employs simulation-based inference to perform highly efficient cosmological inference using normalizing flows. It leverages the predictive power of high-fidelity simulations and robust…
▽ More
We present the first cosmological constraints from analyzing higher-order galaxy clustering on non-linear scales. We use ${\rm S{\scriptsize IM}BIG}$, a forward modeling framework for galaxy clustering analyses that employs simulation-based inference to perform highly efficient cosmological inference using normalizing flows. It leverages the predictive power of high-fidelity simulations and robustly extracts cosmological information from regimes inaccessible with current standard analyses. In this work, we apply ${\rm S{\scriptsize IM}BIG}$ to a subset of the BOSS galaxy sample and analyze the redshift-space bispectrum monopole, $B_0(k_1, k_2, k_3)$, to $k_{\rm max}=0.5\,h/{\rm Mpc}$. We achieve 1$σ$ constraints of $Ω_m=0.293^{+0.027}_{-0.027}$ and $σ_8= 0.783^{+0.040}_{-0.038}$, which are more than 1.2 and 2.4$\times$ tighter than constraints from standard power spectrum analyses of the same dataset. We also derive 1.4, 1.4, 1.7$\times$ tighter constraints on $Ω_b$, $h$, $n_s$. This improvement comes from additional cosmological information in higher-order clustering on non-linear scales and, for $σ_8$, is equivalent to the gain expected from a standard analysis on a $\sim$4$\times$ larger galaxy sample. Even with our BOSS subsample, which only spans 10% of the full BOSS volume, we derive competitive constraints on the growth of structure: $S_8 = 0.774^{+0.056}_{-0.053}$. Our constraint is consistent with results from both cosmic microwave background and weak lensing. Combined with a $ω_b$ prior from Big Bang Nucleosynthesis, we also derive a constraint on $H_0=67.6^{+2.2}_{-1.8}\,{\rm km\,s^{-1}\,Mpc^{-1}}$ that is consistent with early universe constraints.
△ Less
Submitted 23 October, 2023;
originally announced October 2023.
-
Light Thermal Self-Interacting Dark Matter in the Shadow of Non-Standard Cosmology
Authors:
Shu-Yu Ho,
Pyungwon Ko,
Dibyendu Nanda
Abstract:
In this paper, we construct a viable model for a GeV scale self-interacting dark matter (DM), where the DM was thermally produced in the early universe. Here, a new vector-like fermion with a dark charge under the $U(1)_{D}$ gauge symmetry serves as a secluded WIMP DM and it can dominantly annihilate into the light dark gauge boson and singlet scalar through the dark gauge interaction. Also, the s…
▽ More
In this paper, we construct a viable model for a GeV scale self-interacting dark matter (DM), where the DM was thermally produced in the early universe. Here, a new vector-like fermion with a dark charge under the $U(1)_{D}$ gauge symmetry serves as a secluded WIMP DM and it can dominantly annihilate into the light dark gauge boson and singlet scalar through the dark gauge interaction. Also, the self-interaction of DM is induced by the light dark gauge boson via the same gauge interaction. In addition to these particles, we further introduce two Weyl fermions and a doublet scalar, by which the dark gauge boson produced from s-wave DM annihilations can mostly decay into active neutrinos after the dark symmetry breaking such that the CMB bound on the DM with low masses can be eluded. In order to have a common parameter region to explain the observed relic abundance and self-interaction of DM, we also study this model in a non-standard cosmological evolution, where the cosmic expansion driven by a new field species is faster than the standard radiation-dominated universe during the frozen time of DM. Reversely, one can also use the self-interacting nature of light thermal DM to examine the non-standard cosmological history of the universe.
△ Less
Submitted 25 March, 2024; v1 submitted 9 October, 2023;
originally announced October 2023.
-
AstroCLIP: A Cross-Modal Foundation Model for Galaxies
Authors:
Liam Parker,
Francois Lanusse,
Siavash Golkar,
Leopoldo Sarra,
Miles Cranmer,
Alberto Bietti,
Michael Eickenberg,
Geraud Krawezik,
Michael McCabe,
Ruben Ohana,
Mariel Pettee,
Bruno Regaldo-Saint Blancard,
Tiberiu Tesileanu,
Kyunghyun Cho,
Shirley Ho
Abstract:
We present AstroCLIP, a single, versatile model that can embed both galaxy images and spectra into a shared, physically meaningful latent space. These embeddings can then be used - without any model fine-tuning - for a variety of downstream tasks including (1) accurate in-modality and cross-modality semantic similarity search, (2) photometric redshift estimation, (3) galaxy property estimation fro…
▽ More
We present AstroCLIP, a single, versatile model that can embed both galaxy images and spectra into a shared, physically meaningful latent space. These embeddings can then be used - without any model fine-tuning - for a variety of downstream tasks including (1) accurate in-modality and cross-modality semantic similarity search, (2) photometric redshift estimation, (3) galaxy property estimation from both images and spectra, and (4) morphology classification. Our approach to implementing AstroCLIP consists of two parts. First, we embed galaxy images and spectra separately by pretraining separate transformer-based image and spectrum encoders in self-supervised settings. We then align the encoders using a contrastive loss. We apply our method to spectra from the Dark Energy Spectroscopic Instrument and images from its corresponding Legacy Imaging Survey. Overall, we find remarkable performance on all downstream tasks, even relative to supervised baselines. For example, for a task like photometric redshift prediction, we find similar performance to a specifically-trained ResNet18, and for additional tasks like physical property estimation (stellar mass, age, metallicity, and sSFR), we beat this supervised baseline by 19\% in terms of $R^2$. We also compare our results to a state-of-the-art self-supervised single-modal model for galaxy images, and find that our approach outperforms this benchmark by roughly a factor of two on photometric redshift estimation and physical property prediction in terms of $R^2$, while remaining roughly in-line in terms of morphology classification. Ultimately, our approach represents the first cross-modal self-supervised model for galaxies, and the first self-supervised transformer-based architectures for galaxy images and spectra.
△ Less
Submitted 14 June, 2024; v1 submitted 4 October, 2023;
originally announced October 2023.
-
Learnable wavelet neural networks for cosmological inference
Authors:
Christian Pedersen,
Michael Eickenberg,
Shirley Ho
Abstract:
Convolutional neural networks (CNNs) have been shown to both extract more information than the traditional two-point statistics from cosmological fields, and marginalise over astrophysical effects extremely well. However, CNNs require large amounts of training data, which is potentially problematic in the domain of expensive cosmological simulations, and it is difficult to interpret the network. I…
▽ More
Convolutional neural networks (CNNs) have been shown to both extract more information than the traditional two-point statistics from cosmological fields, and marginalise over astrophysical effects extremely well. However, CNNs require large amounts of training data, which is potentially problematic in the domain of expensive cosmological simulations, and it is difficult to interpret the network. In this work we apply the learnable scattering transform, a kind of convolutional neural network that uses trainable wavelets as filters, to the problem of cosmological inference and marginalisation over astrophysical effects. We present two models based on the scattering transform, one constructed for performance, and one constructed for interpretability, and perform a comparison with a CNN. We find that scattering architectures are able to outperform a CNN, significantly in the case of small training data samples. Additionally we present a lightweight scattering network that is highly interpretable.
△ Less
Submitted 24 July, 2023;
originally announced July 2023.
-
Neutrino mass constraint from an Implicit Likelihood Analysis of BOSS voids
Authors:
Leander Thiele,
Elena Massara,
Alice Pisani,
ChangHoon Hahn,
David N. Spergel,
Shirley Ho,
Benjamin Wandelt
Abstract:
Cosmic voids identified in the spatial distribution of galaxies provide complementary information to two-point statistics. In particular, constraints on the neutrino mass sum, $\sum m_ν$, promise to benefit from the inclusion of void statistics. We perform inference on the CMASS NGC sample of SDSS-III/BOSS with the aim of constraining $\sum m_ν$. We utilize the void size function, the void galaxy…
▽ More
Cosmic voids identified in the spatial distribution of galaxies provide complementary information to two-point statistics. In particular, constraints on the neutrino mass sum, $\sum m_ν$, promise to benefit from the inclusion of void statistics. We perform inference on the CMASS NGC sample of SDSS-III/BOSS with the aim of constraining $\sum m_ν$. We utilize the void size function, the void galaxy cross power spectrum, and the galaxy auto power spectrum. To extract constraints from these summary statistics we use a simulation-based approach, specifically implicit likelihood inference. We populate approximate gravity-only, particle neutrino cosmological simulations with an expressive halo occupation distribution model. With a conservative scale cut of kmax=0.15 h/Mpc and a Planck-inspired LCDM prior, we find upper bounds on $\sum m_ν$ of 0.43 and 0.35 eV from the galaxy auto power spectrum and the full data vector, respectively (95% credible interval). We observe hints that the void statistics may be most effective at constraining $\sum m_ν$ from below. We also substantiate the usual assumption that the void size function is Poisson distributed.
△ Less
Submitted 14 July, 2023;
originally announced July 2023.
-
Machine Learning Classification of Repeating FRBs from FRB121102
Authors:
Bjorn Jasper R. Raquel,
Tetsuya Hashimoto,
Tomotsugu Goto,
Bo Han Chen,
Yuri Uno,
Tiger Yu-Yang Hsiao,
Seong Jin Kim,
Simon C. -C. Ho
Abstract:
Fast Radio Bursts (FRBs) are mysterious bursts in the millisecond timescale at radio wavelengths. Currently, there is little understanding about the classification of repeating FRBs, based on difference in physics, which is of great importance in understanding their origin. Recent works from the literature focus on using specific parameters to classify FRBs to draw inferences on the possible physi…
▽ More
Fast Radio Bursts (FRBs) are mysterious bursts in the millisecond timescale at radio wavelengths. Currently, there is little understanding about the classification of repeating FRBs, based on difference in physics, which is of great importance in understanding their origin. Recent works from the literature focus on using specific parameters to classify FRBs to draw inferences on the possible physical mechanisms or properties of these FRB subtypes. In this study, we use publicly available 1652 repeating FRBs from FRB121102 detected with the Five-hundred-meter Aperture Spherical Telescope (FAST), and studied them with an unsupervised machine learning model. By fine-tuning the hyperparameters of the model, we found that there is an indication for four clusters from the bursts of FRB121102 instead of the two clusters ("Classical" and "Atypical") suggested in the literature. Wherein, the "Atypical" cluster can be further classified into three sub-clusters with distinct characteristics. Our findings show that the clustering result we obtained is more comprehensive not only because our study produced results which are consistent with those in the literature but also because our work uses more physical parameters to create these clusters. Overall, our methods and analyses produced a more holistic approach in clustering the repeating FRBs of FRB121102.
△ Less
Submitted 6 July, 2023; v1 submitted 6 July, 2023;
originally announced July 2023.
-
The Atacama Cosmology Telescope: High-resolution component-separated maps across one-third of the sky
Authors:
William R. Coulton,
Mathew S. Madhavacheril,
Adriaan J. Duivenvoorden,
J. Colin Hill,
Irene Abril-Cabezas,
Peter A. R. Ade,
Simone Aiola,
Tommy Alford,
Mandana Amiri,
Stefania Amodeo,
Rui An,
Zachary Atkins,
Jason E. Austermann,
Nicholas Battaglia,
Elia Stefano Battistelli,
James A. Beall,
Rachel Bean,
Benjamin Beringue,
Tanay Bhandarkar,
Emily Biermann,
Boris Bolliet,
J Richard Bond,
Hongbo Cai,
Erminia Calabrese,
Victoria Calafut
, et al. (129 additional authors not shown)
Abstract:
Observations of the millimeter sky contain valuable information on a number of signals, including the blackbody cosmic microwave background (CMB), Galactic emissions, and the Compton-$y$ distortion due to the thermal Sunyaev-Zel'dovich (tSZ) effect. Extracting new insight into cosmological and astrophysical questions often requires combining multi-wavelength observations to spectrally isolate one…
▽ More
Observations of the millimeter sky contain valuable information on a number of signals, including the blackbody cosmic microwave background (CMB), Galactic emissions, and the Compton-$y$ distortion due to the thermal Sunyaev-Zel'dovich (tSZ) effect. Extracting new insight into cosmological and astrophysical questions often requires combining multi-wavelength observations to spectrally isolate one component. In this work, we present a new arcminute-resolution Compton-$y$ map, which traces out the line-of-sight-integrated electron pressure, as well as maps of the CMB in intensity and E-mode polarization, across a third of the sky (around 13,000 sq.~deg.). We produce these through a joint analysis of data from the Atacama Cosmology Telescope (ACT) Data Release 4 and 6 at frequencies of roughly 93, 148, and 225 GHz, together with data from the \textit{Planck} satellite at frequencies between 30 GHz and 545 GHz. We present detailed verification of an internal linear combination pipeline implemented in a needlet frame that allows us to efficiently suppress Galactic contamination and account for spatial variations in the ACT instrument noise. These maps provide a significant advance, in noise levels and resolution, over the existing \textit{Planck} component-separated maps and will enable a host of science goals including studies of cluster and galaxy astrophysics, inferences of the cosmic velocity field, primordial non-Gaussianity searches, and gravitational lensing reconstruction of the CMB.
△ Less
Submitted 3 July, 2023;
originally announced July 2023.
-
A T-Dwarf Candidate from JWST Early Release NIRCam data
Authors:
Po-Ya Wang,
Tomotsugu Goto,
Simon C. -C. Ho,
Yu-Wei Lin,
Cossas K. -W. Wu,
Chih-Teng Ling,
Tetsuya Hashimoto,
Seong Jin Kim,
Tiger Y. -Y. Hsiao
Abstract:
We present a distant T$-$type brown dwarf candidate at $\approx2.55$ kpc discovered in the Cosmic Evolution Early Release Science (CEERS) fields by James Webb Space Telescope (JWST) NIRCam. In addition to the superb sensitivity, we utilised 7 filters from JWST in near-IR and thus is advantageous in finding faint, previously unseen brown dwarfs. From the model spectra in new JWST/NIRCam filter wave…
▽ More
We present a distant T$-$type brown dwarf candidate at $\approx2.55$ kpc discovered in the Cosmic Evolution Early Release Science (CEERS) fields by James Webb Space Telescope (JWST) NIRCam. In addition to the superb sensitivity, we utilised 7 filters from JWST in near-IR and thus is advantageous in finding faint, previously unseen brown dwarfs. From the model spectra in new JWST/NIRCam filter wavelengths, the selection criteria of F115W-F277W$<$-0.8 and F277W-F444W$>$1.1 were chosen to target the spectrum features of brown dwarfs having temperatures from 500K to 1300K. Searching through the data from Early Release Observations (ERO) and Early Release Science (ERS), we find 1 promising candidate in the CEERS field. The result of SED fitting suggested an early T spectral type with a low effective temperature of T$_\text{eff}\approx$1300K, the surface gravity of $\log{g}\approx5.25\text{cm s}^{-2}$, and an eddy diffusion parameter of logK$_{zz}\approx7\text{cm}^2 \text{s}^{-1}$, which indicates an age of $\approx$1.8Gyr and a mass of $\approx0.05$M$_{\odot}$. In contrast to typically found T$-$dwarf within several hundred parsecs, the estimated distance of the source is $\approx2.55$kpc, showing the JWST's power to extend the search to a much larger distance.
△ Less
Submitted 5 June, 2023;
originally announced June 2023.
-
The Atacama Cosmology Telescope: DR6 Gravitational Lensing Map and Cosmological Parameters
Authors:
Mathew S. Madhavacheril,
Frank J. Qu,
Blake D. Sherwin,
Niall MacCrann,
Yaqiong Li,
Irene Abril-Cabezas,
Peter A. R. Ade,
Simone Aiola,
Tommy Alford,
Mandana Amiri,
Stefania Amodeo,
Rui An,
Zachary Atkins,
Jason E. Austermann,
Nicholas Battaglia,
Elia Stefano Battistelli,
James A. Beall,
Rachel Bean,
Benjamin Beringue,
Tanay Bhandarkar,
Emily Biermann,
Boris Bolliet,
J Richard Bond,
Hongbo Cai,
Erminia Calabrese
, et al. (134 additional authors not shown)
Abstract:
We present cosmological constraints from a gravitational lensing mass map covering 9400 sq. deg. reconstructed from CMB measurements made by the Atacama Cosmology Telescope (ACT) from 2017 to 2021. In combination with BAO measurements (from SDSS and 6dF), we obtain the amplitude of matter fluctuations $σ_8 = 0.819 \pm 0.015$ at 1.8% precision, $S_8\equivσ_8({Ω_{\rm m}}/0.3)^{0.5}=0.840\pm0.028$ an…
▽ More
We present cosmological constraints from a gravitational lensing mass map covering 9400 sq. deg. reconstructed from CMB measurements made by the Atacama Cosmology Telescope (ACT) from 2017 to 2021. In combination with BAO measurements (from SDSS and 6dF), we obtain the amplitude of matter fluctuations $σ_8 = 0.819 \pm 0.015$ at 1.8% precision, $S_8\equivσ_8({Ω_{\rm m}}/0.3)^{0.5}=0.840\pm0.028$ and the Hubble constant $H_0= (68.3 \pm 1.1)\, \text{km}\,\text{s}^{-1}\,\text{Mpc}^{-1}$ at 1.6% precision. A joint constraint with CMB lensing measured by the Planck satellite yields even more precise values: $σ_8 = 0.812 \pm 0.013$, $S_8\equivσ_8({Ω_{\rm m}}/0.3)^{0.5}=0.831\pm0.023$ and $H_0= (68.1 \pm 1.0)\, \text{km}\,\text{s}^{-1}\,\text{Mpc}^{-1}$. These measurements agree well with $Λ$CDM-model extrapolations from the CMB anisotropies measured by Planck. To compare these constraints to those from the KiDS, DES, and HSC galaxy surveys, we revisit those data sets with a uniform set of assumptions, and find $S_8$ from all three surveys are lower than that from ACT+Planck lensing by varying levels ranging from 1.7-2.1$σ$. These results motivate further measurements and comparison, not just between the CMB anisotropies and galaxy lensing, but also between CMB lensing probing $z\sim 0.5-5$ on mostly-linear scales and galaxy lensing at $z\sim 0.5$ on smaller scales. We combine our CMB lensing measurements with CMB anisotropies to constrain extensions of $Λ$CDM, limiting the sum of the neutrino masses to $\sum m_ν < 0.13$ eV (95% c.l.), for example. Our results provide independent confirmation that the universe is spatially flat, conforms with general relativity, and is described remarkably well by the $Λ$CDM model, while paving a promising path for neutrino physics with gravitational lensing from upcoming ground-based CMB surveys.
△ Less
Submitted 12 August, 2024; v1 submitted 11 April, 2023;
originally announced April 2023.
-
The Atacama Cosmology Telescope: A Measurement of the DR6 CMB Lensing Power Spectrum and its Implications for Structure Growth
Authors:
Frank J. Qu,
Blake D. Sherwin,
Mathew S. Madhavacheril,
Dongwon Han,
Kevin T. Crowley,
Irene Abril-Cabezas,
Peter A. R. Ade,
Simone Aiola,
Tommy Alford,
Mandana Amiri,
Stefania Amodeo,
Rui An,
Zachary Atkins,
Jason E. Austermann,
Nicholas Battaglia,
Elia Stefano Battistelli,
James A. Beall,
Rachel Bean,
Benjamin Beringue,
Tanay Bhandarkar,
Emily Biermann,
Boris Bolliet,
J Richard Bond,
Hongbo Cai,
Erminia Calabrese
, et al. (133 additional authors not shown)
Abstract:
We present new measurements of cosmic microwave background (CMB) lensing over $9400$ sq. deg. of the sky. These lensing measurements are derived from the Atacama Cosmology Telescope (ACT) Data Release 6 (DR6) CMB dataset, which consists of five seasons of ACT CMB temperature and polarization observations. We determine the amplitude of the CMB lensing power spectrum at $2.3\%$ precision ($43σ$ sign…
▽ More
We present new measurements of cosmic microwave background (CMB) lensing over $9400$ sq. deg. of the sky. These lensing measurements are derived from the Atacama Cosmology Telescope (ACT) Data Release 6 (DR6) CMB dataset, which consists of five seasons of ACT CMB temperature and polarization observations. We determine the amplitude of the CMB lensing power spectrum at $2.3\%$ precision ($43σ$ significance) using a novel pipeline that minimizes sensitivity to foregrounds and to noise properties. To ensure our results are robust, we analyze an extensive set of null tests, consistency tests, and systematic error estimates and employ a blinded analysis framework. The baseline spectrum is well fit by a lensing amplitude of $A_{\mathrm{lens}}=1.013\pm0.023$ relative to the Planck 2018 CMB power spectra best-fit $Λ$CDM model and $A_{\mathrm{lens}}=1.005\pm0.023$ relative to the $\text{ACT DR4} + \text{WMAP}$ best-fit model. From our lensing power spectrum measurement, we derive constraints on the parameter combination $S^{\mathrm{CMBL}}_8 \equiv σ_8 \left({Ω_m}/{0.3}\right)^{0.25}$ of $S^{\mathrm{CMBL}}_8= 0.818\pm0.022$ from ACT DR6 CMB lensing alone and $S^{\mathrm{CMBL}}_8= 0.813\pm0.018$ when combining ACT DR6 and Planck NPIPE CMB lensing power spectra. These results are in excellent agreement with $Λ$CDM model constraints from Planck or $\text{ACT DR4} + \text{WMAP}$ CMB power spectrum measurements. Our lensing measurements from redshifts $z\sim0.5$--$5$ are thus fully consistent with $Λ$CDM structure growth predictions based on CMB anisotropies probing primarily $z\sim1100$. We find no evidence for a suppression of the amplitude of cosmic structure at low redshifts
△ Less
Submitted 28 May, 2024; v1 submitted 11 April, 2023;
originally announced April 2023.
-
Future Constraints on Dark Matter with Gravitationally Lensed Fast Radio Bursts Detected by BURSTT
Authors:
Simon C. -C. Ho,
Tetsuya Hashimoto,
Tomotsugu Goto,
Yu-Wei Lin,
Seong Jin Kim,
Yuri Uno,
Tiger Y. -Y. Hsiao
Abstract:
Understanding dark matter is one of the most urgent questions in modern physics. A very interesting candidate is primordial black holes (PBHs; Carr2016). For the mass ranges of $< 10^{-16} M_{\odot}$ and $> 100 M_{\odot}$, PBHs have been ruled out. However, they are still poorly constrained in the mass ranges of $10^{-16} - 100 M_{\odot}$ (Belotsky et al. 2019). Fast radio bursts (FRBs) are millis…
▽ More
Understanding dark matter is one of the most urgent questions in modern physics. A very interesting candidate is primordial black holes (PBHs; Carr2016). For the mass ranges of $< 10^{-16} M_{\odot}$ and $> 100 M_{\odot}$, PBHs have been ruled out. However, they are still poorly constrained in the mass ranges of $10^{-16} - 100 M_{\odot}$ (Belotsky et al. 2019). Fast radio bursts (FRBs) are millisecond flashes of radio light of unknown origin mostly from outside the Milky Way. Due to their short timescales, gravitationally lensed FRBs, which are yet to be detected, have been proposed as a useful probe for constraining the presence of PBHs in the mass window of $< 100M_{\odot}$ (Mu�oz et al. 2016). Up to now, the most successful project in finding FRBs has been CHIME. Due to its large field of view (FoV), CHIME is detecting at least 600 FRBs since 2018. However, none of them is confirmed to be gravitationally lensed (Leung et al. 2022). Taiwan plans to build a new telescope, BURSTT dedicated to detecting FRBs. Its survey area will be 25 times greater than CHIME. BURSTT can localize all of these FRBs through very-long-baseline interferometry (VLBI). We estimate the probability to find gravitationally lensed FRBs, based on the scaled redshift distribution from the latest CHIME catalog and the lensing probability function from Mu�oz et al. (2016). BURSTT-2048 can detect ~ 24 lensed FRBs out of ~ 1,700 FRBs per annum. With BURSTT's ability to detect nanosecond FRBs, we can constrain PBHs to form a part of dark matter down to $10^{-4}M_{\odot}$.
△ Less
Submitted 11 April, 2023;
originally announced April 2023.
-
Latent Stochastic Differential Equations for Modeling Quasar Variability and Inferring Black Hole Properties
Authors:
Joshua Fagin,
Ji Won Park,
Henry Best,
James Hung-Hsu Chan,
K. E Saavik Ford,
Matthew J. Graham,
V. Ashley Villar,
Shirley Ho,
Matthew O'Dowd
Abstract:
Quasars are bright and unobscured active galactic nuclei (AGN) thought to be powered by the accretion of matter around supermassive black holes at the centers of galaxies. The temporal variability of a quasar's brightness contains valuable information about its physical properties. The UV/optical variability is thought to be a stochastic process, often represented as a damped random walk described…
▽ More
Quasars are bright and unobscured active galactic nuclei (AGN) thought to be powered by the accretion of matter around supermassive black holes at the centers of galaxies. The temporal variability of a quasar's brightness contains valuable information about its physical properties. The UV/optical variability is thought to be a stochastic process, often represented as a damped random walk described by a stochastic differential equation (SDE). Upcoming wide-field telescopes such as the Rubin Observatory Legacy Survey of Space and Time (LSST) are expected to observe tens of millions of AGN in multiple filters over a ten year period, so there is a need for efficient and automated modeling techniques that can handle the large volume of data. Latent SDEs are machine learning models well suited for modeling quasar variability, as they can explicitly capture the underlying stochastic dynamics. In this work, we adapt latent SDEs to jointly reconstruct multivariate quasar light curves and infer their physical properties such as the black hole mass, inclination angle, and temperature slope. Our model is trained on realistic simulations of LSST ten year quasar light curves, and we demonstrate its ability to reconstruct quasar light curves even in the presence of long seasonal gaps and irregular sampling across different bands, outperforming a multioutput Gaussian process regression baseline. Our method has the potential to provide a deeper understanding of the physical properties of quasars and is applicable to a wide range of other multivariate time series with missing data and irregular sampling.
△ Less
Submitted 16 April, 2024; v1 submitted 9 April, 2023;
originally announced April 2023.
-
Posterior Sampling of the Initial Conditions of the Universe from Non-linear Large Scale Structures using Score-Based Generative Models
Authors:
Ronan Legin,
Matthew Ho,
Pablo Lemos,
Laurence Perreault-Levasseur,
Shirley Ho,
Yashar Hezaveh,
Benjamin Wandelt
Abstract:
Reconstructing the initial conditions of the universe is a key problem in cosmology. Methods based on simulating the forward evolution of the universe have provided a way to infer initial conditions consistent with present-day observations. However, due to the high complexity of the inference problem, these methods either fail to sample a distribution of possible initial density fields or require…
▽ More
Reconstructing the initial conditions of the universe is a key problem in cosmology. Methods based on simulating the forward evolution of the universe have provided a way to infer initial conditions consistent with present-day observations. However, due to the high complexity of the inference problem, these methods either fail to sample a distribution of possible initial density fields or require significant approximations in the simulation model to be tractable, potentially leading to biased results. In this work, we propose the use of score-based generative models to sample realizations of the early universe given present-day observations. We infer the initial density field of full high-resolution dark matter N-body simulations from the present-day density field and verify the quality of produced samples compared to the ground truth based on summary statistics. The proposed method is capable of providing plausible realizations of the early universe density field from the initial conditions posterior distribution marginalized over cosmological parameters and can sample orders of magnitude faster than current state-of-the-art methods.
△ Less
Submitted 7 April, 2023;
originally announced April 2023.
-
Upper limits on transmitter rate of extragalactic civilizations placed by Breakthrough Listen observations
Authors:
Yuri Uno,
Tetsuya Hashimoto,
Tomotsugu Goto,
Simon C. -C. Ho,
Tzu-Yin Hsu,
Ross Burns
Abstract:
The Search for Extra-Terrestrial Intelligence (SETI) has been conducted for over sixty years, yet no technosignatures have been identified. Previous studies have focused on stars in our galaxy, with few searches in the extragalactic Universe despite a larger volume being available. Civilizations capable of harvesting energy from a star or a galaxy are classified as KII or KIII on the Kardashev sca…
▽ More
The Search for Extra-Terrestrial Intelligence (SETI) has been conducted for over sixty years, yet no technosignatures have been identified. Previous studies have focused on stars in our galaxy, with few searches in the extragalactic Universe despite a larger volume being available. Civilizations capable of harvesting energy from a star or a galaxy are classified as KII or KIII on the Kardashev scale, respectively. Technosignatures from such advanced civilizations would be extremely luminous and detectable by current radio telescopes, even from distant galaxies. To explore the frontier of extragalactic SETI, we investigate the likely prevalence of extragalactic civilizations possessing a radio transmitter, known as the transmitter rate, based on observational results from the Breakthrough Listen (BL) observations. We calculated the transmitter rate by considering the background galaxies in the field of view of target stars in BL observations. We used a statistical method to derive the total mass of stars in those background galaxies from a galaxy stellar mass function. Our statistical method suggests that less than one in hundreds of trillions of extragalactic civilizations within 969 Mpc possess a radio transmitter above 7.7$\times$10$^{26}$ W of power, assuming one civilization per one-solar-mass stellar system. Additionally, we cross-matched the BL survey fields with the WISE$\times$SuperCOSMOS Photometric Redshift Catalogue and compared with the statistical method. Our result sets the strictest limits to date on the transmitter rate at such high power levels, emphasizing the high efficiency of searching for radio transmitters in galaxies and the rarity of technologically advanced civilizations in our Universe.
△ Less
Submitted 8 April, 2023; v1 submitted 5 April, 2023;
originally announced April 2023.
-
Classifying a frequently repeating fast radio burst, FRB 20201124A, with unsupervised machine learning
Authors:
Bo Han Chen,
Tetsuya Hashimoto,
Tomotsugu Goto,
Bjorn Jasper R. Raquel,
Yuri Uno,
Seong Jin Kim,
Tiger Y. -Y. Hsiao,
Simon C. -C. Ho
Abstract:
Fast radio bursts (FRBs) are astronomical transients with millisecond timescales. Although most of the FRBs are not observed to repeat, a few of them are detected to repeat more than hundreds of times. There exist a large variety of physical properties among these bursts, suggesting heterogeneous mechanisms of FRBs. In this paper, we conduct a categorisation on the extremely frequently repeating F…
▽ More
Fast radio bursts (FRBs) are astronomical transients with millisecond timescales. Although most of the FRBs are not observed to repeat, a few of them are detected to repeat more than hundreds of times. There exist a large variety of physical properties among these bursts, suggesting heterogeneous mechanisms of FRBs. In this paper, we conduct a categorisation on the extremely frequently repeating FRB 20201124A with the assistance of machine learning, as such techniques have the potential to use subtle differences and correlations that humans are unaware of to better classify bursts. The research is carried out by applying the unsupervised Uniform Manifold Approximation and Projection (UMAP) model on the FRB 20201124A data provided by Five-hundred-meter Aperture Spherical radio Telescope (FAST). The algorithm eventually categorises the bursts into three clusters. In addition to the two categories in previous work based on waiting time, a new way for categorisation has been found. The three clusters are either high energy, high frequency, or low frequency, reflecting the distribution of FRB energy and frequency. Importantly, a similar machine learning result is found in another frequently repeating FRB20121102A, implying a common mechanism among this kind of FRB. This work is one of the first steps towards the systematical categorisation of the extremely frequently repeating FRBs.
△ Less
Submitted 31 March, 2023; v1 submitted 29 March, 2023;
originally announced March 2023.
-
Predicting the Initial Conditions of the Universe using a Deterministic Neural Network
Authors:
Vaibhav Jindal,
Albert Liang,
Aarti Singh,
Shirley Ho,
Drew Jamieson
Abstract:
Finding the initial conditions that led to the current state of the universe is challenging because it involves searching over an intractable input space of initial conditions, along with modeling their evolution via tools such as N-body simulations which are computationally expensive. Recently, deep learning has emerged as a surrogate for N-body simulations by directly learning the mapping betwee…
▽ More
Finding the initial conditions that led to the current state of the universe is challenging because it involves searching over an intractable input space of initial conditions, along with modeling their evolution via tools such as N-body simulations which are computationally expensive. Recently, deep learning has emerged as a surrogate for N-body simulations by directly learning the mapping between the linear input of an N-body simulation and the final nonlinear output from the simulation, significantly accelerating the forward modeling. However, this still does not reduce the search space for initial conditions. In this work, we pioneer the use of a deterministic convolutional neural network for learning the reverse mapping and show that it accurately recovers the initial linear displacement field over a wide range of scales ($<1$-$2\%$ error up to nearly $k\simeq0.8$-$0.9 \text{ Mpc}^{-1}h$), despite the one-to-many mapping of the inverse problem (due to the divergent backward trajectories at smaller scales). Specifically, we train a V-Net architecture, which outputs the linear displacement of an N-body simulation, given the nonlinear displacement at redshift $z=0$ and the cosmological parameters. The results of our method suggest that a simple deterministic neural network is sufficient for accurately approximating the initial linear states, potentially obviating the need for the more complex and computationally demanding backward modeling methods that were recently proposed.
△ Less
Submitted 13 December, 2023; v1 submitted 23 March, 2023;
originally announced March 2023.
-
Infrared galaxies detected by the Atacama Cosmology Telescope
Authors:
Ece Kilerci,
Tetsuya Hashimoto,
Tomotsugu Goto,
Ersin Göğüş,
Seong Jin Kim,
Simon C. -C. Ho,
Yi Hang Valerie Wong
Abstract:
We report on 167 infrared (IR) galaxies selected by AKARI and IRAS and detected in the Atacama Cosmology Telescope (ACT) Data Release 5 (DR5) sky maps at the 98, 150 and 220 GHz frequency bands. Of these detections, 134 (80%) of the millimeter counterparts are first-time identifications with ACT. We expand the previous ACT extragalactic source catalogs, by including new 98 GHz detections and measu…
▽ More
We report on 167 infrared (IR) galaxies selected by AKARI and IRAS and detected in the Atacama Cosmology Telescope (ACT) Data Release 5 (DR5) sky maps at the 98, 150 and 220 GHz frequency bands. Of these detections, 134 (80%) of the millimeter counterparts are first-time identifications with ACT. We expand the previous ACT extragalactic source catalogs, by including new 98 GHz detections and measurements from ACT DR5. We also report flux density measurements at the 98, 150, and 220 GHz frequency bands. We compute $α_{98-150}$, $α_{98-220}$, and $α_{150-220}$ millimeter-wave spectral indices and far-IR to millimeter-wave spectral indices between 90 micron and 98, 150, and 220 GHz. We specify the galaxy type, based on $α_{150-220}$. We combine publicly available multiwavelength data-including ultraviolet, optical, near-IR, mid-IR, far-IR, and the millimeter measurements obtained in this work-and perform spectral energy distribution (SED) fitting with CIGALE. With the radio emission decomposition advantage of CIGALE V2022.0, we identify the origins of the millimeter emissions for 69 galaxies in our sample. Our analysis also shows that millimeter data alone indicates the need for a radio synchrotron component in the SEDs that are produced by active galactic nuclei (AGNs) and/or star formation. We present SEDs and measured physical properties of these galaxies, such as the dust luminosity, AGN luminosity, the total IR luminosity, and the ratio of the IR and radio luminosity. We quantify the relationships between the total IR luminosity and the millimeter-band luminosities, which can be used in the absence of SED analysis.
△ Less
Submitted 22 March, 2023;
originally announced March 2023.
-
Pitfalls of AI classification of rare objects: Galaxy Mergers
Authors:
W. J. Pearson,
L. E. Suelves,
S. C. -C. Ho,
N. Oi,
NEP Team,
GAMA Team
Abstract:
Galaxy mergers are hugely important in our current dark matter cosmology. These powerful events cause the disruption of the merging galaxies, pushing the gas, stars and dust of the galaxies resulting in changes to morphologies. This disruption can also cause more extreme events inside the galaxies: periods of extreme star formation rates and the rapid increase in active galactic nuclei activity. H…
▽ More
Galaxy mergers are hugely important in our current dark matter cosmology. These powerful events cause the disruption of the merging galaxies, pushing the gas, stars and dust of the galaxies resulting in changes to morphologies. This disruption can also cause more extreme events inside the galaxies: periods of extreme star formation rates and the rapid increase in active galactic nuclei activity. Hence, to better understand what goes on in these rare events, we need to be able to identify statistically large samples.
In the last few years, the growth of artificial intelligence techniques has seen application to identifying galaxy mergers. These techniques have shown to be highly accurate and their application has grown beyond academic studies of ``can we?'' to deeper scientific use. However, these classifications are not without their problems.
In this proceedings, we will explore how galaxy merger classification can be improved by adding pre-extracted galaxy morphologies alongside the traditional imaging data. This demonstrates that current neural networks are not extracting all the information from the images they are given. It will also explore how the resulting samples of rare objects could be highly contaminated. This has a knock on impact on the upcoming large scale surveys like Euclid and Rubin-LSST.
△ Less
Submitted 22 February, 2023;
originally announced February 2023.
-
Understanding and predicting cadence effects in the characterization of exoplanet transits
Authors:
Julio Hernandez Camero,
Cynthia S. K. Ho,
Vincent Van Eylen
Abstract:
We investigate the effect of observing cadence on the precision of radius ratio values obtained from transit light curves by performing uniform Markov Chain Monte Carlo fits of 46 exoplanets observed by the Transiting Exoplanet Survey Satellite (TESS) in multiple cadences. We find median improvements of almost 50% when comparing fits to 20s and 120s cadence light curves to 1800s cadence light curv…
▽ More
We investigate the effect of observing cadence on the precision of radius ratio values obtained from transit light curves by performing uniform Markov Chain Monte Carlo fits of 46 exoplanets observed by the Transiting Exoplanet Survey Satellite (TESS) in multiple cadences. We find median improvements of almost 50% when comparing fits to 20s and 120s cadence light curves to 1800s cadence light curves, and of 37% when comparing 600s cadence to 1800s cadence. Such improvements in radius precision are important, for example, to precisely constrain the properties of the radius valley or to characterize exoplanet atmospheres. We also implement a numerical Information Analysis to predict the precision of parameter estimates for different observing cadences. We tested this analysis on our sample and found it reliably predicts the effect of shortening observing cadence with errors in the predicted % precision of <0,5% for most cases. We apply this method to 157 TESS object of interest that have only been observed with 1800s cadence to predict the precision improvement that could be obtained by reobservations with shorter cadences and provide the full table of expected improvements. We report the 10 planet candidates that would benefit the most from reobservations at short cadence. Our implementation of the Information Analysis for the prediction of the precision of exoplanet parameters, Prediction of Exoplanet Precisions using Information in Transit Analysis (PEPITA) is made publicly available.
△ Less
Submitted 6 February, 2023;
originally announced February 2023.
-
A deep radius valley revealed by Kepler short cadence observations
Authors:
Cynthia S. K. Ho,
Vincent Van Eylen
Abstract:
The characteristics of the radius valley, i.e., an observed lack of planets between 1.5-2 Earth radii at periods shorter than about 100 days, provide insights into the formation and evolution of close-in planets. We present a novel view of the radius valley by refitting the transits of 431 planets using Kepler 1-minute short cadence observations, the vast majority of which have not been previously…
▽ More
The characteristics of the radius valley, i.e., an observed lack of planets between 1.5-2 Earth radii at periods shorter than about 100 days, provide insights into the formation and evolution of close-in planets. We present a novel view of the radius valley by refitting the transits of 431 planets using Kepler 1-minute short cadence observations, the vast majority of which have not been previously analysed in this way. In some cases, the updated planetary parameters differ significantly from previous studies, resulting in a deeper radius valley than previously observed. This suggests that planets are likely to have a more homogeneous core composition at formation. Furthermore, using support-vector machines, we find that the radius valley location strongly depends on orbital period and stellar mass and weakly depends on stellar age, with $\partial \log {\left(R_{p, \text{valley}} \right)}/ \partial \log{P} = -0.096_{-0.027}^{+0.023}$, $\partial \log {\left(R_{p, \text{valley}} \right)}/ \partial \log{M_{\star}} = 0.231_{-0.064}^{+0.053}$, and $\partial \log {\left(R_{p, \text{valley}} \right)}/ \partial \log{\left( \text{age} \right)} = 0.033_{-0.025}^{+0.017}$. These findings favour thermally-driven mass loss models such as photoevaporation and core-powered mass loss, with a slight preference for the latter scenario. Finally, this work highlights the value of transit observations with short photometric cadence to precisely determine planet radii, and we provide an updated list of precisely and homogeneously determined parameters for the planets in our sample.
△ Less
Submitted 10 January, 2023;
originally announced January 2023.
-
The molecular gas kinematics in the host galaxy of non-repeating FRB 180924B
Authors:
Tzu-Yin Hsu,
Tetsuya Hashimoto,
Bunyo Hatsukade,
Tomotsugu Goto,
Po-Ya Wang,
Chih Teng Ling,
Simon C. -C. Ho,
Yuri Uno
Abstract:
Fast radio bursts (FRBs) are millisecond-duration transients with large dispersion measures. The origin of FRBs is still mysterious. One of the methods to comprehend FRB origin is to probe the physical environments of FRB host galaxies. Mapping molecular-gas kinematics in FRB host galaxies is critical because it results in star formation that is likely connected to the birth of FRB progenitors. Ho…
▽ More
Fast radio bursts (FRBs) are millisecond-duration transients with large dispersion measures. The origin of FRBs is still mysterious. One of the methods to comprehend FRB origin is to probe the physical environments of FRB host galaxies. Mapping molecular-gas kinematics in FRB host galaxies is critical because it results in star formation that is likely connected to the birth of FRB progenitors. However, most previous works of FRB host galaxies have focused on its stellar component. Therefore, we, for the first time, report the molecular gas kinematics in the host galaxy of the non-repeating FRB 180924B at $z= 0.3216$. Two velocity components of the CO (3-2) emission line are detected in its host galaxy with the Atacama Large Millimeter/submillimeter Array (ALMA): the peak of one component ($-155.40$ km s$^{-1}$) is near the centre of the host galaxy, and another ($-7.76$ km s$^{-1}$) is near the FRB position. The CO (3-2) spectrum shows asymmetric profiles with A$_{\rm peak}$ $=2.03\pm 0.39$, where A$_{\rm peak}$ is the peak flux density ratio between the two velocity components. The CO (3-2) velocity map also indicates an asymmetric velocity gradient from $-180$ km s$^{-1}$ to 8 km s$^{-1}$. These results indicate a disturbed kinetic structure of molecular gas in the host galaxy. Such disturbed kinetic structures are reported for repeating FRB host galaxies using HI emission lines in previous works. Our finding indicates that non-repeating and repeating FRBs could commonly appear in disturbed kinetic environments, suggesting a possible link between the gas kinematics and FRB progenitors.
△ Less
Submitted 7 December, 2022;
originally announced December 2022.
-
Supernova Axion Emissivity with $Δ(1232)$ Resonance in Heavy Baryon Chiral Perturbation Theory
Authors:
Shu-Yu Ho,
Jongkuk Kim,
Pyungwon Ko,
Jae-hyeon Park
Abstract:
In this paper, we evaluate the energy loss rate of supernovae induced by the axion emission process $π^- + p \to n + a$ with the $Δ(1232)$ resonance in the heavy baryon chiral perturbation theory for the first time. Given the axion-nucleon-$Δ$ interactions, we include the previously ignored $Δ$-mediated graphs to the $π^- + p \to n + a$ process. In particular, the $Δ^0$-mediated diagram can give a…
▽ More
In this paper, we evaluate the energy loss rate of supernovae induced by the axion emission process $π^- + p \to n + a$ with the $Δ(1232)$ resonance in the heavy baryon chiral perturbation theory for the first time. Given the axion-nucleon-$Δ$ interactions, we include the previously ignored $Δ$-mediated graphs to the $π^- + p \to n + a$ process. In particular, the $Δ^0$-mediated diagram can give a resonance contribution to the supernova axion emission rate when the center-of-mass energy of the pion and proton approaches the $Δ(1232)$ mass. With these new contributions, we find that for the typical supernova temperatures, compared with the earlier work with the axion-nucleon (and axion-pion-nucleon contact) interactions, the supernova axion emissivity can be enhanced by a factor of $\sim$4(2) in the Kim-Shifman-Vainshtein-Zakharov model and up to a factor of $\sim$5(2) in the Dine-Fischler-Srednicki-Zhitnitsky model with small $\tanβ$ values. Remarkably, we notice that the $Δ(1232)$ resonance gives a destructive contribution to the supernova axion emission rate at high supernova temperatures, which is a nontrivial result in this study.
△ Less
Submitted 10 April, 2023; v1 submitted 2 December, 2022;
originally announced December 2022.
-
Recovering Galaxy Cluster Convergence from Lensed CMB with Generative Adversarial Networks
Authors:
Liam Parker,
Dongwon Han,
Pablo Lemos Portela,
Shirley Ho
Abstract:
We present a new method which leverages conditional Generative Adversarial Networks (cGAN) to reconstruct galaxy cluster convergence from lensed CMB temperature maps. Our model is constructed to emphasize structure and high-frequency correctness relative to the Residual U-Net approach presented by Caldeira, et. al. (2019). Ultimately, we demonstrate that while both models perform similarly in the…
▽ More
We present a new method which leverages conditional Generative Adversarial Networks (cGAN) to reconstruct galaxy cluster convergence from lensed CMB temperature maps. Our model is constructed to emphasize structure and high-frequency correctness relative to the Residual U-Net approach presented by Caldeira, et. al. (2019). Ultimately, we demonstrate that while both models perform similarly in the no-noise regime (as well as after random off-centering of the cluster center), cGAN outperforms ResUNet when processing CMB maps noised with 5uK/arcmin white noise or astrophysical foregrounds (tSZ and kSZ); this out-performance is especially pronounced at high l, which is exactly the regime in which the ResUNet under-performs traditional methods.
△ Less
Submitted 16 November, 2022;
originally announced November 2022.
-
Emulating cosmological growth functions with B-Splines
Authors:
Ngai Pok Kwan,
Chirag Modi,
Yin Li,
Shirley Ho
Abstract:
In the light of GPU accelerations, sequential operations such as solving ordinary differential equations can be bottlenecks for gradient evaluations and hinder potential speed gains. In this work, we focus on growth functions and their time derivatives in cosmological particle mesh simulations and show that these are the majority time cost when using gradient based inference algorithms. We propose…
▽ More
In the light of GPU accelerations, sequential operations such as solving ordinary differential equations can be bottlenecks for gradient evaluations and hinder potential speed gains. In this work, we focus on growth functions and their time derivatives in cosmological particle mesh simulations and show that these are the majority time cost when using gradient based inference algorithms. We propose to construct novel conditional B-spline emulators which directly learn an interpolating function for the growth factor as a function of time, conditioned on the cosmology. We demonstrate that these emulators are sufficiently accurate to not bias our results for cosmological inference and can lead to over an order of magnitude gains in time, especially for small to intermediate size simulations.
△ Less
Submitted 11 November, 2022;
originally announced November 2022.
-
A Neural Network Subgrid Model of the Early Stages of Planet Formation
Authors:
Thomas Pfeil,
Miles Cranmer,
Shirley Ho,
Philip J. Armitage,
Tilman Birnstiel,
Hubert Klahr
Abstract:
Planet formation is a multi-scale process in which the coagulation of $\mathrm{μm}$-sized dust grains in protoplanetary disks is strongly influenced by the hydrodynamic processes on scales of astronomical units ($\approx 1.5\times 10^8 \,\mathrm{km}$). Studies are therefore dependent on subgrid models to emulate the micro physics of dust coagulation on top of a large scale hydrodynamic simulation.…
▽ More
Planet formation is a multi-scale process in which the coagulation of $\mathrm{μm}$-sized dust grains in protoplanetary disks is strongly influenced by the hydrodynamic processes on scales of astronomical units ($\approx 1.5\times 10^8 \,\mathrm{km}$). Studies are therefore dependent on subgrid models to emulate the micro physics of dust coagulation on top of a large scale hydrodynamic simulation. Numerical simulations which include the relevant physical effects are complex and computationally expensive. Here, we present a fast and accurate learned effective model for dust coagulation, trained on data from high resolution numerical coagulation simulations. Our model captures details of the dust coagulation process that were so far not tractable with other dust coagulation prescriptions with similar computational efficiency.
△ Less
Submitted 8 November, 2022;
originally announced November 2022.
-
${\rm S{\scriptsize IM}BIG}$: A Forward Modeling Approach To Analyzing Galaxy Clustering
Authors:
ChangHoon Hahn,
Michael Eickenberg,
Shirley Ho,
Jiamin Hou,
Pablo Lemos,
Elena Massara,
Chirag Modi,
Azadeh Moradinezhad Dizgah,
Bruno Régaldo-Saint Blancard,
Muntazir M. Abidi
Abstract:
We present the first-ever cosmological constraints from a simulation-based inference (SBI) analysis of galaxy clustering from the new ${\rm S{\scriptsize IM}BIG}$ forward modeling framework. ${\rm S{\scriptsize IM}BIG}$ leverages the predictive power of high-fidelity simulations and provides an inference framework that can extract cosmological information on small non-linear scales, inaccessible w…
▽ More
We present the first-ever cosmological constraints from a simulation-based inference (SBI) analysis of galaxy clustering from the new ${\rm S{\scriptsize IM}BIG}$ forward modeling framework. ${\rm S{\scriptsize IM}BIG}$ leverages the predictive power of high-fidelity simulations and provides an inference framework that can extract cosmological information on small non-linear scales, inaccessible with standard analyses. In this work, we apply ${\rm S{\scriptsize IM}BIG}$ to the BOSS CMASS galaxy sample and analyze the power spectrum, $P_\ell(k)$, to $k_{\rm max}=0.5\,h/{\rm Mpc}$. We construct 20,000 simulated galaxy samples using our forward model, which is based on high-resolution ${\rm Q{\scriptsize UIJOTE}}$ $N$-body simulations and includes detailed survey realism for a more complete treatment of observational systematics. We then conduct SBI by training normalizing flows using the simulated samples and infer the posterior distribution of $Λ$CDM cosmological parameters: $Ω_m, Ω_b, h, n_s, σ_8$. We derive significant constraints on $Ω_m$ and $σ_8$, which are consistent with previous works. Our constraints on $σ_8$ are $27\%$ more precise than standard analyses. This improvement is equivalent to the statistical gain expected from analyzing a galaxy sample that is $\sim60\%$ larger than CMASS with standard methods. It results from additional cosmological information on non-linear scales beyond the limit of current analytic models, $k > 0.25\,h/{\rm Mpc}$. While we focus on $P_\ell$ in this work for validation and comparison to the literature, ${\rm S{\scriptsize IM}BIG}$ provides a framework for analyzing galaxy clustering using any summary statistic. We expect further improvements on cosmological constraints from subsequent ${\rm S{\scriptsize IM}BIG}$ analyses of summary statistics beyond $P_\ell$.
△ Less
Submitted 1 November, 2022;
originally announced November 2022.
-
${\rm S{\scriptsize IM}BIG}$: Mock Challenge for a Forward Modeling Approach to Galaxy Clustering
Authors:
ChangHoon Hahn,
Michael Eickenberg,
Shirley Ho,
Jiamin Hou,
Pablo Lemos,
Elena Massara,
Chirag Modi,
Azadeh Moradinezhad Dizgah,
Bruno Régaldo-Saint Blancard,
Muntazir M. Abidi
Abstract:
Simulation-Based Inference of Galaxies (${\rm S{\scriptsize IM}BIG}$) is a forward modeling framework for analyzing galaxy clustering using simulation-based inference. In this work, we present the ${\rm S{\scriptsize IM}BIG}$ forward model, which is designed to match the observed SDSS-III BOSS CMASS galaxy sample. The forward model is based on high-resolution ${\rm Q{\scriptsize UIJOTE}}$ $N$-body…
▽ More
Simulation-Based Inference of Galaxies (${\rm S{\scriptsize IM}BIG}$) is a forward modeling framework for analyzing galaxy clustering using simulation-based inference. In this work, we present the ${\rm S{\scriptsize IM}BIG}$ forward model, which is designed to match the observed SDSS-III BOSS CMASS galaxy sample. The forward model is based on high-resolution ${\rm Q{\scriptsize UIJOTE}}$ $N$-body simulations and a flexible halo occupation model. It includes full survey realism and models observational systematics such as angular masking and fiber collisions. We present the "mock challenge" for validating the accuracy of posteriors inferred from ${\rm S{\scriptsize IM}BIG}$ using a suite of 1,500 test simulations constructed using forward models with a different $N$-body simulation, halo finder, and halo occupation prescription. As a demonstration of ${\rm S{\scriptsize IM}BIG}$, we analyze the power spectrum multipoles out to $k_{\rm max} = 0.5\,h/{\rm Mpc}$ and infer the posterior of $Λ$CDM cosmological and halo occupation parameters. Based on the mock challenge, we find that our constraints on $Ω_m$ and $σ_8$ are unbiased, but conservative. Hence, the mock challenge demonstrates that ${\rm S{\scriptsize IM}BIG}$ provides a robust framework for inferring cosmological parameters from galaxy clustering on non-linear scales and a complete framework for handling observational systematics. In subsequent work, we will use ${\rm S{\scriptsize IM}BIG}$ to analyze summary statistics beyond the power spectrum including the bispectrum, marked power spectrum, skew spectrum, wavelet statistics, and field-level statistics.
△ Less
Submitted 1 November, 2022;
originally announced November 2022.
-
$\texttt{Mangrove}$: Learning Galaxy Properties from Merger Trees
Authors:
Christian Kragh Jespersen,
Miles Cranmer,
Peter Melchior,
Shirley Ho,
Rachel S. Somerville,
Austen Gabrielpillai
Abstract:
Efficiently mapping baryonic properties onto dark matter is a major challenge in astrophysics. Although semi-analytic models (SAMs) and hydrodynamical simulations have made impressive advances in reproducing galaxy observables across cosmologically significant volumes, these methods still require significant computation times, representing a barrier to many applications. Graph Neural Networks (GNN…
▽ More
Efficiently mapping baryonic properties onto dark matter is a major challenge in astrophysics. Although semi-analytic models (SAMs) and hydrodynamical simulations have made impressive advances in reproducing galaxy observables across cosmologically significant volumes, these methods still require significant computation times, representing a barrier to many applications. Graph Neural Networks (GNNs) have recently proven to be the natural choice for learning physical relations. Among the most inherently graph-like structures found in astrophysics are the dark matter merger trees that encode the evolution of dark matter halos. In this paper we introduce a new, graph-based emulator framework, $\texttt{Mangrove}$, and show that it emulates the galactic stellar mass, cold gas mass and metallicity, instantaneous and time-averaged star formation rate, and black hole mass -- as predicted by a SAM -- with root mean squared error up to two times lower than other methods across a $(75 Mpc/h)^3$ simulation box in 40 seconds, 4 orders of magnitude faster than the SAM. We show that $\texttt{Mangrove}$ allows for quantification of the dependence of galaxy properties on merger history. We compare our results to the current state of the art in the field and show significant improvements for all target properties. $\texttt{Mangrove}$ is publicly available.
△ Less
Submitted 24 October, 2022;
originally announced October 2022.
-
Detection Rate of Fast Radio Bursts in the Milky Way with BURSTT
Authors:
Decmend Fang-Jie Ling,
Tetsuya Hashimoto,
Shotaro Yamasaki,
Tomotsugu Goto,
Seong Jin Kim,
Simon C. -C. Ho,
Tiger Y. -Y. Hsiao,
Yi Hang Valerie Wong
Abstract:
Fast radio bursts (FRBs) are intense bursts of radio emission with durations of milliseconds. Although researchers have found them happening frequently all over the sky, they are still in the dark to understand what causes the phenomena because the existing radio observatories have encountered certain challenges during the discovery of FRB progenitors. The construction of Bustling Universe Radio S…
▽ More
Fast radio bursts (FRBs) are intense bursts of radio emission with durations of milliseconds. Although researchers have found them happening frequently all over the sky, they are still in the dark to understand what causes the phenomena because the existing radio observatories have encountered certain challenges during the discovery of FRB progenitors. The construction of Bustling Universe Radio Survey Telescope in Taiwan (BURSTT) is being proposed to solve these challenges. We simulate mock Galactic FRB-like events by applying a range of spatial distributions, pulse widths and luminosity functions. The effect of turbulent Interstellar Medium (ISM) on the detectability of FRB-like events within the Milky Way plane is considered to estimate the dispersion measure and pulse scattering of mock events. We evaluate the fraction of FRB-like events in the Milky Way that are detectable by BURSTT and compare the result with those by Survey for Transient Astronomical Radio Emission 2 (STARE2) and Galactic Radio Explorer (GReX). We find that BURSTT could increase the detection rate by more than two orders of magnitude compared with STARE2 and GReX, depending on the slope of luminosity function of the events. We also investigate the influence of the specifications of BURSTT on its detection improvement. This leads to the fact that greatly higher sensitivity and improved coverage of the Milky Way plane have significant effects on the detection improvement of BURSTT. We find that the upgrade version of BURSTT, BURSTT-2048 could increase the detection rate of faint Galactic FRB-like events by a factor of 3.
△ Less
Submitted 20 October, 2022;
originally announced October 2022.
-
Particle clustering in turbulence: Prediction of spatial and statistical properties with deep learning
Authors:
Yan-Mong Chan,
Natascha Manger,
Yin Li,
Chao-Chin Yang,
Zhaohuan Zhu,
Philip J. Armitage,
Shirley Ho
Abstract:
We investigate the utility of deep learning for modeling the clustering of particles that are aerodynamically coupled to turbulent fluids. Using a Lagrangian particle module within the Athena++ hydrodynamics code, we simulate the dynamics of particles in the Epstein drag regime within a periodic domain of isotropic forced hydrodynamic turbulence. This setup is an idealized model relevant to the co…
▽ More
We investigate the utility of deep learning for modeling the clustering of particles that are aerodynamically coupled to turbulent fluids. Using a Lagrangian particle module within the Athena++ hydrodynamics code, we simulate the dynamics of particles in the Epstein drag regime within a periodic domain of isotropic forced hydrodynamic turbulence. This setup is an idealized model relevant to the collisional growth of micron to mm-sized dust particles in early stage planet formation. The simulation data are used to train a U-Net deep learning model to predict gridded three-dimensional representations of the particle density and velocity fields, given as input the corresponding fluid fields. The trained model qualitatively captures the filamentary structure of clustered particles in a highly non-linear regime. We assess model fidelity by calculating metrics of the density field (the radial distribution function) and of the velocity field (the relative velocity and the relative radial velocity between particles). Although trained only on the spatial fields, the model predicts these statistical quantities with errors that are typically <10%. Our results suggest that, given appropriately expanded training data, deep learning could complement direct numerical simulations in predicting particle clustering within turbulent flows.
△ Less
Submitted 6 January, 2024; v1 submitted 5 October, 2022;
originally announced October 2022.
-
Constraining the Hubble constant and its lower limit from the proper motion of extragalactic radio jets
Authors:
Tiger Yu-Yang Hsiao,
Tomotsugu Goto,
Tetsuya Hashimoto,
Daryl Joe D. Santos,
Yi Hang Valerie Wong,
Seong Jin Kim,
Bjorn Jasper R. Raquel,
Simon C. -C. Ho,
Bo-Han Chen,
Ece Kilerci,
Ting-Yi Lu,
Alvina Y. L. On,
Yu-Wei Lin,
Cossas K. -W. Wu
Abstract:
The Hubble constant ($H_{0}$) is a measurement to describe the expansion rate of the Universe in the current era. However, there is a $4.4σ$ discrepancy between the measurements from the early Universe and the late Universe. In this research, we propose a model-free and distance-free method to constrain $H_{0}$. Combining Friedman-Lemaître-Robertson-Walker cosmology with geometrical relation of th…
▽ More
The Hubble constant ($H_{0}$) is a measurement to describe the expansion rate of the Universe in the current era. However, there is a $4.4σ$ discrepancy between the measurements from the early Universe and the late Universe. In this research, we propose a model-free and distance-free method to constrain $H_{0}$. Combining Friedman-Lemaître-Robertson-Walker cosmology with geometrical relation of the proper motion of extragalactic jets, the lower limit ($H_{\rm 0,min}$) of $H_{0}$ can be determined using only three cosmology-free observables: the redshifts of the host galaxies, as well as the approaching and receding angular velocities of radio jets. Using these, we propose to use the Kolmogorov-Smirnov test (K-S test) between cumulative distribution functions of $H_{\rm 0,min}$ to differentiate cosmology. We simulate 100, 200, and 500 extragalactic jets with 3 levels of accuracy of the proper motion ($μ_{a}$ and $μ_{r}$), at $10\%$, $5\%$, and $1\%$, corresponding to the accuracies of the current and future radio interferometers. We perform K-S tests between the simulated samples as theoretical distributions with different $H_{0}$ and power-law index of velocity distribution of jets and mock observational data. Our result suggests increasing sample sizes leads to tighter constraints on both power-law index and the Hubble constant at moderate accuracy (i.e., $10\%$ and $5\%$) while at $1\%$ accuracy, increasing sample sizes leads to tighter constraints on power-law index more. Improving accuracy results in better constraints in the Hubble constant compared with the power-law index in all cases but it alleviates the degeneracy.
△ Less
Submitted 11 September, 2022;
originally announced September 2022.
-
Can luminous Lyman alpha emitters at $z$ $\simeq$ 5.7 and $z$ $\simeq$ 6.6 suppress star formation?
Authors:
Daryl Joe D. Santos,
Tomotsugu Goto,
Tetsuya Hashimoto,
Seong Jin Kim,
Ting-Yi Lu,
Yi-Hang Valerie Wong,
Simon C. -C. Ho,
Tiger Y. -Y. Hsiao
Abstract:
Addressing how strong UV radiation affects galaxy formation is central to understanding their evolution. The quenching of star formation via strong UV radiation (from starbursts or AGN) has been proposed in various scenes to solve certain astrophysical problems. Around luminous sources, some evidence of decreased star formation has been found but is limited to a handful of individual cases. No dir…
▽ More
Addressing how strong UV radiation affects galaxy formation is central to understanding their evolution. The quenching of star formation via strong UV radiation (from starbursts or AGN) has been proposed in various scenes to solve certain astrophysical problems. Around luminous sources, some evidence of decreased star formation has been found but is limited to a handful of individual cases. No direct, conclusive evidence on the actual role of strong UV radiation in quenching star formation has been found. Here we present statistical evidence of decreased number density of faint (AB magnitude $\geq$ 24.75 mag) Lyαemitters (LAEs) around bright (AB magnitude < 24.75 mag) LAEs even when the radius goes up to 10 pMpc for $z$ $\simeq$ 5.7 LAEs. A similar trend is found for z $\simeq$ 6.6 LAEs but only within 1 pMpc radius from the bright LAEs. We use a large sample of 1077 (962) LAEs at $z$ $\simeq$ 5.7 ($z$ $\simeq$ 6.6) selected in total areas of 14 (21) deg$^2$ with Subaru/Hyper Suprime-Cam narrow-band data, and thus, the result is of statistical significance for the first time at these high redshift ranges. A simple analytical calculation indicates that the radiation from the central LAE is not enough to suppress LAEs with AB mag $\geq$ 24.75 mag around them, suggesting additional physical mechanisms we are unaware of are at work. Our results clearly show that the environment is at work for the galaxy formation at $z$ $\sim$ 6 in the Universe.
△ Less
Submitted 6 September, 2022;
originally announced September 2022.
-
The SZ flux-mass ($Y$-$M$) relation at low halo masses: improvements with symbolic regression and strong constraints on baryonic feedback
Authors:
Digvijay Wadekar,
Leander Thiele,
J. Colin Hill,
Shivam Pandey,
Francisco Villaescusa-Navarro,
David N. Spergel,
Miles Cranmer,
Daisuke Nagai,
Daniel Angl�s-Alc�zar,
Shirley Ho,
Lars Hernquist
Abstract:
Feedback from active galactic nuclei (AGN) and supernovae can affect measurements of integrated SZ flux of halos ($Y_\mathrm{SZ}$) from CMB surveys, and cause its relation with the halo mass ($Y_\mathrm{SZ}-M$) to deviate from the self-similar power-law prediction of the virial theorem. We perform a comprehensive study of such deviations using CAMELS, a suite of hydrodynamic simulations with exten…
▽ More
Feedback from active galactic nuclei (AGN) and supernovae can affect measurements of integrated SZ flux of halos ($Y_\mathrm{SZ}$) from CMB surveys, and cause its relation with the halo mass ($Y_\mathrm{SZ}-M$) to deviate from the self-similar power-law prediction of the virial theorem. We perform a comprehensive study of such deviations using CAMELS, a suite of hydrodynamic simulations with extensive variations in feedback prescriptions. We use a combination of two machine learning tools (random forest and symbolic regression) to search for analogues of the $Y-M$ relation which are more robust to feedback processes for low masses ($M\lesssim 10^{14}\, h^{-1} \, M_\odot$); we find that simply replacing $Y\rightarrow Y(1+M_*/M_\mathrm{gas})$ in the relation makes it remarkably self-similar. This could serve as a robust multiwavelength mass proxy for low-mass clusters and galaxy groups. Our methodology can also be generally useful to improve the domain of validity of other astrophysical scaling relations.
We also forecast that measurements of the $Y-M$ relation could provide percent-level constraints on certain combinations of feedback parameters and/or rule out a major part of the parameter space of supernova and AGN feedback models used in current state-of-the-art hydrodynamic simulations. Our results can be useful for using upcoming SZ surveys (e.g., SO, CMB-S4) and galaxy surveys (e.g., DESI and Rubin) to constrain the nature of baryonic feedback. Finally, we find that the an alternative relation, $Y-M_*$, provides complementary information on feedback than $Y-M$
△ Less
Submitted 28 April, 2023; v1 submitted 5 September, 2022;
originally announced September 2022.