Showing 1–50 of 540 results for author: Yu, H

Improved background modeling for dark matter search with COSINE-100

Authors: G. H. Yu, N. Carlin, J. Y. Cho, J. J. Choi, S. Choi, A. C. Ezeribe, L. E. Franca, C. Ha, I. S. Hahn, S. J. Hollick, E. J. Jeon, H. W. Joo, W. G. Kang, M. Kauer, B. H. Kim, H. J. Kim, J. Kim, K. W. Kim, S. H. Kim, S. K. Kim, W. K. Kim, Y. D. Kim, Y. H. Kim, Y. J. Ko, D. H. Lee , et al. (33 additional authors not shown)

Abstract: COSINE-100 aims to conclusively test the claimed dark matter annual modulation signal detected by DAMA/LIBRA collaboration. DAMA/LIBRA has released updated analysis results by lowering the energy threshold to 0.75 keV through various upgrades. They have consistently claimed to have observed the annual modulation. In COSINE-100, it is crucial to lower the energy threshold for a direct comparison wi… ▽ More COSINE-100 aims to conclusively test the claimed dark matter annual modulation signal detected by DAMA/LIBRA collaboration. DAMA/LIBRA has released updated analysis results by lowering the energy threshold to 0.75 keV through various upgrades. They have consistently claimed to have observed the annual modulation. In COSINE-100, it is crucial to lower the energy threshold for a direct comparison with DAMA/LIBRA, which also enhances the sensitivity of the search for low-mass dark matter, enabling COSINE-100 to explore this area. Therefore, it is essential to have a precise and quantitative understanding of the background spectrum across all energy ranges. This study expands the background modeling from 0.7 to 4000 keV using 2.82 years of COSINE-100 data. The modeling has been improved to describe the background spectrum across all energy ranges accurately. Assessments of the background spectrum are presented, considering the nonproportionality of NaI(Tl) crystals at both low and high energies and the characteristic X-rays produced by the interaction of external backgrounds with materials such as copper. Additionally, constraints on the fit parameters obtained from the alpha spectrum modeling fit are integrated into this model. These improvements are detailed in the paper. △ Less

Submitted 19 August, 2024; originally announced August 2024.

EDEN: Exploring Disks Embedded in N-body simulations of Milky-Way-mass halos from Symphony

Authors: Yunchong Wang, Philip Mansfield, Ethan O. Nadler, Elise Darragh-Ford, Risa H. Wechsler, Daneng Yang, Hai-Bo Yu

Abstract: We investigate the impact of galactic disks on the tidal stripping of cold dark matter subhalos within Milky Way (MW)-mass halos ($M_{\rm vir}\sim 10^{12}\mathrm{M_{\odot}}$) using a new simulation suite, EDEN. By re-simulating 45 MW-mass zoom-in halos from the N-body Symphony compilation with embedded disk potentials, which evolve according to star formation histories predicted by the UniverseMac… ▽ More We investigate the impact of galactic disks on the tidal stripping of cold dark matter subhalos within Milky Way (MW)-mass halos ($M_{\rm vir}\sim 10^{12}\mathrm{M_{\odot}}$) using a new simulation suite, EDEN. By re-simulating 45 MW-mass zoom-in halos from the N-body Symphony compilation with embedded disk potentials, which evolve according to star formation histories predicted by the UniverseMachine model, we self-consistently tie disk growth to halo accretion rate and significantly expand the range of disk masses and formation histories studied. We use the particle-tracking-based subhalo finder Symfind to enhance the robustness of subhalo tracking. We find that disks near the median disk-to-halo mass ratio of our sample ($M_{\ast, \rm Disk}/M_{\rm vir, host} = 2\%$) reduce subhalo peak mass functions within 100 kpc by about $10\%$ for peak masses above $ 10^8\mathrm{M_{\odot}}$. Heavier, MW/M31-like disks ($M_{\ast, \rm Disk}/M_{\rm vir, host} \gtrsim 5\%$) lead to a reduction of more than $40\%$. Subhalo abundance suppression is more pronounced near halo centers, with particularly enhanced stripping for subhalos accreted over 8 Gyr ago on orbits with pericenters < 100 kpc. Suppression is further amplified when disk mass is increased within fixed halo and disk assembly histories. In all cases, the suppression we measure should be interpreted as stripping below the mass resolution limit rather than complete subhalo disruption. This study reshapes our understanding of the MW's impact on its satellites, suggesting it strips subhalos more efficiently than typical MW-mass galaxies due to its larger disk-to-halo mass ratio and earlier disk formation. △ Less

Submitted 2 August, 2024; originally announced August 2024.

Comments: 30 pages, 13 figures, 1 table. Submitted to ApJ, comments welcome

Study of Wide-Field-of-View X-ray Observations of the Virgo Cluster Using the Lobster Eye Imager for Astronomy

Authors: Wen-Cheng Feng, Shu-Mei Jia, Hai-Hui Zhao, Heng Yu, Hai-Wu Pan, Cheng-Kui Li, Yu-Lin Cheng, Shan-Shan Weng, Yong Chen, Yuan Liu, Zhi-Xing Ling, Chen Zhang

Abstract: The Lobster Eye Imager for Astronomy (LEIA) is the pathfinder of the wide-field X-ray telescope used in the Einstein Probe mission. In this study, we present an image of the Virgo Cluster taken by LEIA in the 0.5-4.5 keV band with an exposure time of $\sim$17.3 ks in the central region. This extended emission is generally consistent with the results obtained by ROSAT. However, the field is affecte… ▽ More The Lobster Eye Imager for Astronomy (LEIA) is the pathfinder of the wide-field X-ray telescope used in the Einstein Probe mission. In this study, we present an image of the Virgo Cluster taken by LEIA in the 0.5-4.5 keV band with an exposure time of $\sim$17.3 ks in the central region. This extended emission is generally consistent with the results obtained by ROSAT. However, the field is affected by bright point sources due to the instrument's Point Spread Function (PSF) effect. Through fitting of the LEIA spectrum of the Virgo Cluster, we obtained a temperature of $2.1^{+0.3}_{-0.1}$ keV, which is consistent with the XMM-Newton results ($\sim$2.3 keV). Above 1.6 keV, the spectrum is dominated by the X-ray background. In summary, this study validates LEIA's extended source imaging and spectral resolution capabilities for the first time. △ Less

Submitted 31 July, 2024; originally announced August 2024.

Comments: 9 pages, 6 figures, 1 table

Revisiting the fundamental metallicity relation with observation and simulation

Authors: Chengyu Ma, Kai Wang, Enci Wang, Yingjie Peng, Haochen Jiang, Haoran Yu, Cheng Jia, Zeyu Chen, Haixin Li, Xu Kong

Abstract: The gas-phase metallicity of galaxies is regulated by multiple astrophysical processes, which makes it a crucial diagnostic of galaxy formation and evolution. Beyond the fundamental mass-metallicity relation, a debate about the secondary galaxy property to predict the metallicity of galaxies arises. Motivated by this, we systematically examine the relationship between gas-phase metallicity and oth… ▽ More The gas-phase metallicity of galaxies is regulated by multiple astrophysical processes, which makes it a crucial diagnostic of galaxy formation and evolution. Beyond the fundamental mass-metallicity relation, a debate about the secondary galaxy property to predict the metallicity of galaxies arises. Motivated by this, we systematically examine the relationship between gas-phase metallicity and other galaxy properties, i.e. star formation rate (SFR) and galaxy size, in addition to stellar mass in both observation and simulation. We utilize the data from the MaNGA (Mapping Nearby Galaxies at Apache Point Observatory) survey and the TNG50 simulations. We find that the combination of $M_*/R_{\rm e}^β$ with $β\sim 0.6-1$ is in much stronger correlation to the metallicity than stellar mass alone, regardless of whether the SFR is included or not, in both observation and simulation. This indicates that galaxy size plays a more important role in determining gas-phase metallicity of galaxies than SFR. In addition, the TNG simulation predicts that the SFR, although being a subdominant role, becomes increasingly important in high-$z$ universe. Finally, we speculate that SFR modulates metallicity on the temporal dimension, synchronized with time-varying gas inflows, and galaxy size regulates metallicity on the spatial dimension by affecting the gravitational potential and the mass loading factor. △ Less

Submitted 31 July, 2024; originally announced July 2024.

Comments: 11 pages, 6 figures, ApJL accepted

FAST observations of neutral hydrogen in the interacting galaxies NGC 3395/3396

Authors: Nai-Ping Yu, Ming Zhu, Jin-Long Xu, Chuan-Peng Zhang, Hai-Yang Yu, Xiao-Lan Liu, Peng Jiang, Mei Ai

Abstract: We report on high-sensitivity neutral hydrogen observations toward the gas-rich interacting galaxies NGC 3395/3396 with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Compared to previous observations carried out by the Very Large Array (VLA) and the Westerbork Synthesis Radio Telescope (WSRT), a more extended HI envelope around this system has been detected. The total HI gas ma… ▽ More We report on high-sensitivity neutral hydrogen observations toward the gas-rich interacting galaxies NGC 3395/3396 with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Compared to previous observations carried out by the Very Large Array (VLA) and the Westerbork Synthesis Radio Telescope (WSRT), a more extended HI envelope around this system has been detected. The total HI gas mass of the NGC 3395/3396 system is estimated to be 7.8 x 109 M. This value is 2.7 times more than that reported based on the VLA interferometric maps. Previous observations found a large HI tail extending to the south-west and a minor tail emerging from the north of this peculiar galaxy pair. Based on the high-sensitivity observations of FAST, an extended HI plume to the north-west and a gas plume to the north-east have been detected for the first time. Neutral hydrogen of the two smaller galaxies IC 2604 and IC 2608 on the south of the system have also been detected. We discuss the origins of these extra gas and possible tidal interactions between these galaxies. NGC 3395/3396's most prominent tidal feature, the south-west tail combined with the new detected north-west plume behaves like a large ring. We suggest the ring might be formed by the previous fly-by interaction between NGC 3395 and NGC 3396 which happened 500 Myr ago. Our study shows that high-sensitivity HI observations are important in revealing low column density gas, which is crucial to a deeper understanding of this interacting system. △ Less

