Showing 1–50 of 282 results for author: Cen, R

Cloud Crushing and Dissipation of Uniformly-Driven Adiabatic Turbulence in Circumgalactic Media

Authors: Alex Lv, Lile Wang, Renyue Cen, Luis C. Ho

Abstract: The circumgalactic medium (CGM) is responsive to kinetic disruptions generated by nearby astrophysical events. In this work, we study the saturation and dissipation of turbulent hydrodynamics within the CGM through an extensive array of 252 numerical simulations with a large parameter space. These simulations are endowed with proper cooling mechanisms to consistently explore the parameter space sp… ▽ More The circumgalactic medium (CGM) is responsive to kinetic disruptions generated by nearby astrophysical events. In this work, we study the saturation and dissipation of turbulent hydrodynamics within the CGM through an extensive array of 252 numerical simulations with a large parameter space. These simulations are endowed with proper cooling mechanisms to consistently explore the parameter space spanned by the average gas density, metallicity, and turbulence driving strength. A dichotomy emerges in the dynamics dissipation behaviors. Disturbances that are hot and subsonic are characterized by weak compression and slow dissipation, resulting in density fluctuations typically $\lesssim 10^{-2}$. Conversely, warm supersonic turbulence, marked by significant compression shocks and subsequent rapid cooling, is associated with substantial clumping factors $\sim 10^0-10^1$. In the supersonic cases, the kinetic energy decay is divided into a rate-limiting phase of shock dissipation and a comparatively swift phase of thermal dissipation, predominantly occurring within the overdense regions. Upon turbulence driving turnoff, the strong density contrasts decay within a relatively brief timescale of $\sim 30 - 300~{\rm Myr}$, depending on the average gas density. Dense clouds are crushed on similar timescales of $ \sim 30 - 100 ~{\rm Myr} $, depending on turbulence driving strength but independent from average gas density. Results of this work also contribute a novel dataset of dissipation timescales that incorporates an understanding of kinematics and thermodynamics in addition to the traditional cooling rate tables, which may serve as a valuable asset for forthcoming simulations that aim to explore gas dynamics on galactic and cosmological scales. △ Less

Submitted 27 June, 2024; originally announced June 2024.

Comments: 24 pages, 20 figures

Hypergraph Unreliability in Quasi-Polynomial Time

Authors: Ruoxu Cen, Jason Li, Debmalya Panigrahi

Abstract: The hypergraph unreliability problem asks for the probability that a hypergraph gets disconnected when every hyperedge fails independently with a given probability. For graphs, the unreliability problem has been studied over many decades, and multiple fully polynomial-time approximation schemes are known starting with the work of Karger (STOC 1995). In contrast, prior to this work, no non-trivial… ▽ More The hypergraph unreliability problem asks for the probability that a hypergraph gets disconnected when every hyperedge fails independently with a given probability. For graphs, the unreliability problem has been studied over many decades, and multiple fully polynomial-time approximation schemes are known starting with the work of Karger (STOC 1995). In contrast, prior to this work, no non-trivial result was known for hypergraphs (of arbitrary rank). In this paper, we give quasi-polynomial time approximation schemes for the hypergraph unreliability problem. For any fixed $\varepsilon \in (0, 1)$, we first give a $(1+\varepsilon)$-approximation algorithm that runs in $m^{O(\log n)}$ time on an $m$-hyperedge, $n$-vertex hypergraph. Then, we improve the running time to $m\cdot n^{O(\log^2 n)}$ with an additional exponentially small additive term in the approximation. △ Less

Submitted 27 March, 2024; originally announced March 2024.

Comments: To appear in STOC 2024

Cosmological imprints in the filament with DisPerSE

Authors: Ziyong Wu, Yu Luo, Wei Wang, Xi Kang, Renyue Cen

Abstract: In the regime of cosmology and large-scale structure formation, filaments are vital components of the cosmic web. This study employs statistical methods to examine the formation, evolution, and cosmological constraints of filaments identified by DisPerSe. We run large-sample of N-body simulations to study the filament length and mass functions and their evolution. In general, the filament length f… ▽ More In the regime of cosmology and large-scale structure formation, filaments are vital components of the cosmic web. This study employs statistical methods to examine the formation, evolution, and cosmological constraints of filaments identified by DisPerSe. We run large-sample of N-body simulations to study the filament length and mass functions and their evolution. In general, the filament length function can be fitted by a power law with both the normalization and power index dependent on redshift and cosmological parameters. It is discovered that filament properties are influenced by various cosmological parameters, with $σ_8$ and $n_s$ exhibiting slightly stronger dependence than $Ω_m$. We also uncover a three-stage filament formation process from $z \sim 3$ to $z \sim 1$: rapid formation of long and short filaments from $z \sim 3$ to $z \sim 2$, persistence of long filaments from $z \sim 2$ to $z \sim 1$, followed by fragmentation and increased prevalence of shorter filaments below $z \sim 1$. Additionally, both large and small-scale perturbations in power spectrum influence filament growth, with a characteristic scale around $k \sim 0.265$, wherein changes in matter fluctuations larger than this scale result in the formation of long filaments, while changes smaller than this scale have the opposite effect. These insights enhance our understanding of filament evolution and their cosmological relevance and also highlight their potential cosmological applications in observations. △ Less

Submitted 23 February, 2024; originally announced February 2024.

Comments: 23 pages, 37 figures, 1 table

Efficient survey design for finding high-redshift galaxies with JWST

Authors: Luka Vujeva, Charles L. Steinhardt, Christian Kragh Jespersen, Brenda L. Frye, Anton M. Koekemoer, Priyamvada Natarajan, Andreas L. Faisst, Pascale Hibon, Lukas J. Furtak, Hakim Atek, Renyue Cen, Albert Sneppen

Abstract: Several large JWST blank field observing programs have not yet discovered the first galaxies expected to form at $15 \leq z \leq 20$. This has motivated the search for more effective survey strategies that will be able to effectively probe this redshift range. Here, we explore the use of gravitationally lensed cluster fields, that have historically been the most effective discovery tool with HST.… ▽ More Several large JWST blank field observing programs have not yet discovered the first galaxies expected to form at $15 \leq z \leq 20$. This has motivated the search for more effective survey strategies that will be able to effectively probe this redshift range. Here, we explore the use of gravitationally lensed cluster fields, that have historically been the most effective discovery tool with HST. In this paper, we analyze the effectiveness of the most massive galaxy clusters that provide the highest median magnification factor within a single JWST NIRCam module in uncovering this population. The results of exploiting these lensing clusters to break the $z > 15$ barrier are compared against the results from large area, blank field surveys such as JADES and CEERS in order to determine the most effective survey strategy for JWST. We report that the fields containing massive foreground galaxy clusters specifically chosen to occupy the largest fraction of a single NIRCam module with high magnification factors in the source plane, whilst containing all multiple images in the image plane within a single module provide the highest probability of both probing the $15 \leq z \leq 20$ regime, as well as discovering the highest redshift galaxy possible with JWST. We also find that using multiple massive clusters in exchange for shallower survey depths is a more time efficient method of probing the $z > 15$ regime. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Comments: 12 pages, 5 figures

NIHAO-GIZMO: A Comparison of Simulated Disc Galaxies from GASOLINE and GIZMO

Authors: Hou-Zun Chen, Xi Kang, Andrea V. Macciò, Tobias Buck, Renyue Cen

Abstract: We utilize the public GIZMO code to simulate twelve disc galaxies from the NIHAO suite simulated with the GASOLINE code, then compare the corresponding galaxies in the two simulations. We find that while both codes with the same initial conditions and large-scale environments can successfully produce similar disc galaxies, significant differences are still seen in many properties of the galaxies,… ▽ More We utilize the public GIZMO code to simulate twelve disc galaxies from the NIHAO suite simulated with the GASOLINE code, then compare the corresponding galaxies in the two simulations. We find that while both codes with the same initial conditions and large-scale environments can successfully produce similar disc galaxies, significant differences are still seen in many properties of the galaxies, particularly in the circumgalactic medium (CGM) environment they reside. Specifically, the thermal feedback recipe used in GASOLINE results in ubiquitous long-lasting collimated outflows, primarily driven by high-density hot interstellar medium (ISM) from the galaxy center, and inflows of gas not aligned with the outflow cools rapidly and flows towards the galactic center. In contrast, galaxies from GIZMO code do not exhibit large-scale outflows at low redshifts, but instead display quasi-virialized hot gaseous halos that arise from the strong interaction between inflow of gas and feedback driven outflow. Therefore, the origins of mass and angular momentum of the cold disc in the two simulations are quite different, even though the final morphologies of corresponding galaxies are similar at $z\sim0$. The differences in the distribution of CGM gas are mainly due to different feedback models implemented in the two codes, thus future observations of CGM provide valuable insight into the physics governing the baryon cycle in disc galaxies. △ Less

Submitted 19 October, 2023; originally announced October 2023.

Comments: 17 pages, 11 figures, Submitted to ApJ, Comments welcome

Scientific Objectives of the Hot Universe Baryon Surveyor (HUBS) Mission

Authors: Joel Bregman, Renyue Cen, Yang Chen, Wei Cui, Taotao Fang, Fulai Guo, Edmund Hodges-Kluck, Rui Huang, Luis C. Ho, Li Ji, Suoqing Ji, Xi Kang, Xiaoyu Lai, Hui Li, Jiangtao Li, Miao Li, Xiangdong Li, Yuan Li, Zhaosheng Li, Guiyun Liang, Helei Liu, Wenhao Liu, Fangjun Lu, Junjie Mao, Gabriele Ponti , et al. (29 additional authors not shown)

