Showing 1–50 of 657 results for author: Li, B

Advancing Gamma-Ray Burst Identification through Transfer Learning with Convolutional Neural Networks

Authors: Peng Zhang, Bing Li, Ren-zhou Gui, Shao-lin Xiong, Yu Wang, Yan-qiu Zhang, Chen-wei Wang, Jia-cong Liu, Wang-chen Xue, Chao Zheng, Zheng-hang Yu, Wen-long Zhang

Abstract: The Rapid and accurate identification of Gamma-Ray Bursts (GRBs) is crucial for unraveling their origins. However, current burst search algorithms frequently miss low-threshold signals or lack universality for observations. In this study, we propose a novel approach utilizing transfer learning experiment based on convolutional neural network (CNN) to establish a universal GRB identification method… ▽ More The Rapid and accurate identification of Gamma-Ray Bursts (GRBs) is crucial for unraveling their origins. However, current burst search algorithms frequently miss low-threshold signals or lack universality for observations. In this study, we propose a novel approach utilizing transfer learning experiment based on convolutional neural network (CNN) to establish a universal GRB identification method, which validated successfully using GECAM-B data. By employing data augmentation techniques, we enhance the diversity and quantity of the GRB sample. We develop a 1D CNN model with a multi-scale feature cross fusion module (MSCFM) to extract features from samples and perform classification. The comparative results demonstrated significant performance improvements following pre-training and transferring on a large-scale dataset. Our optimal model achieved an impressive accuracy of 96.41% on the source dataset of GECAM-B, and identified three previously undiscovered GRBs by contrast with manual analysis of GECAM-B observations. These innovative transfer learning and data augmentation methods presented in this work hold promise for applications in multi-satellite exploration scenarios characterized by limited data sets and a scarcity of labeled samples in high-energy astronomy. △ Less

Submitted 24 August, 2024; originally announced August 2024.

Comments: 17 pages, 7 figures

Investigating the competition between the deconfinement and chiral phase transitions in light of the multimessenger observations of neutron stars

Authors: Wen-Li Yuan, Bikai Gao, Yan Yan, Bolin Li, Renxin Xu

Abstract: We extend the parity doublet model for hadronic matter and study the possible presence of quark matter inside the cores of neutron stars with the Nambu-Jona-Lasinio (NJL) model. Considering the uncertainties of the QCD phase diagram and the location of the critical endpoint, we aim to explore the competition between the chiral phase transition and the deconfinement phase transition systematically,… ▽ More We extend the parity doublet model for hadronic matter and study the possible presence of quark matter inside the cores of neutron stars with the Nambu-Jona-Lasinio (NJL) model. Considering the uncertainties of the QCD phase diagram and the location of the critical endpoint, we aim to explore the competition between the chiral phase transition and the deconfinement phase transition systematically, regulated by the vacuum pressure $-B$ in the NJL model. Employing a Maxwell construction, a sharp first-order deconfinement phase transition is implemented combining the parity doublet model for the hadronic phase and the NJL model for the high-energy quark phase. The position of the chiral phase transition is obtained from the NJL model self-consistently. We find stable neutron stars with a quark core within a specific parameter space that satisfies current astronomical observations. The observations suggest a relatively large chiral invariant mass $m_0=600$ MeV in the parity doublet model and a larger split between the chiral and deconfinement phase transitions while assuming the first-order deconfinement phase transition. The maximum mass of the hybrid star that we obtain is $\sim 2.2 M_{\odot}$. △ Less

Submitted 12 August, 2024; v1 submitted 10 August, 2024; originally announced August 2024.

Comments: 10pages,7 figures

Emulation of $f(R)$ modified gravity from $Λらむだ$CDM using conditional GANs

Authors: Yash Gondhalekar, Sownak Bose, Baojiu Li, Carolina Cuesta-Lazaro

Abstract: A major aim of current and upcoming cosmological surveys is testing deviations from the standard $Λらむだ$CDM model, but the full scientific value of these surveys will only be realised through efficient simulation methods that keep up with the increasing volume and precision of observational data. $N$-body simulations of modified gravity (MG) theories are computationally expensive since highly non-line… ▽ More A major aim of current and upcoming cosmological surveys is testing deviations from the standard $Λらむだ$CDM model, but the full scientific value of these surveys will only be realised through efficient simulation methods that keep up with the increasing volume and precision of observational data. $N$-body simulations of modified gravity (MG) theories are computationally expensive since highly non-linear equations need to be solved to model the non-linear matter evolution; this represents a significant bottleneck in the path to reach the data volume and resolution attained by equivalent $Λらむだ$CDM simulations. We develop a field-level, neural-network-based emulator that generates density and velocity divergence fields under the $f(R)$ gravity MG model from the corresponding $Λらむだ$CDM simulated fields. Using attention mechanisms and a complementary frequency-based loss function, our model is able to learn this intricate mapping. We further use the idea of latent space extrapolation to generalise our emulator to $f(R)$ models with differing field strengths. The predictions of our emulator agree with the $f(R)$ simulations to within 5% for matter density and to within 10% for velocity divergence power spectra up to $k \sim 2\,h$ $\mathrm{Mpc}^{-1}$. But for a few select cases, higher-order statistics are reproduced with $\lesssim$10% agreement with the $f(R)$ simulations. Latent extrapolation allows our emulator to generalise to different $f(R)$ model variants without explicitly training on those variants. Given a $Λらむだ$CDM simulation, the GPU-based emulator is able to reproduce the equivalent $f(R)$ realisation $\sim$600 times faster than full $N$-body simulations. This lays the foundations for a valuable tool for realistic yet rapid mock field generation and robust cosmological analyses. △ Less

Submitted 22 July, 2024; originally announced July 2024.

Comments: 20 pages, 10 figures. To be submitted to MNRAS

Negative neutrino masses as a mirage of dark energy

Authors: Willem Elbers, Carlos S. Frenk, Adrian Jenkins, Baojiu Li, Silvia Pascoli

Abstract: The latest cosmological constraints on the sum of the neutrino masses depend on prior physical assumptions about the mass spectrum. To test the accordance of cosmological and laboratory constraints in the absence of such priors, we introduce an effective neutrino mass parameter that extends consistently to negative values. For the $Λらむだ$CDM model, we analyze data from Planck, ACT, and DESI and find a… ▽ More The latest cosmological constraints on the sum of the neutrino masses depend on prior physical assumptions about the mass spectrum. To test the accordance of cosmological and laboratory constraints in the absence of such priors, we introduce an effective neutrino mass parameter that extends consistently to negative values. For the $Λらむだ$CDM model, we analyze data from Planck, ACT, and DESI and find a $2.8-3.3σしぐま$ tension with the constraints from oscillation experiments. Motivated by recent hints of evolving dark energy, we analyze the $w_0w_a$ and mirage dark energy models, showing that they favour larger masses consistent with laboratory data, respectively $\sum m_{νにゅー,\mathrm{eff}} = 0.06_{-0.10}^{+0.15}\,\mathrm{eV}$ and $\sum m_{νにゅー,\mathrm{eff}} = 0.04_{-0.11}^{+0.15}\,\mathrm{eV}$ (both at 68\%). △ Less

Submitted 15 July, 2024; originally announced July 2024.

Comments: 7 pages, 5 figures

Bayesian Inference of Fine-Features of Nuclear Equation of State from Future Neutron Star Radius Measurements to 0.1km Accuracy

Authors: Bao-An Li, Xavier Grundler, Wen-Jie Xie, Nai-Bo Zhang

Abstract: To more precisely constrain the Equation of State (EOS) of supradense neutron-rich nuclear matter, future high-precision X-ray and gravitational wave observatories are proposed to measure the radii of neutron stars (NSs) with an accuracy better than about 0.1 km. However, it remains unclear what particular aspects (other than the stiffness generally spoken of in the literature) of the EOS and to w… ▽ More To more precisely constrain the Equation of State (EOS) of supradense neutron-rich nuclear matter, future high-precision X-ray and gravitational wave observatories are proposed to measure the radii of neutron stars (NSs) with an accuracy better than about 0.1 km. However, it remains unclear what particular aspects (other than the stiffness generally spoken of in the literature) of the EOS and to what precision they will be better constrained. In this work, within a Bayesian framework using a meta-model EOS for NSs, we infer the posterior probability distribution functions (PDFs) of incompressibility $K_{0}$ and skewness $J_{0}$ of symmetric nuclear matter (SNM) as well as the slope $L$, curvature $K_{\rm{sym}}$, and skewness $J_{\rm{sym}}$ characterizing the density dependence of nuclear symmetry energy $E_{\rm{sym}}(ρろー)$, respectively, from mean values of NS radii consistent with existing observations and an expected accuracy $ΔでるたR$ ranging from about 1.0 km to 0.1 km. We found that (1) the $ΔでるたR$ has little effect on inferring the stiffness of SNM at suprasaturation densities, (2) smaller $ΔでるたR$ reveals more accurately not only the PDFs but also pairwise correlations among parameters characterizing high-density $E_{\rm{sym}}(ρろー)$, (3) a double-peak feature of the PDF($K_{\rm{sym}}$) corresponding to the strong $K_{\rm{sym}}-J_{\rm{sym}}$ and $K_{\rm{sym}}-L$ anti-correlations is revealed when $ΔでるたR$ is less than about 0.2 km, (4) the high-precision radius measurement for canonical NSs is more useful than that for massive ones for constraining the EOS of nucleonic matter around (2-3) times the saturation density of SNM. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Comments: 21 pages with 14 figures

A new subclass of gamma-ray burst originating from compact binary merger

Authors: Chen-Wei Wang, Wen-Jun Tan, Shao-Lin Xiong, Shu-Xu Yi, Rahim Moradi, Bing Li, Zhen Zhang, Yu Wang, Yan-Zhi Meng, Jia-Cong Liu, Yue Wang, Sheng-Lun Xie, Wang-Chen Xue, Zheng-Hang Yu, Peng Zhang, Wen-Long Zhang, Yan-Qiu Zhang, Chao Zheng

Abstract: Type I gamma-ray bursts (GRBs) are believed to originate from compact binary merger usually with duration less than 2 seconds for the main emission. However, recent observations of GRB 211211A and GRB 230307A indicate that some merger-origin GRBs could last much longer. Since they show strikingly similar properties (indicating a common mechanism) which are different from the classic "long"-short b… ▽ More Type I gamma-ray bursts (GRBs) are believed to originate from compact binary merger usually with duration less than 2 seconds for the main emission. However, recent observations of GRB 211211A and GRB 230307A indicate that some merger-origin GRBs could last much longer. Since they show strikingly similar properties (indicating a common mechanism) which are different from the classic "long"-short burst (e.g. GRB 060614), forming an interesting subclass of type I GRBs, we suggest to name them as type IL GRBs. By identifying the first peak of GRB 230307A as a quasi-thermal precursor, we find that the prompt emission of type IL GRB is composed of three episodes: (1) a precursor followed by a short quiescent (or weak emission) period, (2) a long-duration main emission, and (3) an extended emission. With this burst pattern, a good candidate, GRB 170228A, was found in the Fermi/GBM archive data, and subsequent temporal and spectral analyses indeed show that GRB 170228A falls in the same cluster with GRB 211211A and GRB 230307A in many diagnostic figures. Thus this burst pattern could be a good reference for rapidly identifying type IL GRB and conducting low-latency follow-up observation. We estimated the occurrence rate and discussed the physical origins and implications for the three emission episodes of type IL GRBs. Our analysis suggests the pre-merger precursor model, especially the super flare model, is more favored for type IL GRBs. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Galaxy clustering in modified gravity from full-physics simulations. I: two-point correlation functions

