Source author record

Daniel Anglés-Alcázar

Daniel Anglés-Alcázar appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

16 published item(s)

preprint2026arXiv

CAMELS Environments: The Impact of Local Neighbours on Galaxy Evolution across the SIMBA, IllustrisTNG, ASTRID, and Swift-EAGLE Simulations

Internal feedback from massive stars and active galactic nuclei (AGN) play a key role in galaxy evolution, but external environmental effects can also strongly influence galaxies. We investigate the impact of environment on galaxy evolution, and its dependence on baryonic physics implementation, using Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) spanning a wide range of stellar and AGN feedback implementations in the SIMBA, IllustrisTNG, ASTRID, and Swift-EAGLE galaxy formation models. We show that satellite galaxies are significantly affected by the environment in all simulation models, with their gas fraction and star formation rate (SFR) suppressed in overdense regions compared to similar mass satellites in underdense environments at $z=0$. Central galaxies are less sensitive to environment but tend to show lower gas fraction and SFR in overdense regions at low stellar mass, transitioning to higher gas fraction and SFR for massive galaxies in higher-density environments. Halo baryon fraction ($f_{\rm B}$) and circumgalactic medium mass fraction ($f_{\rm CGM}$) at $z=0$ show clear environmental effects. In SIMBA, low-mass haloes in overdense regions have systematically lower $f_{\rm B}$ and $f_{\rm CGM}$ at fixed halo mass, while Swift-EAGLE haloes in overdense regions have systematically higher $f_{\rm B}$ and $f_{\rm CGM}$ across the full halo mass range, and IllustrisTNG and ASTRID show opposite trends at the low and high mass ends. Environmental effects can flip at higher redshift, with SFR and $f_{\rm B}$ increasing with local density in low-mass haloes before quenching at an increasing overdensity threshold. Our results demonstrate that the impact of environment on galaxy evolution depends significantly on galaxy formation model, and higher-density environments can either suppress or enhance star formation depending on galaxy mass and cosmic epoch.

preprint2026arXiv

Cosmological back-reaction of baryons on dark matter in the CAMELS simulations

Baryonic processes such as radiative cooling and feedback from massive stars and active galactic nuclei (AGN) directly redistribute baryons in the Universe but also indirectly redistribute dark matter due to changes in the gravitational potential. In this work, we investigate this "back-reaction" of baryons on dark matter using thousands of cosmological hydrodynamic simulations from the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project, including parameter variations in the SIMBA, IllustrisTNG, ASTRID, and Swift-EAGLE galaxy formation models. Matching haloes to corresponding N-body (dark matter-only) simulations, we find that virial masses decrease owing to the ejection of baryons by feedback. Relative to N-body simulations, halo profiles show an increased dark matter density in the center (due to radiative cooling) and a decrease in density farther out (due to feedback), with both effects being strongest in SIMBA (> 450% increase at r < 0.01 Rvir). The clustering of dark matter strongly responds to changes in baryonic physics, with dark matter power spectra in some simulations from each model showing as much as 20% suppression or increase in power at k ~ 10 h/Mpc relative to N-body simulations. We find that the dark matter back-reaction depends intrinsically on cosmology (Omega_m and sigma_8) at fixed baryonic physics, and varies strongly with the details of the feedback implementation. These results emphasize the need for marginalizing over uncertainties in baryonic physics to extract cosmological information from weak lensing surveys as well as their potential to constrain feedback models in galaxy evolution.

preprint2022arXiv

Black hole -- galaxy scaling relations in FIRE: the importance of black hole location and mergers