Abstract: The Hot Universe Baryon Surveyor (HUBS) is a proposed space-based X-ray telescope for detecting X-ray emissions from the hot gas content in our universe. With its unprecedented spatially-resolved high-resolution spectroscopy and large field of view, the HUBS mission will be uniquely qualified to measure the physical and chemical properties of the hot gas in the interstellar medium, the circumgalac… ▽ More The Hot Universe Baryon Surveyor (HUBS) is a proposed space-based X-ray telescope for detecting X-ray emissions from the hot gas content in our universe. With its unprecedented spatially-resolved high-resolution spectroscopy and large field of view, the HUBS mission will be uniquely qualified to measure the physical and chemical properties of the hot gas in the interstellar medium, the circumgalactic medium, the intergalactic medium, and the intracluster medium. These measurements will be valuable for two key scientific goals of HUBS, namely to unravel the AGN and stellar feedback physics that governs the formation and evolution of galaxies, and to probe the baryon budget and multi-phase states from galactic to cosmological scales. In addition to these two goals, the HUBS mission will also help us solve some problems in the fields of galaxy clusters, AGNs, diffuse X-ray backgrounds, supernova remnants, and compact objects. This paper discusses the perspective of advancing these fields using the HUBS telescope. △ Less

Submitted 11 July, 2023; originally announced July 2023.

Comments: 52 pages, 22 figures. Accepted for publication in Science China: Physics, Mechanics and Astronomy

Self-consistent Combined HST, K-band, and Spitzer Photometric Catalogs of the BUFFALO Survey Fields

Authors: Amanda Pagul, F. Javier Sánchez, Iary Davidzon, Anton M. Koekemoer, Hakim Atek, Renyue Cen, Lukas J. Furtak, Mathilde Jauzac, Guillaume Mahler, Bahram Mobasher, Mireia Montes, Mario Nonino, Keren Sharon, Charles L. Steinhardt, John R. Weaver

Abstract: This manuscript presents new astronomical source catalogs using data from the BUFFALO Survey. These catalogs contain detailed information for over 100,000 astronomical sources in the 6 BUFFALO clusters: Abell 370, Abell 2744, Abell S1063, MACS 0416, MACS 0717, and MACS 1149 spanning a total 240 arcmin^2. The catalogs include positions and forced photometry measurements of these objects in the F275… ▽ More This manuscript presents new astronomical source catalogs using data from the BUFFALO Survey. These catalogs contain detailed information for over 100,000 astronomical sources in the 6 BUFFALO clusters: Abell 370, Abell 2744, Abell S1063, MACS 0416, MACS 0717, and MACS 1149 spanning a total 240 arcmin^2. The catalogs include positions and forced photometry measurements of these objects in the F275W, F336W, F435W, F606W, F814W, F105W, F125W, F140W, and F160W HST -bands, Keck-NIRC2/VLT-HAWKI Ks band, and IRAC Channel 1 and 2 bands. Additionally, we include photometry measurements in the F475W, F625W, and F110W bands for Abell 370. This catalog also includes photometric redshift estimates computed via template fitting using LePhare. When comparing to spectroscopic reference, we obtain an outlier fraction of 9.2% and scatter, normalized median absolute deviation (NMAD), of 0.062. The catalogs are publicly available for their use by the community. △ Less

Submitted 10 July, 2023; originally announced July 2023.

Comments: 16 pages, to be submitted to ApJS

Inspiraling streams of enriched gas observed around a massive galaxy 11 billion years ago

Authors: Shiwu Zhang, Zheng Cai, Dandan Xu, Rhythm Shimakawa, Fabrizio Arrigoni Battaia, Jason Xavier Prochaska, Renyue Cen, Zheng Zheng, Yunjing Wu, Qiong Li, Liming Dou, Jianfeng Wu, Ann Zabludoff, Xiaohui Fan, Yanli Ai, Emmet Gabriel Golden-Marx, Miao Li, Youjun Lu, Xiangcheng Ma, Sen Wang, Ran Wang, Feng Yuan

Abstract: Stars form in galaxies, from gas that has been accreted from the intergalactic medium. Simulations have shown that recycling of gas-the reaccretion of gas that was previously ejected from a galaxy-could sustain star formation in the early Universe. We observe the gas surrounding a massive galaxy at redshift 2.3 and detect emission lines from neutral hydrogen, helium, and ionized carbon that extend… ▽ More Stars form in galaxies, from gas that has been accreted from the intergalactic medium. Simulations have shown that recycling of gas-the reaccretion of gas that was previously ejected from a galaxy-could sustain star formation in the early Universe. We observe the gas surrounding a massive galaxy at redshift 2.3 and detect emission lines from neutral hydrogen, helium, and ionized carbon that extend 100 kiloparsecs from the galaxy. The kinematics of this circumgalactic gas is consistent with an inspiraling stream. The carbon abundance indicates that the gas had already been enriched with elements heavier than helium, previously ejected from a galaxy. We interpret the results as evidence of gas recycling during high-redshift galaxy assembly. △ Less

Submitted 3 May, 2023; originally announced May 2023.

Comments: Published in Science, 5 May 2023 (accepted version), Main text 20 pages, four figures in the main text, and 13 figures and 4 tables in the supplementary materials;

Beyond the Quadratic Time Barrier for Network Unreliability

Authors: Ruoxu Cen, William He, Jason Li, Debmalya Panigrahi

Abstract: Karger (STOC 1995) gave the first FPTAS for the network (un)reliability problem, setting in motion research over the next three decades that obtained increasingly faster running times, eventually leading to a $\tilde{O}(n^2)$-time algorithm (Karger, STOC 2020). This represented a natural culmination of this line of work because the algorithmic techniques used can enumerate $Θ(n^2)$ (near)-minimum… ▽ More Karger (STOC 1995) gave the first FPTAS for the network (un)reliability problem, setting in motion research over the next three decades that obtained increasingly faster running times, eventually leading to a $\tilde{O}(n^2)$-time algorithm (Karger, STOC 2020). This represented a natural culmination of this line of work because the algorithmic techniques used can enumerate $Θ(n^2)$ (near)-minimum cuts. In this paper, we go beyond this quadratic barrier and obtain a faster FPTAS for the network unreliability problem. Our algorithm runs in $m^{1+o(1)} + \tilde{O}(n^{1.5})$ time. Our main contribution is a new estimator for network unreliability in very reliable graphs. These graphs are usually the bottleneck for network unreliability since the disconnection event is elusive. Our estimator is obtained by defining an appropriate importance sampling subroutine on a dual spanning tree packing of the graph. To complement this estimator for very reliable graphs, we use recursive contraction for moderately reliable graphs. We show that an interleaving of sparsification and contraction can be used to obtain a better parametrization of the recursive contraction algorithm that yields a faster running time matching the one obtained for the very reliable case. △ Less

Submitted 20 July, 2023; v1 submitted 13 April, 2023; originally announced April 2023.

Generative Adversarial Networks for Scintillation Signal Simulation in EXO-200

Authors: S. Li, I. Ostrovskiy, Z. Li, L. Yang, S. Al Kharusi, G. Anton, I. Badhrees, P. S. Barbeau, D. Beck, V. Belov, T. Bhatta, M. Breidenbach, T. Brunner, G. F. Cao, W. R. Cen, C. Chambers, B. Cleveland, M. Coon, A. Craycraft, T. Daniels, L. Darroch, S. J. Daugherty, J. Davis, S. Delaquis, A. Der Mesrobian-Kabakian , et al. (65 additional authors not shown)

Abstract: Generative Adversarial Networks trained on samples of simulated or actual events have been proposed as a way of generating large simulated datasets at a reduced computational cost. In this work, a novel approach to perform the simulation of photodetector signals from the time projection chamber of the EXO-200 experiment is demonstrated. The method is based on a Wasserstein Generative Adversarial N… ▽ More Generative Adversarial Networks trained on samples of simulated or actual events have been proposed as a way of generating large simulated datasets at a reduced computational cost. In this work, a novel approach to perform the simulation of photodetector signals from the time projection chamber of the EXO-200 experiment is demonstrated. The method is based on a Wasserstein Generative Adversarial Network - a deep learning technique allowing for implicit non-parametric estimation of the population distribution for a given set of objects. Our network is trained on real calibration data using raw scintillation waveforms as input. We find that it is able to produce high-quality simulated waveforms an order of magnitude faster than the traditional simulation approach and, importantly, generalize from the training sample and discern salient high-level features of the data. In particular, the network correctly deduces position dependency of scintillation light response in the detector and correctly recognizes dead photodetector channels. The network output is then integrated into the EXO-200 analysis framework to show that the standard EXO-200 reconstruction routine processes the simulated waveforms to produce energy distributions comparable to that of real waveforms. Finally, the remaining discrepancies and potential ways to improve the approach further are highlighted. △ Less

Submitted 8 May, 2023; v1 submitted 11 March, 2023; originally announced March 2023.

Comments: As accepted by JINST

Journal ref: JINST 18 P06005 2023

Search for Two-neutrino Double-Beta Decay of $^{136}\rm Xe$ to the $0^+_1$ excited state of $^{136}\rm Ba$ with the Complete EXO-200 Dataset

