Researcher profile

Yueying Ni

Yueying Ni contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

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

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.

preprint2024arXiv

PRIYA: A New Suite of Lyman-alpha Forest Simulations for Cosmology

We present the PRIYA suite of cosmological simulations, based on the code and hydrodynamic model of the ASTRID simulation, and designed for cosmological analyses of the Lyman-$α$ forest. Our simulation suite spans a $9$-dimensional parameter space, including $4$ cosmological parameters and $5$ astrophysical/thermal parameters. We have run $48$ low fidelity simulations with $1536^3$ particles in a $120$ Mpc/h box and $3$ high fidelity simulations with $3072^3$ particles in a $120$ Mpc/h box. All our simulations include a full physics model for galaxy formation, including supernova and AGN feedback, and thus also contain a realistic population of DLAs. We advance on earlier simulations suites by larger particle loads, by incorporating new physical models for patchy hydrogen and helium reionization, and by self-consistently incorporating a model for AGN feedback. We show that patchy helium reionization imprints an excess in the 1D flux power spectrum on large scales, which may allow future measurements of helium reionization bubble sizes. Simulation parameters are chosen based on a Latin hypercube design and a Gaussian process is used to interpolate to arbitrary parameter combinations. We build a multi-fidelity emulator for the 1D flux power spectrum and the mean IGM temperature. We show that our final interpolation error is $< 1\%$ and that our simulations produce a flux power spectrum converged at the percent level for $z=5.4$ - $2.2$. Our simulation suite will be used to interpret Lyman-$α$ forest 1D flux power spectra from SDSS and future DESI data releases.

preprint2022arXiv

Probing the $z\gtrsim6$ quasars in a universe with IllustrisTNG physics: Impact of gas-based black hole seeding models

We explore implications of a range of black hole (BH) seeding prescriptions on the formation of the brightest $z\gtrsim6$ quasars in cosmological hydrodynamic simulations. The underlying galaxy formation model is the same as in IllustrisTNG. Using constrained initial conditions, we study the growth of BHs in rare overdense regions (forming $\gtrsim10^{12}M_{\odot}/h$ halos by $z=7$) using a $(9~\mathrm{Mpc}/h)^3$ simulated volume. BH growth is maximal within halos that are compact and have a low tidal field. For these halos, we consider an array of gas-based seeding prescriptions wherein $M_{\mathrm{seed}}=10^4-10^6~M_{\odot}/h$ seeds are inserted in halos above critical thresholds for halo mass and dense, metal-poor gas mass (defined as $\tilde{M}_{\mathrm{h}}$ and $\tilde{M}_{\mathrm{sf,mp}}$, respectively, in units of $M_{\mathrm{seed}}$). We find that a seed model with $\tilde{M}_{\mathrm{sf,mp}}=5$ and $\tilde{M}_{\mathrm{h}}=3000$ successfully produces a $z\sim6$ quasar with $\sim10^9~M_{\odot}$ mass and $\sim10^{47}~\mathrm{ergs~s^ {-1}}$ luminosity. BH mergers play a crucial role at $z\gtrsim9$, causing an early boost in BH mass at a time when accretion-driven BH growth is negligible. When more stringent seeding conditions are applied (for e.g., $\tilde{M}_{\mathrm{sf,mp}}=1000$), the relative paucity of BH seeds results in a much lower merger rate. In this case, $z\gtrsim6$ quasars can only be formed if we enhance the maximum allowed BH accretion rates (by factors $\gtrsim10$) compared to the accretion model used in IllustrisTNG. This can be achieved either by allowing for super-Eddington accretion, or by reducing the radiative efficiency. Our results show that progenitors of $z\sim6$ quasars have distinct BH merger histories for different seeding models, which will be distinguishable with LISA observations.

preprint2022arXiv

Robust field-level inference with dark matter halos