The concurrent growth of supermassive black holes (SMBHs) and their host galaxies remains to be fully explored, especially at high redshift. While often understood as a consequence of self-regulation via AGN feedback, it can also be explained by alternative SMBH accretion models. Here, we expand on previous work by studying the growth of SMBHs with the help of a large suite of cosmological zoom-in simulations (MassiveFIRE) that are part of the Feedback in Realistic Environments (FIRE) project. The growth of SMBHs is modelled in post-processing with different black hole accretion models, placements, and merger treatments, and validated by comparing to on-the-fly calculations. Scaling relations predicted by the gravitational torque driven accretion (GTDA) model agree with observations at low redshift without the need for AGN feedback, in contrast to models in which the accretion rate depends strongly on SMBH mass. At high redshift, we find deviations from the local scaling relations in line with previous theoretical results. In particular, SMBHs are under-massive, presumably due to stellar feedback, but start to grow efficiently once their host galaxies reach $M_* \sim 10^{10} M_{\odot}$. We analyse and explain these findings in the context of a simple analytic model. Finally, we show that the predicted scaling relations depend sensitively on the SMBH location and the efficiency of SMBH merging, particularly in low-mass systems. These findings highlight the relevance of understanding the evolution of SMBH-galaxy scaling relations to predict the rate of gravitational wave signals from SMBH mergers across cosmic history.

preprint2022arXiv

Breaking baryon-cosmology degeneracy with the electron density power spectrum

Uncertain feedback processes in galaxies affect the distribution of matter, currently limiting the power of weak lensing surveys. If we can identify cosmological statistics that are robust against these uncertainties, or constrain these effects by other means, then we can enhance the power of current and upcoming observations from weak lensing surveys such as DES, Euclid, the Rubin Observatory, and the Roman Space Telescope. In this work, we investigate the potential of the electron density auto-power spectrum as a robust probe of cosmology and baryonic feedback. We use a suite of (magneto-)hydrodynamic simulations from the CAMELS project and perform an idealized analysis to forecast statistical uncertainties on a limited set of cosmological and physically-motivated astrophysical parameters. We find that the electron number density auto-correlation, measurable through either kinematic Sunyaev-Zel'dovich observations or through Fast Radio Burst dispersion measures, provides tight constraints on $Ω_{m}$ and the mean baryon fraction in intermediate-mass halos, $\bar{f}_{\mathrm{bar}}$. By obtaining an empirical measure for the associated systematic uncertainties, we find these constraints to be largely robust to differences in baryonic feedback models implemented in hydrodynamic simulations. We further discuss the main caveats associated with our analysis, and point out possible directions for future work.

preprint2022arXiv

Cosmology with one galaxy?

Galaxies can be characterized by many internal properties such as stellar mass, gas metallicity, and star-formation rate. We quantify the amount of cosmological and astrophysical information that the internal properties of individual galaxies and their host dark matter halos contain. We train neural networks using hundreds of thousands of galaxies from 2,000 state-of-the-art hydrodynamic simulations with different cosmologies and astrophysical models of the CAMELS project to perform likelihood-free inference on the value of the cosmological and astrophysical parameters. We find that knowing the internal properties of a single galaxy allow our models to infer the value of $Ω_{\rm m}$, at fixed $Ω_{\rm b}$, with a $\sim10\%$ precision, while no constraint can be placed on $σ_8$. Our results hold for any type of galaxy, central or satellite, massive or dwarf, at all considered redshifts, $z\leq3$, and they incorporate uncertainties in astrophysics as modeled in CAMELS. However, our models are not robust to changes in subgrid physics due to the large intrinsic differences the two considered models imprint on galaxy properties. We find that the stellar mass, stellar metallicity, and maximum circular velocity are among the most important galaxy properties to determine the value of $Ω_{\rm m}$. We believe that our results can be explained taking into account that changes in the value of $Ω_{\rm m}$, or potentially $Ω_{\rm b}/Ω_{\rm m}$, affect the dark matter content of galaxies. That effect leaves a distinct signature in galaxy properties to the one induced by galactic processes. Our results suggest that the low-dimensional manifold hosting galaxy properties provides a tight direct link between cosmology and astrophysics.

preprint2022arXiv

From EMBER to FIRE: predicting high resolution baryon fields from dark matter simulations with Deep Learning

