Showing 1–50 of 120 results for author: Zhu, H

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

Data quality control system and long-term performance monitor of the LHAASO-KM2A

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

Abstract: The KM2A is the largest sub-array of the Large High Altitude Air Shower Observatory (LHAASO). It consists of 5216 electromagnetic particle detectors (EDs) and 1188 muon detectors (MDs). The data recorded by the EDs and MDs are used to reconstruct primary information of cosmic ray and gamma-ray showers. This information is used for physical analysis in gamma-ray astronomy and cosmic ray physics. To… ▽ More The KM2A is the largest sub-array of the Large High Altitude Air Shower Observatory (LHAASO). It consists of 5216 electromagnetic particle detectors (EDs) and 1188 muon detectors (MDs). The data recorded by the EDs and MDs are used to reconstruct primary information of cosmic ray and gamma-ray showers. This information is used for physical analysis in gamma-ray astronomy and cosmic ray physics. To ensure the reliability of the LHAASO-KM2A data, a three-level quality control system has been established. It is used to monitor the status of detector units, stability of reconstructed parameters and the performance of the array based on observations of the Crab Nebula and Moon shadow. This paper will introduce the control system and its application on the LHAASO-KM2A data collected from August 2021 to July 2023. During this period, the pointing and angular resolution of the array were stable. From the observations of the Moon shadow and Crab Nebula, the results achieved using the two methods are consistent with each other. According to the observation of the Crab Nebula at energies from 25 TeV to 100 TeV, the time averaged pointing errors are estimated to be $-0.003^{\circ} \pm 0.005^{\circ}$ and $0.001^{\circ} \pm 0.006^{\circ}$ in the R.A. and Dec directions, respectively. △ Less

Submitted 13 June, 2024; v1 submitted 20 May, 2024; originally announced May 2024.

Comments: 15 pages, 9 figures

Energy dependence of the low-frequency quasi-periodic oscillations in Swift J1727.8-1613

Authors: Haifan Zhu, Wei Wang

Abstract: Based on observations from the Insight-Hard X-ray Modulation Telescope (Insight-HXMT), an analysis of Type-C quasi-periodic oscillations (QPOs) observed during the outburst of the new black hole candidate Swift J1727.8-1613 in 2023 was conducted. This analysis scrutinized the QPO's evolution throughout the outburst, particularly noting its rapid frequency escalation during two flare events. Utiliz… ▽ More Based on observations from the Insight-Hard X-ray Modulation Telescope (Insight-HXMT), an analysis of Type-C quasi-periodic oscillations (QPOs) observed during the outburst of the new black hole candidate Swift J1727.8-1613 in 2023 was conducted. This analysis scrutinized the QPO's evolution throughout the outburst, particularly noting its rapid frequency escalation during two flare events. Utilizing the energy range covered by Insight-HXMT, a dependency of the QPO frequency on energy was observed. Below approximately 3 Hz, minimal variations in frequency with energy were noted, whereas clear variations with photon energy were observed when it exceeded approximately 3 Hz. Additionally, a sharp drop in the rate of change was observed when the frequency exceeded approximately 8 Hz. This behavior, similar to several previously reported sources, suggests the presence of a common underlying physical mechanism. Moreover, the QPO rms-frequency relationship can be explained by the Lense-Thirring precession model. The relationship between rms-energy and phase lag with frequency suggests the black hole system as a high-inclination source. △ Less

Submitted 15 May, 2024; originally announced May 2024.

Comments: 12 pages, 9 figures, accept for the publication in ApJ

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

The Local Dark Matter Kinematic Substructure Based on LAMOST K Giants

Authors: Hai Zhu, Rui Guo, Juntai Shen, Jianglai Liu, Chao Liu, Xiang-Xiang Xue, Lan Zhang, Shude Mao

Abstract: Numerical simulations indicate that correlations exist between the velocity distributions of stars and dark matter (DM). We study the local DM velocity distribution based on these correlations. We select K giants from LAMOST DR8 cross-matched with Gaia DR3, which has robust measurements of three-dimensional velocity and metallicity, and separate them into the disk, halo substructure and main halo… ▽ More Numerical simulations indicate that correlations exist between the velocity distributions of stars and dark matter (DM). We study the local DM velocity distribution based on these correlations. We select K giants from LAMOST DR8 cross-matched with Gaia DR3, which has robust measurements of three-dimensional velocity and metallicity, and separate them into the disk, halo substructure and main halo components in the chemo-dynamical space utilizing the Gaussian Mixture Model. The substructure component is highly radially anisotropic, and possibly related to the Gaia-Enceladus-Sausage (GES) merger event, while the halo component is isotropic and accreted from the earliest mergers following the Maxwell-Boltzmann Distribution (Standard Halo Model, SHM). We find that the GES-like substructure contributes $\sim85\%$ of the local non-disk stars in the Solar neighbourhood, which is nearly invariant when applying different volume cuts or additional angular momentum constraints. Utilizing the metallicity-stellar-mass relation and the stellar-mass-halo-mass relation, we find that $\sim25_{-15}^{+24}\%$ of local DM is in the kinematic substructure. Combined with the stellar distributions of non-disk components, we compute the velocity distribution of local DM. The modified heliocentric velocity distribution of local DM shifts to a lower speed and has a sharper peak compared to the SHM, which yields updated detection limits for the DM direct detection experiments. Our work confirms that the local DM velocity distribution deviates from the SHM, and needs to be properly accounted in the DM detection experiments. △ Less

Submitted 30 April, 2024; originally announced April 2024.

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.

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)

On the Physical Nature of Ly$α$ Transmission Spikes in High Redshift Quasar Spectra

Authors: Hanjue Zhu, Nickolay Gnedin, Camille Avestruz

Abstract: We investigate Lyman-alpha (Ly$α$) transmission spikes at $5.2 < z < 6.8$ using synthetic quasar spectra from the ``Cosmic Reionization On Computers" simulations. We focus on understanding the relationship between these spikes and the properties of the intergalactic medium (IGM). Disentangling the complex interplay between IGM physics and the influence of galaxies on the generation of these spikes… ▽ More We investigate Lyman-alpha (Ly$α$) transmission spikes at $5.2 < z < 6.8$ using synthetic quasar spectra from the ``Cosmic Reionization On Computers" simulations. We focus on understanding the relationship between these spikes and the properties of the intergalactic medium (IGM). Disentangling the complex interplay between IGM physics and the influence of galaxies on the generation of these spikes presents a significant challenge. To address this, we employ Explainable Boosting Machines, an interpretable machine learning algorithm, to quantify the relative impact of various IGM properties on the Ly$α$ flux. Our findings reveal that gas density is the primary factor influencing absorption strength, followed by the intensity of background radiation and the temperature of the IGM. Ionizing radiation from local sources (i.e. galaxies) appears to have a minimal effect on Ly$α$ flux. The simulations show that transmission spikes predominantly occur in regions of low gas density. Our results challenge recent observational studies suggesting the origin of these spikes in regions with enhanced radiation. We demonstrate that Ly$α$ transmission spikes are largely a product of the large-scale structure, of which galaxies are biased tracers. △ Less

Submitted 9 January, 2024; originally announced January 2024.

Comments: 10 pages, 8 figures, comments are welcome

Nonlinear Effects In Black Hole Ringdown From Scattering Experiments I: spin and initial data dependence of quadratic mode coupling

Authors: Hengrui Zhu, Justin L. Ripley, Frans Pretorius, Sizheng Ma, Keefe Mitman, Robert Owen, Michael Boyle, Yitian Chen, Nils Deppe, Lawrence E. Kidder, Jordan Moxon, Kyle C. Nelli, Harald P. Pfeiffer, Mark A. Scheel, William Throwe, Nils L. Vu

Abstract: We investigate quadratic quasinormal mode coupling in black hole spacetime through numerical simulations of single perturbed black holes using both numerical relativity and second-order black hole perturbation theory. Focusing on the dominant $\ell=|m|=2$ quadrupolar modes, we find good agreement (within $\sim10\%$) between these approaches, with discrepancies attributed to truncation error and un… ▽ More We investigate quadratic quasinormal mode coupling in black hole spacetime through numerical simulations of single perturbed black holes using both numerical relativity and second-order black hole perturbation theory. Focusing on the dominant $\ell=|m|=2$ quadrupolar modes, we find good agreement (within $\sim10\%$) between these approaches, with discrepancies attributed to truncation error and uncertainties from mode fitting. Our results align with earlier studies extracting the coupling coefficients from select binary black hole merger simulations, showing consistency for the same remnant spins. Notably, the coupling coefficient is insensitive to a diverse range of initial data, including configurations that led to a significant (up to $5\%$) increase in the remnant black hole mass. These findings present opportunities for testing the nonlinear dynamics of general relativity with ground-based gravitational wave observatories. Lastly, we provide evidence of a bifurcation in coupling coefficients between counter-rotating and co-rotating quasinormal modes as black hole spin increases. △ Less