Authors: Michael Collier, Sownak Bose, Baojiu Li

Abstract: We present an in-depth investigation of galaxy clustering based on a new suite of realistic large-box galaxy-formation simulations in $f(R)$ gravity, with a subgrid physics model that has been recalibrated to reproduce various observed stellar and gas properties. We focus on the two-point correlation functions of the luminous red galaxies (LRGs) and emission line galaxies (ELGs), which are primary… ▽ More We present an in-depth investigation of galaxy clustering based on a new suite of realistic large-box galaxy-formation simulations in $f(R)$ gravity, with a subgrid physics model that has been recalibrated to reproduce various observed stellar and gas properties. We focus on the two-point correlation functions of the luminous red galaxies (LRGs) and emission line galaxies (ELGs), which are primary targets of ongoing and future galaxy surveys such as DESI. One surprising result is that, due to several nontrivial effects of modified gravity on matter clustering and the galaxy-halo connection, the clustering signal does not depend monotonically on the fifth-force strength. For LRGs this complicated behaviour poses a challenge to meaningfully constraining this model. For ELGs, in contrast, this can be straightforwardly explained by the time evolution of the fifth force, which means that weaker $f(R)$ models can display nearly the same -- up to $25\%$ -- deviations from $Λらむだ$CDM as the strongest ones, albeit at lower redshifts. This implies that even very weak $f(R)$ models can be strongly constrained, unlike with most other observations. Our results show that galaxy formation acquires a significant environment dependence in $f(R)$ gravity which, if not properly accounted for, may lead to biased constraints on the model. This highlights the essential role of hydrodynamical simulations in future tests of gravity exploring precision galaxy-clustering data from the likes of DESI and Euclid. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Comments: 18 pages, 8 figures + 2 figures in appendix

Timing and Scintillation Studies of Pulsars in Globular Cluster M3 (NGC 5272) with FAST

Authors: Baoda Li, Li-yun Zhang, Jumei Yao, Dejiang Yin, Ralph P. Eatough, Minghui Li, Yifeng Li, Yujie Lian, Yu Pan, Yinfeng Dai, Yaowei Li, Xingnan Zhang, Tianhao Su, Yuxiao Wu, Tong Liu, Kuo Liu, Lin Wang, Lei Qian, Zhichen Pan

Abstract: We present the phase-connected timing solutions of all the five pulsars in globular cluster (GC) M3 (NGC 5272), namely PSRs M3A to F (PSRs J1342+2822A to F), with the exception of PSR M3C, from FAST archival data. In these timing solutions, those of PSRs M3E, and F are obtained for the first time. We find that PSRs M3E and F have low mass companions, and are in circular orbits with periods of 7.1… ▽ More We present the phase-connected timing solutions of all the five pulsars in globular cluster (GC) M3 (NGC 5272), namely PSRs M3A to F (PSRs J1342+2822A to F), with the exception of PSR M3C, from FAST archival data. In these timing solutions, those of PSRs M3E, and F are obtained for the first time. We find that PSRs M3E and F have low mass companions, and are in circular orbits with periods of 7.1 and 3.0 days, respectively. For PSR M3C, we have not detected it in all the 41 observations. We found no X-ray counterparts for these pulsars in archival Chandra images in the band of 0.2-20 keV. We noticed that the pulsars in M3 seem to be native. From the Auto-Correlation Function (ACF) analysis of the M3A's and M3B's dynamic spectra, the scintillation timescale ranges from $7.0\pm0.3$ min to $60.0\pm0.6$ min, and the scintillation bandwidth ranges from $4.6\pm0.2$ MHz to $57.1\pm1.1$ MHz. The measured scintillation bandwidths from the dynamic spectra indicate strong scintillation, and the scattering medium is anisotropic. From the secondary spectra, we captured a scintillation arc only for PSR M3B with a curvature of $649\pm23 {\rm m}^{-1} {\rm mHz}^{-2}$. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: 14 pages, 4 figures, accepted for publication in The Astrophysical Journal

Matter Power Spectra in Modified Gravity: A Comparative Study of Approximations and $N$-Body Simulations

Authors: Benjamin Bose, Ashim Sen Gupta, Bartolomeo Fiorini, Guilherme Brando, Farbod Hassani, Tessa Baker, Lucas Lombriser, Baojiu Li, Cheng-Zong Ruan, Cesar Hernandez-Aguayo, Luis Atayde, Noemi Frusciante

Abstract: Testing gravity and the concordance model of cosmology, $Λらむだ$CDM, at large scales is a key goal of this decade's largest galaxy surveys. Here we present a comparative study of dark matter power spectrum predictions from different numerical codes in the context of three popular theories of gravity that induce scale-independent modifications to the linear growth of structure: nDGP, Cubic Galileon and… ▽ More Testing gravity and the concordance model of cosmology, $Λらむだ$CDM, at large scales is a key goal of this decade's largest galaxy surveys. Here we present a comparative study of dark matter power spectrum predictions from different numerical codes in the context of three popular theories of gravity that induce scale-independent modifications to the linear growth of structure: nDGP, Cubic Galileon and K-mouflage. In particular, we compare the predictions from full $N$-body simulations, two $N$-body codes with approximate time integration schemes, a parametrised modified $N$-body implementation and the analytic halo model reaction approach. We find the modification to the $Λらむだ$CDM spectrum is in $2\%$ agreement for $z\leq1$ and $k\leq 1~h/{\rm Mpc}$ over all gravitational models and codes, in accordance with many previous studies, indicating these modelling approaches are robust enough to be used in forthcoming survey analyses under appropriate scale cuts. We further make public the new code implementations presented, specifically the halo model reaction K-mouflage implementation and the relativistic Cubic Galileon implementation. △ Less

Submitted 19 June, 2024; originally announced June 2024.

Comments: 20 pages, 4 figures, 4 tables

Constraining modified gravity with weak lensing peaks

Authors: Christopher T. Davies, Joachim Harnois-Déraps, Baojiu Li, Benjamin Giblin, César Hernández-Aguayo, Enrique Paillas

Abstract: It is well established that maximizing the information extracted from upcoming and ongoing stage-IV weak-lensing surveys requires higher-order summary statistics that complement the standard two-point statistics. In this work, we focus on weak-lensing peak statistics to test two popular modified gravity models, $f(R)$ and nDGP, using the FORGE and BRIDGE weak-lensing simulations, respectively. Fro… ▽ More It is well established that maximizing the information extracted from upcoming and ongoing stage-IV weak-lensing surveys requires higher-order summary statistics that complement the standard two-point statistics. In this work, we focus on weak-lensing peak statistics to test two popular modified gravity models, $f(R)$ and nDGP, using the FORGE and BRIDGE weak-lensing simulations, respectively. From these simulations we measure the peak statistics as a function of both cosmological and modified gravity parameters simultaneously. Our findings indicate that the peak abundance is sensitive to the strength of modified gravity, while the peak two-point correlation function is sensitive to the nature of the screening mechanism in a modified gravity model. We combine these simulated statistics with a Gaussian Process Regression emulator and a Gaussian likelihood to generate stage-IV forecast posterior distributions for the modified gravity models. We demonstrate that, assuming small scales can be correctly modelled, peak statistics can be used to distinguish GR from $f(R)$ and nDGP models at the two-sigma level with a stage-IV survey area of $300 \, \rm{deg}^2$ and $1000 \, \rm{deg}^2$, respectively. Finally, we show that peak statistics can constrain $\log_{10}\left(|f_{R0}|\right) = -6$ to 2\% precision, and $\log_{10}(H_0 r_c) = 0.5$ to 25\% precision. △ Less

Submitted 1 August, 2024; v1 submitted 17 June, 2024; originally announced June 2024.

Comments: 25 pages, 16 figures

ExaGRyPE: Numerical General Relativity Solvers Based upon the Hyperbolic PDEs Solver Engine ExaHyPE

Authors: Han Zhang, Baojiu Li, Tobias Weinzierl, Cristian Barrera-Hinojosa

Abstract: ExaGRyPE describes a suite of solvers for numerical relativity, based upon ExaHyPE 2, the second generation of our Exascale Hyperbolic PDE Engine. The presented generation of ExaGRyPE solves the Einstein equations in the CCZ4 formulation under a 3+1 foliation and focuses on black hole spacetimes. It employs a block-structured Cartesian grid carrying a higher-order order Finite Difference scheme wi… ▽ More ExaGRyPE describes a suite of solvers for numerical relativity, based upon ExaHyPE 2, the second generation of our Exascale Hyperbolic PDE Engine. The presented generation of ExaGRyPE solves the Einstein equations in the CCZ4 formulation under a 3+1 foliation and focuses on black hole spacetimes. It employs a block-structured Cartesian grid carrying a higher-order order Finite Difference scheme with adaptive mesh refinement, it facilitates massive parallelism combining message passing, domain decomposition and task parallelism, and it supports the injection of particles into the grid as data probes or tracers. We introduce the ExaGRyPE-specific building blocks within ExaHyPE 2, and discuss its software architecture and compute-n-feel: For this, we formalize the creation of any specific simulation with ExaGRyPE as a sequence of lowering operations, where abstract logical tasks are successively broken into finer tasks until we obtain an abstraction level that can be mapped onto a C++ executable. The overall program logic is fully specified via a domain-specific Python interface, we map this logic onto a more detailed set of numerical tasks, subsequently lower this representation onto technical tasks that the underlying ExaHyPE engine uses to parallelize the application, before eventually the technical tasks are mapped onto task graphs including the actual PDE term evaluations, initial conditions, boundary conditions, and so forth. We end up with a rigorous separation of concerns which shields ExaGRyPE users from technical details and hence simplifies the development of novel physical models. We present the simulations and data for the gauge wave, static single black holes and rotating binary black hole systems, demonstrating that the code base is mature and usable. However, we also uncover domain-specific numerical challenges that need further study by the community in future work. △ Less

Submitted 18 June, 2024; v1 submitted 17 June, 2024; originally announced June 2024.

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

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

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

Submitted 12 June, 2024; originally announced June 2024.

Comments: 17 pages, 12 figures, accepted by PRL

Impact of the nuclear equation of state on the formation of twin stars

Authors: Nai-Bo Zhang, Bao-An Li

Abstract: Twin stars-two stable neutron stars (NSs) with the same mass but different radii have long been proposed to appear as a consequence of a possible first-order phase transition in NS matter. Within a meta-model for the EOS of hybrid stars, we revisit the viability of twin stars and its dependence on numerous parameters characterizing the EOS of nuclear matter, quark matter, and the phase transition… ▽ More Twin stars-two stable neutron stars (NSs) with the same mass but different radii have long been proposed to appear as a consequence of a possible first-order phase transition in NS matter. Within a meta-model for the EOS of hybrid stars, we revisit the viability of twin stars and its dependence on numerous parameters characterizing the EOS of nuclear matter, quark matter, and the phase transition between them. While essentially no experimental constraint exists for the last two, parameters characterizing the EOS of neutron-rich nucleonic matter have been constrained within various ranges by terrestrial experiments and astrophysical observations. Within these ranges, the impact of nuclear EOS on the formation of twin stars is studied. It is found that the symmetry energy of neutron-rich nucleonic matter notably influences the formation of twin stars, particularly through its slope $L$ and curvature $K_{\rm sym}$. Conversely, the EOS of symmetric nuclear matter shows minimal influence on the formation of twin stars, partially owing to its relatively well-constrained parameters. △ Less

