Showing 1–50 of 334 results for author: Ho, S

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.

Comments: 15 pages, 20 figures, 4 tables. Accepted for publication in MNRAS. The 3 min summary: https://www.youtube.com/watch?v=mWUebbgUOh8

{\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.

Comments: 15 pages, 6 figures

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.

Comments: 15 pages, 10 figures, 8 tables, 2 appendices, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 686, A94 (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.

Comments: Accepted for publication in MNRAS

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.

Comments: 22 pages, 22 figures, 7 tables. Accepted for publication in MNRAS. A summary video can be found at https://youtu.be/TRb6bjmGfOU

${\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.

Comments: 23 pages, 12 figures, 2 tables

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.

Comments: 19 pages, accepted for publication in MNRAS

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.

Comments: 12 pages, 20 figures, 4 tables. Accepted by MNRAS. A summary video is at https://www.youtube.com/watch?v=UtPaVTFM4f8&ab_channel=NTHUCosmology

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.

Comments: 15 pages, 12 figures, published in MNRAS, https://doi.org/10.1093/mnras/stad3499. A summary video is https://youtu.be/Md6wragrYyM

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.

Comments: 18 pages, 7 figures, 3 tables, 4 pages of appendices. Submitted to ApJ

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.

Comments: 11 pages, 9 figures, Accepted for the NeurIPS 2023 AI4Science Workshop

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.

Comments: 14 pages, 4 figures. A previous version of the paper was published in the ICML 2023 Workshop on Machine Learning for Astrophysics

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.

Comments: 11+5 pages, 8+2 figures, published in Physical Review D

${\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.

Comments: 13 pages, 5 figures, submitted to Nature Astronomy, comments welcome

${\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.

Comments: 13 pages, 7 figures, submitted to PRD, comments welcome

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.

Comments: 18 pages, 5 figures, 1 table, matches the version accepted for publication in JCAP

Report number: KIAS-P23045

Journal ref: JCAP 04(2024)055

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.

Comments: 18 pages, accepted in Monthly Notices of the Royal Astronomical Society, Presented at the NeurIPS 2023 AI4Science Workshop

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.

Comments: Accepted at ICML 2022 Workshop on Machine Learning for Astrophysics, Baltimore, Maryland, USA, 2022

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.

Comments: 10+8 pages, 11+7 figures

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.

Comments: 24 pages, 14 figures, accepted for publication in MNRAS. For summary video, please see https://www.youtube.com/watch?v=wYx6t2G__84&list=PLOpYDs2PkYlYIiKDjDz6r6aKXcXdJZXYb&index=13&ab_channel=NCHUAstronomy

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

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

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

Submitted 3 July, 2023; originally announced July 2023.

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

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.

Comments: 5 pages, 6 figures and 1 table; accepted for publication in MNRAS; A summary video is available at https://youtu.be/PQW79tuS0mI

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

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

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

Submitted 11 April, 2023; originally announced April 2023.

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

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

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

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

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

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

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

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.

Comments: Accepted for publication in ApJ. A summary video is available at this https://youtu.be/yivrtvuMDHE

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.

Comments: 22 pages, 12 figures

Journal ref: The Astrophysical Journal, Volume 965, Number 2, April 2024

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.

Comments: 8 pages, 7 figures

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.

Comments: 5 pages, 1 figure, accepted for publication in MNRAS. A summary video is available: https://www.youtube.com/watch?v=_xaRCnLMjsY&t=21s&ab_channel=NCHUAstronomy

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.

Comments: 8 pages, 9 figures, accepted for publication in MNRAS. For summary video, please see https://www.youtube.com/watch?v=1PPCMKcrIf4&list=PLOpYDs2PkYlYIiKDjDz6r6aKXcXdJZXYb&index=14&ab_channel=NCHUAstronomy

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.

Comments: Camera ready version for NeurIPS 2023 AI for Science workshop https://ai4sciencecommunity.github.io/neurips23.html

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.

Comments: 29 pages, 6 tables, 11 figures, 1 appendix with 1 further figure. Accepted for publication in ApJS

Journal ref: The Astrophysical Journal Supplement Series, Volume 265, 2023, Number 2

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.

Comments: 4 pages, 3 figures, proceedings for EAS 2022 S11, to be published in Memorie della SAIt

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.

Comments: 15 pages, 9 figures. Accepted for publication in MNRAS

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.

Comments: 19 pages, 17 figures, accepted in MNRAS

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.

Comments: 5 pages, 4 figures, Accepted for publication in MNRAS, https://www.youtube.com/watch?v=CldxLE7Pdwk&t=1s

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.

Comments: 17 pages, 4 figures, typo is fixed, references updated. Matches version published in Physical Review D

Report number: KIAS-P22082

Journal ref: Phys.Rev.D 107 (2023) 7, 075002

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.

Comments: 6 pages, 4 figures, accepted at the Machine Learning and the Physical Sciences workshop, NeurIPS 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.

Comments: 9 pages, 5 figures

${\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.

Comments: 28 pages, 6 figures

$\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.

Comments: 15 pages, 9 figures, 3 tables, 10 pages of Appendices. Accepted for publication in ApJ

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.

Comments: 9 pages, 9 figures, 2 tables, accepted for publication in the Monthly Notices of the Royal Astronomical Society

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.

Comments: 20 pages, 13 figures; Accepted to ApJ on 2023 October 30

Journal ref: ApJ 960 19 (2024)

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.

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

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.

Comments: Accepted for publication at MNRAS

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.

Comments: Version appearing in MNRAS. Minor change to Fig.6 and added Fig. A5 compared to the previous version. 7+5 pages. The code and data associated with this paper are available at https://github.com/JayWadekar/ScalingRelations_ML

Journal ref: MNRAS, 522, 2628 (2023)

Source counts at 7.7 to 21 $μ$m in CEERS field with James Webb Space Telescope

Authors: Cossas K. -W. Wu, Chih-Teng Ling, Tomotsugu Goto, Ece Kilerci, Seong Jin Kim, Tetsuya Hashimoto, Yu-Wei Lin, Po-Ya Wang, Yuri Uno, Simon C. -C. Ho, Tiger Yu-Yang Hsiao

Abstract: Source counts -- the number density of sources as a function of flux density -- represent one of the fundamental metrics in observational cosmology due to their straightforward and simple nature. It is an important tool that provides information on galaxy formation and evolution. Source counting is a direct measurement. Compared to advanced analyzes that require more observational input such as lu… ▽ More Source counts -- the number density of sources as a function of flux density -- represent one of the fundamental metrics in observational cosmology due to their straightforward and simple nature. It is an important tool that provides information on galaxy formation and evolution. Source counting is a direct measurement. Compared to advanced analyzes that require more observational input such as luminosity/mass functions, it is less affected by any cosmological parameter assumptions or any errors propagated from luminosities. In this study, we present source counts at the six mid-infrared bands, i.e., 7.7, 10, 12.8, 15, 18, and 21 $μ$m from the MIR instrument of the James Webb Space Telescope (JWST). Contrasted with the infrared source counts achieved by prior generations of infrared space telescopes, our source counts delve up to $\sim$100 times deeper, showcasing the exceptional sensitivity of the JWST, and aligning with the model predictions based on preceding observations. In a follow-up study, we utilize our source counts to establish a new IR galaxy population evolutionary model that provides a physical interpretation. △ Less

Submitted 13 June, 2023; v1 submitted 5 September, 2022; originally announced September 2022.

Comments: 10 pages, 24 figures, accepted for publication in MNRAS, https://youtu.be/xlX25O8T89g

Galaxy source counts at 7.7 $μ$m, 10 $μ$m and 15 $μ$m with the James Webb Space Telescope

Authors: Chih-Teng Ling, Seong Jin Kim, Cossas K. -W. Wu, Tomotsugu Goto, Ece Kilerci, Tetsuya Hashimoto, Yu-Wei Lin, Po-Ya Wang, Simon C. -C. Ho, Tiger Yu-Yang Hsiao

Abstract: We present mid-infrared galaxy number counts based on the Early Release Observations obtained by the James Webb Space Telescope (JWST) at 7.7-, 10- and 15-$μ$m (F770W, F1000W and F1500W, respectively) bands of the Mid-Infrared Instrument (MIRI). Due to the superior sensitivity of JWST, the 80 percent completeness limits reach 0.32, 0.79 and 2.0 $μ$Jy in F770W, F1000W and F1500W filters, respective… ▽ More We present mid-infrared galaxy number counts based on the Early Release Observations obtained by the James Webb Space Telescope (JWST) at 7.7-, 10- and 15-$μ$m (F770W, F1000W and F1500W, respectively) bands of the Mid-Infrared Instrument (MIRI). Due to the superior sensitivity of JWST, the 80 percent completeness limits reach 0.32, 0.79 and 2.0 $μ$Jy in F770W, F1000W and F1500W filters, respectively, i.e., $\sim$100 times deeper than previous space infrared telescopes such as Spitzer or AKARI. The number counts reach much deeper than the broad bump around $0.05\sim0.5$ mJy due to polycyclic aromatic hydrocarbon (PAH) emissions. An extrapolation towards fainter flux from the evolutionary models in the literature agrees amazingly well with the new data, where the extrapolated faint-end of infrared luminosity functions combined with the cosmic star-formation history to higher redshifts can reproduce the deeper number counts by JWST. Our understanding of the faint infrared sources has been confirmed by the observed data due to the superb sensitivity of JWST. △ Less

Submitted 21 September, 2022; v1 submitted 8 August, 2022; originally announced August 2022.

Comments: 6 pages, 8 figures. Accepted for publication in MNRAS

Development of digital sideband separating down-conversion for Yuan-Tseh Lee Array

Authors: Chao-Te Li, Derek Kubo, Jen-Chieh Cheng, John Kuroda, Ranjani Srinivasan, Solomon Ho, Kim Guzzino, Ming-Tang Chen

Abstract: This report presents a down-conversion method involving digital sideband separation for the Yuan Tseh Lee Array to double the processing bandwidth. The receiver consists of a MMIC HEMT LNA frontend operating at a wavelength of 3 mm, and sub-harmonic mixers that output signals at intermediate frequencies of 2 - 18 GHz. The sideband separation scheme involves an analog 90 degree hybrid followed by t… ▽ More This report presents a down-conversion method involving digital sideband separation for the Yuan Tseh Lee Array to double the processing bandwidth. The receiver consists of a MMIC HEMT LNA frontend operating at a wavelength of 3 mm, and sub-harmonic mixers that output signals at intermediate frequencies of 2 - 18 GHz. The sideband separation scheme involves an analog 90 degree hybrid followed by two mixers that provide down conversion of the IF signal to a pair of in phase (I) and quadrature (Q) signals in baseband. The I and Q baseband signals are digitized using 5 Giga sample per second analog to digital converters. A second hybrid is digitally implemented using field programmable gate arrays to produce two sidebands, each with a bandwidth of 1.6 GHz. The 2 x 1.6 GHz band can be tuned to cover any 3.2 GHz window within the aforementioned IF range of the array. Sideband rejection ratios (SRRs) above 20 dB can be obtained across the 3.2 GHz bandwidth by equalizing the power and delay between the I and Q baseband signals. Furthermore, SRRs above 30 dB can be achieved when calibration is applied. △ Less

Submitted 4 August, 2022; originally announced August 2022.

Comments: 8 pages, 12 figures

Journal ref: Proc. SPIE 9914, Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy VIII, 99140E (19 July 2016)