Researcher profile

Aaron Smith

Aaron Smith contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
7topics
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

11 published item(s)

preprint2025arXiv

RIGEL: Feedback-regulated cloud-scale star formation efficiency in a simulated dwarf galaxy merger

Major mergers of galaxies are likely to trigger bursty star formation activities. The accumulation of dense gas and the boost of star formation efficiency (SFE) are considered to be the two main drivers of starbursts. However, it remains unclear how each process operates on the scale of individual star-forming clouds. Here, we present a high-resolution (2 Msun) RHD simulation of a gas-rich dwarf galaxy merger using the RIGEL model to investigate how mergers affect the properties of the structure of dense star-forming gas and the cloud-scale SFE. We tracked the evolution of sub-virial dense clouds in the simulation by mapping them across successive snapshots taken at intervals of 0.2 Myr. We find that the merger triggers a 130 fold increase in the SFR and shortens the galaxy-wide gas-depletion time by two orders of magnitude compared to those in two isolated galaxies. However, the depletion time of individual clouds and their lifetime distribution remained unchanged over the simulation period. The cloud life cycles and cloud-scale SFE are determined by the local stellar feedback rather than such environmental factors as tidal fields regardless of the merger process, and the integrated SFE ($ε_{\rm int}$) of clouds in complex environments remains well-described by an $ε_{\rm int}-Σ_{\rm tot}$ relation found in idealized isolated-cloud experiments. During the peak of the starburst, the media SFE was lower by only 0.17-0.33 dex compared to the value when the galaxies were not interacting. The merger boosts the SFR through the accumulation and compression of dense gas fueling star formation. Strong tidal torques assemble $>10^5$ Msun clouds, which seed massive star clusters. The average separation between star-forming clouds decreases during the merger, which in turn decreases the cloud--cluster spatial de-correlation from >1 kpc to 0.1 kpc depicted in tuning fork diagrams.

preprint2025arXiv

The THESAN-ZOOM project: Star-formation efficiencies in high-redshift galaxies

Recent JWST observations hint at unexpectedly intense cosmic star-formation in the early Universe, often attributed to enhanced star-formation efficiencies (SFEs). Here, we analyze the SFE in THESAN-ZOOM, a novel zoom-in radiation-hydrodynamic simulation campaign of high-redshift ($z \gtrsim 3$) galaxies employing a state-of-the-art galaxy formation model resolving the multiphase interstellar medium (ISM). The halo-scale SFE ($ε^{\ast}_{\rm halo}$) - the fraction of baryons accreted by a halo that are converted to stars - follows a double power-law dependence on halo mass, with a mild redshift evolution above $M_{\rm halo} \gtrsim 10^{9.5}\,{\rm M}_{\odot}$. The power-law slope is roughly $1/3$ at large halo masses, consistent with expectations when gas outflows are momentum-driven. At lower masses, the slope is roughly $2/3$ and is more aligned with the energy-driven outflow scenario. $ε^{\ast}_{\rm halo}$ is a factor of $2-3$ larger than commonly assumed in empirical galaxy-formation models at $M_{\rm halo} \lesssim 10^{11}\,{\rm M}_{\odot}$. On galactic (kpc) scales, the Kennicutt-Schmidt (KS) relation of neutral gas is universal in THESAN-ZOOM, following $Σ_{\rm SFR} \propto Σ_{\rm gas}^2$, indicative of a turbulent energy balance in the ISM maintained by stellar feedback. The rise of $ε^{\ast}_{\rm halo}$ with halo mass can be traced primarily to increasing gas surface densities in massive galaxies, while the underlying KS relation and neutral, star-forming gas fraction remain unchanged. Although the increase in $ε^{\ast}_{\rm halo}$ with redshift is relatively modest, it is sufficient to explain the large observed number density of UV-bright galaxies at $z \gtrsim 12$. However, reproducing the brightest sources at $M_{\rm UV} \lesssim -21$ may require extrapolating the SFE beyond the halo mass range directly covered by THESAN-ZOOM.