Hydrodynamic simulations provide a powerful, but computationally expensive, approach to study the interplay of dark matter and baryons in cosmological structure formation. Here we introduce the EMulating Baryonic EnRichment (EMBER) Deep Learning framework to predict baryon fields based on dark-matter-only simulations thereby reducing computational cost. EMBER comprises two network architectures, U-Net and Wasserstein Generative Adversarial Networks (WGANs), to predict two-dimensional gas and HI densities from dark matter fields. We design the conditional WGANs as stochastic emulators, such that multiple target fields can be sampled from the same dark matter input. For training we combine cosmological volume and zoom-in hydrodynamical simulations from the Feedback in Realistic Environments (FIRE) project to represent a large range of scales. Our fiducial WGAN model reproduces the gas and HI power spectra within 10% accuracy down to ~10 kpc scales. Furthermore, we investigate the capability of EMBER to predict high resolution baryon fields from low resolution dark matter inputs through upsampling techniques. As a practical application, we use this methodology to emulate high-resolution HI maps for a dark matter simulation of a L=100 Mpc/h comoving cosmological box. The gas content of dark matter haloes and the HI column density distributions predicted by EMBER agree well with results of large volume cosmological simulations and abundance matching models. Our method provides a computationally efficient, stochastic emulator for augmenting dark matter only simulations with physically consistent maps of baryon fields.

preprint2022arXiv

Percent-level constraints on baryonic feedback with spectral distortion measurements

High-significance measurements of the monopole thermal Sunyaev-Zel'dovich CMB spectral distortions have the potential to tightly constrain poorly understood baryonic feedback processes. The sky-averaged Compton-y distortion and its relativistic correction are measures of the total thermal energy in electrons in the observable universe and their mean temperature. We use the CAMELS suite of hydrodynamic simulations to explore possible constraints on parameters describing the subgrid implementation of feedback from active galactic nuclei and supernovae, assuming a PIXIE-like measurement. The small 25 Mpc/h CAMELS boxes present challenges due to the significant cosmic variance. We utilize machine learning to construct interpolators through the noisy simulation data. Using the halo model, we translate the simulation halo mass functions into correction factors to reduce cosmic variance where required. Our results depend on the subgrid model. In the case of IllustrisTNG, we find that the best-determined parameter combination can be measured to ~2% and corresponds to a product of AGN and SN feedback. In the case of SIMBA, the tightest constraint is ~0.2% on a ratio between AGN and SN feedback. A second orthogonal parameter combination can be measured to ~8%. Our results demonstrate the significant constraining power a measurement of the late-time spectral distortion monopoles would have for baryonic feedback models.

preprint2022arXiv

Reionization with Simba: How much does astrophysics matter in modeling cosmic reionization?

Traditional large-scale models of reionization usually employ simple deterministic relations between halo mass and luminosity to predict how reionization proceeds. We here examine the impact on modelling reionization of using more detailed models for the ionizing sources as identified within the $100~{\rm Mpc/h}$ cosmological hydrodynamic simulation {\sc Simba}, coupled with post-processed radiative transfer. Comparing with simple (one-to-one) models, the main difference of using {\sc Simba} sources is the scatter in the relation between dark matter halos and star formation, and hence ionizing emissivity. We find that, at the power spectrum level, the ionization morphology remains mostly unchanged, regardless of the variability in the number of sources or escape fraction. In particular, the power spectrum shape remains unaffected and its amplitude changes slightly by less than 5-10\%, throughout reionization, depending on the scale and neutral fraction. Our results show that simplified models of ionizing sources remain viable to efficiently model the structure of reionization on cosmological scales, although the precise progress of reionization requires accounting for the scatter induced by astrophysical effects.

preprint2022arXiv

Running Late: Testing Delayed Supermassive Black Hole Growth Models Against the Quasar Luminosity Function

