Showing 1–50 of 141 results for author: Gao, L

Baryon Acoustic Oscillations analyses with Density-Split Statistics

Authors: Tengpeng Xu, Yan-Chuan Cai, Yun Chen, Mark Neyrinck, Liang Gao, Qiao Wang

Abstract: Accurate modeling for the evolution of the Baryon Acoustic Oscillations (BAO) is essential for using it as a standard ruler to probe cosmology. We explore the non-linearity of the BAO in different environments using the density-split statistics and compare them to the case of the conventional two-point correlation function (2PCF). We detect density-dependent shifts for the position of the BAO with… ▽ More Accurate modeling for the evolution of the Baryon Acoustic Oscillations (BAO) is essential for using it as a standard ruler to probe cosmology. We explore the non-linearity of the BAO in different environments using the density-split statistics and compare them to the case of the conventional two-point correlation function (2PCF). We detect density-dependent shifts for the position of the BAO with respect to its linear version using halos from N-body simulations. Around low/high-densities, the scale of the BAO expands/contracts due to non-linear peculiar velocities. As the simulation evolves from redshift 1 to 0, the difference in the magnitude of the shifts between high- and low-density regions increases from the sub-percent to the percent level. In contrast, the scale of the BAO does not evolve in the total 2PCF in the same redshift range. The width of the BAO around high density regions increases as the universe evolves, similar to the known broadening of the BAO in the 2PCF due to non-linear evolution. In contrast, the width is smaller and stable for low density regions. We discuss possible implications for the reconstructions of the BAO in light of our results. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: 16 pages, 10 figures

Hyper-sampling imaging

Authors: Ze Zhang, Hemeng Xue, Mingtao Shang, Hongfei Yu, Jinchao Liang, Meiling Guan, Chengming Sun, Huahua Wang, Shufeng Wang, Zhengyu Ye, Feng Gao, Lu Gao

Abstract: In our research, we have developed a novel mechanism that allows for a significant reduction in the smallest sampling unit of digital image sensors (DIS) to as small as 1/16th of a pixel, through measuring the intra-pixel quantum efficiency for the first time and recomputing the image. Employing our method, the physical sampling resolution of DIS can be enhanced by 16 times. The method has undergo… ▽ More In our research, we have developed a novel mechanism that allows for a significant reduction in the smallest sampling unit of digital image sensors (DIS) to as small as 1/16th of a pixel, through measuring the intra-pixel quantum efficiency for the first time and recomputing the image. Employing our method, the physical sampling resolution of DIS can be enhanced by 16 times. The method has undergone rigorous testing in real-world imaging scenarios. △ Less

Submitted 27 June, 2024; originally announced June 2024.

Constraints on Ultra Heavy Dark Matter Properties from Dwarf Spheroidal Galaxies with LHAASO Observations

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (255 additional authors not shown)

Abstract: In this work we try to search for signals generated by ultra-heavy dark matter at the Large High Altitude Air Shower Observatory (LHAASO) data. We look for possible gamma-ray by dark matter annihilation or decay from 16 dwarf spheroidal galaxies in the field of view of LHAASO. Dwarf spheroidal galaxies are among the most promising targets for indirect detection of dark matter which have low fluxes… ▽ More In this work we try to search for signals generated by ultra-heavy dark matter at the Large High Altitude Air Shower Observatory (LHAASO) data. We look for possible gamma-ray by dark matter annihilation or decay from 16 dwarf spheroidal galaxies in the field of view of LHAASO. Dwarf spheroidal galaxies are among the most promising targets for indirect detection of dark matter which have low fluxes of astrophysical $γ$-ray background while large amount of dark matter. By analyzing more than 700 days observational data at LHAASO, no significant dark matter signal from 1 TeV to 1 EeV is detected. Accordingly we derive the most stringent constraints on the ultra-heavy dark matter annihilation cross-section up to EeV. The constraints on the lifetime of dark matter in decay mode are also derived. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: 17 pages, 12 figures, accepted by PRL

Discovery of Very-high-energy Gamma-ray Emissions from the Low Luminosity AGN NGC 4278 by LHAASO

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (255 additional authors not shown)

Abstract: The first source catalog of Large High Altitude Air Shower Observatory reported the detection of a very-high-energy gamma ray source, 1LHAASO J1219+2915. In this paper a further detailed study of the spectral and temporal behavior of this point-like source have been carried. The best-fit position of the TeV source ($\rm{RA}=185.05^{\circ}\pm0.04^{\circ}$, $\rm{Dec}=29.25^{\circ}\pm0.03^{\circ}$) i… ▽ More The first source catalog of Large High Altitude Air Shower Observatory reported the detection of a very-high-energy gamma ray source, 1LHAASO J1219+2915. In this paper a further detailed study of the spectral and temporal behavior of this point-like source have been carried. The best-fit position of the TeV source ($\rm{RA}=185.05^{\circ}\pm0.04^{\circ}$, $\rm{Dec}=29.25^{\circ}\pm0.03^{\circ}$) is compatible with NGC 4278 within $\sim0.03$ degree. Variation analysis shows an indication of the variability at a few months level in the TeV band, which is consistent with low frequency observations. Based on these observations, we report the detection of TeV $γ$-ray emissions from this low-luminosity AGN NGC 4278. The observations by LHAASO-WCDA during active period has a significance level of 8.8\,$σ$ with best-fit photon spectral index $\varGamma=2.56\pm0.14$ and a flux $f_{1-10\,\rm{TeV}}=(7.0\pm1.1_{\rm{sta}}\pm0.35_{\rm{syst}})\times10^{-13}\,\rm{photons\,cm^{-2}\,s^{-1}}$, or approximately $5\%$ of the Crab Nebula. The discovery of VHE from NGC 4278 indicates that the compact, weak radio jet can efficiently accelerate particles and emit TeV photons. △ Less

Submitted 13 May, 2024; originally announced May 2024.

Comments: 11 pages, 5 figures

Kinematic morphology of low-mass galaxies in IllustrisTNG

Authors: Guangquan Zeng, Lan Wang, Liang Gao, Hang Yang

Abstract: The origin of diverse kinematic morphologies observed in low-mass galaxies is unclear. In this study, we investigate the kinematic morphologies of central galaxies with stellar mass $10^{8.5-9.0} M_\odot$ at $z=0$ in the TNG50-1 cosmological simulation. The majority of the low-mass galaxies in TNG50-1 are dispersion-dominated, consistent with observations. By tracing the evolutionary histories of… ▽ More The origin of diverse kinematic morphologies observed in low-mass galaxies is unclear. In this study, we investigate the kinematic morphologies of central galaxies with stellar mass $10^{8.5-9.0} M_\odot$ at $z=0$ in the TNG50-1 cosmological simulation. The majority of the low-mass galaxies in TNG50-1 are dispersion-dominated, consistent with observations. By tracing the evolutionary histories of simulated low-mass galaxies, we find that while most stars form in rotating cold gas discs, the orientation of the star-forming discs relative to the galaxies may evolve with cosmic time. If the cold gas disc remains aligning with the galaxy during its evolution, stars formed at different times share the same rotational direction, leading to a rotation-dominated system. On the contrary, frequent misalignment of cold gas disc would result in a dispersion-dominated system. In addition, we also find that the two-body scattering can have a non-negligible numerical heating effect on the simulated galaxy morphology, especially at central regions of galaxies and for relatively low-mass galaxies. By comparing results of simulations with different resolutions, our results suggest that the simulated morphology of galaxies is roughly reliable when their number of stellar particles exceeds about $10^{4}$, and bulge morphology of galaxies can not be resolved robustly at the resolution level of TNG50-1. △ Less

Submitted 3 July, 2024; v1 submitted 22 April, 2024; originally announced April 2024.

Comments: 12 pages, 10 figures (additional 1 figure in the appendix). Accepted for publication in MNRAS

LHAASO-KM2A detector simulation using Geant4

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (254 additional authors not shown)

Abstract: KM2A is one of the main sub-arrays of LHAASO, working on gamma ray astronomy and cosmic ray physics at energies above 10 TeV. Detector simulation is the important foundation for estimating detector performance and data analysis. It is a big challenge to simulate the KM2A detector in the framework of Geant4 due to the need to track numerous photons from a large number of detector units (>6000) with… ▽ More KM2A is one of the main sub-arrays of LHAASO, working on gamma ray astronomy and cosmic ray physics at energies above 10 TeV. Detector simulation is the important foundation for estimating detector performance and data analysis. It is a big challenge to simulate the KM2A detector in the framework of Geant4 due to the need to track numerous photons from a large number of detector units (>6000) with large altitude difference (30 m) and huge coverage (1.3 km^2). In this paper, the design of the KM2A simulation code G4KM2A based on Geant4 is introduced. The process of G4KM2A is optimized mainly in memory consumption to avoid memory overffow. Some simpliffcations are used to signiffcantly speed up the execution of G4KM2A. The running time is reduced by at least 30 times compared to full detector simulation. The particle distributions and the core/angle resolution comparison between simulation and experimental data of the full KM2A array are also presented, which show good agreement. △ Less

Submitted 7 April, 2024; originally announced April 2024.

Simulating emission line galaxies for the next generation of large-scale structure surveys

Authors: Wenxiang Pei, Qi Guo, Ming Li, Qiao Wang, Jiaxin Han, Jia Hu, Tong Su, Liang Gao, Jie Wang, Yu Luo, Chengliang Wei

Abstract: We investigate emission line galaxies across cosmic time by combining the modified L-Galaxies semi-analytical galaxy formation model with the JiuTian cosmological simulation. We improve the tidal disruption model of satellite galaxies in L-Galaxies to address the time dependence problem. We utilise the public code CLOUDY to compute emission line ratios for a grid of HII region models. The emission… ▽ More We investigate emission line galaxies across cosmic time by combining the modified L-Galaxies semi-analytical galaxy formation model with the JiuTian cosmological simulation. We improve the tidal disruption model of satellite galaxies in L-Galaxies to address the time dependence problem. We utilise the public code CLOUDY to compute emission line ratios for a grid of HII region models. The emission line models assume the same initial mass function as that used to generate the spectral energy distribution of semi-analytical galaxies, ensuring a coherent treatment for modelling the full galaxy spectrum. By incorporating these emission line ratios with galaxy properties, we reproduce observed luminosity functions for H$α$, H$β$, [OII], and [OIII] in the local Universe and at high redshifts. We also find good agreement between model predictions and observations for auto-correlation and cross-correlation functions of [OII]-selected galaxies, as well as their luminosity dependence. The bias of emission line galaxies depends on both luminosity and redshift. At lower redshifts, it remains constant with increasing luminosity up to around $\sim 10^{42.5}\rm \, erg\,s^{-1}$ and then rises steeply for higher luminosities. The transition luminosity increases with redshift and becomes insignificant above $z$=1.5. Generally, galaxy bias shows an increasing trend with redshift. However, for luminous galaxies, the bias is higher at low redshifts, as the strong luminosity dependence observed at low redshifts diminishes at higher redshifts. We provide a fitting formula for the bias of emission line galaxies as a function of luminosity and redshift, which can be utilised for large-scale structure studies with future galaxy surveys. △ Less

Submitted 29 March, 2024; originally announced April 2024.

Comments: 22 pages, 18 figures, 5 tables, including 3 Appendix; Accepted by MNRAS

The Influence of Baryons on Low-mass Haloes

Authors: Haonan Zheng, Sownak Bose, Carlos S. Frenk, Liang Gao, Adrian Jenkins, Shihong Liao, Volker Springel, Jie Wang, Simon D. M. White