Authors: EXO-200 Collaboration, :, S. Al Kharusi, G. Anton, I. Badhrees, P. S. Barbeau, D. Beck, V. Belov, T. Bhatta, M. Breidenbach, T. Brunner, G. F. Cao, W. R. Cen, C. Chambers, B. Cleveland, M. Coon, A. Craycraft, T. Daniels, L. Darroch, S. J. Daugherty, J. Davis, S. Delaquis, A. Der Mesrobian-Kabakian, R. DeVoe, J. Dilling , et al. (83 additional authors not shown)

Abstract: A new search for two-neutrino double-beta ($2νββ$) decay of $^{136}\rm Xe$ to the $0^+_1$ excited state of $^{136}\rm Ba$ is performed with the full EXO-200 dataset. A deep learning-based convolutional neural network is used to discriminate signal from background events. Signal detection efficiency is increased relative to previous searches by EXO-200 by more than a factor of two. With the additio… ▽ More A new search for two-neutrino double-beta ($2νββ$) decay of $^{136}\rm Xe$ to the $0^+_1$ excited state of $^{136}\rm Ba$ is performed with the full EXO-200 dataset. A deep learning-based convolutional neural network is used to discriminate signal from background events. Signal detection efficiency is increased relative to previous searches by EXO-200 by more than a factor of two. With the addition of the Phase II dataset taken with an upgraded detector, the median 90$\%$ confidence level half-life sensitivity of $2νββ$ decay to the $0^+_1$ state of $^{136}\rm Ba$ is $2.9 \times 10^{24}~\rm yr$ using a total $^{136}\rm Xe$ exposure of $234.1~\rm kg~yr$. No statistically significant evidence for $2νββ$ decay to the $0^+_1$ state is observed, leading to a lower limit of $T^{2ν}_{1/2}(0^+ \rightarrow 0^+_1) > 1.4\times10^{24}~\rm yr$ at 90$\%$ confidence level, improved by 70$\%$ relative to the current world's best constraint. △ Less

Submitted 16 October, 2023; v1 submitted 2 March, 2023; originally announced March 2023.

Comments: 9 pages, 7 figures, 2 tables

Journal ref: Chinese Physics C, Volume 47, Number 10, 2023

UV & Ly$α$ halos of Ly$α$ emitters across environments at z=2.84

Authors: Satoshi Kikuta, Yuichi Matsuda, Shigeki Inoue, Charles C. Steidel, Renyue Cen, Zheng Zheng, Hidenobu Yajima, Rieko Momose, Masatoshi Imanishi, Yutaka Komiyama

Abstract: We present UV and Ly$α$ radial surface brightness (SB) profiles of Ly$α$ emitters (LAEs) at $z=2.84$ detected with the Hyper Suprime-Cam (HSC) on the Subaru Telescope. The depth of our data, together with the wide field coverage including a protocluster, enable us to study the dependence of Ly$α$ halos (LAHs) on various galaxy properties, including Mpc-scale environments. UV and Ly$α$ images of 34… ▽ More We present UV and Ly$α$ radial surface brightness (SB) profiles of Ly$α$ emitters (LAEs) at $z=2.84$ detected with the Hyper Suprime-Cam (HSC) on the Subaru Telescope. The depth of our data, together with the wide field coverage including a protocluster, enable us to study the dependence of Ly$α$ halos (LAHs) on various galaxy properties, including Mpc-scale environments. UV and Ly$α$ images of 3490 LAEs are extracted, and stacking the images yields SB sensitivity of $\sim1\times10^{-20}\mathrm{~erg~s^{-1}~cm^{-2}~arcsec^{-2}}$ in Ly$α$, reaching the expected level of optically thick gas illuminated by the UV background at $z\sim3$. Fitting of the two-component exponential function gives the scale-lengths of $1.56\pm0.01$ and $10.4\pm0.3$ pkpc. Dividing the sample according to their photometric properties, we find that while the dependence of halo scale-length on environment outside of the protocluster core is not clear, LAEs in the central regions of protoclusters appear to have very large LAHs which could be caused by combined effects of source overlapping and diffuse Ly$α$ emission from cool intergalactic gas permeating the forming protocluster core irradiated by active members. For the first time, we identify ``UV halos'' around bright LAEs which are probably due to a few lower-mass satellite galaxies. Through comparison with recent numerical simulations, we conclude that, while scattered Ly$α$ photons from the host galaxies are dominant, star formation in satellites evidently contributes to LAHs, and that fluorescent Ly$α$ emission may be boosted within protocluster cores at cosmic noon and/or near bright QSOs. △ Less

Submitted 24 February, 2023; originally announced February 2023.

Comments: 29 pages, 20 Figures, Accepted for publication in ApJ

Steiner Connectivity Augmentation and Splitting-off in Poly-logarithmic Maximum Flows

Authors: Ruoxu Cen, William He, Jason Li, Debmalya Panigrahi

Abstract: We give an almost-linear time algorithm for the Steiner connectivity augmentation problem: given an undirected graph, find a smallest (or minimum weight) set of edges whose addition makes a given set of terminals $τ$-connected (for any given $τ> 0$). The running time of our algorithm is dominated by polylogarithmic calls to any maximum flow subroutine; using the recent almost-linear time maximum f… ▽ More We give an almost-linear time algorithm for the Steiner connectivity augmentation problem: given an undirected graph, find a smallest (or minimum weight) set of edges whose addition makes a given set of terminals $τ$-connected (for any given $τ> 0$). The running time of our algorithm is dominated by polylogarithmic calls to any maximum flow subroutine; using the recent almost-linear time maximum flow algorithm (Chen et al., FOCS 2022), we get an almost-linear running time for our algorithm as well. This is tight up to the polylogarithmic factor even for just two terminals. Prior to our work, an almost-linear (in fact, near-linear) running time was known only for the special case of global connectivity augmentation, i.e., when all vertices are terminals (Cen et al., STOC 2022). We also extend our algorithm to the closely related Steiner splitting-off problem, where the edges incident on a vertex have to be {\em split-off} while maintaining the (Steiner) connectivity of a given set of terminals. Prior to our work, a nearly-linear time algorithm was known only for the special case of global connectivity (Cen et al., STOC 2022). The only known generalization beyond global connectivity was to preserve all pairwise connectivities using a much slower algorithm that makes $n$ calls to an all-pairs maximum flow (or Gomory-Hu tree) subroutine (Lau and Yung, SICOMP 2013), as against polylog(n) calls to a (single-pair) maximum flow subroutine in this work. △ Less

Submitted 10 November, 2022; originally announced November 2022.

Comments: To appear in SODA 2023

Bridging the Gap between Cosmic Dawn and Reionization favors Faint Galaxies-dominated Models

Authors: Ankita Bera, Sultan Hassan, Aaron Smith, Renyue Cen, Enrico Garaldi, Rahul Kannan, Mark Vogelsberger

Abstract: It has been claimed that traditional models struggle to explain the tentative detection of the 21\,cm absorption trough centered at $z\sim17$ measured by the EDGES collaboration. On the other hand, it has been shown that the EDGES results are consistent with an extrapolation of a declining UV luminosity density, following a simple power-law of deep Hubble Space Telescope observations of… ▽ More It has been claimed that traditional models struggle to explain the tentative detection of the 21\,cm absorption trough centered at $z\sim17$ measured by the EDGES collaboration. On the other hand, it has been shown that the EDGES results are consistent with an extrapolation of a declining UV luminosity density, following a simple power-law of deep Hubble Space Telescope observations of $4 < z < 9$ galaxies. We here explore the conditions by which the EDGES detection is consistent with current reionization and post-reionization observations, including the neutral hydrogen fraction at $z\sim6$--$8$, Thomson scattering optical depth, and ionizing emissivity at $z\sim5$. By coupling a physically motivated source model derived from radiative transfer hydrodynamic simulations of reionization to a Markov Chain Monte Carlo sampler, we find that it is entirely possible to reconcile the high-redshift (cosmic dawn) and low-redshift (reionization) existing constraints. In particular, we find that high contribution from low-mass halos along with high photon escape fractions are required to simultaneously reproduce cosmic dawn and reionization constraints. Our analysis further confirms that low-mass galaxies produce a flatter emissivity evolution, which leads to an earlier onset of reionization with gradual and longer duration, resulting in a higher optical depth. While our faint-galaxies dominated models successfully reproduce the measured globally averaged quantities over the first one billion years, they underestimate the late redshift-instantaneous measurements in efficiently star-forming and massive systems. We show that our (simple) physically-motivated semi-analytical prescription produces consistent results with the (sophisticated) state-of-the-art \thesan radiation-magneto-hydrodynamic simulation of reionization. △ Less

Submitted 27 October, 2023; v1 submitted 28 September, 2022; originally announced September 2022.

Comments: 14 pages, 6 figures. Accepted for publication in ApJ. Comments are welcome

Search for MeV Electron Recoils from Dark Matter in EXO-200

Authors: EXO-200 Collaboration, :, S. Al Kharusi, G. Anton, I. Badhrees, P. S. Barbeau, D. Beck, V. Belov, T. Bhatta, M. Breidenbach, T. Brunner, G. F. Cao, W. R. Cen, C. Chambers, B. Cleveland, M. Coon, A. Craycraft, T. Daniels, L. Darroch, S. J. Daugherty, J. Davis, S. Delaquis, A. Der Mesrobian-Kabakian, R. DeVoe, J. Dilling , et al. (83 additional authors not shown)