We train graph neural networks on halo catalogues from Gadget N-body simulations to perform field-level likelihood-free inference of cosmological parameters. The catalogues contain $\lesssim$5,000 halos with masses $\gtrsim 10^{10}~h^{-1}M_\odot$ in a periodic volume of $(25~h^{-1}{\rm Mpc})^3$; every halo in the catalogue is characterized by several properties such as position, mass, velocity, concentration, and maximum circular velocity. Our models, built to be permutationally, translationally, and rotationally invariant, do not impose a minimum scale on which to extract information and are able to infer the values of $Ω_{\rm m}$ and $σ_8$ with a mean relative error of $\sim6\%$, when using positions plus velocities and positions plus masses, respectively. More importantly, we find that our models are very robust: they can infer the value of $Ω_{\rm m}$ and $σ_8$ when tested using halo catalogues from thousands of N-body simulations run with five different N-body codes: Abacus, CUBEP$^3$M, Enzo, PKDGrav3, and Ramses. Surprisingly, the model trained to infer $Ω_{\rm m}$ also works when tested on thousands of state-of-the-art CAMELS hydrodynamic simulations run with four different codes and subgrid physics implementations. Using halo properties such as concentration and maximum circular velocity allow our models to extract more information, at the expense of breaking the robustness of the models. This may happen because the different N-body codes are not converged on the relevant scales corresponding to these parameters.

preprint2022arXiv

The ASTRID Simulation: Galaxy Formation and Reionization

We introduce the ASTRID simulation, a large-scale cosmological hydrodynamic simulation in a $250$ Mpc/h box with $2\times 5500^3$ particles. ASTRID contains a large number of high redshift galaxies, which can be compared to future survey data, and resolves galaxies in halos more massive than $2\times 10^9 M_\odot$. ASTRID has been run from $z=99$ to $z=3$. As a particular focus is modelling the high redshift Universe, it contains models for inhomogeneous hydrogen and helium reionization, baryon relative velocities and massive neutrinos, as well as supernova and AGN feedback. The black hole model includes mergers driven by dynamical friction rather than repositioning. We briefly summarise the implemented models, and the technical choices we took when developing the simulation code. We validate the model, showing good agreement with observed UV luminosity functions, galaxy stellar mass functions and specific star formation rates. We show that the redshift at which a given galaxy underwent hydrogen reionization has a large effect on the halo gas fraction. Finally, at $z=6$, halos with $M \sim 2\times 10^9 M_\odot$ which have been reionized have a star formation rate $1.5$ times greater than those which have not yet been reionized.

preprint2022arXiv

The ASTRID simulation: the evolution of Supermassive Black Holes

We present the evolution of black holes (BHs) and their relationship with their host galaxies in Astrid, a large-volume cosmological hydrodynamical simulation with box size 250 $h^{-1} \rm Mpc$ containing $2\times5500^3$ particles evolved to z=3. Astrid statistically models BH gas accretion and AGN feedback to their environments, applies a power-law distribution for BH seed mass $M_{\rm sd}$, uses a dynamical friction model for BH dynamics and executes a physical treatment of BH mergers. The BH population is broadly consistent with empirical constraints on the BH mass function, the bright end of the luminosity functions, and the time evolution of BH mass and accretion rate density. The BH mass and accretion exhibit a tight correlation with host stellar mass and star formation rate. We trace BHs seeded before z>10 down to z=3, finding that BHs carry virtually no imprint of the initial $M_{\rm sd}$ except those with the smallest $M_{\rm sd}$, where less than 50\% of them have doubled in mass. Gas accretion is the dominant channel for BH growth compared to BH mergers. With dynamical friction, Astrid predicts a significant delay for BH mergers after the first encounter of a BH pair, with a typical elapse time of about 200 Myrs. There are in total $4.5 \times 10^5$ BH mergers in Astrid at z>3, $\sim 10^3$ of which have X-ray detectable EM counterparts: a bright kpc scale dual AGN with $L_X>10^{43}$ erg/s. BHs with $M_{\rm BH} \sim 10^{7-8} M_{\odot}$ experience the most frequent mergers. Galaxies that host BH mergers are unbiased tracers of the overall $M_{\rm BH} - M_{*}$ relation. Massive ($>10^{11} M_{\odot}$) galaxies have a high occupation number (>10) of BHs, and hence host the majority of BH mergers.

preprint2022arXiv

The BlueTides Mock Image Catalogue: Simulated observations of high-redshift galaxies and predictions for JWST imaging surveys