Abstract: The Voids-within-Voids-within-Voids (VVV) project used dark-matter-only simulations to study the abundance and structure of dark matter haloes over the full mass range populated in the standard $Λ\mathrm{CDM}$ cosmology. Here we explore how baryonic effects modify these results for $z=0$ halo masses in the range $10^4$ to $10^7~\mathrm{M_\odot}$, below the threshold for galaxy formation. Our main… ▽ More The Voids-within-Voids-within-Voids (VVV) project used dark-matter-only simulations to study the abundance and structure of dark matter haloes over the full mass range populated in the standard $Λ\mathrm{CDM}$ cosmology. Here we explore how baryonic effects modify these results for $z=0$ halo masses in the range $10^4$ to $10^7~\mathrm{M_\odot}$, below the threshold for galaxy formation. Our main study focuses on three simulations from identical initial conditions at $z=127$, one following dark matter only, one including non-radiative gas, and one additionally including the baryonic physics relevant in this halo mass range (cooling and photoheating). In the non-radiative simulation, above $10^{5.5}~\mathrm{M_\odot}$, halo abundance and internal structure are very similar to the dark-matter-only simulation, and the baryon to dark matter ratio is everywhere close to the cosmic value. At lower mass, this ratio drops and haloes are less concentrated and less massive in the non-radiative case. Test simulations at higher resolution show this to be mainly a resolution effect; the expected drop in baryon content due to residual pressure effects only becomes substantial for $z=0$ haloes below $\sim 10^{2.7}~\mathrm{M_\odot}$. However, gas is heated by reionization at $z=6$ in our ``full physics'' run, and this results in almost complete expulsion of gas from all haloes in our simulated mass range. This suppresses the halo mass function by $\sim 30 \%$, lowers halo concentration, and consequently weakens the dark matter annihilation signal by $\sim 40-60 \%$. △ Less

Submitted 25 March, 2024; originally announced March 2024.

Comments: 12+2 pages, 9 figures

Measurements of All-Particle Energy Spectrum and Mean Logarithmic Mass of Cosmic Rays from 0.3 to 30 PeV with LHAASO-KM2A

Authors: The LHAASO Collaboration, Zhen Cao, F. Aharonian, Q. An, A. Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen , et al. (256 additional authors not shown)

Abstract: We present the measurements of all-particle energy spectrum and mean logarithmic mass of cosmic rays in the energy range of 0.3-30 PeV using data collected from LHAASO-KM2A between September 2021 and December 2022, which is based on a nearly composition-independent energy reconstruction method, achieving unprecedented accuracy. Our analysis reveals the position of the knee at… ▽ More We present the measurements of all-particle energy spectrum and mean logarithmic mass of cosmic rays in the energy range of 0.3-30 PeV using data collected from LHAASO-KM2A between September 2021 and December 2022, which is based on a nearly composition-independent energy reconstruction method, achieving unprecedented accuracy. Our analysis reveals the position of the knee at $3.67 \pm 0.05 \pm 0.15$ PeV. Below the knee, the spectral index is found to be -$2.7413 \pm 0.0004 \pm 0.0050$, while above the knee, it is -$3.128 \pm 0.005 \pm 0.027$, with the sharpness of the transition measured with a statistical error of 2%. The mean logarithmic mass of cosmic rays is almost heavier than helium in the whole measured energy range. It decreases from 1.7 at 0.3 PeV to 1.3 at 3 PeV, representing a 24% decline following a power law with an index of -$0.1200 \pm 0.0003 \pm 0.0341$. This is equivalent to an increase in abundance of light components. Above the knee, the mean logarithmic mass exhibits a power law trend towards heavier components, which is reversal to the behavior observed in the all-particle energy spectrum. Additionally, the knee position and the change in power-law index are approximately the same. These findings suggest that the knee observed in the all-particle spectrum corresponds to the knee of the light component, rather than the medium-heavy components. △ Less

Submitted 26 March, 2024; v1 submitted 15 March, 2024; originally announced March 2024.

Comments: 8 pages, 3 figures

Journal ref: Physical Review Letters 132, 131002 (2024)

HiFAST: an HI data calibration and imaging pipeline for FAST

Authors: Yingjie Jing, Jie Wang, Chen Xu, Ziming Liu, Qingze Chen, Tiantian Liang, Jinlong Xu, Yixian Cao, Jing Wang, Huijie Hu, Chuan-Peng Zhang, Qi Guo, Liang Gao, Mei Ai, Hengqian Gan, Xuyang Gao, Jinlin Han, Ligang Hou, Zhipeng Hou, Peng Jiang, Xu Kong, Fujia Li, Zerui Liu, Li Shao, Hengxing Pan , et al. (8 additional authors not shown)