Submitted 11 June, 2024; originally announced June 2024.

Comments: 9 pages including 7 figures

Unraveling Trace Anomaly of Supradense Matter via Neutron Star Compactness Scaling

Authors: Bao-Jun Cai, Bao-An Li

Abstract: The trace anomaly $Δでるた\equiv 1/3-P/\varepsilon$ quantifies the possibly broken conformal symmetry in supradense matter under pressure $P$ at energy density $\varepsilon$. Perturbative QCD (pQCD) predicts a vanishing $Δでるた$ at extremely high energy or baryon densities when the conformal symmetry is realized but its behavior at intermediate densities reachable in neutron stars (NSs) are still very uncert… ▽ More The trace anomaly $Δでるた\equiv 1/3-P/\varepsilon$ quantifies the possibly broken conformal symmetry in supradense matter under pressure $P$ at energy density $\varepsilon$. Perturbative QCD (pQCD) predicts a vanishing $Δでるた$ at extremely high energy or baryon densities when the conformal symmetry is realized but its behavior at intermediate densities reachable in neutron stars (NSs) are still very uncertain. The extraction of $Δでるた$ from NS observations strongly depends on the employed model for nuclear Equation of State (EOS). Based on the analytical results from perturbatively analyzing the dimensionless Tolman-Oppenheimer-Volkoff (TOV) equations that are further verified numerically by using $10^5$ EOSs generated randomly with a meta-model in a very broad EOS parameter space constrained by terrestrial nuclear experiments and astrophysical observations, here we first show that the compactness $ξくしー\equiv M_{\rm{NS}}/R$ of a NS with mass $M_{\rm{NS}}$ and radius $R$ scales very accurately with $Πぱい_{\rm{c}}\equiv\mathrm{X}/(1+3\mathrm{X}^2+4\mathrm{X})$ where $\mathrm{X}\equiv P_{\rm{c}}/\varepsilon_{\rm{c}}$ is the ratio of pressure over energy density at NS centers. The scaling of NS compactness thus enables one to readily read off the central trace anomaly $Δでるた_{\rm{c}}=1/3-\mathrm{X}$ directly from the observational data of either the mass-radius or red-shift measurements. We then demonstrate indeed that the available NS data themselves from recent X-ray and gravitational wave observations can determine model-independently the trace anomaly as a function of energy density in NS cores, providing a stringent test of existing NS models and a clear guidance in a new direction for further understanding the nature and EOS of supradense matter. △ Less

Submitted 7 June, 2024; originally announced June 2024.

Comments: 5 pages of main text with 4 pages of supplementary materials

Modelling the redshift-space cluster-galaxy correlation function on Mpc scales with emulation of the pairwise velocity distribution

Authors: Andrew Robertson, Eric Huff, Katarina Markovic, Baojiu Li

Abstract: We present a method for modelling the cluster-galaxy correlation function in redshift-space, down to ~ Mpc scales. The method builds upon the so-called Galaxy Infall Kinematics (GIK) model, a parametric model for the pairwise velocities of galaxies with respect to nearby galaxy clusters. We fit the parameters of the GIK model to a suite of simulations run with different cosmologies, and use Gaussi… ▽ More We present a method for modelling the cluster-galaxy correlation function in redshift-space, down to ~ Mpc scales. The method builds upon the so-called Galaxy Infall Kinematics (GIK) model, a parametric model for the pairwise velocities of galaxies with respect to nearby galaxy clusters. We fit the parameters of the GIK model to a suite of simulations run with different cosmologies, and use Gaussian Processes to emulate how the GIK parameters depend upon cosmology. This emulator can then be combined with knowledge of the real-space clustering of clusters and galaxies, to predict the cluster-galaxy correlation function in redshift space. Fitting this model to an observed correlation function enables the extraction of cosmological parameter constraints, and we present forecasts for a DESI-like survey. We also perform tests of the robustness of our constraints from fitting to mock data extracted from N-body simulations, finding that fitting to scales < 3 Mpc/h leads to a biased inference on cosmology, due to model misspecification on these scales. Finally, we discuss what steps will need to be taken in order to apply our method to real data. △ Less

Submitted 3 June, 2024; originally announced June 2024.

Comments: 24 pages, 11 figures, under review at MNRAS

FAST Discovery of Eight Isolated Millisecond Pulsars in NGC 6517

Authors: Dejiang Yin, Li-yun Zhang, Lei Qian, Ralph P. Eatough, Baoda Li, Duncan R. Lorimer, Yinfeng Dai, Yaowei Li, Xingnan Zhang, Minghui Li, Tianhao Su, Yuxiao Wu, Yu Pan, Yujie Lian, Tong Liu, Zhen Yan, Zhichen Pan

Abstract: We present the discovery of 8 isolated millisecond pulsars in Globular Cluster (GC) NGC 6517 using the Five-Hundred-meter Aperture Spherical radio Telescope (FAST). The spin periods of those pulsars (namely PSR J1801-0857K to R, or, NGC 6517K to R) are all shorter than 10 ms. With these discoveries, NGC 6517 is currently the GC with the most known pulsars in the FAST sky. The largest difference in… ▽ More We present the discovery of 8 isolated millisecond pulsars in Globular Cluster (GC) NGC 6517 using the Five-Hundred-meter Aperture Spherical radio Telescope (FAST). The spin periods of those pulsars (namely PSR J1801-0857K to R, or, NGC 6517K to R) are all shorter than 10 ms. With these discoveries, NGC 6517 is currently the GC with the most known pulsars in the FAST sky. The largest difference in dispersion measure of the pulsars in NGC 6517 is 11.2 cm$^{-3}$ pc, the second among all GCs. The fraction of isolated pulsars in this GC (16 of 17, 94$\%$) is consistent with previous studies indicating an overabundance of isolated pulsars in the densest GCs, especially in those undergoing cluster core collapse. Considering the FAST GC pulsar discoveries, we modeled the GC pulsar population using the empirical Bayesian method described by Turk and Lorimer with the recent counts. Using this approach, we find that the expected number of potential pulsars in GCs seems to be correlated with the central escape velocity, hence, the GCs Liller 1, NGC 6441, M54 (NGC 6715), and $ωおめが$-Cen (NGC 5139) are expected to host the largest numbers of pulsars. △ Less

Submitted 28 May, 2024; originally announced May 2024.

Comments: 21 pages, 2 figures, accepted for publication in The Astrophysical Journal Letters

Gliese 12 b, A Temperate Earth-sized Planet at 12 Parsecs Discovered with TESS and CHEOPS

Authors: Shishir Dholakia, Larissa Palethorpe, Alexander Venner, Annelies Mortier, Thomas G. Wilson, Chelsea X. Huang, Ken Rice, Vincent Van Eylen, Emma Nabbie, Ryan Cloutier, Walter Boschin, David Ciardi, Laetitia Delrez, Georgina Dransfield, Elsa Ducrot, Zahra Essack, Mark E. Everett, Michaël Gillon, Matthew J. Hooton, Michelle Kunimoto, David W. Latham, Mercedes López-Morales, Bin Li, Fan Li, Scott McDermott , et al. (11 additional authors not shown)

Abstract: We report on the discovery of Gliese 12 b, the nearest transiting temperate, Earth-sized planet found to date. Gliese 12 is a bright ($V=12.6$ mag, $K=7.8$ mag) metal-poor M4V star only $12.162\pm0.005$ pc away from the Solar System with one of the lowest stellar activity levels known for an M-dwarf. A planet candidate was detected by TESS based on only 3 transits in sectors 42, 43, and 57, with a… ▽ More We report on the discovery of Gliese 12 b, the nearest transiting temperate, Earth-sized planet found to date. Gliese 12 is a bright ($V=12.6$ mag, $K=7.8$ mag) metal-poor M4V star only $12.162\pm0.005$ pc away from the Solar System with one of the lowest stellar activity levels known for an M-dwarf. A planet candidate was detected by TESS based on only 3 transits in sectors 42, 43, and 57, with an ambiguity in the orbital period due to observational gaps. We performed follow-up transit observations with CHEOPS and ground-based photometry with MINERVA-Australis, SPECULOOS, and Purple Mountain Observatory, as well as further TESS observations in sector 70. We statistically validate Gliese 12 b as a planet with an orbital period of $12.76144\pm0.00006$ days and a radius of $1.0\pm{0.1}$ R$_\oplus$, resulting in an equilibrium temperature of $\sim$315K. Gliese 12 b has excellent future prospects for precise mass measurement, which may inform how planetary internal structure is affected by the stellar compositional environment. Gliese 12 b also represents one of the best targets to study whether Earth-like planets orbiting cool stars can retain their atmospheres, a crucial step to advance our understanding of habitability on Earth and across the Galaxy. △ Less

Submitted 21 May, 2024; originally announced May 2024.

Comments: 19 pages, 7 figures, Accepted for publication in MNRAS, Authors Shishir Dholakia and Larissa Palethorpe contributed equally

Constraining regularization ambiguities in Loop Quantum Cosmology via CMB

Authors: Bao-Fei Li, Meysam Motaharfar, Parampreet Singh

Abstract: In order to investigate the potential observational signals of different regularization ambiguities in loop quantum cosmological models, we systematically compute and compare the primordial scalar power spectra and the resulting angular power spectra in the standard loop quantum cosmology (LQC) and its Thiemann regularized versions -- modified LQC-I/II (mLQC-I/II), using both the dressed metric an… ▽ More In order to investigate the potential observational signals of different regularization ambiguities in loop quantum cosmological models, we systematically compute and compare the primordial scalar power spectra and the resulting angular power spectra in the standard loop quantum cosmology (LQC) and its Thiemann regularized versions -- modified LQC-I/II (mLQC-I/II), using both the dressed metric and the hybrid approaches. All three loop quantum cosmological models yield a non-singular bounce with a post-bounce physics that converges rapidly in a few Planck seconds. Using Starobinsky potential and the initial conditions for the background dynamics chosen to yield the same inflationary e-foldings, which are fixed to be $65$ in all three LQC models, we require that all three models result in the same scale-invariant regime for the primordial power spectrum with a relative difference of less than one percent. This permits us to explore the differences resulting from the deep Planck regime in the angular power spectrum. For the adiabatic states, our results demonstrate that the angular power spectrum predicted by the hybrid approach has a smaller deviation from the angular power spectrum predicted by the standard $Λらむだ$CDM cosmological model at large angles in comparison with the dressed metric approach for all three models. The angular power spectrum predicted by mLQC-I in both the hybrid and the dressed metric approaches shows the smallest deviation from the one predicted by the standard $Λらむだ$CDM cosmological model at large angular scales, except for the case of fourth order adiabatic initial states in the hybrid approach. On the contrary, mLQC-II results in the largest deviations for the amplitude of the angular power spectrum at large angles and is most disfavored. △ Less

Submitted 28 August, 2024; v1 submitted 20 May, 2024; originally announced May 2024.

Comments: 29 pages, 12 figures. Minor revision to match published version in Phys Rev D

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

Reassessment of the dipole in the distribution of quasars on the sky

Authors: Arefe Abghari, Emory F. Bunn, Lukas T. Hergt, Boris Li, Douglas Scott, Raelyn M. Sullivan, Dingchen Wei