We present a mock image catalogue of ~100,000 MUV=-22.5 to -19.6 mag galaxies at z=7-12 from the BlueTides cosmological simulation. We create mock images of each galaxy with the James Webb (JWST), Hubble, Roman, and Euclid Space Telescopes, as well as Subaru, and VISTA, with a range of near- and mid-infrared filters. We perform photometry on the mock images to estimate the success of these instruments for detecting high-z galaxies. We predict that JWST will have unprecedented power in detecting high-z galaxies, with a 95% completeness limit at least 2.5 magnitudes fainter than VISTA and Subaru, 1.1 magnitudes fainter than Hubble, and 0.9 magnitudes fainter than Roman, for the same wavelength and exposure time. Focusing on JWST, we consider a range of exposure times and filters, and find that the NIRCam F356W and F277W filters will detect the faintest galaxies, with 95% completeness at m=27.4 mag in 10ks exposures. We also predict the number of high-z galaxies that will be discovered by upcoming JWST imaging surveys. We predict that the COSMOS-Web survey will detect ~1000 MUV<-20.1 mag galaxies at 6.5<z<7.5, by virtue of its large survey area. JADES-Medium will detect almost 100% of MUV<-20 mag galaxies at z<8.5 due to its significant depth, however with its smaller survey area it will detect only ~100 of these galaxies at 6.5<z<7.5. Cosmic variance results in a large range in the number of predicted galaxies each survey will detect, which is more evident in smaller surveys such as CEERS and the PEARLS NEP and GOODS-S fields.

preprint2022arXiv

The Impact of Dust on the Sizes of Galaxies in the Epoch of Reionization

We study the sizes of galaxies in the Epoch of Reionization using a sample of ~100,000 galaxies from the BlueTides cosmological hydrodynamical simulation from z=7 to 11. We measure the galaxy sizes from stellar mass and luminosity maps, defining the effective radius as the minimum radius which could enclose the pixels containing 50% of the total mass/light in the image. We find an inverse relationship between stellar mass and effective half-mass radius, suggesting that the most massive galaxies are more compact and dense than lower mass galaxies, which have flatter mass distributions. We find a mildly negative relation between intrinsic far-ultraviolet luminosity and size, while we find a positive size-luminosity relation when measured from dust-attenuated images. This suggests that dust is the predominant cause of the observed positive size-luminosity relation, with dust preferentially attenuating bright sight lines resulting in a flatter emission profile and thus larger measured effective radii. We study the size-luminosity relation across the rest-frame ultraviolet and optical, and find that the slope decreases at longer wavelengths; this is a consequence of the relation being caused by dust, which produces less attenuation at longer wavelengths. We find that the far-ultraviolet size-luminosity relation shows mild evolution from z=7 to 11, and galaxy size evolves with redshift as $R\propto(1+z)^{-m}$, where $m=0.662\pm0.009$. Finally, we investigate the sizes of z=7 quasar host galaxies, and find that while the intrinsic sizes of quasar hosts are small relative to the overall galaxy sample, they have comparable sizes when measured from dust-attenuated images.

preprint2021arXiv

Massive Black Hole Mergers with Orbital Information: Predictions from the ASTRID Simulation

We examine massive black hole (MBH) mergers and their associated gravitational wave signals from the large-volume cosmological simulation Astrid. Astrid includes galaxy formation and black hole models recently updated with a MBH seed population between $3\times 10^4M_{\odot}/h$ and $3\times 10^5M_{\odot}/h$ and a sub-grid dynamical friction (DF) model to follow the MBH dynamics down to $1.5\;\text{ckpc}/h$. We calculate initial eccentricities of MBH orbits directly from the simulation at kpc-scales, and find orbital eccentricities above $0.7$ for most MBH pairs before the numerical merger. After approximating unresolved evolution on scales below ${\sim 200\,\text{pc}}$, we find that the in-simulation DF on large scales accounts for more than half of the total orbital decay time ($\sim 500\,\text{Myrs}$) due to DF. The binary hardening time is an order of magnitude longer than the DF time, especially for the seed-mass binaries ($M_\text{BH}<2M_\text{seed}$). As a result, only $\lesssim20\%$ of seed MBH pairs merge at $z>3$ after considering both unresolved DF evolution and binary hardening. These $z>3$ seed-mass mergers are hosted in a biased population of galaxies with the highest stellar masses of $>10^9\,M_\odot$. With the higher initial eccentricity prediction from Astrid, we estimate an expected merger rate of $0.3-0.7$ per year from the $z>3$ MBH population. This is a factor of $\sim 7$ higher than the prediction using the circular orbit assumption. The LISA events are expected at a similar rate, and comprise $\gtrsim 60\%$ seed-seed mergers, $\sim 30\%$ involving only one seed-mass MBH, and $\sim 10\%$ mergers of non-seed MBHs.