Submitted 1 January, 2024; originally announced January 2024.

Non-perturbative study of quantum many-body correlation effects in neutron stars: Equation of state

Authors: Hao-Fu Zhu, Xufen Wu, Guo-Zhu Liu

Abstract: Although neutron stars have been studied for decades, their internal structure remains enigmatic, mainly due to large uncertainties in the equation of state. In neutron stars, nucleons are strongly interacting by exchanging mesons, which leads to significant quantum many-body correlation effects. Previous mean-field calculations failed to capture these effects. Here, we develop a non-perturbative… ▽ More Although neutron stars have been studied for decades, their internal structure remains enigmatic, mainly due to large uncertainties in the equation of state. In neutron stars, nucleons are strongly interacting by exchanging mesons, which leads to significant quantum many-body correlation effects. Previous mean-field calculations failed to capture these effects. Here, we develop a non-perturbative quantum field-theoretic approach to handle strongly correlated nuclear matters. Based on a linear $σ$-$ω$-$ρ$ model with a density-dependent isovector coupling $Γ_ρ$, we reproduce all the six frequently-used nuclear matter observables by tuning six parameters. After including quantum many-body effects into the equation of state of realistic neutron-star matters, we obtain a mass-radius relation that is in good agreement with recent astrophysical observations. Our results show that the maximal neutron-star mass supported by neutron-dominated matters is $\sim 2.67M_{\odot}$, where $M_{\odot}$ is solar mass. △ Less

Submitted 15 April, 2024; v1 submitted 19 December, 2023; originally announced December 2023.

Comments: 5+15 pages, 8 figure, 2 tables, significantly revised version

Black Hole Spectroscopy for Precessing Binary Black Hole Coalescences

Authors: Hengrui Zhu, Harrison Siegel, Keefe Mitman, Maximiliano Isi, Will M. Farr, Michael Boyle, Nils Deppe, Lawrence E. Kidder, Sizheng Ma, Jordan Moxon, Kyle C. Nelli, Harald P. Pfeiffer, Mark A. Scheel, Saul A. Teukolsky, William Throwe, Vijay Varma, Nils L. Vu

Abstract: To accurately perform black hole spectroscopy, it is essential to know which quasinormal modes dominate astrophysical ringdown signals. In this Letter, we present a phenomenological description of the quasinormal modes that are excited in the ringdowns of comparable mass, quasi-circular precessing binary black hole coalescences. By analyzing an exhaustive catalog of numerical relativity simulation… ▽ More To accurately perform black hole spectroscopy, it is essential to know which quasinormal modes dominate astrophysical ringdown signals. In this Letter, we present a phenomenological description of the quasinormal modes that are excited in the ringdowns of comparable mass, quasi-circular precessing binary black hole coalescences. By analyzing an exhaustive catalog of numerical relativity simulations, we confirm that the relative fundamental quasinormal mode amplitudes of precessing systems are related to those of non-precessing systems by a simple rotation, and that additional structure in the spectrum is connected to the system's kick velocity and other asymmetries in the orbital dynamics. We find that the ringdowns of precessing systems need not be dominated by the ${(\ell,m)=(2,\pm 2)}$ quasinormal modes. These results build upon previous works on waveform modeling, and are consistent with a recent ringdown analysis of the LIGO-Virgo gravitational wave signal GW190521. △ Less

Submitted 13 December, 2023; originally announced December 2023.

Comments: Data Release and Analysis Scripts: https://github.com/HengruiPrinceton/precession_ringdown

Modulating binary dynamics via the termination of black hole superradiance

Authors: Kaiyuan Fan, Xi Tong, Yi Wang, Hui-Yu Zhu

Abstract: A superradiant cloud of ultralight bosons near a rotating black hole provides a smoking gun for particle physics in the infrared. However, tidal perturbations from a nearby binary companion can destabilise the boson cloud and even terminate superradiance. In this work, we consider the backreaction of superradiance termination to the dynamics of general binary orbits parametrised by their semi-latu… ▽ More A superradiant cloud of ultralight bosons near a rotating black hole provides a smoking gun for particle physics in the infrared. However, tidal perturbations from a nearby binary companion can destabilise the boson cloud and even terminate superradiance. In this work, we consider the backreaction of superradiance termination to the dynamics of general binary orbits parametrised by their semi-latus rectum, eccentricity and inclination angle. Our analysis focuses on Extreme Mass Ratio Inspiral (EMRI) systems and employs the period-average approximation to derive evolution equations of these binary parameters in the Newtonian limit. We find that the binary evolution history can be significantly modulated by the backreaction towards large circular equatorial orbits with reduced termination rate. This process can generically happen even away from the resonance bands. Our work therefore serves as a first step towards probing ultralight bosons through the statistics of EMRI binary parameters in the future. △ Less

Submitted 11 January, 2024; v1 submitted 28 November, 2023; originally announced November 2023.

Comments: 13 pages, 5 figures

Constraining Neutrino Cosmologies with Nonlinear Reconstruction

Authors: Shi-Hui Zang, Hong-Ming Zhu

Abstract: Nonlinear gravitational evolution induces strong nonlinearities in the observed cosmological density fields, leading to positive off-diagonal correlations in the power spectrum covariance. This has caused the information saturation in the power spectrum, e.g., the neutrino mass constraints from the nonlinear power spectra are lower than their linear counterparts by a factor of $\sim2$ at $z=0$. In… ▽ More Nonlinear gravitational evolution induces strong nonlinearities in the observed cosmological density fields, leading to positive off-diagonal correlations in the power spectrum covariance. This has caused the information saturation in the power spectrum, e.g., the neutrino mass constraints from the nonlinear power spectra are lower than their linear counterparts by a factor of $\sim2$ at $z=0$. In this paper, we explore how nonlinear reconstruction methods improve the cosmological information from nonlinear cosmic fields. By applying nonlinear reconstruction to cold dark matter fields from the Quijote simulations, we find that nonlinear reconstruction can improve the constraints on cosmological parameters significantly, nearly reaching the linear theory limit. For neutrino mass, the result is only $12\%$ lower than the linear power spectrum, i.e., the theoretical best result. This makes nonlinear reconstruction an efficient and useful method to extract neutrino information from current and upcoming galaxy surveys. △ Less

Submitted 27 November, 2023; originally announced November 2023.

Comments: 14 pages, 8 figures, accepted version

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

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

Single and coupled cavity mode sensing schemes using a diagnostic field

Authors: Aaron W. Goodwin-Jones, Haochen Zhu, Carl Blair, Daniel D. Brown, Joris van Heijningen, Li Ju, Chunnong Zhao

Abstract: Precise optical mode matching is of critical importance in experiments using squeezed-vacuum states. Automatic spatial-mode matching schemes have the potential to reduce losses and improve loss stability. However, in quantum-enhanced coupled-cavity experiments, such as gravitational-wave detectors, one must also ensure that the sub-cavities are also mode matched. We propose a new mode sensing sche… ▽ More Precise optical mode matching is of critical importance in experiments using squeezed-vacuum states. Automatic spatial-mode matching schemes have the potential to reduce losses and improve loss stability. However, in quantum-enhanced coupled-cavity experiments, such as gravitational-wave detectors, one must also ensure that the sub-cavities are also mode matched. We propose a new mode sensing scheme, which works for simple and coupled cavities. The scheme requires no moving parts, nor tuning of Gouy phases. Instead a diagnostic field tuned to the HG20/LG10 mode frequency is used. The error signals are derived to be proportional to the difference in waist position, and difference in Rayleigh ranges, between the sub-cavity eigenmodes. The two error signals are separable by 90 degrees of demodulation phase. We demonstrate reasonable error signals for a simplified Einstein Telescope optical design. This work will facilitate routine use of extremely high levels of squeezing in current and future gravitational-wave detectors. △ Less

Submitted 29 August, 2023; originally announced August 2023.

Report number: LIGO-P2300010

A Joint Fermi-GBM and Swift-BAT Analysis of Gravitational-Wave Candidates from the Third Gravitational-wave Observing Run

Authors: C. Fletcher, J. Wood, R. Hamburg, P. Veres, C. M. Hui, E. Bissaldi, M. S. Briggs, E. Burns, W. H. Cleveland, M. M. Giles, A. Goldstein, B. A. Hristov, D. Kocevski, S. Lesage, B. Mailyan, C. Malacaria, S. Poolakkil, A. von Kienlin, C. A. Wilson-Hodge, The Fermi Gamma-ray Burst Monitor Team, M. Crnogorčević, J. DeLaunay, A. Tohuvavohu, R. Caputo, S. B. Cenko , et al. (1674 additional authors not shown)