Observations of massive galaxies at low redshift have revealed approximately linear scaling relations between the mass of a supermassive black hole (SMBH) and properties of its host galaxy. How these scaling relations evolve with redshift and whether they extend to lower-mass galaxies however remain open questions. Recent galaxy formation simulations predict a delayed, or "two-phase", growth of SMBHs: slow, highly intermittent BH growth due to repeated gas ejection by stellar feedback in low-mass galaxies, followed by more sustained gas accretion that eventually brings BHs onto the local scaling relations. The predicted two-phase growth implies a steep increase, or "kink", in BH-galaxy scaling relations at a stellar mass $M_{*}\sim5\times10^{10} M_{\odot}$. We develop a parametric, semi-analytic model to compare different SMBH growth models against observations of the quasar luminosity function (QLF) at $z\sim0.5-4$. We compare models in which the relation between SMBH mass and galaxy mass is purely linear versus two-phase models. The models are anchored to the observed galaxy stellar mass function, and the BH mass functions at different redshifts are consistently connected by the accretion rates contributing to the QLF. The best fits suggest that two-phase evolution is significantly preferred by the QLF data over a purely linear scaling relation. Moreover, when the model parameters are left free, the two-phase model fits imply a transition mass consistent with that predicted by simulations. Our analysis motivates further observational tests, including measurements of BH masses and AGN activity at the low-mass end, which could more directly test two-phase SMBH growth.

preprint2021arXiv

IQ Collaboratory III: The Empirical Dust Attenuation Framework -- Taking Hydrodynamical Simulations with a Grain of Dust

We present the Empirical Dust Attenuation (EDA) framework -- a flexible prescription for assigning realistic dust attenuation to simulated galaxies based on their physical properties. We use the EDA to forward model synthetic observations for three state-of-the-art large-scale cosmological hydrodynamical simulations: SIMBA, IllustrisTNG, and EAGLE. We then compare the optical and UV color-magnitude relations, $(g-r) - M_r$ and $(FUV-NUV)-M_r$, of the simulations to a $M_r < -20$ and UV complete SDSS galaxy sample using likelihood-free inference. Without dust, none of the simulations match observations, as expected. With the EDA, however, we can reproduce the observed color-magnitude with all three simulations. Furthermore, the attenuation curves predicted by our dust prescription are in good agreement with the observed attenuation-slope relations and attenuation curves of star-forming galaxies. However, the EDA does not predict star-forming galaxies with low $A_V$ since simulated star-forming galaxies are intrinsically much brighter than observations. Additionally, the EDA provides, for the first time, predictions on the attenuation curves of quiescent galaxies, which are challenging to measure observationally. Simulated quiescent galaxies require shallower attenuation curves with lower amplitude than star-forming galaxies. The EDA, combined with forward modeling, provides an effective approach for shedding light on dust in galaxies and probing hydrodynamical simulations. This work also illustrates a major limitation in comparing galaxy formation models: by adjusting dust attenuation, simulations that predict significantly different galaxy populations can reproduce the same UV and optical observations.

preprint2020arXiv

Measuring dynamical masses from gas kinematics in simulated high-redshift galaxies

Advances in instrumentation have recently extended detailed measurements of gas kinematics to large samples of high-redshift galaxies. Relative to most nearby, thin disk galaxies, in which gas rotation accurately traces the gravitational potential, the interstellar medium (ISM) of z>1 galaxies is typically more dynamic and exhibits elevated turbulence. If not properly modeled, these effects can strongly bias dynamical mass measurements. We use high-resolution FIRE-2 cosmological zoom-in simulations to analyze the physical effects that must be considered to correctly infer dynamical masses from gas kinematics. Our analysis covers a range of galaxy properties from low-redshift Milky-Way-mass galaxies to massive high-redshift galaxies (M_* > 10^11 M_sun at z=1). Selecting only snapshots where a disk is present, we calculate the rotational profile v_phi(r) of the cool (10^3.5 K < T < 10^4.5 K) gas and compare it to the circular velocity v_c=sqrt(GM/r). In the simulated galaxies, the gas rotation traces the circular velocity at intermediate radii, but the two quantities diverge significantly in the center and in the outer disk. Our simulations appear to over-predict observed rotational velocities in the centers of massive galaxies (likely from a lack of black hole feedback), so we focus on larger radii. Gradients in the turbulent pressure at these radii can provide additional radial support and bias dynamical mass measurements low by up to 40%. In both the interior and exterior, the gas' motion can be significantly non-circular due to e.g. bars, satellites, and inflows/outflows. We discuss the accuracy of commonly-used analytic models for pressure gradients (or "asymmetric drift") in the ISM of high-redshift galaxies.