preprint2023arXiv

There and back again: understanding the critical properties of backsplash galaxies

Backsplash galaxies are galaxies that once resided inside a cluster, and have migrated back oustide as they move towards the apocentre of their orbit. The kinematic properties of these galaxies are well understood, thanks to the significant study of backsplashers in dark matter-only simulations, but their intrinsic properties are not well constrained due to modelling uncertainties in sub-grid physics, ram pressure stripping, dynamical friction, and tidal forces. In this paper, we use the IllustrisTNG300-1 simulation, with a baryonic resolution of $M_{\rm b} \approx 1.1\times 10^7$ M$_\odot$, to study backsplash galaxies around 1302 isolated galaxy clusters with mass $10^{13.0} < M_{\rm 200,mean} / {\rm M}_\odot< 10^{15.5}$. We employ a decision tree classifier to extract features of galaxies that make them likely to be backsplash galaxies, compared to nearby field galaxies, and find that backsplash galaxies have low gas fractions, high mass-to-light ratios, large stellar sizes, and low black hole occupation fractions. We investigate in detail the origins of these large sizes, and hypothesise their origins are linked to the tidal environments in the cluster. We show that the black hole recentreing scheme employed in many cosmological simulations leads to the loss of black holes from galaxies accreted into clusters, and suggest improvements to these models. Generally, we find that backsplash galaxies are a useful population to test and understand numerical galaxy formation models due to their challenging environments and evolutionary pathways that interact with poorly constrained physics.

preprint2022arXiv

H-alpha emission in local galaxies: star formation, time variability and the diffuse ionized gas

The nebular recombination line H$α$ is widely used as a star-formation rate (SFR) indicator in the local and high-redshift Universe. We present a detailed H$α$ radiative transfer study of high-resolution isolated Milky-Way and Large Magellanic Cloud simulations that include radiative transfer, non-equilibrium thermochemistry, and dust evolution. We focus on the spatial morphology and temporal variability of the H$α$ emission, and its connection to the underlying gas and star formation properties. The H$α$ and H$β$ radial and vertical surface brightness profiles are in excellent agreement with observations of nearby galaxies. We find that the fraction of H$α$ emission from collisional excitation amounts to $f_{\rm col}\sim5-10\%$, only weakly dependent on radius and vertical height, and that scattering boosts the H$α$ luminosity by $\sim40\%$. The dust correction via the Balmer decrement works well (intrinsic H$α$ emission recoverable within $25\%$), though the dust attenuation law depends on the amount of attenuation itself both on spatially resolved and integrated scales. Important for the understanding of the H$α$-SFR connection is the dust and helium absorption of ionizing radiation (Lyman continuum [LyC] photons), which are about $f_{\rm abs}\approx28\%$ and $f_{\rm He}\approx9\%$, respectively. Together with an escape fraction of $f_{\rm esc}\approx6\%$, this reduces the available budget for hydrogen line emission by nearly half ($f_{\rm H}\approx57\%$). We discuss the impact of the diffuse ionized gas, showing - among other things - that the extraplanar H$α$ emission is powered by LyC photons escaping the disc. Future applications of this framework to cosmological (zoom-in) simulations will assist in the interpretation of spectroscopy of high-redshift galaxies with the upcoming James Webb Space Telescope.

preprint2022arXiv

Rapid Convergence of Informed Importance Tempering

Informed Markov chain Monte Carlo (MCMC) methods have been proposed as scalable solutions to Bayesian posterior computation on high-dimensional discrete state spaces, but theoretical results about their convergence behavior in general settings are lacking. In this article, we propose a class of MCMC schemes called informed importance tempering (IIT), which combine importance sampling and informed local proposals, and derive generally applicable spectral gap bounds for IIT estimators. Our theory shows that IIT samplers have remarkable scalability when the target posterior distribution concentrates on a small set. Further, both our theory and numerical experiments demonstrate that the informed proposal should be chosen with caution: the performance of some proposals may be very sensitive to the shape of the target distribution. We find that the &#34;square-root proposal weighting&#34; scheme tends to perform well in most settings.