Abstract: We present Fermi Gamma-ray Burst Monitor (Fermi-GBM) and Swift Burst Alert Telescope (Swift-BAT) searches for gamma-ray/X-ray counterparts to gravitational wave (GW) candidate events identified during the third observing run of the Advanced LIGO and Advanced Virgo detectors. Using Fermi-GBM on-board triggers and sub-threshold gamma-ray burst (GRB) candidates found in the Fermi-GBM ground analyses,… ▽ More We present Fermi Gamma-ray Burst Monitor (Fermi-GBM) and Swift Burst Alert Telescope (Swift-BAT) searches for gamma-ray/X-ray counterparts to gravitational wave (GW) candidate events identified during the third observing run of the Advanced LIGO and Advanced Virgo detectors. Using Fermi-GBM on-board triggers and sub-threshold gamma-ray burst (GRB) candidates found in the Fermi-GBM ground analyses, the Targeted Search and the Untargeted Search, we investigate whether there are any coincident GRBs associated with the GWs. We also search the Swift-BAT rate data around the GW times to determine whether a GRB counterpart is present. No counterparts are found. Using both the Fermi-GBM Targeted Search and the Swift-BAT search, we calculate flux upper limits and present joint upper limits on the gamma-ray luminosity of each GW. Given these limits, we constrain theoretical models for the emission of gamma-rays from binary black hole mergers. △ Less

Submitted 25 August, 2023; originally announced August 2023.

Timing Analysis of a New Black Hole Candidate MAXI J1803-298 with Insight-HXMT and NICER

Authors: Haifan Zhu, Xiao Chen, Wei Wang

Abstract: We present the timing analysis results of MAXI J1803$-$298, a black hole candidate, during its 2021 outburst using data obtained from the \textit{Insight-}HXMT and \textit{NICER} telescopes. Our analysis reveals that the source undergoes a state transition from the low/hard state to the hard intermediate state, followed by the soft intermediate state, and ultimately reaching the high/soft state. W… ▽ More We present the timing analysis results of MAXI J1803$-$298, a black hole candidate, during its 2021 outburst using data obtained from the \textit{Insight-}HXMT and \textit{NICER} telescopes. Our analysis reveals that the source undergoes a state transition from the low/hard state to the hard intermediate state, followed by the soft intermediate state, and ultimately reaching the high/soft state. We searched for the quasi-periodic oscillations (QPOs) and studied the characteristics of the outburst. At the beginning of the outburst, the source was in the hard state, many type-C QPOs were seen in the \textit{Insight-}HXMT data, and the frequency of these QPOs increased from $\sim 0.16$ to $2.6$ Hz. Our analysis of the type-C QPOs' rms-frequency relationship indicates a turning point in the frequency. We also analyzed the phase lag versus frequency and energy relationship and deduced that the source likely has a high inclination angle, consistent with previous research. The observed rms and phase lag features in type-C QPOs could be explained by the Lense-Thirring precession model, while the alternatives would be still viable. The lag spectrum of type-B QPO exhibits a "\textit{U-shaped}" pattern similar to many other sources, and the type-B QPOs' rms increase as the energy rises. This phenomenon can be explained by the dual-corona model. △ Less

Submitted 24 June, 2023; originally announced June 2023.

Comments: 11 pages, reference: 2023, MNRAS, 523, 4394-4404

Spatially Resolved Kinematics of Extraplanar Diffuse Ionized Gas in NGC$\,$3511 and NGC$\,$3513

Authors: Hanjue Zhu, Erin Boettcher, Hsiao-Wen Chen

Abstract: Gaseous, disk-halo interfaces are shaped by physical processes that are critical to disk galaxy evolution, including gas accretion and galactic outflows. However, observations indicate that extraplanar diffuse ionized gas (eDIG) layers have scale heights several times higher than their thermal scale heights. This discrepancy poses a challenge to our current understanding of the disk-halo interface… ▽ More Gaseous, disk-halo interfaces are shaped by physical processes that are critical to disk galaxy evolution, including gas accretion and galactic outflows. However, observations indicate that extraplanar diffuse ionized gas (eDIG) layers have scale heights several times higher than their thermal scale heights. This discrepancy poses a challenge to our current understanding of the disk-halo interface. In this paper, we present a spatially-resolved case study of the eDIG layers in a nearby pair of sub-$L_*$ disk galaxies NGC$\,$3511/3513 using long-slit spectroscopy. We decompose optical nebular lines from the warm interstellar medium and disk-halo interfaces into narrow and broad velocity components. We show that in NGC$\,$3511, the broad component has three distinctive characteristics in comparison to the narrow component: (1) elevated [NII]$λ6583/Hα$ and [SII]$λ6716/Hα$ line ratios, (2) significantly higher velocity dispersions (a median $\langleσ\rangle_{\text{Broad}} = 24\,$km/s compared to $\langleσ\rangle_{\text{Narrow}} = 13\,$km/s), and (3) a rotational velocity lag. Together, these characteristics support an origin in an extraplanar, thick gaseous disk. In NGC$\,$3513, the broad component is consistent with localized outflows making their way out of the galactic disk. Our findings demonstrate that gas circulation at the disk-halo interface is present in both galaxies. Additionally, we test a dynamic equilibrium model with pressure support supplied by thermal and turbulent motions. Based on measurements of the eDIG velocity dispersion in NGC$\,$3511, we demonstrate that turbulent motions increase the scale height by at least a factor of a few above the thermal scale height, with $h_{z} \gtrsim 0.2 - 0.4$ kpc at $R = 3 - 5$ kpc. This highlights the importance of turbulent motions to the vertical structure of the gaseous, disk-halo interface. △ Less

Submitted 22 June, 2023; originally announced June 2023.

Comments: 14 pages, 14 figures, comments are welcome

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

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)

A moderate spin for the black hole in X-ray binary MAXI J1348-630 revealed by Insight-HXMT

Authors: Hanji Wu, Wei Wang, Na Sai, Haifan Zhu, Jiashi Chen

Abstract: MAXI J1348-630 is a low-mass X-ray black hole binary located in the Galaxy and undergone the X-ray outburst in 2019. We analyzed the observation data in very soft state during the outburst between MJD 58588 and MJD 58596 based on the Insight-HXMT observations from 2 -- 20 keV via the continuum fitting method to measure the spin of the stellar-mass black hole in MAXI J1348-630. The inner disk tempe… ▽ More MAXI J1348-630 is a low-mass X-ray black hole binary located in the Galaxy and undergone the X-ray outburst in 2019. We analyzed the observation data in very soft state during the outburst between MJD 58588 and MJD 58596 based on the Insight-HXMT observations from 2 -- 20 keV via the continuum fitting method to measure the spin of the stellar-mass black hole in MAXI J1348-630. The inner disk temperature and the apparent inner disk radius were found to be $0.47\pm 0.01 \rm keV$ and $5.33\pm 0.10 \ R_{g}$ from the observation data modeled by the multicolor disc blackbody model. Assuming the distance of the source $D\sim 3.4 \rm kpc$, the mass of the black hole $M\sim 11 \ M_{\odot}$, and the inclination of the system $i\sim 29.2^{\circ}$, the spin is determined to be $a_{\star}=0.41\pm 0.03$ for fixing hardening factor at 1.6 and $n_{H}=8.6\times 10^{21} \rm cm^{-2}$. Besides, considering the uncertainty of the parameters $D, M, i$ of this system, with the Monte Carlo analysis, we still confirm the moderate spin of the black hole as $a_{\star}=0.42^{+0.13}_{-0.50}$. Some spectral parameters (e.g., column density and hardening factor) which could affect the measurements of the BH spin are also briefly discussed. △ Less

Submitted 26 April, 2023; originally announced April 2023.

Comments: 10 pages, 14 figures, 5 tables, accept for publication in MNRAS

Cosmic Reionization On Computers: Statistics, Physical Properties and Environment of Lyman Limit Systems at $z\sim6$

Authors: Jiawen Fan, Hanjue Zhu, Camille Avestruz, Nickolay Y. Gnedin

Abstract: Lyman limit systems (LLSs) are dense hydrogen clouds with high enough HI column densities to absorb Lyman continuum photons emitted from distant quasars. Their high column densities imply an origin in dense environments; however, the statistics and distribution of LLSs at high redshifts still remain uncertain. In this paper, we use self-consistent radiative transfer cosmological simulations from t… ▽ More Lyman limit systems (LLSs) are dense hydrogen clouds with high enough HI column densities to absorb Lyman continuum photons emitted from distant quasars. Their high column densities imply an origin in dense environments; however, the statistics and distribution of LLSs at high redshifts still remain uncertain. In this paper, we use self-consistent radiative transfer cosmological simulations from the "Cosmic Reionization On Computers" (CROC) project to study the physical properties of LLSs at the tail end of cosmic reionization at $z\sim6$. We generate 3000 synthetic quasar sightlines to obtain a large number of LLS samples in the simulations. In addition, with the high physical fidelity and resolution of CROC, we are able to quantify the association between these LLS samples and nearby galaxies. Our results show that the fraction LLSs spatially associated with nearby galaxies is increasing with the HI column density. Moreover, we find that LLSs that are not near any galaxy typically reside in filamentary structures connecting neighboring galaxies in the intergalactic medium (IGM). This quantification of the distribution and associations of LLSs to large scale structures informs our understanding of the IGM-galaxy connection during the Epoch of Reionization, and provides a theoretical basis for interpreting future observations. △ Less

