Showing 1–50 of 81 results for author: Ni, Y

Discrepancies Between JWST Observations and Simulations of Quenched Massive Galaxies at $z > 3$: A Comparative Study With IllustrisTNG and ASTRID

Authors: Emma Jane Weller, Fabio Pacucci, Yueying Ni, Lars Hernquist, Minjung Park

Abstract: Recent JWST observations have uncovered an unexpectedly large population of massive quiescent galaxies at $z>3$. Using the cosmological simulations IllustrisTNG and ASTRID, we identify analogous galaxies and investigate their abundance, formation, quenching mechanisms, and post-quenching evolution for stellar masses $9.5 < \log_{10}{(M_\star/{\rm M}_\odot)} < 12$. We apply three different quenchin… ▽ More Recent JWST observations have uncovered an unexpectedly large population of massive quiescent galaxies at $z>3$. Using the cosmological simulations IllustrisTNG and ASTRID, we identify analogous galaxies and investigate their abundance, formation, quenching mechanisms, and post-quenching evolution for stellar masses $9.5 < \log_{10}{(M_\star/{\rm M}_\odot)} < 12$. We apply three different quenching definitions and find that both simulations significantly underestimate the comoving number density of quenched massive galaxies at $z \gtrsim 3$ compared to JWST observations by up to $\sim 2$ dex. This fact highlights the necessity for improved physical models of AGN feedback in galaxy formation simulations. In both simulations, the high-$z$ quenched massive galaxies often host overmassive central black holes above the standard $M_{BH}-M_\star$ relation, implying that the AGN feedback plays a crucial role in quenching galaxies in the early Universe. The typical quenching timescales for these galaxies are $\sim 200-600$ Myr. IllustrisTNG primarily employs AGN kinetic feedback, while ASTRID relies on AGN thermal feedback, which is less effective and has a longer quenching timescale. We also study the post-quenching evolution of the high-$z$ massive quiescent galaxies and find that many experience subsequent reactivation of star formation, evolving into primary progenitors of $z=0$ brightest cluster galaxies. △ Less

Submitted 4 June, 2024; originally announced June 2024.

Comments: Submitted to The Astrophysical Journal. 13 pages, 12 figures

Debiasing with Diffusion: Probabilistic reconstruction of Dark Matter fields from galaxies with CAMELS

Authors: Victoria Ono, Core Francisco Park, Nayantara Mudur, Yueying Ni, Carolina Cuesta-Lazaro, Francisco Villaescusa-Navarro

Abstract: Galaxies are biased tracers of the underlying cosmic web, which is dominated by dark matter components that cannot be directly observed. Galaxy formation simulations can be used to study the relationship between dark matter density fields and galaxy distributions. However, this relationship can be sensitive to assumptions in cosmology and astrophysical processes embedded in the galaxy formation mo… ▽ More Galaxies are biased tracers of the underlying cosmic web, which is dominated by dark matter components that cannot be directly observed. Galaxy formation simulations can be used to study the relationship between dark matter density fields and galaxy distributions. However, this relationship can be sensitive to assumptions in cosmology and astrophysical processes embedded in the galaxy formation models, that remain uncertain in many aspects. In this work, we develop a diffusion generative model to reconstruct dark matter fields from galaxies. The diffusion model is trained on the CAMELS simulation suite that contains thousands of state-of-the-art galaxy formation simulations with varying cosmological parameters and sub-grid astrophysics. We demonstrate that the diffusion model can predict the unbiased posterior distribution of the underlying dark matter fields from the given stellar mass fields, while being able to marginalize over uncertainties in cosmological and astrophysical models. Interestingly, the model generalizes to simulation volumes approximately 500 times larger than those it was trained on, and across different galaxy formation models. Code for reproducing these results can be found at https://github.com/victoriaono/variational-diffusion-cdm △ Less

Submitted 15 March, 2024; originally announced March 2024.

Helium-deficient ER UMa-type dwarf nova below the period minimum with a hot secondary

Authors: Youngdae Lee, Dae-Sik Moon, Sang Chul Kim, Hong Soo Park, Yuan Qi Ni

Abstract: We present the discovery of a peculiar dwarf nova KSP-OT-201712a using high-cadence, multi-color observations made with the Korea Microlensing Telescope Network. KSP-OT-201712a exhibits a rare presence of outbursts during standstills as well as strong H$α$ emission for a dwarf nova below the period minimum with an orbital period of 58.75 $\pm$ 0.02 minutes. The outburst cycles are ~ 6.6 days withi… ▽ More We present the discovery of a peculiar dwarf nova KSP-OT-201712a using high-cadence, multi-color observations made with the Korea Microlensing Telescope Network. KSP-OT-201712a exhibits a rare presence of outbursts during standstills as well as strong H$α$ emission for a dwarf nova below the period minimum with an orbital period of 58.75 $\pm$ 0.02 minutes. The outburst cycles are ~ 6.6 days within standstills but increase to ~ 15 days outside of them. Both B-V and V-I colors become bluer and redder as the outburst luminosities increase and decrease, respectively, for the outburst within standstill, while they evolve in the opposite directions outside of the standstills. The presence of strong double-peaked H$α$ and weak He I emission lines with He/H flux ratio of 0.27, together with absorption lines of Mg b and Na D in the source, leads to the estimation Teff ~ 4570 $\pm$ 40 K, [Fe/H] ~ 0.06 $\pm$ 0.15 dex, and log g ~ 4.5 $\pm$ 0.1 for its secondary. KSP-OT-201712a is the second He-deficient dwarf nova below the period minimum, while the temperature of the secondary is measured for the first time in such objects. We identify it to be an ER UMa type dwarf nova suggesting that the evolution of dwarf novae across the period minimum is accompanied by large mass transfers. The high temperature of the secondary indicates that the system started its mass transfer when the secondary was about 93$\%$ of its main sequence age. The system will evolve to a helium cataclysmic variable or to AM CVn once its hydrogen envelope is exhausted before it explodes as a Type Ia supernova. △ Less

Submitted 5 February, 2024; originally announced February 2024.

Comments: 12 pages, 5 figures, accepted for publication in ApJ

The evolution of galaxy morphology from redshift z=6 to 3: Mock JWST observations of galaxies in the ASTRID simulation

Authors: Patrick LaChance, Rupert Croft, Yueying Ni, Nianyi Chen, Tiziana Di Matteo, Simeon Bird

Abstract: We present mock JWST observations for more than 215,000 different galaxies from the Astrid simulation with $3 \leq z \leq 6$. The mock observations are made using the BPASS stellar SED model, and a simple dust model. They are then viewed through NIRCam filters, convolved with a PSF, have noise added, and are drizzled together to emulate the Cosmic Evolution Early Release Science (CEERS) survey. We… ▽ More We present mock JWST observations for more than 215,000 different galaxies from the Astrid simulation with $3 \leq z \leq 6$. The mock observations are made using the BPASS stellar SED model, and a simple dust model. They are then viewed through NIRCam filters, convolved with a PSF, have noise added, and are drizzled together to emulate the Cosmic Evolution Early Release Science (CEERS) survey. We analyse this dataset by computing a number of morphological measures and find our catalog to have comparable statistics to similar mock catalogs, and the first release of CEERS data. We find that most of the Sersic indices of galaxies in our redshift range are lower than observed, with most having n less than one. Additionally, we observe the sizes of galaxies of all masses to increase from redshift z=6 to redshift z=3 consistent with other results. The number of galaxies in our catalog allows us to examine how relationships like the mass-size relation evolve with redshift, and compare the accuracy of a variety of traditional galaxy classification techniques (Sersic fit, Asymmetry-Concentration, and Gini-$M_{20}$) within our redshift range. We find the mass-size relation to be nearly flat at redshift z=6, and consistently increases as redshift decreases, and find the galaxy classification methods have minimal correlation with each other in our redshift range. We also investigate the impact that different stages of our imaging pipeline have on these morphological measures to determine how robust mock catalogs are to different choices at each step. Finally, we test the addition of incorporating light from AGNs into our pipeline and find that while the population of galaxies that have significant AGN luminosity is low, those galaxies do tend to have higher Sersic indices once the AGN luminosity is added, rectifying some of the systematic bias towards lower Sersic indices present in our dataset. △ Less

Submitted 29 January, 2024; originally announced January 2024.

Comments: 17 pages, 14 figures

Simulated host galaxy analogs of high-z quasars observed with JWST

Authors: Sabrina Berger, Madeline A. Marshall, J. Stuart B. Wyithe, Tiziana di Matteo, Yueying Ni, Stephen M. Wilkins

Abstract: The hosts of two low-luminosity high-z quasars, J2255+0251 and J2236+0032, were recently detected using JWST's NIRCam instrument. These represent the first high-z quasar host galaxy stellar detections and open a new window into studying high-z quasars. We examine the implications of the measured properties of J2255+0251 and J2236+0032 within the context of the hydrodynamic simulation BlueTides at… ▽ More The hosts of two low-luminosity high-z quasars, J2255+0251 and J2236+0032, were recently detected using JWST's NIRCam instrument. These represent the first high-z quasar host galaxy stellar detections and open a new window into studying high-z quasars. We examine the implications of the measured properties of J2255+0251 and J2236+0032 within the context of the hydrodynamic simulation BlueTides at z = 6.5. We find that these observed quasars fall on the BlueTides stellar to black hole mass relation and have similar luminosities to the brightest simulated quasars. We predict their star formation rates, estimating approximately $10^{2-3}$ $M_{\odot}/ \rm yr$ for both quasar hosts. J2255+0251 and J2236+0032's host galaxy radii also fall within estimates of the radii of the simulated host galaxies of similar luminosity quasars. We generate mock JWST NIRCam images of analogs to the observed quasars within BlueTides and perform a point source removal to illustrate both a qualitative and quantitative comparison of the measured and simulated radii and magnitudes. The quasar subtraction works well for similar luminosity quasars, and the recovered host images are consistent with what was observed for J2255+0251 and J2236+0032, further supporting the success of those observations. We also use our mock imaging pipeline to make predictions for the detection of J2255+0251 and J2236+0032's hosts in upcoming JWST observations. We anticipate that the simulation analogs of future high-z quasar host discoveries will allow us to make accurate predictions of their properties beyond the capabilities of JWST. △ Less

Submitted 18 April, 2024; v1 submitted 16 January, 2024; originally announced January 2024.

Comments: Accepted to MNRAS. 15 pages, 11 figures

z~2 dual AGN host galaxies are disky: stellar kinematics in the ASTRID Simulation

Authors: Ekaterine Dadiani, Tiziana Di Matteo, Nianyi Chen, Patrick Lachance, Yue Shen, Yu-Ching Chen, Rupert Croft, Yueying Ni, Simeon Bird

Abstract: We study dual AGN host galaxy morphologies at $z=2$ using the ASTRID simulation, selecting black hole (BH) pairs with small separation ($Δr<30\rm{kpc}$), high mass ($M_{\text{BH,12}}>10^7M_\odot$), and luminosity ($L_{\text{bol,12}}>10^{43}\rm{erg/s}$). We kinematically decompose (using MORDOR) $\sim1000$ dual AGN hosts into standard components - a `disk' (thin and thick disk, pseudo-bulge) and 'b… ▽ More We study dual AGN host galaxy morphologies at $z=2$ using the ASTRID simulation, selecting black hole (BH) pairs with small separation ($Δr<30\rm{kpc}$), high mass ($M_{\text{BH,12}}>10^7M_\odot$), and luminosity ($L_{\text{bol,12}}>10^{43}\rm{erg/s}$). We kinematically decompose (using MORDOR) $\sim1000$ dual AGN hosts into standard components - a `disk' (thin and thick disk, pseudo-bulge) and 'bulge' (bulge and halo) and define disk-dominated galaxies by the disk-to-total $D/T\geq0.5$. In ASTRID, $60.9\pm2.1\%$ of dual AGN hosts (independent of separation) are disk-dominated, with the $D/T$ distribution peaking at $\sim0.7$. Notably, hosts of BH pairs have similar morphologies (most either both disk or bulge-dominated). In dual-AGN hosts, the $D/T$ increases from $\sim17\% $ at $M_{\rm *}\sim 10^{9} M_{\odot}$ to $ 64\% $ for $M_{\rm *} \sim 10^{11.5} M_{\odot}$, and the pseudo-bulge is the dominant component of the disk fraction at the high mass end. Moreover, dual AGN hosts exhibit a higher fraction of disk/large pseudo-bulge than single-AGN hosts. The Disk-to-Total ratio is approximately constant with BH mass or AGN luminosity. We also create mock images of dual AGN host galaxies, employing morphological fitting software Statmorph to calculate morphological parameters and compare them with our kinematic decomposition results. Around $83.3\pm2.4\%$ of galaxies display disk-like profiles, of which $\sim60.7\pm2.2\%$ are kinematically confirmed as disks. Seŕsic indices and half-mass radii of dual AGN host galaxies align with observational measurements from HST at $z\sim2$. Around $34\%$ are identified as mergers from the $\text{Gini}-M_{20}$ relation. We find two dual AGN hosted by galaxies that exhibit disk-like seŕsic index $n_{12}<1$ and $(D/T)_{12}>0.5$, which are in remarkable agreement with properties of recently discovered dual quasars in disk galaxies at $z\sim 2$. △ Less