Abstract: We present a search for electron-recoil signatures from the charged-current absorption of fermionic dark matter using the EXO-200 detector. We report an average electron recoil background rate of $6.8 \times 10^{-4}\, \mathrm{cts}\,\mathrm{kg}^{-1}\mathrm{yr}^{-1}\mathrm{keV}^{-1}$ above $4\,\mathrm{MeV}$ and find no statistically significant excess over our background projection. Using a total… ▽ More We present a search for electron-recoil signatures from the charged-current absorption of fermionic dark matter using the EXO-200 detector. We report an average electron recoil background rate of $6.8 \times 10^{-4}\, \mathrm{cts}\,\mathrm{kg}^{-1}\mathrm{yr}^{-1}\mathrm{keV}^{-1}$ above $4\,\mathrm{MeV}$ and find no statistically significant excess over our background projection. Using a total ${}^{136}\mathrm{Xe}$ exposure of $234.1\,\mathrm{kg}\,\mathrm{yr}$ we exclude new parameter space for the charged-current absorption cross-section for dark matter masses between $m_χ= 2.6\,\mathrm{MeV} - 11.6\,\mathrm{MeV}$ with a minimum of $6\times 10^{-51}\,\mathrm{cm}^2$ at $8.3\,\mathrm{MeV}$ at the $90\%$ confidence level. △ Less

Submitted 20 February, 2023; v1 submitted 2 July, 2022; originally announced July 2022.

Comments: 9 pages, 6 figures, 1 table

Journal ref: Phys. Rev. D 107, 012007 (2023)

Edge Connectivity Augmentation in Near-Linear Time

Authors: Ruoxu Cen, Jason Li, Debmalya Panigrahi

Abstract: We give an $\tilde{O}(m)$-time algorithm for the edge connectivity augmentation problem and the closely related edge splitting-off problem. This is optimal up to lower order terms and closes the long line of work on these problems. We give an $\tilde{O}(m)$-time algorithm for the edge connectivity augmentation problem and the closely related edge splitting-off problem. This is optimal up to lower order terms and closes the long line of work on these problems. △ Less

Submitted 9 May, 2022; originally announced May 2022.

Comments: To be published in STOC 2022

Patchy Kinetic Sunyaev-Zel'dovich Effect with Controlled Reionization History and Morphology

Authors: Nianyi Chen, Hy Trac, Suvodip Mukherjee, Renyue Cen

Abstract: Using the novel semi-numerical code for reionization AMBER, we model the patchy kinetic Sunyaev-Zel'dovich (kSZ) effect by directly specifying the reionization history with the redshift midpoint $z_\mathrm{mid}$, duration $Δ_\mathrm{z}$, and asymmetry $A_\mathrm{z}$. We further control the ionizing sources and radiation through the minimum halo mass $M_\mathrm{h}$ and the radiation mean free path… ▽ More Using the novel semi-numerical code for reionization AMBER, we model the patchy kinetic Sunyaev-Zel'dovich (kSZ) effect by directly specifying the reionization history with the redshift midpoint $z_\mathrm{mid}$, duration $Δ_\mathrm{z}$, and asymmetry $A_\mathrm{z}$. We further control the ionizing sources and radiation through the minimum halo mass $M_\mathrm{h}$ and the radiation mean free path $λ_\mathrm{mfp}$. AMBER reproduces the free electron number density and the patchy kSZ power spectrum of radiation-hydrodynamic simulations at the target resolution ($1\,{\rm Mpc}/h$) with matched reionization parameters. With a suite of $(2\,{\rm Gpc}/h)^3$ simulations using AMBER, we first constrain the redshift midpoint $6.0<z_{\rm mid}<8.9$ using the Planck2018 Thomson optical depth result (95\% CL). Then, assuming $z_{\rm mid}=8$, we find that the amplitude of $D^{\rm pkSZ}_{\ell=3000}$ scales linearly with the duration of reionization $Δ_z$, and is consistent with the $1σ$ upper limit from the South Pole Telescope (SPT) results up to $Δ_z<5.1$ ($Δ_z$ encloses $5\%$ to $95\%$ ionization). Moreover, a shorter $λ_{\rm mfp}$ can lead to a $\sim 10\%$ lower $D^{\rm pkSZ}_{\ell=3000}$ and a flatter slope in the $Δ_z-D^{\rm pkSZ}_{\ell=3000}$ scaling relation, thereby affecting the constraints on $Δ_z$ at $\ell=3000$. Allowing $z_{\rm mid}$ and $λ_{\rm mfp}$ to vary simultaneously, we get spectra consistent with the SPT result ($95\%$ CL) up to $Δ_z=12.8$ (but $A_z>8$ is needed to ensure an end of reionization before $z=5.5$). We show that constraints on the asymmetry require $\sim 0.1\,μk^2$ measurement accuracy at multipoles other than $\ell=3000$. Finally, we find that the amplitude and shape of the kSZ spectrum are only weakly sensitive to $M_h$ under a fixed reionization history and radiation mean-free path. △ Less

Submitted 8 March, 2022; originally announced March 2022.

Comments: 20 pages, 17 Figures; comments welcome

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

Authors: David J. Lagattuta, Johan Richard, Franz Erik Bauer, Catherine Cerny, Adélaïde Claeyssens, Lucia Guaita, Mathilde Jauzac, Alexandre Jeanneau, Anton M. Koekemoer, Guillaume Mahler, Gonzalo Prieto Lyon, Matteo Bianconi, Thomas Connor, Renyue Cen, Alastair Edge, Andreas L. Faisst, Marceau Limousin, Richard Massey, Mauro Sereno, Keren Sharon, John R. Weaver

Abstract: 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 n… ▽ More 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. △ Less

Submitted 9 February, 2022; originally announced February 2022.

Comments: 20 pages, 12 figures, 3 tables (including 1 in an appendix). Accepted in MNRAS. Data Release products available at https://astro.dur.ac.uk/~hbpn39/pilot-wings.html

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

Authors: Yuri Oku, Kengo Tomida, Kentaro Nagamine, Ikkoh Shimizu, Renyue Cen

Abstract: 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 scale… ▽ More 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. △ Less

Submitted 26 July, 2022; v1 submitted 3 January, 2022; originally announced January 2022.

Comments: 36 pages, 22 figures, accepted for publication in The Astrophysical Journal Supplement Series. Movies available at https://www.yurioku.com/research/

Minimum Cuts in Directed Graphs via Partial Sparsification

Authors: Ruoxu Cen, Jason Li, Danupon Nanongkai, Debmalya Panigrahi, Kent Quanrud, Thatchaphol Saranurak

Abstract: We give an algorithm to find a minimum cut in an edge-weighted directed graph with $n$ vertices and $m$ edges in $\tilde O(n\cdot \max(m^{2/3}, n))$ time. This improves on the 30 year old bound of $\tilde O(nm)$ obtained by Hao and Orlin for this problem. Our main technique is to reduce the directed mincut problem to $\tilde O(\min(n/m^{1/3}, \sqrt{n}))$ calls of {\em any} maxflow subroutine. Usin… ▽ More We give an algorithm to find a minimum cut in an edge-weighted directed graph with $n$ vertices and $m$ edges in $\tilde O(n\cdot \max(m^{2/3}, n))$ time. This improves on the 30 year old bound of $\tilde O(nm)$ obtained by Hao and Orlin for this problem. Our main technique is to reduce the directed mincut problem to $\tilde O(\min(n/m^{1/3}, \sqrt{n}))$ calls of {\em any} maxflow subroutine. Using state-of-the-art maxflow algorithms, this yields the above running time. Our techniques also yield fast {\em approximation} algorithms for finding minimum cuts in directed graphs. For both edge and vertex weighted graphs, we give $(1+ε)$-approximation algorithms that run in $\tilde O(n^2 / ε^2)$ time. △ Less

Submitted 17 November, 2021; originally announced November 2021.

Comments: To appear in FOCS 2021. This paper subsumes arXiv:2104.06933 and arXiv:2104.07898

Augmenting Edge Connectivity via Isolating Cuts

Authors: Ruoxu Cen, Jason Li, Debmalya Panigrahi

Abstract: We give an algorithm for augmenting the edge connectivity of an undirected graph by using the isolating cuts framework (Li and Panigrahi, FOCS '20). Our algorithm uses poly-logarithmic calls to any max-flow algorithm, which yields a running time of $\tilde O(m + n^{3/2})$ and improves on the previous best time of $\tilde O(n^2)$ (Benczúr and Karger, SODA '98) for this problem. We also obtain an id… ▽ More We give an algorithm for augmenting the edge connectivity of an undirected graph by using the isolating cuts framework (Li and Panigrahi, FOCS '20). Our algorithm uses poly-logarithmic calls to any max-flow algorithm, which yields a running time of $\tilde O(m + n^{3/2})$ and improves on the previous best time of $\tilde O(n^2)$ (Benczúr and Karger, SODA '98) for this problem. We also obtain an identical improvement in the running time of the closely related edge splitting off problem in undirected graphs. △ Less

Submitted 3 November, 2021; originally announced November 2021.

Comments: To be published in SODA 2022

AMBER: A Semi-numerical Abundance Matching Box for the Epoch of Reionization

Authors: Hy Trac, Nianyi Chen, Ian Holst, Marcelo A. Alvarez, Renyue Cen

Abstract: 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.… ▽ More 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. △ Less

Submitted 10 May, 2022; v1 submitted 21 September, 2021; originally announced September 2021.

Comments: AMBER is publicly available at https://github.com/hytrac/amber

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

