Researcher profile

Hy Trac

Hy Trac contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

10 published item(s)

preprint2026arXiv

Cross-correlating the patchy screening and kinetic Sunyaev-Zel'dovich effects as a new probe of reionization

The kinetic Sunyaev-Zel'dovich effect (kSZ) and patchy screening effect are two complementary cosmic microwave background (CMB) probes of the reionization era. The kSZ effect is a relatively strong signal, but is difficult to disentangle from other sources of temperature anisotropy, whereas patchy screening is weaker but can be reconstructed using the cleaner polarization channel. Here, we explore the potential of using upcoming CMB surveys to correlate a reconstructed map of patchy screening with (the square of) the kSZ map, and what a detection of this cross-correlation would mean for reionization science. To do this, we use simulations and theory to quantify the contributions to this signal from different redshifts. We then use the expected survey properties for CMB-S4 and CMB-HD to make detection forecasts. We find that, for or our fiducial reionization scenario, CMB-S4 will obtain a hint of this signal at up to 1.8$σ$, and CMB-HD will detect it at up to 14$σ$. We explore the physical interpretation of the signal and find that it is uniquely sensitive to the first half of reionization and to the bispectrum of the ionized gas distribution.

preprint2022arXiv

AMBER: A Semi-numerical Abundance Matching Box for the Epoch of Reionization

The Abundance Matching Box for the Epoch of Reionization (AMBER) is a semi-numerical code for modeling the cosmic dawn. The new algorithm is not based on the excursion set formalism for reionization, but takes the novel approach of calculating the reionization-redshift field $z_\mathrm{re}(\boldsymbol{x})$ assuming that hydrogen gas encountering higher radiation intensity are photoionized earlier. Redshift values are assigned while matching the abundance of ionized mass according to a given mass-weighted ionization fraction $\bar{x}_\mathrm{i}(z)$. The code has the unique advantage of allowing users to directly specify the reionization history through the redshift midpoint $z_\mathrm{mid}$, duration $Δ_\mathrm{z}$, and asymmetry $A_\mathrm{z}$ input parameters. The reionization process is further controlled through the minimum halo mass $M_\mathrm{min}$ for galaxy formation and the radiation mean free path $l_\mathrm{mfp}$ for radiative transfer. We implement improved methods for constructing density, velocity, halo, and radiation fields, which are essential components for modeling reionization observables. We compare AMBER with two other semi-numerical methods and find that our code more accurately reproduces the results from radiation-hydrodynamic simulations. The parallelized code is over four orders of magnitude faster than radiative transfer simulations and will efficiently enable large-volume models, full-sky mock observations, and parameter-space studies. AMBER will be made publicly available to facilitate and transform studies of the EoR.

preprint2022arXiv

The Dynamical Mass of the Coma Cluster from Deep Learning

In 1933, Fritz Zwicky's famous investigations of the mass of the Coma cluster led him to infer the existence of dark matter \cite{1933AcHPh...6..110Z}. His fundamental discoveries have proven to be foundational to modern cosmology; as we now know such dark matter makes up 85\% of the matter and 25\% of the mass-energy content in the universe. Galaxy clusters like Coma are massive, complex systems of dark matter in addition to hot ionized gas and thousands of galaxies, and serve as excellent probes of the dark matter distribution. However, empirical studies show that the total mass of such systems remains elusive and difficult to precisely constrain. Here, we present new estimates for the dynamical mass of the Coma cluster based on Bayesian deep learning methodologies developed in recent years. Using our novel data-driven approach, we predict Coma's $\mthc$ mass to be $10^{15.10 \pm 0.15}\ \hmsun$ within a radius of $1.78 \pm 0.03\ h^{-1}\mathrm{Mpc}$ of its center. We show that our predictions are rigorous across multiple training datasets and statistically consistent with historical estimates of Coma's mass. This measurement reinforces our understanding of the dynamical state of the Coma cluster and advances rigorous analyses and verification methods for empirical applications of machine learning in astronomy.

preprint2021arXiv

A Hydro-Particle-Mesh Code for Efficient and Rapid Simulations of the Intracluster Medium