Submitted 21 December, 2023; originally announced December 2023.

Comments: 15 pages, 12 figures, submitted to the Open Journal of Astrophysics

MAGICS I. The First Few Orbits Encode the Fate of Seed Massive Black Hole Pairs

Authors: Nianyi Chen, Diptajyoti Mukherjee, Tiziana Di Matteo, Yueying Ni, Simeon Bird, Rupert Croft

Abstract: The elusive massive black hole (MBH) seeds stand to be revealed by the Laser Space Antenna Interferometer through mergers. As an aftermath of galaxy mergers, MBH coalescence is a vastly multi-scale process connected to galaxy formation. We introduce the "Massive black hole Assembly in Galaxies Informed by Cosmological Simulations" (MAGICS) suite, with galaxy/MBH properties and orbits recovered fro… ▽ More The elusive massive black hole (MBH) seeds stand to be revealed by the Laser Space Antenna Interferometer through mergers. As an aftermath of galaxy mergers, MBH coalescence is a vastly multi-scale process connected to galaxy formation. We introduce the "Massive black hole Assembly in Galaxies Informed by Cosmological Simulations" (MAGICS) suite, with galaxy/MBH properties and orbits recovered from large-volume cosmological simulation ASTRID. The simulations include subgrid star formation, supernovae feedback, and MBH accretion/feedback. In this first suite, we extract fifteen representative galaxy mergers with seed MBHs to examine their dynamics at an improved mass and spatial resolution (by $\sim2000$ and $\sim20$) and follow MBH orbits down to $\sim10\,\text{pc}$. We find that the seed MBH energy loss and orbital decay are largely governed by global torques induced by the galaxy merger process on scales resolvable by cosmological simulations. Specifically, pairs sink quickly if their orbits shrink rapidly below $1\,\text{kpc}$ during the first $\sim200\,\text{Myr}$ of pairing due to effective energy loss in major galaxy mergers, whereas MBHs gaining energy in minor galaxy mergers with head-on collisions are likely to stall. High initial eccentricities ($e_\text{init}>0.5$) and high stellar densities at kpc scales ($ρ_\text{star}>0.05\,M_\odot/\text{pc}^3$) also lead to most efficient decays. $\sim50\%$ high-redshift seed MBH pairs experience consecutive galaxy mergers and are more likely to stall at $\sim1\,\text{kpc}$. For a subset of systems, we carry out N-Body re-simulations until binary formation and find that some stalled systems merge at high-z when embedded in sufficient nuclear star clusters. △ Less

Submitted 25 April, 2024; v1 submitted 14 December, 2023; originally announced December 2023.

Comments: 21 pages, 20 Figures. Submitted to The Open Journal of Astrophysics. Comments welcome!

SN2023ixf in Messier 101: the twilight years of the progenitor as seen by Pan-STARRS

Authors: Conor L. Ransome, V. Ashley Villar, Anna Tartaglia, Sebastian Javier Gonzalez, Wynn V. Jacobson-Galán, Charles D. Kilpatrick, Raffaella Margutti, Ryan J. Foley, Matthew Grayling, Yuan Qi Ni, Ricardo Yarza, Christine Ye, Katie Auchettl, Thomas de Boer, Kenneth C. Chambers, David A. Coulter, Maria R. Drout, Diego Farias, Christa Gall, Hua Gao, Mark E. Huber, Adaeze L. Ibik, David O. Jones, Nandita Khetan, Chien-Cheng Lin , et al. (6 additional authors not shown)

Abstract: The nearby type II supernova, SN2023ixf in M101 exhibits signatures of early-time interaction with circumstellar material in the first week post-explosion. This material may be the consequence of prior mass loss suffered by the progenitor which possibly manifested in the form of a detectable pre-supernova outburst. We present an analysis of the long-baseline pre-explosion photometric data in $g$,… ▽ More The nearby type II supernova, SN2023ixf in M101 exhibits signatures of early-time interaction with circumstellar material in the first week post-explosion. This material may be the consequence of prior mass loss suffered by the progenitor which possibly manifested in the form of a detectable pre-supernova outburst. We present an analysis of the long-baseline pre-explosion photometric data in $g$, $w$, $r$, $i$, $z$ and $y$ filters from Pan-STARRS as part of the Young Supernova Experiment, spanning $\sim$5,000 days. We find no significant detections in the Pan-STARRS pre-explosion light curve. We train a multilayer perceptron neural network to classify pre-supernova outbursts. We find no evidence of eruptive pre-supernova activity to a limiting absolute magnitude of $-7$. The limiting magnitudes from the full set of $gwrizy$ (average absolute magnitude $\approx$-8) data are consistent with previous pre-explosion studies. We use deep photometry from the literature to constrain the progenitor of SN2023ixf, finding that these data are consistent with a dusty red supergiant (RSG) progenitor with luminosity $\log\left(L/L_\odot\right)$$\approx$5.12 and temperature $\approx$3950K, corresponding to a mass of 14-20 M$_\odot$ △ Less

Submitted 7 December, 2023; originally announced December 2023.

Comments: 19 pages, 8 figures, 1 table

Formation of a long filament through the connection of two filament segments observed by CHASE

Authors: H. T. Li, X. Cheng, Y. W. Ni, C. Li, S. H. Rao, J. H. Guo, M. D. Ding, P. F. Chen

Abstract: We present imaging and spectroscopic diagnostics of a long filament during its formation with the observations from the Chinese H$α$ Solar Explorer and Solar Dynamics Observatory. The seed filament first appeared at about 05:00 UT on 2022 September 13. Afterwards, it grew gradually and connected to another filament segment nearby, building up a long filament at about 20:00 UT on the same day. The… ▽ More We present imaging and spectroscopic diagnostics of a long filament during its formation with the observations from the Chinese H$α$ Solar Explorer and Solar Dynamics Observatory. The seed filament first appeared at about 05:00 UT on 2022 September 13. Afterwards, it grew gradually and connected to another filament segment nearby, building up a long filament at about 20:00 UT on the same day. The CHASE H$α$ spectra show an obvious centroid absorption with mild broadening at the main spine of the long filament, which is interpreted as the evidence of filament material accumulation. More interestingly, near the footpoints of the filament, persistent redshifts have been detected in the H$α$ spectra during the filament formation, indicating continuous drainage of filament materials. Furthermore, through inspecting the extreme ultraviolet images and magnetograms, it is found that EUV jets and brightenings appeared repeatedly at the junction of the two filament segments, where opposite magnetic polarities converged and canceled to each other continuously. These results suggest the occurrence of intermittent magnetic reconnection that not only connects magnetic structures of the two filament segments but also supplies cold materials for the filament channel likely by the condensation of injected hot plasma, even though a part of cold materials fall down to the filament footpoints at the same time. △ Less

Submitted 24 November, 2023; originally announced November 2023.

Comments: 11 pages, 6 figures, Accepted for publication in ApJL

Modelling the propagation of coronal mass ejections with COCONUT: implementation of the Regularized Biot-Savart Laws flux rope model

Authors: Jinhan Guo, L. Linan, S. Poedts, Y. Guo, A. Lani, B. Schmieder, M. Brchnelova, B. Perri, T. Baratashvili, Y. W. Ni, P. F. Chen

Abstract: Context: Coronal mass ejections (CMEs) are rapid eruptions of magnetized plasma that occur on the Sun, which are known as the main drivers of adverse space weather. Accurately tracking their evolution in the heliosphere in numerical models is of utmost importance for space weather forecasting. Aims: The main objective of this paper is to implement the Regularized Biot-Savart Laws (RBSL) method in… ▽ More Context: Coronal mass ejections (CMEs) are rapid eruptions of magnetized plasma that occur on the Sun, which are known as the main drivers of adverse space weather. Accurately tracking their evolution in the heliosphere in numerical models is of utmost importance for space weather forecasting. Aims: The main objective of this paper is to implement the Regularized Biot-Savart Laws (RBSL) method in a new global corona model COCONUT. This approach has the capability to construct the magnetic flux rope with an axis of arbitrary shape. Methods: We present the implementation process of the RBSL flux rope model in COCONUT, which is superposed onto a realistic solar wind reconstructed from the observed magnetogram around the minimum of solar activity. Based on this, we simulate the propagation of an S-shaped flux rope from the solar surface to a distance of 25 solar radii. Results: Our simulation successfully reproduces the birth process of a CME originating from a sigmoid in a self-consistent way. The model effectively captures various physical processes and retrieves the prominent features of the CMEs in observations. In addition, the simulation results indicate that the magnetic topology of the CME flux rope at around 20 solar radii deviates from a coherent structure, and manifests as a mix of open and closed field lines with diverse footpoints. Conclusions: This work demonstrates the potential of the RBSL flux rope model in reproducing CME events that are more consistent with observations. Moreover, our findings strongly suggest that magnetic reconnection during the CME propagation plays a critical role in destroying the coherent characteristic of a CME flux rope. △ Less

Submitted 22 November, 2023; originally announced November 2023.

Comments: 14 pages, 8 figures, accepted for publication in A&A

Probabilistic reconstruction of Dark Matter fields from biased tracers using diffusion models

Authors: Core Francisco Park, Victoria Ono, Nayantara Mudur, Yueying Ni, Carolina Cuesta-Lazaro

Abstract: Galaxies are biased tracers of the underlying cosmic web, which is dominated by dark matter components that cannot be directly observed. The relationship between dark matter density fields and galaxy distributions can be sensitive to assumptions in cosmology and astrophysical processes embedded in the galaxy formation models, that remain uncertain in many aspects. Based on state-of-the-art galaxy… ▽ More Galaxies are biased tracers of the underlying cosmic web, which is dominated by dark matter components that cannot be directly observed. The relationship between dark matter density fields and galaxy distributions can be sensitive to assumptions in cosmology and astrophysical processes embedded in the galaxy formation models, that remain uncertain in many aspects. Based on state-of-the-art galaxy formation simulation suites with varied cosmological parameters and sub-grid astrophysics, we develop a diffusion generative model to predict the unbiased posterior distribution of the underlying dark matter fields from the given stellar mass fields, while being able to marginalize over the uncertainties in cosmology and galaxy formation. △ Less

Submitted 14 November, 2023; originally announced November 2023.

Detecting Wandering Intermediate-Mass Black Holes with AXIS in the Milky Way and Local Massive Galaxies

Authors: Fabio Pacucci, Bryan Seepaul, Yueying Ni, Nico Cappelluti, Adi Foord

