BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Renyue Cen

Renyue Cen contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.GAastro-ph.COastro-ph.HEastro-ph.IMastro-ph.SR

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
16works
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

16 published item(s)

preprint2025arXiv

Rapid formation of a very massive star >50000 $M_\odot$ and subsequently an IMBH from runaway collisions. Direct N-body and Monte Carlo simulations of dense star clusters

Context. We present simulations of a massive young star cluster using \textsc{Nbody6++GPU} and \textsc{MOCCA}. The cluster is initially more compact than previously published models, with one million stars, a total mass of $5.86 \times 10^5~\mathrm{M}_{\odot}$, and a half-mass radius of $0.1~\mathrm{pc}$. Aims. We analyse the formation and growth of a very massive star (VMS) through successive stellar collisions and investigate the subsequent formation of an intermediate-mass black hole (IMBH) in the core of a dense star cluster. Methods. We use both direct \textit{N}-body and Monte Carlo simulations, incorporating updated stellar evolution prescriptions (SSE/BSE) tailored to massive stars and VMSs. These include revised treatments of stellar radii, rejuvenation, and mass loss during collisions. While the prescriptions represent reasonable extrapolations into the VMS regime, the internal structure and thermal state of VMSs formed through stellar collisions remain uncertain, and future work may require further refinement. Results. We find that runaway stellar collisions in the cluster core produce a VMS exceeding $5 \times 10^4~\mathrm{M}_{\odot}$ within 5 Myr, which subsequently collapses into an IMBH. Conclusions. Our model suggests that dense stellar environments may enable the formation of very massive stars and massive black hole seeds through runaway stellar collisions. These results provide a potential pathway for early black hole growth in star clusters and offer theoretical context for interpreting recent JWST observations of young, compact clusters at high redshift.

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

Resolving the Origins and Pathways of Ionizing Radiation Escape with UV Integral Field Spectroscopy

The Epoch of Reionization marks the last major phase transition in the early Universe, during which the majority of neutral hydrogen once filling the intergalactic medium was ionized by the first galaxies. The James Webb Space Telescope (JWST) is now identifying promising galaxy candidates capable of producing sufficient ionizing photons to drive this transformation. However, the fraction of these photons that escape into intergalactic space--the escape fraction--remains highly uncertain. Stellar feedback is thought to play a critical role in carving low-density channels that allow ionizing radiation to escape, but the dominant mechanisms, their operation, and their connection to observable signatures are not well understood. Local analogs of high-redshift galaxies offer a powerful alternative for studying these processes, since ionizing radiation is unobservable at high redshift due to intergalactic absorption. However, current UV space-based instrumentation lacks the spatial resolution and sensitivity required to fully address this problem. The core challenge lies in the multiscale nature of LyC escape: ionizing photons are generated on scales of 1--100 pc in super star clusters but must traverse the circumgalactic medium which can extend beyond 100 kpc. The proposed Habitable Worlds Observatory (HWO) will provide a platform for future UV instruments capable of resolving these scales. In this article, we present a science case for understanding how LyC photons escape from star-forming galaxies and define the observational requirements for future instruments aboard HWO, including a UV integral field spectrograph (IFS).

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Detecting preheating in proto-clusters with Lyman-$α$ Forest Tomography

Studies of low redshift galaxy clusters suggest the intra-cluster medium (ICM) has experienced non-gravitational heating during the formation phase of the clusters. Using simple phenomenological heating prescriptions, we simulate the effect of this preheating of the nascent ICM in galaxy proto-clusters and examine its effect on Lyman-$α$ (Ly$α$) forest tomographic maps. We analyse a series of cosmological zoom-in simulations of proto-clusters within the framework of the Ly$α$ transmission-dark matter (DM) density distribution. We find that the more energy is injected into the proto-ICM at $z$ = 3, the more the distribution at high DM density tilts towards higher Ly$α$ transmission. This effect has been confirmed in both low-resolution simulations adopting a preheating scheme based on entropy floors, as well as in higher-resolution simulations with another scheme based on energy floors. The evolution of the slope of this distribution is shown to vary with redshift. The methodology developed here can be applied to current and upcoming Ly$α$ forest tomographic survey data to help constrain feedback models in galaxy proto-clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Osaka Feedback Model II: Modeling Supernova Feedback Based on High-Resolution Simulations