preprint2020arXiv

The impact of quenching on galaxy profiles in the Simba simulation

We study specific star formation rate (sSFR) and gas profiles of star forming and green valley galaxies in the Simba cosmological hydrodynamic simulation. Star-forming galaxy half-light radii ($R_{\rm half}$) at $z=0$ and their evolution ($\propto(1+z)^{-0.78}$) agree with observations. Passive galaxy $R_{\rm half}$ agree with observations at high redshift, but by $z=0$ are too large, owing to numerical heating. We compare Simba $z=0$ sSFR radial profiles for star forming and green valley galaxies to observations. Simba shows strong central depressions in star formation rate (SFR), sSFR, and gas fraction in green valley galaxies and massive star-forming systems, qualitatively as observed, owing to black hole X-ray feedback, which pushes central gas outwards. Turning off X-ray feedback leads to centrally peaked sSFR profiles as in other simulations. In conflict with observations, Simba yields green valley galaxies with strongly dropping sSFR profiles beyond $\gtrsim R_{\rm half}$, regardless of AGN feedback. The central depression owes to lowering molecular gas content; the drop in the outskirts owes to reduced star formation efficiency. Simba's satellites have higher central sSFR and lower outskirts sSFR than centrals, in qualitative agreement with observations. At $z=2$ Simba does not show central depressions in massive star-forming galaxies, suggesting Simba's X-ray feedback should be more active at high-$z$. High resolution tests indicate central sSFR suppression is not sensitive to numerical resolution. Reproducing the central sSFR depression in $z=0$ green valley galaxies represents a unique success of Simba. The remaining discrepancies highlight the importance of SFR and gas profiles in constraining quenching mechanisms.

preprint2016arXiv

Gravitational torque-driven black hole growth and feedback in cosmological simulations

We investigate black hole-host galaxy scaling relations in cosmological simulations with a self-consistent black hole growth and feedback model. The sub-grid accretion model captures the key scalings governing angular momentum transport from galactic scales down to parsec scales, while our kinetic feedback implementation enables the injection of outflows with properties chosen to match observed nuclear outflows. We show that "quasar mode" feedback can have a large impact on the thermal properties of the intergalactic medium and the growth of galaxies and massive black holes for kinetic feedback efficiencies as low as 0.1% relative to the bolometric luminosity. Nonetheless, our simulations suggest that the black hole-host scaling relations are only weakly dependent on the effects of black hole feedback on galactic scales, owing to feedback suppressing the growth of galaxies and massive black holes by a similar amount. In contrast, the rate at which gravitational torques feed the central black hole relative to the host galaxy star formation rate governs the slope and normalization of the black hole-host correlations. Our results suggest that a common gas supply regulated by gravitational torques is the primary driver of the observed co-evolution of black holes and galaxies.

preprint2015arXiv

Torque-Limited Growth of Massive Black Holes in Galaxies Across Cosmic Time