Abstract: This white paper explores the detectability of intermediate-mass black holes (IMBHs) wandering in the Milky Way (MW) and massive local galaxies, with a particular emphasis on the role of AXIS. IMBHs, ranging within $10^{3-6} \, M_\odot$, are commonly found at the centers of dwarf galaxies and may exist, yet undiscovered, in the MW. By using model spectra for advection-dominated accretion flows (AD… ▽ More This white paper explores the detectability of intermediate-mass black holes (IMBHs) wandering in the Milky Way (MW) and massive local galaxies, with a particular emphasis on the role of AXIS. IMBHs, ranging within $10^{3-6} \, M_\odot$, are commonly found at the centers of dwarf galaxies and may exist, yet undiscovered, in the MW. By using model spectra for advection-dominated accretion flows (ADAFs), we calculated the expected fluxes emitted by a population of wandering IMBHs with a mass of $10^5 \, M_\odot$ in various MW environments and extrapolated our results to massive local galaxies. Around $40\%$ of the potential population of wandering IMBHs in the MW can be detected in an AXIS deep field. We proposed criteria to aid in selecting IMBH candidates using already available optical surveys. We also showed that IMBHs wandering in $>200$ galaxies within $10$ Mpc can be easily detected with AXIS when passing within dense galactic environments (e.g., molecular clouds and cold neutral medium). In summary, we highlighted the potential X-ray detectability of wandering IMBHs in local galaxies and provided insights for guiding future surveys. Detecting wandering IMBHs is crucial for understanding their demographics, evolution, and the merging history of galaxies. △ Less

Submitted 14 November, 2023; originally announced November 2023.

Comments: This White Paper is part of a series commissioned for the AXIS Probe Concept Mission; additional AXIS White Papers can be found at the AXIS website: http://axis.astro.umd.edu/ with a mission overview here: arXiv:2311.00780. Review article, 7 pages, 3 figures

Data-driven Modeling of a Coronal Magnetic Flux Rope: from Birth to Death

Authors: J. H. Guo, Y. W. Ni, Y. Guo, C. Xia, B. Schmieder, S. Poedts, Z. Zhong, Y. H. Zhou, F. Yu, P. F. Chen

Abstract: Magnetic flux ropes are a bundle of twisted magnetic field lines produced by internal electric currents, which are responsible for solar eruptions and are the major drivers of geomagnetic storms. As such, it is crucial to develop a numerical model that can capture the entire evolution of a flux rope, from its birth to death, in order to predict whether adverse space weather events might occur or n… ▽ More Magnetic flux ropes are a bundle of twisted magnetic field lines produced by internal electric currents, which are responsible for solar eruptions and are the major drivers of geomagnetic storms. As such, it is crucial to develop a numerical model that can capture the entire evolution of a flux rope, from its birth to death, in order to predict whether adverse space weather events might occur or not. In this paper, we develop a data-driven modeling that combines a time-dependent magneto-frictional approach with a thermodynamic magnetohydrodynamic model. Our numerical modeling successfully reproduces the formation and confined eruption of an observed flux rope, and unveils the physical details behind the observations. Regarding the long-term evolution of the active region, our simulation results indicate that the flux cancellation due to collisional shearing plays a critical role in the formation of the flux rope, corresponding to a substantial increase in magnetic free energy and helicity. Regarding the eruption stage, the deformation of the flux rope during its eruption can cause an increase in the downward tension force, which suppresses it from further rising. This finding may shed light on why some torus-unstable flux ropes lead to failed eruptions after large-angle rotations. Moreover, we find that twisted fluxes can accumulate during the confined eruptions, which would breed the subsequent eruptive flares. △ Less

Submitted 30 October, 2023; originally announced October 2023.

Comments: 30 pages, 10 figures, Accepted for ApJ

Field-level simulation-based inference with galaxy catalogs: the impact of systematic effects

Authors: Natalí S. M. de Santi, Francisco Villaescusa-Navarro, L. Raul Abramo, Helen Shao, Lucia A. Perez, Tiago Castro, Yueying Ni, Christopher C. Lovell, Elena Hernandez-Martinez, Federico Marinacci, David N. Spergel, Klaus Dolag, Lars Hernquist, Mark Vogelsberger

Abstract: It has been recently shown that a powerful way to constrain cosmological parameters from galaxy redshift surveys is to train graph neural networks to perform field-level likelihood-free inference without imposing cuts on scale. In particular, de Santi et al. (2023) developed models that could accurately infer the value of $Ω_{\rm m}$ from catalogs that only contain the positions and radial velocit… ▽ More It has been recently shown that a powerful way to constrain cosmological parameters from galaxy redshift surveys is to train graph neural networks to perform field-level likelihood-free inference without imposing cuts on scale. In particular, de Santi et al. (2023) developed models that could accurately infer the value of $Ω_{\rm m}$ from catalogs that only contain the positions and radial velocities of galaxies that are robust to uncertainties in astrophysics and subgrid models. However, observations are affected by many effects, including 1) masking, 2) uncertainties in peculiar velocities and radial distances, and 3) different galaxy selections. Moreover, observations only allow us to measure redshift, intertwining galaxies' radial positions and velocities. In this paper we train and test our models on galaxy catalogs, created from thousands of state-of-the-art hydrodynamic simulations run with different codes from the CAMELS project, that incorporate these observational effects. We find that, although the presence of these effects degrades the precision and accuracy of the models, and increases the fraction of catalogs where the model breaks down, the fraction of galaxy catalogs where the model performs well is over 90 %, demonstrating the potential of these models to constrain cosmological parameters even when applied to real data. △ Less

Submitted 9 May, 2024; v1 submitted 23 October, 2023; originally announced October 2023.

Comments: 39 pages, 25 figures. For the reference in the abstract (de Santi et al. 2023) see arXiv:2302.14101

Evaluating residual acceleration noise for TianQin gravitational waves observatory with an empirical magnetic field model

Authors: Wei Su, Ze-Bing Zhou, Yan Wang, Chen Zhou, P. F. Chen, Wei Hong, J. H. Peng, Yun Yang, Y. W. Ni

Abstract: TianQin (TQ) project plans to deploy three satellites in space around the Earth to measure the displacement change of test masses caused by gravitational waves via laser interferometry. The requirement of the acceleration noise of the test mass is on the order of $10^{-15}~\,{\rm m}\,{\rm s}^{-2}\,{\rm Hz}^{-1/2}$ in the sensitive frequency range of TQ, %the extremely precise acceleration measurem… ▽ More TianQin (TQ) project plans to deploy three satellites in space around the Earth to measure the displacement change of test masses caused by gravitational waves via laser interferometry. The requirement of the acceleration noise of the test mass is on the order of $10^{-15}~\,{\rm m}\,{\rm s}^{-2}\,{\rm Hz}^{-1/2}$ in the sensitive frequency range of TQ, %the extremely precise acceleration measurement requirements make it necessary to investigate acceleration noise due to space magnetic fields. which is so stringent that the acceleration noise caused by the interaction of the space magnetic field with the test mass needs to be investigated. In this work, by using the Tsyganenko model, a data-based empirical space magnetic field model, we obtain the magnetic field distribution around TQ's orbit spanning two solar cycles in 23 years from 1998 to 2020. With the obtained space magnetic field, we derive the distribution and amplitude spectral densities (ASDs) of the acceleration noise of TQ in 23 years. Our results reveal that the average values of the ratio of the acceleration noise cauesd by the space magnetic field to the requirements of TQ at 1 mHz ($R_{\rm 1mHz}$) and 6 mHz ($R_{\rm 6mHz}$) are 0.123$\pm$0.052 and 0.027$\pm$0.013, respectively. The occurence probabilities of $R_{\rm 1mHz}>0.2$ and $>0.3$ are only 7.9% and 1.2%, respectively, and $R_{\rm 6mHz}$ never exceeds 0.2. △ Less

Submitted 30 November, 2023; v1 submitted 15 October, 2023; originally announced October 2023.

Comments: 13 pages, 10 figures, Published in PRD

Cosmology with multiple galaxies

Authors: Chaitanya Chawak, Francisco Villaescusa-Navarro, Nicolas Echeverri Rojas, Yueying Ni, ChangHoon Hahn, Daniel Angles-Alcazar

Abstract: Recent works have discovered a relatively tight correlation between $Ω_{\rm m}$ and properties of individual simulated galaxies. Because of this, it has been shown that constraints on $Ω_{\rm m}$ can be placed using the properties of individual galaxies while accounting for uncertainties on astrophysical processes such as feedback from supernova and active galactic nuclei. In this work, we quantif… ▽ More Recent works have discovered a relatively tight correlation between $Ω_{\rm m}$ and properties of individual simulated galaxies. Because of this, it has been shown that constraints on $Ω_{\rm m}$ can be placed using the properties of individual galaxies while accounting for uncertainties on astrophysical processes such as feedback from supernova and active galactic nuclei. In this work, we quantify whether using the properties of multiple galaxies simultaneously can tighten those constraints. For this, we train neural networks to perform likelihood-free inference on the value of two cosmological parameters ($Ω_{\rm m}$ and $σ_8$) and four astrophysical parameters using the properties of several galaxies from thousands of hydrodynamic simulations of the CAMELS project. We find that using properties of more than one galaxy increases the precision of the $Ω_{\rm m}$ inference. Furthermore, using multiple galaxies enables the inference of other parameters that were poorly constrained with one single galaxy. We show that the same subset of galaxy properties are responsible for the constraints on $Ω_{\rm m}$ from one and multiple galaxies. Finally, we quantify the robustness of the model and find that without identifying the model range of validity, the model does not perform well when tested on galaxies from other galaxy formation models. △ Less

Submitted 21 September, 2023; originally announced September 2023.

Comments: 13 pages, 7 figures

Modeling Quasar Proximity Zones in a Realistic Cosmological Environment with a Self-consistent Light Curve

Authors: Yihao Zhou, Huanqing Chen, Tiziana Di Matteo, Yueying Ni, Rupert A. C. Croft, Simeon Bird

Abstract: We study quasar proximity zones in a simulation that includes a self-consistent quasar formation model and realistic IGM environments. The quasar host halo is $10^{13}\ M_{\mathrm{\odot}}$ at $z=6$, more massive than typical halos studied in previous work. Between $6<z<7.5$, the quasar luminosity varies rapidly, with a mean magnitude of $M_{UV,mean}=-24.8$ and the fluctuation reaching up to two or… ▽ More We study quasar proximity zones in a simulation that includes a self-consistent quasar formation model and realistic IGM environments. The quasar host halo is $10^{13}\ M_{\mathrm{\odot}}$ at $z=6$, more massive than typical halos studied in previous work. Between $6<z<7.5$, the quasar luminosity varies rapidly, with a mean magnitude of $M_{UV,mean}=-24.8$ and the fluctuation reaching up to two orders of magnitude. Using this light curve to post-process the dense environment around the quasar, we find that the proximity zone size ($R_{p}$) ranges between $0.5-5$ pMpc. We show that the light curve variability causes a similar degree of scatter in $R_{p}$ as does the density fluctuation, both of which result in a standard deviation of $\sim 0.3$ pMpc). The $R_{p}$ traces the light curve fluctuations closely but with a time delay of $\sim 10^4\ \mathrm{yr}$, breaking the correspondence between the $R_{p}$ and the contemporaneous $M_{UV}$. This also indicates that we can only infer quasar activity within the past $\sim 10^4$ years instead of the integrated lifetime from $R_{p}$ in the later part of cosmic reionization. Compared with the variable light curve, a constant light curve underestimates the $R_{p}$ by 13% at the dim end ($M_{UV}\sim -23.5$), and overestimates the $R_{p}$ by 30% at the bright end ($M_{UV}\sim -26$). By calculating the $R_{p}$ generated by a number of quasars, we show that variable light curves predict a wider $R_{p}$ distribution than lightbulb models, and readily explain the extremely small $R_{p}$ values that have been observed. △ Less

Submitted 20 September, 2023; originally announced September 2023.

Comments: 15 pages, 12 figures, comments welcome

Extreme Tidal Stripping May Explain the Overmassive Black Hole in Leo I: a Proof of Concept

Authors: Fabio Pacucci, Yueying Ni, Abraham Loeb

Abstract: A recent study found dynamical evidence of a supermassive black hole of $\sim 3 \times 10^{6} \, \rm M_\odot$ at the center of Leo I, the most distant dwarf spheroidal galaxy of the Milky Way. This black hole, comparable in mass to the Milky Way's Sgr A*, places the system $>2$ orders of magnitude above the standard $M_\bullet-M_{\star}$ relation. We investigate the possibility, from a dynamical s… ▽ More A recent study found dynamical evidence of a supermassive black hole of $\sim 3 \times 10^{6} \, \rm M_\odot$ at the center of Leo I, the most distant dwarf spheroidal galaxy of the Milky Way. This black hole, comparable in mass to the Milky Way's Sgr A*, places the system $>2$ orders of magnitude above the standard $M_\bullet-M_{\star}$ relation. We investigate the possibility, from a dynamical standpoint, that Leo I's stellar system was originally much more massive and, thus, closer to the relation. Extreme tidal disruption from one or two close passages within the Milky Way's virial radius could have removed most of its stellar mass. A simple analytical model suggests that the progenitor of Leo I could have experienced a mass loss in the range of $32\% - 57\%$ from a single pericenter passage, depending on the stellar velocity dispersion estimate. This mass loss percentage increases to the range $66\% - 78\%$ if the pericenter occurs at the minimum distance allowed by current orbital reconstructions. Detailed N-body simulations show that the mass loss could reach $\sim 90\%$ with up to two passages, again with pericenter distances compatible with the minimum value allowed by Gaia data. Despite very significant uncertainties in the properties of Leo I, we reproduce its current position and velocity dispersion, as well as the final stellar mass enclosed in $1$ kpc ($\sim 5\times 10^6 \, \rm M_\odot$) within a factor $< 2$. The most recent tidal stream is directed along our line of sight toward Leo I, making it difficult to detect. Evidence from this extreme tidal disruption event could be present in current Gaia data in the form of extended tidal streams. △ Less

Submitted 26 September, 2023; v1 submitted 5 September, 2023; originally announced September 2023.