Submitted 30 July, 2024; originally announced July 2024.

Revisiting the dynamical masses of the transiting planets in the young AU Mic system: Potential AU Mic b inflation at $\sim$20 Myr

Authors: M. Mallorquín, V. J. S. Béjar, N. Lodieu, M. R. Zapatero Osorio, H. Yu, A. Suárez Mascareño, M. Damasso, J. Sanz-Forcada, I. Ribas, A. Reiners, A. Quirrenbach, P. J. Amado, J. A. Caballero, S. Aigrain, O. Barragán, S. Dreizler, A. Fernández-Martín, E. Goffo, Th. Henning, A. Kaminski, B. Klein, R. Luque, D. Montes, J. C. Morales, E. Nagel , et al. (4 additional authors not shown)

Abstract: Understanding planet formation is important in the context of the origin of planetary systems in general and of the Solar System in particular, as well as to predict the likelihood of finding Jupiter, Neptune, and Earth analogues around other stars. We aim to precisely determine the radii and dynamical masses of transiting planets orbiting the young M star AU Mic using public photometric and spect… ▽ More Understanding planet formation is important in the context of the origin of planetary systems in general and of the Solar System in particular, as well as to predict the likelihood of finding Jupiter, Neptune, and Earth analogues around other stars. We aim to precisely determine the radii and dynamical masses of transiting planets orbiting the young M star AU Mic using public photometric and spectroscopic datasets. We characterise the stellar activity and physical properties (radius, mass, density) of the transiting planets in the young AU Mic system through joint transit and radial velocity fits with Gaussian processes. We determine a radius of $R^{b}$= 4.79 +/- 0.29 R$_\oplus$, a mass of $M^{b}$= 9.0 +/- 2.7 M$_\oplus$, and a bulk density of $ρ^{b}$ = 0.49 +/- 0.16 g cm$^{-3}$ for the innermost transiting planet AU Mic b. For the second known transiting planet, AU Mic c, we infer a radius of $R^{c}$= 2.79 +/- 0.18 R$_\oplus$, a mass of $M^{c}$= 14.5 +/- 3.4 M$_\oplus$, and a bulk density of $ρ^{c}$ = 3.90 +/- 1.17 g cm$^{-3}$. According to theoretical models, AU Mic b may harbour an H2 envelope larger than 5\% by mass, with a fraction of rock and a fraction of water. AU Mic c could be made of rock and/or water and may have an H2 atmosphere comprising at most 5\% of its mass. AU Mic b has retained most of its atmosphere but might lose it over tens of millions of years due to the strong stellar radiation, while AU Mic c likely suffers much less photo-evaporation because it lies at a larger separation from its host. Using all the datasets in hand, we determine a 3$σ$ upper mass limit of $M^{[d]}\sin{i}$ = 8.6 M$_{\oplus}$ for the AU Mic 'd' TTV-candidate. In addition, we do not confirm the recently proposed existence of the planet candidate AU Mic 'e' with an orbital period of 33.4 days. △ Less

Submitted 23 July, 2024; originally announced July 2024.

Comments: Accepted in A&A. 28 pages, 15 figures

Gravothermal collapse and the diversity of galactic rotation curves

Authors: M. Grant Roberts, Manoj Kaplinghat, Mauro Valli, Hai-Bo Yu

Abstract: The rotation curves of spiral galaxies exhibit a great diversity that challenge our understanding of galaxy formation and the nature of dark matter. Previous studies showed that in self-interacting dark matter (SIDM) models with a cross section per unit mass of $σ/m\approx{\cal O}(1)~{\rm cm^2/g}$, the predicted dark matter central densities are a good match to the observed densities in galaxies.… ▽ More The rotation curves of spiral galaxies exhibit a great diversity that challenge our understanding of galaxy formation and the nature of dark matter. Previous studies showed that in self-interacting dark matter (SIDM) models with a cross section per unit mass of $σ/m\approx{\cal O}(1)~{\rm cm^2/g}$, the predicted dark matter central densities are a good match to the observed densities in galaxies. In this work, we explore a regime with a larger cross section of $σ/m\approx20-40~{\rm cm^2/g}$ in dwarf galactic halos. We will show that such strong dark matter self-interactions can further amplify the diversity of halo densities inherited from their assembly history. High concentration halos can enter the gravothermal collapse phase within $10~{\rm Gyr}$, resulting in a high density, while low concentration ones remain in the expansion phase and have a low density. We fit the rotation curves of $14$ representative low surface brightness galaxies and demonstrate how the large range of observed central densities are naturally accommodated in the strong SIDM regime of $σ/m\approx20-40~{\rm cm^2/g}$. Galaxies that are outliers in the previous studies due to their high halo central densities, are no longer outliers in this SIDM regime as their halos would be in the collapse phase. For galaxies with a low density, the SIDM fits are robust to the variation of the cross section. Our findings open up a new window for testing gravothermal collapse, the unique signature of strong dark matter self-interactions, and exploring broad SIDM model space. △ Less

Submitted 20 July, 2024; originally announced July 2024.

Comments: 21 pages, 14 figures

Swift-BAT GUANO follow-up of gravitational-wave triggers in the third LIGO-Virgo-KAGRA observing run

Authors: Gayathri Raman, Samuele Ronchini, James Delaunay, Aaron Tohuvavohu, Jamie A. Kennea, Tyler Parsotan, Elena Ambrosi, Maria Grazia Bernardini, Sergio Campana, Giancarlo Cusumano, Antonino D'Ai, Paolo D'Avanzo, Valerio D'Elia, Massimiliano De Pasquale, Simone Dichiara, Phil Evans, Dieter Hartmann, Paul Kuin, Andrea Melandri, Paul O'Brien, Julian P. Osborne, Kim Page, David M. Palmer, Boris Sbarufatti, Gianpiero Tagliaferri , et al. (1797 additional authors not shown)

Abstract: We present results from a search for X-ray/gamma-ray counterparts of gravitational-wave (GW) candidates from the third observing run (O3) of the LIGO-Virgo-KAGRA (LVK) network using the Swift Burst Alert Telescope (Swift-BAT). The search includes 636 GW candidates received in low latency, 86 of which have been confirmed by the offline analysis and included in the third cumulative Gravitational-Wav… ▽ More We present results from a search for X-ray/gamma-ray counterparts of gravitational-wave (GW) candidates from the third observing run (O3) of the LIGO-Virgo-KAGRA (LVK) network using the Swift Burst Alert Telescope (Swift-BAT). The search includes 636 GW candidates received in low latency, 86 of which have been confirmed by the offline analysis and included in the third cumulative Gravitational-Wave Transient Catalogs (GWTC-3). Targeted searches were carried out on the entire GW sample using the maximum--likelihood NITRATES pipeline on the BAT data made available via the GUANO infrastructure. We do not detect any significant electromagnetic emission that is temporally and spatially coincident with any of the GW candidates. We report flux upper limits in the 15-350 keV band as a function of sky position for all the catalog candidates. For GW candidates where the Swift-BAT false alarm rate is less than 10$^{-3}$ Hz, we compute the GW--BAT joint false alarm rate. Finally, the derived Swift-BAT upper limits are used to infer constraints on the putative electromagnetic emission associated with binary black hole mergers. △ Less

Submitted 13 July, 2024; originally announced July 2024.

Comments: 50 pages, 10 figures, 4 tables

Supernova Pointing Capabilities of DUNE

Authors: DUNE Collaboration, A. Abed Abud, B. Abi, R. Acciarri, M. A. Acero, M. R. Adames, G. Adamov, M. Adamowski, D. Adams, M. Adinolfi, C. Adriano, A. Aduszkiewicz, J. Aguilar, B. Aimard, F. Akbar, K. Allison, S. Alonso Monsalve, M. Alrashed, A. Alton, R. Alvarez, T. Alves, H. Amar, P. Amedo, J. Anderson, D. A. Andrade , et al. (1340 additional authors not shown)

Abstract: The determination of the direction of a stellar core collapse via its neutrino emission is crucial for the identification of the progenitor for a multimessenger follow-up. A highly effective method of reconstructing supernova directions within the Deep Underground Neutrino Experiment (DUNE) is introduced. The supernova neutrino pointing resolution is studied by simulating and reconstructing electr… ▽ More The determination of the direction of a stellar core collapse via its neutrino emission is crucial for the identification of the progenitor for a multimessenger follow-up. A highly effective method of reconstructing supernova directions within the Deep Underground Neutrino Experiment (DUNE) is introduced. The supernova neutrino pointing resolution is studied by simulating and reconstructing electron-neutrino charged-current absorption on $^{40}$Ar and elastic scattering of neutrinos on electrons. Procedures to reconstruct individual interactions, including a newly developed technique called ``brems flipping'', as well as the burst direction from an ensemble of interactions are described. Performance of the burst direction reconstruction is evaluated for supernovae happening at a distance of 10 kpc for a specific supernova burst flux model. The pointing resolution is found to be 3.4 degrees at 68% coverage for a perfect interaction-channel classification and a fiducial mass of 40 kton, and 6.6 degrees for a 10 kton fiducial mass respectively. Assuming a 4% rate of charged-current interactions being misidentified as elastic scattering, DUNE's burst pointing resolution is found to be 4.3 degrees (8.7 degrees) at 68% coverage. △ Less

Submitted 14 July, 2024; originally announced July 2024.

Comments: 25 pages, 16 figures

Report number: FERMILAB-PUB-24-0319-LBNF

Hyper-sampling imaging

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

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

Submitted 27 June, 2024; originally announced June 2024.

Testing the parametric model for self-interacting dark matter using matched halos in cosmological simulations

Authors: Daneng Yang, Ethan O. Nadler, Hai-Bo Yu