Submitted 1 March, 2024; v1 submitted 14 December, 2022; originally announced December 2022.

Comments: 11 pages, 10 figures, submitted to ApJ

Journal ref: 2024 ApJ 963 45

Cosmic tidal reconstruction in redshift space

Authors: Shi-Hui Zang, Hong-Ming Zhu, Marcel Schmittfull, Ue-Li Pen

Abstract: The gravitational coupling between large- and small-scale density perturbations leads to anisotropic distortions to the local small-scale matter fluctuations. Such local anisotropic distortions can be used to reconstruct the large-scale matter distribution, known as tidal reconstruction. In this paper, we apply the tidal reconstruction methods to simulated galaxies in redshift space. We find that… ▽ More The gravitational coupling between large- and small-scale density perturbations leads to anisotropic distortions to the local small-scale matter fluctuations. Such local anisotropic distortions can be used to reconstruct the large-scale matter distribution, known as tidal reconstruction. In this paper, we apply the tidal reconstruction methods to simulated galaxies in redshift space. We find that redshift-space distortions lead to anisotropic reconstruction results. While the reconstructed radial modes are more noisy mainly due to the small-scale velocity dispersion, the transverse modes are still reconstructed with high fidelity, well correlated with the original large-scale density modes. The bias of the reconstructed field at large scales shows a simple angular dependence which can be described by a form similar to the linear redshift-space distortion. The noise power spectrum is nearly isotropic and scale-independent on large scales. This makes the reconstructed tides field an ideal tracer for cosmic variance cancellation and multi-tracer analysis and has profound implications for future 21~cm intensity mapping surveys. △ Less

Submitted 8 December, 2022; originally announced December 2022.

Comments: 18 pages, 11 figures

Magneto hydrodynamic simulations of the supernova remnant G1.9+0.3

Authors: Shaobo Zhang, Wenwu Tian, Mengfei Zhang, Hui Zhu, Xiaohong Cui

Abstract: The youngest Galactic supernova remnant G1.9+0.3 shows a discrete feature between its radio and X-ray morphologies. The observed radio morphology features a single maximum in the north, while the X-ray observation shows two opposite 'ears' on the east and west sides. Using 3D magneto hydrodynamical simulations, we investigate the formation of the discrete feature of the remnant. We have tested dif… ▽ More The youngest Galactic supernova remnant G1.9+0.3 shows a discrete feature between its radio and X-ray morphologies. The observed radio morphology features a single maximum in the north, while the X-ray observation shows two opposite 'ears' on the east and west sides. Using 3D magneto hydrodynamical simulations, we investigate the formation of the discrete feature of the remnant. We have tested different parameters for better simulation and reproduced similar discrete features under an environment with density gradient and an environment with clump, which provides a possible explanation of the observation. △ Less

Submitted 19 January, 2023; v1 submitted 20 November, 2022; originally announced November 2022.

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

Journal ref: The Astrophysical Journal, 942:94 (5pp),(2023)

Adaptive Analytical Ray Tracing of Black Hole Photon Rings

Authors: Alejandro Cárdenas-Avendaño, Alexandru Lupsasca, Hengrui Zhu

Abstract: Recent interferometric observations by the Event Horizon Telescope have resolved the horizon-scale emission from sources in the vicinity of nearby supermassive black holes. Future space-based interferometers promise to measure the "photon ring"--a narrow, ring-shaped, lensed feature predicted by general relativity, but not yet observed--and thereby open a new window into strong gravity. Here we pr… ▽ More Recent interferometric observations by the Event Horizon Telescope have resolved the horizon-scale emission from sources in the vicinity of nearby supermassive black holes. Future space-based interferometers promise to measure the "photon ring"--a narrow, ring-shaped, lensed feature predicted by general relativity, but not yet observed--and thereby open a new window into strong gravity. Here we present AART: an Adaptive Analytical Ray-Tracing code that exploits the integrability of light propagation in the Kerr spacetime to rapidly compute high-resolution simulated black hole images, together with the corresponding radio visibility accessible on very long space-ground baselines. The code samples images on a nonuniform adaptive grid that is specially tailored to the lensing behavior of the Kerr geometry and is therefore particularly well-suited to studying photon rings. This numerical approach guarantees that interferometric signatures are correctly computed on long baselines, and the modularity of the code allows for detailed studies of equatorial sources with complex emission profiles and time variability. To demonstrate its capabilities, we use AART to simulate a black hole movie of a stochastic, non-stationary, non-axisymmetric equatorial source; by time-averaging the visibility amplitude of each snapshot, we are able to extract the projected diameter of the photon ring and recover the shape predicted by general relativity. △ Less

Submitted 4 June, 2024; v1 submitted 14 November, 2022; originally announced November 2022.

Comments: 39 pages, 16 figures. Updated to match published version. V4: corrected two very minor typos (in Eqs. 50 and 70)

Journal ref: Phys. Rev. D 107, 043030, 2023

A Study of Radio Knots within Supernova Remnant Cassiopeia A

Authors: Xianhuan Lei, Hui Zhu, Haiyan Zhang, Wenwu Tian

Abstract: The study on the dynamic evolution of young supernova remnants (SNRs) is an important way to understand the density structure of the progenitor's circumstellar medium. We have reported the acceleration or deceleration, proper motion and brightness changes of 260 compact radio features in the second youngest known SNR Cas A at 5\,GHz based on the VLA data of five epochs from 1987 to 2004. The radio… ▽ More The study on the dynamic evolution of young supernova remnants (SNRs) is an important way to understand the density structure of the progenitor's circumstellar medium. We have reported the acceleration or deceleration, proper motion and brightness changes of 260 compact radio features in the second youngest known SNR Cas A at 5\,GHz based on the VLA data of five epochs from 1987 to 2004. The radio expansion center locates at $α(1950)=23^{\rm h}21^{\rm m}9^{\rm s}_{\cdot}7 \pm 0^{\rm s}_{\cdot}29, δ(1950)=+58^{\circ}32^{\prime}25^{\prime\prime}_{\cdot}2 \pm 2^{\prime\prime}_{\cdot}2$. Three-quarters of the compact knots are decelerating, this suggests that there are significant density fluctuations in the stellar winds of the remnant's progenitor. We have verified that the acceleration or deceleration of compact knots are not related with the distribution of brightness. The brightening, fading, disappearing or new appearing of compact radio features in Cas A suggests that the magnetic field in the remnant is changing rapidly. △ Less

Submitted 9 November, 2022; originally announced November 2022.

Comments: 20 pages,9 figures, ApJ accepted

Missing one-loop contributions in secondary gravitational waves

Authors: Chao Chen, Atsuhisa Ota, Hui-Yu Zhu, Yuhang Zhu

Abstract: We find several missing one-loop-order contributions in previous considerations about secondary gravitational waves induced at nonlinear order in cosmological perturbations. We consider a consistent perturbative expansion to third-order in cosmological perturbations, including higher-order interactions and iterative solutions ignored in the previous literature. Tensor fluctuations induced by the s… ▽ More We find several missing one-loop-order contributions in previous considerations about secondary gravitational waves induced at nonlinear order in cosmological perturbations. We consider a consistent perturbative expansion to third-order in cosmological perturbations, including higher-order interactions and iterative solutions ignored in the previous literature. Tensor fluctuations induced by the source with two scalar and one tensor perturbations are correlated with the first-order tensor fluctuation and thus give a one-loop-order correction to the tensor power spectrum. The missing loop correction is \textit{scale-invariant} and \textit{negative} in the superhorion region, which secondarily reduces the initial primordial tensor power spectrum prior to the horizon re-entry. Such an IR behavior is very different from the auto-spectrum of second-order induced tensor modes discussed in the previous literature and can be important for the actual gravitational wave measurements. For a sharp peak of scalar fluctuations with $A_ζ=10^{-2}$ at $k_*=10^{5}h/{\rm Mpc}$ motivated by the LIGO/Virgo events, we show that the tensor power spectrum at the cosmic microwave background scale reduces by at most 35%. Hence, the polarization B-mode might not be seen because of the reduction of the original tensor spectrum due to the secondary effect of primordial black hole formation. △ Less

Submitted 31 October, 2022; originally announced October 2022.

Comments: 13 pages, 6 figures

Mass-Metallicity Relation during the Epoch of Reionization in the CROC Simulations

Authors: Isaac Noel, Hanjue Zhu, Nickolay Gnedin

