Researcher profile

Xuejian Shen

Xuejian Shen contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

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

8 published item(s)

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.

preprint2024arXiv

Gas Morphology of Milky Way-like Galaxies in the TNG50 Simulation: Signals of Twisting and Stretching

We present an in-depth analysis of gas morphologies for a sample of 25 Milky Way-like galaxies from the IllustrisTNG TNG50 simulation. We constrain the morphology of cold, warm, hot gas, and gas particles as a whole using a Local Shell Iterative Method (LSIM) and explore its observational implications by computing the hard-to-soft X-ray ratio, which ranges between $10^{-3}$-$10^{-2}$ in the inner $\sim 50 \rm kpc$ of the distribution and $10^{-5}$-$10^{-4}$ at the outer portion of the hot gas distribution. We group galaxies into three main categories: simple, stretched, and twisted. These categories are based on the radial reorientation of the principal axes of the reduced inertia tensor. We find that a vast majority ($77\%$) of the galaxies in our sample exhibit twisting patterns in their radial profiles. Additionally, we present detailed comparisons between 1) the gaseous distributions belonging to individual temperature regimes, 2) the cold gas distributions and stellar distributions, and 3) the gaseous distributions and dark matter (DM) halos. We find a strong correlation between the morphological properties of the cold gas and stellar distributions. Furthermore, we find a correlation between gaseous distributions with DM halo that increases with gas temperature, implying that we may use the warm-hot gaseous morphology as a tracer to probe the DM morphology. Finally, we show gaseous distributions exhibit significantly more prolate morphologies than the stellar distributions and DM halos, which we hypothesize is due to stellar and AGN feedback.

preprint2022arXiv

Dissipative Dark Matter on FIRE: II. Observational signatures and constraints from local dwarf galaxies

We analyze the first set of cosmological baryonic zoom-in simulations of galaxies in dissipative self-interacting dark matter (dSIDM). The simulations utilize the FIRE-2 galaxy formation physics with the inclusion of dissipative dark matter self-interactions modelled as a constant fractional energy dissipation ($f_{\rm diss}=0.5$). In this paper, we examine the properties of dwarf galaxies with $M_{\ast} \sim 10^{5}\operatorname{-}10^{9}\,{\rm M}_{\odot}$ in both isolation and within Milky Way-mass hosts. For isolated dwarfs, we find more compact galaxy sizes and promotion of stellar/neutral gas disk formation in dSIDM with $(σ/m)\leq 1\,{\rm cm^2\,g^{-1}}$ but they are still consistent with observed galaxy sizes and masses. In addition, as a result of the steeper central density profiles developed in dSIDM, the sub-kpc circular velocities of isolated dwarfs in models with $(σ/m)\geq 0.1\,{\rm cm^2\,g^{-1}}$ are enhanced by about a factor of two, which are still consistent with the measured stellar velocity dispersions of Local Group dwarfs but in tension with the HI rotation curves of more massive field dwarfs. Meanwhile, for satellites of the simulated Milky Way-mass hosts, the median circular velocity profiles are marginally affected by dSIDM physics, but dSIDM may help address the missing compact dwarf satellites in CDM. The number of satellites is slightly enhanced in dSIDM, but the differences are small compared with the large host-to-host variations. In conclusion, the dSIDM models with constant cross-section $(σ/m) \gtrsim 0.1\,{\rm cm^2\,g^{-1}}$ (assuming $f_{\rm diss}=0.5$) are effectively ruled out in bright dwarfs ($M_{\rm halo}\sim 10^{11}\,{\rm M}_{\odot}$) by circular velocity constraints. However, models with lower effective cross-sections (at this halo mass/velocity scale) are still viable and can give rise to non-trivial observable signatures.

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.

preprint2020arXiv

High redshift JWST predictions from IllustrisTNG: II. Galaxy line and continuum spectral indices and dust attenuation curves

We present predictions for high redshift ($z=2-10$) galaxy populations based on the IllustrisTNG simulation suite and a full Monte Carlo dust radiative transfer post-processing. Specifically, we discuss the ${\rm H}_α$ and ${\rm H}_β$ + $[\rm O \,III]$ luminosity functions up to $z=8$. The predicted ${\rm H}_β$ + $[\rm O \,III]$ luminosity functions are consistent with present observations at $z\lesssim 3$ with $\lesssim 0.1\,{\rm dex}$ differences in luminosities. However, the predicted ${\rm H}_α$ luminosity function is $\sim 0.3\,{\rm dex}$ dimmer than the observed one at $z\simeq 2$. Furthermore, we explore continuum spectral indices, the Balmer break at $4000$Å(D4000) and the UV continuum slope $β$. The median D4000 versus sSFR relation predicted at $z=2$ is in agreement with the local calibration despite a different distribution pattern of galaxies in this plane. In addition, we reproduce the observed $A_{\rm UV}$ versus $β$ relation and explore its dependence on galaxy stellar mass, providing an explanation for the observed complexity of this relation. We also find a deficiency in heavily attenuated, UV red galaxies in the simulations. Finally, we provide predictions for the dust attenuation curves of galaxies at $z=2-6$ and investigate their dependence on galaxy colors and stellar masses. The attenuation curves are steeper in galaxies at higher redshifts, with bluer colors, or with lower stellar masses. We attribute these predicted trends to dust geometry. Overall, our results are consistent with present observations of high redshift galaxies. Future JWST observations will further test these predictions.