Abstract: We systemically evaluate the performance of the self-interacting dark matter (SIDM) halo model proposed in arXiv:2305.16176 with matched halos from high-resolution cosmological CDM and SIDM simulations. The model incorporates SIDM effects along mass evolution histories of CDM halos and it is applicable to both isolated halos and suhbhalos. We focus on the accuracy of the model in predicting halo d… ▽ More We systemically evaluate the performance of the self-interacting dark matter (SIDM) halo model proposed in arXiv:2305.16176 with matched halos from high-resolution cosmological CDM and SIDM simulations. The model incorporates SIDM effects along mass evolution histories of CDM halos and it is applicable to both isolated halos and suhbhalos. We focus on the accuracy of the model in predicting halo density profiles at $z=0$ and the evolution of maximum circular velocity. We find the model predictions agree with the simulations within $10\%-50\%$ for most of the simulated (sub)halos, $50\%-100\%$ for extreme cases. This indicates that the model effectively captures the gravothermal evolution of the halos with very strong, velocity-dependent self-interactions. For an example application, we apply the model to study the impact of various SIDM scenarios on strong lensing perturber systems, demonstrating its utility in predicting SIDM effects for small-scale structure analyses. Our findings confirm that the model is an effective tool for mapping CDM halos into their SIDM counterparts. △ Less

Submitted 15 June, 2024; originally announced June 2024.

Comments: 20 pages, 19 figures

Interstellar Nitrogen Isotope Ratios: Measurements on tracers of C$^{14}$N and C$^{15}$N

Authors: J. L. Chen, J. S. Zhang, C. Henkel, Y. T. Yan, H. Z. Yu, Y. X. Wang, Y. P. Zou, J. Y. Zhao, X. Y. Wang

Abstract: The nitrogen isotope ratio 14N/15N is a powerful tool to trace Galactic stellar nucleosynthesis and constraining Galactic chemical evolution. Previous observations have found lower 14N/15N ratios in the Galactic center and higher values in the Galactic disk. This is consistent with the inside-out formation scenario of our Milky Way. However, previous studies mostly utilized double isotope ratios a… ▽ More The nitrogen isotope ratio 14N/15N is a powerful tool to trace Galactic stellar nucleosynthesis and constraining Galactic chemical evolution. Previous observations have found lower 14N/15N ratios in the Galactic center and higher values in the Galactic disk. This is consistent with the inside-out formation scenario of our Milky Way. However, previous studies mostly utilized double isotope ratios also including 12C/13C, which introduces additional uncertainties. Here we therefore present observations of C14N and its rare isotopologue, C15N, toward a sample of star forming regions, measured by the IRAM 30 m and/or the ARO 12 m telescope at $λ$ ~3 mm wavelength. For those 35 sources detected in both isotopologues, physical parameters are determined. Furthermore we have obtained nitrogen isotope ratios using the strongest hyperfine components of CN and C15N. For those sources showing small deviations from Local Thermodynamical Equilibrium and/or self-absorption, the weakest hyperfine component, likely free of the latter effect, was used to obtain reliable 14N/15N values. Our measured 14N/15N isotope ratios from C14N and C15N measurements are compatible with those from our earlier measurements of NH3 and 15NH3 (Paper I), i.e., increasing ratios to a Galacticentric distance of ~9 kpc. The unweighted second order polynomial fit yields $\frac{{\rm C^{14}N}}{{\rm C^{15}N}} = (-4.85 \pm 1.89)\;{\rm kpc^{-2}} \times R_{\rm GC}^{2} + (82.11 \pm 31.93) \;{\rm kpc^{-1}} \times R_{\rm GC} - (28.12 \pm 126.62)$. Toward the outer galaxy, the isotope ratio tends to decrease, supporting an earlier finding by H13CN/HC15N. Galactic chemical evolution models are consistent with our measurements of the 14N/15N isotope ratio, i.e. a rising trend from the Galactic center region to approximately 9 kpc, followed by a decreasing trend with increasing $R_{\rm GC}$ toward the outer Galaxy. △ Less

Submitted 13 June, 2024; originally announced June 2024.

Comments: 34 pages, 9 figures, 6 tables

Journal ref: The Astrophysical Journal (2004)

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

The destiny of open cluster NGC 6530: past and future

Authors: Delong Jia, Heng Yu, Zhengyi Shao, Lu Li

Abstract: Studying the structures of open clusters is crucial for understanding stellar evolution and galactic dynamics. Based on Gaia DR3 data, we apply the hierarchical clustering algorithm to a young open cluster NGC 6530 and group its members into 5 substructures. By linear tracing with the kinematic information of their members, we find that: Sub 1 is the core of the cluster. It is expanding slowly. Su… ▽ More Studying the structures of open clusters is crucial for understanding stellar evolution and galactic dynamics. Based on Gaia DR3 data, we apply the hierarchical clustering algorithm to a young open cluster NGC 6530 and group its members into 5 substructures. By linear tracing with the kinematic information of their members, we find that: Sub 1 is the core of the cluster. It is expanding slowly. Sub 2 consists of less bound members, which began escaping from the core about 0.78 Myr ago. Sub 3 is associated with a young star forming region. It will merge with the core after 0.72 Myr; Sub 4, as an outskirt group, is also moving towards the core, but won't end up falling in. While Sub 5 is composed of less-bound members with field contamination. This work reveals the complex internal structure and evolutionary trends of the cluster NGC 6530. It also shows the potential of the hierarchical clustering algorithm in star cluster structure analysis. △ Less

Submitted 14 July, 2024; v1 submitted 11 June, 2024; originally announced June 2024.

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

Probing vector chirality in the early Universe

Authors: Junsup Shim, Ue-Li Pen, Hao-Ran Yu, Teppei Okumura

Abstract: We explore the potential of detecting parity violation in primordial vector fossils using late-time galaxy spins. Utilizing $N$-body simulations, we use halo spins as a reliable proxy for galaxy spins to investigate how effectively such primordial vectorial parity asymmetry remains in galaxy spins at low redshifts. We develop a novel approach to generate initial conditions with substantial parity… ▽ More We explore the potential of detecting parity violation in primordial vector fossils using late-time galaxy spins. Utilizing $N$-body simulations, we use halo spins as a reliable proxy for galaxy spins to investigate how effectively such primordial vectorial parity asymmetry remains in galaxy spins at low redshifts. We develop a novel approach to generate initial conditions with substantial parity asymmetry, while maintaining the initial matter power spectrum unchanged. From the parity broken initial condition and halos evolved from it, we construct the initial spin and halo spin fields, respectively. Focusing on the helicity of these vector fields, we detect substantial asymmetry in the initial spin field as a consequence of parity violation in the primordial vector fossil. In addition, we discover that over $50\%$ of the primordial asymmetry in the initial spin field remains in the late-time halo spin field on a range of scales. Given the tight correlation between halo spins and observable galaxy spins, we expect to detect the current amplitude of vectorial parity asymmetry potentially up to $16σ$-level in observation, when utilizing galaxy samples from DESI BGS. Our findings demonstrate that the primordial imprints of vectorial parity violation persist through non-linear gravitational evolution, highlighting the reliability of galaxy spin as a sensitive probe for testing the vectorial parity-invariance in the early Universe. △ Less

Submitted 10 June, 2024; originally announced June 2024.

Comments: 5 pages, 3 figures. Submitted to PRL

Alleviating the Hubble-constant tension and the growth tension via a transition of absolute magnitude favored by the Pantheon+ sample

Authors: Yang Liu, Hongwei Yu, Puxun Wu

Abstract: We establish a cosmological-model-independent method to extract the apparent magnitude and its derivative at different redshifts from the Pantheon+ type Ia supernova sample, and find that the obtained values deviate clearly from the prediction of the $Λ$CDM model at the lowest redshift. This deviation can be explained as a result of a transition of the absolute magnitude $M$ in the low redshift re… ▽ More We establish a cosmological-model-independent method to extract the apparent magnitude and its derivative at different redshifts from the Pantheon+ type Ia supernova sample, and find that the obtained values deviate clearly from the prediction of the $Λ$CDM model at the lowest redshift. This deviation can be explained as a result of a transition of the absolute magnitude $M$ in the low redshift region. The observations seem to favor this transition since the minimum values of $χ^2$ for two ansatzes of a varying $M$ are less than that of a constant $M$. The Hubble constant tension is alleviated from larger than $5σ$ to be about $1$ to $2σ$ for a varying $M$, and the growth tension can be resolved after attributing the variation of $M$ to a modification of the effective Newton's constant. △ Less

Submitted 24 July, 2024; v1 submitted 5 June, 2024; originally announced June 2024.

Comments: 20 pages, 5 figures, 4 tables. Published in Physical Review D (Letter)

Are WASP-107-like Systems Consistent with High-eccentricity Migration?

Authors: Hang Yu, Fei Dai

Abstract: WASP-107 b seems to be a poster child of the long-suspected high-eccentricity migration scenario. It is on a 5.7-day, polar orbit. The planet is Jupiter-like in radius but Neptune-like in mass with exceptionally low density. WASP-107 c is on a 1100-day, $e=0.28$ orbit with at least Saturn mass. Planet b may still have a residual eccentricity of $0.06\pm 0.04$: the ongoing tidal dissipation leads t… ▽ More WASP-107 b seems to be a poster child of the long-suspected high-eccentricity migration scenario. It is on a 5.7-day, polar orbit. The planet is Jupiter-like in radius but Neptune-like in mass with exceptionally low density. WASP-107 c is on a 1100-day, $e=0.28$ orbit with at least Saturn mass. Planet b may still have a residual eccentricity of $0.06\pm 0.04$: the ongoing tidal dissipation leads to the observed internally heated atmosphere and hydrodynamic atmospheric erosion. We present a population synthesis study coupling octopole Lidov-Kozai oscillations with various short-range forces, while simultaneously accounting for the radius inflation and tidal disruption of the planet. We find that a high-eccentricity migration scenario can successfully explain nearly all observed system properties. Our simulations further suggest that the initial location of WASP-107 b at the onset of migration is likely within the snowline ($<0.5\,{\rm AU}$). More distant initial orbits usually lead to tidal disruption or orbit crossing. WASP-107 b most likely lost no more than 20% of its mass during the high-eccentricity migration, i.e. it did not form as a Jupiter-mass object. More vigorous tidally-induced mass loss leads to disruption of the planet during migration. We predict that the current-day mutual inclination between the planets b and c is substantial: at least 25-55$^\circ$ which may be tested with future Gaia astrometric observations. Knowing the current-day mutual inclination may further constrain the initial orbit of planet b. We suggest that the proposed high-eccentricity migration scenario of WASP-107 may be applicable to HAT-P-11, GJ-3470, HAT-P-18, and GJ-436 which have similar orbital architectures. △ Less