Abstract: The low-redshift mass-metallicity relation (MZR) is well studied, but the high-redshift MZR remains difficult to observe. To study the early MZR further, we analyze the Cosmic Reionization on Computers (CROC) simulations with a focus on the MZR from redshifts 5 to 10. We find that, across all redshifts, CROC galaxies exhibit similar stellar-phase and gas-phase MZRs that flatten with higher stellar… ▽ More The low-redshift mass-metallicity relation (MZR) is well studied, but the high-redshift MZR remains difficult to observe. To study the early MZR further, we analyze the Cosmic Reionization on Computers (CROC) simulations with a focus on the MZR from redshifts 5 to 10. We find that, across all redshifts, CROC galaxies exhibit similar stellar-phase and gas-phase MZRs that flatten with higher stellar mass. We attribute this flattening to the inaccurate star formation and feedback modeling in CROC (star formation is overly suppressed in massive CROC galaxies). In addition, we show that the ratio between stellar metallicity and gas metallicity ($Z_*/Z_{gas}$) decreases as stellar age increases, meaning that in CROC galaxies, gas accretion rate is lower than metal production rate. With JWST we will be able to compare our predictions to observations of the Epoch of Reionization and understand better early galaxy formation. △ Less

Submitted 30 October, 2022; originally announced October 2022.

Comments: 5 pages, 5 figures, submitted to ApJ

Galaxy Spin Classification I: Z-wise vs S-wise Spirals With Chirality Equivariant Residual Network

Authors: He Jia, Hong-Ming Zhu, Ue-Li Pen

Abstract: The angular momentum of galaxies (galaxy spin) contains rich information about the initial condition of the Universe, yet it is challenging to efficiently measure the spin direction for the tremendous amount of galaxies that are being mapped by the ongoing and forthcoming cosmological surveys. We present a machine learning based classifier for the Z-wise vs S-wise spirals, which can help to break… ▽ More The angular momentum of galaxies (galaxy spin) contains rich information about the initial condition of the Universe, yet it is challenging to efficiently measure the spin direction for the tremendous amount of galaxies that are being mapped by the ongoing and forthcoming cosmological surveys. We present a machine learning based classifier for the Z-wise vs S-wise spirals, which can help to break the degeneracy in the galaxy spin direction measurement. The proposed Chirality Equivariant Residual Network (CE-ResNet) is manifestly equivariant under a reflection of the input image, which guarantees that there is no inherent asymmetry between the Z-wise and S-wise probability estimators. We train the model with Sloan Digital Sky Survey (SDSS) images, with the training labels given by the Galaxy Zoo 1 (GZ1) project. A combination of data augmentation tricks are used during the training, making the model more robust to be applied to other surveys. We find a $\sim\!30\%$ increase of both types of spirals when Dark Energy Spectroscopic Instrument (DESI) images are used for classification, due to the better imaging quality of DESI. We verify that the $\sim\!7σ$ difference between the numbers of Z-wise and S-wise spirals is due to human bias, since the discrepancy drops to $<\!1.8σ$ with our CE-ResNet classification results. We discuss the potential systematics that are relevant to the future cosmological applications. △ Less

Submitted 6 December, 2022; v1 submitted 9 October, 2022; originally announced October 2022.

Comments: 13+4 pages, 11 figures, 2 tables, accepted by ApJ

Cosmic Reionization On Computers: Baryonic Effects on Halo Concentrations During the Epoch of Reionization

Authors: Hanjue Zhu, Nickolay Y. Gnedin

Abstract: Baryons both increase halo concentration through adiabatic contraction and expel mass through feedback processes. However, it is not well understood how the radiation fields prevalent during the epoch of reionization affect the evolution of concentration in dark matter halos. We investigate how baryonic physics during the epoch of reionization modify the structure of dark matter halos in the Cosmi… ▽ More Baryons both increase halo concentration through adiabatic contraction and expel mass through feedback processes. However, it is not well understood how the radiation fields prevalent during the epoch of reionization affect the evolution of concentration in dark matter halos. We investigate how baryonic physics during the epoch of reionization modify the structure of dark matter halos in the Cosmic Reionization On Computers (CROC) simulations. We use two different measures of halo concentration to quantify the effects. We compare concentrations of halos matched between full physics simulations and dark-matter-only simulations with identical initial conditions between $5 \leq z \leq 9$. Baryons in full physics simulations do pull matter towards the center, increasing the maximum circular velocity compared to dark-matter-only simulations. However, their overall effects are much less than if all the baryons were simply centrally concentrated indicating that heating processes efficiently counteract cooling effects. Finally, we show that the baryonic effects on halo concentrations at $z\approx5$ are relatively insensitive to environmental variations of reionization history. These results are pertinent to models of galaxy-halo connection during the epoch of reionization. △ Less

Submitted 29 September, 2022; originally announced September 2022.

Comments: 9 pages, 7 figures, comments are welcome

Complementary astrometry of Cassini Imaging Science Subsystem Images of Phoebe

Authors: Q. F. Zhang, W. H. Qin, Y. L. Ma, V. Lainey, N. J. Cooper, N. Rambaux, Y. Li, W. H. Zhu

Abstract: Phoebe is the only major satellite of Saturn with a retrograde orbit. The Cassini Imaging Science Subsystem (ISS) took a lot of Phoebe images between 2004 and 2017, but only a selection of them has been reduced. In this paper, we reduced the remaining ISS images of Phoebe. In the reduction, the Gaia EDR3 catalogue was used to provide the reference stars' positions, and the modified moment was used… ▽ More Phoebe is the only major satellite of Saturn with a retrograde orbit. The Cassini Imaging Science Subsystem (ISS) took a lot of Phoebe images between 2004 and 2017, but only a selection of them has been reduced. In this paper, we reduced the remaining ISS images of Phoebe. In the reduction, the Gaia EDR3 catalogue was used to provide the reference stars' positions, and the modified moment was used to measure the centre of image stars and Phoebe. Finally, a total of 834 ISS images of Phoebe have been reduced successfully. Compared with the JPL ephemeris SAT375, Phoebe's positions are consistent. The average residuals in the right ascension and declination are 0.08" and $-0.05$", and the standard deviations of the residuals are about 0.2". In terms of residuals in linear units, the means in the right ascension and declination are about 5 km and $-2$ km, respectively; The standard deviations are about 11 km. Compared with the JPL ephemeris SAT427 and IMCCE ephemeris PH20, our measurements show a strong bias and a large dispersion. △ Less

Submitted 22 August, 2022; originally announced August 2022.

Comments: 8 Pages, 5 figures, Journal article

Pulsar Double-lensing Sheds Light on the Origin of Extreme Scattering Events

Authors: Hengrui Zhu, Daniel Baker, Ue-Li Pen, Dan R. Stinebring, Marten H. van Kerkwijk

Abstract: In extreme scattering events, the brightness of a compact radio source drops significantly, as light is refracted out of the line of sight by foreground plasma lenses. Despite recent efforts, the nature of these lenses has remained a puzzle, because any roughly round lens would be so highly overpressurized relative to the interstellar medium that it could only exist for about a year. This, combine… ▽ More In extreme scattering events, the brightness of a compact radio source drops significantly, as light is refracted out of the line of sight by foreground plasma lenses. Despite recent efforts, the nature of these lenses has remained a puzzle, because any roughly round lens would be so highly overpressurized relative to the interstellar medium that it could only exist for about a year. This, combined with a lack of constraints on distances and velocities, has led to a plethora of theoretical models. We present observations of a dramatic double-lensing event in pulsar PSR~B0834+06 and use a novel phase-retrieval technique to show that the data can be reproduced remarkably well with a two-screen model: one screen with many small lenses and another with a single, strong one. We further show that the latter lens is so strong that it would inevitably cause extreme scattering events. Our observations show that the lens moves slowly and is highly elongated on the sky. If similarly elongated along the line of sight, as would arise naturally from a sheet of plasma viewed nearly edge-on, no large over-pressure is required and hence the lens could be long-lived. Our new technique opens up the possibility of probing interstellar plasma structures in detail, leading to understanding crucial for high-precision pulsar timing and the subsequent detection of gravitational waves. △ Less

Submitted 29 December, 2022; v1 submitted 14 August, 2022; originally announced August 2022.

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)

A multi-cubic-kilometre neutrino telescope in the western Pacific Ocean

Authors: Z. P. Ye, F. Hu, W. Tian, Q. C. Chang, Y. L. Chang, Z. S. Cheng, J. Gao, T. Ge, G. H. Gong, J. Guo, X. X. Guo, X. G. He, J. T. Huang, K. Jiang, P. K. Jiang, Y. P. Jing, H. L. Li, J. L. Li, L. Li, W. L. Li, Z. Li, N. Y. Liao, Q. Lin, F. Liu, J. L. Liu , et al. (33 additional authors not shown)