Authors: Robin Kooistra, Shigeki Inoue, Khee-Gan Lee, Renyue Cen, Naoki Yoshida

Abstract: 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 cosm… ▽ More 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. △ Less

Submitted 25 January, 2022; v1 submitted 21 September, 2021; originally announced September 2021.

Comments: Accepted for publication in ApJ, 24 pages, 16 figures, 1 table

Search for Majoron-emitting modes of $^{136}$Xe double beta decay with the complete EXO-200 dataset

Authors: S. Al Kharusi, G. Anton, I. Badhrees, P. S. Barbeau, D. Beck, V. Belov, T. Bhatta, M. Breidenbach, T. Brunner, G. F. Cao, W. R. Cen, C. Chambers, B. Cleveland, M. Coon, A. Craycraft, T. Daniels, L. Darroch, S. J. Daugherty, J. Davis, S. Delaquis, A. Der Mesrobian-Kabakian, R. DeVoe, J. Dilling, A. Dolgolenko, M. J. Dolinski , et al. (81 additional authors not shown)

Abstract: A search for Majoron-emitting modes of the neutrinoless double-beta decay of $^{136}$Xe is performed with the full EXO-200 dataset. This dataset consists of a total $^{136}$Xe exposure of 234.1 kg$\cdot$yr, and includes data with detector upgrades that have improved the energy threshold relative to previous searches. A lower limit of T$_{1/2}^{\rm{^{136}Xe}}>$4.3$\cdot$10$^{24}$ yr at 90\% C.L. on… ▽ More A search for Majoron-emitting modes of the neutrinoless double-beta decay of $^{136}$Xe is performed with the full EXO-200 dataset. This dataset consists of a total $^{136}$Xe exposure of 234.1 kg$\cdot$yr, and includes data with detector upgrades that have improved the energy threshold relative to previous searches. A lower limit of T$_{1/2}^{\rm{^{136}Xe}}>$4.3$\cdot$10$^{24}$ yr at 90\% C.L. on the half-life of the spectral index $n=1$ Majoron decay was obtained, a factor of 3.6 more stringent than the previous limit from EXO-200, corresponding to a constraint on the Majoron-neutrino coupling constant of $|\langle g_{ee}^{M}\rangle|$$<(0.4$-$0.9)\cdot10^{-5}$. The lower threshold and the additional data taken resulted in a factor 8.4 improvement for the $n=7$ mode compared to the previous EXO search. This search provides the most stringent limits to-date on the Majoron-emitting decays of $^{136}$Xe with spectral indices $n=1,2,3,$ and 7. △ Less

Submitted 17 November, 2021; v1 submitted 3 September, 2021; originally announced September 2021.

The EXO-200 detector, part II: Auxiliary Systems

Authors: N. Ackerman, J. Albert, M. Auger, D. J. Auty, I. Badhrees, P. S. Barbeau, L. Bartoszek, E. Baussan, V. Belov, C. Benitez-Medina, T. Bhatta, M. Breidenbach, T. Brunner, G. F. Cao, W. R. Cen, C. Chambers, B. Cleveland, R. Conley, S. Cook, M. Coon, W. Craddock, A. Craycraft, W. Cree, T. Daniels, L. Darroch , et al. (135 additional authors not shown)

Abstract: The EXO-200 experiment searched for neutrinoless double-beta decay of $^{136}$Xe with a single-phase liquid xenon detector. It used an active mass of 110 kg of 80.6%-enriched liquid xenon in an ultra-low background time projection chamber with ionization and scintillation detection and readout. This paper describes the design and performance of the various support systems necessary for detector op… ▽ More The EXO-200 experiment searched for neutrinoless double-beta decay of $^{136}$Xe with a single-phase liquid xenon detector. It used an active mass of 110 kg of 80.6%-enriched liquid xenon in an ultra-low background time projection chamber with ionization and scintillation detection and readout. This paper describes the design and performance of the various support systems necessary for detector operation, including cryogenics, xenon handling, and controls. Novel features of the system were driven by the need to protect the thin-walled detector chamber containing the liquid xenon, to achieve high chemical purity of the Xe, and to maintain thermal uniformity across the detector. △ Less

Submitted 22 October, 2021; v1 submitted 13 July, 2021; originally announced July 2021.

Comments: Manuscript updated in response to JINST reviewer comments

Minimum Cuts in Directed Graphs via $\sqrt{n}$ Max-Flows

Authors: Ruoxu Cen, Jason Li, Danupon Nanongkai, Debmalya Panigrahi, Thatchaphol Saranurak

Abstract: We give an algorithm to find a mincut in an $n$-vertex, $m$-edge weighted directed graph using $\tilde O(\sqrt{n})$ calls to any maxflow subroutine. Using state of the art maxflow algorithms, this yields a directed mincut algorithm that runs in $\tilde O(m\sqrt{n} + n^2)$ time. This improves on the 30 year old bound of $\tilde O(mn)$ obtained by Hao and Orlin for this problem. We give an algorithm to find a mincut in an $n$-vertex, $m$-edge weighted directed graph using $\tilde O(\sqrt{n})$ calls to any maxflow subroutine. Using state of the art maxflow algorithms, this yields a directed mincut algorithm that runs in $\tilde O(m\sqrt{n} + n^2)$ time. This improves on the 30 year old bound of $\tilde O(mn)$ obtained by Hao and Orlin for this problem. △ Less

Submitted 16 April, 2021; originally announced April 2021.

Spectral Signatures of Population III and Envelope-stripped Stars in Galaxies at the Epoch of Reionization

Authors: Elizabeth Berzin, Amy Secunda, Renyue Cen, Alexander Menegas, Ylva Götberg

Abstract: While most simulations of the epoch of reionization have focused on single-stellar populations in star-forming dwarf galaxies, products of binary evolution are expected to significantly contribute to emissions of hydrogen-ionizing photons. Among these products are stripped stars (or helium stars), which have their envelopes stripped from interactions with binary companions, leaving an exposed heli… ▽ More While most simulations of the epoch of reionization have focused on single-stellar populations in star-forming dwarf galaxies, products of binary evolution are expected to significantly contribute to emissions of hydrogen-ionizing photons. Among these products are stripped stars (or helium stars), which have their envelopes stripped from interactions with binary companions, leaving an exposed helium core. Previous work has suggested these stripped stars can dominate the LyC photon output of high-redshift low luminosity galaxies. Other sources of hard radiation in the early universe include zero-metallicity Population III stars, which may have similar SED properties to galaxies with radiation dominated by stripped star emissions. Here, we use two metrics (the power-law exponent over wavelength intervals 240-500 Å, 600-900 Å, and 1200-2000 Å, and the ratio of total luminosity in FUV wavelengths to LyC wavelengths) to compare the SEDs of simulated galaxies with only single-stellar evolution, galaxies containing stripped stars, and galaxies containing Population III stars, with four different IMFs. We find that stripped stars significantly alter the SEDs in the LyC range of galaxies at the epoch of reionization. SEDs in galaxies with stripped stars present have lower power-law indices in the LyC range and lower FUV to LyC luminosity ratios. These differences in SEDs are present at all considered luminosities ($M_{UV} > -15$, AB system), and are most pronounced for lower luminosity galaxies. We also find that SEDs of galaxies with stripped stars and Pop III stars are distinct from each other for all tested IMFs. △ Less

Submitted 10 August, 2021; v1 submitted 16 February, 2021; originally announced February 2021.

Comments: 16 pages, 9 figures. Accepted to ApJ

Physics of Non-Universal Larson's Relation

Authors: Renyue Cen

Abstract: 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's Relation relating velocity dispersion ($σ_R$) to cloud size ($R$). We report on two key new findings. First, the correct form of the Larson's Relation is… ▽ More 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's Relation relating velocity dispersion ($σ_R$) to cloud size ($R$). We report on two key new findings. First, the correct form of the Larson'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'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. △ Less

Submitted 3 December, 2020; originally announced December 2020.

Comments: 12 pages, 2 figures, accepted to ApJ Letters

Delayed Photons from Binary Evolution Help Reionize the Universe

Authors: Amy Secunda, Renyue Cen, Taysun Kimm, Ylva Gotberg, Selma E. de Mink

Abstract: High-resolution numerical simulations including feedback and aimed at calculating the escape fraction (fesc) of hydrogen-ionizing photons often assume stellar radiation based on single-stellar population synthesis models. However, strong evidence suggests the binary fraction of massive stars is 70%. Moreover, simulations so far yield values of fesc falling only on the lower end of the roughly 10-2… ▽ More High-resolution numerical simulations including feedback and aimed at calculating the escape fraction (fesc) of hydrogen-ionizing photons often assume stellar radiation based on single-stellar population synthesis models. However, strong evidence suggests the binary fraction of massive stars is 70%. Moreover, simulations so far yield values of fesc falling only on the lower end of the roughly 10-20% range, the amount presumed necessary to reionize the Universe. Analyzing a high-resolution (4 pc) cosmological radiation hydrodynamic simulation we study how fesc changes when we include two different products of binary stellar evolution - stars stripped of their hydrogen envelopes and massive blue stragglers. Both produce significant amounts of ionizing photons 10-200 Myr after each starburst. We find the relative importance of these photons are amplified with respect to escaped ionizing photons, because peaks in star formation rates (SFRs) and fesc are often out of phase by this 10-200 Myr. Additionally, low mass, bursty galaxies emit Lyman continuum radiation primarily from binary products when SFRs are low. Observations of these galaxies by the James Webb Space Telescope could provide crucial information on the evolution of binary stars as a function of redshift. Overall, including stripped stars and massive blue stragglers increases our photon-weighted mean escape fraction by around 13% and 10%, respectively, resulting in a mean fesc of 17%. Our results emphasize that using updated stellar population synthesis models with binary stellar evolution provides a more sound physical basis for stellar reionization. △ Less