Feedback from supernovae (SNe) is an essential mechanism that self-regulates the growth of galaxies, and a better model of SN feedback is still needed in galaxy formation simulations. In the first part of this paper, using an Eulerian hydrodynamic code Athena++, we find universal scaling relations for the time evolution of momentum and radius for a superbubble, when the momentum and time are scaled by those at the shell-formation time. In the second part of this paper, we develop an SN feedback model based on the Athena++ simulation results utilizing Voronoi tessellation around each star particle, and implement it into the GADGET3-Osaka smoothed particle hydrodynamic code. Our feedback model was demonstrated to be isotropic and conservative in terms of energy and momentum. We examined the mass/energy/metal loading factors and find that our stochastic thermal feedback model produced galactic outflow that carries metals high above the galactic plane but with weak suppression of star formation. Additional mechanical feedback further suppressed star formation and brought the simulation results in better agreement with the observations of the Kennicutt--Schmidt relation, with all the results being within the uncertainties of observed data. We argue that both thermal and mechanical feedback are necessary for the SN feedback model of galaxy evolution when an individual SN bubble is unresolved.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Pilot-WINGS: An extended MUSE view of the structure of Abell 370

We investigate the strong-lensing cluster Abell 370 (A370) using a wide Integral Field Unit (IFU) spectroscopic mosaic from the Multi-Unit Spectroscopic Explorer (MUSE). IFU spectroscopy provides significant insight into the structure and mass content of galaxy clusters, yet IFU-based cluster studies focus almost exclusively on the central Einstein-radius region. Covering over 14 arcmin$^2$, the new MUSE mosaic extends significantly beyond the A370 Einstein radius, providing, for the first time, a detailed look at the cluster outskirts. Combining these data with wide-field, multi-band Hubble Space Telescope (HST) imaging from the BUFFALO project, we analyse the distribution of objects within the cluster and along the line of sight. Identifying 416 cluster galaxies, we use kinematics to trace the radial mass profile of the halo, providing a mass estimate independent from the lens model. We also measure radially-averaged properties of the cluster members, tracking their evolution as a function of infall. Thanks to the high spatial resolution of our data, we identify six cluster members acting as galaxy-galaxy lenses, which constrain localized mass distributions beyond the Einstein radius. Finally, taking advantage of MUSE's 3D capabilities, we detect and analyse multiple spatially extended overdensities outside of the cluster that influence lensing-derived halo mass estimates. We stress that much of this work is only possible thanks to the robust, extended IFU coverage, highlighting its importance even in less optically dense cluster regions. Overall, this work showcases the power of combining HST+MUSE, and serves as the initial step towards a larger and wider program targeting several clusters.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Absorption line abundances in the SMC-like galaxy UGC 5282: evidence for ISM dilution from inflows on kiloparsec scales

We present an HST Cosmic Origins Spectrograph (COS) spectrum of the QSO SDSS J095109.12+330745.8 ($z_{em} = 0.645$) whose sightline passes through the SMC-like dwarf galaxy UGC 5282 ($M_B = -16.0$, $cz = 1577$ km s$^{-1}$ ), 1.2 kpc in projection from the central HII region of the galaxy. Damped Ly$α$ (DLA) absorption is detected at the redshift of UGC 5282 with $\log$[$N$(HI) cm$^{-2}\:$] = 20.89$^{+0.12}_{-0.21}\:\:$. Analysis of the accompanying SII, PII and OI metal lines yields a neutral gas metallicity, $Z_{HI}$, of [S/H] $\simeq$ [P/H] $= -0.80\pm0.24$. The metallicity of ionized gas from the central HII region, $Z_{HII}\:$, measured from its emission lines is [O/H]$= -0.37\pm0.10$, a difference of $+0.43\pm 0.26$ from $Z_{HI}\:$. This difference $δ$ is consistent with that seen towards HII regions in other star-forming galaxies and supports the idea that ionized gas near star forming regions shows systematically higher metallicities than exist in the rest of a galaxy&#39;s neutral interstellar medium (ISM). The positive values of $δ$ found in UGC 5282 (and the other star forming galaxies) is likely due to infalling low metallicity gas from the intergalactic medium that mixes with the galaxy&#39;s ISM on kpc scales. This model is also consistent with broad Ly$α$ emission detected at the bottom of the DLA absorption, offset by $\sim 125$ km s$^{-1}$ from the absorption velocity. Models of galaxy evolution that attempt to replicate population characteristics, such as the mass-metallicity relation, may need to start with a galaxy metallicity represented by $Z_{HI}$ rather than that measured traditionally from $Z_{HII}\:$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Assembly Conformity of Structure Growth: Fossil versus Normal Groups of Galaxies