Abstract: Next-generation neutrino telescopes with significantly improved sensitivity are required to pinpoint the sources of the diffuse astrophysical neutrino flux detected by IceCube and uncover the century-old puzzle of cosmic ray origins. A detector near the equator will provide a unique viewpoint of the neutrino sky, complementing IceCube and other neutrino telescopes in the Northern Hemisphere. Here… ▽ More Next-generation neutrino telescopes with significantly improved sensitivity are required to pinpoint the sources of the diffuse astrophysical neutrino flux detected by IceCube and uncover the century-old puzzle of cosmic ray origins. A detector near the equator will provide a unique viewpoint of the neutrino sky, complementing IceCube and other neutrino telescopes in the Northern Hemisphere. Here we present results from an expedition to the north-eastern region of the South China Sea, in the western Pacific Ocean. A favorable neutrino telescope site was found on an abyssal plain at a depth of $\sim$ 3.5km. At depths below 3km, the sea current speed, water absorption and scattering lengths for Cherenkov light, were measured to be $v_{\mathrm{c}}<$10cm/s, $λ_{\mathrm{abs} }\simeq$ 27m and $λ_{\mathrm{sca} }\simeq$ 63m, respectively. Accounting for these measurements, we present the design and expected performance of a next-generation neutrino telescope, TRopIcal DEep-sea Neutrino Telescope (TRIDENT). With its advanced photon-detection technology and large dimensions, TRIDENT expects to observe the IceCube steady source candidate NGC 1068 with 5$σ$ significance within 1 year of operation. This level of sensitivity will open a new arena for diagnosing the origin of cosmic rays and probing fundamental physics over astronomical baselines. △ Less

Submitted 13 May, 2024; v1 submitted 10 July, 2022; originally announced July 2022.

Comments: 34 pages,12 figures. Correspondence should be addressed to D. L. Xu: donglianxu@sjtu.edu.cn

BURSTT: Bustling Universe Radio Survey Telescope in Taiwan

Authors: Hsiu-Hsien Lin, Kai-yang Lin, Chao-Te Li, Yao-Huan Tseng, Homin Jiang, Jen-Hung Wang, Jen-Chieh Cheng, Ue-Li Pen, Ming-Tang Chen, Pisin Chen, Yaocheng Chen, Tomotsugu Goto, Tetsuya Hashimoto, Yuh-Jing Hwang, Sun-Kun King, Derek Kubo, Chung-Yun Kuo, Adam Mills, Jiwoo Nam, Peter Oshiro, Chang-Shao Shen, Hsien-Chun Tseng, Shih-Hao Wang, Vigo Feng-Shun Wu, Geoffrey Bower , et al. (22 additional authors not shown)

Abstract: Fast Radio Bursts (FRBs) are bright millisecond-duration radio transients that appear about 1,000 times per day, all-sky, for a fluence threshold 5 Jy ms at 600 MHz. The FRB radio-emission physics and the compact objects involved in these events are subjects of intense active debate. To better constrain source models, the Bustling Universe Radio Survey Telescope in Taiwan (BURSTT) is optimized to… ▽ More Fast Radio Bursts (FRBs) are bright millisecond-duration radio transients that appear about 1,000 times per day, all-sky, for a fluence threshold 5 Jy ms at 600 MHz. The FRB radio-emission physics and the compact objects involved in these events are subjects of intense active debate. To better constrain source models, the Bustling Universe Radio Survey Telescope in Taiwan (BURSTT) is optimized to discover and localize a large sample of rare, high-fluence, nearby FRBs. This is the population most amenable to multi-messenger, multi-wavelength follow-up, allowing deeper understanding of source mechanisms. BURSTT will provide horizon-to-horizon sky coverage with a half power field-of-view (FoV) of $\sim$10$^{4}$ deg$^{2}$, a 400 MHz effective bandwidth between 300-800 MHz, and sub-arcsecond localization, made possible using outrigger stations hundreds to thousands of km from the main array. Initially, BURSTT will employ 256 antennas. After tests of various antenna designs and optimization of system performance we plan to expand to 2048 antennas. We estimate that BURSTT-256 will detect and localize $\sim$100 bright ($\geq$100 Jy ms) FRBs per year. Another advantage of BURSTT's large FoV and continuous operation will be greatly enhanced monitoring of FRBs for repetition. The current lack of sensitive all-sky observations likely means that many repeating FRBs are currently cataloged as single-event FRBs. △ Less

Submitted 26 September, 2022; v1 submitted 17 June, 2022; originally announced June 2022.

Comments: Published version

Journal ref: PASP 134 094106 (2022)

LAMOST MRS-N Observations of the W80 Region

Authors: Yao Li, Chao-Jian Wu, Yong-Qiang Yao, Wei Zhang, Jia Yin, Juan-Juan Ren, Chih-Hao Hsia, Rui Zhuang, Jian-Jun Chen, Yu-Zhong Wu, Hui Zhu, Bin Li, Yong-Hui Hou, Meng-Yuan Yao, Hong Wu

Abstract: The spectral observations and analysis for the W80 Region are presented by using the data of Medium-Resolution Spectroscopic Survey of Nebulae (MRS-N) with the Large Sky Area Multi-Object Fiber Spectroscopy Telescope (LAMOST). A total of 2982 high-quality nebular spectra have been obtained in the 20 square degree field of view (FoV) which covers the W80 complex, and the largest sample of spectral… ▽ More The spectral observations and analysis for the W80 Region are presented by using the data of Medium-Resolution Spectroscopic Survey of Nebulae (MRS-N) with the Large Sky Area Multi-Object Fiber Spectroscopy Telescope (LAMOST). A total of 2982 high-quality nebular spectra have been obtained in the 20 square degree field of view (FoV) which covers the W80 complex, and the largest sample of spectral data have been established for the first time. The relative intensities, radial velocities (RVs), and Full Widths at Half Maximum (FWHMs) are measured with the high spectral resolution of LAMOST MRS, for H$α$ $λ$ 6563 Å, [\ion{N}{ii}] $λ$$λ$ 6548 Å, 6584 Å\ , and [\ion{S}{ii}] $λ$$λ$ 6716 Å, 6731 Å\ emission lines. In the field of view of whole W80 Region, the strongest line emissions are found to be consistent with the bright nebulae, NGC 7000, IC 5070, and LBN 391, and weak line emissions also truly exist in the Middle Region, where no bright nebulae are detected by the wide-band optical observations. The large-scale spectral observations to the W80 Region reveal the systematic spatial variations of RVs and FWHMs, and several unique structural features. A 'curved feature' to the east of the NGC 7000, and a 'jet feature' to the west of the LBN 391 are detected to be showing with larger radial velocities. A 'wider FWHM region' is identified in the eastern part of the NGC 7000. The variations of [\ion{S}{ii}] / H$α$ ratios display a gradient from southwest to northeast in the NGC 7000 region, and manifest a ring shape around the 'W80 bubble' ionized by an O-type star in the L935. Further spectral and multi-band observations are guaranteed to investigate in detail the structural features. △ Less

Submitted 29 May, 2022; originally announced May 2022.

Comments: 18 pages, 14 figures, accepted for publication on RAA

Termination of Superradiance from a Binary Companion

Authors: Xi Tong, Yi Wang, Hui-Yu Zhu

Abstract: We study the impact of a binary companion on black hole superradiance at orbital frequencies away from the gravitational-collider-physics (GCP) resonance bands. A superradiant state can couple to a strongly absorptive state via the tidal perturbation of the companion, thereby acquiring a suppressed superradiance rate. Below a critical binary separation, this superradiance rate becomes negative, an… ▽ More We study the impact of a binary companion on black hole superradiance at orbital frequencies away from the gravitational-collider-physics (GCP) resonance bands. A superradiant state can couple to a strongly absorptive state via the tidal perturbation of the companion, thereby acquiring a suppressed superradiance rate. Below a critical binary separation, this superradiance rate becomes negative, and the boson cloud gets absorbed by the black hole. This critical binary separation leads to tight constraints on GCP. Especially, a companion with mass ratio $q>10^{-3}$ invalidates all GCP fine structure transitions, as well as almost all Bohr transitions except those from the $|ψ_{211}\rangle$ state. Meanwhile, the backreaction on the companion manifests itself as a torque acting on the binary, producing floating/sinking orbits that can be verified via pulsar timing. In addition, the possible termination of cloud growth may help to alleviate the current bounds on the ultralight boson mass from various null detections. △ Less

Submitted 21 May, 2022; originally announced May 2022.

Comments: 10 pages, 5 figures

Revealing the Galaxy-Halo Connection Through Machine Learning

Authors: Ryan Hausen, Brant E. Robertson, Hanjue Zhu, Nickolay Y. Gnedin, Piero Madau, Evan E. Schneider, Bruno Villasenor, Nicole E. Drakos