Submitted 29 July, 2020; originally announced July 2020.

Comments: 20 pages, 7 figures, submitted to ApJ

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

Authors: Zack Li, Renyue Cen

Abstract: 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 be… ▽ More 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 "Assembly Conformity" 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 "baryonic environmental" 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. △ Less

Submitted 21 July, 2020; originally announced July 2020.

Comments: 15 pages, 15 figures

Journal ref: The Astrophysical Journal (2020), Volume 898, Page 39

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

Authors: Takuma Izumi, John D. Silverman, Knud Jahnke, Andreas Schulze, Renyue Cen, Malte Schramm, Tohru Nagao, Lutz Wisotzki, Wiphu Rujopakarn

Abstract: 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… ▽ More 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. △ Less

Submitted 6 July, 2020; v1 submitted 4 June, 2020; originally announced June 2020.

Comments: 17 pages, 8 figures, 4 tables (including 1 table and 1 figure in the appendix). Accepted for publication in ApJ; Minor updates after the acceptance

Sparsification of Directed Graphs via Cut Balance

Authors: Ruoxu Cen, Yu Cheng, Debmalya Panigrahi, Kevin Sun

Abstract: In this paper, we consider the problem of designing cut sparsifiers and sketches for directed graphs. To bypass known lower bounds, we allow the sparsifier/sketch to depend on the balance of the input graph, which smoothly interpolates between undirected and directed graphs. We give nearly matching upper and lower bounds for both for-all (cf. Benczúr and Karger, STOC 1996) and for-each (Andoni et… ▽ More In this paper, we consider the problem of designing cut sparsifiers and sketches for directed graphs. To bypass known lower bounds, we allow the sparsifier/sketch to depend on the balance of the input graph, which smoothly interpolates between undirected and directed graphs. We give nearly matching upper and lower bounds for both for-all (cf. Benczúr and Karger, STOC 1996) and for-each (Andoni et al., ITCS 2016) cut sparsifiers/sketches as a function of cut balance, defined the maximum ratio of the cut value in the two directions of a directed graph (Ene et al., STOC 2016). We also show an interesting application of digraph sparsification via cut balance by using it to give a very short proof of a celebrated maximum flow result of Karger and Levine (STOC 2002). △ Less

Submitted 11 May, 2021; v1 submitted 2 June, 2020; originally announced June 2020.

SCORCH. III. Analytical Models of Reionization with Varying Clumping Factors

Authors: Nianyi Chen, Aristide Doussot, Hy Trac, Renyue Cen

Abstract: 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… ▽ More 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. △ Less

Submitted 26 January, 2021; v1 submitted 16 April, 2020; originally announced April 2020.

Comments: 14 pages, 8 figures

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

Authors: David V. Bowen, Todd M. Tripp, Edward B. Jenkins, Max Pettini, Doron Chelouche, Renyue Cen, Donald G. York

Abstract: 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… ▽ More 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'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'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}\:$. △ Less

Submitted 30 March, 2020; v1 submitted 13 March, 2020; originally announced March 2020.

Comments: Accepted for publication in ApJ. This version: typos corrected, references updated

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

Authors: Yiping Ao, Zheng Zheng, Christian Henkel, Shiyu Nie, Alexandre Beelen, Renyue Cen, Mark Dijkstra, Paul J. Francis, James E. Geach, Kotaro Kohno, Matthew D. Lehnert, Karl M. Menten, Junzhi Wang, Axel Weiss

Abstract: 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/ou… ▽ More 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. △ Less

Submitted 13 March, 2020; originally announced March 2020.

Comments: 30 pages, 7 figures, submitted to Nature Astronomy. The full-text access to a view-only version of the published paper is available by using the following SharedIt link: https://rdcu.be/b2JkM

Measurement of the Spectral Shape of the beta-decay of 137Xe to the Ground State of 137Cs in EXO-200 and Comparison with Theory

Authors: S. Al Kharusi, G. Anton, I. Badhrees, P. S. Barbeau, D. Beck, V. Belov, T. Bhatta, M. Breidenbach, T. Brunner, G. F. Cao, W. R. Cen, C. Chambers, B. Cleveland, M. Coon, A. Craycraft, T. Daniels, L. Darroch, S. J. Daugherty, J. Davis, S. Delaquis, A. Der Mesrobian-Kabakian, R. DeVoe, J. Dilling, A. Dolgolenko, M. J. Dolinski , et al. (83 additional authors not shown)

Abstract: We report on a comparison between the theoretically predicted and experimentally measured spectra of the first-forbidden non-unique $β$-decay transition $^{137}\textrm{Xe}(7/2^-)\to\,^{137}\textrm{Cs}(7/2^+)$. The experimental data were acquired by the EXO-200 experiment during a deployment of an AmBe neutron source. The ultra-low background environment of EXO-200, together with dedicated source d… ▽ More We report on a comparison between the theoretically predicted and experimentally measured spectra of the first-forbidden non-unique $β$-decay transition $^{137}\textrm{Xe}(7/2^-)\to\,^{137}\textrm{Cs}(7/2^+)$. The experimental data were acquired by the EXO-200 experiment during a deployment of an AmBe neutron source. The ultra-low background environment of EXO-200, together with dedicated source deployment and analysis procedures, allowed for collection of a pure sample of the decays, with an estimated signal-to-background ratio of more than 99-to-1 in the energy range from 1075 to 4175 keV. In addition to providing a rare and accurate measurement of the first-forbidden non-unique $β$-decay shape, this work constitutes a novel test of the calculated electron spectral shapes in the context of the reactor antineutrino anomaly and spectral bump. △ Less

Submitted 7 May, 2020; v1 submitted 31 January, 2020; originally announced February 2020.

Comments: Version as accepted by PRL

Journal ref: Phys. Rev. Lett. 124, 232502 (2020)

Physics of Prodigious Lyman Continuum Leakers

Authors: Renyue Cen

Abstract: 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 ra… ▽ More 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. △ Less

Submitted 29 January, 2020; originally announced January 2020.

Comments: 15 pages, 4 figures, in press of ApJL

The BUFFALO HST Survey

Authors: Charles L. Steinhardt, Mathilde Jauzac, Ana Acebron, Hakim Atek, Peter Capak, Iary Davidzon, Dominique Eckert, David Harvey, Anton M. Koekemoer, Claudia D. P. Lagos, Guillaume Mahler, Mireia Montes, Anna Niemiec, Mario Nonino, P. A. Oesch, Johan Richard, Steven A. Rodney, Matthieu Schaller, Keren Sharon, Louis-Gregory Strolger, Joseph Allingham, Adam Amara, Yannick Bah'e, Celine Boehm, Sownak Bose , et al. (70 additional authors not shown)

Abstract: 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 no… ▽ More 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. △ Less

Submitted 13 February, 2020; v1 submitted 27 January, 2020; originally announced January 2020.

Comments: Accepted ApJS; MAST archive will be live concurrent with publication

Indoor Localization System of ROS mobile robot based on Visible Light Communication

Authors: Weipeng Guan, Shihuan Chen, Shangsheng Wen, Wenyuan Hou, Zequn Tan, Ruihong Cen

Abstract: In this paper, an indoor robot localization system based on Robot Operating System (ROS) and visible light communication (VLC) is presented. On the basis of our previous work, we innovatively designed a VLC localization and navigation package based on Robot Operating System (ROS), which contains the LED-ID detection and recognition method, the video target tracking algorithm and the double-lamp po… ▽ More In this paper, an indoor robot localization system based on Robot Operating System (ROS) and visible light communication (VLC) is presented. On the basis of our previous work, we innovatively designed a VLC localization and navigation package based on Robot Operating System (ROS), which contains the LED-ID detection and recognition method, the video target tracking algorithm and the double-lamp positioning algorithm. This package exploited the principle of double-lamp positioning and the loose coupling characteristics of the ROS system, which is implemented by loosely coupled ROS nodes. Consequently, this paper combines ROS and VLC, aiming at promoting the application of VLC positioning in mature robotic systems. Moreover, it pushed forward the development of localization application based on VLC technology and lays a foundation for transplanting to other ROS robot platforms. Experimental results show that the proposed system can provide indoor localization within 1 cm and possesses a good real-time performance which takes only 0.4 seconds for one-time positioning. And if a high-performance laptop is used, the single positioning time can be reduced to 0.08 seconds. Therefore, this study confirms the practical application and the superior performance of VLC technology in ROS robot, showing the great potential of VLC localization. T he video demo of the proposed robot positioning system based on VLC can be seen in * △ Less

Submitted 6 January, 2020; originally announced January 2020.

On Post-Starburst Galaxies Dominating Tidal Disruption Events

Authors: Renyue Cen