preprint2020arXiv

The early growth of supermassive black holes in cosmological hydrodynamic simulations with constrained Gaussian realizations

The paper examines the early growth of supermassive black holes (SMBHs) in cosmological hydrodynamic simulations with different BH seeding scenarios. Employing the constrained Gaussian realization, we reconstruct the initial conditions in the large-volume BlueTides simulation and run them to $z=6$ to cross-validate that the method reproduces the first quasars and their environments. Our constrained simulations in a volume of $(15\, h^{-1}{\rm Mpc})^3$ successfully recover the evolution of large-scale structure and the stellar and BH masses in the vicinity of a $\sim10^{12}\, M_{\odot}$ halo which we identified in BlueTides at $z\sim7$ hosting a $\sim10^9\, M_{\odot}$ SMBH. Among our constrained simulations, only the ones with a low-tidal field and high-density peak in the initial conditions induce the fastest BH growth required to explain the $z>6$ quasars. We run two sets of simulations with different BH seed masses of $5\times10^3$, $5\times10^4$, and $5\times10^5\, h^{-1}M_{\odot}$, (a) with the same ratio of halo to BH seed mass and (b) with the same halo threshold mass. At $z=6$, all the SMBHs converge in mass to $\sim10^9\, M_{\odot}$ except for the one with the smallest seed in (b) undergoing critical BH growth and reaching $10^8$ -- $10^9\, M_{\odot}$, albeit with most of the growth in (b) delayed compared to set (a). The finding of eight BH mergers in the small-seed scenario (four with masses $10^4$ -- $10^6\, M_{\odot}$ at $z>12$), six in the intermediate-seed scenario, and zero in the large-seed scenario suggests that the vast BHs in the small-seed scenario merge frequently during the early phases of the growth of SMBHs. The increased BH merger rate for the low-mass BH seed and halo threshold scenario provides an exciting prospect for discriminating BH formation mechanisms with the advent of multi-messenger astrophysics and next-generation gravitational wave facilities.

preprint2019arXiv

QSO obscuration at high redshift ($z \gtrsim 7$): Predictions from the BlueTides simulation

High-$z$ AGNs hosted in gas rich galaxies are expected to grow through significantly obscured accretion phases. This may limit or bias their observability. In this work, we use \textsc{BlueTides}, a large volume cosmological simulation of galaxy formation to examine quasar obscuration for the highest-redshift ($z \geq 7$) supermassive black holes residing in the center of galaxies. We find that for the bright quasars, most of the high column density gas ($>90\%$) resides in the innermost regions of the host galaxy, (typically within $< 10$ ckpc), while the gas in the outskirts is a minor contributor to the $N_\mathrm H$. The brightest quasars can have large angular variations in galactic obscuration, over 2 orders of magnitude, where the lines of sight with the lowest obscuration are those formed via strong gas outflows driven by AGN feedback. We find that for the overall AGN population, the mean $N_\mathrm H$ is generally larger for high luminosity and BH mass, while the $N_\mathrm H$ distribution is significantly broadened, developing a low $N_\mathrm H $ wing due to the angular variations driven by the AGN outflows/feedback. The obscured fraction P($N_{\rm H} > 10^{23} {\rm cm}^{-2}$) typically range from 0.6 to 1.0 for increasing $L_{X}$ (with $L_X > 10^{43} \rm{ergs/s}$), with no clear trend of redshift evolution. With respect to the galaxy host property, we find a linear relation between $N_{\rm H}$, $M_*$ and $M_{\rm H_2}$ with $\log N_{\rm H} = (0.24 \pm 0.03) \log M_{*} + (20.7 \pm 0.3)$ and $\log N_{\rm H} = (0.47 \pm 0.03) \log M_{\rm H_2} + (18.4 \pm 0.3)$. The dust optical depth in the UV band $τ_{\mathrm UV}$ has tight positive correlation with $N_{\rm H}$. Our dust extincted UVLF is about 1.5 dex lower than the intrinsic UVLF, implying that more than 99\% of the $z \sim 7$ AGNs are heavily dust extincted and therefore would be missed by the UV band observation.