preprint2020arXiv

The Bolometric Quasar Luminosity Function at z = 0-7

In this paper, we provide updated constraints on the bolometric quasar luminosity function (QLF) from $z=0$ to $z=7$. The constraints are based on an observational compilation that includes observations in the rest-frame IR, B band, UV, soft and hard X-ray in past decades. Our method follows Hopkins et al. 2007 with an updated quasar SED model and bolometric and extinction corrections. The new best-fit bolometric quasar luminosity function behaves qualitatively different from the Hopkins et al. 2007 model at high redshift. Compared with the old model, the number density normalization decreases towards higher redshift and the bright-end slope is steeper at $z\gtrsim 2$. Due to the paucity of measurements at the faint end, the faint end slope at $z\gtrsim 5$ is quite uncertain. We present two models, one featuring a progressively steeper faint-end slope at higher redshift and the other featuring a shallow faint-end slope at $z\gtrsim 5$. Further multi-band observations of the faint-end QLF are needed to distinguish between these models. The evolutionary pattern of the bolometric QLF can be interpreted as an early phase likely dominated by the hierarchical assembly of structures and a late phase likely dominated by the quenching of galaxies. We explore the implications of this model on the ionizing photon production by quasars, the CXB spectrum, the SMBH mass density and mass functions. The predicted hydrogen photoionization rate contributed by quasars is subdominant during the epoch of reionization and only becomes important at $z\lesssim 3$. The predicted CXB spectrum, cosmic SMBH mass density and SMBH mass function are generally consistent with existing observations.

preprint2019arXiv

Galaxy Formation with BECDM -- II. Cosmic Filaments and First Galaxies

Bose-Einstein Condensate Dark Matter (BECDM; also known as Fuzzy Dark Matter) is motivated by fundamental physics and has recently received significant attention as a serious alternative to the established Cold Dark Matter (CDM) model. We perform cosmological simulations of BECDM gravitationally coupled to baryons and investigate structure formation at high redshifts ($z \gtrsim 5$) for a boson mass $m=2.5\cdot 10^{-22}~{\rm eV}$, exploring the dynamical effects of its wavelike nature on the cosmic web and the formation of first galaxies. Our BECDM simulations are directly compared to CDM as well as to simulations where the dynamical quantum potential is ignored and only the initial suppression of the power spectrum is considered -- a Warm Dark Matter-like ("WDM") model often used as a proxy for BECDM. Our simulations confirm that "WDM" is a good approximation to BECDM on large cosmological scales even in the presence of the baryonic feedback. Similarities also exist on small scales, with primordial star formation happening both in isolated haloes and continuously along cosmic filaments; the latter effect is not present in CDM. Global star formation and metal enrichment in these first galaxies are delayed in BECDM/"WDM" compared to the CDM case: in BECDM/"WDM" first stars form at $z\sim 13$/$13.5$ while in CDM star formation starts at $z\sim 35$. The signature of BECDM interference, not present in "WDM", is seen in the evolved dark matter power spectrum: although the small scale structure is initially suppressed, power on kpc scales is added at lower redshifts. Our simulations lay the groundwork for realistic simulations of galaxy formation in BECDM.

preprint2019arXiv

High redshift JWST predictions from IllustrisTNG: Dust modelling and galaxy luminosity functions

The James Webb Space Telescop (JWST) promises to revolutionise our understanding of the early Universe, and contrasting its upcoming observations with predictions of the $Λ$CDM model requires detailed theoretical forecasts. Here, we exploit the large dynamic range of the IllustrisTNG simulation suite, TNG50, TNG100, and TNG300, to derive multi-band galaxy luminosity functions from $z=2$ to $z=10$. We put particular emphasis on the exploration of different dust attenuation models to determine galaxy luminosity functions for the rest-frame ultraviolet (UV), and apparent wide NIRCam bands. Our most detailed dust model is based on continuum Monte Carlo radiative transfer calculations employing observationally calibrated dust properties. This calibration results in constraints on the redshift evolution of the dust attenuation normalisation and dust-to-metal ratios yielding a stronger redshift evolution of the attenuation normalisation compared to most previous theoretical studies. Overall we find good agreement between the rest-frame UV luminosity functions and observational data for all redshifts, also beyond the regimes used for the dust-model calibrations. Furthermore, we also recover the observed high redshift ($z=4-6$) UV luminosity versus stellar mass relation, the H$α$ versus star formation rate relation, and the H$α$ luminosity function at $z=2$. The bright end ($M_{\rm UV}>-19.5$) cumulative galaxy number densities are consistent with observational data. For the F200W NIRCam band, we predict that JWST will detect $\sim 80$ ($\sim 200$) galaxies with a signal-to-noise ratio of $10$ ($\sim 5$) within the NIRCam field of view, $2.2\times2.2 \,{\rm arcmin}^{2}$, for a total exposure time of $10^5{\rm s}$ in the redshift range $z=8 \pm 0.5$. These numbers drop to $\sim 10$ ($\sim 40$) for an exposure time of $10^4{\rm s}$.