Abstract: A starburst induced by a galaxy merger may create a relatively thin central stellar disk at radius $\le 100$pc. We calculate the rate of tidal disruption events (TDEs) by the inspiraling secondary supermassive black (SMBH) through the disk. With a small enough stellar velocity dispersion ($σ/v_c \le 0.1$) in the disk, it is shown that $10^5-10^6$ TDEs of solar-type main sequence stars per post-sta… ▽ More A starburst induced by a galaxy merger may create a relatively thin central stellar disk at radius $\le 100$pc. We calculate the rate of tidal disruption events (TDEs) by the inspiraling secondary supermassive black (SMBH) through the disk. With a small enough stellar velocity dispersion ($σ/v_c \le 0.1$) in the disk, it is shown that $10^5-10^6$ TDEs of solar-type main sequence stars per post-starburst galaxy (PSB) can be produced to explain their dominance in producing observed TDEs. Although the time it takes to bring the secondary SMBH to the disk apparently varies in the range of $\sim 0.1-1$Gyr since the starburst, depending on its landing location and subsequently due to dynamical friction with stars exterior to the central stellar disk in question, the vast majority of TDEs by the secondary SMBH in any individual PSB occurs within a space of time shorter than $\sim 30$Myr. Five unique testable predictions of this model are suggested. △ Less

Submitted 9 December, 2019; originally announced December 2019.

Comments: 21 pages, 4 figures, to appear in ApJ Letters

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

Authors: Yutaka Fujita, Renyue Cen, Irina Zhuravleva

Abstract: 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… ▽ More 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. △ Less

Submitted 20 May, 2020; v1 submitted 2 December, 2019; originally announced December 2019.

Journal ref: MNRAS, 494, 5507-5519 (2020)

Roundtrip Spanners with $(2k-1)$ Stretch

Authors: Ruoxu Cen, Ran Duan, Yong Gu