Comments: Accepted for publication in The Astrophysical Journal Letters. 9 pages, 4 figures. A video of the N-body simulation can be viewed here: https://www.fabiopacucci.com/resources/simulations/

Journal ref: The Astrophysical Journal Letters, 2023, Volume 956, Number 2

Data-constrained Magnetohydrodynamic Simulation of an Intermediate Solar Filament Eruption

Authors: Yang Guo, Jinhan Guo, Yiwei Ni, M. D. Ding, P. F. Chen, Chun Xia, Rony Keppens, Kai E. Yang

Abstract: Solar eruptive activities could occur in weak magnetic field environments and over large spatial scales, especially relevant to eruptions involving intermediate or quiescent solar filaments. To handle the large scales, we implement and apply a flux rope embedding method using regularized Biot-Savart laws in the spherical coordinate system. Combined with a potential field source surface model and a… ▽ More Solar eruptive activities could occur in weak magnetic field environments and over large spatial scales, especially relevant to eruptions involving intermediate or quiescent solar filaments. To handle the large scales, we implement and apply a flux rope embedding method using regularized Biot-Savart laws in the spherical coordinate system. Combined with a potential field source surface model and a magneto-frictional method, a nonlinear force-free field comprising a flux rope embedded in a potential field is constructed. Using the combined nonlinear force-free field as the initial condition, we then perform a zero-$β$ data-constrained magnetohydrodynamic (MHD) simulation for an M8.7 flare at 03:38 UT on 2012 January 23. The MHD model reproduces the eruption process, flare ribbon evolution (represented by the quasi-separatrix layer evolution) and kinematics of the flux rope. This approach could potentially model global-scale eruptions from weak field regions. △ Less

Submitted 3 September, 2023; originally announced September 2023.

Comments: 23 pages, 7 figures, accepted for publicaiton in ApJ

Understanding the Lateral Drifting of an Erupting Filament with a Data-constrained Magnetohydrodynamic Simulation

Authors: Jinhan Guo, Ye Qiu, Yiwei Ni, Yang Guo, Chuan Li, Yuhang Gao, Brigitte Schmieder, Stefaan Poedts, Pengfei Chen

Abstract: Solar filaments often exhibit rotation and deflection during eruptions, which would significantly affect the geoeffectiveness of the corresponding coronal mass ejections (CMEs). Therefore, understanding the mechanisms that lead to such rotation and lateral displacement of filaments is a great concern to space weather forecasting. In this paper, we examine an intriguing filament eruption event obse… ▽ More Solar filaments often exhibit rotation and deflection during eruptions, which would significantly affect the geoeffectiveness of the corresponding coronal mass ejections (CMEs). Therefore, understanding the mechanisms that lead to such rotation and lateral displacement of filaments is a great concern to space weather forecasting. In this paper, we examine an intriguing filament eruption event observed by the Chinese Hα Solar Explorer (CHASE) and the Solar Dynamics Observatory (SDO). The filament, which eventually evolves into a CME, exhibits significant lateral drifting during its rising. Moreover, the orientation of the CME flux rope axis deviates from that of the pre-eruptive filament observed in the source region. To investigate the physical processes behind these observations, we perform a data-constrained magnetohydrodynamic (MHD) simulation. Many prominent observational features in the eruption are reproduced by our numerical model, including the morphology of the eruptive filament, eruption path, and flare ribbons. The simulation results reveal that the magnetic reconnection between the flux-rope leg and neighboring low-lying sheared arcades may be the primary mechanism responsible for the lateral drifting of the filament material. Such a reconnection geometry leads to flux-rope footpoint migration and a reconfiguration of its morphology. As a consequence, the filament material hosted in the flux rope drifts laterally, and the CME flux rope deviates from the pre-eruptive filament. This finding underscores the importance of external magnetic reconnection in influencing the orientation of a flux rope axis during eruption. △ Less

Submitted 17 August, 2023; originally announced August 2023.

Comments: 25 pages, 14 figures, Accepted for publication in ApJ

Cosmological baryon spread and impact on matter clustering in CAMELS

Authors: Matthew Gebhardt, Daniel Anglés-Alcázar, Josh Borrow, Shy Genel, Francisco Villaescusa-Navarro, Yueying Ni, Christopher Lovell, Daisuke Nagai, Romeel Davé, Federico Marinacci, Mark Vogelsberger, Lars Hernquist

Abstract: We quantify the cosmological spread of baryons relative to their initial neighboring dark matter distribution using thousands of state-of-the-art simulations from the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project. We show that dark matter particles spread relative to their initial neighboring distribution owing to chaotic gravitational dynamics on spatial scales com… ▽ More We quantify the cosmological spread of baryons relative to their initial neighboring dark matter distribution using thousands of state-of-the-art simulations from the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project. We show that dark matter particles spread relative to their initial neighboring distribution owing to chaotic gravitational dynamics on spatial scales comparable to their host dark matter halo. In contrast, gas in hydrodynamic simulations spreads much further from the initial neighboring dark matter owing to feedback from supernovae (SNe) and Active Galactic Nuclei (AGN). We show that large-scale baryon spread is very sensitive to model implementation details, with the fiducial \textsc{SIMBA} model spreading $\sim$40\% of baryons $>$1\,Mpc away compared to $\sim$10\% for the IllustrisTNG and \textsc{ASTRID} models. Increasing the efficiency of AGN-driven outflows greatly increases baryon spread while increasing the strength of SNe-driven winds can decrease spreading due to non-linear coupling of stellar and AGN feedback. We compare total matter power spectra between hydrodynamic and paired $N$-body simulations and demonstrate that the baryonic spread metric broadly captures the global impact of feedback on matter clustering over variations of cosmological and astrophysical parameters, initial conditions, and galaxy formation models. Using symbolic regression, we find a function that reproduces the suppression of power by feedback as a function of wave number ($k$) and baryonic spread up to $k \sim 10\,h$\,Mpc$^{-1}$ while highlighting the challenge of developing models robust to variations in galaxy formation physics implementation. △ Less

Submitted 21 July, 2023; originally announced July 2023.

Comments: 17 pages, 15 figures

Fly-by galaxy encounters with multiple black holes produce star-forming linear wakes

Authors: Nianyi Chen, Patrick LaChance, Yueying Ni, Tiziana Di Matteo, Rupert Croft, Priyamvada Natarajan, Simeon Bird