Abstract: Understanding the connections between galaxy stellar mass, star formation rate, and dark matter halo mass represents a key goal of the theory of galaxy formation. Cosmological simulations that include hydrodynamics, physical treatments of star formation, feedback from supernovae, and the radiative transfer of ionizing photons can capture the processes relevant for establishing these connections. T… ▽ More Understanding the connections between galaxy stellar mass, star formation rate, and dark matter halo mass represents a key goal of the theory of galaxy formation. Cosmological simulations that include hydrodynamics, physical treatments of star formation, feedback from supernovae, and the radiative transfer of ionizing photons can capture the processes relevant for establishing these connections. The complexity of these physics can prove difficult to disentangle and obfuscate how mass-dependent trends in the galaxy population originate. Here, we train a machine learning method called Explainable Boosting Machines (EBMs) to infer how the stellar mass and star formation rate of nearly 6 million galaxies simulated by the Cosmic Reionization on Computers (CROC) project depend on the physical properties of halo mass, the peak circular velocity of the galaxy during its formation history $v_\mathrm{peak}$, cosmic environment, and redshift. The resulting EBM models reveal the relative importance of these properties in setting galaxy stellar mass and star formation rate, with $v_\mathrm{peak}$ providing the most dominant contribution. Environmental properties provide substantial improvements for modeling the stellar mass and star formation rate in only $\lesssim10\%$ of the simulated galaxies. We also provide alternative formulations of EBM models that enable low-resolution simulations, which cannot track the interior structure of dark matter halos, to predict the stellar mass and star formation rate of galaxies computed by high-resolution simulations with detailed baryonic physics. △ Less

Submitted 21 April, 2022; originally announced April 2022.

Comments: 27 pages, 18 figures, Submitted to AAS Journals

Unusually high HCO+/CO ratios in and outside supernova remnant W49B

Authors: Ping Zhou, Gao-Yuan Zhang, Xin Zhou, Maria Arias, Bon-Chul Koo, Jacco Vink, Zhi-Yu Zhang, Lei Sun, Fu-Jun Du, Hui Zhu, Yang Chen, Stefano Bovino, Yong-Hyun Lee

Abstract: Galactic supernova remnants (SNRs) and their environments provide the nearest laboratories to study SN feedback. We performed molecular observations toward SNR W49B, the most luminous Galactic SNR in the X-ray band, aiming to explore signs of multiple feedback channels of SNRs on nearby molecular clouds (MCs). We found very broad HCO+ lines with widths of dv = 48--75 km/s in the SNR southwest, pro… ▽ More Galactic supernova remnants (SNRs) and their environments provide the nearest laboratories to study SN feedback. We performed molecular observations toward SNR W49B, the most luminous Galactic SNR in the X-ray band, aiming to explore signs of multiple feedback channels of SNRs on nearby molecular clouds (MCs). We found very broad HCO+ lines with widths of dv = 48--75 km/s in the SNR southwest, providing strong evidence that W49B is perturbing MCs at a systemic velocity of $V_{LSR}=61$--65 km/s, and placing W49B at a distance of $7.9\pm 0.6$ kpc. We observed unusually high-intensity ratios of HCO+ J=1-0/CO J=1-0 not only at shocked regions ($1.1\pm 0.4$ and $0.70\pm 0.16$), but also in quiescent clouds over 1 pc away from the SNR's eastern boundary (> 0.2). By comparing with the magnetohydrodynamics shock models, we interpret that the high ratio in the broad-line regions can result from a cosmic-ray (CR) induced chemistry in shocked MCs, where the CR ionization rate is enhanced to around 10--100 times of the Galactic level. The high HCO+/CO ratio outside the SNR is probably caused by the radiation precursor, while the luminous X-ray emission of W49B can explain a few properties in this region. The above results provide observational evidence that SNRs can strongly influence the molecular chemistry in and outside the shock boundary via their shocks, CRs, and radiation. We propose that the HCO+/CO ratio is a potentially useful tool to probe an SNR's multichannel influence on MCs. △ Less

Submitted 25 May, 2022; v1 submitted 24 March, 2022; originally announced March 2022.

Comments: 19 pages, 11 figures, ApJ in press

Extracting high-order cosmological information in galaxy surveys with power spectra

Authors: Yuting Wang, Gong-Bo Zhao, Kazuya Koyama, Will J. Percival, Ryuichi Takahashi, Chiaki Hikage, Héctor Gil-Marín, ChangHoon Hahn, Ruiyang Zhao, Weibing Zhang, Xiaoyong Mu, Yu Yu, Hong-Ming Zhu, Fei Ge

Abstract: The reconstruction method was proposed more than a decade ago to boost the signal of baryonic acoustic oscillations measured in galaxy redshift surveys, which is one of key probes for dark energy. After moving the observed overdensities in galaxy surveys back to their initial position, the reconstructed density field is closer to a linear Gaussian field, with higher-order information moved back in… ▽ More The reconstruction method was proposed more than a decade ago to boost the signal of baryonic acoustic oscillations measured in galaxy redshift surveys, which is one of key probes for dark energy. After moving the observed overdensities in galaxy surveys back to their initial position, the reconstructed density field is closer to a linear Gaussian field, with higher-order information moved back into the power spectrum. We find that by jointly analysing power spectra measured from the pre- and post-reconstructed galaxy samples, higher-order information beyond the $2$-point power spectrum can be efficiently extracted, which generally yields an information gain upon the analysis using the pre- or post-reconstructed galaxy sample alone. This opens a window to easily use higher-order information when constraining cosmological models. △ Less

Submitted 9 April, 2024; v1 submitted 10 February, 2022; originally announced February 2022.

Comments: 38 pages including 4 figures and 1 table, with a Supplement Information document attached; Version accepted to Communications Physics

Journal ref: Communications Physics, 7, 130 (2024)

Design, Uncertainty Analysis and Measurement of a Silicon-based Platelet THz Corrugated Horn

Authors: Jie Hu, Wei-Tao Lv, Hao-tian Zhu, Zheng Lou, Dong Liu, Jiang-Qiao Ding, Sheng-Cai Shi

Abstract: Platelets corrugated horn is a promising technology for their scalability to a large corrugated horn array. In this paper, we present the design, fabrication, measurement and uncertainty analysis of a wideband 170-320 GHz platelet corrugated horn that features with low sidelobe across the band (<-30 dB). We also propose an accurate and universal method to analyze the axial misalignment of the plat… ▽ More Platelets corrugated horn is a promising technology for their scalability to a large corrugated horn array. In this paper, we present the design, fabrication, measurement and uncertainty analysis of a wideband 170-320 GHz platelet corrugated horn that features with low sidelobe across the band (<-30 dB). We also propose an accurate and universal method to analyze the axial misalignment of the platelets for the first time. It is based on the mode matching (MM) method with a closed-form solution to off-axis circular waveguide discontinuities obtained by using Graf addition theorem for the Bessel functions. The uncertainties introduced in the fabrication have been quantitatively analyzed using the Monte Carlo method. The analysis shows the cross-polarization of the corrugated horn degrades significantly with the axial misalignment. It well explains the discrepancy between the designed and the measured cross-polarization of platelets corrugated horn fabricated in THz band. The method can be used to determine the fabrication tolerance needed for other THz corrugated horns and evaluate the impact of the corrugated horn for astronomical observations. △ Less

Submitted 3 January, 2022; originally announced January 2022.

Comments: Accepted by IEEE Transaction on Antenna and Propagation

Probing the local interstellar medium with scintillometry of the bright pulsar B1133+16

Authors: James W. McKee, Hengrui Zhu, Daniel R. Stinebring, James M. Cordes

Abstract: The interstellar medium hosts a population of scattering screens, most of unknown origin. Scintillation studies of pulsars provides a sensitive tool for resolving these scattering screens and a means of measuring their properties. In this paper, we report our analysis of 34 years of Arecibo observations of PSR B1133+16, from which we have obtained high-quality dynamic spectra and their associated… ▽ More The interstellar medium hosts a population of scattering screens, most of unknown origin. Scintillation studies of pulsars provides a sensitive tool for resolving these scattering screens and a means of measuring their properties. In this paper, we report our analysis of 34 years of Arecibo observations of PSR B1133+16, from which we have obtained high-quality dynamic spectra and their associated scintillation arcs, arising from the scattering screens located along the line of sight to the pulsar. We have identified six individual scattering screens that are responsible for the observed scintillation arcs, which persist for decades. Using the assumption that the scattering screens have not changed significantly in this time, we have modeled the variations in arc curvature throughout the Earth's orbit and extracted information about the placement, orientation, and velocity of five of the six screens, with the highest-precision distance measurement placing a screen at just $5.46^{+0.54}_{-0.59}$ pc from the Earth. We associate the more distant of these screens with an under-dense region of the Local Bubble. △ Less

Submitted 22 December, 2021; originally announced December 2021.

Comments: 16 pages, 9 figures, accepted for publication in ApJ

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

Fermi-LAT detection of extended gamma-ray emission in the vicinity of SNR G045.7$-$00.4: evidence for escaping cosmic rays interacting with the surrounding molecular clouds

Authors: Hai-Ming Zhang, Ruo-Yu Liu, Yang Su, Hui Zhu, Shao-Qiang Xi, Xiang-Yu Wang