Abstract: The Five-hundred-meter Aperture Spherical radio Telescope (FAST) has the largest aperture and a 19-beam L-band receiver, making it powerful for investigating the neutral hydrogen atomic gas (HI) in the universe. We present HiFAST (https://hifast.readthedocs.io), a dedicated, modular, and self-contained calibration and imaging pipeline for processing the HI data of FAST. The pipeline consists of fr… ▽ More The Five-hundred-meter Aperture Spherical radio Telescope (FAST) has the largest aperture and a 19-beam L-band receiver, making it powerful for investigating the neutral hydrogen atomic gas (HI) in the universe. We present HiFAST (https://hifast.readthedocs.io), a dedicated, modular, and self-contained calibration and imaging pipeline for processing the HI data of FAST. The pipeline consists of frequency-dependent noise diode calibration, baseline fitting, standing wave removal using an FFT-based method, flux density calibration, stray radiation correction, and gridding to produce data cubes. These modules can be combined as needed to process the data from most FAST observation modes: tracking, drift scanning, On-The-Fly mapping, and most of their variants. With HiFAST, the RMS noises of the calibrated spectra from all 19 beams were only slightly (~ 5%) higher than the theoretical expectation. The results for the extended source M33 and the point sources are consistent with the results from Arecibo. The moment maps (0,1 and 2) of M33 agree well with the results from the Arecibo Galaxy Environment Survey (AGES) with a fractional difference of less than 10%. For a common sample of 221 sources with signal-to-noise ratio S/N >10 from the Arecibo Legacy Fast ALFA (ALFALFA) survey, the mean value of fractional difference in the integrated flux density, $S_{\mathrm{int}}$, between the two datasets is approximately 0.005 %, with a dispersion of 15.4%. Further checks on the integrated flux density of 23 sources with seven observations indicate that the variance in the flux density of the source with luminous objects ($S_\mathrm{int}$ $ > 2.5$ Jy km s$^{-1}$) is less than 5%. Our tests suggest that the FAST telescope, with the efficient, precise, and user-friendly pipeline HiFAST, will yield numerous significant scientific findings in the investigation of the HI in the universe. △ Less

Submitted 30 January, 2024; originally announced January 2024.

Comments: Accepted by SCPMA. 21 pages, 14 figures. The pipeline is accessible at https://hifast.readthedocs.io

A younger Universe implied by satellite pair correlations from SDSS observations of massive galaxy groups

Authors: Qing Gu, Qi Guo, Marius Cautun, Shi Shao, Wenxiang Pei, Wenting Wang, Liang Gao, Jie Wang

Abstract: Many of the satellites of galactic-mass systems such as the Miky Way, Andromeda and Centaurus A show evidence of coherent motions to a larger extent than most of the systems predicted by the standard cosmological model. It is an open question if correlations in satellite orbits are present in systems of different masses. Here , we report an analysis of the kinematics of satellite galaxies around m… ▽ More Many of the satellites of galactic-mass systems such as the Miky Way, Andromeda and Centaurus A show evidence of coherent motions to a larger extent than most of the systems predicted by the standard cosmological model. It is an open question if correlations in satellite orbits are present in systems of different masses. Here , we report an analysis of the kinematics of satellite galaxies around massive galaxy groups. Unlike what is seen in Milky Way analogues, we find an excess of diametrically opposed pairs of satellites that have line-of-sight velocity offsets from the central galaxy of the same sign. This corresponds to a $\pmb{6.0σ}$ ($\pmb{p}$-value $\pmb{=\ 9.9\times10^{-10}}$) detection of non-random satellite motions. Such excess is predicted by up-to-date cosmological simulations but the magnitude of the effect is considerably lower than in observations. The observational data is discrepant at the $\pmb{4.1σ}$ and $\pmb{3.6σ}$ level with the expectations of the Millennium and the Illustris TNG300 cosmological simulations, potentially indicating that massive galaxy groups assembled later in the real Universe. The detection of velocity correlations of satellite galaxies and tension with theoretical predictions is robust against changes in sample selection. Using the largest sample to date, our findings demonstrate that the motions of satellite galaxies represent a challenge to the current cosmological model. △ Less

Submitted 18 January, 2024; originally announced January 2024.

Comments: 28 pages, 9 figures, accepted for publication in Nature Astronomy

Constraints on Axion-like Particles from the Observation of Galactic Sources by LHAASO

Authors: Jun Li, Xiao-Jun Bi, Lin-Qing Gao, Xiaoyuan Huang, Run-Min Yao, Peng-Fei Yin

Abstract: High-energy photons may oscillate with axion-like particles (ALPs) when they propagate through the Milky Way's magnetic field, resulting in an alteration in the observed photon energy spectrum. The ultra-high energy gamma-ray spectra, measured by the Large High Altitude Air Shower Observatory (LHAASO) up to $\mathcal{O}(1)~\mathrm{PeV}$, provide a promising opportunity to investigate the ALP-photo… ▽ More High-energy photons may oscillate with axion-like particles (ALPs) when they propagate through the Milky Way's magnetic field, resulting in an alteration in the observed photon energy spectrum. The ultra-high energy gamma-ray spectra, measured by the Large High Altitude Air Shower Observatory (LHAASO) up to $\mathcal{O}(1)~\mathrm{PeV}$, provide a promising opportunity to investigate the ALP-photon oscillation effect. In this study, we utilize the gamma-ray spectra of four Galactic sources measured by LHAASO, including the Crab Nebula, LHAASO J2226+6057, LHAASO J1908+0621, and LHAASO J1825-1326, to explore this effect. We employ the $\rm CL_s$ method to set constraints on the ALP parameters. Combing the observations of the four sources, our analysis reveals that the ALP-photon coupling $g_{aγ}$ is constrained to be smaller than $1.4\times10^{-10}$ ${\rm GeV}^{-1}$ for the ALP mass of $\sim 4\times10^{-7} ~\mathrm{eV}$ at the 95\% C.L. By combing the observations of the Crab Nebula from LHAASO and other experiments, we find that the ALP-photon coupling could be set to be about $7.2\times10^{-11}$ ${\rm GeV}^{-1}$ for the ALP mass $\sim 4 \times10^{-7}~\mathrm{eV}$ , which is in close proximity to the CAST constraint. △ Less

Submitted 3 January, 2024; originally announced January 2024.

Does or did the supernova remnant Cassiopeia A operate as a PeVatron?

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (255 additional authors not shown)

Abstract: For decades, supernova remnants (SNRs) have been considered the prime sources of Galactic Cosmic rays (CRs). But whether SNRs can accelerate CR protons to PeV energies and thus dominate CR flux up to the knee is currently under intensive theoretical and phenomenological debate. The direct test of the ability of SNRs to operate as CR PeVatrons can be provided by ultrahigh-energy (UHE;… ▽ More For decades, supernova remnants (SNRs) have been considered the prime sources of Galactic Cosmic rays (CRs). But whether SNRs can accelerate CR protons to PeV energies and thus dominate CR flux up to the knee is currently under intensive theoretical and phenomenological debate. The direct test of the ability of SNRs to operate as CR PeVatrons can be provided by ultrahigh-energy (UHE; $E_γ\geq 100$~TeV) $γ$-rays. In this context, the historical SNR Cassiopeia A (Cas A) is considered one of the most promising target for UHE observations. This paper presents the observation of Cas A and its vicinity by the LHAASO KM2A detector. The exceptional sensitivity of LHAASO KM2A in the UHE band, combined with the young age of Cas A, enabled us to derive stringent model-independent limits on the energy budget of UHE protons and nuclei accelerated by Cas A at any epoch after the explosion. The results challenge the prevailing paradigm that Cas A-type SNRs are major suppliers of PeV CRs in the Milky Way. △ Less

Submitted 25 October, 2023; originally announced October 2023.

Comments: 11 pages, 3 figures, Accepted by the APJL

The abundance of dark matter haloes down to Earth mass

Authors: Haonan Zheng, Sownak Bose, Carlos S. Frenk, Liang Gao, Adrian Jenkins, Shihong Liao, Yizhou Liu, Jie Wang

Abstract: We use the Voids-within-Voids-within-Voids (VVV) simulations, a suite of successive nested N-body simulations with extremely high resolution (denoted, from low to high resolution, by L0 to L7), to test the Press-Schechter (PS), Sheth-Tormen (ST), and extended Press-Schechter (EPS) formulae for the halo abundance over the entire mass range, from mini-haloes of $10^{-6}\ \mathrm{M_\odot}$, to cluste… ▽ More We use the Voids-within-Voids-within-Voids (VVV) simulations, a suite of successive nested N-body simulations with extremely high resolution (denoted, from low to high resolution, by L0 to L7), to test the Press-Schechter (PS), Sheth-Tormen (ST), and extended Press-Schechter (EPS) formulae for the halo abundance over the entire mass range, from mini-haloes of $10^{-6}\ \mathrm{M_\odot}$, to cluster haloes of $10^{15}\ \mathrm{M_\odot}$, at different redshifts, from $z=30$ to the present. We find that at $z=0$ and $z=2$, ST best reproduces the results of L0, which has the mean cosmic density (overdensity $δ=0$), at $10^{11-15} ~\mathrm{M_\odot}$. The higher resolution levels (L1-L7) are biased underdense regions ($δ<-0.6$). The EPS formalism takes this into account since it gives the mass function of a region conditioned, in this case, on having a given underdensity. EPS provides good matches to these higher levels, with deviations $\lesssim 20\%$, at $10^{-6-12.5} ~\mathrm{M_\odot}$. At $z \sim 7-15$, the ST predictions for L0 and the EPS for L1-L7 show somewhat larger deviations from the simulation results. However, at even higher redshifts, $z \sim 30$, EPS fits the simulations well again. We confirm our results by picking more subvolumes from the L0 simulation, finding that our conclusions depend only weakly on the size and overdensity of the region. The good agreement of EPS with the higher-level simulations implies that PS (or ST) gives an accurate description of the total halo mass function in representative regions of the universe. △ Less

Submitted 10 April, 2024; v1 submitted 24 October, 2023; originally announced October 2023.

Comments: 10 pages, 5 figures (additional 2 figures in the appendix)

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 528, Issue 4, March 2024, Pages 7300-7309

MaNGA DynPop -- VI. Matter density slopes from dynamical models of 6000 galaxies versus cosmological simulations: the interplay between baryonic and dark matter

Authors: Shubo Li, Ran Li, Kai Zhu, Shengdong Lu, Michele Cappellari, Shude Mao, Chunxiang Wang, Liang Gao

Abstract: We try to understand the trends in the mass density slopes as a function of galaxy properties. We use the results from the best Jeans Anisotropic Modelling (JAM) of the integral-field stellar kinematics for near 6000 galaxies from the MaNGA DynPop project, with stellar masses of $10^{9-12}\ {\rm M_{\odot}}$, including both early-type and late-type galaxies. We use the mass-weighted density slopes… ▽ More We try to understand the trends in the mass density slopes as a function of galaxy properties. We use the results from the best Jeans Anisotropic Modelling (JAM) of the integral-field stellar kinematics for near 6000 galaxies from the MaNGA DynPop project, with stellar masses of $10^{9-12}\ {\rm M_{\odot}}$, including both early-type and late-type galaxies. We use the mass-weighted density slopes for the stellar $\overlineγ_*$, dark $\overlineγ_{\rm DM}$, and total $\overlineγ_{\rm T}$ mass from the MaNGA DynPop project. The $\overlineγ_{\rm T}$ approaches a constant value of 2.2 for high $σ_{\rm e}$ galaxies, and flattens for lg$(σ_{\rm e}/{\rm km\ s^{-1}})\lesssim2.3$, reaching 1.5 for lg$(σ_{\rm e}/{\rm km\ s^{-1}})\approx1.8$. The total and stellar slopes track each other tightly, with $\overlineγ_{\rm T}\approx\overlineγ_*-0.174$ over the full $σ_{\rm e}$ range. This confirms the dominance of stellar matter within $R_{\rm e}$. We also show that there is no perfect conspiracy between baryonic and dark matter, as $\overlineγ_*$ and $\overlineγ_{\rm DM}$ do not vary inversely within the $σ_{\rm e}$ range. We find that the central galaxies from TNG50 and TNG100 simulations do not reproduce the observed galaxy mass distribution, which we attribute to the overestimated dark matter fraction, possibly due to a constant IMF and excessive adiabatic contraction effects in the simulations. Finally, we present the stacked dark matter density profiles and show that they are slightly steeper than the pure dark matter simulation prediction of $\overlineγ_{\rm DM}\approx1$, suggesting moderate adiabatic contraction in the central region of galaxies. Our work demonstrate the power of stellar dynamics modelling for probing the interaction between stellar and dark matter and testing galaxy formation theories. △ Less

Submitted 1 April, 2024; v1 submitted 20 October, 2023; originally announced October 2023.

Comments: 17 pages, 14 figures; accepted by MNRAS

Constraints on Axion-like Particles from the Observation of GRB 221009A by LHAASO

Authors: Lin-Qing Gao, Xiao-Jun Bi, Jun Li, Run-Min Yao, Peng-Fei Yin

Abstract: The LHAASO collaboration recently reported the measurement of the gamma-ray spectra of GRB 221009A, which is the brightest burst ever, covering an energy range from 0.3 $\mathrm{TeV}$ to about 10 $\mathrm{TeV}$. Based on the observation, we investigate the ALP-photon oscillation effect in the host galaxy of GRB 221009A and the Milky Way. The ${\rm CL_s}$ method is applied to set constraints on the… ▽ More The LHAASO collaboration recently reported the measurement of the gamma-ray spectra of GRB 221009A, which is the brightest burst ever, covering an energy range from 0.3 $\mathrm{TeV}$ to about 10 $\mathrm{TeV}$. Based on the observation, we investigate the ALP-photon oscillation effect in the host galaxy of GRB 221009A and the Milky Way. The ${\rm CL_s}$ method is applied to set constraints on the ALP parameters in this study. Given the uncertain magnetic field configuration in the host galaxy, we use three different models: a homogeneous magnetic field model, a magnetic field model identical to that of the Milky Way, and a model constructed from the HST observations of the host galaxy. We find that the constraints derived using these three host galaxy magnetic field models are comparable. Our results are complementary in the small ALP mass regions compared with other experiments. △ Less

Submitted 15 January, 2024; v1 submitted 17 October, 2023; originally announced October 2023.

Comments: 8 pages, 13 figures

Very high energy gamma-ray emission beyond 10 TeV from GRB 221009A

Authors: Zhen Cao, F. Aharonian, Q. An, A. Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (255 additional authors not shown)

Abstract: The highest energy gamma-rays from gamma-ray bursts (GRBs) have important implications for their radiation mechanism. Here we report for the first time the detection of gamma-rays up to 13 TeV from the brightest GRB 221009A by the Large High Altitude Air-shower Observatory (LHAASO). The LHAASO-KM2A detector registered more than 140 gamma-rays with energies above 3 TeV during 230$-$900s after the t… ▽ More The highest energy gamma-rays from gamma-ray bursts (GRBs) have important implications for their radiation mechanism. Here we report for the first time the detection of gamma-rays up to 13 TeV from the brightest GRB 221009A by the Large High Altitude Air-shower Observatory (LHAASO). The LHAASO-KM2A detector registered more than 140 gamma-rays with energies above 3 TeV during 230$-$900s after the trigger. The intrinsic energy spectrum of gamma-rays can be described by a power-law after correcting for extragalactic background light (EBL) absorption. Such a hard spectrum challenges the synchrotron self-Compton (SSC) scenario of relativistic electrons for the afterglow emission above several TeV. Observations of gamma-rays up to 13 TeV from a source with a measured redshift of z=0.151 hints more transparency in intergalactic space than previously expected. Alternatively, one may invoke new physics such as Lorentz Invariance Violation (LIV) or an axion origin of very high energy (VHE) signals. △ Less

Submitted 22 November, 2023; v1 submitted 13 October, 2023; originally announced October 2023.

Comments: 49pages, 11figures

Journal ref: Science Advances, 9, eadj2778 (2023) 15 November 2023

Constraints on Axion-like Particles from Observations of Mrk 421 using the ${\rm CL_s}$ Method

Authors: Lin-Qing Gao, Xiao-Jun Bi, Jun-Guang Guo, Wenbin Lin, Peng-Fei Yin

Abstract: Axion-like particles (ALPs) may undergo mixing with photons in the presence of astrophysical magnetic fields, leading to alterations in the observed high energy $γ$-ray spectra. In this study, we investigate the ALP-photon oscillation effect using the spectra of the blazar Mrk 421 over 15 observation periods measured by Major Atmospheric Gamma Imaging Cherenkov Telescopes (MAGIC) and Fermi Large A… ▽ More Axion-like particles (ALPs) may undergo mixing with photons in the presence of astrophysical magnetic fields, leading to alterations in the observed high energy $γ$-ray spectra. In this study, we investigate the ALP-photon oscillation effect using the spectra of the blazar Mrk 421 over 15 observation periods measured by Major Atmospheric Gamma Imaging Cherenkov Telescopes (MAGIC) and Fermi Large Area Telescope (Fermi-LAT). Compared with previous studies, we generate mock data under the ALP hypothesis and employ the ${\rm CL_s}$ method to set constraints on the ALP parameters. This method is widely utilized in high energy experiments and avoids the exclusion of specific parameter regions where distinguishing between the null and ALP hypotheses is challenging. We find that the ALP-photon coupling $g_{aγ}$ is constrained to be smaller than $\sim 2\times10^{-11}$ GeV$^{-1}$ for ALP masses ranging from $10^{-9}$ eV to $10^{-7}$ eV at the 95\% confidence level. We also present the constraints derived from the TS distribution under the null hypothesis, which is commonly utilized in previous astrophysical ALP studies. Our results reveal that the combined constraints of all the periods obtained from both methods are consistent. However, the ${\rm CL_s}$ method remains effective in cases where the latter method fails to provide constraints for specific observation periods. △ Less

Submitted 5 March, 2024; v1 submitted 5 September, 2023; originally announced September 2023.

Comments: 10 pages, 27 figures

The First LHAASO Catalog of Gamma-Ray Sources

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (255 additional authors not shown)

Abstract: We present the first catalog of very-high energy and ultra-high energy gamma-ray sources detected by the Large High Altitude Air Shower Observatory (LHAASO). The catalog was compiled using 508 days of data collected by the Water Cherenkov Detector Array (WCDA) from March 2021 to September 2022 and 933 days of data recorded by the Kilometer Squared Array (KM2A) from January 2020 to September 2022.… ▽ More We present the first catalog of very-high energy and ultra-high energy gamma-ray sources detected by the Large High Altitude Air Shower Observatory (LHAASO). The catalog was compiled using 508 days of data collected by the Water Cherenkov Detector Array (WCDA) from March 2021 to September 2022 and 933 days of data recorded by the Kilometer Squared Array (KM2A) from January 2020 to September 2022. This catalog represents the main result from the most sensitive large coverage gamma-ray survey of the sky above 1 TeV, covering declination from $-$20$^{\circ}$ to 80$^{\circ}$. In total, the catalog contains 90 sources with an extended size smaller than $2^\circ$ and a significance of detection at $> 5σ$. Based on our source association criteria, 32 new TeV sources are proposed in this study. Among the 90 sources, 43 sources are detected with ultra-high energy ($E > 100$ TeV) emission at $> 4σ$ significance level. We provide the position, extension, and spectral characteristics of all the sources in this catalog. △ Less

Submitted 27 November, 2023; v1 submitted 26 May, 2023; originally announced May 2023.

Comments: 40 pages, 13 figures, 4 tables

Journal ref: The Astrophysical Journal Supplement Series, 271 (2024) 25

The Lobster Eye Imager for Astronomy Onboard the SATech-01 Satellite

Authors: Z. X. Ling, X. J. Sun, C. Zhang, S. L. Sun, G. Jin, S. N. Zhang, X. F. Zhang, J. B. Chang, F. S. Chen, Y. F. Chen, Z. W. Cheng, W. Fu, Y. X. Han, H. Li, J. F. Li, Y. Li, Z. D. Li, P. R. Liu, Y. H. Lv, X. H. Ma, Y. J. Tang, C. B. Wang, R. J. Xie, Y. L. Xue, A. L. Yan , et al. (101 additional authors not shown)

Abstract: The Lobster Eye Imager for Astronomy (LEIA), a pathfinder of the Wide-field X-ray Telescope of the Einstein Probe (EP) mission, was successfully launched onboard the SATech-01 satellite of the Chinese Academy of Sciences on 27 July 2022. In this paper, we introduce the design and on-ground test results of the LEIA instrument. Using state-of-the-art Micro-Pore Optics (MPO), a wide field-of-view (Fo… ▽ More The Lobster Eye Imager for Astronomy (LEIA), a pathfinder of the Wide-field X-ray Telescope of the Einstein Probe (EP) mission, was successfully launched onboard the SATech-01 satellite of the Chinese Academy of Sciences on 27 July 2022. In this paper, we introduce the design and on-ground test results of the LEIA instrument. Using state-of-the-art Micro-Pore Optics (MPO), a wide field-of-view (FoV) of 346 square degrees (18.6 degrees * 18.6 degrees) of the X-ray imager is realized. An optical assembly composed of 36 MPO chips is used to focus incident X-ray photons, and four large-format complementary metal-oxide semiconductor (CMOS) sensors, each of 6 cm * 6 cm, are used as the focal plane detectors. The instrument has an angular resolution of 4 - 8 arcmin (in FWHM) for the central focal spot of the point spread function, and an effective area of 2 - 3 cm2 at 1 keV in essentially all the directions within the field of view. The detection passband is 0.5 - 4 keV in the soft X-rays and the sensitivity is 2 - 3 * 10-11 erg s-1 cm-2 (about 1 mini-Crab) at 1,000 second observation. The total weight of LEIA is 56 kg and the power is 85 W. The satellite, with a design lifetime of 2 years, operates in a Sun-synchronous orbit of 500 km with an orbital period of 95 minutes. LEIA is paving the way for future missions by verifying in flight the technologies of both novel focusing imaging optics and CMOS sensors for X-ray observation, and by optimizing the working setups of the instrumental parameters. In addition, LEIA is able to carry out scientific observations to find new transients and to monitor known sources in the soft X-ray band, albeit limited useful observing time available. △ Less

Submitted 24 May, 2023; originally announced May 2023.

Comments: Accepted by RAA

Assessing Mass Loss and Stellar-to-Halo Mass Ratio of Satellite Galaxies: A Galaxy-Galaxy Lensing Approach Utilizing DECaLS DR8 Data

Authors: Chunxiang Wang, Ran Li, Huanyuan Shan, Weiwei Xu, Ji Yao, Yingjie Jing, Liang Gao, Nan Li, Yushan Xie, Kai Zhu, Hang Yang, Qingze Chen

Abstract: The galaxy-galaxy lensing technique allows us to measure the subhalo mass of satellite galaxies, studying their mass loss and evolution within galaxy clusters and providing direct observational validation for theories of galaxy formation. In this study, we use the weak gravitational lensing observations from DECaLS DR8, in combination with the redMaPPer galaxy cluster catalog from Sloan Digital Sk… ▽ More The galaxy-galaxy lensing technique allows us to measure the subhalo mass of satellite galaxies, studying their mass loss and evolution within galaxy clusters and providing direct observational validation for theories of galaxy formation. In this study, we use the weak gravitational lensing observations from DECaLS DR8, in combination with the redMaPPer galaxy cluster catalog from Sloan Digital Sky Survey data (SDSS) DR8 to accurately measure the dark matter halo mass of satellite galaxies. We confirm a significant increase in the stellar-to-halo mass ratio of satellite galaxies with their halo-centric radius, indicating clear evidence of mass loss due to tidal stripping. Additionally, we find that this mass loss is strongly dependent on the mass of the satellite galaxies, with satellite galaxies above $10^{11}~{\rm M_{\odot}/h}$ experiencing more pronounced mass loss compared to lower mass satellites, reaching 86\% at projected halo-centric radius $0.5R_{\rm 200c}$. The average mass loss rate, when not considering halo-centric radius, displays a U-shaped variation with stellar mass, with galaxies of approximately $4\times10^{10}~{\rm M_{\odot}/h}$ exhibiting the least mass loss, around 60\%. We compare our results with state-of-the-art hydrodynamical numerical simulations and find that the satellite galaxy stellar-to-halo mass ratio in the outskirts of galaxy clusters is higher compared to the predictions of the Illustris-TNG project about factor 5. Furthermore, the Illustris-TNG project's numerical simulations did not predict the observed dependence of satellite galaxy mass loss rate on satellite galaxy mass. △ Less

Submitted 31 January, 2024; v1 submitted 23 May, 2023; originally announced May 2023.

Comments: 14 pages, 9 figures

Measurement of ultra-high-energy diffuse gamma-ray emission of the Galactic plane from 10 TeV to 1 PeV with LHAASO-KM2A

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (255 additional authors not shown)

Abstract: The diffuse Galactic $γ$-ray emission, mainly produced via interactions between cosmic rays and the interstellar medium and/or radiation field, is a very important probe of the distribution, propagation, and interaction of cosmic rays in the Milky Way. In this work we report the measurements of diffuse $γ$-rays from the Galactic plane between 10 TeV and 1 PeV energies, with the square kilometer ar… ▽ More The diffuse Galactic $γ$-ray emission, mainly produced via interactions between cosmic rays and the interstellar medium and/or radiation field, is a very important probe of the distribution, propagation, and interaction of cosmic rays in the Milky Way. In this work we report the measurements of diffuse $γ$-rays from the Galactic plane between 10 TeV and 1 PeV energies, with the square kilometer array of the Large High Altitude Air Shower Observatory (LHAASO). Diffuse emissions from the inner ($15^{\circ}<l<125^{\circ}$, $|b|<5^{\circ}$) and outer ($125^{\circ}<l<235^{\circ}$, $|b|<5^{\circ}$) Galactic plane are detected with $29.1σ$ and $12.7σ$ significance, respectively. The outer Galactic plane diffuse emission is detected for the first time in the very- to ultra-high-energy domain ($E>10$~TeV). The energy spectrum in the inner Galaxy regions can be described by a power-law function with an index of $-2.99\pm0.04$, which is different from the curved spectrum as expected from hadronic interactions between locally measured cosmic rays and the line-of-sight integrated gas content. Furthermore, the measured flux is higher by a factor of $\sim3$ than the prediction. A similar spectrum with an index of $-2.99\pm0.07$ is found in the outer Galaxy region, and the absolute flux for $10\lesssim E\lesssim60$ TeV is again higher than the prediction for hadronic cosmic ray interactions. The latitude distributions of the diffuse emission are consistent with the gas distribution, while the longitude distributions show clear deviation from the gas distribution. The LHAASO measurements imply that either additional emission sources exist or cosmic ray intensities have spatial variations. △ Less

Submitted 19 August, 2023; v1 submitted 9 May, 2023; originally announced May 2023.

Comments: 12 pages, 8 figures, 5 tables; accepted for publication in Physical Review Letters; source mask file provided as ancillary file

Journal ref: Phys. Rev. Lett. 131, 151001 (2023)

MaNGA DynPop -- IV. Stacked total density profile of galaxy groups and clusters from combining dynamical models of integral-field stellar kinematics and galaxy-galaxy lensing

Authors: Chunxiang Wang, Ran Li, Kai Zhu, Huanyuan Shan, Weiwei Xu, Michele Cappellari, Liang Gao, Nan Li, Shengdong Lu, Shude Mao, Ji Yao, Yushan Xie

Abstract: We present the measurement of total and stellar/dark matter decomposed mass density profile around a sample of galaxy groups and clusters with dynamical masses derived from integral-field stellar kinematics from the MaNGA survey in Paper I and weak lensing derived from the DECaLS imaging survey. Combining the two data sets enables accurate measurement of the radial density distribution from severa… ▽ More We present the measurement of total and stellar/dark matter decomposed mass density profile around a sample of galaxy groups and clusters with dynamical masses derived from integral-field stellar kinematics from the MaNGA survey in Paper I and weak lensing derived from the DECaLS imaging survey. Combining the two data sets enables accurate measurement of the radial density distribution from several kpc to Mpc scales. Intriguingly, we find that the excess surface density derived from stellar kinematics in the inner region cannot be explained by simply adding an NFW dark matter halo extrapolated from lensing measurement at a larger scale to a stellar mass component derived from the NASA-Sloan Atlas (NSA) catalogue. We find that a good fit to both data sets requires a stellar mass normalization about 3 times higher than that derived from the NSA catalogue, which would require an unrealistically too-heavy initial mass function for stellar mass estimation. If we keep the stellar mass normalization to that of the NSA catalogue but allow a varying inner dark matter density profile, we obtain an asymptotic slope of $γ_{\rm gnfw}$= $1.82_{-0.25}^{+0.15}$ and $γ_{\rm gnfw}$= $1.48_{-0.41}^{+0.20}$ for the group bin and the cluster bin respectively, significantly steeper than the NFW case. We also compare the total mass inner density slopes with those from TNG300 and find that the values from the simulation are lower than the observation by about $2σ$ level. △ Less

Submitted 18 December, 2023; v1 submitted 23 April, 2023; originally announced April 2023.

Comments: 12 pages, 9 figures, accepted by MNRAS

MaNGA DynPop -- III. Stellar dynamics versus stellar population relations in 6000 early-type and spiral galaxies: Fundamental Plane, mass-to-light ratios, total density slopes, and dark matter fractions

Authors: Kai Zhu, Shengdong Lu, Michele Cappellari, Ran Li, Shude Mao, Liang Gao, Junqiang Ge

Abstract: We present dynamical scaling relations, combined with the stellar population properties, for a subsample of about 6000 nearby galaxies with the most reliable dynamical models extracted from the full MaNGA sample of 10K galaxies. We show that the inclination-corrected mass plane (MP) for both early-type galaxies (ETGs) and late-type galaxies (LTGs), which links dynamical mass, projected half-light… ▽ More We present dynamical scaling relations, combined with the stellar population properties, for a subsample of about 6000 nearby galaxies with the most reliable dynamical models extracted from the full MaNGA sample of 10K galaxies. We show that the inclination-corrected mass plane (MP) for both early-type galaxies (ETGs) and late-type galaxies (LTGs), which links dynamical mass, projected half-light radius $R_{\rm e}$, and the second stellar velocity moment $σ_{\rm e}$ within $R_{\rm e}$, satisfies the virial theorem and is even tighter than the uncorrected one. We find a clear parabolic relation between $\lg(M/L)_{\rm e}$, the total mass-to-light ratio within a sphere of radius $R_{\rm e}$, and $\lgσ_{\rm e}$, with the $M/L$ increasing with $σ_{\rm e}$ and for older stellar populations. However, the relation for ETGs is linear and the one for the youngest galaxies is constant. We confirm and improve the relation between average logarithmic total density slopes $\overline{γ_{\rm T}}$ and $σ_{\rm e}$: $\overline{γ_{\rm T}}$ become steeper with increasing $σ_{\rm e}$ until $\lg(σ_{\rm e}/{\rm km\,s^{-1}})\approx 2.2$ and then remain constant around $\overline{γ_{\rm T}}\approx -2.2$. The $\overline{γ_{\rm T}}-σ_{\rm e}$ variation is larger for LTGs than ETGs. At fixed $σ_{\rm e}$ the total density profiles steepen with galaxy age and for ETGs. We find generally low dark matter fractions, median $f_{\rm DM}(<R_{\rm e})=8$ per cent, within a sphere of radius $R_{\rm e}$. However, we find that $f_{\rm DM}(<R_{\rm e})$ depends on $σ_{\rm e}$ better than stellar mass: dark matter increases to a median $f_{\rm DM}(<R_{\rm e})=33$ percent for galaxies with $σ_{\rm e}\lesssim100{\rm km\,s^{-1}}$. The increased $f_{\rm DM}(<R_{\rm e})$ at low $σ_{\rm e}$ explains the parabolic $\lg(M/L)_{\rm e}-\lgσ_{\rm e}$ relation. △ Less

Submitted 4 November, 2023; v1 submitted 23 April, 2023; originally announced April 2023.

Comments: 26 pages, 21 figures, 1 table; Accepted by MNRAS

MaNGA DynPop -- I. Quality-assessed stellar dynamical modelling from integral-field spectroscopy of 10K nearby galaxies: a catalogue of masses, mass-to-light ratios, density profiles and dark matter

Authors: Kai Zhu, Shengdong Lu, Michele Cappellari, Ran Li, Shude Mao, Liang Gao

Abstract: This is the first paper in our series on the combined analysis of the Dynamics and stellar Population (DynPop) for the MaNGA survey in the final SDSS Data Release 17 (DR17). Here we present a catalogue of dynamically-determined quantities for over 10000 nearby galaxies based on integral-field stellar kinematics from the MaNGA survey. The dynamical properties are extracted using the axisymmetric Je… ▽ More This is the first paper in our series on the combined analysis of the Dynamics and stellar Population (DynPop) for the MaNGA survey in the final SDSS Data Release 17 (DR17). Here we present a catalogue of dynamically-determined quantities for over 10000 nearby galaxies based on integral-field stellar kinematics from the MaNGA survey. The dynamical properties are extracted using the axisymmetric Jeans Anisotropic Modelling (JAM) method, which was previously shown to be the most accurate for this kind of study. We assess systematic uncertainties using eight dynamical models with different assumptions. We use two orientations of the velocity ellipsoid: either cylindrically-aligned JAM$_{\rm cyl}$ or spherically-aligned JAM$_{\rm sph}$. We also make four assumptions for the models' dark vs. luminous matter distributions: (1) mass-follows-light, (2) free NFW dark halo, (3) cosmologically-constrained NFW halo, (4) generalized NFW dark halo, i.e. with free inner slope. In this catalogue, we provide the quantities related to the mass distributions (e.g. the density slopes and enclosed mass within a sphere of a given radius for total mass, stellar mass, and dark matter mass components). We also provide the complete models which can be used to compute the full luminous and mass distribution of each galaxy. Additionally, we visually assess the qualities of the models to help with model selections. We estimate the observed scatter in the measured quantities which decreases as expected with improvements in quality. For the best data quality, we find a remarkable consistency of measured quantities between different models, highlighting the robustness of the results. △ Less

Submitted 4 June, 2023; v1 submitted 23 April, 2023; originally announced April 2023.

Comments: 30 pages, 15 figures, 4 tables, accepted by MNRAS, the catalogue is available here https://manga-dynpop.github.io/

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 522, Issue 4, pp.6326-6353, (2023)

Exploring the multiband gravitational wave background with a semi-analytic galaxy formation model

Authors: Zhencheng Li, Zhen Jiang, Xi-Long Fan, Yun Chen, Liang Gao, Qi Guo, Shenghua Yu

Abstract: An enormous number of compact binary systems, spanning from stellar to supermassive levels, emit substantial gravitational waves during their final evolutionary stages, thereby creating a stochastic gravitational wave background (SGWB). We calculate the merger rates of stellar compact binaries and massive black hole binaries using a semi-analytic galaxy formation model -- Galaxy Assembly with Bina… ▽ More An enormous number of compact binary systems, spanning from stellar to supermassive levels, emit substantial gravitational waves during their final evolutionary stages, thereby creating a stochastic gravitational wave background (SGWB). We calculate the merger rates of stellar compact binaries and massive black hole binaries using a semi-analytic galaxy formation model -- Galaxy Assembly with Binary Evolution (GABE) in a unified and self-consistent approach, followed by an estimation of the multi-band SGWB contributed by those systems. We find that the amplitudes of the principal peaks of the SGWB energy density are within one order of magnitude $Ω_{GW} \sim 10^{-9}- 10^{-8}$. This SGWB could easily be detected by the Square Kilometre Array (SKA), as well as planned interferometric detectors, such as the Einstein Telescope (ET) and the Laser Interferometer Space Antenna (LISA). The energy density of this background varies as $Ω_{GW} \propto f^{2/3}$ in SKA band. The shape of the SGWB spectrum in the frequency range $\sim[10^{-4}$,$1]$Hz could allow the LISA to distinguish the black hole seed models. The amplitude of the SGWB from merging stellar binary black holes (BBHs) at $\sim 100$ Hz is approximately 10 and 100 times greater than those from merging binary neutron stars (BNSs) and neutron-star-black-hole (NSBH) mergers, respectively. Note that, since the cosmic star formation rate density predicted by GABE is somewhat lower than observational results by $\sim 0.2$ dex at z < $\sim 2$, the amplitude of the SGWB in the frequency range $\sim[1$, $10^{4}]$ Hz may be underestimated by a similar factor at most. △ Less

Submitted 4 December, 2023; v1 submitted 17 April, 2023; originally announced April 2023.

Comments: 12 pages, 5 figures

Journal ref: MNRAS, 527, 3, 2024, 5616-5626

The mass accretion history of dark matter haloes down to Earth mass

Authors: Yizhou Liu, Liang Gao, Sownak Bose, Carlos S. Frenk, Adrian Jenkins, Volker Springel, Jie Wang, Simon D. M. White, Haonan Zheng

Abstract: We take advantage of the unprecedented dynamical range provided by the "Cosmic-Zoom" project to study the mass accretion history (MAH) of present-day dark matter haloes over the entire mass range present in the $Λ$CDM paradigm when the dark matter is made of weakly interacting massive particles of mass $100\ \mathrm{GeV}$. In particular, we complement previous studies by exploring the MAHs of halo… ▽ More We take advantage of the unprecedented dynamical range provided by the "Cosmic-Zoom" project to study the mass accretion history (MAH) of present-day dark matter haloes over the entire mass range present in the $Λ$CDM paradigm when the dark matter is made of weakly interacting massive particles of mass $100\ \mathrm{GeV}$. In particular, we complement previous studies by exploring the MAHs of haloes with mass from $10^8\ h^{-1}\mathrm{M_{\odot}}$ down to Earth mass, $10^{-6}\ h^{-1}\mathrm{M_{\odot}}$. The formation redshift of low-mass haloes anti-correlates weakly with mass, peaking at $z=3$ for haloes of mass $10^{-4}\ h^{-1}\mathrm{M_{\odot}}$. Even lower masses are affected by the free-streaming cutoff in the primordial spectrum of density fluctuations and form at lower redshift. We compare MAHs in our simulations with predictions from two analytical models based on the extended Press-Schechter theory (EPS), and three empirical models derived by fitting and extrapolating either results from cosmological $N$-body simulations or Monte Carlo realizations of halo growth. All models fit our simulations reasonably well over the mass range for which they were calibrated. While the empirical models match better for more massive haloes, $M>10^{10}\ h^{-1}\mathrm{M_{\odot}}$, the analytical models do better when extrapolated down to Earth mass. At the higher masses, we explore the correlation between local environment density and MAH, finding that biases are relatively weak, with typical MAHs for haloes in extremely low-density and in typical regions differing by less than $20$ percent at high redshift. We conclude that EPS theory predicts the hierarchical build-up of dark matter haloes quite well over the entire mass range. △ Less

Submitted 28 March, 2023; originally announced March 2023.

SKA Science Data Challenge 2: analysis and results

Authors: P. Hartley, A. Bonaldi, R. Braun, J. N. H. S. Aditya, S. Aicardi, L. Alegre, A. Chakraborty, X. Chen, S. Choudhuri, A. O. Clarke, J. Coles, J. S. Collinson, D. Cornu, L. Darriba, M. Delli Veneri, J. Forbrich, B. Fraga, A. Galan, J. Garrido, F. Gubanov, H. Håkansson, M. J. Hardcastle, C. Heneka, D. Herranz, K. M. Hess , et al. (83 additional authors not shown)

Abstract: The Square Kilometre Array Observatory (SKAO) will explore the radio sky to new depths in order to conduct transformational science. SKAO data products made available to astronomers will be correspondingly large and complex, requiring the application of advanced analysis techniques to extract key science findings. To this end, SKAO is conducting a series of Science Data Challenges, each designed t… ▽ More The Square Kilometre Array Observatory (SKAO) will explore the radio sky to new depths in order to conduct transformational science. SKAO data products made available to astronomers will be correspondingly large and complex, requiring the application of advanced analysis techniques to extract key science findings. To this end, SKAO is conducting a series of Science Data Challenges, each designed to familiarise the scientific community with SKAO data and to drive the development of new analysis techniques. We present the results from Science Data Challenge 2 (SDC2), which invited participants to find and characterise 233245 neutral hydrogen (Hi) sources in a simulated data product representing a 2000~h SKA MID spectral line observation from redshifts 0.25 to 0.5. Through the generous support of eight international supercomputing facilities, participants were able to undertake the Challenge using dedicated computational resources. Alongside the main challenge, `reproducibility awards' were made in recognition of those pipelines which demonstrated Open Science best practice. The Challenge saw over 100 participants develop a range of new and existing techniques, with results that highlight the strengths of multidisciplinary and collaborative effort. The winning strategy -- which combined predictions from two independent machine learning techniques to yield a 20 percent improvement in overall performance -- underscores one of the main Challenge outcomes: that of method complementarity. It is likely that the combination of methods in a so-called ensemble approach will be key to exploiting very large astronomical datasets. △ Less

Submitted 14 March, 2023; originally announced March 2023.

Comments: Under review by MNRAS; 28 pages, 16 figures

Dark energy and matter interacting scenario to relieve $H_0$ and $S_8$ tensions

Authors: Li-Yang Gao, She-Sheng Xue, Xin Zhang

Abstract: We consider a new cosmological model (called $\tildeΛ$CDM) in which the vacuum energy interacts with matter and radiation, and test this model using the current cosmological observations. Using the CMB+BAO+SN (CBS) dataset to constrain the model, we find that $H_0$ and $S_8$ tensions are relieved to $2.87σ$ and $2.77σ$, respectively. However, in this case, the $\tildeΛ$CDM model is not favored by… ▽ More We consider a new cosmological model (called $\tildeΛ$CDM) in which the vacuum energy interacts with matter and radiation, and test this model using the current cosmological observations. Using the CMB+BAO+SN (CBS) dataset to constrain the model, we find that $H_0$ and $S_8$ tensions are relieved to $2.87σ$ and $2.77σ$, respectively. However, in this case, the $\tildeΛ$CDM model is not favored by the data, compared with $Λ$CDM. We find that when the $H_0$ and $S_8$ data are added to the data combination, the situation is significantly improved. In the CBS+$H_0$ case, the model relieves the $H_0$ tension to $0.47σ$, and the model is favored over $Λ$CDM. In the CBS+$H_0$+$S_8$ case, we obtain a synthetically best situation, in which the $H_0$ and $S_8$ tensions are relieved to $0.72σ$ and $2.11σ$, respectively. In this case, the model is most favored by the data. Therefore, this cosmological model can greatly relieve the $H_0$ tension and simultaneously effectively alleviate the $S_8$ tension. △ Less

Submitted 28 March, 2024; v1 submitted 26 December, 2022; originally announced December 2022.

Comments: 7 pages, 4 figures

Journal ref: Chinese Physics C Vol. 48, No. 6 (2024) 051001

The halo concentration and mass relation traced by satellite galaxies

Authors: Qing Gu, Qi Guo, Tianchi Zhang, Wenting Wang, Quan Guo, Liang Gao

Abstract: We study the relation between halo concentration and mass (c-M relation) using the Seventh and Eighth Data Release of the Sloan Digital Sky Survey (SDSS DR7 and DR8) galaxy catalogue. Assuming that the satellite galaxies follow the distribution of dark matter, we derive the halo concentration by fitting the satellite radial profile with a Nararro Frank and White (NFW) format. The derived c-M relat… ▽ More We study the relation between halo concentration and mass (c-M relation) using the Seventh and Eighth Data Release of the Sloan Digital Sky Survey (SDSS DR7 and DR8) galaxy catalogue. Assuming that the satellite galaxies follow the distribution of dark matter, we derive the halo concentration by fitting the satellite radial profile with a Nararro Frank and White (NFW) format. The derived c-M relation covers a wide halo mass range from $10^{11.6}$ to $10^{14.1} \rm\ M_\odot$. We confirm the anti-correlation between the halo mass and concentration as predicted in cosmological simulations. Our results are in good agreement with those derived using galaxy dynamics and gravitational lensing for halos of $10^{11.6}-10^{12.9} \rm\ M_\odot$, while they are slightly lower for halos of $10^{12.9}-10^{14.1}\rm\ M_\odot$. It is because blue satellite galaxies are less concentrated, especially in the inner regions. Instead of using all satellite galaxies, red satellites could be better tracers of the underlying dark matter distribution in galaxy groups. △ Less

Submitted 20 December, 2022; originally announced December 2022.

Comments: 8 pages, 4 figures, accepted for publication in RAA

Eliminating polarization leakage effect for neutral hydrogen intensity mapping with deep learning

Authors: Li-Yang Gao, Yichao Li, Shulei Ni, Xin Zhang

Abstract: The neutral hydrogen (HI) intensity mapping (IM) survey is regarded as a promising approach for cosmic large-scale structure (LSS) studies. A major issue for the HI IM survey is to remove the bright foreground contamination. A key to successfully remove the bright foreground is to well control or eliminate the instrumental effects. In this work, we consider the instrumental effect of polarization… ▽ More The neutral hydrogen (HI) intensity mapping (IM) survey is regarded as a promising approach for cosmic large-scale structure (LSS) studies. A major issue for the HI IM survey is to remove the bright foreground contamination. A key to successfully remove the bright foreground is to well control or eliminate the instrumental effects. In this work, we consider the instrumental effect of polarization leakage and use the U-Net approach, a deep learning-based foreground removal technique, to eliminate the polarization leakage effect. The thermal noise is assumed to be a subdominant factor compared with the polarization leakage for future HI IM surveys and ignored in this analysis. In this method, the principal component analysis (PCA) foreground subtraction is used as a preprocessing step for the U-Net foreground subtraction. Our results show that the additional U-Net processing could either remove the foreground residual after the conservative PCA subtraction or compensate for the signal loss caused by the aggressive PCA preprocessing. Finally, we test the robustness of the U-Net foreground subtraction technique and show that it is still reliable in the case of existing constraint error on HI fluctuation amplitude. △ Less

Submitted 4 September, 2023; v1 submitted 16 December, 2022; originally announced December 2022.

Comments: 13 pages, 13 figures; accepted for publication in MNRAS

Journal ref: MNRAS 525, 5278-5290 (2023)

Flux Variations of Cosmic Ray Air Showers Detected by LHAASO-KM2A During a Thunderstorm on 10 June 2021

Authors: LHAASO Collaboration, F. Aharonian, Q. An, Axikegu, L. X. Bai, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Zhe Cao, Zhen Cao, J. Chang, J. F. Chang, E. S. Chen, Liang Chen, Liang Chen, Long Chen, M. J. Chen, M. L. Chen, S. H. Chen, S. Z. Chen, T. L. Chen, X. J. Chen , et al. (248 additional authors not shown)

Abstract: The Large High Altitude Air Shower Observatory (LHAASO) has three sub-arrays, KM2A, WCDA and WFCTA. The flux variations of cosmic ray air showers were studied by analyzing the KM2A data during the thunderstorm on 10 June 2021. The number of shower events that meet the trigger conditions increases significantly in atmospheric electric fields, with maximum fractional increase of 20%. The variations… ▽ More The Large High Altitude Air Shower Observatory (LHAASO) has three sub-arrays, KM2A, WCDA and WFCTA. The flux variations of cosmic ray air showers were studied by analyzing the KM2A data during the thunderstorm on 10 June 2021. The number of shower events that meet the trigger conditions increases significantly in atmospheric electric fields, with maximum fractional increase of 20%. The variations of trigger rates (increases or decreases) are found to be strongly dependent on the primary zenith angle. The flux of secondary particles increases significantly, following a similar trend with that of the shower events. To better understand the observed behavior, Monte Carlo simulations are performed with CORSIKA and G4KM2A (a code based on GEANT4). We find that the experimental data (in saturated negative fields) are in good agreement with simulations, assuming the presence of a uniform upward electric field of 700 V/cm with a thickness of 1500 m in the atmosphere above the observation level. Due to the acceleration/deceleration and deflection by the atmospheric electric field, the number of secondary particles with energy above the detector threshold is modified, resulting in the changes in shower detection rate. △ Less

Submitted 6 December, 2022; v1 submitted 25 July, 2022; originally announced July 2022.

Comments: 18 pages, 11 figures

Journal ref: Chinese Phys. C 47 015001 (2023)

The Ultramarine Simulation: properties of dark matter haloes before redshift 5.5

Authors: Qiao Wang, Liang Gao, Cheng Meng

Abstract: We introduce the Ultramarine simulation, an extremely large $N$-body simulation of the structure formation and evolution to redshift 5.5 at which cosmic reionization was just completed. The simulation evolves 2.1 trillion particles within a $512~h^{-1}$Mpc cube and has an unprecedented mass and force resolution for large volume simulations of this kind, $5.6\times 10^6 h^{-1}$M$_\odot$ and… ▽ More We introduce the Ultramarine simulation, an extremely large $N$-body simulation of the structure formation and evolution to redshift 5.5 at which cosmic reionization was just completed. The simulation evolves 2.1 trillion particles within a $512~h^{-1}$Mpc cube and has an unprecedented mass and force resolution for large volume simulations of this kind, $5.6\times 10^6 h^{-1}$M$_\odot$ and $1.2~h^{-1}$kpc, respectively. We present some basic statistical results of the simulation, including the halo mass function, halo bias parameter as well as halo mass-concentration relation at high redshifts, and compare them with some existing representative models. We find excellent agreement with some models on the high redshift halo mass functions, but neither the halo bias factor nor halo mass-concentration relation. All halo bias models for comparison over-predicate high redshift halo bias by large factors, an accurate fit to our simulation is given. High redshift dark matter haloes still can be reasonably described with NFW model, the halo mass-concentration relations are monotonic, with more massive haloes having lower concentration, in disfavor of the upturn feature reported by some studies. The mass concentration relation has little evolution between $z=5.5$ to $z=10$, in contrast to strong evolution predicted by most existing models. In addition, concentration parameters of high redshift dark matter haloes are much lower than most model predictions. △ Less

Submitted 23 October, 2022; v1 submitted 13 June, 2022; originally announced June 2022.

Comments: 9 pages, 8 figures. Accepted for publication in MNRAS

Eliminating Primary Beam Effect in Foreground Subtraction of Neutral Hydrogen Intensity Mapping Survey with Deep Learning

Authors: Shulei Ni, Yichao Li, Li-Yang Gao, Xin Zhang

Abstract: In the neutral hydrogen (HI) intensity mapping (IM) survey, the foreground contamination on the cosmological signals is extremely severe, and the systematic effects caused by radio telescopes themselves further aggravate the difficulties in subtracting foreground. In this work, we investigate whether the deep learning method, concretely the 3D U-Net algorithm here, can play a crucial role in foreg… ▽ More In the neutral hydrogen (HI) intensity mapping (IM) survey, the foreground contamination on the cosmological signals is extremely severe, and the systematic effects caused by radio telescopes themselves further aggravate the difficulties in subtracting foreground. In this work, we investigate whether the deep learning method, concretely the 3D U-Net algorithm here, can play a crucial role in foreground subtraction when considering the systematic effect caused by the telescope's primary beam. We consider two beam models, i.e., the Gaussian beam model as a simple case and the Cosine beam model as a sophisticated case. The traditional principal component analysis (PCA) method is employed as a comparison and, more importantly, as the preprocessing step for the U-Net method to reduce the sky map dynamic range. We find that in the case of the Gaussian beam, the PCA method can effectively clean the foreground. However, the PCA method cannot handle the systematic effect induced by the Cosine beam, and the additional U-Net method can improve the result significantly. In order to show how well the PCA and U-Net methods can recover the HI signals, we also derive the HI angular power spectra, as well as the HI 2D power spectra, after performing the foreground subtractions. It is found that, in the case of Gaussian beam, the concordance with the original HI map using U-Net is better than that using PCA by $27.4\%$, and in the case of Cosine beam, the concordance using U-Net is better than that using PCA by $144.8\%$. Therefore, the U-Net based foreground subtraction can efficiently eliminate the telescope primary beam effect and shed new light on recovering the HI power spectrum for future HI IM experiments. △ Less

Submitted 31 May, 2022; v1 submitted 5 April, 2022; originally announced April 2022.

Comments: 19 pages, 13 figures

Journal ref: The Astrophysical Journal 934, 83 (2022)

Meshless Methods for Magnetohydrodynamics with Vector Potential

Authors: Xiongbiao Tu, Qiao Wang, Haonan Zheng, Liang Gao

Abstract: We present a meshless method for magnetohydrodynamics by evolving the vector potential of magnetic fields. A novel scheme and numerical techniques are developed to restrict the divergence of magnetic field, based on the Meshless Finite Mass/Volume with HLLD Riemann solver for conservative flux calculation. We found the magnetic field could be stabilized by a proper smoothing process and so the lon… ▽ More We present a meshless method for magnetohydrodynamics by evolving the vector potential of magnetic fields. A novel scheme and numerical techniques are developed to restrict the divergence of magnetic field, based on the Meshless Finite Mass/Volume with HLLD Riemann solver for conservative flux calculation. We found the magnetic field could be stabilized by a proper smoothing process and so the long-term evolution becomes available. To verify the new scheme, we perform the Brio-Wu shock tube, 2D and 3D Orszag-Tang vortex and Magnetorotational Instability test problems. Our results suggest that our method is robust and has better precision on central offset, amplitude and detailed pattern than an existing meshless code$-$GIZMO. △ Less

Submitted 3 September, 2022; v1 submitted 8 February, 2022; originally announced February 2022.

Comments: 17 pages, 11 figures, accepted for publication in Journal of Computational Physics

Direct measurement of non-thermal electron acceleration from magnetically driven reconnection in a laboratory plasma

Authors: Abraham Chien, Lan Gao, Shu Zhang, Hantao Ji, Eric G. Blackman, William Daughton, Adam Stanier, Ari Le, Fan Guo, Russ Follett, Hui Chen, Gennady Fiksel, Gabriel Bleotu, Robert C. Cauble, Sophia N. Chen, Alice Fazzini, Kirk Flippo, Omar French, Dustin H. Froula, Julien Fuchs, Shinsuke Fujioka, Kenneth Hill, Sallee Klein, Carolyn Kuranz, Philip Nilson , et al. (2 additional authors not shown)

Abstract: Magnetic reconnection is a ubiquitous astrophysical process that rapidly converts magnetic energy into some combination of plasma flow energy, thermal energy, and non-thermal energetic particles, including energetic electrons. Various reconnection acceleration mechanisms in different low-$β$ (plasma-to-magnetic pressure ratio) and collisionless environments have been proposed theoretically and stu… ▽ More Magnetic reconnection is a ubiquitous astrophysical process that rapidly converts magnetic energy into some combination of plasma flow energy, thermal energy, and non-thermal energetic particles, including energetic electrons. Various reconnection acceleration mechanisms in different low-$β$ (plasma-to-magnetic pressure ratio) and collisionless environments have been proposed theoretically and studied numerically, including first- and second-order Fermi acceleration, betatron acceleration, parallel electric field acceleration along magnetic fields, and direct acceleration by the reconnection electric field. However, none of them have been heretofore confirmed experimentally, as the direct observation of non-thermal particle acceleration in laboratory experiments has been difficult due to short Debye lengths for \textit{in-situ} measurements and short mean free paths for \textit{ex-situ} measurements. Here we report the direct measurement of accelerated non-thermal electrons from low-$β$ magnetically driven reconnection in experiments using a laser-powered capacitor coil platform. We use kiloJoule lasers to drive parallel currents to reconnect MegaGauss-level magnetic fields in a quasi-axisymmetric geometry. The angular dependence of the measured electron energy spectrum and the resulting accelerated energies, supported by particle-in-cell simulations, indicate that the mechanism of direct electric field acceleration by the out-of-plane reconnection electric field is at work. Scaled energies using this mechanism show direct relevance to astrophysical observations. Our results therefore validate one of the proposed acceleration mechanisms by reconnection, and establish a new approach to study reconnection particle acceleration with laboratory experiments in relevant regimes. △ Less

Submitted 24 January, 2022; originally announced January 2022.

Comments: 13 pages, 6 figures

Peta-electron volt gamma-ray emission from the Crab Nebula

Authors: The LHAASO Collaboration, Zhen Cao, F. Aharonian, Q. An, Axikegu, L. X. Bai, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, H. Cai, J. T. Cai, Zhe Cao, J. Chang, J. F. Chang, B. M. Chen, E. S. Chen, J. Chen, Liang Chen, Liang Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen , et al. (250 additional authors not shown)

Abstract: The Crab pulsar and the surrounding nebula powered by the pulsar's rotational energy through the formation and termination of a relativistic electron-positron wind is a bright source of gamma-rays carrying crucial information about this complex conglomerate. We report the detection of $γ$-rays with a spectrum showing gradual steepening over three energy decades, from $5\times 10^{-4}$ to $1.1$ pet… ▽ More The Crab pulsar and the surrounding nebula powered by the pulsar's rotational energy through the formation and termination of a relativistic electron-positron wind is a bright source of gamma-rays carrying crucial information about this complex conglomerate. We report the detection of $γ$-rays with a spectrum showing gradual steepening over three energy decades, from $5\times 10^{-4}$ to $1.1$ petaelectronvolt (PeV). The ultra-high-energy photons exhibit the presence of a PeV electron accelerator (a pevatron) with an acceleration rate exceeding 15% of the absolute theoretical limit. Assuming that unpulsed $γ$-rays are produced at the termination of the pulsar's wind, we constrain the pevatron's size, between $0.025$ and $0.1$ pc, and the magnetic field $\approx 110 μ$G. The production rate of PeV electrons, $2.5 \times 10^{36}$ erg $\rm s^{-1}$, constitutes 0.5% of the pulsar's spin-down luminosity, although we do not exclude a non-negligible contribution of PeV protons to the production of the highest energy $γ$-rays. △ Less

Submitted 11 November, 2021; originally announced November 2021.

Comments: 43 pages, 13 figures, 2 tables; Published in Science

Journal ref: Science, 2021, Vol 373, Issue 6553, pp. 425-430

The Formation and Evolution of Massive Galaxies

Authors: Yingjie Jing, Yu Rong, Jie Wang, Qi Guo, Liang Gao

Abstract: The discovery of massive galaxies at high redshifts, especially the passive ones, poses a big challenge for the current standard galaxy formation models. Here we use the semi-analytic galaxy formation model developed by Henriques et al. to explore the formation and evolution of massive galaxies (MGs, stellar-mass $M_{*}> 10^{11}$ M$_{\odot}$). Different from previous works, we focus on the ones ju… ▽ More The discovery of massive galaxies at high redshifts, especially the passive ones, poses a big challenge for the current standard galaxy formation models. Here we use the semi-analytic galaxy formation model developed by Henriques et al. to explore the formation and evolution of massive galaxies (MGs, stellar-mass $M_{*}> 10^{11}$ M$_{\odot}$). Different from previous works, we focus on the ones just formed (e.g. just reach $\simeq 10^{11}$ M$_{\odot}$). We find that most of the MGs are formed around $z=0.6$, with the earliest formation at $z>4$. Interestingly, although most of the MGs in the local Universe are passive, we find that only $13\%$ of the MGs are quenched at the formation time. Most of the quenched MGs at formation already hosts a very massive supermassive black hole (SMBH) which could power the very effective AGN feedback. For the star-forming MGs, the ones with more massive SMBH prefer to quench in shorter timescales; in particular, those with $M_{\textrm{SMBH}} > 10^{7.5}$ M$_{\odot}$ have a quenching timescale of $\sim 0.5$ Gyr and the characteristic $M_{\textrm{SMBH}}$ depends on the chosen stellar mass threshold in the definition of MGs as a result of their co-evolution. We also find that the "in-situ" star formation dominates the stellar mass growth of MGs until they are formed. Over the whole redshift range, we find the quiescent MGs prefer to stay in more massive dark matter halos, and have more massive SMBH and less cold gas masses. Our results provide a new angle on the whole life of the growth of MGs in the Universe. △ Less

Submitted 1 November, 2021; originally announced November 2021.

Comments: 11 pages, 11 figures; published in RAA

The galaxy size to halo spin relation of disk galaxies in cosmological hydrodynamical simulations

Authors: Hang Yang, Liang Gao, Carlos S. Frenk, Robert J. J. Grand, Qi Guo, Shihong Liao, Shi Shao

Abstract: In the standard disk galaxy formation model, the sizes of galactic disks are tightly related to the spin parameters $λ$ of their dark matter haloes. The model has been wildly adopted by various semi-analytic galaxy formation models which have been extremely successful to interpret a large body of observational data. However, the size-$λ$ correlation was rarely seen in most modern hydrodynamical si… ▽ More In the standard disk galaxy formation model, the sizes of galactic disks are tightly related to the spin parameters $λ$ of their dark matter haloes. The model has been wildly adopted by various semi-analytic galaxy formation models which have been extremely successful to interpret a large body of observational data. However, the size-$λ$ correlation was rarely seen in most modern hydrodynamical simulations of galaxy formation. In this short paper, we make use of 4 sets of large hydrodynamical simulations to explore the size-spin parameter relation with a large sample of simulated disk galaxies, and compare it with a popular disk galaxy formation model of Mo et al. (1998). Intriguingly, galactic sizes correlate strongly with spin parameters of their dark matter haloes in the simulations developed by the IllustrisTNG collaborations, albeit the relation does not always agree with prediction of MMW98 model over all stellar mass range we examined. There is also a size-spin correlation for the Milky way analogies in the EAGLE simulations, while it is relatively weaker than that of the IllustrisTNG counterparts. For the dwarfs in the simulations from the EAGLE collaboration, there is NULL correlation. We conclude that either the detailed subgrid physics or hydrodynamics solvers account for the size-spin relation, which will be explored in our future work. △ Less

Submitted 10 March, 2023; v1 submitted 8 October, 2021; originally announced October 2021.

Journal ref: 2023MNRAS.518.5253Y

About One-point Statistics of the Ratio of Two Fourier-transformed Cosmic Fields and an Application

Authors: Ming Li, Jun Pan, Pengjie Zhang, Jie Wang, Longlong Feng, Liang Gao, Xi Kang, Guoliang Li, Weipeng Lin, Haihui Wang

Abstract: The Fourier transformation is an effective and efficient operation of Gaussianization at the one-point level. Using a set of N-body simulation data, we verified that the one-point distribution functions of the dark matter momentum divergence and density fields closely follow complex Gaussian distributions. The one-point distribution function of the quotient of two complex Gaussian variables is int… ▽ More The Fourier transformation is an effective and efficient operation of Gaussianization at the one-point level. Using a set of N-body simulation data, we verified that the one-point distribution functions of the dark matter momentum divergence and density fields closely follow complex Gaussian distributions. The one-point distribution function of the quotient of two complex Gaussian variables is introduced and studied. Statistical theories are then applied to model one-point statistics about the growth of individual Fourier mode of the dark matter density field, which can be obtained by the ratio of two Fourier transformed cosmic fields. Our simulation results proved that the models based on the Gaussian approximation are impressively accurate, and our analysis revealed many interesting aspects about the growth of dark matter's density fluctuation in Fourier space. △ Less

Submitted 3 July, 2022; v1 submitted 30 September, 2021; originally announced September 2021.

Comments: 15 pages, 13 figures. Accepted by ApJ

The impact of filaments on dwarf galaxy properties in the Auriga simulations

Authors: Haonan Zheng, Shihong Liao, Jia Hu, Liang Gao, Robert J. J. Grand, Qing Gu, Qi Guo

Abstract: With a hydrodynamical simulation using a simple galaxy formation model without taking into account feedback, our previous work has shown that dense and massive filaments at high redshift can provide potential wells to trap and compress gas, and hence affect galaxy formation in their resident low-mass haloes. In this paper, we make use of the Auriga simulations, a suite of high-resolution zoom-in h… ▽ More With a hydrodynamical simulation using a simple galaxy formation model without taking into account feedback, our previous work has shown that dense and massive filaments at high redshift can provide potential wells to trap and compress gas, and hence affect galaxy formation in their resident low-mass haloes. In this paper, we make use of the Auriga simulations, a suite of high-resolution zoom-in hydrodynamical simulations of Milky Way-like galaxies, to study whether the conclusion still holds in the simulations with a sophisticated galaxy formation model. In agreement with the results of our previous work, we find that, comparing to their counterparts with similar halo masses in field, dwarf galaxies residing in filaments tend to have higher baryonic and stellar fractions. At the fixed parent halo mass, the filament dwarfs tend to have slightly higher star formation rates than those of field ones. But overall we do not find a clear difference in galaxy g - r colours between the filament and field populations. We also show that at high redshifts, the gas components in dwarf galaxies tend to have their spins aligned with the filaments in which they reside. Our results support a picture in which massive filaments at high redshift assist gas accretion and enhance star formation in their resident dwarf sized dark matter haloes. △ Less

Submitted 27 May, 2022; v1 submitted 17 September, 2021; originally announced September 2021.

Comments: 9 pages, 7 figures, accepted for publication in MNRAS

Journal ref: Monthly Notice of the Royal Astronomical Society, Volume 514, Issue 2, August 2022, Pages 2488-2496

The role of mergers and gas accretion in black hole growth and galaxy evolution

Authors: TianChi Zhang, Qi Guo, Yan Qu, Liang Gao

Abstract: We use a semi-analytic galaxy formation model to study the co-evolution of supermassive black holes (SMBHs) with their host galaxies. Although the coalescence of SMBHs is not important, the quasar-mode accretion induced by mergers plays a dominant role in the growth of SMBHs. Mergers play a more important role in the growth of SMBH host galaxies than in the SMBH growth. It is the combined contribu… ▽ More We use a semi-analytic galaxy formation model to study the co-evolution of supermassive black holes (SMBHs) with their host galaxies. Although the coalescence of SMBHs is not important, the quasar-mode accretion induced by mergers plays a dominant role in the growth of SMBHs. Mergers play a more important role in the growth of SMBH host galaxies than in the SMBH growth. It is the combined contribution from quasar mode accretion and mergers to the SMBH growth and the combined contribution from starburst and mergers to their host galaxy growth that determine the observed scaling relation between the SMBH masses and their host galaxy masses. We also find that mergers are more important in the growth of SMBH host galaxies compared to normal galaxies which share the same stellar mass range as the SMBH host galaxies. △ Less

Submitted 7 September, 2021; originally announced September 2021.

Comments: 10 pages, 5 figures, version published by RAA

Journal ref: RAA 2021 Vol. 21 No. 8, 212(10pp)

Formation of massive disk galaxies in the IllustrisTNG simulation

Authors: Guangquan Zeng, Lan Wang, Liang Gao

Abstract: We investigate the formation history of massive disk galaxies in hydro-dynamical simulation--the IllustrisTNG, to study why massive disk galaxies survive through cosmic time. 83 galaxies in the simulation are selected with M$_{*,z=0}$ $>8\times10^{10}$ M$_\odot$ and kinematic bulge-to-total ratio less than $0.3$. We find that 8.4 percent of these massive disk galaxies have quiet merger histories a… ▽ More We investigate the formation history of massive disk galaxies in hydro-dynamical simulation--the IllustrisTNG, to study why massive disk galaxies survive through cosmic time. 83 galaxies in the simulation are selected with M$_{*,z=0}$ $>8\times10^{10}$ M$_\odot$ and kinematic bulge-to-total ratio less than $0.3$. We find that 8.4 percent of these massive disk galaxies have quiet merger histories and preserve disk morphology since formed. 54.2 percent have a significant increase in bulge components in history, then become disks again till present time. The rest 37.3 percent experience prominent mergers but survive to remain disky. While mergers and even major mergers do not always turn disk galaxies into ellipticals, we study the relations between various properties of mergers and the morphology of merger remnants. We find a strong dependence of remnant morphology on the orbit type of major mergers. Specifically, major mergers with a spiral-in falling orbit mostly lead to disk-dominant remnants, and major mergers of head-on galaxy-galaxy collision mostly form ellipticals. This dependence of remnant morphology on orbit type is much stronger than the dependence on cold gas fraction or orbital configuration of merger system as previously studied. △ Less

Submitted 8 September, 2021; v1 submitted 20 May, 2021; originally announced May 2021.

Comments: 12 pages, 1 table, 8 figures. Accepted for publication in MNRAS

Relieving the $H_0$ tension with a new interacting dark energy model

Authors: Li-Yang Gao, Ze-Wei Zhao, She-Sheng Xue, Xin Zhang

Abstract: We investigate an extended cosmological model motivated by the asymptotic safety of gravitational field theory, in which the matter and radiation densities and the cosmological constant receive a correction parametrized by the parameters $δ_G$ and $δ_Λ$, leading to that both the evolutions of the matter and radiation densities and the cosmological constant slightly deviate from the standard forms.… ▽ More We investigate an extended cosmological model motivated by the asymptotic safety of gravitational field theory, in which the matter and radiation densities and the cosmological constant receive a correction parametrized by the parameters $δ_G$ and $δ_Λ$, leading to that both the evolutions of the matter and radiation densities and the cosmological constant slightly deviate from the standard forms. Here we explain this model as a scenario of vacuum energy interacting with matter and radiation. We consider two cases of the model: {(i) ${\tildeΛ}$CDM with one additional free parameter $δ_G$, with $δ_{\rm G}$ and $δ_Λ$ related by a low-redshift limit relation and (ii) e${\tildeΛ}$CDM with two additional free parameters $δ_G$ and $δ_Λ$ that are independent of each other.} We use two data combinations, CMB+BAO+SN (CBS) and CMB+BAO+SN+$H_0$ (CBSH), to constrain the models. We find that, in the case of using the CBS data, neither ${\tildeΛ}$CDM nor e${\tildeΛ}$CDM can effectively alleviate the $H_0$ tension. However, it is found that using the CBSH data the $H_0$ tension can be greatly relieved by the models. In particular, in the case of e${\tildeΛ}$CDM, the $H_0$ tension can be resolved to 0.71$σ$. We conclude that as an interacting dark energy model, ${\tildeΛ}$CDM is much better than $Λ(t)$CDM in the sense of both relieving the $H_0$ tension and fitting to the current observational data. △ Less

Submitted 26 May, 2021; v1 submitted 26 January, 2021; originally announced January 2021.

Comments: 17 pages, 6 figures; accepted for publication in JCAP

Journal ref: JCAP 07 (2021) 005

Major Scientific Challenges and Opportunities in Understanding Magnetic Reconnection and Related Explosive Phenomena in Solar and Heliospheric Plasmas

Authors: H. Ji, J. Karpen, A. Alt, S. Antiochos, S. Baalrud, S. Bale, P. M. Bellan, M. Begelman, A. Beresnyak, A. Bhattacharjee, E. G. Blackman, D. Brennan, M. Brown, J. Buechner, J. Burch, P. Cassak, B. Chen, L. -J. Chen, Y. Chen, A. Chien, L. Comisso, D. Craig, J. Dahlin, W. Daughton, E. DeLuca , et al. (83 additional authors not shown)

Abstract: Magnetic reconnection underlies many explosive phenomena in the heliosphere and in laboratory plasmas. The new research capabilities in theory/simulations, observations, and laboratory experiments provide the opportunity to solve the grand scientific challenges summarized in this whitepaper. Success will require enhanced and sustained investments from relevant funding agencies, increased interagen… ▽ More Magnetic reconnection underlies many explosive phenomena in the heliosphere and in laboratory plasmas. The new research capabilities in theory/simulations, observations, and laboratory experiments provide the opportunity to solve the grand scientific challenges summarized in this whitepaper. Success will require enhanced and sustained investments from relevant funding agencies, increased interagency/international partnerships, and close collaborations of the solar, heliospheric, and laboratory plasma communities. These investments will deliver transformative progress in understanding magnetic reconnection and related explosive phenomena including space weather events. △ Less

Submitted 16 September, 2020; originally announced September 2020.

Comments: 4 pages (including a title page), white paper submitted to Helio2050 workshop at https://www.hou.usra.edu/meetings/helio2050/. arXiv admin note: substantial text overlap with arXiv:2004.00079

Galaxy properties in the cosmic web of EAGLE simulation

Authors: Wenxiao Xu, Qi Guo, Haonan Zheng, Liang Gao, Cedric Lacey, Qing Gu, Shihong Liao, Shi Shao, Tianxiang Mao, Tianchi Zhang, Xuelei Chen

Abstract: We investigate the dependence of the galaxy properties on cosmic web environments using the most up-to-date hydrodynamic simulation: Evolution and Assembly of Galaxies and their Environments (EAGLE). The baryon fractions in haloes and the amplitudes of the galaxy luminosity function decrease going from knots to filaments to sheets to voids. Interestingly, the value of L$^*$ varies dramatically in… ▽ More We investigate the dependence of the galaxy properties on cosmic web environments using the most up-to-date hydrodynamic simulation: Evolution and Assembly of Galaxies and their Environments (EAGLE). The baryon fractions in haloes and the amplitudes of the galaxy luminosity function decrease going from knots to filaments to sheets to voids. Interestingly, the value of L$^*$ varies dramatically in different cosmic web environments. At z = 0, we find a characteristic halo mass of $10^{12} h^{-1}\rm M_{\odot}$, below which the stellar-to-halo mass ratio is higher in knots while above which it reverses. This particular halo mass corresponds to a characteristic stellar mass of $1.8\times 10^{10} h^{-1}\rm M_{\odot}$. Below the characteristic stellar mass central galaxies have redder colors, lower sSFRs and higher metallicities in knots than those in filaments, sheets and voids, while above this characteristic stellar mass, the cosmic web environmental dependences either reverse or vanish. Such dependences can be attributed to the fact that the active galaxy fraction decreases along voids, sheets, filaments and knots. The cosmic web dependences get weaker towards higher redshifts for most of the explored galaxy properties and scaling relations, except for the gas metallicity vs. stellar mass relation. △ Less

Submitted 15 September, 2020; originally announced September 2020.

Comments: 14 pages, 14 figures, 1 table. Accepted for publication in MNRAS

LESSER: A Catalogue of Spectroscopically selected sample of Lyman-$α$ Emitters Lensed By Galaxies

Authors: Xiaoyue Cao, Ran Li, Yiping Shu, Shude Mao, Jean-Paul Kneib, Liang Gao

Abstract: We introduce the LEnSed laeS in the Eboss suRvey (LESSER) project, which aims to search for lensed Lyman-$α$ Emitters (LAEs) in the Extended Baryon Oscillation Spectroscopic Survey (eBOSS). The final catalog contains 361 candidate lensing systems. The lens galaxies are luminous red galaxies (LRGs) at redshift $0.4 < z < 0.8$, and the source galaxies are LAEs at redshift $2 < z < 3$. The spectral r… ▽ More We introduce the LEnSed laeS in the Eboss suRvey (LESSER) project, which aims to search for lensed Lyman-$α$ Emitters (LAEs) in the Extended Baryon Oscillation Spectroscopic Survey (eBOSS). The final catalog contains 361 candidate lensing systems. The lens galaxies are luminous red galaxies (LRGs) at redshift $0.4 < z < 0.8$, and the source galaxies are LAEs at redshift $2 < z < 3$. The spectral resolution of eBOSS ($\sim$2000) allows us to further identify the fine structures of Lyman-$α$ ($\rm Lyα$) emissions. Among our lensed LAE candidates, 281 systems present single-peaked line profiles while 80 systems show double-peaked features. Future spectroscopic/imaging follow-up observations of the catalog may shed light on the origin of diverse $\rm Lyα$ line morphology, and provide promising labs for studying low mass dark matter haloes/subhaloes. △ Less

Submitted 18 October, 2020; v1 submitted 4 August, 2020; originally announced August 2020.

Comments: 11 Pages, 9 figures, 2 tables, accepted by MNRAS. For catalog data, see https://github.com/caoxiaoyue/search_lensed_laes/tree/master/paper_material/catalog ; For spectra images, see https://github.com/caoxiaoyue/search_lensed_laes/tree/master/paper_material/combine_spectra ; For Lens images, see https://github.com/caoxiaoyue/search_lensed_laes/tree/master/paper_material/image_stamp

A Fast Radio Burst discovered in FAST drift scan survey

Authors: Weiwei Zhu, Di Li, Rui Luo, Chenchen Miao, Bing Zhang, Laura Spitler, Duncan Lorimer, Michael Kramer, David Champion, Youling Yue, Andrew Cameron, Marilyn Cruces, Ran Duan, Yi Feng, Jun Han, George Hobbs, Chenhui Niu, Jiarui Niu, Zhichen Pan, Lei Qian, Dai Shi, Ningyu Tang, Pei Wang, Hongfeng Wang, Mao Yuan , et al. (45 additional authors not shown)

Abstract: We report the discovery of a highly dispersed fast radio burst, FRB~181123, from an analysis of $\sim$1500~hr of drift-scan survey data taken using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The pulse has three distinct emission components, which vary with frequency across our 1.0--1.5~GHz observing band. We measure the peak flux density to be $>0.065$~Jy and the correspondi… ▽ More We report the discovery of a highly dispersed fast radio burst, FRB~181123, from an analysis of $\sim$1500~hr of drift-scan survey data taken using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The pulse has three distinct emission components, which vary with frequency across our 1.0--1.5~GHz observing band. We measure the peak flux density to be $>0.065$~Jy and the corresponding fluence $>0.2$~Jy~ms. Based on the observed dispersion measure of 1812~cm$^{-3}$~pc, we infer a redshift of $\sim 1.9$. From this, we estimate the peak luminosity and isotropic energy to be $\lesssim 2\times10^{43}$~erg~s$^{-1}$ and $\lesssim 2\times10^{40}$~erg, respectively. With only one FRB from the survey detected so far, our constraints on the event rate are limited. We derive a 95\% confidence lower limit for the event rate of 900 FRBs per day for FRBs with fluences $>0.025$~Jy~ms. We performed follow-up observations of the source with FAST for four hours and have not found a repeated burst. We discuss the implications of this discovery for our understanding of the physical mechanisms of FRBs. △ Less

Submitted 29 April, 2020; originally announced April 2020.

Comments: 7 pages, 3 figures, accepted for publication in ApJL

The effect of hydrodynamics alone on the subhalo population in a $Λ$CDM rich cluster sized dark matter halo

Authors: Junyi Jia, Liang Gao, Yan Qu

Abstract: We perform a set of non-radiative hydro-dynamical (NHD) simulations of a rich cluster sized dark matter halo from the Phoenix project with 3 different numerical resolutions, to investigate the effect of hydrodynamics alone on the subhalo population in the halo. Compared to dark matter only (DMO) simulations of the same halo, subhaloes are less abundant for relatively massive subhaloes (… ▽ More We perform a set of non-radiative hydro-dynamical (NHD) simulations of a rich cluster sized dark matter halo from the Phoenix project with 3 different numerical resolutions, to investigate the effect of hydrodynamics alone on the subhalo population in the halo. Compared to dark matter only (DMO) simulations of the same halo, subhaloes are less abundant for relatively massive subhaloes ($M_{sub} > 2.5 \times 10^9h^{-1}M_{\odot}$, or $V_{max} > 70 kms^{-1}$) but more abundant for less massive subhaloes in the NHD simulations. This results in different shapes in the subhalo mass/$V_{max}$ function in two different sets of simulations. At given subhalo mass, the subhaloes less massive than $10^{10} h^{-1}M_{\odot}$ have larger $V_{max}$ in the NHD than DMO simulations, while $V_{max}$ is similar for the subhaloes more massive than the mass value. This is mainly because the progenitors of present day low mass subhaloes have larger concentration parameters in the NHD than DMO simulations. The survival number fraction of the accreted low mass progenitors of the main halo at redshift 2 is about 50 percent higher in the NHD than DMO simulations. △ Less

Submitted 29 December, 2020; v1 submitted 2 April, 2020; originally announced April 2020.

Comments: 8 pages, 6 figres

Major Scientific Challenges and Opportunities in Understanding Magnetic Reconnection and Related Explosive Phenomena throughout the Universe

Authors: H. Ji, A. Alt, S. Antiochos, S. Baalrud, S. Bale, P. M. Bellan, M. Begelman, A. Beresnyak, E. G. Blackman, D. Brennan, M. Brown, J. Buechner, J. Burch, P. Cassak, L. -J. Chen, Y. Chen, A. Chien, D. Craig, J. Dahlin, W. Daughton, E. DeLuca, C. F. Dong, S. Dorfman, J. Drake, F. Ebrahimi , et al. (75 additional authors not shown)

Abstract: This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process. This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process. △ Less

Submitted 31 March, 2020; originally announced April 2020.

Comments: 6 pages, 1 figure, white paper submitted to both Plasma 2020 and Astro 2020 Decadal Surveys