We introduce the cosmological HYPER code based on an innovative hydro-particle-mesh (HPM) algorithm for efficient and rapid simulations of gas and dark matter. For the HPM algorithm, we update the approach of Gnedin & Hui (1998) to expand the scope of its application from the lower-density intergalactic medium (IGM) to the higher-density intracluster medium (ICM). While the original algorithm tracks only one effective particle species, the updated version separately tracks the gas and dark matter particles as they do not exactly trace each other on small scales. For the approximate hydrodynamics solver, the pressure term in the gas equations of motion is calculated using robust physical models. In particular, we use a dark matter halo model, ICM pressure profile, and IGM temperature-density relation, all of which can be systematically varied for parameter-space studies. We show that the HYPER simulation results are in good agreement with the halo model expectations for the density, temperature, and pressure radial profiles. Simulated galaxy cluster scaling relations for Sunyaev-Zel'dovich (SZ) and X-ray observables are also in good agreement with mean predictions, with scatter comparable to that found in hydrodynamic simulations. HYPER also produces lightcone catalogs of dark matter halos and full-sky tomographic maps of the lensing convergence, SZ effect, and X-ray emission. These simulation products are useful for testing data analysis pipelines, generating training data for machine learning, understanding selection and systematic effects, and for interpreting astrophysical and cosmological constraints.

preprint2021arXiv

A Machine Learning Approach for Dynamical Mass Measurements of Galaxy Clusters

We present a modern machine learning approach for cluster dynamical mass measurements that is a factor of two improvement over using a conventional scaling relation. Different methods are tested against a mock cluster catalog constructed using halos with mass >= 10^14 Msolar/h from Multidark's publicly-available N-body MDPL halo catalog. In the conventional method, we use a standard M(sigma_v) power law scaling relation to infer cluster mass, M, from line-of-sight (LOS) galaxy velocity dispersion, sigma_v. The resulting fractional mass error distribution is broad, with width=0.87 (68% scatter), and has extended high-error tails. The standard scaling relation can be simply enhanced by including higher-order moments of the LOS velocity distribution. Applying the kurtosis as a correction term to log(sigma_v) reduces the width of the error distribution to 0.74 (16% improvement). Machine learning can be used to take full advantage of all the information in the velocity distribution. We employ the Support Distribution Machines (SDMs) algorithm that learns from distributions of data to predict single values. SDMs trained and tested on the distribution of LOS velocities yield width=0.46 (47% improvement). Furthermore, the problematic tails of the mass error distribution are effectively eliminated. Decreasing cluster mass errors will improve measurements of the growth of structure and lead to tighter constraints on cosmological parameters.

preprint2021arXiv

Hydrodynamic Response of the Intergalactic Medium to Reionization

The intergalactic medium is expected to clump on scales down to $10^4-10^8$ M$_{\odot}$ before the onset of reionization. The impact of these small-scale structures on reionization is poorly understood despite the modern understanding that gas clumpiness limits the growth of H II regions. We use a suite of radiation-hydrodynamics simulations that capture the $\sim 10^4$ $M_\odot$ Jeans mass of unheated gas to study density fluctuations during reionization. Our simulations track the complex ionization and hydrodynamical response of gas in the wake of ionization fronts. The clumping factor of ionized gas (proportional to the recombination rate) rises to a peak value of $5-20$ approximately $Δt = 10$ Myr after ionization front passage, depending on the incident intensity, redshift, and degree to which the gas had been pre-heated by the first X-ray sources. The clumping factor reaches its relaxed value of $\approx 3$ by $Δt = 300$ Myr. The mean free path of Lyman-limit photons evolves in unison, being up to several times shorter in un-relaxed, recently reionized regions compared to those that were reionized much earlier. Assessing the impact of this response on the global reionizaton process, we find that un-relaxed gaseous structures boost the total number of recombinations by $\approx 50$ % and lead to spatial fluctuations in the mean free path that persist appreciably for several hundred million years after the completion of reionization.

preprint2021arXiv

SCORCH. III. Analytical Models of Reionization with Varying Clumping Factors