Abstract: We look for simulated star-forming linear wakes such as the one recently discovered by van Dokkum et al. (2023) in the cosmological hydrodynamical simulation ASTRID. Amongst the runaway black holes in ASTRID, none are able to produce clear star-forming wakes. Meanwhile, fly-by encounters, typically involving a compact galaxy (with a central black hole) and a star-forming galaxy (with a duo of blac… ▽ More We look for simulated star-forming linear wakes such as the one recently discovered by van Dokkum et al. (2023) in the cosmological hydrodynamical simulation ASTRID. Amongst the runaway black holes in ASTRID, none are able to produce clear star-forming wakes. Meanwhile, fly-by encounters, typically involving a compact galaxy (with a central black hole) and a star-forming galaxy (with a duo of black holes) reproduce remarkably well many of the key properties (its length and linearity; recent star formation, etc.) of the observed star-forming linear feature. We predict the feature to persist for approximately 100 Myr in such a system and hence constitute a rare event. The feature contains a partly stripped galaxy (with $M_{\rm gal}=10^9 \sim 10^{10}M_\odot$) and a dual BH system ($M_{\rm BH}=10^5 \sim 10^7\,M_\odot$) in its brightest knot. X-ray emission from AGN in the knot should be detectable in such systems. After $100\sim 200\,{\rm Myrs}$ from the first fly-by, the galaxies merge leaving behind a triple black hole system in a (still) actively star-forming early-type remnant of mass $\sim 5\times 10^{10}\,M_\odot$. Follow-up JWST observations may be key for revealing the nature of these linear features by potentially detecting the older stellar populations constituting the bright knot. Confirmation of such detections may therefore help discriminate a fly-by encounter from a massive BH wake to reveal the origin of such features. △ Less

Submitted 3 July, 2023; originally announced July 2023.

Comments: 8 pages, 5 figures, comments welcome

Efficient Reionization in a Large Hydrodynamic Galaxy Formation Simulation

Authors: James E. Davies, Simeon Bird, Simon Mutch, Yueying Ni, Yu Feng, Rupert Croft, Tiziana Di Matteo, J. Stuart B. Wyithe

Abstract: Accuracy in the topology and statistics of a simulated Epoch of Reionization (EoR) are vital to draw connections between observations and physical processes. While full radiative transfer models produce the most accurate reionization models, they are highly computationally expensive, and are infeasible for the largest cosmological simulations. Instead, large simulations often include EoR models th… ▽ More Accuracy in the topology and statistics of a simulated Epoch of Reionization (EoR) are vital to draw connections between observations and physical processes. While full radiative transfer models produce the most accurate reionization models, they are highly computationally expensive, and are infeasible for the largest cosmological simulations. Instead, large simulations often include EoR models that are pre-computed via the initial density field, or post-processed where feedback effects are ignored. We introduce Astrid-ES, a resimulation of the Astrid epoch of reionisation $20 > z > 5.5$ which includes an on-the-fly excursion-set reionization algorithm. Astrid-ES produces more accurate reionization histories without significantly impacting the computational time. This model directly utilises the star particles produced in the simulation to calculate the EoR history and includes a UV background which heats the gas particles after their reionization. We contrast the reionization topology and statistics in Astrid-ES with the previously employed parametric reionisation model, finding that in Astrid-ES, ionised regions are more correlated with galaxies, and the 21cm power-spectrum shows an increase in large scale power. We calculate the relation between the size of HII regions and the UV luminosity of the brightest galaxy within them. Prior to the overlap phase, we find a power-law fit of $\mathrm{log} (R) = -0.314 M_\mathrm{UV} - 2.550 \mathrm{log}(1+z) + 7.408$ with a standard deviation $σ_R < 0.15 \mathrm{dex}$ across all mass bins. We also examine the properties of halos throughout reionization, finding that while the properties of halos in the simulation are correlated with the redshift of reionisation, they are not greatly affected by reionisation itself. △ Less

Submitted 8 August, 2023; v1 submitted 13 June, 2023; originally announced June 2023.

Comments: 12 pages, 11 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society

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

Authors: Simeon Bird, Martin Fernandez, Ming-Feng Ho, Mahdi Qezlou, Reza Monadi, Yueying Ni, Nianyi Chen, Rupert Croft, Tiziana Di Matteo

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

Submitted 2 January, 2024; v1 submitted 8 June, 2023; originally announced June 2023.

Comments: 25 pages, 11 figures, matches version published in JCAP

Journal ref: JCAP 10 (2023) 037

AI-assisted super-resolution cosmological simulations III: Time evolution

Authors: Xiaowen Zhang, Patrick Lachance, Yueying Ni, Yin Li, Rupert A. C. Croft, Tiziana Di Matteo, Simeon Bird, Yu Feng

Abstract: In this work, we extend our recently developed super-resolution (SR) model for cosmological simulations to produce fully time consistent evolving representations of the particle phase-space distribution. We employ a style-based constrained generative adversarial network (Style-GAN) where the changing cosmic time is an input style parameter to the network. The matter power spectrum and halo mass fu… ▽ More In this work, we extend our recently developed super-resolution (SR) model for cosmological simulations to produce fully time consistent evolving representations of the particle phase-space distribution. We employ a style-based constrained generative adversarial network (Style-GAN) where the changing cosmic time is an input style parameter to the network. The matter power spectrum and halo mass function agree well with results from high-resolution N-body simulations over the full trained redshift range ($10 \le z \le 0$). Furthermore, we assess the temporal consistency of our SR model by constructing halo merger trees. We examine progenitors, descendants and mass growth along the tree branches. All statistical indicators demonstrate the ability of our SR model to generate satisfactory high-resolution simulations based on low-resolution inputs. △ Less

Submitted 20 May, 2023; originally announced May 2023.

Comments: 12 pages, 11 figures, code and movie available in https://github.com/sagasv5-xw/map2map on styled srsgan branch

Cosmology with one galaxy? -- The ASTRID model and robustness

Authors: Nicolas Echeverri, Francisco Villaescusa-Navarro, Chaitanya Chawak, Yueying Ni, ChangHoon Hahn, Elena Hernandez-Martinez, Romain Teyssier, Daniel Angles-Alcazar, Klaus Dolag, Tiago Castro

Abstract: Recent work has pointed out the potential existence of a tight relation between the cosmological parameter $Ω_{\rm m}$, at fixed $Ω_{\rm b}$, and the properties of individual galaxies in state-of-the-art cosmological hydrodynamic simulations. In this paper, we investigate whether such a relation also holds for galaxies from simulations run with a different code that made use of a distinct subgrid… ▽ More Recent work has pointed out the potential existence of a tight relation between the cosmological parameter $Ω_{\rm m}$, at fixed $Ω_{\rm b}$, and the properties of individual galaxies in state-of-the-art cosmological hydrodynamic simulations. In this paper, we investigate whether such a relation also holds for galaxies from simulations run with a different code that made use of a distinct subgrid physics: Astrid. We find that also in this case, neural networks are able to infer the value of $Ω_{\rm m}$ with a $\sim10\%$ precision from the properties of individual galaxies while accounting for astrophysics uncertainties as modeled in CAMELS. This tight relationship is present at all considered redshifts, $z\leq3$, and the stellar mass, the stellar metallicity, and the maximum circular velocity are among the most important galaxy properties behind the relation. In order to use this method with real galaxies, one needs to quantify its robustness: the accuracy of the model when tested on galaxies generated by codes different from the one used for training. We quantify the robustness of the models by testing them on galaxies from four different codes: IllustrisTNG, SIMBA, Astrid, and Magneticum. We show that the models perform well on a large fraction of the galaxies, but fail dramatically on a small fraction of them. Removing these outliers significantly improves the accuracy of the models across simulation codes. △ Less

Submitted 12 April, 2023; originally announced April 2023.

Comments: 16 pages, 12 figures

The CAMELS project: Expanding the galaxy formation model space with new ASTRID and 28-parameter TNG and SIMBA suites

Authors: Yueying Ni, Shy Genel, Daniel Anglés-Alcázar, Francisco Villaescusa-Navarro, Yongseok Jo, Simeon Bird, Tiziana Di Matteo, Rupert Croft, Nianyi Chen, Natalí S. M. de Santi, Matthew Gebhardt, Helen Shao, Shivam Pandey, Lars Hernquist, Romeel Dave

Abstract: We present CAMELS-ASTRID, the third suite of hydrodynamical simulations in the Cosmology and Astrophysics with MachinE Learning (CAMELS) project, along with new simulation sets that extend the model parameter space based on the previous frameworks of CAMELS-TNG and CAMELS-SIMBA, to provide broader training sets and testing grounds for machine-learning algorithms designed for cosmological studies.… ▽ More We present CAMELS-ASTRID, the third suite of hydrodynamical simulations in the Cosmology and Astrophysics with MachinE Learning (CAMELS) project, along with new simulation sets that extend the model parameter space based on the previous frameworks of CAMELS-TNG and CAMELS-SIMBA, to provide broader training sets and testing grounds for machine-learning algorithms designed for cosmological studies. CAMELS-ASTRID employs the galaxy formation model following the ASTRID simulation and contains 2,124 hydrodynamic simulation runs that vary 3 cosmological parameters ($Ω_m$, $σ_8$, $Ω_b$) and 4 parameters controlling stellar and AGN feedback. Compared to the existing TNG and SIMBA simulation suites in CAMELS, the fiducial model of ASTRID features the mildest AGN feedback and predicts the least baryonic effect on the matter power spectrum. The training set of ASTRID covers a broader variation in the galaxy populations and the baryonic impact on the matter power spectrum compared to its TNG and SIMBA counterparts, which can make machine-learning models trained on the ASTRID suite exhibit better extrapolation performance when tested on other hydrodynamic simulation sets. We also introduce extension simulation sets in CAMELS that widely explore 28 parameters in the TNG and SIMBA models, demonstrating the enormity of the overall galaxy formation model parameter space and the complex non-linear interplay between cosmology and astrophysical processes. With the new simulation suites, we show that building robust machine-learning models favors training and testing on the largest possible diversity of galaxy formation models. We also demonstrate that it is possible to train accurate neural networks to infer cosmological parameters using the high-dimensional TNG-SB28 simulation set. △ Less

Submitted 4 April, 2023; originally announced April 2023.

Origin of high-velocity ejecta and early red excess emission in the infant Type Ia supernova 2021aefx

Authors: Yuan Qi Ni, Dae-Sik Moon, Maria R. Drout, Christopher D. Matzner, Kelvin C. C. Leong, Sang Chul Kim, Hong Soo Park, Youngdae Lee

Abstract: \object{SN 2021aefx} is a normal Type Ia Supernova (SN) with red excess emission over the first $\sim$ 2 days. We present detailed analysis of this SN using our high-cadence KMTNet multi-band photometry, spectroscopy, and publicly available data. We provide the first measurements of its epochs of explosion (MJD 59529.32 $\pm$ 0.16) as well as ``first light'' (MJD 59529.85 $\pm$ 0.55) associated wi… ▽ More \object{SN 2021aefx} is a normal Type Ia Supernova (SN) with red excess emission over the first $\sim$ 2 days. We present detailed analysis of this SN using our high-cadence KMTNet multi-band photometry, spectroscopy, and publicly available data. We provide the first measurements of its epochs of explosion (MJD 59529.32 $\pm$ 0.16) as well as ``first light'' (MJD 59529.85 $\pm$ 0.55) associated with the main ejecta ${\rm{^{56}Ni}}$ distribution. This places our first detection of SN 2021aefx at $\sim -$0.5 hours since ``first light'', indicating the presence of additional power sources. Our peak-spectrum confirms its Type Ia sub-classification as intermediate between Core-Normal and Broad-Line, and we estimate the ejecta mass to be $\sim$ 1.34 $M_{\odot}$. The pre-peak spectral evolution identifies fast-expanding material reaching $>$ 40,000 km s$^{-1}$ (the fastest ever observed in Type Ia SNe) and at least two distinct homologously-expanding ejecta components: (1) a normal-velocity (12,400 km s$^{-1}$) component consistent with the typical photospheric evolution of Chandrasekhar-mass ejecta; and (2) a high-velocity (23,500 km s$^{-1}$) component visible during the first $\sim$ 3.6 days post-explosion, which locates the component within the outer $<$ 16\% of the ejecta mass. Asymmetric, subsonic explosion processes producing a non-spherical photosphere provide an explanation for the simultaneous presence of the two components, as well as the red excess emission via a slight ${\rm{^{56}Ni}}$ enrichment in the outer $\sim$ 0.5\% of the ejecta mass. Our spectrum from 300 days post-peak advances the constraint against non-degenerate companions and further supports a near-Chandrasekhar-mass explosion origin. Off-center ignited delayed-detonations of Chandrasekhar-mass white dwarfs may be responsible for the observed features of SN 2021aefx in some normal Type Ia SNe. △ Less

Submitted 5 April, 2023; v1 submitted 2 April, 2023; originally announced April 2023.

Comments: Submitted for publication in ApJ. 29 pages, 14 figures, 4 tables

Boosting Line Intensity Map Signal-to-Noise with the Ly-$α$ Forest Cross-Correlation

Authors: Mahdi Qezlou, Simeon Bird, Adam Lidz, Guochao Sun, Andrew B. Newman, Gwen C. Rudie, Yueying Ni, Rupert Croft, Tiziana Di Matteo

Abstract: We forecast the prospects for cross-correlating future line intensity mapping (LIM) surveys with the current and future Ly-$α$ forest data. We use large cosmological hydrodynamic simulations to model the expected emission signal for the CO rotational transition in the COMAP LIM experiment at the 5-year benchmark and the Ly-$α$ forest absorption signal for various surveys, including eBOSS, DESI, an… ▽ More We forecast the prospects for cross-correlating future line intensity mapping (LIM) surveys with the current and future Ly-$α$ forest data. We use large cosmological hydrodynamic simulations to model the expected emission signal for the CO rotational transition in the COMAP LIM experiment at the 5-year benchmark and the Ly-$α$ forest absorption signal for various surveys, including eBOSS, DESI, and PFS. We show that CO$\times$Ly-$α$ forest can significantly enhance the detection signal-to-noise ratio of CO, with a $200$ to $300 \%$ improvement when cross-correlated with the forest observed in the Prime Focus Spectrograph (PFS) survey and a $50$ to $75\%$ enhancement for the currently available eBOSS or the upcoming DESI observations. We compare to the signal-to-noise improvements expected for a galaxy survey and show that CO$\times$Ly-$α$ is competitive with even a spectroscopic galaxy survey in raw signal-to-noise. Furthermore, our study suggests that the clustering of CO emission is tightly constrained by CO$\times$Ly-$α$ forest, due to the increased signal-to-noise ratio and the simplicity of Ly-$α$ absorption power spectrum modeling. Any foreground contamination or systematics are expected not to be shared between LIM surveys and Ly-$α$ forest observations; this provides an unbiased inference. Our findings highlight the potential benefits of utilizing the Ly-$α$ forest to aid in the initial detection of signals in line intensity experiments. For example, we also estimate that [CII]$\times$Ly-$α$ forest measurements from EXCLAIM and DESI/eBOSS, respectively, should have a larger signal-to-noise ratio than planned [CII]$\times$quasar observations by about an order of magnitude. Our results can be readily applied to actual data thanks to the observed quasar spectra in eBOSS Stripe 82, which overlaps with several LIM surveys. △ Less

Submitted 1 November, 2023; v1 submitted 30 March, 2023; originally announced March 2023.

Comments: Codes and the produced data are available at https://github.com/qezlou/lali

Dwarf Galaxy Discoveries from the KMTNet Supernova Program III. the Milky-Way Analog NGC~2997 Group

Authors: Tony Junjing Fan, Dae-Sik Moon, Hong Soo Park, Dennis Zaritsky, Sang Chul Kim, Youngdae Lee, Ting S. Li, Yuan Qi Ni, Jeehye Shin, Sang-Mok Cha, Yongseok Lee

Abstract: We present the discovery of 48 new and the analysis of 55, including 7 previously discovered, dwarf galaxy candidates around the giant spiral galaxy NGC~2997 using deep $BVI$ images from the KMTNet Supernova Program. Their $V$-band central surface brightness and total absolute magnitudes are in the range of 20.3--26.7 mag arcsec$^{-2}$ and --(8.02--17.69) mag, respectively, while the $I$-band effe… ▽ More We present the discovery of 48 new and the analysis of 55, including 7 previously discovered, dwarf galaxy candidates around the giant spiral galaxy NGC~2997 using deep $BVI$ images from the KMTNet Supernova Program. Their $V$-band central surface brightness and total absolute magnitudes are in the range of 20.3--26.7 mag arcsec$^{-2}$ and --(8.02--17.69) mag, respectively, while the $I$-band effective radii are 0.14--2.97 kpc. We obtain $α$ $\simeq$ --1.43 $\pm$ 0.02 for the faint-end slope of their luminosity function, comparable to previously measured values but shallower than theoretical predictions based on $Λ$CDM models. The distance-independent distributions of their mass and color suggest that the group could have recently accreted new massive members from the surrounding fields. The systematically bluer colors of the brighter members indicate younger stellar population and higher star formation activities in them, which appears to be consistent with similar findings from the SAGA or ELVES survey. We suggest that the massive and bluer dwarf galaxies in the group have experienced less environmental quenching due to their recent accretion, while environmental quenching is more effective for the low-mass members. The interpretation of NGC~2997 being populationally young with recent accretion of massive members is also consistent with the overall morphological distribution of the dwarf galaxies showing a lack of morphologically evolved candidates but a plethora of irregularly shaped ones. Our detection rate of dwarf galaxy candidates in the NGC~2997 group and their inferred star formation activities are comparable to those found in Milky Way analog systems from recent surveys within the magnitude limit M$_{V}$ $\lesssim$ --13 mag. △ Less

Submitted 26 July, 2023; v1 submitted 29 March, 2023; originally announced March 2023.

Comments: 24 pages, 14 figures

Thermodynamic and Magnetic Topology Evolution of the X1.0 Flare on 2021 October 28 Simulated by a Data-driven Radiative Magnetohydrodynamic Model

Authors: Jin-han Guo, Yi-wei Ni, Ze Zhong, Yang Guo, Chun Xia, Hai-tang Li, Stefaan Poedts, Brigitte Schmieder, Peng-fei Chen

Abstract: Solar filament eruptions, flares and coronal mass ejections (CMEs) are manifestations of drastic release of energy in the magnetic field, which are related to many eruptive phenomena from the Earth magnetosphere to black hole accretion disks. With the availability of high-resolution magnetograms on the solar surface, observational data-based modelling is a promising way to quantitatively study the… ▽ More Solar filament eruptions, flares and coronal mass ejections (CMEs) are manifestations of drastic release of energy in the magnetic field, which are related to many eruptive phenomena from the Earth magnetosphere to black hole accretion disks. With the availability of high-resolution magnetograms on the solar surface, observational data-based modelling is a promising way to quantitatively study the underlying physical mechanisms behind observations. By incorporating thermal conduction and radiation losses in the energy equation, we develop a new data-driven radiative magnetohydrodynamic (MHD) model, which has the capability to capture the thermodynamic evolution compared to our previous zero-\b{eta} model. Our numerical results reproduce major observational characteristics of the X1.0 flare on 2021 October 28 in NOAA active region (AR) 12887, including the morphology of the eruption, kinematic of flare ribbons, extreme-ultraviolet (EUV) radiations, and two components of the EUV waves predicted by the magnetic stretching model, i.e., a fast-mode shock wave and a slower apparent wave due to successive stretching of magnetic field lines. Moreover, some intriguing phenomena are revealed in the simulation. We find that flare ribbons separate initially and ultimately stop at the outer stationary quasi-separatrix layers (QSLs). Such outer QSLs correspond to the border of the filament channel and determine the final positions of flare ribbons, which can be used to predict the size and the lifetime of a flare before it occurs. In addition, the side view of the synthesized EUV and white-light images exhibit typical three-part structures of CMEs, where the bright leading front is roughly cospatial with the non-wave component of the EUV wave, reinforcing the magnetic stretching model for the slow component of EUV waves. △ Less

Submitted 24 March, 2023; originally announced March 2023.

Comments: 46 pages, 15 figures, accepted for publication in The Astrophysical Journal Supplement

Triple and Quadruple Black Holes in the ASTRID Simulation at $z \sim 2$

Authors: Calvin Hoffman, Nianyi Chen, Tiziana Di Matteo, Yueying Ni, Simeon Bird, Rupert Croft, Abraham Loeb

Abstract: We use the ASTRID cosmological hydrodynamic simulation to investigate the properties and evolution of triple and quadruple Massive Black Hole (MBH) systems at $z = 2-3$. Only a handful of MBH tuple systems have been detected to date. In ASTRID, we find $4\%$ of the $M_{\rm BH}>10^7\,M_\odot$ are in tuples with $Δr_{\rm max} < 200\,{\rm kpc}$. The tuple systems span a range of separations with the… ▽ More We use the ASTRID cosmological hydrodynamic simulation to investigate the properties and evolution of triple and quadruple Massive Black Hole (MBH) systems at $z = 2-3$. Only a handful of MBH tuple systems have been detected to date. In ASTRID, we find $4\%$ of the $M_{\rm BH}>10^7\,M_\odot$ are in tuples with $Δr_{\rm max} < 200\,{\rm kpc}$. The tuple systems span a range of separations with the majority of the observable AGN systems at $Δr \sim 50-100$ kpc. They include some of the most massive BHs (up to $10^{10} \,M_\odot$) but with at least one of the components of $M_{\rm BH} \sim 10^7 \,M_\odot$. Tuples' host galaxies are typically massive with $M_* \sim 10^{10-11} \,M_\odot$. We find that $>10\%$ massive halos with $M_{\rm halo} > 10^{13} M_\odot$ host MBH tuples. Following the subsequent interactions between MBHs in tuples, we found that in $\sim 5\%$ of the triplets all three MBHs merge within a Gyr, and $15\%$ go through one merger. As a by-product of the complex multi-galaxy interaction of these systems, we also find that up to $\sim 5\%$ of tuples lead to runaway MBHs. In ASTRID, virtually all of the ultramassive black holes ($>10^{10} \,M_\odot $) have undergone a triple quasar phase while for BHs with $M_{\rm BH} \sim 10^9 \,M_\odot$ this fraction drops to $50\%$. △ Less

Submitted 8 March, 2023; originally announced March 2023.

Comments: 10 pages, 9 figures; comments welcome

A universal equation to predict $Ω_{\rm m}$ from halo and galaxy catalogues

Authors: Helen Shao, Natalí S. M de Santi, Francisco Villaescusa-Navarro, Romain Teyssier, Yueying Ni, Daniel Angles-Alcazar, Shy Genel, Lars Hernquist, Ulrich P. Steinwandel, Tiago Castro, Elena Hernandez-Martınez, Klaus Dolag, Christopher C. Lovell, Eli Visbal, Lehman H. Garrison, Mihir Kulkarni

Abstract: We discover analytic equations that can infer the value of $Ω_{\rm m}$ from the positions and velocity moduli of halo and galaxy catalogues. The equations are derived by combining a tailored graph neural network (GNN) architecture with symbolic regression. We first train the GNN on dark matter halos from Gadget N-body simulations to perform field-level likelihood-free inference, and show that our… ▽ More We discover analytic equations that can infer the value of $Ω_{\rm m}$ from the positions and velocity moduli of halo and galaxy catalogues. The equations are derived by combining a tailored graph neural network (GNN) architecture with symbolic regression. We first train the GNN on dark matter halos from Gadget N-body simulations to perform field-level likelihood-free inference, and show that our model can infer $Ω_{\rm m}$ with $\sim6\%$ accuracy from halo catalogues of thousands of N-body simulations run with six different codes: Abacus, CUBEP$^3$M, Gadget, Enzo, PKDGrav3, and Ramses. By applying symbolic regression to the different parts comprising the GNN, we derive equations that can predict $Ω_{\rm m}$ from halo catalogues of simulations run with all of the above codes with accuracies similar to those of the GNN. We show that by tuning a single free parameter, our equations can also infer the value of $Ω_{\rm m}$ from galaxy catalogues of thousands of state-of-the-art hydrodynamic simulations of the CAMELS project, each with a different astrophysics model, run with five distinct codes that employ different subgrid physics: IllustrisTNG, SIMBA, Astrid, Magneticum, SWIFT-EAGLE. Furthermore, the equations also perform well when tested on galaxy catalogues from simulations covering a vast region in parameter space that samples variations in 5 cosmological and 23 astrophysical parameters. We speculate that the equations may reflect the existence of a fundamental physics relation between the phase-space distribution of generic tracers and $Ω_{\rm m}$, one that is not affected by galaxy formation physics down to scales as small as $10~h^{-1}{\rm kpc}$. △ Less

Submitted 28 February, 2023; originally announced February 2023.

Comments: 32 pages, 13 figures, summary video: https://youtu.be/STZHvDHkVgo

Robust Field-level Likelihood-free Inference with Galaxies

Authors: Natalí S. M. de Santi, Helen Shao, Francisco Villaescusa-Navarro, L. Raul Abramo, Romain Teyssier, Pablo Villanueva-Domingo, Yueying Ni, Daniel Anglés-Alcázar, Shy Genel, Elena Hernandez-Martinez, Ulrich P. Steinwandel, Christopher C. Lovell, Klaus Dolag, Tiago Castro, Mark Vogelsberger

Abstract: We train graph neural networks to perform field-level likelihood-free inference using galaxy catalogs from state-of-the-art hydrodynamic simulations of the CAMELS project. Our models are rotational, translational, and permutation invariant and do not impose any cut on scale. From galaxy catalogs that only contain $3$D positions and radial velocities of $\sim 1, 000$ galaxies in tiny… ▽ More We train graph neural networks to perform field-level likelihood-free inference using galaxy catalogs from state-of-the-art hydrodynamic simulations of the CAMELS project. Our models are rotational, translational, and permutation invariant and do not impose any cut on scale. From galaxy catalogs that only contain $3$D positions and radial velocities of $\sim 1, 000$ galaxies in tiny $(25~h^{-1}{\rm Mpc})^3$ volumes our models can infer the value of $Ω_{\rm m}$ with approximately $12$ % precision. More importantly, by testing the models on galaxy catalogs from thousands of hydrodynamic simulations, each having a different efficiency of supernova and AGN feedback, run with five different codes and subgrid models - IllustrisTNG, SIMBA, Astrid, Magneticum, SWIFT-EAGLE -, we find that our models are robust to changes in astrophysics, subgrid physics, and subhalo/galaxy finder. Furthermore, we test our models on $1,024$ simulations that cover a vast region in parameter space - variations in $5$ cosmological and $23$ astrophysical parameters - finding that the model extrapolates really well. Our results indicate that the key to building a robust model is the use of both galaxy positions and velocities, suggesting that the network have likely learned an underlying physical relation that does not depend on galaxy formation and is valid on scales larger than $\sim10~h^{-1}{\rm kpc}$. △ Less

Submitted 18 July, 2023; v1 submitted 27 February, 2023; originally announced February 2023.

Comments: 34 pages, 12 figures. For a video summarizing the results, see https://youtu.be/b59ep7cyPOs

Journal ref: Volume 952, Number 1, Year 2023, Pages 69

High-redshift supermassive black hole mergers in simulations with dynamical friction modelling

Authors: Colin DeGraf, Nianyi Chen, Yueying Ni, Tiziana Di Matteo, Simeon Bird, Michael Tremmel, Rupert Croft

Abstract: In the near future, projects like LISA and Pulsar Timing Arrays are expected to detect gravitational waves from mergers between supermassive black holes, and it is crucial to precisely model the underlying merger populations now to maximize what we can learn from this new data. Here we characterize expected high-redshift (z > 2) black hole mergers using the very large volume Astrid cosmological si… ▽ More In the near future, projects like LISA and Pulsar Timing Arrays are expected to detect gravitational waves from mergers between supermassive black holes, and it is crucial to precisely model the underlying merger populations now to maximize what we can learn from this new data. Here we characterize expected high-redshift (z > 2) black hole mergers using the very large volume Astrid cosmological simulation, which uses a range of seed masses to probe down to low-mass BHs, and directly incorporates dynamical friction so as to accurately model the dynamical processes which bring black holes to the galaxy center where binary formation and coalescence will occur. The black hole populations in Astrid include black holes down to 10$^{4.5}$ M$_\odot$, and remain broadly consistent with the TNG simulations at scales > 10$^6$ M$_\odot$ (the seed mass used in TNG). By resolving lower-mass black holes, the overall merger rate is ~5x higher than in TNG. However, incorporating dynamical friction delays mergers compared to a recentering scheme, reducing the high-z merger rate mass-matched mergers by a factor of ~2x. We also calculate the expected LISA Signal-to-Noise values, and show that the distribution peaks at high SNR (>100), emphasizing the importance of implementing a seed mass well below LISA's peak sensitivity (10$^6$ M$_\odot$) to resolve the majority of LISA's GW detections. △ Less

Submitted 1 February, 2023; originally announced February 2023.

Comments: 11 pages, 10 figures, submitted to MNRAS

Inferring the impact of feedback on the matter distribution using the Sunyaev Zel'dovich effect: Insights from CAMELS simulations and ACT+DES data

Authors: Shivam Pandey, Kai Lehman, Eric J. Baxter, Yueying Ni, Daniel Anglés-Alcázar, Shy Genel, Francisco Villaescusa-Navarro, Ana Maria Delgado, Tiziana di Matteo

Abstract: Feedback from active galactic nuclei and stellar processes changes the matter distribution on small scales, leading to significant systematic uncertainty in weak lensing constraints on cosmology. We investigate how the observable properties of group-scale halos can constrain feedback's impact on the matter distribution using Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS). Ex… ▽ More Feedback from active galactic nuclei and stellar processes changes the matter distribution on small scales, leading to significant systematic uncertainty in weak lensing constraints on cosmology. We investigate how the observable properties of group-scale halos can constrain feedback's impact on the matter distribution using Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS). Extending the results of previous work to smaller halo masses and higher wavenumber, $k$, we find that the baryon fraction in halos contains significant information about the impact of feedback on the matter power spectrum. We explore how the thermal Sunyaev Zel'dovich (tSZ) signal from group-scale halos contains similar information. Using recent Dark Energy Survey (DES) weak lensing and Atacama Cosmology Telescope (ACT) tSZ cross-correlation measurements and models trained on CAMELS, we obtain $10\%$ constraints on feedback effects on the power spectrum at $k \sim 5\, h/{\rm Mpc}$. We show that with future surveys, it will be possible to constrain baryonic effects on the power spectrum to $\mathcal{O}(<1\%)$ at $k = 1\, h/{\rm Mpc}$ and $\mathcal{O}(3\%)$ at $k = 5\, h/{\rm Mpc}$ using the methods that we introduce here. Finally, we investigate the impact of feedback on the matter bispectrum, finding that tSZ observables are highly informative in this case. △ Less