Submitted 9 July, 2024; v1 submitted 31 May, 2024; originally announced June 2024.

Comments: 19 pages, 10 figures. Accepted by ApJ

Nonthermal Acceleration of Electrons, Positrons and Protons at a Nonrelativistic Quasiparallel Collisionless Shock

Authors: Huan Yu, Qi Xia, Jun Fang

Abstract: Energetic positrons have been observed in the interstellar medium, and high-energy positrons with relativistic energies up to approximately 1 TeV have been detected in Galactic cosmic rays. We conducted a study on the acceleration of particles, specifically positrons, in a nonrelativistic quasiparallel collisionless shock induced by a plasma consisting of protons, electrons, and positrons. The pos… ▽ More Energetic positrons have been observed in the interstellar medium, and high-energy positrons with relativistic energies up to approximately 1 TeV have been detected in Galactic cosmic rays. We conducted a study on the acceleration of particles, specifically positrons, in a nonrelativistic quasiparallel collisionless shock induced by a plasma consisting of protons, electrons, and positrons. The positron-to-proton number density ratio in the plasma is 0.1. We focused on a representative shock with a sonic Mach number of 17.1 and an Alfvénic Mach number of 16.8 in the rest frame of the shock. To investigate the acceleration mechanisms of particles including positrons in the shock, we utilized one-dimensional particle-in-cell (PIC) simulations. It was found that all three species of particles in the shock can be accelerated and exhibit power law spectra. At the shock front, a significant portion of incoming upstream particles are reflected and undergo significant energy increase, and these reflected particles can be efficiently injected into the process of diffusive shock acceleration (DSA). Moveover, the reflected positrons can be further accelerated by an electric field parallel to the magnetic field when they move along the magnetic field upstream of the shock. As a result, positrons can be preferentially accelerated to be injected in the DSA process compared to electrons. △ Less

Submitted 29 May, 2024; originally announced May 2024.

Comments: 7 pages, 7 figures, Accepted for publication in ApJ

Investigating stellar activity through eight years of Sun-as-a-star observations

Authors: Baptiste Klein, Suzanne Aigrain, Michael Cretignier, Khaled Al Moulla, Xavier Dumusque, Oscar Barragán, Haochuan Yu, Annelies Mortier, Federica Rescigno, Andrew Collier Cameron, Mercedes López-Morales, Nadège Meunier, Alessandro Sozzetti, Niamh K. O'Sullivan

Abstract: Stellar magnetic activity induces both distortions and Doppler-shifts in the absorption line profiles of Sun-like stars. Those effects produce apparent radial velocity (RV) signals which greatly hamper the search for potentially habitable, Earth-like planets. In this work, we investigate these distortions in the Sun using cross-correlation functions (CCFs), derived from intensive monitoring with t… ▽ More Stellar magnetic activity induces both distortions and Doppler-shifts in the absorption line profiles of Sun-like stars. Those effects produce apparent radial velocity (RV) signals which greatly hamper the search for potentially habitable, Earth-like planets. In this work, we investigate these distortions in the Sun using cross-correlation functions (CCFs), derived from intensive monitoring with the high-precision spectrograph HARPS-N. We show that the RV signal arising from line-shape variations on time-scales associated with the solar rotation and activity cycle can be robustly extracted from the data, reducing the RV dispersion by half. Once these have been corrected, activity-induced Doppler-shifts remain, that are modulated at the solar rotation period, and that are most effectively modelled in the time domain, using Gaussian Processes (GPs). Planet signatures are still best retrieved with multi-dimensonal GPs, when activity is jointly modelled from the raw RVs and indicators of the line width or of the Ca II H and K emission. After GP modelling, the residual RVs exhibit a dispersion of 0.6-0.8 m/s, likely to be dominated by signals induced by super-granulation. Finally, we find that the statistical properties of the RVs evolve significantly over time, and that this evolution is primarily driven by sunspots, which control the smoothness of the signal. Such evolution, which reduces the sensitivity to long-period planet signatures, is no longer seen in the activity-induced Doppler-shifts, which is promising for long term RV monitoring surveys such as the Terra Hunting Experiment or the PLATO follow-up campaign. △ Less

Submitted 22 May, 2024; v1 submitted 20 May, 2024; originally announced May 2024.

Comments: Accepted by MNRAS on the 2024 May 20. 17 pages, 14 figures (plus 8 pages of Appendix, 9 Appendix figures)

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

Uncovering Stealth Bias in LISA observations of Double White Dwarf Binaries due to Tidal Coupling

Authors: Grace Fiacco, Neil J. Cornish, Hang Yu

Abstract: Double white dwarfs are important gravitational wave sources for LISA, as they are some of the most numerous compact systems in our universe. Here we consider finite-sized effects due to tidal interactions, as they are expected to have a measurable impact on these systems. Previous studies suggested that tidal effects would allow the individual masses to be measured, but there was a subtle error i… ▽ More Double white dwarfs are important gravitational wave sources for LISA, as they are some of the most numerous compact systems in our universe. Here we consider finite-sized effects due to tidal interactions, as they are expected to have a measurable impact on these systems. Previous studies suggested that tidal effects would allow the individual masses to be measured, but there was a subtle error in those analyses. Using a fully Bayesian analysis we find that while tidal effects do not allow us to constrain the individual masses, they do yield informative lower bounds on the total mass of the system. Including tidal effects is crucial to the accuracy of our estimation of the chirp and total mass. Neglecting tidal effects leads to significant biases towards higher chirp masses, and we see that the lower bound of the total masses is biased towards a higher value as well. For many systems observed by LISA, tidal effects can lead to a "stealth" bias, since only the first derivative of the frequency can be measured. To separate tidal effects from the usual point particle decay we need to be able to measure the change in the second derivative of the frequency cause by the tides. This can only be done for high frequency systems observed with high signal-to-noise. The bias, if not accounted for, can have significant astrophysical implications; for example, it could lead to an incorrect estimation of the population of potential Type IA supernovae progenitors. △ Less

Submitted 16 May, 2024; originally announced May 2024.

Comments: 19 pages, 17 figures, to be submitted to PRD

Enhancing Morphological Measurements of Cosmic Web with Delaunay Tessellation Field Estimation

Authors: Yu Liu, Yu Yu, Pengjie Zhang, Hao-Ran Yu

Abstract: The density fields constructed by traditional mass assignment methods are susceptible to irritating discreteness, which hinders morphological measurements of cosmic large-scale structure (LSS) through Minkowski functionals (MFs). For alleviating this issue, fixed-kernel smoothing methods are commonly used in literatures, at the expense of losing substantial structural information. In this work, we… ▽ More The density fields constructed by traditional mass assignment methods are susceptible to irritating discreteness, which hinders morphological measurements of cosmic large-scale structure (LSS) through Minkowski functionals (MFs). For alleviating this issue, fixed-kernel smoothing methods are commonly used in literatures, at the expense of losing substantial structural information. In this work, we propose to measure MFs with Delaunay tessellation field estimation (DTFE) technique, with the goal to maximize extractions of morphological information from sparse tracers. We perform our analyses starting from matter fields and progressively extending to halo fields. At matter field level, we elucidate how discreteness affects the morphological measurements of LSS. Then, by comparing with traditional Gaussian smoothing scheme, we preliminarily showcase the advantages of DTFE for enhancing measurements of MFs from sparse tracers. At halo field level, we first numerically investigate various systematic effects on MFs of DTFE fields, which are induced by finite voxel sizes, halo number densities, halo weightings, and redshift space distortions (RSDs), respectively. Then, we explore the statistical power of MFs measured with DTFE for extracting cosmological information encoded in RSDs. We find that MFs measured with DTFE exhibit improvements by $\sim$ $2$ orders of magnitude in discriminative power for RSD effects and by a factor of $\sim$ $3$-$5$ in constraining power on structure growth rate over the MFs measured with Gaussian smoothing. These findings demonstrate the remarkable enhancements in statistical power of MFs achieved by DTFE, showing enormous application potentials of our method in extracting various key cosmological information from galaxy surveys. △ Less

Submitted 19 July, 2024; v1 submitted 14 May, 2024; originally announced May 2024.

Comments: 19 pages, 15 figures, and 1 table; accepted for publication in ApJS

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

Lagrangian space remapping and the angular momentum reconstruction from cosmic structures

Authors: Sijia Li, Ming-Jie Sheng, Haikun Li, Hao-Ran Yu

Abstract: Large scale structures provide valuable information of the primordial perturbations that encode the secrets of the origin of the Universe. It is an essential step to map between observables and their initial coordinates, called Lagrangian space, from which primordial perturbations transfer their information to structures via linear theory. By using numerical simulations and state-of-the-art recons… ▽ More Large scale structures provide valuable information of the primordial perturbations that encode the secrets of the origin of the Universe. It is an essential step to map between observables and their initial coordinates, called Lagrangian space, from which primordial perturbations transfer their information to structures via linear theory. By using numerical simulations and state-of-the-art reconstruction techniques, we report the accuracy of estimating the Lagrangian coordinates of galaxies and galaxy clusters, represented by dark matter halos in various ranges of mass, and study the accuracy of this remapping on the angular momentum (spin) reconstruction. Our work shows that galaxy groups and clusters, represented by halos with mass $\gtrsim 10^{13}M_\odot$, can be accurately remapped to Lagrangian space, and their spin reconstruction errors are dominated by the reconstructed initial gravitational potential. For all mass ranges, the errors of Lagrangian remapping, as well as redshift space distortions, play subdominant roles in estimating their angular momenta. This study explains the low correlation level between observed galaxy spins and reconstructed cosmic initial conditions and illustrates the potential of using angular momenta of cosmic structures to improve the reconstruction of primordial perturbations. △ Less