Abstract: A roundtrip spanner of a directed graph $G$ is a subgraph of $G$ preserving roundtrip distances approximately for all pairs of vertices. Despite extensive research, there is still a small stretch gap between roundtrip spanners in directed graphs and undirected graphs. For a directed graph with real edge weights in $[1,W]$, we first propose a new deterministic algorithm that constructs a roundtrip… ▽ More A roundtrip spanner of a directed graph $G$ is a subgraph of $G$ preserving roundtrip distances approximately for all pairs of vertices. Despite extensive research, there is still a small stretch gap between roundtrip spanners in directed graphs and undirected graphs. For a directed graph with real edge weights in $[1,W]$, we first propose a new deterministic algorithm that constructs a roundtrip spanner with $(2k-1)$ stretch and $O(k n^{1+1/k}\log (nW))$ edges for every integer $k> 1$, then remove the dependence of size on $W$ to give a roundtrip spanner with $(2k-1)$ stretch and $O(k n^{1+1/k}\log n)$ edges. While keeping the edge size small, our result improves the previous $2k+ε$ stretch roundtrip spanners in directed graphs [Roditty, Thorup, Zwick'02; Zhu, Lam'18], and almost matches the undirected $(2k-1)$-spanner with $O(n^{1+1/k})$ edges [Althöfer et al. '93] when $k$ is a constant, which is optimal under Erdös conjecture. △ Less

Submitted 30 April, 2020; v1 submitted 27 November, 2019; originally announced November 2019.

Comments: 12 pages

Measurements of electron transport in liquid and gas Xenon using a laser-driven photocathode

Authors: O. Njoya, T. Tsang, M. Tarka, W. Fairbank, K. S. Kumar, T. Rao, T. Wager, S. Al Kharusi, G. Anton, I. J. Arnquist, I. Badhrees, P. S. Barbeau, D. Beck, V. Belov, T. Bhatta, J. P. Brodsky, E. Brown, T. Brunner, E. Caden, G. F. Cao, L. Cao, W. R. Cen, C. Chambers, B. Chana, S. A. Charlebois , et al. (131 additional authors not shown)

Abstract: Measurements of electron drift properties in liquid and gaseous xenon are reported. The electrons are generated by the photoelectric effect in a semi-transparent gold photocathode driven in transmission mode with a pulsed ultraviolet laser. The charges drift and diffuse in a small chamber at various electric fields and a fixed drift distance of 2.0 cm. At an electric field of 0.5 kV/cm, the measur… ▽ More Measurements of electron drift properties in liquid and gaseous xenon are reported. The electrons are generated by the photoelectric effect in a semi-transparent gold photocathode driven in transmission mode with a pulsed ultraviolet laser. The charges drift and diffuse in a small chamber at various electric fields and a fixed drift distance of 2.0 cm. At an electric field of 0.5 kV/cm, the measured drift velocities and corresponding temperature coefficients respectively are $1.97 \pm 0.04$ mm/$μ$s and $(-0.69\pm0.05)$\%/K for liquid xenon, and $1.42 \pm 0.03$ mm/$μ$s and $(+0.11\pm0.01)$\%/K for gaseous xenon at 1.5 bar. In addition, we measure longitudinal diffusion coefficients of $25.7 \pm 4.6$ cm$^2$/s and $149 \pm 23$ cm$^2$/s, for liquid and gas, respectively. The quantum efficiency of the gold photocathode is studied at the photon energy of 4.73 eV in liquid and gaseous xenon, and vacuum. These charge transport properties and the behavior of photocathodes in a xenon environment are important in designing and calibrating future large scale noble liquid detectors. △ Less

Submitted 24 November, 2019; originally announced November 2019.

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

Authors: John D. Silverman, Tommaso Treu, Xuheng Ding, Knud Jahnke, Vardha N. Bennert, Simon Birrer, Malte Schramm, Andreas Schulze, Jeyhan S. Kartaltepe, David B. Sanders, Renyue Cen

Abstract: 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… ▽ More 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. △ Less

Submitted 30 October, 2019; originally announced October 2019.

Comments: 4 pages; 1 figure; To be resubmitted to ApJ Letters imminently; minor corrections based on comments from the referee have been included; Further comments are welcome

The Physical Origins of the Identified and Still Missing Components of the Warm-Hot Intergalactic Medium: Insights from Deep Surveys in the Field of Blazar 1ES1553+113

Authors: Sean D. Johnson, John S. Mulchaey, Hsiao-Wen Chen, Nastasha A. Wijers, Thomas Connor, Sowgat Muzahid, Joop Schaye, Renyue Cen, Scott G. Carlsten, Jane Charlton, Maria R. Drout, Andy D. Goulding, Terese T. Hansen, Gregory L. Walth

Abstract: The relationship between galaxies and the state/chemical enrichment of the warm-hot intergalactic medium (WHIM) expected to dominate the baryon budget at low-z provides sensitive constraints on structure formation and galaxy evolution models. We present a deep redshift survey in the field of 1ES1553+113, a blazar with a unique combination of UV+X-ray spectra for surveys of the circum-/intergalacti… ▽ More The relationship between galaxies and the state/chemical enrichment of the warm-hot intergalactic medium (WHIM) expected to dominate the baryon budget at low-z provides sensitive constraints on structure formation and galaxy evolution models. We present a deep redshift survey in the field of 1ES1553+113, a blazar with a unique combination of UV+X-ray spectra for surveys of the circum-/intergalactic medium (CGM/IGM). Nicastro et al. 2018 reported the detection of two O VII WHIM absorbers at $z=0.4339$ and $0.3551$ in its spectrum, suggesting that the WHIM is metal-rich and sufficient to close the missing baryons problem. Our survey indicates that the blazar is a member of a $z=0.433$ group and that the higher-$z$ O VII candidate arises from its intragroup medium. The resulting bias precludes its use in baryon censuses. The $z=0.3551$ candidate occurs in an isolated environment 630 kpc from the nearest galaxy (with stellar mass $\log M_*/M_\odot \approx 9.7$) which we show is unexpected for the WHIM. Finally, we characterize the galactic environments of broad H I Ly$α$ absorbers (Doppler widths of $b=40-80$ \kms; $T\lesssim4\times10^5$ K) which provide metallicity independent WHIM probes. On average, broad Ly$α$, absorbers are ${\approx}2\times$ closer to the nearest luminous ($L>0.25 L_*$) galaxy (700 kpc) than narrow ($b<30$ \kms; $T\lesssim4\times10^5$ K) ones (1300 kpc) but ${\approx}2\times$ further than O\,VI absorbers (350 kpc). These observations suggest that gravitational collapse heats portions of the IGM to form the WHIM but with feedback that does not enrich the IGM far beyond galaxy/group halos to levels currently observable in UV/X-ray metal lines. △ Less

Submitted 29 September, 2019; v1 submitted 25 September, 2019; originally announced September 2019.

Comments: Accepted to ApJ Letters; 10 pages, 4 figures, 1 table with the full version available on the journal webpage (or by request); v2 corrects one reference

Measurement of the scintillation and ionization response of liquid xenon at MeV energies in the EXO-200 experiment

Authors: EXO-200 Collaboration, :, G. Anton, I. Badhrees, P. S. Barbeau, D. Beck, V. Belov, T. Bhatta, M. Breidenbach, T. Brunner, G. F. Cao, W. R. Cen, C. Chambers, B. Cleveland, M. Coon, A. Craycraft, T. Daniels, L. Darroch, S. J. Daugherty, J. Davis, S. Delaquis, A. Der Mesrobian-Kabakian, R. DeVoe, J. Dilling, A. Dolgolenko , et al. (78 additional authors not shown)

Abstract: Liquid xenon (LXe) is employed in a number of current and future detectors for rare event searches. We use the EXO-200 experimental data to measure the absolute scintillation and ionization yields generated by $γ$ interactions from $^{228}$Th (2615~keV), $^{226}$Ra (1764~keV) and $^{60}$Co (1332~keV and 1173~keV) calibration sources, over a range of electric fields. The $W$-value that defines the… ▽ More Liquid xenon (LXe) is employed in a number of current and future detectors for rare event searches. We use the EXO-200 experimental data to measure the absolute scintillation and ionization yields generated by $γ$ interactions from $^{228}$Th (2615~keV), $^{226}$Ra (1764~keV) and $^{60}$Co (1332~keV and 1173~keV) calibration sources, over a range of electric fields. The $W$-value that defines the recombination-independent energy scale is measured to be $11.5~\pm~0.5$~(syst.)~$\pm~0.1$~(stat.) eV. These data are also used to measure the recombination fluctuations in the number of electrons and photons produced by the calibration sources at the MeV-scale, which deviate from extrapolations of lower-energy data. Additionally, a semi-empirical model for the energy resolution of the detector is developed, which is used to constrain the recombination efficiency, i.e., the fraction of recombined electrons that result in the emission of a detectable photon. Detailed measurements of the absolute charge and light yields for MeV-scale electron recoils are important for predicting the performance of future neutrinoless double beta decay detectors. △ Less

Submitted 15 June, 2020; v1 submitted 12 August, 2019; originally announced August 2019.

Journal ref: Phys. Rev. C 101, 065501 (2020)

Simulation of charge readout with segmented tiles in nEXO

Authors: Z. Li, W. R. Cen, A. Robinson, D. C. Moore, L. J. Wen, A. Odian, S. Al Kharusi, G. Anton, I. J. Arnquist, I. Badhrees, P. S. Barbeau, D. Beck, V. Belov, T. Bhatta, J. P. Brodsky, E. Brown, T. Brunner, E. Caden, G. F. Cao, L. Cao, C. Chambers, B. Chana, S. A. Charlebois, M. Chiu, B. Cleveland , et al. (128 additional authors not shown)

Abstract: nEXO is a proposed experiment to search for the neutrino-less double beta decay ($0νββ$) of $^{136}$Xe in a tonne-scale liquid xenon time projection chamber (TPC). The nEXO TPC will be equipped with charge collection tiles to form the anode. In this work, the charge reconstruction performance of this anode design is studied with a dedicated simulation package. A multi-variate method and a deep neu… ▽ More nEXO is a proposed experiment to search for the neutrino-less double beta decay ($0νββ$) of $^{136}$Xe in a tonne-scale liquid xenon time projection chamber (TPC). The nEXO TPC will be equipped with charge collection tiles to form the anode. In this work, the charge reconstruction performance of this anode design is studied with a dedicated simulation package. A multi-variate method and a deep neural network are developed to distinguish simulated $0νββ$ signals from backgrounds arising from trace levels of natural radioactivity in the detector materials. These simulations indicate that the nEXO TPC with charge-collection tiles shows promising capability to discriminate the $0νββ$ signal from backgrounds. The estimated half-life sensitivity for $0νββ$ decay is improved by $\sim$20$~(32)\%$ with the multi-variate~(deep neural network) methods considered here, relative to the sensitivity estimated in the nEXO pre-conceptual design report. △ Less

Submitted 11 October, 2019; v1 submitted 17 July, 2019; originally announced July 2019.

Journal ref: JINST 14 (2019) no.09, P09020

Search for Neutrinoless Double-Beta Decay with the Complete EXO-200 Dataset

Authors: G. Anton, I. Badhrees, P. S. Barbeau, D. Beck, V. Belov, T. Bhatta, M. Breidenbach, T. Brunner, G. F. Cao, W. R. Cen, C. Chambers, B. Cleveland, M. Coon, A. Craycraft, T. Daniels, M. Danilov, L. Darroch, S. J. Daugherty, J. Davis, S. Delaquis, A. Der Mesrobian-Kabakian, R. DeVoe, J. Dilling, A. Dolgolenko, M. J. Dolinski , et al. (77 additional authors not shown)

Abstract: A search for neutrinoless double-beta decay ($0νββ$) in $^{136}$Xe is performed with the full EXO-200 dataset using a deep neural network to discriminate between $0νββ$ and background events. Relative to previous analyses, the signal detection efficiency has been raised from 80.8% to 96.4$\pm$3.0% and the energy resolution of the detector at the Q-value of $^{136}$Xe $0νββ$ has been improved from… ▽ More A search for neutrinoless double-beta decay ($0νββ$) in $^{136}$Xe is performed with the full EXO-200 dataset using a deep neural network to discriminate between $0νββ$ and background events. Relative to previous analyses, the signal detection efficiency has been raised from 80.8% to 96.4$\pm$3.0% and the energy resolution of the detector at the Q-value of $^{136}$Xe $0νββ$ has been improved from $σ/E=1.23\%$ to $1.15\pm0.02\%$ with the upgraded detector. Accounting for the new data, the median 90% confidence level $0νββ$ half-life sensitivity for this analysis is $5.0 \cdot 10^{25}$ yr with a total $^{136}$Xe exposure of 234.1 kg$\cdot$yr. No statistically significant evidence for $0νββ$ is observed, leading to a lower limit on the $0νββ$ half-life of $3.5\cdot10^{25}$ yr at the 90% confidence level. △ Less

Submitted 18 October, 2019; v1 submitted 6 June, 2019; originally announced June 2019.

Comments: v1, 7 pages, 5 figures; v2, fix references; v3, update to accepted version

Journal ref: Phys. Rev. Lett. 123, 161802 (2019)

Ly$α$ view around a z=2.84 hyperluminous QSO at a node of the cosmic web

Authors: Satoshi Kikuta, Yuichi Matsuda, Renyue Cen, Charles C. Steidel, Masafumi Yagi, Tomoki Hayashino, Masatoshi Imanishi, Yutaka Komiyama, Rieko Momose, Tomoki Saito

Abstract: We report on the results of deep and wide-field (1.1 deg$^2$) narrow-band observations with Subaru/Hyper Suprime-Cam (HSC) of a field around a hyperluminous QSO (HLQSO), HS1549+1919, residing in a protocluster at $z=2.84$, to map the large-scale structure of Ly$α$ emitters (LAEs). One HSC pointing enables us to detect 3490 LAEs and 76 extended Ly$α$ blobs (LABs), probing diverse environments from… ▽ More We report on the results of deep and wide-field (1.1 deg$^2$) narrow-band observations with Subaru/Hyper Suprime-Cam (HSC) of a field around a hyperluminous QSO (HLQSO), HS1549+1919, residing in a protocluster at $z=2.84$, to map the large-scale structure of Ly$α$ emitters (LAEs). One HSC pointing enables us to detect 3490 LAEs and 76 extended Ly$α$ blobs (LABs), probing diverse environments from voids to protoclusters. The HLQSO is found to be near the center of the protocluster, which corresponds to the intersection of $\sim$100 cMpc-scale structures of LAEs. LABs are basically distributed along the large-scale structure, with larger ones particularly clustered around the HLQSO, confirming a previously noted tendency of LABs to prefer denser environments. Moreover, the shapes of LABs near the HLQSO appear to be aligned with the large-scale structure. Finally, a deep Ly$α$ image reveals a diffuse Ly$α$ nebula along a filamentary structure with no luminous UV/sub-mm counterpart. We suggest that the diffuse nebula is due to a cold filament with high clumping factor illuminated by the QSO, with a required high clumpiness provided by unresolved residing halos of mass $\leq 10^{9-10}M_\odot$. △ Less

Submitted 16 April, 2019; originally announced April 2019.

Comments: 6 pages, 3 figures, submitted to PASJ

On the Assembly Bias of Cool Core Clusters Traced by H$α$ Nebulae

Authors: Elinor Medezinski, Michael McDonald, Surhud More, Hironao Miyatake, Nicholas Battaglia, Massimo Gaspari, David Spergel, Renyue Cen

Abstract: Do cool-core (CC) and noncool-core (NCC) clusters live in different environments? We make novel use of H$α$ emission lines in the central galaxies of redMaPPer clusters as proxies to construct large (1,000's) samples of CC and NCC clusters, and measure their relative assembly bias using both clustering and weak lensing. We increase the statistical significance of the bias measurements from cluster… ▽ More Do cool-core (CC) and noncool-core (NCC) clusters live in different environments? We make novel use of H$α$ emission lines in the central galaxies of redMaPPer clusters as proxies to construct large (1,000's) samples of CC and NCC clusters, and measure their relative assembly bias using both clustering and weak lensing. We increase the statistical significance of the bias measurements from clustering by cross-correlating the clusters with an external galaxy redshift catalog from the Sloan Digital Sky Survey III, the LOWZ sample. Our cross-correlations can constrain assembly bias up to a statistical uncertainty of 6%. Given our H$α$ criteria for CC and NCC, we find no significant differences in their clustering amplitude. Interpreting this difference as the absence of halo assembly bias, our results rule out the possibility of having different large-scale (tens of Mpc) environments as the source of diversity observed in cluster cores. Combined with recent observations of the overall mild evolution of CC and NCC properties, such as central density and CC fraction, this would suggest that either the cooling properties of the cluster core are determined early on solely by the local (<200 kpc) gas properties at formation or that local merging leads to stochastic CC relaxation and disruption in a periodic way, preserving the average population properties over time. Studying the small-scale clustering in clusters at high redshift would help shed light on the exact scenario. △ Less

Submitted 12 March, 2019; originally announced March 2019.

Comments: 17 pages, 9 figures, 2 tables, to be submitted to ApJ; comments welcome