Abstract: We investigate claims of an anomalously large amplitude of the dipole in the distribution of quasars on the sky. Two main issues indicate that the systematic uncertainties in the derived quasar-density dipole are underestimated. Firstly, the spatial distribution of the quasars is not a pure dipole, possessing low-order multipoles of comparable size to the dipole. These multipoles are unexpected an… ▽ More We investigate claims of an anomalously large amplitude of the dipole in the distribution of quasars on the sky. Two main issues indicate that the systematic uncertainties in the derived quasar-density dipole are underestimated. Firstly, the spatial distribution of the quasars is not a pure dipole, possessing low-order multipoles of comparable size to the dipole. These multipoles are unexpected and presumably caused by unknown systematic effects; we cannot be confident that the dipole amplitude is not also affected by the same systematics until the origin of these fluctuations is understood. Secondly, the 50 percent sky cut associated with the quasar catalogue strongly couples the multipoles, meaning that the power estimate at ell=1 contains significant contributions from ell>1. In particular, the dominant quadrupole mode in the Galactic mask strongly couples the dipole with the octupole, leading to a large uncertainty in the dipole amplitude. Together these issues mean that the dipole in the quasar catalogue has an uncertainty large enough that consistency with the cosmic microwave background (CMB) dipole cannot be ruled out. More generally, current data sets are insufficiently clean to robustly measure the quasar dipole and future studies will require samples that are larger (preferably covering more of the sky) and free of systematic effects to make strong claims regarding their consistency with the CMB dipole. △ Less

Submitted 15 May, 2024; originally announced May 2024.

Comments: 19 pages, 9 figures

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

A short review on the pulsar magnetic inclination angles

Authors: Biao-Peng Li, Zhi-Fu Gao

Abstract: The inclination angle $χかい$ between magnetic and rotation axes of pulsars is an important parameter in pulsar physics. The changes in the inclination angle of a pulsar would lead to observable effects, such as changes in the pulse beam width and braking index of the star. In this paper, we perform a short review on the evolution of pulsar's magnetic inclination angle, as well as the latest research… ▽ More The inclination angle $χかい$ between magnetic and rotation axes of pulsars is an important parameter in pulsar physics. The changes in the inclination angle of a pulsar would lead to observable effects, such as changes in the pulse beam width and braking index of the star. In this paper, we perform a short review on the evolution of pulsar's magnetic inclination angle, as well as the latest research progress. Using an alignment rotator model in vacuum, we investigate the magnetic inclination angle change rates for 12 high-braking index pulsars without glitch, whose timing observations are obtained using the Nanshan 25-m Radio Telescope at Xinjiang Astronomical Observatory. For our purpose, three representative pulsars J0157+6212, J1743-3150 and J1857+0526 are chosen and their rotation and inclination angle evolutions are further investigated. In the future, radio and X-ray polarimetric observations will provide more information about the inclination angles of pulsars, which could help us understand the origin of the variations in $χかい$ of pulsars and shed light on the range of possibilities of pulsar magnetic field configuration. A continuous study of the pulsar inclination angle will provide an important window into additional physical processes at work in the young and highly magnetized pulsars. △ Less

Submitted 23 April, 2024; v1 submitted 22 April, 2024; originally announced April 2024.

Comments: 9 pages, 4 figures, published in AN

Journal ref: 2023, Astron. Nachr., e20220111

A short review of the pulsar magnetic inclination angles (II)

Authors: Biao-Peng Li, Wen-Qi Ma, Zhi-Fu Gao

Abstract: The pulsar magnetic inclination angle is a key parameter for pulsar physics. It influences the observable properties of pulsars, such as the pulse beam width, braking index, polarisation, and emission geometry. In this study, we give a brief overview of the current state of knowledge and research on this parameter and its implications for the internal physics of pulsars. We use the observed pulsar… ▽ More The pulsar magnetic inclination angle is a key parameter for pulsar physics. It influences the observable properties of pulsars, such as the pulse beam width, braking index, polarisation, and emission geometry. In this study, we give a brief overview of the current state of knowledge and research on this parameter and its implications for the internal physics of pulsars. We use the observed pulsar data of magnetic inclination angle and braking index to constrain the star's number of precession cycles, $ξくしー$, which reflects the interaction between superfluid neutrons and other particles inside a neutron star\,(NS). We apply the method proposed by Cheng et al. (2019) to analyse the data of PSR J2013+3845 and obtain the constraints for $ξくしー$ ranging from $2.393\times 10^{5}$ to $1.268\times10^{6}$. And further analysis suggests that the internal magnetic field structure of PSR J2013+3845 is likely dominated by toroidal component. This study may help us understand the process of internal viscous dissipation and the related evolution of the inclination angles of pulsars, and may have important implications for the study of continuous gravitational wave emissions from NS. △ Less

Submitted 23 April, 2024; v1 submitted 21 April, 2024; originally announced April 2024.

Comments: 7 pages, 2 figures, published in AN

Journal ref: Astron.Nachr., e20230167 (2023)

Non-linear power spectrum and forecasts for Generalized Cubic Covariant Galileon

Authors: Luís Atayde, Noemi Frusciante, Benjamin Bose, Santiago Casas, Baojiu Li

Abstract: To fully exploit the data from next generation surveys, we need an accurate modelling of the matter power spectrum up to non-linear scales. Therefore in this work we present the halo model reaction framework for the Generalized Cubic Covariant Galileon (GCCG) model, a modified gravity model within the Horndeski class of theories which extends the cubic covariant Galileon (G3) by including power la… ▽ More To fully exploit the data from next generation surveys, we need an accurate modelling of the matter power spectrum up to non-linear scales. Therefore in this work we present the halo model reaction framework for the Generalized Cubic Covariant Galileon (GCCG) model, a modified gravity model within the Horndeski class of theories which extends the cubic covariant Galileon (G3) by including power laws of the derivatives of the scalar field in the K-essence and cubic terms. We modify the publicly available software ReACT for the GCCG in order to obtain an accurate prediction of the non-linear power spectrum. In the limit of the G3 model we compare the modified ReACT code to $N$-body simulations and we find agreement within 5\% for a wide range of scales and redshifts. We then study the relevant effects of the modifications introduced by the GCCG on the non-linear matter power spectrum. Finally, we provide forecasts from spectroscopic and photometric primary probes by next generation surveys using a Fisher matrix method. We show that future data will be able to constrain at 1$σしぐま$ the two additional parameters of the model at the percent level and that considering non-linear corrections to the matter power spectrum beyond the linear regime is crucial to obtain this result. △ Less

Submitted 17 April, 2024; originally announced April 2024.

Comments: 17 pages, 7 figures, code available at https://github.com/nebblu/ACTio-ReACTio

Horizontally and vertically polarized kink oscillations in curved solar coronal loops

Authors: Mingzhe Guo, Tom Van Doorsselaere, Bo Li, Marcel Goossens

Abstract: Kink oscillations are frequently observed in coronal loops. This work aims to numerically clarify the influence of loop curvature on horizontally and vertically polarized kink oscillations. Working within the framework of ideal MHD, we conduct 3D simulations of axial fundamental kink oscillations in curved density-enhanced loops embedded in a potential magnetic field. Both horizontal and vertical… ▽ More Kink oscillations are frequently observed in coronal loops. This work aims to numerically clarify the influence of loop curvature on horizontally and vertically polarized kink oscillations. Working within the framework of ideal MHD, we conduct 3D simulations of axial fundamental kink oscillations in curved density-enhanced loops embedded in a potential magnetic field. Both horizontal and vertical polarizations are examined, and their oscillation frequencies are compared with WKB expectations. We discriminate two different density specifications. In the first (dubbed"uniform-density"), the density is axially uniform and varies continuously in the transverse direction toward a uniform ambient corona. Some further stratification is implemented in the second specification (dubbed"stratified"), allowing us to address the effect of evanescent barriers. Examining the oscillating profiles of the initially perturbed uniform-density loops, we found that the frequencies for both polarizations deviate from the WKB expectation by $\sim 10\%$. In the stratified loop, however, the frequency of the horizontal polarization deviates to a larger extent ($\sim 25\%$). We illustrate the lateral leakage of kink modes through wave tunnelling in 3D simulations, for the first time. Despite this, in both loops, the damping time-to-period ratios are similar and close to the analytical predictions for straight configurations under the thin-tube-thin-boundary (TTTB) assumption. The WKB expectation for straight configurations can reasonably describe the eigenfrequency of kink oscillations only in loops without an asymmetrical cross-loop density profile perpendicular to the oscillating direction. Lateral leakage via wave tunnelling is found to be less efficient than resonant absorption, meaning that the latter remains a robust damping mechanism for kink motions even when loop curvature is included. △ Less

Submitted 8 April, 2024; originally announced April 2024.

Comments: accepted for publication in A&A

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.

Impact of The Newly Revised Gravitational Redshift of X-ray Burster GS 1826-24 on The Equation of State of Supradense Neutron-Rich Matter

Authors: Wen-Jie Xie, Bao-An Li, Nai-Bo Zhang