preprint2022arXiv

Rate-optimal refinement strategies for local approximation MCMC

Many Bayesian inference problems involve target distributions whose density functions are computationally expensive to evaluate. Replacing the target density with a local approximation based on a small number of carefully chosen density evaluations can significantly reduce the computational expense of Markov chain Monte Carlo (MCMC) sampling. Moreover, continual refinement of the local approximation can guarantee asymptotically exact sampling. We devise a new strategy for balancing the decay rate of the bias due to the approximation with that of the MCMC variance. We prove that the error of the resulting local approximation MCMC (LA-MCMC) algorithm decays at roughly the expected $1/\sqrt{T}$ rate, and we demonstrate this rate numerically. We also introduce an algorithmic parameter that guarantees convergence given very weak tail bounds, significantly strengthening previous convergence results. Finally, we apply LA-MCMC to a computationally intensive Bayesian inverse problem arising in groundwater hydrology.

preprint2022arXiv

The THESAN project: predictions for multi-tracer line intensity mapping in the epoch of reionization

Line intensity mapping (LIM) is rapidly emerging as a powerful technique to study galaxy formation and cosmology in the high-redshift Universe. We present LIM estimates of select spectral lines originating from the interstellar medium (ISM) of galaxies and 21 cm emission from neutral hydrogen gas in the Universe using the large volume, high resolution THESAN reionization simulations. A combination of sub-resolution photo-ionization modelling for HII regions and Monte Carlo radiative transfer calculations is employed to estimate the dust-attenuated spectral energy distributions (SEDs) of high-redshift galaxies ($z\gtrsim5.5$). We show that the derived photometric properties such as the ultraviolet (UV) luminosity function and the UV continuum slopes match observationally inferred values, demonstrating the accuracy of the SED modelling. We provide fits to the luminosity--star formation rate relation (L-SFR) for the brightest emission lines and find that important differences exist between the derived scaling relations and the widely used low-$z$ ones because the interstellar medium of reionization era galaxies is generally less metal-enriched than in their low redshift counterparts. We use these relations to construct line intensity maps of nebular emission lines and cross correlate with the 21 cm emission. Interestingly, the wavenumber at which the correlation switches sign ($k_\mathrm{transition}$) depends heavily on the reionization model and to a lesser extent on the targeted emission line, which is consistent with the picture that $k_\mathrm{transition}$ probes the typical sizes of ionized regions. The derived scaling relations and intensity maps represent a timely state-of-the-art framework for forecasting and interpreting results from current and upcoming LIM experiments.

preprint2022arXiv

The THESAN project: properties of the intergalactic medium and its connection to Reionization-era galaxies