Using a semi-analytic method calibrated to the global star formation history and the stellar mass function at $z=0$, we attempt to understand the most stellar deficient galaxy groups. We argue such groups are a kind of fossil group (FGs) -- in comparison to the normal groups of galaxies, they assemble both halo and stellar mass earlier. We find there is a central galaxy and satellite conformity between these FGs and normal groups: centrals and satellites in the former form earlier and more stellar deficient than their counterparts of the latter. We term this effect &#34;Assembly Conformity&#34; of dark matter halos. This effect accounts for about 70 percent of the difference in stellar content between FGs and normal groups. When split by the peak redshift for the star formation rate of a group, the mass functions of satellite halos on either side of the peak redshift are found to be indistinguishable between FGs and normal groups, indicating a self-similarity of halo assembly with respect to the peak. The &#34;baryonic environmental&#34; effect due to ram-pressure and gas heating accounts for about 30 percent of the difference in stellar content. While the total stellar mass of FGs is lower than that of normal groups, we predict that the mass of the brightest central galaxy of FGs is, on average, higher than that of normal groups. We also predict that in the central galaxies of FGs, there is a negative stellar age gradient from the center outward, where the opposite is expected for those in normal groups.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Circumnuclear Molecular Gas in Low-redshift Quasars and Matched Star-forming Galaxies

A series of gravitational instabilities in a circumnuclear gas disk (CND) are required to trigger gas transport to a central supermassive black hole (SMBH) and ignite Active Galactic Nuclei (AGNs). A test of this scenario is to investigate whether an enhanced molecular gas mass surface density ($Σ_{\rm mol}$) is found in the CND-scale of quasars relative to a comparison sample of inactive galaxies. Here we performed sub-kpc resolution CO(2-1) observations with ALMA of four low-redshift ($z \sim 0.06$), luminous ($\sim 10^{45}$ erg s$^{-1}$) quasars with each matched to a different star-forming galaxy, having similar redshift, stellar mass, and star-formation rate. We detected CO(2-1) emission from all quasars, which show diverse morphologies. Contrary to expectations, $Σ_{\rm mol}$ of the quasar sample, computed from the CO(2-1) luminosity, tends to be smaller than the comparison sample at $r < 500$ pc; there is no systematic enhancement of $Σ_{\rm mol}$ in our quasars. We discuss four possible scenarios that would explain the lower molecular gas content (or CO(2-1) luminosity as an actual observable) at the CND-scale of quasars, i.e., AGN-driven outflows, gas-rich minor mergers, time-delay between the onsets of a starburst-phase and a quasar-phase, and X-ray-dominated region (XDR) effects on the gas chemical abundance and excitation. While not extensively discussed in the literature, XDR effects can have an impact on molecular mass measurements particularly in the vicinity of luminous quasar nuclei; therefore higher resolution molecular gas observations, which are now viable using ALMA, need to be considered.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Evidence for Infalling Gas in a Lyman-$α$ Blob

Lyman-$α$ blobs (LABs) are spatially extended nebulae of emission in the Ly$α$ line of hydrogen, seen at high redshifts$^{1,2}$, and most commonly found in the dense environment of star-forming galaxies$^{3,4}$. The origin of Ly$α$ emission in the LABs is still unclear and under debate$^{5}$. Proposed powering sources generally fall into two categories: (1) photoionization, galactic super-winds/outflows, resonant scattering of Ly$α$ photons from starbursts or active galactic nuclei (AGNs)$^{6,7,8,9,10}$ and (2) cooling radiation from cold streams of gas accreting onto galaxies$^{12}$. Here we analyze the gas kinematics within a LAB providing rare observational evidence for infalling gas. This is consistent with the release of gravitational accretion energy as cold streams radiate Ly$α$ photons. It also provides direct evidence for possible cold streams feeding the central galaxies. The infalling gas is not important by mass but hints at more than one mechanism to explain the origin of the extended Ly$α$ emission around young galaxies. It is also possible that the infalling gas may represent material falling back to the galaxy from where it originated, forming a galactic fountain.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Non-steady heating of cool cores of galaxy clusters by ubiquitous turbulence and AGN