Submitted 8 July, 2024; v1 submitted 7 May, 2024; originally announced May 2024.

Comments: 7 pages, 4 figures, 1 table. Matches the accepted version in Physical Review D

Journal ref: Physical Review D, 110, 023512, 2024

Illuminating Black Hole Shadow with Dark Matter Annihilation

Authors: Yifan Chen, Ran Ding, Yuxin Liu, Yosuke Mizuno, Jing Shu, Haiyue Yu, Yanjie Zeng

Abstract: The Event Horizon Telescope (EHT) has revolutionized our ability to study black holes by providing unprecedented spatial resolution and unveiling horizon-scale details. With advancements leading to the next-generation EHT, there is potential to probe even deeper into the black hole's dark region, especially the inner shadow characterized by low-intensity foreground emissions from the jet, thanks t… ▽ More The Event Horizon Telescope (EHT) has revolutionized our ability to study black holes by providing unprecedented spatial resolution and unveiling horizon-scale details. With advancements leading to the next-generation EHT, there is potential to probe even deeper into the black hole's dark region, especially the inner shadow characterized by low-intensity foreground emissions from the jet, thanks to a significant enhancement in dynamic range by two orders of magnitude. We demonstrate how such enhanced observations could transform supermassive black holes into powerful probes for detecting annihilating dark matter, which can form a dense profile in the vicinity of supermassive black holes, by examining the morphology of the black hole image. △ Less

Submitted 2 May, 2024; v1 submitted 25 April, 2024; originally announced April 2024.

Comments: 19 pages, 6 figures, references updated

Constraints on Cosmological Models from Quasars Calibrated with Type Ia Supernova by a Gaussian Process

Authors: Haixiang Zhang, Yang Liu, Hongwei Yu, Xiaodong Nong, Nan Liang, Puxun Wu

Abstract: In this paper, we use quasars calibrated from type Ia supernova (SN Ia) to constrain cosmological models. We consider three different X-ray luminosity ($L_{X}$) - ultraviolet luminosity ($L_{UV}$) relations of quasars, i.e., the standard $L_{X}$-$L_{UV}$ relation and two redshift-evolutionary relations (Type I and Type II) respectively constructed from copula and considering a redshift correction… ▽ More In this paper, we use quasars calibrated from type Ia supernova (SN Ia) to constrain cosmological models. We consider three different X-ray luminosity ($L_{X}$) - ultraviolet luminosity ($L_{UV}$) relations of quasars, i.e., the standard $L_{X}$-$L_{UV}$ relation and two redshift-evolutionary relations (Type I and Type II) respectively constructed from copula and considering a redshift correction to the luminosity of quasars. Only in the case of the Type I relation, quasars can always provide effective constraints on the $Λ$CDM model. Furthermore, we show that, when the observational Hubble data (OHD) are added, the constraints on the absolute magnitude $M$ of SN Ia and the Hubble constant $H_0$ can be obtained. In the $Λ$CDM model, the OHD measurements plus quasars with the Type I relation yields $M$ =$-19.321^{+0.085}_{-0.076}$, which is in good agreement with the measurement from SH0ES ($M=-19.253\pm{0.027}$), and $H_0$ = $70.80\pm3.6~\mathrm{km~s^{-1}Mpc^{-1}}$, falling between the measurements from SH0ES and the Planck cosmic microwave background radiation data. △ Less

Submitted 13 June, 2024; v1 submitted 24 April, 2024; originally announced April 2024.

Comments: 9 pages,7 figures,4 tables. A new reference added. Published in MNRAS

Journal ref: MNRAS 530, 4493 (2024)

LIGO operates with quantum noise below the Standard Quantum Limit

Authors: Wenxuan Jia, Victoria Xu, Kevin Kuns, Masayuki Nakano, Lisa Barsotti, Matthew Evans, Nergis Mavalvala, Rich Abbott, Ibrahim Abouelfettouh, Rana Adhikari, Alena Ananyeva, Stephen Appert, Koji Arai, Naoki Aritomi, Stuart Aston, Matthew Ball, Stefan Ballmer, David Barker, Beverly Berger, Joseph Betzwieser, Dripta Bhattacharjee, Garilynn Billingsley, Nina Bode, Edgard Bonilla, Vladimir Bossilkov , et al. (146 additional authors not shown)

Abstract: Precision measurements of space and time, like those made by the detectors of the Laser Interferometer Gravitational-wave Observatory (LIGO), are often confronted with fundamental limitations imposed by quantum mechanics. The Heisenberg uncertainty principle dictates that the position and momentum of an object cannot both be precisely measured, giving rise to an apparent limitation called the Stan… ▽ More Precision measurements of space and time, like those made by the detectors of the Laser Interferometer Gravitational-wave Observatory (LIGO), are often confronted with fundamental limitations imposed by quantum mechanics. The Heisenberg uncertainty principle dictates that the position and momentum of an object cannot both be precisely measured, giving rise to an apparent limitation called the Standard Quantum Limit (SQL). Reducing quantum noise below the SQL in gravitational-wave detectors, where photons are used to continuously measure the positions of freely falling mirrors, has been an active area of research for decades. Here we show how the LIGO A+ upgrade reduced the detectors' quantum noise below the SQL by up to 3 dB while achieving a broadband sensitivity improvement, more than two decades after this possibility was first presented. △ Less

Submitted 22 April, 2024; originally announced April 2024.

Report number: LIGO-P2400059

TOI-837b is a Young Saturn-sized Exoplanet with a Massive 70 $M_{\oplus}$ Core

Authors: Oscar Barragán, Haochuan Yu, Alix Violet Freckelton, Annabella Meech, Michael Cretignier, Annelies Mortier, Suzanne Aigrain, Baptiste Klein, Niamh K. O'Sullivan, Edward Gillen, Louise Dyregaard Nielsen, Manuel Mallorquín, Norbert Zicher

Abstract: We present an exhaustive photometric and spectroscopic analysis of TOI-837, a F9/G0 35 Myr young star, hosting a transiting exoplanet, TOI-837b, with an orbital period of 8.32d. Utilising data from TESS and ground-based observations, we determine a planetary radius of 0.82 R_J for TOI-837b. Through detailed HARPS spectroscopic time series analysis, we derive a Doppler semi-amplitude of 35 m/s, cor… ▽ More We present an exhaustive photometric and spectroscopic analysis of TOI-837, a F9/G0 35 Myr young star, hosting a transiting exoplanet, TOI-837b, with an orbital period of 8.32d. Utilising data from TESS and ground-based observations, we determine a planetary radius of 0.82 R_J for TOI-837b. Through detailed HARPS spectroscopic time series analysis, we derive a Doppler semi-amplitude of 35 m/s, corresponding to a planetary mass of 0.39 M_J. The derived planetary properties suggest a substantial core of approximately 70 M_E, constituting about 60% of the planet's total mass. This finding poses a significant challenge to existing theoretical models of core formation. We propose that future atmospheric observations with JWST could provide insights into resolving ambiguities of TOI-837b, offering new perspectives on its composition, formation, and evolution. △ Less

Submitted 24 June, 2024; v1 submitted 21 April, 2024; originally announced April 2024.

Comments: MNRAS published version, MNRAS, 531, 4275

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.

Observation of Gravitational Waves from the Coalescence of a $2.5\text{-}4.5~M_\odot$ Compact Object and a Neutron Star

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, R. Abbott, I. Abouelfettouh, F. Acernese, K. Ackley, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, D. Agarwal, M. Agathos, M. Aghaei Abchouyeh, O. D. Aguiar, I. Aguilar, L. Aiello, A. Ain, P. Ajith, S. Akçay, T. Akutsu, S. Albanesi, R. A. Alfaidi, A. Al-Jodah , et al. (1771 additional authors not shown)

Abstract: We report the observation of a coalescing compact binary with component masses $2.5\text{-}4.5~M_\odot$ and $1.2\text{-}2.0~M_\odot$ (all measurements quoted at the 90% credible level). The gravitational-wave signal GW230529_181500 was observed during the fourth observing run of the LIGO-Virgo-KAGRA detector network on 2023 May 29 by the LIGO Livingston Observatory. The primary component of the so… ▽ More We report the observation of a coalescing compact binary with component masses $2.5\text{-}4.5~M_\odot$ and $1.2\text{-}2.0~M_\odot$ (all measurements quoted at the 90% credible level). The gravitational-wave signal GW230529_181500 was observed during the fourth observing run of the LIGO-Virgo-KAGRA detector network on 2023 May 29 by the LIGO Livingston Observatory. The primary component of the source has a mass less than $5~M_\odot$ at 99% credibility. We cannot definitively determine from gravitational-wave data alone whether either component of the source is a neutron star or a black hole. However, given existing estimates of the maximum neutron star mass, we find the most probable interpretation of the source to be the coalescence of a neutron star with a black hole that has a mass between the most massive neutron stars and the least massive black holes observed in the Galaxy. We provisionally estimate a merger rate density of $55^{+127}_{-47}~\text{Gpc}^{-3}\,\text{yr}^{-1}$ for compact binary coalescences with properties similar to the source of GW230529_181500; assuming that the source is a neutron star-black hole merger, GW230529_181500-like sources constitute about 60% of the total merger rate inferred for neutron star-black hole coalescences. The discovery of this system implies an increase in the expected rate of neutron star-black hole mergers with electromagnetic counterparts and provides further evidence for compact objects existing within the purported lower mass gap. △ Less