Submitted 5 January, 2023; originally announced January 2023.

Comments: 18 pages, 15 figures. Comments are welcome

Orbital and Radiative Properties of Wandering Intermediate-Mass Black Holes in the ASTRID Simulation

Authors: Emma Jane Weller, Fabio Pacucci, Yueying Ni, Nianyi Chen, Tiziana Di Matteo, Magdalena Siwek, Lars Hernquist

Abstract: Intermediate-Mass Black Holes (IMBHs) of $10^3-10^6 \, M_\odot$ are commonly found at the center of dwarf galaxies. Simulations and observations convincingly show that a sizable population of IMBHs could wander off-center in galaxies. We use the cosmological simulation ASTRID to study the orbital and radiative properties of wandering IMBHs in massive galaxies at $z\sim3$. We find that this populat… ▽ More Intermediate-Mass Black Holes (IMBHs) of $10^3-10^6 \, M_\odot$ are commonly found at the center of dwarf galaxies. Simulations and observations convincingly show that a sizable population of IMBHs could wander off-center in galaxies. We use the cosmological simulation ASTRID to study the orbital and radiative properties of wandering IMBHs in massive galaxies at $z\sim3$. We find that this population of black holes has large orbital inclinations ($60^\circ\pm22^\circ$) with respect to the principal plane of the host. The eccentricity of their orbits is also significant ($0.6\pm0.2$) and decreases with time. Wandering IMBHs undergo spikes of accretion activity around the pericenter of their orbits, with rates $10^{-3}-10^{-5}$ times the Eddington rate and a median accretion duty cycle of $\sim 12\%$. Their typical spectral energy distribution peaks in the infrared at $\sim 11 \, μ\rm m$ rest-frame. Assuming a standard value of $10\%$ for the matter-to-energy radiative efficiency, IMBHs reach $2-10$ keV X-ray luminosities $>10^{37} \, \mathrm{erg\,s^{-1}}$ for $\sim10\%$ of the time. This luminosity corresponds to fluxes $>10^{-15} \, \mathrm{erg \, s^{-1} \, cm^{-2}}$ within $10$ Mpc. They could be challenging to detect because of competing emissions from X-ray binaries and the interstellar medium. X-ray luminosities $> 10^{41} \, \mathrm{erg \, s^{-1}}$, in the hyper-luminous X-ray sources (HLXs) regime, are reached by $\sim 7\%$ of the IMBHs. These findings suggest that HLXs are a small subset of the wandering IMBH population, which is characterized by luminosities $10^3-10^4$ times fainter. Dedicated surveys are needed to assess the demographics of this missing population of black holes. △ Less