In the Simulations and Constructions of the Reionization of Cosmic Hydrogen (SCORCH) project, we compare analytical models of the hydrogen ionization fraction with radiation-hydrodynamic simulations. We derive analytical models of the mass-weighted hydrogen ionization fraction from the local ionization balance equations as a more accurate alternative to the widely adopted model based on the volume filling factor. In particular, our model has a recombination term quadratic in the ionization fraction, which is consistent with the two-body interaction nature of recombination. Then, we use the radiation-hydrodynamic simulations to study the clumping factors needed to solve the analytical equations, and provide accurate fitting functions. We find that the ionized hydrogen clumping factors from our radiative transfer simulations are significantly different than those from other simulations that use a uniform photoionization background. In addition to redshift dependence, we also see the dependence of ionized hydrogen clumping factor on ionization fraction, and we incorporate this into our fits. We calculate the reionization histories using our analytical models and clumping factors and compare with widely adopted models, and all of our models achieve $<7\%$ difference from simulation results while the other models have $>20\%$ deviations. The Thomson optical depths from reionization calculated from our analytical models result in $<5\%$ deviation from simulations, while the previous analytical models have $>20\%$ difference in and could result in biased conclusions of the IGM reionization.

preprint2021arXiv

The Role of Machine Learning in the Next Decade of Cosmology

In recent years, machine learning (ML) methods have remarkably improved how cosmologists can interpret data. The next decade will bring new opportunities for data-driven cosmological discovery, but will also present new challenges for adopting ML methodologies and understanding the results. ML could transform our field, but this transformation will require the astronomy community to both foster and promote interdisciplinary research endeavors.

preprint2020arXiv

A Model-Insensitive Baryon Acoustic Oscillation Feature in the 21 cm Signal from Reionization

We examine the impact of baryon-dark matter relative velocities on intergalactic small-scale structure and the 21 cm signal during reionization. Streaming velocities reduced clumping in the intergalactic medium (IGM) on mass scales of $\sim 10^4 - 10^8$ M$_{\odot}$. This effect produced a distinct baryon acoustic oscillation (BAO) feature in the 21 cm power spectrum at wave numbers $k\sim 0.1$ h/Mpc, near which forthcoming surveys will be most sensitive. In contrast to the highly uncertain impact of streaming velocities on star formation, the effect on clumping is better constrained because it is set mainly by cosmology and straightforward gas dynamics. We quantify the latter using coupled radiation-hydrodynamic simulations that capture the Jeans scale of pre-reionization gas. The clumping factor of ionized gas is reduced by 5-10\% in regions with RMS streaming velocities. The suppression peaks $\approx 5$ Myr after a region is reionized, but disappears within 200 Myr due to pressure smoothing. We model the corresponding impact on the 21 cm signal and find that the BAO feature is most likely to appear at $\approx$ 10 \% ionization. During this phase, the feature may appear at the 1 \% (5 \%) level at $k \sim 0.1 (0.06)$ h/Mpc with an amplitude that varies by a factor of $< 10$ across a range of reionization histories. We also provide a model for the signal originating from streaming velocity&#39;s impact on ionizing sources, which can vary by 4 orders of magnitude depending on highly uncertain source properties. We find that the clumping signal probably dominates the source one unless Population III star formation in $10^6 - 10^8$ M$_{\odot}$ halos contributed significantly to the first 10\% of reionization.

preprint2020arXiv

Aging Halos: Implications of the Magnitude Gap on Conditional Statistics of Stellar and Gas Properties of Massive Halos

Cold dark matter model predicts that the large-scale structure grows hierarchically. Small dark matter halos form first. Then, they grow gradually via continuous merger and accretion. These halos host the majority of baryonic matter in the Universe in the form of hot gas and cold stellar phase. Determining how baryons are partitioned into these phases requires detailed modeling of galaxy formation and their assembly history. It is speculated that formation time of the same mass halos might be correlated with their baryonic content. To evaluate this hypothesis, we employ halos of mass above $10^{14}\,M_{\odot}$ realized by TNG300 solution of the IllustrisTNG project. Formation time is not directly observable. Hence, we rely on the magnitude gap between the brightest and the fourth brightest halo galaxy member, which is shown that traces formation time of the host halo. We compute the conditional statistics of the stellar and gas content of halos conditioned on their total mass and magnitude gap. We find a strong correlation between magnitude gap and gas mass, BCG stellar mass, and satellite galaxies stellar mass, but not the total stellar mass of halo. Conditioning on the magnitude gap can reduce the scatter about halo property--halo mass relation and has a significant impact on the conditional covariance. Reduction in the scatter can be as significant as 30%, which implies more accurate halo mass prediction. Incorporating the magnitude gap has a potential to improve cosmological constraints using halo abundance and allows us to gain insight into the baryon evolution within these systems.