Submitted 26 July, 2024; v1 submitted 5 April, 2024; originally announced April 2024.

Comments: 45 pages (10 pages author list, 13 pages main text, 1 page acknowledgements, 13 pages appendices, 8 pages bibliography), 17 figures, 16 tables. Update to match version published in The Astrophysical Journal Letters. Data products available from https://zenodo.org/records/10845779

Report number: LIGO-P2300352

Journal ref: ApJL 970, L34 (2024)

SASHIMI-SIDM: Semi-analytical subhalo modelling for self-interacting dark matter at sub-galactic scales

Authors: Shin'ichiro Ando, Shunichi Horigome, Ethan O. Nadler, Daneng Yang, Hai-Bo Yu

Abstract: We combine the semi-analytical structure formation model, SASHIMI, which predicts subhalo populations in collisionless, cold dark matter (CDM), with a parametric model that maps CDM halos to self-interacting dark matter (SIDM) halos. The resulting model, SASHIMI-SIDM, generates SIDM subhalo populations down to sub-galactic mass scales, for an arbitrary input cross section, in minutes. We show that… ▽ More We combine the semi-analytical structure formation model, SASHIMI, which predicts subhalo populations in collisionless, cold dark matter (CDM), with a parametric model that maps CDM halos to self-interacting dark matter (SIDM) halos. The resulting model, SASHIMI-SIDM, generates SIDM subhalo populations down to sub-galactic mass scales, for an arbitrary input cross section, in minutes. We show that SASHIMI-SIDM agrees with SIDM subhalo populations from high-resolution cosmological zoom-in simulations in resolved regimes. Crucially, we predict that the fraction of core-collapsed subhalos peaks at a mass scale determined by the input SIDM cross section and decreases toward higher halo masses, consistent with the predictions of gravothermal models and cosmological simulations. For the first time, we also show that the core-collapsed fraction decreases toward lower halo masses; this result is uniquely enabled by our semi-analytical approach. As a proof of principle, we apply SASHIMI-SIDM to predict the boost to the local dark matter density and annihilation rate from core-collapsed SIDM subhalos, which can be enhanced relative to CDM by an order of magnitude for viable SIDM models. Thus, SASHIMI-SIDM provides an efficient and reliable tool for scanning SIDM parameter space and testing it with astrophysical observations. The code is publicly available at https://github.com/shinichiroando/sashimi-si. △ Less

Submitted 25 March, 2024; originally announced March 2024.

Comments: 21 pages, 7 figures. The code accompanying this submission is available at https://github.com/shinichiroando/sashimi-si

Almost Optically Dark Galaxies in DECaLS (I): Detection, Optical Properties and Possible Origins

Authors: Lin Du, Wei Du, Cheng Cheng, Ming Zhu, Haiyang Yu, Hong Wu

Abstract: We report the discovery of eight optical counterparts of ALFALFA extragalactic objects from DECaLS, five of which are discovered for the first time. These objects were flagged as HI emission sources with no optical counterparts in SDSS before. Multi-band data reveal their unusual physical properties. They are faint and blue ($g-r=-0.35\sim0.55$), with quite low surface brightness (… ▽ More We report the discovery of eight optical counterparts of ALFALFA extragalactic objects from DECaLS, five of which are discovered for the first time. These objects were flagged as HI emission sources with no optical counterparts in SDSS before. Multi-band data reveal their unusual physical properties. They are faint and blue ($g-r=-0.35\sim0.55$), with quite low surface brightness ($μ_{\rm g,peak}=24.88\sim26.41\,{\rm mag}/{\rm arcsec}^2$), irregular morphologies, low stellar masses ($log_{10}(M_{*}/M_\odot)=5.27\sim7.15$), low star formation rates ($SFR=0.21\sim9.24\times10^{-3}\,{M_\odot}\,{\rm yr}^{-1}$), and remarkably high HI-to-stellar mass ratios ($log_{10}(M_{\rm HI}/M_{*}) = 1.72\sim3.22$, except AGC\,215415). They deviate from the scaling relations between HI and optical properties defined by the ALFALFA sample and the baryonic Tully-Fisher relation. They agree well with the main sequence of star-forming galaxies but exhibit low star-forming efficiency. Based on their physical properties and environments, we speculate that six of these objects may have originated from tidal processes, while the remaining two appear to have isolated origins. They may have had a relatively calm evolutionary history and only begun to form stars recently. △ Less

Submitted 18 March, 2024; originally announced March 2024.

Comments: 32 pages, 11 figures, accepted by the Astrophysical Journal

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)

Ultralight vector dark matter search using data from the KAGRA O3GK run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, R. Abbott, H. Abe, I. Abouelfettouh, F. Acernese, K. Ackley, C. Adamcewicz, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, I. Aguilar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi , et al. (1778 additional authors not shown)

Abstract: Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we prese… ▽ More Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we present the result of a search for $U(1)_{B-L}$ gauge boson DM using the KAGRA data from auxiliary length channels during the first joint observation run together with GEO600. By applying our search pipeline, which takes into account the stochastic nature of ultralight DM, upper bounds on the coupling strength between the $U(1)_{B-L}$ gauge boson and ordinary matter are obtained for a range of DM masses. While our constraints are less stringent than those derived from previous experiments, this study demonstrates the applicability of our method to the lower-mass vector DM search, which is made difficult in this measurement by the short observation time compared to the auto-correlation time scale of DM. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Comments: 20 pages, 5 figures

Report number: LIGO-P2300250

Numerical challenges for energy conservation in N-body simulations of collapsing self-interacting dark matter haloes

Authors: Moritz S. Fischer, Klaus Dolag, Hai-Bo Yu

Abstract: Dark matter (DM) haloes can be subject to gravothermal collapse if DM is not collisionless but has strong self-interactions. When the scattering can efficiently transfer heat from the centre to the outskirts, the central region of the halo collapses and reaches densities much higher than those for collisionless DM. This phenomenon is potentially observable in studies of strong lensing. Current the… ▽ More Dark matter (DM) haloes can be subject to gravothermal collapse if DM is not collisionless but has strong self-interactions. When the scattering can efficiently transfer heat from the centre to the outskirts, the central region of the halo collapses and reaches densities much higher than those for collisionless DM. This phenomenon is potentially observable in studies of strong lensing. Current theoretical efforts are motivated by observations of surprisingly dense substructures. A comparison with them requires accurate predictions. One method to obtain such predictions is to use N-body simulations. The collapsed haloes are extreme systems that pose severe challenges to state-of-the-art codes used to model self-interacting dark matter (SIDM). We investigate the root of such problems with a focus on energy non-conservation. We run N-body simulations with and without DM self-interactions of an isolated DM-only halo and change numerical parameters relevant to the accuracy of the simulation. We find that not only the numerical scheme for the DM self-interactions can lead to energy non-conservation but also the modelling of gravitational interaction and the time integration are problematic. The issues we found are: (a) particles changing their time step in a non-time-reversible manner; (b) the asymmetry in the tree-based gravitational force evaluation; (c) SIDM velocity kicks break the symplectic nature and time symmetry. Tuning the parameters of the simulation to achieve a high accuracy allows for conserving energy not only at early stages of the evolution but also at late ones. However, the cost of the simulations becomes prohibitively large. Some problems making the simulations of the gravothermal collapse phase inaccurate, can be overcome by choosing appropriate numerical schemes. However, others remain challenging. Our findings motivate further work on addressing these challenges. △ Less

Submitted 1 March, 2024; originally announced March 2024.

Comments: 10 pages, 6 figures

CDM and SIDM Interpretations of the Strong Gravitational Lensing Object JWST-ER1

Authors: Demao Kong, Daneng Yang, Hai-Bo Yu

Abstract: van Dokkum et al. (arXiv:2309.07969) reported the discovery of JWST-ER1, a strong lensing object at redshift $z\approx2$, using data from the James Webb Space Telescope. The lens mass within the Einstein ring is $5.9$ times higher than the expected stellar mass from a Chabrier initial mass function, indicating a high dark matter density. In this work, we show that a cold dark matter halo, influenc… ▽ More van Dokkum et al. (arXiv:2309.07969) reported the discovery of JWST-ER1, a strong lensing object at redshift $z\approx2$, using data from the James Webb Space Telescope. The lens mass within the Einstein ring is $5.9$ times higher than the expected stellar mass from a Chabrier initial mass function, indicating a high dark matter density. In this work, we show that a cold dark matter halo, influenced by gas-driven adiabatic contraction, can account for the observed lens mass. We interpret the measurement of JWST-ER1 in the self-interacting dark matter scenario and show that the cross section per particle mass $σ/m\approx0.1~{\rm cm^2/g}$ is generally favored. Intriguingly, $σ/m\approx0.1~{\rm cm^2/g}$ can also be consistent with the strong lensing observations of early-type galaxies at redshift $z\approx0.2$, where adiabatic contraction is not observed overall. △ Less

Submitted 1 April, 2024; v1 submitted 24 February, 2024; originally announced February 2024.

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

The boundary of cosmic filaments

Authors: Wei Wang, Peng Wang, Hong Guo, Xi Kang, Noam I. Libeskind, Daniela Galarraga-Espinosa, Volker Springel, Rahul Kannan, Lars Hernquist, Rudiger Pakmor, Haoran Yu, Sownak Bose, Quan Guo, Luo Yu, Cesar Hernandez-Aguayo