Submitted 2 April, 2023; v1 submitted 28 October, 2022; originally announced October 2022.

Comments: Accepted for publication in MNRAS. This is the final version of the manuscript. 9 pages, 7 figures

A vast population of wandering and merging IMBHs at cosmic noon

Authors: Tiziana Di Matteo, Yueying Ni, Nianyi Chen, Rupert Croft, Simeon Bird, Fabio Pacucci, Angelo Ricarte, Michael Tremmel

Abstract: Massive black holes in the centers of galaxies today must have grown by several orders of magnitude from seed black holes formed at early times. Detecting a population of intermediate mass black holes (IMBHs) can provide constraints on these elusive BH seeds. Here we use the large volume, cosmological hydrodynamical simulation Astrid, which includes IMBH seeds and dynamical friction to investigate… ▽ More Massive black holes in the centers of galaxies today must have grown by several orders of magnitude from seed black holes formed at early times. Detecting a population of intermediate mass black holes (IMBHs) can provide constraints on these elusive BH seeds. Here we use the large volume, cosmological hydrodynamical simulation Astrid, which includes IMBH seeds and dynamical friction to investigate the population of IMBH seeds. Dynamical friction is largely inefficient at sinking and merging seed IMBHs at high-z. This leads to an extensive population (several hundred per galaxy) of wandering IMBHs in large halos at z~2. A small fraction of these IMBHs are detectable as HLXs, Hyper Luminous X-ray sources. Importantly, at z ~ 2, IMBHs mergers produce the peak of GW events. We find close to a million GW events in Astrid between z=2-3 involving seed IMBH mergers. These GW events (almost all detectable by LISA) at cosmic noon should provide strong constraints on IMBH seed models and their formation mechanisms. At the center of massive galaxies, where the number of IMBHs can be as high as 10-100, SMBH-IMBH pairs can form. These Intermediate mass ratio inspirals (IMRIs) and extreme mass ratio inspirals (EMRIs), will require the next generation of milli-muHz space-based GW interferometers to be detected. Large populations of IMBHs around massive black holes will probe their environments and MBH causal structure. △ Less

Submitted 26 October, 2022; originally announced October 2022.

Super-resolution simulation of the Fuzzy Dark Matter cosmological model

Authors: Meris Sipp, Patrick LaChance, Rupert Croft, Yueying Ni, Tiziana Di Matteo

Abstract: AI super-resolution, combining deep learning and N-body simulations has been shown to successfully reproduce the large scale structure and halo abundances in the Lambda Cold Dark Matter cosmological model. Here, we extend its use to models with a different dark matter content, in this case Fuzzy Dark Matter (FDM), in the approximation that the difference is encoded in the initial power spectrum. W… ▽ More AI super-resolution, combining deep learning and N-body simulations has been shown to successfully reproduce the large scale structure and halo abundances in the Lambda Cold Dark Matter cosmological model. Here, we extend its use to models with a different dark matter content, in this case Fuzzy Dark Matter (FDM), in the approximation that the difference is encoded in the initial power spectrum. We focus on redshift z = 2, with simulations that model smaller scales and lower masses, the latter by two orders of magnitude, than has been done in previous AI super-resolution work. We find that the super-resolution technique can reproduce the power spectrum and halo mass function to within a few percent of full high resolution calculations. We also find that halo artifacts, caused by spurious numerical fragmentation of filaments, are equally present in the super-resolution outputs. Although we have not trained the super-resolution algorithm using full quantum pressure FDM simulations, the fact that it performs well at the relevant length and mass scales means that it has promise as technique which could avoid the very high computational cost of the latter, in some contexts. We conclude that AI super-resolution can become a useful tool to extend the range of dark matter models covered in mock catalogs. △ Less

Submitted 23 October, 2022; originally announced October 2022.

Comments: 7 pages, 4 figures

Robust field-level inference with dark matter halos

Authors: Helen Shao, Francisco Villaescusa-Navarro, Pablo Villanueva-Domingo, Romain Teyssier, Lehman H. Garrison, Marco Gatti, Derek Inman, Yueying Ni, Ulrich P. Steinwandel, Mihir Kulkarni, Eli Visbal, Greg L. Bryan, Daniel Angles-Alcazar, Tiago Castro, Elena Hernandez-Martinez, Klaus Dolag

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

Submitted 14 September, 2022; originally announced September 2022.

Comments: 25 pages, 11 figures, summary video: https://youtu.be/qkw92Z6owJU

Prominence fine structures in weakly twisted and highly twisted magnetic flux ropes

Authors: J. H. Guo, Y. W. Ni, Y. H. Zhou, Y. Guo, B. Schmieder, P. F. Chen

Abstract: Many prominences are supported by magnetic flux ropes. One important question is how we can determine whether the flux rope is weakly-twisted or strongly-twisted. In this paper, we attempted to check whether prominences supported by weakly-twisted and strongly-twisted flux ropes can manifest different features so that we might distinguish the two types of magnetic structures by their appearance. W… ▽ More Many prominences are supported by magnetic flux ropes. One important question is how we can determine whether the flux rope is weakly-twisted or strongly-twisted. In this paper, we attempted to check whether prominences supported by weakly-twisted and strongly-twisted flux ropes can manifest different features so that we might distinguish the two types of magnetic structures by their appearance. We performed pseudo three-dimensional simulations of two magnetic flux ropes with different twists. We found that the resulting two prominences differ in many aspects. The prominence supported by a weakly-twisted flux rope is composed mainly of transient threads, forming high-speed flows inside the prominence. Its horns are evident. Conversely, the one supported by a highly-twisted flux rope consists mainly of stable quasi-stationary threads, including longer independently trapped threads and shorter magnetically connected threads. It is also revealed that the prominence spine deviates from the flux rope axis in the vertical direction and from the photospheric polarity inversion line projected on the solar surface, especially for the weakly-twisted magnetic flux rope. The two types of prominences differ significantly in appearance. It is also suggested that a piling-up of short threads in highly-twisted flux ropes might account for the vertical-like threads in some prominences. △ Less

Submitted 14 September, 2022; originally announced September 2022.

Comments: 10 pages, 7 figures, accepted for publication in A&A

Journal ref: A&A 667, A89 (2022)

Ultramassive black holes formed by triple quasar mergers at $z\sim 2$

Authors: Yueying Ni, Tiziana DiMatteo, Nianyi Chen, Rupert Croft, Simeon Bird

Abstract: The origin of rare and elusive ultramassive black holes (UMBH, with MBH > 1e10 Msun) is an open question. Using the large volume cosmological hydrodynamic simulation ASTRID, we report on the formation of an extremely massive UMBH with MBH ~ 1e11 Msun at z~2. The UMBH is assembled as a result of two successive mergers of massive galaxies each with stellar mass M* > 3e11 Msun that also produces a br… ▽ More The origin of rare and elusive ultramassive black holes (UMBH, with MBH > 1e10 Msun) is an open question. Using the large volume cosmological hydrodynamic simulation ASTRID, we report on the formation of an extremely massive UMBH with MBH ~ 1e11 Msun at z~2. The UMBH is assembled as a result of two successive mergers of massive galaxies each with stellar mass M* > 3e11 Msun that also produces a bright, rare triple quasar system powered by three ~10^9 Msun black holes. The second merger of supermassive black holes (SMBHs) follows the first after 150 Myrs. The merger events lead to sustained Eddington accretion onto the central SMBH, forming an UMBH in the center of a massive compact stellar core with M* > 2e12 Msun. The strong feedback of the UMBH quenches the surrounding star formation to < 10 Msun/yr in the inner 50 kpc/h region. There are two more UMBHs with MBH > 5e10 Msun at z>2 in ASTRID which are also produced by major mergers of galaxies, and their progenitors can be observed as quasar triplets of lower luminosity. The rarely observed quasar multiples can be the cradle of UMBHs at high redshift, and likely end up in the center of the most massive clusters. △ Less

Submitted 6 September, 2022; v1 submitted 2 September, 2022; originally announced September 2022.

Comments: 7 pages, 4 figures, comments welcome

First detection of transverse vertical oscillation during the expansion of coronal loops

Authors: Qingmin Zhang, Chuan Li, Dong Li, Ye Qiu, Yanjie Zhang, Yiwei Ni

Abstract: In this Letter, we perform a detailed analysis of the M5.5-class eruptive flare occurring in active region 12929 on 2022 January 20. The eruption of a hot channel generates a fast coronal mass ejection (CME) and a dome-shaped extreme-ultraviolet (EUV) wave at speeds of 740$-$860 km s$^{-1}$. The CME is associated with a type II radio burst, implying that the EUV wave is a fast-mode shock wave. Dur… ▽ More In this Letter, we perform a detailed analysis of the M5.5-class eruptive flare occurring in active region 12929 on 2022 January 20. The eruption of a hot channel generates a fast coronal mass ejection (CME) and a dome-shaped extreme-ultraviolet (EUV) wave at speeds of 740$-$860 km s$^{-1}$. The CME is associated with a type II radio burst, implying that the EUV wave is a fast-mode shock wave. During the impulsive phase, the flare shows quasi-periodic pulsations (QPPs) in EUV, hard X-ray, and radio wavelengths. The periods of QPPs range from 18 s to 113 s, indicating that flare energy is released and nonthermal electrons are accelerated intermittently with multiple time scales. The interaction between the EUV wave and low-lying adjacent coronal loops (ACLs) results in contraction, expansion, and transverse vertical oscillation of ACLs. The speed of contraction in 171, 193, and 211 Å is higher than that in 304 Å. The periods of oscillation are 253 s and 275 s in 304 Å and 171 Å, respectively. A new scenario is proposed to explain the interaction. The equation that interprets the contraction and oscillation of the overlying coronal loops above a flare core can also interpret the expansion and oscillation of ACLs, suggesting that the two phenomena are the same in essence. △ Less

Submitted 31 August, 2022; originally announced September 2022.

Comments: 11 pages, 12 figures, accepted for publication in ApJL

Statistics of Galactic-Scale Quasar Pairs at Cosmic Noon

Authors: Yue Shen, Hsiang-Chih Hwang, Masamune Oguri, Nianyi Chen, Tiziana Di Matteo, Yueying Ni, Simeon Bird, Nadia Zakamska, Xin Liu, Yu-Ching Chen, Kaitlin M. Kratter