We combine cosmological hydrodynamic simulations with analytic models to evaluate the role of galaxy-scale gravitational torques on the evolution of massive black holes at the centers of star-forming galaxies. We confirm and extend our earlier results to show that torque-limited growth yields black holes and host galaxies evolving on average along the Mbh-Mbulge relation from early times down to z = 0 and that convergence onto the scaling relation occurs independent of the initial conditions and with no need for mass averaging through mergers or additional self-regulation processes. Smooth accretion dominates the long-term evolution, with black hole mergers with mass ratios >1:5 representing typically a small fraction of the total growth. Winds from the accretion disk are required to eject significant mass to suppress black hole growth, but there is no need for coupling this wind to galactic-scale gas to regulate black holes in a non-linear feedback loop. Torque-limited growth yields a close-to-linear relation for the star formation rate and the black hole accretion rate averaged over galaxy evolution time scales. However, the SFR-AGN connection has significant scatter owing to strong variability of black hole accretion at all resolved time scales. Eddington ratios can be described by a broad lognormal distribution with median value evolving roughly as (1 + z)^1.9, suggesting a main sequence for black hole growth similar to the cosmic evolution of specific SFRs. Our results offer an attractive scenario consistent with available observations in which cosmological gas infall and transport of angular momentum in the galaxy by gravitational instabilities regulate the long-term co-evolution of black holes and star-forming galaxies.

preprint2014arXiv

Black Hole-Galaxy Correlations without Self-Regulation

Recent models of black hole growth in a cosmological context have forwarded a paradigm in which the growth is self-regulated by feedback from the black hole itself. Here we use cosmological zoom simulations of galaxy formation down to z = 2 to show that such strong self-regulation is required in the popular spherical Bondi accretion model, but that a plausible alternative model in which black hole growth is limited by galaxy-scale torques does not require self-regulation. Instead, this torque-limited accretion model yields black holes and galaxies evolving on average along the observed scaling relations by relying only on a fixed, 5% mass retention rate onto the black hole from the radius at which the accretion flow is fed. Feedback from the black hole may (and likely does) occur, but does not need to couple to galaxy-scale gas in order to regulate black hole growth. We show that this result is insensitive to variations in the initial black hole mass, stellar feedback, or other implementation details. The torque-limited model allows for high accretion rates at very early epochs (unlike the Bondi case), which if viable can help explain the rapid early growth of black holes, while by z = 2 it yields Eddington factors of 1%-10%. This model also yields a less direct correspondence between major merger events and rapid phases of black hole growth. Instead, growth is more closely tied to cosmological disk feeding, which may help explain observational studies showing that, at least at z > 1, active galaxies do not preferentially show merger signatures.

preprint2014arXiv

Cosmological Zoom Simulations of z = 2 Galaxies: The Impact of Galactic Outflows

We use high-resolution cosmological zoom simulations with ~200 pc resolution at z = 2 and various prescriptions for galactic outflows in order to explore the impact of winds on the morphological, dynamical, and structural properties of eight individual galaxies with halo masses ~ 10^11--2x10^12 Msun at z = 2. We present a detailed comparison to spatially and spectrally resolved Hα and other observations of z ~ 2 galaxies. We find that simulations without winds produce massive, compact galaxies with low gas fractions, super-solar metallicities, high bulge fractions, and much of the star formation concentrated within the inner kpc. Strong winds are required to maintain high gas fractions, redistribute star-forming gas over larger scales, and increase the velocity dispersion of simulated galaxies, more in agreement with the large, extended, turbulent disks typical of high-redshift star-forming galaxies. Winds also suppress early star formation to produce high-redshift cosmic star formation efficiencies in better agreement with observations. Sizes, rotation velocities, and velocity dispersions all scale with stellar mass in accord with observations. Our simulations produce a diversity of morphological characteristics - among our three most massive galaxies, we find a quiescent grand-design spiral, a very compact star-forming galaxy, and a clumpy disk undergoing a minor merger; the clumps are evident in Hα but not in the stars. Rotation curves are generally slowly rising, particularly when calculated using azimuthal velocities rather than enclosed mass. Our results are broadly resolution-converged. These results show that cosmological simulations including outflows can produce disk galaxies similar to those observed during the peak epoch of cosmic galaxy growth.