Abstract: For decades, the boundary of cosmic filaments have been a subject of debate. In this work, we determine the physically-motivated radii of filaments by constructing stacked galaxy number density profiles around the filament spines. We find that the slope of the profile changes with distance to the filament spine, reaching its minimum at approximately 1 Mpc at z = 0 in both state-of-the-art hydrodyn… ▽ More For decades, the boundary of cosmic filaments have been a subject of debate. In this work, we determine the physically-motivated radii of filaments by constructing stacked galaxy number density profiles around the filament spines. We find that the slope of the profile changes with distance to the filament spine, reaching its minimum at approximately 1 Mpc at z = 0 in both state-of-the-art hydrodynamical simulations and observational data. This can be taken as the average value of the filament radius. Furthermore, we note that the average filament radius rapidly decreases from z = 4 to z = 1, and then slightly increases. Moreover, we find that the filament radius depends on the filament length, the distance from connected clusters, and the masses of the clusters. These results suggest a two-phase formation scenario of cosmic filaments. The filaments experience rapid contraction before z = 1, but their density distribution has remained roughly stable since then. The subsequent mass transport along the filaments to the connected clusters is likely to have contributed to the formation of the clusters themselves. △ Less

Submitted 3 August, 2024; v1 submitted 18 February, 2024; originally announced February 2024.

Comments: 12 pages, 11 figures, accepted by mnras

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

Resonant amplification of curvature perturbations in inflation model with periodical derivative coupling

Authors: Li-Yang Chen, Hongwei Yu, Puxun Wu

Abstract: In this paper, we introduce a weak, transient and periodical derivative coupling between the inflaton field and gravity, and find that the square of the sound speed of the curvature perturbations becomes a periodic function, which results in that the equation of the curvature perturbations can be transformed into the form of the Mathieu equation in the sub-horizon limit. Thus, the parametric reson… ▽ More In this paper, we introduce a weak, transient and periodical derivative coupling between the inflaton field and gravity, and find that the square of the sound speed of the curvature perturbations becomes a periodic function, which results in that the equation of the curvature perturbations can be transformed into the form of the Mathieu equation in the sub-horizon limit. Thus, the parametric resonance will amplify the curvature perturbations so as to generate a formation of abundant primordial black holes (PBHs). We show that the generated PBHs can make up most of dark matter. Associated with the generation of PBHs, the large scalar perturbations will give rise to the scalar induced gravitational waves which may be detected by future gravitational wave projects. △ Less

Submitted 27 January, 2024; v1 submitted 15 January, 2024; originally announced January 2024.

Comments: 19 pages, 4 figures, refs added

Journal ref: Phys. Lett. B 849 (2024) 138457

Modelling stellar variability in archival HARPS data: I -- Rotation and activity properties with multi-dimensional Gaussian Processes

Authors: Haochuan Yu, Suzanne Aigrain, Baptiste Klein, Oscar Barragán, Annelies Mortier, Niamh K. O'Sullivan, Michael Cretignier

Abstract: Although instruments for measuring the radial velocities (RVs) of stars now routinely reach sub-meter per second accuracy, the detection of low-mass planets is still very challenging. The rotational modulation and evolution of spots and/or faculae can induce variations in the RVs at the level of a few m/s in Sun-like stars. To overcome this, a multi-dimensional Gaussian Process framework has been… ▽ More Although instruments for measuring the radial velocities (RVs) of stars now routinely reach sub-meter per second accuracy, the detection of low-mass planets is still very challenging. The rotational modulation and evolution of spots and/or faculae can induce variations in the RVs at the level of a few m/s in Sun-like stars. To overcome this, a multi-dimensional Gaussian Process framework has been developed to model the stellar activity signal using spectroscopic activity indicators together with the RVs. A recently published computationally efficient implementation of this framework, S+LEAF 2, enables the rapid analysis of large samples of targets with sizeable data sets. In this work, we apply this framework to HARPS observations of 268 well-observed targets with precisely determined stellar parameters. Our long-term goal is to quantify the effectiveness of this framework to model and mitigate activity signals for stars of different spectral types and activity levels. In this first paper in the series, we initially focus on the activity indicators (S-index and Bisector Inverse Slope), and use them to a) measure rotation periods for 49 slow rotators in our sample, b) explore the impact of these results on the spin-down of middle-aged late F, G & K stars, and c) explore indirectly how the spot to facular ratio varies across our sample. Our results should provide valuable clues for planning future RV planet surveys such as the Terra Hunting Experiment or the PLATO ground-based follow-up observations program, and help fine-tune current stellar structure and evolution models. △ Less

Submitted 10 January, 2024; originally announced January 2024.

Comments: Accepted for publication in MNRAS

Self-interacting dark matter interpretation of Crater II

Authors: Xingyu Zhang, Hai-Bo Yu, Daneng Yang, Haipeng An

Abstract: The satellite galaxy Crater II of the Milky Way is extremely cold and exceptionally diffuse. These unusual properties are challenging to understand in the standard model of cold dark matter. We investigate the formation of Crater II in self-interacting dark matter (SIDM), where dark matter particles can scatter and thermalize. We conduct a series of controlled N-body simulations to model the tidal… ▽ More The satellite galaxy Crater II of the Milky Way is extremely cold and exceptionally diffuse. These unusual properties are challenging to understand in the standard model of cold dark matter. We investigate the formation of Crater II in self-interacting dark matter (SIDM), where dark matter particles can scatter and thermalize. We conduct a series of controlled N-body simulations to model the tidal evolution of Crater II, varying the self-interacting cross section, orbit parameters, and initial stellar distribution. Dark matter self-interactions lead to halo core formation and the distribution of stars expands accordingly. A cored SIDM halo also boosts tidal mass loss, allowing for a high orbit. Our simulations show that SIDM halos with a $1~{\rm kpc}$ core can simultaneously explain the low stellar velocity dispersion and the large half-light radius of Crater II, remaining robust to the initial distribution of stars. For the orbit motivated by the measurements from Gaia Early Data Release 3, the favored self-interacting cross section is approximately $60~{\rm cm^2/g}$ on the mass scale of Crater II. △ Less

Submitted 5 June, 2024; v1 submitted 10 January, 2024; originally announced January 2024.

Comments: 9 pages, 5 figures; Revised for publication in The Astrophysical Journal Letters

Observations favor the redshift-evolutionary $L_X$-$L_{UV}$ relation of quasars from copula

Authors: Bao Wang, Yang Liu, Hongwei Yu, Puxun Wu

Abstract: We compare, with data from the quasars, the Hubble parameter measurements, and the Pantheon+ type Ia supernova, three different relations between X-ray luminosity ($L_X$) and ultraviolet luminosity ($L_{UV}$) of quasars. These three relations consist of the standard and two redshift-evolutionary $L_X$-$L_{UV}$ relations which are constructed respectively by considering a redshift dependent correct… ▽ More We compare, with data from the quasars, the Hubble parameter measurements, and the Pantheon+ type Ia supernova, three different relations between X-ray luminosity ($L_X$) and ultraviolet luminosity ($L_{UV}$) of quasars. These three relations consist of the standard and two redshift-evolutionary $L_X$-$L_{UV}$ relations which are constructed respectively by considering a redshift dependent correction to the luminosities of quasars and using the statistical tool called copula. By employing the PAge approximation for a cosmological-model-independent description of the cosmic background evolution and dividing the quasar data into the low-redshift and high-redshift parts, we find that the constraints on the PAge parameters from the low-redshift and high-redshift data, which are obtained with the redshift-evolutionary relations, are consistent with each other, while they are not when the standard relation is considered. If the data are used to constrain the coefficients of the relations and the PAge parameters simultaneously, then the observations support the redshift-evolutionary relations at more than $3σ$. The Akaike and Bayes information criteria indicate that there is strong evidence against the standard relation and mild evidence against the redshift-evolutionary relation constructed by considering a redshift dependent correction to the luminosities of quasars. This suggests that the redshift-evolutionary $L_X$-$L_{UV}$ relation of quasars constructed from copula is favored by the observations. △ Less

Submitted 2 January, 2024; originally announced January 2024.

Comments: 15 pages, 5 figures. Accepted for publication in ApJ

The FAST all sky HI survey (FASHI): The first release of catalog

Authors: Chuan-Peng Zhang, M. Zhu, P. Jiang, C. Cheng, J. Wang, J. Wang, J. -L. Xu, X. -L. Liu, N. -P. Yu, L. Qian, H. Yu, M. Ai, Y. Jing, C. Xu, Z. Liu, X. Guan, C. Sun, Q. Yang, M. Huang, Q. Hao, FAST Collaboration

Abstract: The FAST All Sky HI survey (FASHI) was designed to cover the entire sky observable by the Five-hundred-meter Aperture Spherical radio Telescope (FAST), spanning approximately 22000 square degrees of declination between -14 deg and +66 deg, and in the frequency range of 1050-1450 MHz, with the expectation of eventually detecting more than 100000 HI sources. Between August 2020 and June 2023, FASHI… ▽ More The FAST All Sky HI survey (FASHI) was designed to cover the entire sky observable by the Five-hundred-meter Aperture Spherical radio Telescope (FAST), spanning approximately 22000 square degrees of declination between -14 deg and +66 deg, and in the frequency range of 1050-1450 MHz, with the expectation of eventually detecting more than 100000 HI sources. Between August 2020 and June 2023, FASHI had covered more than 7600 square degrees, which is approximately 35% of the total sky observable by FAST. It has a median detection sensitivity of around 0.76 mJy/beam and a spectral line velocity resolution of ~6.4 km/s at a frequency of ~1.4 GHz. As of now, a total of 41741 extragalactic HI sources have been detected in the frequency range 1305.5-1419.5 MHz, corresponding to a redshift limit of z<0.09. By cross-matching FASHI sources with the Siena Galaxy Atlas (SGA) and the Sloan Digital Sky Survey (SDSS) catalogs, we found that 16972 (40.7%) sources have spectroscopic redshifts and 10975 (26.3%) sources have only photometric redshifts. Most of the remaining 13794 (33.0%) HI sources are located in the direction of the Galactic plane, making their optical counterparts difficult to identify due to high extinction or high contamination of Galactic stellar sources. Based on current survey results, the FASHI survey is an unprecedented blind extragalactic HI survey. It has higher spectral and spatial resolution and broader coverage than the Arecibo Legacy Fast ALFA Survey (ALFALFA). When completed, FASHI will provide the largest extragalactic HI catalog and an objective view of HI content and large-scale structure in the local universe. △ Less