Abstract: Thanks to the recent advancement in producing rare isotopes and measuring their masses with unprecedented precision, the updated nuclear masses around the waiting-point nucleus $^{64}$Ge in the rapid-proton capture process have led to a significant revision of the surface gravitational redshift of the neutron star (NS) in GS 1826-24 by re-fitting its X-ray burst light curve ({\it X. Zhou et al., N… ▽ More Thanks to the recent advancement in producing rare isotopes and measuring their masses with unprecedented precision, the updated nuclear masses around the waiting-point nucleus $^{64}$Ge in the rapid-proton capture process have led to a significant revision of the surface gravitational redshift of the neutron star (NS) in GS 1826-24 by re-fitting its X-ray burst light curve ({\it X. Zhou et al., Nature Physics {\bf 19}, 1091 (2023)}) using Modules for Experiments in Stellar Astrophysics (MESA). The resulting NS compactness $ξくしー$ is between 0.244 and 0.342 at 95\% confidence level and its upper boundary is significantly smaller than the maximum $ξくしー$ previously known. Incorporating this new data within a comprehensive Bayesian statistical framework, we investigate its impact on the Equation of State (EOS) of supradense neutron-rich matter and the required spin frequency for GW190814's minor $m_2$ with mass $2.59\pm 0.05$M$_{\odot}$ to be a rotationally stable pulsar. We found that the EOS of high-density symmetric nuclear matter (SNM) has to be softened significantly while the symmetry energy at supersaturation densities stiffened compared to our prior knowledge from earlier analyses using data from both astrophysical observations and terrestrial nuclear experiments. In particular, the skewness $J_0$ characterizing the stiffness of high-density SNM decreases significantly, while the slope $L$, curvature $K_{\rm{sym}}$, and skewness $J_{\rm{sym}}$ of nuclear symmetry energy all increase appreciably compared to their fiducial values. We also found that the most probable spin rate for the $m_2$ to be a stable pulsar is very close to its mass-shedding limit once the revised redshift data from GS 1826-24 is considered, making the $m_2$ unlikely the most massive NS observed so far. △ Less

Submitted 15 July, 2024; v1 submitted 2 April, 2024; originally announced April 2024.

Comments: Additional results and discussions added. Phys. Rev. D in press

Journal ref: Phys. Rev. D 110, 043025 (2024)

The FLAMINGO project: the coupling between baryonic feedback and cosmology in light of the $S_8$ tension

Authors: Willem Elbers, Carlos S. Frenk, Adrian Jenkins, Baojiu Li, John C. Helly, Roi Kugel, Matthieu Schaller, Joop Schaye, Joey Braspenning, Juliana Kwan, Ian G. McCarthy, Jaime Salcido, Marcel P. van Daalen, Bert Vandenbroucke, Silvia Pascoli

Abstract: Large-scale structure surveys have reported measurements of the density of matter, $Ωおめが_\mathrm{m}$, and the amplitude of clustering, $σしぐま_8$, that are in tension with the values inferred from observations of the cosmic microwave background. While this may be a sign of new physics that slows the growth of structure at late times, strong astrophysical feedback processes could also be responsible. In th… ▽ More Large-scale structure surveys have reported measurements of the density of matter, $Ωおめが_\mathrm{m}$, and the amplitude of clustering, $σしぐま_8$, that are in tension with the values inferred from observations of the cosmic microwave background. While this may be a sign of new physics that slows the growth of structure at late times, strong astrophysical feedback processes could also be responsible. In this work, we argue that astrophysical processes are not independent of cosmology and that their coupling naturally leads to stronger baryonic feedback in cosmological models with suppressed structure formation or when combined with a mechanism that removes dark matter from halos. We illustrate this with two well-motivated extensions of the Standard Model known to suppress structure formation: massive neutrinos and decaying dark matter. Our results, based on the FLAMINGO suite of hydrodynamical simulations, show that the combined effect of baryonic and non-baryonic suppression mechanisms is greater than the sum of its parts, particularly for decaying dark matter. We also show that the dependence of baryonic feedback on cosmology can be modelled as a function of the ratio $f_\mathrm{b}/c^2_\mathrm{v}\sim f_\mathrm{b}/(Ωおめが_\mathrm{m}σしぐま_8)^{1/4}$ of the universal baryon fraction, $f_\mathrm{b}$, to a velocity-based definition of halo concentration, $c^2_\mathrm{v}$, giving an accurate fitting formula for the baryonic suppression of the matter power spectrum. Although the combination of baryonic and non-baryonic suppression mechanisms can resolve the tension, the models with neutrinos and decaying dark matter are challenged by constraints on the expansion history. △ Less

Submitted 19 March, 2024; originally announced March 2024.

Comments: 16 pages, 9 figures, submitted to MNRAS

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)

AT2023lli: A Tidal Disruption Event with Prominent Optical Early Bump and Delayed Episodic X-ray Emission

Authors: Shifeng Huang, Ning Jiang, Jiazheng Zhu, Yibo Wang, Tinggui Wang, Shan-Qin Wang, Wen-Pei Gan, En-Wei Liang, Yu-Jing Qin, Zheyu Lin, Lin-Na Xu, Min-Xuan Cai, Ji-An Jiang, Xu Kong, Jiaxun Li, Long Li, Jian-Guo Wang, Ze-Lin Xu, Yongquan Xue, Ye-Fei Yuan, Jingquan Cheng, Lulu Fan, Jie Gao, Lei Hu, Weida Hu , et al. (20 additional authors not shown)

Abstract: High-cadence, multiwavelength observations have continuously revealed the diversity of tidal disruption events (TDEs), thus greatly advancing our knowledge and understanding of TDEs. In this work, we conducted an intensive optical-UV and X-ray follow-up campaign of TDE AT2023lli, and found a remarkable month-long bump in its UV/optical light curve nearly two months prior to maximum brightness. The… ▽ More High-cadence, multiwavelength observations have continuously revealed the diversity of tidal disruption events (TDEs), thus greatly advancing our knowledge and understanding of TDEs. In this work, we conducted an intensive optical-UV and X-ray follow-up campaign of TDE AT2023lli, and found a remarkable month-long bump in its UV/optical light curve nearly two months prior to maximum brightness. The bump represents the longest separation time from the main peak among known TDEs to date. The main UV/optical outburst declines as $t^{-4.10}$, making it one of the fastest decaying optically selected TDEs. Furthermore, we detected sporadic X-ray emission 30 days after the UV/optical peak, accompanied by a reduction in the period of inactivity. It is proposed that the UV/optical bump could be caused by the self-intersection of the stream debris, whereas the primary peak is generated by the reprocessed emission of the accretion process. In addition, our results suggest that episodic X-ray radiation during the initial phase of decline may be due to the patched obscurer surrounding the accretion disk, a phenomenon associated with the inhomogeneous reprocessing process. The double TDE scenario, in which two stars are disrupted in sequence, is also a possible explanation for producing the observed early bump and main peak. We anticipate that the multicolor light curves of TDEs, especially in the very early stages, and the underlying physics can be better understood in the near future with the assistance of dedicated surveys such as the deep high-cadence survey of the 2.5-meter Wide Field Survey Telescope (WFST). △ Less

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

Comments: 14 pages, 8 figures,accepted for publication by ApJL

A Search for Radio Pulsars in Supernova Remnants Using FAST with One Pulsar Discovered

Authors: Zhen Zhang, Wen-Ming Yan, Jian-Ping Yuan, Na Wang, Jun-Tao Bai, Zhi-Gang Wen, Bao-Da Li, Jin-Tao Xie, De Zhao, Yu-Bin Wang, Nan-Nan Zhai

Abstract: We report on the results of a search for radio pulsars in five supernova remnants (SNRs) with FAST. The observations were made using the 19-beam receiver in the Snapshot mode. The integration time for each pointing is 10 min. We discovered a new pulsar PSR J1845$-$0306 which has a spin period of 983.6 ms and a dispersion measure of 444.6$\pm$2.0 cm$^{-3}$ pc in observations of SNR G29.6+0.1. To ju… ▽ More We report on the results of a search for radio pulsars in five supernova remnants (SNRs) with FAST. The observations were made using the 19-beam receiver in the Snapshot mode. The integration time for each pointing is 10 min. We discovered a new pulsar PSR J1845$-$0306 which has a spin period of 983.6 ms and a dispersion measure of 444.6$\pm$2.0 cm$^{-3}$ pc in observations of SNR G29.6+0.1. To judge the association between the pulsar and the SNR, further verification is needed. We also re-detected some known pulsars in the data from SNRs G29.6+0.1 and G29.7$-$0.3. No pulsars were detected in observations of other three SNRs. △ Less

Submitted 27 February, 2024; originally announced February 2024.

Comments: 6 pages, 2 figures, 2 tables published in CPL

Journal ref: Chin. Phys. Lett. 2024, 41 (2): 029701 February 2024

Horizontally Polarized Kink Oscillations Supported by Solar Coronal Loops in an Asymmetric Environment

Authors: Mijie Shi, Bo Li, Shengju Yuan

Abstract: Kink oscillations are ubiquitously observed in solar coronal loops, their understanding being crucial in the contexts of coronal seismology and atmospheric heating. We study kink modes supported by a straight coronal loop embeded in an asymmetric environment using three-dimensional magnetohydrodynamic (MHD) simulations. We implement the asymmetric effect by setting different exterior densities bel… ▽ More Kink oscillations are ubiquitously observed in solar coronal loops, their understanding being crucial in the contexts of coronal seismology and atmospheric heating. We study kink modes supported by a straight coronal loop embeded in an asymmetric environment using three-dimensional magnetohydrodynamic (MHD) simulations. We implement the asymmetric effect by setting different exterior densities below and above the loop interior, and initiate the simulation using a kink-like velocity perturbation perpendicular to the loop plane, mimicking the frequently measured horizontally polarized kink modes. We find that the external velocity fields show fan blade structures propagating in the azimuthal direction as a result of the successive excitation of higher azimuthal Fourier modes. Resonant absorption and phase mixing can still occur despite an asymmetric environment, leading to the development of small scales at loop boundaries. These small scales nonetheless develop asymmetrically at the upper and lower boundaries due to the different gradients of the Alfven speed. These findings enrich our understanding of kink modes in coronal loops embedded within an asymmetric environment, providing insights helpful for future high-resolution observations. △ Less

Submitted 16 February, 2024; originally announced February 2024.

Comments: Accepted for publication in A&A

Damped kink motions in a system of two solar coronal tubes with elliptic cross-sections

Authors: Mijie Shi, Bo Li, Shaoxia Chen, Hui Yu, Mingzhe Guo

Abstract: This study is motivated by observations of coordinated transverse displacements in neighboring solar active region loops, addressing specifically how the behavior of kink motions in straight two-tube equilibria is impacted by tube interactions and tube cross-sectional shapes.We work with linear, ideal, pressureless magnetohydrodynamics. Axially standing kink motions are examined as an initial valu… ▽ More This study is motivated by observations of coordinated transverse displacements in neighboring solar active region loops, addressing specifically how the behavior of kink motions in straight two-tube equilibria is impacted by tube interactions and tube cross-sectional shapes.We work with linear, ideal, pressureless magnetohydrodynamics. Axially standing kink motions are examined as an initial value problem for transversely structured equilibria involving two identical, field-aligned, density-enhanced tubes with elliptic cross-sections (elliptic tubes). Continuously nonuniform layers are implemented around both tube boundaries. We numerically follow the system response to external velocity drivers, largely focusing on the quasi-mode stage of internal flows to derive the pertinent periods and damping times. The periods and damping times we derive for two-circular-tube setups justify available modal results found with the T-matrix approach. Regardless of cross-sectional shapes, our nonuniform layers feature the development of small-scale shears and energy accumulation around Alfven resonances, indicative of resonant absorption and phase-mixing. As with two-circular-tube systems, our configurational symmetries make it still possible to classify lower-order kink motions by the polarization and symmetric properties of the internal flows; hence such mode labels as $S_x$ and $A_x$. However, the periods and damping times for two-elliptic-tube setups further depend on cross-sectional aspect ratios, with $A_x$ motions occasionally damped less rapidly than $S_x$ motions. We find uncertainties up to $\sim 20\%$ ($\sim 50\%$) for the axial Alfven time (the inhomogeneity lengthscale) if the periods (damping times) computed for two-elliptic-tube setups are seismologically inverted with canonical theories for isolated circular tubes. △ Less

Submitted 3 March, 2024; v1 submitted 23 January, 2024; originally announced January 2024.

Comments: Accepted for publication in A&A

Research on the knee region of cosmic ray by using a novel type of electron-neutron detector array

Authors: Bing-Bing Li, Xin-Hua Ma, Shu-Wang Cui, Hao-Kun Chen, Tian-Lu Chen, Danzengluobu, Wei Gao, Hai-Bing Hu, Denis Kuleshov, Kirill Kurinov, Hu Liu, Mao-Yuan Liu, Ye Liu, Da-Yu Peng, Yao-Hui Qi, Oleg Shchegolev, Yuri Stenkin, Li-Qiao Yin, Heng-Yu Zhang, Liang-Wei Zhang

Abstract: By accurately measuring composition and energy spectrum of cosmic ray, the origin problem of so called "keen" region (energy > 1 PeV) can be solved. However, up to the present, the results of the spectrum in the knee region obtained by several previous experiments have shown obvious differences, so they cannot give effective evidence for judging the theoretical models on the origin of the knee. Re… ▽ More By accurately measuring composition and energy spectrum of cosmic ray, the origin problem of so called "keen" region (energy > 1 PeV) can be solved. However, up to the present, the results of the spectrum in the knee region obtained by several previous experiments have shown obvious differences, so they cannot give effective evidence for judging the theoretical models on the origin of the knee. Recently, the Large High Altitude Air Shower Observatory (LHAASO) has reported several major breakthroughs and important results in astro-particle physics field. Relying on its advantages of wide-sky survey, high altitude location and large area detector arrays, the research content of LHAASO experiment mainly includes ultra high-energy gamma-ray astronomy, measurement of cosmic ray spectra in the knee region, searching for dark matter and new phenomena of particle physics at higher energy. The electron and Thermal Neutron detector (EN-Detector) is a new scintillator detector which applies thermal neutron detection technology to measure cosmic ray extensive air shower (EAS). This technology is an extension of LHAASO. The EN-Detector Array (ENDA) can highly efficiently measure thermal neutrons generated by secondary hadrons so called "skeleton" of EAS. In this paper, we perform the optimization of ENDA configuration, and obtain expectations on the ENDA results, including thermal neutron distribution, trigger efficiency and capability of cosmic ray composition separation. The obtained real data results are consistent with those by the Monte Carlo simulation. △ Less

Submitted 23 January, 2024; originally announced January 2024.

fkPT: Constraining scale-dependent modified gravity with the full-shape galaxy power spectrum

Authors: Mario A. Rodriguez-Meza, Alejandro Aviles, Hernan E. Noriega, Cheng-Zong Ruan, Baojiu Li, Mariana Vargas-Magaña, Jorge L. Cervantes-Cota

Abstract: Modified gravity models with scale-dependent linear growth typically exhibit an enhancement in the power spectrum beyond a certain scale. The conventional methods for extracting cosmological information usually involve inferring modified gravity effects via Redshift Space Distortions (RSD), particularly through the time evolution of $fσしぐま_8$. However, classical galaxy RSD clustering analyses encount… ▽ More Modified gravity models with scale-dependent linear growth typically exhibit an enhancement in the power spectrum beyond a certain scale. The conventional methods for extracting cosmological information usually involve inferring modified gravity effects via Redshift Space Distortions (RSD), particularly through the time evolution of $fσしぐま_8$. However, classical galaxy RSD clustering analyses encounter difficulties in accurately capturing the spectrum's enhanced power, which is better obtained from the broad-band power spectrum. In this sense, full-shape analyses aim to consider survey data using comprehensive and precise models of the whole power spectrum. Yet, a major challenge in this approach is the slow computation of non-linear loop integrals for scale-dependent modified gravity, precluding the estimation of cosmological parameters using Markov Chain Monte Carlo methods. Based on recent studies, in this work we develop a perturbation theory tailored for Modified Gravity, or analogous scenarios introducing additional scales, such as in the presence of massive neutrinos. Our approach only needs the calculation of the scale-dependent growth rate $f(k,t)$ and the limit of the perturbative kernels at large scales. We called this approximate technique as fk-Perturbation Theory and implemented it into the code fkpt, capable of computing the redshift space galaxy power spectrum in a fraction of a second. We validate our modeling and code with the $f(R)$ theory MG-GLAM and General Relativity NSeries sets of simulations. The code is available at https://github.com/alejandroaviles/fkpt △ Less

Submitted 4 May, 2024; v1 submitted 16 December, 2023; originally announced December 2023.

Comments: 34 pages, 9 figures. v2

Journal ref: JCAP03(2024)049

The bright black hole X-ray binary 4U 1543--47 during 2021 outburst: a thick accretion disk inflated by high luminosity

Authors: S. J. Zhao, L. Tao, P. P. Li, R. Soria, H. Feng, Y. X. Zhang, R. C. Ma, W. D. Zhang, E. L. Qiao, Q. Q. Yin, S. N. Zhang, L. Zhang, Q. C. Bu, X. Ma, Y. Huang, M. Y. Ge, X. B. Li, Q. C. Zhao, J. Q. Peng, Y. X. Xiao

Abstract: The black hole X-ray binary source 4U 1543--47 experienced a super-Eddington outburst in 2021, reaching a peak flux of up to $\sim1.96\times10^{-7}\rm erg\ \rm cm^{-2}\ \rm s^{-1}$ ($\sim 8.2$ Crab) in the 2--10\,keV band. Soon after the outburst began, it rapidly transitioned into the soft state. Our goal is to understand how the accretion disk structure deviates from a standard thin disk when th… ▽ More The black hole X-ray binary source 4U 1543--47 experienced a super-Eddington outburst in 2021, reaching a peak flux of up to $\sim1.96\times10^{-7}\rm erg\ \rm cm^{-2}\ \rm s^{-1}$ ($\sim 8.2$ Crab) in the 2--10\,keV band. Soon after the outburst began, it rapidly transitioned into the soft state. Our goal is to understand how the accretion disk structure deviates from a standard thin disk when the accretion rate is near Eddington. To do so, we analyzed spectra obtained from quasi-simultaneous observations conducted by the Hard X-ray Modulation Telescope (Insight-HXMT), the Nuclear Spectroscopic Telescope Array (NuSTAR), and the Neil Gehrels Swift Observatory (Swift). These spectra are well-fitted by a model comprising a disk, a weak corona, and a reflection component. We suggest that the reflection component is caused by disk self-irradiation, that is by photons emitted from the inner disk which return to the accretion disk surface, as their trajectories are bent by the strong gravity field. In this scenario, the best-fitting parameters imply that the reflected flux represents more than half of the total flux. Using general relativistic ray-tracing simulations, we show that this scenario is viable when the disk becomes geometrically thick, with a funnel-like shape, as the accretion rate is near or above the Eddington limit. In the specific case of 4U 1543--47, an angle $\gtrsim$ 45 deg between the disk surface and the equatorial plane can explain the required amount of self-irradiation. △ Less

Submitted 15 December, 2023; originally announced December 2023.

Comments: Accepted for publication in Astronomy and Astrophysics. 15 pages, 4 tables, 12 figures

Journal ref: A&A 685, A42 (2024)

Three Pulsars Discovered in Globular Cluster M15 (NGC 7078) with FAST

Authors: Yuxiao Wu, Zhichen Pan, Lei Qian, Scott Ransom, BoJun Wang, Zhen Yan, Jintao Luo, Liyun Zhang, Minghui Li, Dejiang Yin, Baoda Li, Yifeng Li, Yinfeng Dai, Yaowei Li, Xinnan Zhang, Tong Liu, Yu Pan

Abstract: We present the discovery of three pulsars in Globular Cluster M15 (NGC 7078) by the Five-hundred-meter Aperture Spherical radio Telescope (FAST). In the three pulsars, PSR~J2129+1210J (M15J) is a millisecond pulsar with a spinning period of 11.84 ms and a dispersion measure of 66.68 pc cm$^{-3}$. Both PSR~J2129+1210K and L (M15K and L) are long period pulsars with spinning periods of 1928 ms and 3… ▽ More We present the discovery of three pulsars in Globular Cluster M15 (NGC 7078) by the Five-hundred-meter Aperture Spherical radio Telescope (FAST). In the three pulsars, PSR~J2129+1210J (M15J) is a millisecond pulsar with a spinning period of 11.84 ms and a dispersion measure of 66.68 pc cm$^{-3}$. Both PSR~J2129+1210K and L (M15K and L) are long period pulsars with spinning periods of 1928 ms and 3961 ms , respectively, while M15L is the GC pulsar with the longest spinning period till now. The discoveries of M15K and L support the theory that core-collapsed Globular Clusters may contain partially recycled long period pulsars. With the same dataset, the timing solutions of M15A to H were updated, and the timing parameter P1 of M15F is different from the previous results, which is approximately 0.027$\times 10^{-18} ss^{-1}$ from our work and $0.032 \times 10^{-18} ss^{-1}$ from Anderson's\citep{anderson-1993}. As predicted by Rodolfi et al. , the luminosity of M15C kept decreasing and the latest detection in our dataset is on December 20$^{\rm th}$, 2022. We also detected M15I for one more time. The different barycentric spin periods indicate that this pulsar should locate in a binary system, manifesting itself as the exceptional one in such a core-collapsing GC. △ Less

Submitted 10 December, 2023; originally announced December 2023.

Comments: 10 pages, 4 figures, 2 tables, submitted to ApJ Letter

Estimation of line-of-sight velocities of individual galaxies using neural networks I. Modelling redshift-space distortions at large scales

Authors: Hongxiang Chen, Jie Wang, Tianxiang Mao, Juntao Ma, Yuxi Meng, Baojiu Li, Yan-Chuan Cai, Mark Neyrinck, Bridget Falck, Alexander S. Szalay

Abstract: We present a scheme based on artificial neural networks (ANN) to estimate the line-of-sight velocities of individual galaxies from an observed redshift-space galaxy distribution. We find an estimate of the peculiar velocity at a galaxy based on galaxy counts and barycenters in shells around it. By training the network with environmental characteristics, such as the total mass and mass center withi… ▽ More We present a scheme based on artificial neural networks (ANN) to estimate the line-of-sight velocities of individual galaxies from an observed redshift-space galaxy distribution. We find an estimate of the peculiar velocity at a galaxy based on galaxy counts and barycenters in shells around it. By training the network with environmental characteristics, such as the total mass and mass center within each shell surrounding every galaxy in redshift space, our ANN model can accurately predict the line-of-sight velocity of each individual galaxy. When this velocity is used to eliminate the RSD effect, the two-point correlation function (TPCF) in real space can be recovered with an accuracy better than 1% at $s$ > 8 $h^{-1}\mathrm{Mpc}$, and 4% on all scales compared to ground truth. The real-space power spectrum can be recovered within 3% on $k$< 0.5 $\mathrm{Mpc}^{-1}h$, and less than 5% for all $k$ modes. The quadrupole moment of the TPCF or power spectrum is almost zero down to $s$ = 10 $h^{-1}\mathrm{Mpc}$ or all $k$ modes, indicating an effective correction of the spatial anisotropy caused by the RSD effect. We demonstrate that on large scales, without additional training with new data, our network is adaptable to different galaxy formation models, different cosmological models, and mock galaxy samples at high redshifts and high biases, achieving less than 10% error for scales greater than 15 $h^{-1}\mathrm{Mpc}$. As it is sensitive to large-scale densities, it does not manage to remove Fingers of God in large clusters, but works remarkably well at recovering real-space galaxy positions elsewhere. Our scheme provides a novel way to predict the peculiar velocity of individual galaxies, to eliminate the RSD effect directly in future large galaxy surveys, and to reconstruct the 3-D cosmic velocity field accurately. △ Less

Submitted 8 July, 2024; v1 submitted 6 December, 2023; originally announced December 2023.

Comments: 15 pages,15 figures, accepted for publication in MNRAS

The Uchuu-GLAM BOSS and eBOSS LRG lightcones: Exploring clustering and covariance errors

Authors: Julia Ereza, Francisco Prada, Anatoly Klypin, Tomoaki Ishiyama, Alex Smith, Carlton M. Baugh, Baojiu Li, César Hernández-Aguayo, José Ruedas

Abstract: This study investigates the clustering and bias of Luminous Red Galaxies (LRG) in the BOSS-LOWZ, -CMASS, -COMB, and eBOSS samples, using two types of simulated lightcones: (i) high-fidelity lightcones from Uchuu $N$-body simulation, employing SHAM technique to assign LRG to (sub)halos, and (ii) 16000 covariance lightcones from GLAM-Uchuu $N$-body simulations, including LRG using HOD data from Uchu… ▽ More This study investigates the clustering and bias of Luminous Red Galaxies (LRG) in the BOSS-LOWZ, -CMASS, -COMB, and eBOSS samples, using two types of simulated lightcones: (i) high-fidelity lightcones from Uchuu $N$-body simulation, employing SHAM technique to assign LRG to (sub)halos, and (ii) 16000 covariance lightcones from GLAM-Uchuu $N$-body simulations, including LRG using HOD data from Uchuu. Our results indicate that Uchuu and GLAM lightcones closely replicate BOSS/eBOSS data, reproducing correlation function and power spectrum across scales from redshifts $0.2$ to $1.0$, from $2$ to $150\,h^{-1}\mathrm{Mpc}$ in configuration space, from $0.005$ to $0.7\,h\mathrm{Mpc}^{-1}$ in Fourier space, and across different LRG stellar masses. Furthermore, comparing with existing MD-Patchy and EZmock BOSS/eBOSS lightcones based on approximate methods, our GLAM-Uchuu lightcones provide more precise clustering estimates. We identify significant deviations from observations within $20\,h^{-1}\mathrm{Mpc}$ scales in MD-Patchy and EZmock, with our covariance matrices indicating that these methods underestimate errors by between $10\%$ and $60\%$. Lastly, we explore the impact of cosmology on galaxy clustering. Our findings suggest that, given the current level of uncertainties in BOSS/eBOSS data, distinguishing models with and without massive neutrino effects on LSS is challenging. This paper highlights the Uchuu and GLAM-Uchuu simulations' robustness in verifying the accuracy of Planck cosmological parameters, providing a strong foundation for enhancing lightcone construction in future LSS surveys. We also demonstrate that generating thousands of galaxy lightcones is feasible using $N$-body simulations with adequate mass and force resolution. △ Less

Submitted 26 June, 2024; v1 submitted 24 November, 2023; originally announced November 2023.

Journal ref: Monthly Notices of the Royal Astronomical Society, stae1543 (2024)

Strong Gravity Extruding Peaks in Speed of Sound Profiles of Massive Neutron Stars

Authors: Bao-Jun Cai, Bao-An Li

Abstract: The speed of sound squared (SSS) $s^2$ in massive neutron stars (NSs) characterizes not only the stiffness of supradense neutron-rich matter within but also equivalently properties of the curved geometry due to the strong-field gravity and matter-geometry coupling. A peaked density or radius profile of $s^2$ has been predicted for massive NSs using various NS Equation of State (EOS) models. Howeve… ▽ More The speed of sound squared (SSS) $s^2$ in massive neutron stars (NSs) characterizes not only the stiffness of supradense neutron-rich matter within but also equivalently properties of the curved geometry due to the strong-field gravity and matter-geometry coupling. A peaked density or radius profile of $s^2$ has been predicted for massive NSs using various NS Equation of State (EOS) models. However, the nature, cause, location and size of the peak in $s^2$ profiles are still very EOS model dependent. In this work, we investigate systematically $s^2$ profiles in massive NSs in a new approach that is independent of the nuclear EOS model and without any presumption about the NS structure and/or composition. In terms of the small quantities (reduced radius, the energy density and pressure scaled by their central values), we perform double-element perturbative expansions in solving perturbatively the scaled Tolman--Oppenheimer--Volkoff (TOV) equations and analyzing $s^2$ profiles from the Newtonian limit to the general relativistic (GR) case. The GR term in the TOV equations plays a twofold role: it compresses NS matter and modifies the pressure/energy density ratio from small values in Newtonian stars showing no $s^2$ peak to large ones for massive NSs possessing a peak in their $s^2$ profiles, and eventually takes away the peak in extremely compact/massive NSs approaching the causality limit. {In particular, the peaked behavior in $s^2$ is expected to emerge near the center of massive NSs like PSR J0740+6620, while a sharp phase transition is unlikely to occur there.} These features revealed from our analyses are universal as they are intrinsic properties of the GR stellar structure equations independent of the still very uncertain EOS of supradense neutron-rich matter in NSs. △ Less

Submitted 15 March, 2024; v1 submitted 21 November, 2023; originally announced November 2023.

Comments: Added/revised some discussions. Phys. Rev. D in press

Journal ref: Phys. Rev. D 109, 083015 (2024)

The first Ka-band (26.1-35 GHz) blind line survey towards Orion KL

Authors: Xunchuan Liu, Tie Liu, Zhiqiang Shen, Sheng-Li Qin, Qiuyi Luo, Yan Gong, Yu Cheng, Christian Henkel, Qilao Gu, Fengyao Zhu, Tianwei Zhang, Rongbing Zhao, Yajun Wu, Bin Li, Juan Li, Zhang Zhao, Jinqing Wang, Weiye Zhong, Qinghui Liu, Bo Xia, Li Fu, Zhen Yan, Chao Zhang, Lingling Wang, Qian Ye , et al. (9 additional authors not shown)

Abstract: We conducted a Ka-band (26.1--35 GHz) line survey towards Orion KL using the TianMa 65-m Radio Telescope (TMRT). It is the first blind line survey in the Ka band, and achieves a sensitivity of mK level (1--3 mK at a spectral resolution of $\sim$1 km s$^{-1}$). In total, 592 Gaussian features are extracted. Among them, 257 radio recombination lines (RRLs) are identified. The maximum $Δでるたn$ of RRLs of… ▽ More We conducted a Ka-band (26.1--35 GHz) line survey towards Orion KL using the TianMa 65-m Radio Telescope (TMRT). It is the first blind line survey in the Ka band, and achieves a sensitivity of mK level (1--3 mK at a spectral resolution of $\sim$1 km s$^{-1}$). In total, 592 Gaussian features are extracted. Among them, 257 radio recombination lines (RRLs) are identified. The maximum $Δでるたn$ of RRLs of H, He and C are 20, 15, and 5, respectively. Through stacking, we have detected the $βべーた$ lines of ion RRLs (RRLs of C$^+$ with possible contribution of other ions like O$^+$) for the first time, and tentative signal of the $γがんま$ lines of ion RRLs can also be seen on the stacked spectrum. Besides, 318 other line features were assigned to 37 molecular species, and ten of these species were not detected in the Q-band survey of TMRT. The vibrationally excited states of nine species were also detected. Emission of most species can be modeled under LTE. A number of transitions of E-CH3OH ($J_2-J_1$) display maser effects, which are confirmed by our modeling, and besides the bumping peak at $J\sim 6$ there is another peak at $J\sim 13$. Methylcyanoacetylene (CH$_3$C$_3$N) is detected in Orion KL for the first time. This work emphasizes that the Ka band, which was long-ignored for spectral line surveys, is very useful for surveying RRLs and molecular lines simultaneously. △ Less

Submitted 20 November, 2023; originally announced November 2023.

Comments: accepted by ApJS

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

Observation of GRB 221009A early afterglow in X/$γがんま$-ray energy band

Authors: Chao Zheng, Yan-Qiu Zhang, Shao-Lin Xiong, Cheng-Kui Li, He Gao, Wang-Chen Xue, Jia-Cong Liu, Chen-Wei Wang, Wen-Jun Tan, Wen-Xi Peng, Zheng-Hua An, Ce Cai, Ming-Yu Ge, Dong-Ya Guo, Yue Huang, Bing Li, Ti-Pei Li, Xiao-Bo Li, Xin-Qiao Li, Xu-Fang Li, Jin-Yuan Liao, Cong-Zhan Liu, Fang-Jun Lu, Xiang Ma, Rui Qiao , et al. (23 additional authors not shown)

Abstract: The early afterglow of a Gamma-ray burst (GRB) can provide critical information on the jet and progenitor of the GRB. The extreme brightness of GRB 221009A allows us to probe its early afterglow in unprecedented detail. In this letter, we report comprehensive observation results of the early afterglow of GRB 221009A (from $T_0$+660 s to $T_0$+1860 s, where $T_0$ is the \textit{Insight}-HXMT/HE tri… ▽ More The early afterglow of a Gamma-ray burst (GRB) can provide critical information on the jet and progenitor of the GRB. The extreme brightness of GRB 221009A allows us to probe its early afterglow in unprecedented detail. In this letter, we report comprehensive observation results of the early afterglow of GRB 221009A (from $T_0$+660 s to $T_0$+1860 s, where $T_0$ is the \textit{Insight}-HXMT/HE trigger time) in X/$γがんま$-ray energy band (from 20 keV to 20 MeV) by \textit{Insight}-HXMT/HE, GECAM-C and \textit{Fermi}/GBM. We find that the spectrum of the early afterglow in 20 keV-20 MeV could be well described by a cutoff power-law with an extra power-law which dominates the low and high energy bands respectively. The cutoff power-law $E_{\rm peak}$ is $\sim$ 30 keV and the power-law photon index is $\sim$ 1.8 throughout the early afterglow phase. By fitting the light curves in different energy bands, we find that a significant achromatic break (from keV to TeV) is required at $T_0$ + 1246$^{+27}_{-26}$ s (i.e. 1021 s since the afterglow starting time $T_{\rm AG}$=$T_0$+225 s), providing compelling evidence of a jet break. Interestingly, both the pre-break and post-break decay slopes vary with energy, and these two slopes become closer in the lower energy band, making the break less identifiable. Intriguingly, the spectrum of the early afterglow experienced a slight hardening before the break and a softening after the break. These results provide new insights into the understanding of this remarkable GRB. △ Less

Submitted 19 January, 2024; v1 submitted 16 October, 2023; originally announced October 2023.

Comments: Accepted for publication in ApJ Letters on 19-Jan-2024, 11 pages, 7 figures and 2 tables

Searching for X-Ray Counterparts of Degree Wide TeV Halos Around Middle-Aged Pulsars with SRG/eROSITA

Authors: A. Khokhriakova, W. Becker, G. Ponti, M. Sasaki, B. Li, R. -Y. Liu

Abstract: Context. Extended gamma-ray TeV emission (TeV halos) around middle-aged pulsars has been detected. A proposed model to explain these TeV halos is that electrons from a degree-wide Pulsar Wind Nebula (PWN) get up-scattered by cosmic microwave background photons through inverse Compton processes. However, no X-ray degree-wide faint diffuse PWNe have been found around these middle-aged pulsars in pre… ▽ More Context. Extended gamma-ray TeV emission (TeV halos) around middle-aged pulsars has been detected. A proposed model to explain these TeV halos is that electrons from a degree-wide Pulsar Wind Nebula (PWN) get up-scattered by cosmic microwave background photons through inverse Compton processes. However, no X-ray degree-wide faint diffuse PWNe have been found around these middle-aged pulsars in previous X-ray observations. Aims. We have performed a search for degree wide PWNe around Geminga, PSR B0656+14, B0540+23, J0633+0632, and J0631+1036, using data from the first four consecutive Spectrum Roentgen Gamma/eROSITA all-sky surveys. In order to better understand the mechanisms underlying the formation of TeV halos, we investigated the magnetic field strength in the degree wide neighbourhood of those pulsars. Results. We did not detect degree-wide diffuse emission around Geminga, PSR B0656+14, B0540+23, J0633+0632, and J0631+1036, which can be attributed to being powered by the rotation-powered pulsars. Indeed, a close inspection of the data shows that the pulsars of interest are all embedded in diffuse emission from supernova remnants like the Monogem Ring or the Rosetta Nebula, while PSR B0540+23 is located ~2.5 degrees away from the bright Crab pulsar, which shines out the eROSITA point-spread function up to the position of PSR B0540+23 and thus reduced the sensitivity to search for degree wide bright diffuse X-ray emission strongly. Conclusions. Despite the non-detection of any degree-wide PWN surrounding the analysed pulsars, we set flux upper limits to provide useful information on magnetic field strength and its spatial distribution around those pulsars, providing additional constraints to the proposed theory for the formation of TeV halos around pulsars. △ Less

Submitted 16 October, 2023; originally announced October 2023.

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

Precessing jet nozzle connecting to a spinning black hole in M87

Authors: Yuzhu Cui, Kazuhiro Hada, Tomohisa Kawashima, Motoki Kino, Weikang Lin, Yosuke Mizuno, Hyunwook Ro, Mareki Honma, Kunwoo Yi, Jintao Yu, Jongho Park, Wu Jiang, Zhiqiang Shen, Evgeniya Kravchenko, Juan-Carlos Algaba, Xiaopeng Cheng, Ilje Cho, Gabriele Giovannini, Marcello Giroletti, Taehyun Jung, Ru-Sen Lu, Kotaro Niinuma, Junghwan Oh, Ken Ohsuga, Satoko Sawada-Satoh , et al. (54 additional authors not shown)

Abstract: The nearby radio galaxy M87 offers a unique opportunity to explore the connections between the central supermassive black hole and relativistic jets. Previous studies of the inner region of M87 revealed a wide opening angle for the jet originating near the black hole. The Event Horizon Telescope resolved the central radio source and found an asymmetric ring structure consistent with expectations f… ▽ More The nearby radio galaxy M87 offers a unique opportunity to explore the connections between the central supermassive black hole and relativistic jets. Previous studies of the inner region of M87 revealed a wide opening angle for the jet originating near the black hole. The Event Horizon Telescope resolved the central radio source and found an asymmetric ring structure consistent with expectations from General Relativity. With a baseline of 17 years of observations, there was a shift in the jet's transverse position, possibly arising from an eight to ten-year quasi-periodicity. However, the origin of this sideways shift remains unclear. Here we report an analysis of radio observations over 22 years that suggests a period of about 11 years in the position angle variation of the jet. We infer that we are seeing a spinning black hole that induces the Lense-Thirring precession of a misaligned accretion disk. Similar jet precession may commonly occur in other active galactic nuclei but has been challenging to detect owing to the small magnitude and long period of the variation. △ Less

Submitted 13 October, 2023; originally announced October 2023.

Comments: 41 pages, 7 figures, 7 tables

Journal ref: 2023, Nature, 621, 711-715

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

Evidence of mini-jet emission in a large emission zone from a magnetically-dominated gamma-ray burst jet

Authors: S. -X. Yi, C. -W. Wang, X. -Y. Shao, R. Moradi, H. Gao, B. Zhang, S. -L. Xiong, S. -N. Zhang, W. -J. Tan, J. -C. Liu, W. -C. Xue, Y. -Q. Zhang, C. Zheng, Y. Wang, P. Zhang, Z. -H. An, C. Cai, P. -Y. Feng, K. Gong, D. -Y. Guo, Y. Huang, B. Li, X. -B. Li, X. -Q. Li, X. -J. Liu , et al. (21 additional authors not shown)

Abstract: The second brightest GRB in history, GRB230307A provides an ideal laboratory to study the details of GRB prompt emission thanks to its extraordinarily high photon statistics and its single broad pulse overall shape characterized by an energy-dependent fast-rise-exponential-decay (FRED) profile. Here we demonstrate that its broad pulse is composed of many rapidly variable short pulses, rather than… ▽ More The second brightest GRB in history, GRB230307A provides an ideal laboratory to study the details of GRB prompt emission thanks to its extraordinarily high photon statistics and its single broad pulse overall shape characterized by an energy-dependent fast-rise-exponential-decay (FRED) profile. Here we demonstrate that its broad pulse is composed of many rapidly variable short pulses, rather than being the superposition of many short pulses on top of a slow component. Such a feature is consistent with the picture of many mini-jets due to local magnetic reconnection events in a large emission zone far from the GRB central engine, as envisaged in the internal-collision-induced magnetic reconnection and turbulence (ICMART) model, but raises a great challenge to the internal shock models that attribute all variability components to collisions among different shells. Since relativistic mini-jets demand strong magnetization in the outflow, this work provides strong evidence for a Poynting-flux-dominated jet composition of this bright GRB. △ Less

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

Comments: 7 pages and 2 figures in the main text. 27 pages and 9 figures in total

Galaxy evolution in modified gravity simulations: using passive galaxies to constrain gravity with upcoming surveys

Authors: Diego Pallero, Facundo A. Gómez, Nelson D. Padilla, Yara L. Jaffé, Carlton M. Baugh, Baojiu Li, César Hernández-Aguayo, Christian Arnold

Abstract: We present a quantitative analysis of the properties of galaxies and structures evolving in universes dominated by different modified gravitational models, including two variants of the f(R)-gravity (F) and two of the Dvali-Gabdadze-Poratti (N) braneworld model, which respectively feature the chameleon and Vainshtein screening mechanisms. Using the Simulation HYdrodynamics BeyONd Einstein (SHYBONE… ▽ More We present a quantitative analysis of the properties of galaxies and structures evolving in universes dominated by different modified gravitational models, including two variants of the f(R)-gravity (F) and two of the Dvali-Gabdadze-Poratti (N) braneworld model, which respectively feature the chameleon and Vainshtein screening mechanisms. Using the Simulation HYdrodynamics BeyONd Einstein (SHYBONE) cosmological hydrodynamical full-physics simulations suite, we study the departures in the properties of galaxies residing in different environments with respect to the standard model (GR). Using two different criteria to compare, we find that structures formed within modified gravity tend to show a denser gas density profile than their GR counterparts. Within the different modified gravity models, N1 and F5 gravity models show greater departures from the standard model, with gas density profiles $ρろー_{\rm IGM} \geq 30\%$ denser in the outskirts for the N1 model, and in the inner parts for the F5 model. Additionally, we find that haloes evolving in MG universes show, in general, larger quenched fractions than GR, reaching up to $20\%$ larger quenching fractions in F5 regardless of the stellar mass of the galaxy. With respect to the other models, F6, N1 and N5 show slightly larger quenched fractions, but no strong differences can be found. These results directly impact the colour distribution of galaxies, making them in MG models redder and older than their GR counterparts. Like GR, once the environment starts to play a role, galaxies rapidly get quenched and the differences between models vanish. △ Less

Submitted 3 October, 2023; originally announced October 2023.

Comments: 17 pages, 12 figures, submitted to MNRAS

Euclid: The search for primordial features

Authors: M. Ballardini, Y. Akrami, F. Finelli, D. Karagiannis, B. Li, Y. Li, Z. Sakr, D. Sapone, A. Achúcarro, M. Baldi, N. Bartolo, G. Cañas-Herrera, S. Casas, R. Murgia, H. A. Winther, M. Viel, A. Andrews, J. Jasche, G. Lavaux, D. K. Hazra, D. Paoletti, J. Valiviita, A. Amara, S. Andreon, N. Auricchio , et al. (104 additional authors not shown)

Abstract: Primordial features, in particular oscillatory signals, imprinted in the primordial power spectrum of density perturbations represent a clear window of opportunity for detecting new physics at high-energy scales. Future spectroscopic and photometric measurements from the $Euclid$ space mission will provide unique constraints on the primordial power spectrum, thanks to the redshift coverage and hig… ▽ More Primordial features, in particular oscillatory signals, imprinted in the primordial power spectrum of density perturbations represent a clear window of opportunity for detecting new physics at high-energy scales. Future spectroscopic and photometric measurements from the $Euclid$ space mission will provide unique constraints on the primordial power spectrum, thanks to the redshift coverage and high-accuracy measurement of nonlinear scales, thus allowing us to investigate deviations from the standard power-law primordial power spectrum. We consider two models with primordial undamped oscillations superimposed on the matter power spectrum, one linearly spaced in $k$-space the other logarithmically spaced in $k$-space. We forecast uncertainties applying a Fisher matrix method to spectroscopic galaxy clustering, weak lensing, photometric galaxy clustering, cross correlation between photometric probes, spectroscopic galaxy clustering bispectrum, CMB temperature and $E$-mode polarization, temperature-polarization cross correlation, and CMB weak lensing. We also study a nonlinear density reconstruction method to retrieve the oscillatory signals in the primordial power spectrum. We find the following percentage relative errors in the feature amplitude with $Euclid$ primary probes for the linear (logarithmic) feature model: 21% (22%) in the pessimistic settings and 18% (18%) in the optimistic settings at 68.3% confidence level (CL) using GC$_{\rm sp}$+WL+GC$_{\rm ph}$+XC. Combining all the sources of information explored expected from $Euclid$ in combination with future SO-like CMB experiment, we forecast ${\cal A}_{\rm lin} \simeq 0.010 \pm 0.001$ at 68.3% CL and ${\cal A}_{\rm log} \simeq 0.010 \pm 0.001$ for GC$_{\rm sp}$(PS rec + BS)+WL+GC$_{\rm ph}$+XC+SO-like both for the optimistic and pessimistic settings over the frequency range $(1,\,10^{2.1})$. △ Less

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

Comments: 23 pages, 9 figures, 4 tables

Journal ref: A&A 683, A220 (2024)

Peculiar orbital characteristics of Earth quasi-satellite 469219 Kamo`oalewa: implications for the Yarkovsky detection and orbital uncertainty propagation

Authors: Shoucun Hu, Bin Li, Haoxuan Jiang, Gang Bao, Jianghui Ji

Abstract: 469219 Kamo`oalewa is selected as one of the primary targets of Tianwen-2 mission, which is currently believed to be the most stable quasi-satellite of Earth. Here we derive a weak detection of the Yarkovsky effect for Kamo`oalewa, giving $A_2 = -1.075\pm0.447\times 10^{-13} \rm{au/d}^2$, with the available ground-based optical observations from Minor Planet Center and a relatively conservative we… ▽ More 469219 Kamo`oalewa is selected as one of the primary targets of Tianwen-2 mission, which is currently believed to be the most stable quasi-satellite of Earth. Here we derive a weak detection of the Yarkovsky effect for Kamo`oalewa, giving $A_2 = -1.075\pm0.447\times 10^{-13} \rm{au/d}^2$, with the available ground-based optical observations from Minor Planet Center and a relatively conservative weighting scheme. Due to the quasi-satellite resonance with Earth, we show that the detection of Yarkovsky effect by orbital fitting with astrometric observations becomes difficult as its orbital drift shows a slow oscillatory growth resulting from the Yarkovsky effect. In addition, we extensively explore the characteristics of orbital uncertainty propagation and find that the positional uncertainty mainly arises from the geocentric radial direction in 2010-2020, and then concentrates in the heliocentric transverse direction in 2020-2030. Furthermore, the heliocentric transverse uncertainty is clearly monthly dependent, which can arrive at a minimum around January and a maximum around July as the orbit moves towards the leading and trailing edges, respectively, in 2025-2027. Finally, we investigate a long-term uncertainty propagation in the quasi-satellite regime, implying that the quasi-satellite resonance with Earth may play a crucial role in constraining the increase of uncertainty over time. Such interesting feature further implies that the orbital precision of Kamo`oalewa is relatively stable at its quasi-satellite phase, which may also be true for other quasi-satellites of Earth. △ Less

Submitted 31 August, 2023; originally announced August 2023.

Comments: 15 pages, 11 figures, accepted for publication in AJ