Abstract: We present the analysis of Fermi Large Area Telecope (LAT) data of the gamma-ray emission in the vicinity of a radio supernova remnant (SNR), G045.7$-$00.4. To study the origin of the gamma-ray emission, we also make use of the CO survey data of Milky Way Imaging Scroll Painting to study the massive molecular gas complex that surrounds the SNR. The whole size of the GeV emission is significantly l… ▽ More We present the analysis of Fermi Large Area Telecope (LAT) data of the gamma-ray emission in the vicinity of a radio supernova remnant (SNR), G045.7$-$00.4. To study the origin of the gamma-ray emission, we also make use of the CO survey data of Milky Way Imaging Scroll Painting to study the massive molecular gas complex that surrounds the SNR. The whole size of the GeV emission is significantly larger than that of the radio morphology. Above 3 GeV, the GeV emission is resolved into two sources: one is spatially consistent with the position of the SNR with a size comparable to that of the radio emission, and the other is located outside of the western boundary of the SNR and spatially coincident with the densest region of the surrounding molecular cloud. We suggest that the GeV emission of the western source may arise from cosmic rays (CRs) which have escaped the SNR and illuminated the surrounding molecular cloud. We find that the gamma-ray spectra of the western source can be consistently explained by this scenario with a total energy of $\sim 10^{50}{\rm erg}$ in escaping CRs assuming the escape is isotropic. △ Less

Submitted 15 November, 2021; v1 submitted 5 November, 2021; originally announced November 2021.

Comments: 11 pages, 7 figures and 1 table -Accepted for publication in ApJ

The Data Processing of the LAMOST Medium-Resolution Spectral Survey of Galactic Nebulae (LAMOST MRS-N Pipeline)

Authors: Chao-Jian Wu, Hong Wu, Wei Zhang, Yao Li, Juan-Juan Ren, Jian-Jun Chen, Chih-Hao Hsia, Yu-Zhong Wu, Hui Zhu, Bin Li, Yong-Hui Hou

Abstract: The Large sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) medium-resolution spectral survey of Galactic Nebulae (MRS-N) has conducted for three years since Sep. 2018 and observed more than 190 thousands nebular spectra and 20 thousands stellar spectra. However, there is not yet a data processing pipeline for nebular data. To significantly improve the accuracy of nebulae classification… ▽ More The Large sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) medium-resolution spectral survey of Galactic Nebulae (MRS-N) has conducted for three years since Sep. 2018 and observed more than 190 thousands nebular spectra and 20 thousands stellar spectra. However, there is not yet a data processing pipeline for nebular data. To significantly improve the accuracy of nebulae classification and their physical parameters, we developed the MRS-N Pipeline. This article presented in detail each data processing step of the MRS-N Pipeline, such as removing cosmic rays, merging single exposure, fitting sky light emission lines, subtracting skylight, wavelength recalibration, measuring nebular parameters, creating catalogs and packing spectra. Finally, a description of the data products, including nebular spectra files and parameter catalogs, is provided. △ Less

Submitted 19 May, 2022; v1 submitted 9 October, 2021; originally announced October 2021.

Comments: Accepted for publication in RAA. 17 pages, 9 figures, 1 table

LAMOST Medium-Resolution Spectral Survey of Galactic Nebulae (LAMOST-MRS-N): Subtraction of Geocoronal Halpha Emission

Authors: Wei Zhang, Hong Wu, Chaojian Wu, Juanjuan Ren, Jianjun Chen, Chih-Hao Hsia, Yuzhong Wu, Hui Zhu, Jianrong Shi, Zhongrui Bai, Zhaoxiang Qi, Yongheng Zhao, Yonghui Hou

Abstract: We introduce a method of subtracting geocoronal Halpha emissions from the spectra of LAMOST medium-resolution spectral survey of Galactic nebulae (LAMOST-MRS-N). The flux ratios of the Halpha sky line to the adjacent OH lambda6554 single line do not show a pattern or gradient distribution in a plate. More interestingly, the ratio is well correlated to solar altitude, which is the angle of the sun… ▽ More We introduce a method of subtracting geocoronal Halpha emissions from the spectra of LAMOST medium-resolution spectral survey of Galactic nebulae (LAMOST-MRS-N). The flux ratios of the Halpha sky line to the adjacent OH lambda6554 single line do not show a pattern or gradient distribution in a plate. More interestingly, the ratio is well correlated to solar altitude, which is the angle of the sun relative to the Earth's horizon. It is found that the ratio decreases from 0.8 to 0.2 with the decreasing solar altitude from -17 to -73 degree. Based on this relation, which is described by a linear function, we can construct the Halpha sky component and subtract it from the science spectrum. This method has been applied to the LAMOST-MRS-N data, and the contamination level of the Halpha sky to nebula is reduced from 40% to less than 10%. The new generated spectra will significantly improve the accuracy of the classifications and the measurements of physical parameters of Galactic nebulae. △ Less

Submitted 18 August, 2021; originally announced August 2021.

Comments: 12 pages, 10 figures, accepted for publication in Research in Astronomy and Astrophysics

Journal ref: RAA 2021 Vol. 21 No. 11, 280

Cosmic Tidal Reconstruction with Halo Fields

Authors: Hong-Ming Zhu, Tian-Xiang Mao, Ue-Li Pen

Abstract: The gravitational coupling between large-scale perturbations and small-scale perturbations leads to anisotropic distortions of the small-scale matter distribution. The measured local small-scale power spectrum can thus be used to infer the large-scale matter distribution. In this paper, we present a new tidal reconstruction algorithm for reconstructing large-scale modes using the full three-dimens… ▽ More The gravitational coupling between large-scale perturbations and small-scale perturbations leads to anisotropic distortions of the small-scale matter distribution. The measured local small-scale power spectrum can thus be used to infer the large-scale matter distribution. In this paper, we present a new tidal reconstruction algorithm for reconstructing large-scale modes using the full three-dimensional tidal shear information. We apply it to simulated dark matter halo fields and the reconstructed large-scale density field correlates well with the original matter density field on large scales, improving upon the previous tidal reconstruction method which only uses two transverse shear fields. This has profound implications for recovering lost 21~cm radial modes due to foreground subtraction and constraining primordial non-Gaussianity using the multi-tracer method with future cosmological surveys. △ Less

Submitted 3 August, 2021; originally announced August 2021.

Comments: 10 pages, 5 figures

Gravitational Collider Physics via Pulsar-Black Hole Binaries II: Fine and Hyperfine Structures are Favored

Authors: Xi Tong, Yi Wang, Hui-Yu Zhu

Abstract: A rotating black hole can be clouded by light bosons via superradiance, and thus acquire an atom-like structure. If such a gravitational atom system is companioned with a pulsar, the pulsar can trigger transitions between energy levels of the gravitational atom, and these transitions can be detected by pulsar timing. We show that in such pulsar-black hole systems, fine and hyperfine structure tran… ▽ More A rotating black hole can be clouded by light bosons via superradiance, and thus acquire an atom-like structure. If such a gravitational atom system is companioned with a pulsar, the pulsar can trigger transitions between energy levels of the gravitational atom, and these transitions can be detected by pulsar timing. We show that in such pulsar-black hole systems, fine and hyperfine structure transitions are more likely to be probed than the Bohr transition. Also, the calculation of these fine and hyperfine structure transitions are under better analytic control. Thus, these fine and hyperfine structure transitions are more ideal probes in the search for gravitational collider signals in pulsar-black hole systems. △ Less

Submitted 11 July, 2021; v1 submitted 25 June, 2021; originally announced June 2021.

Comments: 9 pages, 4 figures; v2: references added

Element abundance analysis of the metal-rich stellar halo and high-velocitythick disk in the galaxy

Authors: Haifan Zhu, Cuihua Du, Yepeng Yan, Jianrong Shi, Jun Ma, Heidi Jo Newberg

Abstract: Based on the second Gaia data release (DR2) and APOGEE (DR16) spectroscopic surveys, wedefined two kinds of star sample: high-velocity thick disk (HVTD) with $vφ>90km/s$ and metal-richstellar halo (MRSH) with $vφ<90km/s$. Due to high resolution spectra data from APOGEE (DR16),we can analyze accurately the element abundance distribution of HVTD and MRSH. These elementsabundance constituted a multid… ▽ More Based on the second Gaia data release (DR2) and APOGEE (DR16) spectroscopic surveys, wedefined two kinds of star sample: high-velocity thick disk (HVTD) with $vφ>90km/s$ and metal-richstellar halo (MRSH) with $vφ<90km/s$. Due to high resolution spectra data from APOGEE (DR16),we can analyze accurately the element abundance distribution of HVTD and MRSH. These elementsabundance constituted a multidimensional data space, and we introduced an algorithm method forprocessing multi-dimensional data to give the result of dimensionality reduction clustering. Accordingto chemical properties analysis, we derived that some HVTD stars could origin from the thin disk,and some MRSH stars from dwarf galaxies, but those stars which have similar chemical abundancecharacteristics in both sample may form in-situ. △ Less

Submitted 7 May, 2021; originally announced May 2021.

Comments: 11 pages, 8 figures, accepted for publication in ApJ