Submitted 10 December, 2023; originally announced December 2023.

Comments: 22 pages, 12 figures, published in SCPMA. All catalogs are available at https://zcp521.github.io/fashi and https://fast.bao.ac.cn/cms/article/271/

Journal ref: Sci. China-Phys. Mech. Astron. 67, 219511 (2024)

NANOGrav hints for first-order confinement-deconfinement phase transition in different QCD-matter scenarios

Authors: Zu-Cheng Chen, Shou-Long Li, Puxun Wu, Hongwei Yu

Abstract: Recent observations from several pulsar timing array (PTA) collaborations have unveiled compelling evidence for a stochastic signal in the nanohertz band. This signal aligns remarkably with a gravitational wave (GW) background, potentially originating from the first-order color charge confinement phase transition. Distinct quantum chromodynamics (QCD) matters, such as quarks or gluons, and diverse… ▽ More Recent observations from several pulsar timing array (PTA) collaborations have unveiled compelling evidence for a stochastic signal in the nanohertz band. This signal aligns remarkably with a gravitational wave (GW) background, potentially originating from the first-order color charge confinement phase transition. Distinct quantum chromodynamics (QCD) matters, such as quarks or gluons, and diverse phase transition processes thereof can yield disparate GW energy density spectra. In this paper, employing the Bayesian analysis on the NANOGrav 15-year data set, we explore the compatibility with the observed PTA signal of the GW from phase transitions of various QCD matter scenarios in the framework of the holographic QCD. We find that the PTA signal can be effectively explained by the GW from the confinement-deconfinement phase transition of pure quark systems in a hard wall model of the holographic QCD where the bubble dynamics, one important source of the GWs, is of the Jouguet detonations. Notably, our analysis decisively rules out the plausibility of the pure gluon QCD-matter scenario and the non-runaway bubble dynamics model for the phase transition in explaining the observed PTA signal. △ Less

Submitted 22 January, 2024; v1 submitted 4 December, 2023; originally announced December 2023.

Comments: 7 pages, 2 figures; accepted by PRD

Constraints on cosmic star formation history via a new modeling of the radio luminosity function of star-forming galaxies

Authors: Wenjie Wang, Zunli Yuan, Hongwei Yu, Jirong Mao

Abstract: Radio wavelengths offer a unique possibility to trace the total star-formation rate (SFR) in galaxies, both obscured and unobscured. To probe the dust-unbiased star-formation history, an accurate measurement of the radio luminosity function (LF) for star-forming galaxies (SFGs) is crucial. We make use of an SFG sample (5900 sources) from the Very Large Array (VLA) COSMOS 3 GHz data to perform a ne… ▽ More Radio wavelengths offer a unique possibility to trace the total star-formation rate (SFR) in galaxies, both obscured and unobscured. To probe the dust-unbiased star-formation history, an accurate measurement of the radio luminosity function (LF) for star-forming galaxies (SFGs) is crucial. We make use of an SFG sample (5900 sources) from the Very Large Array (VLA) COSMOS 3 GHz data to perform a new modeling of the radio LF. By integrating the analytical LF, we aim to calculate the history of the cosmic SFR density (SFRD) from $z\sim5$ onwards. For the first time, we use both models of the pure luminosity evolution (PLE) and joint luminosity+density evolution (LADE) to fit the LFs directly to the radio data using a full maximum-likelihood analysis, considering the sample completeness correction. We also incorporate updated observations of local radio LFs and radio source counts into the fitting process to obtain additional constraints. We find that the PLE model cannot be used to describe the evolution of the radio LF at high redshift ($z>2$). By construct, our LADE models can successfully fit a large amount of data on radio LFs and source counts of SFGs from recent observations. We therefore conclude that density evolution is genuinely indispensable in modeling the evolution of SFG radio LFs. Our SFRD curve shows a good fit to the SFRD points derived by previous radio estimates. In view of the fact that our radio LFs are not biased, as opposed those of previous studies performed by fitting the $1/V_{\rm max}$ LF points, our SFRD results should be an improvement on these previous estimates. Below $z\sim1.5$, our SFRD matches a published multiwavelength compilation, while our SFRD turns over at a slightly higher redshift ($2<z<2.5$) and falls more rapidly out to high redshift. △ Less

Submitted 22 January, 2024; v1 submitted 15 November, 2023; originally announced November 2023.

Comments: 14 pages, 10 figures, final version, to be published in A&A

Spin speed correlations and the evolution of galaxy-halo systems

Authors: Ming-Jie Sheng, Lin Zhu, Hao-Ran Yu, Hongchuan Ma, Haikun Li, Peng Wang, Xi Kang

Abstract: Galaxy angular momenta (spins) contain valuable cosmological information, complementing their positions and velocities. The baryonic spin direction of galaxies has been probed as a reliable tracer of their host halos and the primordial spin modes. Here we use the TNG100 simulation of the IllustrisTNG project to study the spin magnitude correlations between dark matter, gas, and stellar components… ▽ More Galaxy angular momenta (spins) contain valuable cosmological information, complementing their positions and velocities. The baryonic spin direction of galaxies has been probed as a reliable tracer of their host halos and the primordial spin modes. Here we use the TNG100 simulation of the IllustrisTNG project to study the spin magnitude correlations between dark matter, gas, and stellar components of galaxy-halo systems and their evolutions across cosmic history. We find that these components generate similar initial spin magnitudes from the same tidal torque in Lagrangian space. At low redshifts, the gas component still traces the spin magnitude of the dark matter halo and the primordial spin magnitude. However, the traceability of the stellar component depends on the $ex$ $situ$ stellar mass fraction, $f_{\rm acc}$. Our results suggest that the galaxy baryonic spin magnitude can also serve as a tracer of their host halo and the initial perturbations, and the galaxy-halo correlations are affected by the similarity of their evolution histories. △ Less

Submitted 2 July, 2024; v1 submitted 14 November, 2023; originally announced November 2023.

Comments: 10 pages, 10 figures. Matches the accepted version in Physical Review D

Journal ref: Physical Review D, 109, 123548, 2024

Alpha backgrounds in NaI(Tl) crystals of COSINE-100

Authors: G. Adhikari, N. Carlin, D. F. F. S. Cavalcante, J. Y. Cho, J. J. Choi, S. Choi, A. C. Ezeribe, L. E. Franca, C. Ha, I. S. Hahn, S. J. Hollick, E. J. Jeon, H. W. Joo, W. G. Kang, M. Kauer, B. H. Kim, H. J. Kim, J. Kim, K. W. Kim, S. H. Kim, S. K. Kim, S. W. Kim, W. K. Kim, Y. D. Kim, Y. H. Kim , et al. (38 additional authors not shown)

Abstract: COSINE-100 is a dark matter direct detection experiment with 106 kg NaI(Tl) as the target material. 210Pb and daughter isotopes are a dominant background in the WIMP region of interest and are detected via beta decay and alpha decay. Analysis of the alpha channel complements the background model as observed in the beta/gamma channel. We present the measurement of the quenching factors and Monte Ca… ▽ More COSINE-100 is a dark matter direct detection experiment with 106 kg NaI(Tl) as the target material. 210Pb and daughter isotopes are a dominant background in the WIMP region of interest and are detected via beta decay and alpha decay. Analysis of the alpha channel complements the background model as observed in the beta/gamma channel. We present the measurement of the quenching factors and Monte Carlo simulation results and activity quantification of the alpha decay components of the COSINE-100 NaI(Tl) crystals. The data strongly indicate that the alpha decays probabilistically undergo two possible quenching factors but require further investigation. The fitted results are consistent with independent measurements and improve the overall understanding of the COSINE-100 backgrounds. Furthermore, the half-life of 216Po has been measured to be 143.4 +/- 1.2 ms, which is consistent with and more precise than recent measurements. △ Less

Submitted 30 January, 2024; v1 submitted 8 November, 2023; originally announced November 2023.

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

Measuring Supermassive Black Hole Properties via Gravitational Radiation from Eccentrically Orbiting Stellar Mass Black Hole Binaries

Authors: Andrew Laeuger, Brian Seymour, Yanbei Chen, Hang Yu

Abstract: There may exist stellar-mass binary black holes (BBH) which merge while orbiting nearby a supermassive black hole (SMBH). In such a triple system, the SMBH will modulate the gravitational waveform of the BBH through orbital Doppler shift and de Sitter precession of the angular momentum. Future space-based GW observatories focused on the milli- and decihertz band will be uniquely poised to observe… ▽ More There may exist stellar-mass binary black holes (BBH) which merge while orbiting nearby a supermassive black hole (SMBH). In such a triple system, the SMBH will modulate the gravitational waveform of the BBH through orbital Doppler shift and de Sitter precession of the angular momentum. Future space-based GW observatories focused on the milli- and decihertz band will be uniquely poised to observe these waveform modulations, as the GW frequency from stellar-mass BBHs varies slowly in this band while modulation effects accumulate. In this work, we apply the Fisher information matrix formalism to estimate how well space-borne GW detectors can measure properties of BBH+SMBH hierarchical triples using the GW from orbiting BBH. We extend previous work by considering the more realistic case of an eccentric orbit around the SMBH, and notably include the effects of orbital pericenter precession. We find that for detector concepts such as LISA, B-DECIGO, and TianGO, we can extract the SMBH mass and semimajor axis of the orbit with a fractional uncertainty below the 0.1% level over a wide range of triple system parameters. Furthermore, we find that the effects of pericenter precession and orbital eccentricity significantly improve our ability to measure this system. We also find that while LISA could measure these systems, the decihertz detector concepts B-DECIGO and TianGO would enable better sensitivity to the triple's parameters. △ Less

Submitted 25 October, 2023; originally announced October 2023.

Comments: 12 pages (main text excluding references and appendices), 11 figures, submitted to PRD

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