The high-redshift intergalactic medium (IGM) and the primeval galaxy population are rapidly becoming the new frontier of extra-galactic astronomy. We investigate the IGM properties and their connection to galaxies at $z\geq5.5$ under different assumptions for the ionizing photon escape and the nature of dark matter, employing our novel THESAN radiation-hydrodynamical simulation suite, designed to provide a comprehensive picture of the emergence of galaxies in a full reionization context. Our simulations have realistic `late&#39; reionization histories, match available constraints on global IGM properties and reproduce the recently-observed rapid evolution of the mean free path of ionizing photons. We additionally examine high-z Lyman-$α$ transmission. The optical depth evolution is consistent with data, and its distribution suggests an even-later reionization than simulated, although with a strong sensitivity to the source model. We show that the effects of these two unknowns can be disentangled by characterising the spectral shape and separation of Lyman-$α$ transmission regions, opening up the possibility to observationally constrain both. For the first time in simulations, THESAN reproduces the modulation of the Lyman-$α$ flux as a function of galaxy distance, demonstrating the power of coupling a realistic galaxy formation model with proper radiation-hydrodynamics. We find this feature to be extremely sensitive on the timing of reionization, while being relatively insensitive to the source model. Overall, THESAN produces a realistic IGM and galaxy population, providing a robust framework for future analysis of the high-z Universe.

preprint2020arXiv

Drift, Minorization, and Hitting Times

The &#34;drift-and-minorization&#34; method, introduced and popularized in (Rosenthal, 1995; Meyn and Tweedie, 1994; Meyn and Tweedie, 2012), remains the most popular approach for bounding the convergence rates of Markov chains used in statistical computation. This approach requires estimates of two quantities: the rate at which a single copy of the Markov chain &#34;drifts&#34; towards a fixed &#34;small set&#34;, and a &#34;minorization condition&#34; which gives the worst-case time for two Markov chains started within the small set to couple with moderately large probability. In this paper, we build on (Oliveira, 2012; Peres and Sousi, 2015) and our work (Anderson, Duanmu, Smith, 2019a; Anderson, Duanmu, Smith, 2019b) to replace the &#34;minorization condition&#34; with an alternative &#34;hitting condition&#34; that is stated in terms of only one Markov chain, and illustrate how this can be used to obtain similar bounds that can be easier to use.

preprint2020arXiv

Fast and memory-optimal dimension reduction using Kac&#39;s walk

In this work, we analyze dimension reduction algorithms based on the Kac walk and discrete variants. (1) For $n$ points in $\mathbb{R}^{d}$, we design an optimal Johnson-Lindenstrauss (JL) transform based on the Kac walk which can be applied to any vector in time $O(d\log{d})$ for essentially the same restriction on $n$ as in the best-known transforms due to Ailon and Liberty [SODA, 2008], and Bamberger and Krahmer [arXiv, 2017]. Our algorithm is memory-optimal, and outperforms existing algorithms in regimes when $n$ is sufficiently large and the distortion parameter is sufficiently small. In particular, this confirms a conjecture of Ailon and Chazelle [STOC, 2006] in a stronger form. (2) The same construction gives a simple transform with optimal Restricted Isometry Property (RIP) which can be applied in time $O(d\log{d})$ for essentially the same range of sparsity as in the best-known such transform due to Ailon and Rauhut [Discrete Comput. Geom., 2014]. (3) We show that by fixing the angle in the Kac walk to be $π/4$ throughout, one obtains optimal JL and RIP transforms with almost the same running time, thereby confirming -- up to a $\log\log{d}$ factor -- a conjecture of Avron, Maymounkov, and Toledo [SIAM J. Sci. Comput., 2010]. Our moment-based analysis of this modification of the Kac walk may also be of independent interest.

preprint2020arXiv

Resonant-line radiative transfer within power-law density profiles

Star-forming regions in galaxies are surrounded by vast reservoirs of gas capable of both emitting and absorbing Lyman-alpha (Lya) radiation. Observations of Lya emitters and spatially extended Lya haloes indeed provide insights into the formation and evolution of galaxies. However, due to the complexity of resonant scattering, only a few analytic solutions are known in the literature. We discuss several idealized but physically motivated scenarios to extend the existing formalism to new analytic solutions, enabling quantitative predictions about the transport and diffusion of Lya photons. This includes a closed form solution for the radiation field and derived quantities including the emergent flux, peak locations, energy density, average internal spectrum, number of scatters, outward force multiplier, trapping time, and characteristic radius. To verify our predictions, we employ a robust gridless Monte Carlo radiative transfer (GMCRT) method, which is straightforward to incorporate into existing ray-tracing codes but requires modifications to opacity-based calculations, including dynamical core-skipping acceleration schemes. We primarily focus on power-law density and emissivity profiles, however both the analytic and numerical methods can be generalized to other cases. Such studies provide additional intuition and understanding regarding the connection between the physical environments and observational signatures of galaxies throughout the Universe.