Abstract: The statistics of galactic-scale quasar pairs can elucidate our understanding of the dynamical evolution of supermassive black hole (SMBH) pairs, the duty cycles of quasar activity in mergers, or even the nature of dark matter, but have been challenging to measure at cosmic noon, the prime epoch of massive galaxy and SMBH formation. Here we measure a double quasar fraction of… ▽ More The statistics of galactic-scale quasar pairs can elucidate our understanding of the dynamical evolution of supermassive black hole (SMBH) pairs, the duty cycles of quasar activity in mergers, or even the nature of dark matter, but have been challenging to measure at cosmic noon, the prime epoch of massive galaxy and SMBH formation. Here we measure a double quasar fraction of $\sim 6.2\pm0.5\times 10^{-4}$ integrated over $\sim 0.3-3$ arcsec separations (projected physical separations of $\sim 3-30\,{\rm kpc}$ at $z\sim 2$) in luminous ($L_{\rm bol}>10^{45.8}\,{\rm erg\,s^{-1}}$) unobscured quasars at $1.5<z<3.5$, using Gaia EDR3-resolved pairs around SDSS DR16 quasars. The measurement was based on a sample of 60 Gaia-resolved double quasars (out of 487 Gaia pairs dominated by quasar+star superpositions) at these separations, corrected for pair completeness in Gaia, which we quantify as functions of pair separation, magnitude of the primary, and magnitude contrast. The double quasar fraction increases towards smaller separations by a factor of $\sim 5$ over these scales. The division between physical quasar pairs and lensed quasars in our sample is currently unknown, requiring dedicated follow-up observations (in particular, deep, sub-arcsec-resolution IR imaging for the closest pairs). Intriguingly, at this point the observed pair statistics are in rough agreement with theoretical predictions both for the lensed quasar population in mock catalogs and for dual quasars in cosmological hydrodynamic simulations. Upcoming wide-field imaging/spectroscopic space missions such as Euclid, CSST and Roman, combined with targeted follow-up observations, will conclusively measure the abundances and host galaxy properties of galactic-scale quasar pairs, offset AGNs, and sub-arcsec lensed quasars across cosmic time. △ Less

Submitted 9 August, 2022; originally announced August 2022.

Comments: 19 pages, 9 figures; submitted to ApJ

Journal ref: ApJ, 943, 38 (2023)

Properties and Evolution of Dual and Offset AGN in the ASTRID Simulation at $z \sim 2$

Authors: Nianyi Chen, Tiziana Di Matteo, Yueying Ni, Michael Tremmel, Colin DeGraf, Yue Shen, A. Miguel Holgado, Simeon Bird, Rupert Croft, Yu Feng

Abstract: We examine the dual (both BHs active) and offset (one BH active) AGN population (comprising $\sim$ 2000 pairs at $0.5\,\text{kpc}\lesssimΔr<30\,\text{kpc}$) at $z=2\sim3$ in the ASTRID simulation covering (360 cMpc)${^3}$. The dual (offset) AGN make up $3.0(2.2)\%$ of all AGN at $z=2$. The dual fraction is roughly constant while the offset fraction increases by a factor of ten from $z=4\sim2$. Com… ▽ More We examine the dual (both BHs active) and offset (one BH active) AGN population (comprising $\sim$ 2000 pairs at $0.5\,\text{kpc}\lesssimΔr<30\,\text{kpc}$) at $z=2\sim3$ in the ASTRID simulation covering (360 cMpc)${^3}$. The dual (offset) AGN make up $3.0(2.2)\%$ of all AGN at $z=2$. The dual fraction is roughly constant while the offset fraction increases by a factor of ten from $z=4\sim2$. Compared with the full AGN population, duals are characterized by a low $M_\text{BH}/M_*$ ratio, a high specific star-formation rate (sSFR) of $\sim 1\,\text{Gyr}^{-1}$, and a high Eddington ratio ($\sim 0.05$, double that of single AGN). The dual AGN are formed in major galaxy mergers (typically involving $M_\text{halo}<10^{13}\,M_\odot$), with BHs that have similar masses. At small separations (when their host galaxies are in the late phase of the merger) duals become $2\sim8$ times brighter (albeit more obscured) than at larger separations. $80\%$ of these bright, close duals merge in the simulation within $\sim500\,\text{Myrs}$. Notably, the initially less-massive BH in duals frequently becomes the brighter AGN during the galaxy merger. In offset AGN, the active BH is typically $\gtrsim 10$ times more massive than its non-active counterpart and than most BHs in duals. Offsets are predominantly formed in minor galaxy mergers with the active BH residing in the center of massive halos ($ M_\text{ halo}\sim 10^{13-14}\,M_\odot$). In these deep potentials, gas stripping is common and the secondary quickly deactivates. The stripping also leads to inefficient orbital decay amongst offsets, which stall at $Δr\sim5\,\text{kpc}$ for a few hundred Myrs. △ Less

Submitted 9 August, 2022; originally announced August 2022.

Comments: 17 Pages, 18 Figures; comments welcome!

The origin and evolution of the normal Type Ia SN 2018aoz with infant-phase reddening and excess emission

Authors: Yuan Qi Ni, Dae-Sik Moon, Maria R. Drout, Abigail Polin, David J. Sand, Santiago GonzÁlez-GaitÁn, Sang Chul Kim, Youngdae Lee, Hong Soo Park, D. Andrew Howell, Peter E. Nugent, Anthony L. Piro, Peter J. Brown, LluÍs Galbany, Jamison Burke, Daichi Hiramatsu, Griffin Hosseinzadeh, Stefano Valenti, Niloufar Afsariardchi, Jennifer E. Andrews, John Antoniadis, Rachael L. Beaton, K. Azalee Bostroem, Raymond G. Carlberg, S. Bradley Cenko , et al. (18 additional authors not shown)

Abstract: SN~2018aoz is a Type Ia SN with a $B$-band plateau and excess emission in the infant-phase light curves $\lesssim$ 1 day after first light, evidencing an over-density of surface iron-peak elements as shown in our previous study. Here, we advance the constraints on the nature and origin of SN~2018aoz based on its evolution until the nebular phase. Near-peak spectroscopic features show the SN is int… ▽ More SN~2018aoz is a Type Ia SN with a $B$-band plateau and excess emission in the infant-phase light curves $\lesssim$ 1 day after first light, evidencing an over-density of surface iron-peak elements as shown in our previous study. Here, we advance the constraints on the nature and origin of SN~2018aoz based on its evolution until the nebular phase. Near-peak spectroscopic features show the SN is intermediate between two subtypes of normal Type Ia: Core-Normal and Broad-Line. The excess emission could have contributions from the radioactive decay of surface iron-peak elements as well as ejecta interaction with either the binary companion or a small torus of circumstellar material. Nebular-phase limits on H$α$ and He~I favour a white dwarf companion, consistent with the small companion size constrained by the low early SN luminosity, while the absence of [O~I] and He~I disfavours a violent merger of the progenitor. Of the two main explosion mechanisms proposed to explain the distribution of surface iron-peak elements in SN~2018aoz, the asymmetric Chandrasekhar-mass explosion is less consistent with the progenitor constraints and the observed blueshifts of nebular-phase [Fe~II] and [Ni~II]. The helium-shell double-detonation explosion is compatible with the observed lack of C spectral features, but current 1-D models are incompatible with the infant-phase excess emission, $B_{\rm max}-V_{\rm max}$ color, and absence of nebular-phase [Ca~II]. Although the explosion processes of SN~2018aoz still need to be more precisely understood, the same processes could produce a significant fraction of Type Ia SNe that appear normal after $\sim$ 1 day. △ Less

Submitted 24 June, 2022; originally announced June 2022.

Comments: Submitted for publication in ApJ. 35 pages, 16 figures, 7 tables

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

Authors: Madeline A. Marshall, Katelyn Watts, Stephen Wilkins, Tiziana Di Matteo, Jussi K. Kuusisto, William J. Roper, Aswin P. Vijayan, Yueying Ni, Yu Feng, Rupert A. C. Croft

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

Submitted 25 August, 2022; v1 submitted 17 June, 2022; originally announced June 2022.

Comments: Published in MNRAS. The BlueTides Mock Image Catalogue is publicly available as a High Level Science Product (HLSP) via the Mikulski Archive for Space Telescopes (MAST), at https://doi.org/10.17909/er09-4527. Example codes for using the catalogue data are available at https://github.com/madelinemarshall/BlueTidesMockImageCatalogue

Calibration procedures for the CHASE/HIS science data

Authors: Ye Qiu, ShiHao Rao, Chuan Li, Cheng Fang, MingDe Ding, Zhen Li, YiWei Ni, WenBo Wang, Jie Hong, Qi Hao, Yu Dai, PengFei Chen, XiaoSheng Wan, Zhi Xu, Wei You, Yuan Yuan, HongJiang Tao, XianSheng Li, YuKun He, Qiang Liu

Abstract: The Hα line is an important optical line in solar observations containing the information from the photosphere to the chromosphere. To study the mechanisms of solar eruptions and the plasma dynamics in the lower atmosphere, the Chinese Hα Solar Explorer (CHASE) was launched into a Sun-synchronous orbit on October 14, 2021. The scientific payload of the CHASE satellite is the Hα Imaging Spectrograp… ▽ More The Hα line is an important optical line in solar observations containing the information from the photosphere to the chromosphere. To study the mechanisms of solar eruptions and the plasma dynamics in the lower atmosphere, the Chinese Hα Solar Explorer (CHASE) was launched into a Sun-synchronous orbit on October 14, 2021. The scientific payload of the CHASE satellite is the Hα Imaging Spectrograph (HIS). The CHASE/HIS acquires, for the first time, seeing-free Hα spectroscopic observations with high spectral and temporal resolutions. It consists of two observational modes. The raster scanning mode provides full-Sun or region-of-interest spectra at Hα (6559.7-6565.9 Å) and Fe I (6567.8-6570.6 Å) wavebands. The continuum imaging mode obtains full-Sun photospheric images at around 6689 Å. In this paper, we present detailed calibration procedures for the CHASE/HIS science data, including the dark-field and flat-field correction, slit image curvature correction, wavelength and intensity calibration, and coordinate transformation. The higher-level data products can be directly used for scientific research. △ Less

Submitted 12 May, 2022; originally announced May 2022.

Comments: 9 pages, 7 figures

Journal ref: Sci. China-Phys. Mech. Astron. 65, 289603 (2022)

The Chinese Hα Solar Explorer (CHASE) mission: An overview

Authors: Chuan Li, Cheng Fang, Zhen Li, MingDe Ding, PengFei Chen, Ye Qiu, Wei You, Yuan Yuan, MinJie An, HongJiang Tao, XianSheng Li, Zhe Chen, Qiang Liu, Gui Mei, Liang Yang, Wei Zhang, WeiQiang Cheng, JianXin Chen, ChangYa Chen, Qiang Gu, QingLong Huang, MingXing Liu, ChengShan Han, HongWei Xin, ChangZheng Chen , et al. (10 additional authors not shown)

Abstract: The Chinese Hα Solar Explorer (CHASE), dubbed "Xihe" - Goddess of the Sun, was launched on October 14, 2021 as the first solar space mission of China National Space Administration (CNSA). The CHASE mission is designed to test a newly developed satellite platform and to acquire the spectroscopic observations in the Hα waveband. The Hα Imaging Spectrograph (HIS) is the scientific payload of the CHAS… ▽ More The Chinese Hα Solar Explorer (CHASE), dubbed "Xihe" - Goddess of the Sun, was launched on October 14, 2021 as the first solar space mission of China National Space Administration (CNSA). The CHASE mission is designed to test a newly developed satellite platform and to acquire the spectroscopic observations in the Hα waveband. The Hα Imaging Spectrograph (HIS) is the scientific payload of the CHASE satellite. It consists of two observational modes: raster scanning mode and continuum imaging mode. The raster scanning mode obtains full-Sun or region-of-interest spectral images from 6559.7 to 6565.9 Å and from 6567.8 to 6570.6 Å with 0.024 Å pixel spectral resolution and 1 minute temporal resolution. The continuum imaging mode obtains photospheric images in continuum around 6689 Å with the full width at half maximum of 13.4 Å. The CHASE mission will advance our understanding of the dynamics of solar activity in the photosphere and chromosphere. In this paper, we present an overview of the CHASE mission including the scientific objectives, HIS instrument overview, data calibration flow, and first results of on-orbit observations. △ Less

Submitted 12 May, 2022; originally announced May 2022.

Comments: 9 pages, 6 figures

Journal ref: Sci. China-Phys. Mech. Astron. 65, 289602 (2022)

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

Authors: Aklant Kumar Bhowmick, Laura Blecha, Yueying Ni, Tiziana Di Matteo, Paul Torrey, Luke Zoltan Kelley, Mark Vogelsberger, Rainer Weinberger, Lars Hernquist

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

Submitted 11 May, 2022; originally announced May 2022.

Comments: 21 pages, 15 figures