Recent cosmological simulations have shown that turbulence should be generally prevailing in clusters because clusters are continuously growing through matter accretion. Using one-dimensional hydrodynamic simulations, we study the heating of cool-core clusters by the ubiquitous turbulence as well as feedback from the central active galactic nuclei (AGNs) for a wide range of cluster and turbulence parameters, focusing on the global stability of the core. We find that the AGN shows intermittent activities in the presence of moderate turbulence similar to the one observed with Hitomi. The cluster core maintains a quasi-equilibrium state for most of the time because the heating through turbulent diffusion is nearly balanced with radiative cooling. The balance is gradually lost because of slight dominance of the radiative cooling, and the AGN is ignited by increased gas inflow. Finally, when the AGN bursts, the core is heated almost instantaneously. Thanks to the pre-existing turbulence, the heated gas is distributed throughout the core without becoming globally unstable and causing catastrophic cooling, and the core recovers the quasi-equilibrium state. The AGN bursts can be stronger in lower-mass clusters. Predictions of our model can be easily checked with future X-ray missions like XRISM and Athena.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Physics of Non-Universal Larson&#39;s Relation

From a new perspective, we re-examine self-gravity and turbulence jointly, in hopes of understanding the physical basis for one of the most important empirical relations governing clouds in the interstellar medium (ISM), the Larson&#39;s Relation relating velocity dispersion ($σ_R$) to cloud size ($R$). We report on two key new findings. First, the correct form of the Larson&#39;s Relation is $σ_R=α_v^{1/5}σ_{pc}(R/1pc)^{3/5}$, where $α_v$ is the virial parameter of clouds and $σ_{pc}$ is the strength of the turbulence, if the turbulence has the Kolmogorov spectrum. Second, the amplitude of the Larson&#39;s Relation, $σ_{pc}$, is not universal, differing by a factor of about two between clouds on the Galactic disk and those at the Galactic center, evidenced by observational data.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Physics of Prodigious Lyman Continuum Leakers

An analysis of the dynamics of a star formation event is performed. It is shown that galaxies able to drive leftover gas to sufficient altitudes in a few million years are characterized by two basic properties: small sizes (<1kpc) and high star formation rate surface densities (Sigma_SFR > 10 Msun/yr/kpc2). For the parameter space of relevance, the outflow is primarily driven by supernovae with radiation pressure being significant but subdominant. Our analysis provides the unifying physical origin for a diverse set of observed LyC leakers, including the green-peas galaxies, [SII]-weak galaxies, Lyman-alpha emitters, with these two characteristics as the common denominator. Among verifiable physical properties of LyC leakers, we predict that (1) the newly formed stellar masses are are typically in the range of 1e8-1e10 Msun, except perhaps ULIRGs, (2) the outflow velocities are typically in the range typically of 100-600km/s, but may exceed 1e3 km/s in ULIRGs, with a strong positive correlation between the stellar masses formed and the outflow velocities, (3) the overall escape fraction of galaxies is expected to increase with increasing redshift, given the cosmological trend that galaxies become denser and more compact with increasing redshift. In addition, two interesting by-product predictions are also borne out. First, ULIRGs appear to be in a parameter region where they should be prodigious LyC leakers, unless there is a large ram-pressure. Second, Lyman break galaxies (LBGs) are not supposed to be prodigious LyC leakers in our model, given their claimed effective radii exceeding 1kpc.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The BUFFALO HST Survey

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Where do quasar hosts lie with respect to the size-mass relation of galaxies?

The evolution of the galaxy size - stellar mass (Mstellar) relation has been a puzzle for over a decade. High redshift galaxies are significantly more compact than galaxies observed today, at an equivalent mass, but how much of this apparent growth is driven by progenitor bias, minor mergers, secular processes, or feedback from AGN is unclear. To help disentangle the physical mechanisms at work by addressing the latter, we study the galaxy size - Mstellar relation of 32 carefully-selected broad-line AGN hosts at 1.2 < z < 1.7 (7.5 < log M_BH < 8.5; L_bol/L_Edd > 0.1). Using HST with multi-band photometry and state-of-the-art modeling techniques, we measure half-light radii while accounting for uncertainties from subtracting bright central point sources. We find AGN hosts to have sizes ranging from 1 to 6 kpc at Mstellar ~ 0.3 - 1 x 10^11 Msun. Thus, many hosts have intermediate sizes as compared to equal-mass star-forming and quiescent galaxies. While inconsistent with the idea that AGN feedback may induce an increase in galaxy sizes, this finding is consistent with hypotheses in which AGNs preferentially occur in systems with prior concentrated gas reservoirs, or are involved in secular compaction processes perhaps responsible for simultaneously building bulges and shutting down star formation. If driven by minor mergers, which do not grow central black holes as fast as they do bulge-like stellar structures, such a process would explain both the galaxy size - mass relation observed here and the evolution in the black hole, bulge mass relation described in a companion paper.

0 citations0 reviewsNo code linkedOpen access