Showing 1–50 of 110 results for author: Hu, Q

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

Prospects for the detection of very-high-energy pulsars with LHAASO and SWGO

Authors: Quan Hu, Yi Zhang, Kaikai Duan, Houdun Zeng

Abstract: Pulsations from the Crab pulsar have been detected by the MAGIC telescopes at energies up to 1.5 TeV, and the pulsed emission from the Vela pulsar was detected by H.E.S.S., reaching tens of TeV. These discoveries, along with the proposed additional emission due to inverse Compton scattering at TeV energies, lead us to consider suitable candidates for detection with current and future extensive air… ▽ More Pulsations from the Crab pulsar have been detected by the MAGIC telescopes at energies up to 1.5 TeV, and the pulsed emission from the Vela pulsar was detected by H.E.S.S., reaching tens of TeV. These discoveries, along with the proposed additional emission due to inverse Compton scattering at TeV energies, lead us to consider suitable candidates for detection with current and future extensive air show (EAS) experiments at very-high-energy (VHE; 0.1 $-$ 100 TeV) ranges. Leveraging energy spectrum data from pulsars as observed by Fermi and Imaging Atmospheric Cherenkov Telescopes (IACTs) and considering the sensitivities of both LHAASO and SWGO, this study evaluates their detectability and estimates the time required for their significant detection. Our results indicate that LHAASO could detect the Crab's pulsed signal within six years, while SWGO might detect Vela's signal within one year. Observations of the most energetic Fermi pulsars with EAS experiments will provide insight into the nature of VHE pulsar emissions, helping to clarify the primary characteristics of VHE pulsars. △ Less

Submitted 28 June, 2024; originally announced July 2024.

Comments: 6 pages, 3 figures, 2 Tables and accepted for publication in MNRAS

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

GWnext 2024: Meeting Summary

Authors: Alejandro Torres-Orjuela, Veronica Vazquez-Aceves, Rui Xu, Jin-Hong Chen, Andrea Derdzinski, Matthias U. Kruckow, Stefano Rinaldi, Lorenzo Speri, Ziming Wang, Garvin Yim, Xue-Ting Zhang, Qian Hu, Miaoxin Liu, Xiangyu Lyu, Zheng Wu, Cong Zhou, Manuel Arca Sedda, Yan-Chen Bi, Hong-Yu Chen, Xian Chen, Jiageng Jiao, Yu-Mei Wu

Abstract: GWnext 2024 was a meeting held in the Kavli Institute for Astronomy and Astrophysics at Peking University in March $4^\text{th} - 8^\text{th}$, 2024. In the meeting researchers at different career stages -- with a particular focus on early career scientists -- working on the different aspects of gravitational wave (GW) astronomy gathered to discuss the current status as well as prospects of the fi… ▽ More GWnext 2024 was a meeting held in the Kavli Institute for Astronomy and Astrophysics at Peking University in March $4^\text{th} - 8^\text{th}$, 2024. In the meeting researchers at different career stages -- with a particular focus on early career scientists -- working on the different aspects of gravitational wave (GW) astronomy gathered to discuss the current status as well as prospects of the field. The meeting was divided into three core sessions: Astrophysics, GW Theory, and Detection. Each session consisted of introductory talks and extended discussion sessions. Moreover, there was a poster session where students could present their results. In this paper, we summarize the results presented during the meeting and present the most important outcomes. △ Less

Submitted 27 May, 2024; originally announced June 2024.

Single Aperture Large Telescope for Universe Studies (SALTUS): Science Overview

Authors: Gordon Chin, Carrie M. Anderson, Jennifer Bergner, Nicolas Biver, Gordon L. Bjoraker, Thibault Cavalie, Michael DiSanti, Jian-Rong Gao, Paul Hartogh, Leon K. Harding, Qing Hu, Daewook Kim, Craig Kulesa, Gert de Lange, David T. Leisawitz, Rebecca C. Levy, Arthur Lichtenberger, Daniel P. Marronh, Joan Najita, Trent Newswander, George H. Rieke, Dimitra Rigopoulou, Peter Roefsema, Nathan X. Roth, Kamber Schwarz , et al. (11 additional authors not shown)

Abstract: The SALTUS Probe mission will provide a powerful far-infrared (far-IR) pointed space observatory to explore our cosmic origins and the possibility of life elsewhere. The observatory employs an innovative deployable 14-m aperture, with a sunshield that will radiatively cool the off-axis primary to <45K. This cooled primary reflector works in tandem with cryogenic coherent and incoherent instruments… ▽ More The SALTUS Probe mission will provide a powerful far-infrared (far-IR) pointed space observatory to explore our cosmic origins and the possibility of life elsewhere. The observatory employs an innovative deployable 14-m aperture, with a sunshield that will radiatively cool the off-axis primary to <45K. This cooled primary reflector works in tandem with cryogenic coherent and incoherent instruments that span the 34 to 660 micron far-IR range at both high and moderate spectral resolutions. △ Less

Submitted 21 May, 2024; originally announced May 2024.

Comments: 49 pages, 10 figures, 3 tables, submitted to SPIE JATIS

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

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

Exploration of morphological coherence in open clusters with a "core-shell'' structure

Authors: Qingshun Hu, Yu Zhang, Songmei Qin, Jing Zhong, Li Chen, Yangping Luo

Abstract: The study of their morphological coherence allows us to obtain a better understanding of the morphological evolution of open clusters. We employed the ellipsoid fitting method to delineate the 3D spatial structure of the sample clusters while using the morphological dislocation (MD) defined in our previous work and the ellipticity ratio (ER) of the clusters' inner and outer structures to character… ▽ More The study of their morphological coherence allows us to obtain a better understanding of the morphological evolution of open clusters. We employed the ellipsoid fitting method to delineate the 3D spatial structure of the sample clusters while using the morphological dislocation (MD) defined in our previous work and the ellipticity ratio (ER) of the clusters' inner and outer structures to characterize the morphological coherence of the sample clusters. The results show an inverse correlation between the ER of the sample clusters and the number of their members, indicating that sample clusters with a much more elliptical external morphology than internal shape generally tend to host a large number of members. Meanwhile, a slight shrinking of the MD of the sample clusters with their members' number may shed light on the significant role of the gravitational binding of the sample clusters in maintaining their morphological stability. Moreover, there are no correlations between the MD and ER of the sample clusters and their age. They are also not significantly correlated with the X-axis, the Y-axis, their orbital eccentricities, and the radial and vertical forces on them. However, the ER of the sample clusters displays some fluctuations in the distributions between it and the above covariates, implying that the morphologies of the sample clusters are sensitive to the external environment if sample effects are not taken into account. Finally, the analysis of the 3D spatial shapes of sample clusters with a small ER or a large ER demonstrates that the number of members lays an important foundation for forming a dense internal system for sample clusters. At the same time, the MD of the sample clusters can serve well as an indicator of their morphological stability, which is built upon a certain amount of member stars. △ Less

Submitted 11 May, 2024; v1 submitted 9 May, 2024; originally announced May 2024.

Comments: 17 pages, 13 figures, 1 table. Accepted by Astronomy & Astrophysics

Journal ref: A&A 687, A291 (2024)

LHAASO-KM2A detector simulation using Geant4

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

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

Submitted 7 April, 2024; originally announced April 2024.

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

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

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

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

Comments: 8 pages, 3 figures

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

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

The PFS view of TOI-677 b: A spin-orbit aligned warm Jupiter in a dynamically hot system

Authors: Qingru Hu, Malena Rice, Xian-Yu Wang, Songhu Wang, Avi Shporer, Johanna K. Teske, Samuel W. Yee, R. Paul Butler, Stephen Shectman, Jeffrey D. Crane, Karen A. Collins, Kevin I. Collins

Abstract: TOI-677 b is part of an emerging class of ``tidally-detached'' gas giants ($a/R_\star \gtrsim 11$) that exhibit large orbital eccentricities and yet low stellar obliquities. Such sources pose a challenge for models of giant planet formation, which must account for the excitation of high eccentricities without large changes in the orbital inclination. In this work, we present a new Rossiter-McLaugh… ▽ More TOI-677 b is part of an emerging class of ``tidally-detached'' gas giants ($a/R_\star \gtrsim 11$) that exhibit large orbital eccentricities and yet low stellar obliquities. Such sources pose a challenge for models of giant planet formation, which must account for the excitation of high eccentricities without large changes in the orbital inclination. In this work, we present a new Rossiter-McLaughlin (RM) measurement for the tidally-detached warm Jupiter TOI-677 b, obtained using high-precision radial velocity observations from the PFS/Magellan spectrograph. Combined with previously published observations from the ESPRESSO/VLT spectrograph, we derive one of the most precisely constrained sky-projected spin-orbit angle measurements to date for an exoplanet. The combined fit offers a refined set of self-consistent parameters, including a low sky-projected stellar obliquity of $λ=3.2^{+1.6}_{-1.5}$ deg and a moderately high eccentricity of $e=0.460^{+0.019}_{-0.018}$, that further constrains the puzzling architecture of this system. We examine several potential scenarios that may have produced the current TOI-677 orbital configuration, ultimately concluding that TOI-677 b most likely had its eccentricity excited through disk-planet interactions. This system adds to a growing population of aligned warm Jupiters on eccentric orbits around hot ($T_{\rm eff}>6100$ K) stars. △ Less

Submitted 11 February, 2024; originally announced February 2024.

Comments: 16 pages, 4 figures, accepted for publication in the Astronomical Journal

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

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

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

Submitted 25 October, 2023; originally announced October 2023.

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

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

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

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

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

Comments: 49pages, 11figures

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

Characterization of Magnetic Flux Contents for Flux Transfer Events and its Implications for Flux Rope Formation at the Earth's Magnetopause

Authors: Shuo Wang, Ying Zou, Qiang Hu, Xueling Shi, Hiroshi Hasegawa

Abstract: Flux transfer events (FTEs) are a type of magnetospheric phenomena that exhibit distinctive observational signatures from the in-situ spacecraft measurements across the Earth's magnetopause. They are generally believed to possess a magnetic field configuration of a magnetic flux rope and formed through magnetic reconnection at the dayside magnetopause, sometimes accompanied with enhanced plasma co… ▽ More Flux transfer events (FTEs) are a type of magnetospheric phenomena that exhibit distinctive observational signatures from the in-situ spacecraft measurements across the Earth's magnetopause. They are generally believed to possess a magnetic field configuration of a magnetic flux rope and formed through magnetic reconnection at the dayside magnetopause, sometimes accompanied with enhanced plasma convection in the ionosphere. We examine two FTE events under the condition of southward interplanetary magnetic field (IMF) with a dawn-dusk component at the magnetopause by applying the Grad-Shafranov (GS) reconstruction method to the in-situ measurements by the Magnetospheric Multiscale (MMS) spacecraft to derive the magnetic flux contents associated with the FTE flux ropes. In particular, given a cylindrical magnetic flux rope configuration derived from the GS reconstruction, the magnetic flux content can be characterized by both the toroidal (axial) and poloidal fluxes. We then estimate the amount of magnetic flux (i.e., the reconnection flux) encompassed by the area ``opened" in the ionosphere, based on the ground-based Super Dual Auroral Radar Network (SuperDARN) observations. We find that for event 1, the FTE flux rope is oriented in the approximate dawn-dusk direction, and the amount of its poloidal magnetic flux agrees with the corresponding reconnection flux. For event 2, the agreement among the estimates of the magnetic fluxes is uncertain. We provide a detailed description about our interpretation for the topological features of the FTE flux ropes, based on a formation scenario of sequential magnetic field reconnection between adjacent field lines, consistent with our results. △ Less

Submitted 16 September, 2023; originally announced September 2023.

Comments: Submitted to JGR-Space Physics, and also to ESS Open Archive

Rapid pre-merger localization of binary neutron stars in third generation gravitational wave detectors

Authors: Qian Hu, John Veitch

Abstract: Pre-merger localization of binary neutron stars (BNSs) is one of the most important scientific goals for the third generation (3G) gravitational wave (GW) detectors. It will enable the electromagnetic observation of the whole process of BNS coalescence, especially for the pre-merger and merger phases which have not been observed yet, opening a window for deeper understandings of compact objects. T… ▽ More Pre-merger localization of binary neutron stars (BNSs) is one of the most important scientific goals for the third generation (3G) gravitational wave (GW) detectors. It will enable the electromagnetic observation of the whole process of BNS coalescence, especially for the pre-merger and merger phases which have not been observed yet, opening a window for deeper understandings of compact objects. To reach this goal, we describe a novel combination of multi-band matched filtering and semi-analytical localization algorithms to achieve early-warning localization of long BNS signals in 3G detectors. Using our method we are able to efficiently simulate one month of observations with a three-detector 3G network, and show that it is possible to provide accurate sky localizations more than 30 minutes before the merger. Our simulation shows that there could be ~ 10 (~ 100) BNS events localized within 100 deg2, 20 (6) minutes before merger, per month of observation. △ Less

Submitted 21 November, 2023; v1 submitted 2 September, 2023; originally announced September 2023.

Comments: 8 pages, 5 figures. Accepted version. Data release available

Report number: ET-0289A-23

Journal ref: Astrophys.J.Lett. 958 (2023) 2, L43

Magnetic Reconnection as the Key Mechanism in Sunspot Rotation Leading to Solar Eruption

Authors: Chaowei Jiang, Xueshang Feng, Xinkai Bian, Peng Zou, Aiying Duan, Xiaoli Yan, Qiang Hu, Wen He, Xinyi Wang, Pingbing Zuo, Yi Wang

Abstract: The rotation of sunspots around their umbral center has long been considered as an important process in leading to solar eruptions, but the underlying mechanism remains unclear. A prevailing physical picture on how sunspot rotation leads to eruption is that, by twisting the coronal magnetic field lines from their footpoints, the rotation can build up a magnetic flux rope and drive it into some kin… ▽ More The rotation of sunspots around their umbral center has long been considered as an important process in leading to solar eruptions, but the underlying mechanism remains unclear. A prevailing physical picture on how sunspot rotation leads to eruption is that, by twisting the coronal magnetic field lines from their footpoints, the rotation can build up a magnetic flux rope and drive it into some kinds of ideal magnetohydrodynamics (MHD) instabilities which initiate eruptions. Here with a data-inspired MHD simulation we studied the rotation of a large sunspot in solar active region NOAA 12158 leading to a major eruption, and found that it is distinct from prevailing theories based on ideal instabilities of twisted flux rope. The simulation suggests that, through successive rotation of the sunspot, the coronal magnetic field is sheared with a central current sheet created progressively within the sheared arcade before the eruption, but without forming a flux rope. Then the eruption is instantly triggered once fast reconnection sets in at the current sheet, while a highly twisted flux rope is created during the eruption. Furthermore, the simulation reveals an intermediate evolution stage between the quasi-static energy-storage phase and the impulsive eruption-acceleration phase. This stage may correspond to the slow-rise phase in observation and it enhances building up of the current sheet. △ Less

Submitted 30 September, 2023; v1 submitted 19 August, 2023; originally announced August 2023.

Comments: Updated from the initial version and text overlap with arXiv:2308.06977 is removed

Data-driven MHD simulation of a sunspot rotating active region leading to solar eruption

Authors: Chaowei Jiang, Xueshang Feng, Xinkai Bian, Peng Zou, Aiying Duan, Xiaoli Yan, Qiang Hu, Wen He, Xinyi Wang, Pingbing Zuo, Yi Wang

Abstract: Solar eruptions are the leading driver of space weather, and it is vital for space weather forecast to understand in what conditions the solar eruptions can be produced and how they are initiated. The rotation of sunspots around their umbral center has long been considered as an important condition in causing solar eruptions. To unveil the underlying mechanisms, here we carried out a data-driven m… ▽ More Solar eruptions are the leading driver of space weather, and it is vital for space weather forecast to understand in what conditions the solar eruptions can be produced and how they are initiated. The rotation of sunspots around their umbral center has long been considered as an important condition in causing solar eruptions. To unveil the underlying mechanisms, here we carried out a data-driven magnetohydrodynamics simulation for the event of a large sunspot with rotation for days in solar active region NOAA 12158 leading to a major eruption. The photospheric velocity as recovered from the time sequence of vector magnetograms are inputted directly at the bottom boundary of the numerical model as the driving flow. Our simulation successfully follows the long-term quasi-static evolution of the active region until the fast eruption, with magnetic field structure consistent with the observed coronal emission and onset time of simulated eruption matches rather well with the observations. Analysis of the process suggests that through the successive rotation of the sunspot the coronal magnetic field is sheared with a vertical current sheet created progressively, and once fast reconnection sets in at the current sheet, the eruption is instantly triggered, with a highly twisted flux rope originating from the eruption. This data-driven simulation stresses magnetic reconnection as the key mechanism in sunspot rotation leading to eruption. △ Less

Submitted 14 August, 2023; originally announced August 2023.

Comments: Accept by A&A

Search for Eccentric Black Hole Coalescences during the Third Observing Run of LIGO and Virgo

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, R. Abbott, H. Abe, 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, R. A. Alfaidi , et al. (1750 additional authors not shown)

Abstract: Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effect… ▽ More Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effects of eccentricity. Here, we present observational results for a waveform-independent search sensitive to eccentric black hole coalescences, covering the third observing run (O3) of the LIGO and Virgo detectors. We identified no new high-significance candidates beyond those that were already identified with searches focusing on quasi-circular binaries. We determine the sensitivity of our search to high-mass (total mass $M>70$ $M_\odot$) binaries covering eccentricities up to 0.3 at 15 Hz orbital frequency, and use this to compare model predictions to search results. Assuming all detections are indeed quasi-circular, for our fiducial population model, we place an upper limit for the merger rate density of high-mass binaries with eccentricities $0 < e \leq 0.3$ at $0.33$ Gpc$^{-3}$ yr$^{-1}$ at 90\% confidence level. △ Less

Submitted 7 August, 2023; originally announced August 2023.

Comments: 24 pages, 5 figures

Report number: LIGO-P2300080

Sequential Small Coronal Mass Ejections Observed In~situ and in White-Light Images by Parker Solar Probe

Authors: Brian E. Wood, Phillip Hess, Yu Chen, Qiang Hu

Abstract: We reconstruct the morphology and kinematics of a series of small transients that erupt from the Sun on 2021 April 24 using observations primarily from Parker Solar Probe (PSP). These sequential small coronal mass ejections (CMEs) may be the product of continuous reconnection at a current sheet, a macroscopic example of the more microscopic reconnection activity that has been proposed to accelerat… ▽ More We reconstruct the morphology and kinematics of a series of small transients that erupt from the Sun on 2021 April 24 using observations primarily from Parker Solar Probe (PSP). These sequential small coronal mass ejections (CMEs) may be the product of continuous reconnection at a current sheet, a macroscopic example of the more microscopic reconnection activity that has been proposed to accelerate the solar wind more generally. These particular CMEs are of interest because they are the first CMEs to hit PSP and be simultaneously imaged by it, using the Wide-field Imager for Solar Probe (WISPR) instrument. Based on imaging from WISPR and STEREO-A, we identify and model six discrete transients, and determine that it is the second of them (CME2) that first hits PSP, although PSP later more obliquely encounters the third transient as well. Signatures of these encounters are seen in the PSP in situ data. Within these data, we identify six candidate magnetic flux ropes (MFRs), all but one of which are associated with the second transient. The five CME2 MFRs have orientations roughly consistent with PSP encountering the right sides of roughly E-W oriented MFRs, which are sloping back towards the Sun. △ Less

Submitted 2 August, 2023; originally announced August 2023.

Comments: 17 pages, 11 figures, to appear in The Astrophysical Journal

Formation of an observed eruptive flux rope above the torus instability threshold through tether-cutting magnetic reconnection

Authors: Avijeet Prasad, Sanjay Kumar, Alphonse C. Sterling, Ronald L. Moore, Guillaume Aulanier, R. Bhattacharyya, Qiang Hu

Abstract: Erupting magnetic flux ropes (MFRs) play a crucial role in producing solar flares. However, the formation of erupting MFRs in complex coronal magnetic configurations and their subsequent evolution in the flaring events are not fully understood. We performed an MHD simulation of active region NOAA 12241 to understand the formation of a rising MFR during the onset of an M6.9 flare on 2014 December 1… ▽ More Erupting magnetic flux ropes (MFRs) play a crucial role in producing solar flares. However, the formation of erupting MFRs in complex coronal magnetic configurations and their subsequent evolution in the flaring events are not fully understood. We performed an MHD simulation of active region NOAA 12241 to understand the formation of a rising MFR during the onset of an M6.9 flare on 2014 December 18, around 21:41 UT. The MHD simulation was initialised with an extrapolated non-force-free magnetic field generated from the photospheric vector magnetogram of the active region taken a few minutes before the flare. The initial magnetic field topology displays a pre-existing sheared arcade enveloping the polarity inversion line. The simulated dynamics exhibit the movement of the oppositely directed legs of the sheared arcade field lines towards each other due to the converging Lorentz force, resulting in the onset of tether-cutting magnetic reconnection that produces an underlying flare arcade and flare ribbons. Concurrently, an MFR above the flare arcade develops inside the sheared arcade and shows a rising motion. The MFR is found to be formed in a torus-unstable region, thereby explaining its eruptive nature. Interestingly, the location and rise of the rope are in good agreement with the corresponding observations seen in EUV channels. Furthermore, the foot points of the simulation's flare arcade match well with the location of the observed parallel ribbons of the flare. The presented simulation supports the development of the MFR by the tether-cutting magnetic reconnection inside the sheared coronal arcade during flare onset. The MFR is then found to extend along the polarity inversion line (PIL) through slip-running reconnection. The MFR's eruptive nature is ascribed both to its formation in the torus-unstable region and also to the runaway tether-cutting reconnection. △ Less

Submitted 13 July, 2023; originally announced July 2023.

Comments: 12 pages, 8 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 677, A43 (2023)

Coronal Magnetic Field Extrapolation and Topological Analysis of Fine-Scale Structures during Solar Flare Precursors

Authors: Wen He, Qiang Hu, Ju Jing, Haimin Wang, Chaowei Jiang, Sushree S. Nayak, Avijeet Prasad

Abstract: Magnetic field plays an important role in various solar eruptions like flares, coronal mass ejections, etc. The formation and evolution of characteristic magnetic field topology in solar eruptions are critical problems that will ultimately help us understand the origination of these eruptions in the solar source regions. With the development of advanced techniques and instruments, observations wit… ▽ More Magnetic field plays an important role in various solar eruptions like flares, coronal mass ejections, etc. The formation and evolution of characteristic magnetic field topology in solar eruptions are critical problems that will ultimately help us understand the origination of these eruptions in the solar source regions. With the development of advanced techniques and instruments, observations with higher resolutions in different wavelengths and fields of view have provided more quantitative information for finer structures. So it is essential to improve our method to study the magnetic field topology in the solar source regions by taking advantage of high-resolution observations. In this study, we employ a nonlinear force-free field (NLFFF) extrapolation method based on a nonuniform grid setting for an M-class flare eruption event (SOL2015-06-22T17:39) with embedded magnetograms from the Solar Dynamics Observatory (SDO) and the Goode Solar Telescope (GST). The extrapolation results employing the embedded magnetogram for the bottom boundary are obtained by maintaining the native resolutions of the corresponding GST and SDO magnetograms. We compare the field line connectivity with the simultaneous GST/H$α$ and SDO/AIA observations for fine-scale structures associated with precursor brightenings. Then we perform a topological analysis of the field line connectivity corresponding to fine-scale magnetic field structures based on the extrapolation results. The results indicate that by combining the high-resolution GST magnetogram with a larger HMI magnetogram, the derived magnetic field topology is consistent with a scenario of magnetic reconnection among sheared field lines across the main polarity inversion line during solar flare precursors. △ Less

Submitted 5 June, 2023; originally announced June 2023.

Comments: 17 pages, 9 figures, 2 tables, submitted to ApJ

The First LHAASO Catalog of Gamma-Ray Sources

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

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

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

Comments: 40 pages, 13 figures, 4 tables

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

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

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

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

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

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

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

Two-injection scenario for the hard X-ray excess observed in Mrk 421

Authors: Wen Hu, Da-hai Yan, Qiang-lin Hu

Abstract: An interesting result recently reported for Mrk 421 is the detection of a significant excess at hard X-ray energies, which could provide useful information for investigating particle acceleration and emission mechanisms in the relativistic jet. Considering a two-injection scenario, we develop a self-consistent one-zone leptonic model to understand the origin of the hard X-ray excess in Mrk 421 dur… ▽ More An interesting result recently reported for Mrk 421 is the detection of a significant excess at hard X-ray energies, which could provide useful information for investigating particle acceleration and emission mechanisms in the relativistic jet. Considering a two-injection scenario, we develop a self-consistent one-zone leptonic model to understand the origin of the hard X-ray excess in Mrk 421 during the period of extremely low X-ray and very high energy (VHE) flux in 2013 January. In the model, two populations of mono-energetic ultrarelativistic electrons are injected into the emission region that is a magnetized plasmoid propagating along the blazar jet. We numerically calculate the emitting electron energy distribution by solving a kinetic equation that incorporates both shock acceleration and stochastic acceleration processes. Moreover, we infer analytic expressions relating electrons acceleration, cooling, escape and injection to the observed spectra and variability. For the injection luminosity in particular, we derive a new approximate analytical expression for the case of continual injection with a mono-energetic distribution. Based on a comparison between the theoretical predictions and the observed SED, we conclude that the hard X-ray excess observed in Mrk 421 may be due to the synchrotron radiation emitted by an additional electrons population, which is co-spatial with an electron population producing simultaneous Optical/UV, soft X-ray, and $γ$-ray emissions. The stochastic acceleration may play a major role in producing the observed X-ray spectrum. △ Less

Submitted 18 April, 2023; originally announced April 2023.

Comments: 12 pages, 4 figures, 2 tables, Accepted for publication in ApJ

Cluster aggregates surrounding Pismis 5 in the Vela Molecular Ridge

Authors: Ming Feng Qin, Yu Zhang, Jinzhong Liu, Fangfang Song, Qingshun Hu, Haozhi Wang, Shuo Ma, Guoliang Lv

Abstract: Context. In the Gaia era, the precision of astrometric data is unprecedented. High-quality data make it easier to find more cluster aggregates and support further confirmation of these open clusters. Aims. We use Gaia DR3 to redetermine the open clusters surrounding Pismis 5 in the Vela Molecular Ridge. We also investigate the basic properties of these clusters. Methods. We apply two clustering al… ▽ More Context. In the Gaia era, the precision of astrometric data is unprecedented. High-quality data make it easier to find more cluster aggregates and support further confirmation of these open clusters. Aims. We use Gaia DR3 to redetermine the open clusters surrounding Pismis 5 in the Vela Molecular Ridge. We also investigate the basic properties of these clusters. Methods. We apply two clustering algorithms (StarGO and pyUPMASK) to identify the open cluster members in a five-dimensional space with Gaia DR3. Results. We identify eight open clusters surrounding Pismis 5 in the Vela Molecular Ridge. The open cluster QZ 1 is newly discovered. Through investigating the comprehensive properties of the clusters, one open binary cluster candidate (Alessi 43 and Collinder 197) and one triple open cluster candidate (Pismis 5, Pismis 5A, and Pismis 5B) are discussed. Conclusions. Binary and triple open cluster candidates have been identified as potential primordial aggregates based on their similar age, position, and motion. According to kinematic speculations, the two aggregate candidates will gradually separate, and their interiors will slowly disintegrate. △ Less

Submitted 26 August, 2023; v1 submitted 1 April, 2023; originally announced April 2023.

Comments: 10 pages, 7 figures

Journal ref: A&A 675, A67 (2023)

A detection of the layered structure of nearby open clusters

Authors: Qingshun Hu, Yu Zhang, Ali Esamdin, Hong Wang, Mingfeng Qin

Abstract: We applied the newly developed rose diagram overlay method to detect the layered structure of 88 nearby open clusters ($\leq$500~pc) on the three projections after the distance correction of their member stars, based on the catalog in literature. The results show that with the rose diagram overlay method, a total of 74 clusters in our sample have a layered structure, while the remaining clusters a… ▽ More We applied the newly developed rose diagram overlay method to detect the layered structure of 88 nearby open clusters ($\leq$500~pc) on the three projections after the distance correction of their member stars, based on the catalog in literature. The results show that with the rose diagram overlay method, a total of 74 clusters in our sample have a layered structure, while the remaining clusters are without a clear layered structure. We for the first time defined the layered structure parameters for the sample clusters. Meanwhile, we found that the layered circle core area ($s$) has a strong positive correlation with the number of cluster members, while the kernel instability index ($η$) has a strong negative correlation with the number of cluster members. Our study provides a novel perspective for the detection of the layered structure of open clusters. △ Less

Submitted 17 May, 2023; v1 submitted 6 February, 2023; originally announced February 2023.

Comments: 11 pages, 9 figures, 6 tables, accepted for publication in A&A

Journal ref: A&A 672, A12 (2023)

The Structure and Origin of Switchbacks: Parker Solar Probe Observations

Authors: Jia Huang, J. C. Kasper, L. A. Fisk, Davin E. Larson, Michael D. McManus, C. H. K. Chen, Mihailo M. Martinović, K. G. Klein, Luke Thomas, Mingzhe Liu, Bennett A. Maruca, Lingling Zhao, Yu Chen, Qiang Hu, Lan K. Jian, J. L. Verniero, Marco Velli, Roberto Livi, P. Whittlesey, Ali Rahmati, Orlando Romeo, Tatiana Niembro, Kristoff Paulson, M. Stevens, A. W. Case , et al. (3 additional authors not shown)

Abstract: Switchbacks are rapid magnetic field reversals that last from seconds to hours. Current Parker Solar Probe (PSP) observations pose many open questions in regard to the nature of switchbacks. For example, are they stable as they propagate through the inner heliosphere, and how are they formed? In this work, we aim to investigate the structure and origin of switchbacks. In order to study the stabili… ▽ More Switchbacks are rapid magnetic field reversals that last from seconds to hours. Current Parker Solar Probe (PSP) observations pose many open questions in regard to the nature of switchbacks. For example, are they stable as they propagate through the inner heliosphere, and how are they formed? In this work, we aim to investigate the structure and origin of switchbacks. In order to study the stability of switchbacks, we suppose the small-scale current sheets therein are generated by magnetic braiding, and they should work to stabilize the switchbacks. With more than one thousand switchbacks identified with PSP observations in seven encounters, we find many more current sheets inside than outside switchbacks, indicating that these microstructures should work to stabilize the S-shaped structures of switchbacks. Additionally, we study the helium variations to trace the switchbacks to their origins. We find both helium-rich and helium-poor populations in switchbacks, implying that the switchbacks could originate from both closed and open magnetic field regions in the Sun. Moreover, we observe that the alpha-proton differential speeds also show complex variations as compared to the local Alfvén speed. The joint distributions of both parameters show that low helium abundance together with low differential speed is the dominant state in switchbacks. The presence of small-scale current sheets in switchbacks along with the helium features are in line with the hypothesis that switchbacks could originate from the Sun via interchange reconnection process. However, other formation mechanisms are not excluded. △ Less

Submitted 22 May, 2023; v1 submitted 24 January, 2023; originally announced January 2023.

Comments: accepted by ApJ

Investigating pre-eruptive magnetic properties at the footprints of erupting magnetic flux ropes

Authors: Wensi Wang, Jiong Qiu, Rui Liu, Chunming Zhu, Kai E Yang, Qiang Hu, Yuming Wang

Abstract: It is well established that solar eruptions are powered by free magnetic energy stored in current-carrying magnetic field in the corona. It has also been generally accepted that magnetic flux ropes (MFRs) are a critical component of many coronal mass ejections (CMEs). What remains controversial is whether MFRs are present well before the eruption. Our aim is to identify progenitors of MFRs, and in… ▽ More It is well established that solar eruptions are powered by free magnetic energy stored in current-carrying magnetic field in the corona. It has also been generally accepted that magnetic flux ropes (MFRs) are a critical component of many coronal mass ejections (CMEs). What remains controversial is whether MFRs are present well before the eruption. Our aim is to identify progenitors of MFRs, and investigate pre-eruptive magnetic properties associated with these progenitors. Here we analyze 28 MFRs erupting within 45 deg from the disk center from 2010 to 2015. All MFRs'feet are well identified by conjugate coronal dimmings. We then calculate magnetic properties at the feet of the MFRs, prior to their eruptions, using Helioseismic and Magnetic Imager (HMI) vector magnetograms. Our results show that only 8 erupting MFRs are associated with significant non-neutralized electric currents, 4 of which also exhibit pre-eruptive dimmings at the foot-prints. Twist and current distributions are asymmetric at the two feet of these MFRs. The presence of pre-eruption dimmings associated with non-neutralized currents suggests the pre-existing MFRs. Furthermore, evolution of conjugate dimmings and electric currents within the foot-prints can provide clues about the internal structure of MFRs and their formation mechanism. △ Less

Submitted 28 November, 2022; originally announced November 2022.

Comments: 40 pages, 13 figures, 6 tables, Accepted for publication in ApJ

Searching for pulsars with phase characteristics

Authors: Bo Peng, Qian-Chen Hu, Qiang Li, Lei Qian, Xiao-Bo Dong, Shi-Lin Peng, Ze-Lin Wang

Abstract: We present a method by using the phase characteristics of radio observation data for pulsar search and candidate identification. The phase characteristics are relations between the pulsar signal and the phase correction in the frequency-domain, and we regard it as a new search diagnostic characteristic. Based on the phase characteristics, a search method is presented: calculating DM (dispersion me… ▽ More We present a method by using the phase characteristics of radio observation data for pulsar search and candidate identification. The phase characteristics are relations between the pulsar signal and the phase correction in the frequency-domain, and we regard it as a new search diagnostic characteristic. Based on the phase characteristics, a search method is presented: calculating DM (dispersion measure) -- frequency data to select candidate frequencies, and then confirming of candidates by using the broadband characteristics of pulsar signals. Based on this method, we performed a search test on short observation data of M15 and M71, which were observed by Five-hundred-meter Aperture spherical radio Telescope (FAST), and some of the Galactic Plane Pulsar Snapshot survey (GPPS) data. Results show that it can get similar search results to PRESTO (PulsaR Exploration and Search TOolkit) while having a faster processing speed. △ Less

Submitted 16 November, 2022; originally announced November 2022.

Accumulating errors in tests of general relativity with gravitational waves: overlapping signals and inaccurate waveforms

Authors: Qian Hu, John Veitch

Abstract: Observations of gravitational waves (GWs) from compact binary coalescences provide powerful tests of general relativity (GR), but systematic errors in data analysis could lead to incorrect scientific conclusions. This issue is especially serious in the third-generation GW detectors in which the signal-to-noise ratio (SNR) is high and the number of detections is large. In this work, we investigate… ▽ More Observations of gravitational waves (GWs) from compact binary coalescences provide powerful tests of general relativity (GR), but systematic errors in data analysis could lead to incorrect scientific conclusions. This issue is especially serious in the third-generation GW detectors in which the signal-to-noise ratio (SNR) is high and the number of detections is large. In this work, we investigate the impacts of overlapping signals and inaccurate waveform models on tests of GR. We simulate mock catalogs for Einstein Telescope and Cosmic Explorer and perform parametric tests of GR using waveform models with different levels of inaccuracy. We find the systematic error in non-GR parameter estimates could accumulate toward a false deviation from GR when combining results from multiple events, although a bayesian model selection analysis may not favour a deviation. Waveform inaccuracies contribute most to the systematic errors, but multiple overlapping signals could magnify the effects of systematics due to the incorrect removal of signals. We also point out that testing GR using selected ''golden binaries'' with high SNR is even more vulnerable to false deviations from GR. The problem of error accumulation is universal; we emphasize that it must be addressed to fully exploit the data from third-generation GW detectors, and that further investigations, particularly in waveform accuracy, will be essential. △ Less

Submitted 13 February, 2023; v1 submitted 10 October, 2022; originally announced October 2022.

Comments: 13 pages, 6 figures

Report number: ET-0211A-22

Journal ref: The Astrophysical Journal, 945, 103 (2023)

Binary open clusters in the Gaia data

Authors: Fangfang Song, Ali Esamdin, Qingshun Hu, Mengfan Zhang

Abstract: Context. Observations indicate that the fraction of potential binary star clusters in the Milky Way is either the same or lower than that of the Magellanic Clouds. The unprecedented precision in the parallax measurements by Gaia has allowed for the discovery of a growing number of new binary open clusters (OCs). Aims. We aim to survey the candidates of truly binary open clusters that are formed si… ▽ More Context. Observations indicate that the fraction of potential binary star clusters in the Milky Way is either the same or lower than that of the Magellanic Clouds. The unprecedented precision in the parallax measurements by Gaia has allowed for the discovery of a growing number of new binary open clusters (OCs). Aims. We aim to survey the candidates of truly binary open clusters that are formed simultaneously, using information from the Gaia database. Methods. Based on the most recent catalog of open clusters, we investigated the interactions of adjacent binary open clusters in our Galaxy within separations of 50 pc. We compared their coordinates, proper motions, parallaxes, and color-magnitude diagrams(CMDs) via binary plots for all candidate pairs. The candidates of truly binary open clusters are selected on the basis of their common proper motions and consistent behaviors in the CMDs of different clusters that are limited to a separation of 50 pc. Results. About ten pairs of the selected binary open clusters appear to be the same clusters, based on evidence that almost half ofthe cluster members are shared. Fourteen pairs are possibly true binaries, implying that they may come from the same clouds, among which five pairs are newly discovered. In addition, two clusters, UBC 46 and UBC 192, were found to be part of the stellar complex LISCA I. Our results confirm that OCs born in groups are usually composed of young open clusters. △ Less

Submitted 27 August, 2022; originally announced August 2022.

Comments: 9 pages, 3 figures

Journal ref: A&A 666, A75 (2022)

A Low-speed Intruder Star in Hyades: A Temporary Residence

Authors: Qingshun Hu, Yu Zhang, Ali Esamdin, Dengkai Jiang, Mingfeng Qin, Ning Chang, Haozhi Wang

Abstract: We hereby report a low-speed (about~21~km$\cdot$~s$^{-1}$ with respect to the Sun) intruder member in the Hyades cluster based on the data in the literature. The results show that the star is a non-native member star for the Hyades, with its radial velocity being smaller than the radial velocity of the Hyades cluster, even exceeding the standard deviation of the radial velocity of the cluster by a… ▽ More We hereby report a low-speed (about~21~km$\cdot$~s$^{-1}$ with respect to the Sun) intruder member in the Hyades cluster based on the data in the literature. The results show that the star is a non-native member star for the Hyades, with its radial velocity being smaller than the radial velocity of the Hyades cluster, even exceeding the standard deviation of the radial velocity of the cluster by a factor of 9. Furthermore, by analyzing and comparing the orbits of this star and its host, it may have intruded into its host in the past 2~Myr. If the star's current motion orbit remains unchanged, it may leave its host in the next 2~Myr. This implies that the intruder star may be temporarily residing in the cluster. This study presents the first observational evidence of a star intrusion into a cluster, which suggests that more evidence may be found. △ Less

Submitted 4 August, 2022; v1 submitted 26 July, 2022; originally announced July 2022.

Comments: 9 pages, 7 figures; accepted for publication in ApJ

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

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

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

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

Comments: 18 pages, 11 figures

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

A photometric study of two contact binaries: CRTS J025408.1+265957 and CRTS J012111.1+272933

Authors: S. Ma, J. Z. Liu, Y. Zhang, Q. S. Hu, G. L. Lv

Abstract: We performed new photometric observations for two contact binaries (i.e., CRTS J025408.1+265957 and CRTS J012111.1+272933), which were observed by the 1.0-m telescope at Xingjiang Astronomical Observatory. From our light curves and several survey data, we derived several sets of photometric solutions. We found that CRTS J025408.1+265957 and CRTS J012111.1+272933 were Aand W-type W UMa, respectivel… ▽ More We performed new photometric observations for two contact binaries (i.e., CRTS J025408.1+265957 and CRTS J012111.1+272933), which were observed by the 1.0-m telescope at Xingjiang Astronomical Observatory. From our light curves and several survey data, we derived several sets of photometric solutions. We found that CRTS J025408.1+265957 and CRTS J012111.1+272933 were Aand W-type W UMa, respectively. The results imply that the spot migrates or disappears in two contact binaries, which were identified by chromospheric activity emissions (e.g. Hα emission) from LAMOST spectra. From the O-C curves, the orbital periods of two contact binaries may be increasing, which is interpreted by the mass transfer from the less massive component to the more massive one. With mass transferring, two contact binaries may evolve from the contact configurations to semi-detected ones as predicted by the theory of thermal relaxation oscillation. △ Less

Submitted 13 July, 2022; originally announced July 2022.

Comments: 16 pages, 10 figures

Searching for Variable Stars in the Open Cluster NGC 2355 and Its Surrounding Region

Authors: Hong Wang, Yu Zhang, Xiangyun Zeng, Qingshun Hu, Jinzhong Liu, Mingfeng Qin, Guoliang Lü

Abstract: We have investigated the variable stars in the field surrounding NGC 2355 based on the time-series photometric observation data. More than 3000 CCD frames were obtained in the V band spread over 13 nights with the Nanshan One-meter Wide-field Telescope. We have detected 88 variable stars, containing 72 new variable stars and 16 known variable stars. By analyzing these light curves, we classified t… ▽ More We have investigated the variable stars in the field surrounding NGC 2355 based on the time-series photometric observation data. More than 3000 CCD frames were obtained in the V band spread over 13 nights with the Nanshan One-meter Wide-field Telescope. We have detected 88 variable stars, containing 72 new variable stars and 16 known variable stars. By analyzing these light curves, we classified the variable stars as follows: 26 eclipsing binaries, 52 pulsating stars, 4 rotating variables, and 6 unclear type variable stars for which their periods are much longer than the time baseline chosen. Employing Gaia DR2 parallax, kinematics, and photometry, the cluster membership of these variable stars were also analyzed for NGC 2355. In addition to the 11 variable members reported by Cantat-Gaudin et al. (2018), we identify 4 more variable member candidates located at the outer region of NGC 2355 and showed homogeneity in space positions and kinematic properties with the cluster members. The main physical parameters of NGC 2355 estimated from the two-color and color-magnitude diagrams are log(age/yr) = 8.9, E(B - V) = 0.24 mag, and [Fe/H] = - 0.07 dex. △ Less

Submitted 13 June, 2022; originally announced June 2022.

Comments: 15 pages, 11 figures. 6 tables,Accepted for publication in AJ

A Magnetic Flux Rope Configuration Derived by Optimization of Two-Spacecraft In-situ Measurements

Authors: Qiang Hu, Wen He, Yu Chen

Abstract: Increasingly one interplanetary coronal mass ejection (ICME) structure can propagate across more than one spacecraft in the solar wind. This usually happens when two or more spacecraft are nearly radially aligned with a relatively small longitudinal separation angle from one another. This provides multi-point measurements of the same structure and enables better characterization and validation of… ▽ More Increasingly one interplanetary coronal mass ejection (ICME) structure can propagate across more than one spacecraft in the solar wind. This usually happens when two or more spacecraft are nearly radially aligned with a relatively small longitudinal separation angle from one another. This provides multi-point measurements of the same structure and enables better characterization and validation of modeling results of the structures embedded in these ICMEs. We report such an event during October 13-14, 2019 when the Solar TErrestrial RElations Observatory Ahead (STA) spacecraft and the Parker Solar Probe (PSP) crossed one ICME structure at two different locations with nominal separations in both heliocentric distances and the longitudinal angles. We first perform an optimal fitting to the STA in-situ measurements, based on an analytic quasi-three dimensional (3D) model, yielding a minimum reduced $χ^2=0.468$. Then we further apply the optimization approach by combining the magnetic field measurements from both spacecraft along their separate paths across the ICME structure. We find that the output based on the optimization (with the minimum reduced $χ^2=3.15$) of the combined two-spacecraft dataset yields a more consistent result, given the much improved agreement of the model output with PSP data. The result demonstrates a magnetic flux rope configuration with clear 3D spatial variations. △ Less

Submitted 2 June, 2022; originally announced June 2022.

Comments: Submitted to Frontiers in Physics (Stellar and Solar Physics)

Validation and interpretation of three-dimensional configuration of a magnetic cloud flux rope

Authors: Qiang Hu, Chunming Zhu, Wen He, Jiong Qiu, Lan K. Jian, Avijeet Prasad

Abstract: One "strong" magnetic cloud (MC) with the magnetic field magnitude reaching $\sim$ 40 nT at 1 au during 2012 June 16-17 is examined in association with a pre-existing magnetic flux rope (MFR) identified on the Sun. The MC is characterized by a quasi-three dimensional (3D) flux rope model based on in situ measurements from the Wind spacecraft. The magnetic flux contents and other parameters are qua… ▽ More One "strong" magnetic cloud (MC) with the magnetic field magnitude reaching $\sim$ 40 nT at 1 au during 2012 June 16-17 is examined in association with a pre-existing magnetic flux rope (MFR) identified on the Sun. The MC is characterized by a quasi-three dimensional (3D) flux rope model based on in situ measurements from the Wind spacecraft. The magnetic flux contents and other parameters are quantified. In addition, a correlative study with the corresponding measurements of the same structure crossed by the Venus Express (VEX) spacecraft at a heliocentric distance 0.7 au and with an angular separation $\sim 6^\circ$ in longitude is performed to validate the MC modeling results. The spatial variation between the Wind and VEX magnetic field measurements is attributed to the 3D configuration of the structure as featured by a knotted bundle of flux. The comparison of the magnetic flux contents between the MC and the source region on the Sun indicates that the 3D reconnection process accompanying an M1.9 flare may correspond to the magnetic reconnection between the field lines of the pre-existing MFR rooted in the opposite polarity footpoints. Such a process reduces the amount of the axial magnetic flux in the erupted flux rope, by approximately 50\%, in this case. △ Less

Submitted 7 April, 2022; originally announced April 2022.

Comments: Submitted to ApJ

The magnetic topology of the inverse Evershed flow

Authors: A. Prasad, M. Ranganathan, C. Beck, D. P. Choudhary, Q. Hu

Abstract: The inverse Evershed flow (IEF) is a mass motion towards sunspots at chromospheric heights. We combined high-resolution observations of NOAA 12418 from the Dunn Solar Telescope and vector magnetic field measurements from the Helioseismic and Magnetic Imager (HMI) to determine the driver of the IEF. We derived chromospheric line-of-sight (LOS) velocities from spectra of H$α$ and Ca II IR. The HMI d… ▽ More The inverse Evershed flow (IEF) is a mass motion towards sunspots at chromospheric heights. We combined high-resolution observations of NOAA 12418 from the Dunn Solar Telescope and vector magnetic field measurements from the Helioseismic and Magnetic Imager (HMI) to determine the driver of the IEF. We derived chromospheric line-of-sight (LOS) velocities from spectra of H$α$ and Ca II IR. The HMI data were used in a non-force-free magnetic field extrapolation to track closed field lines near the sunspot in the active region. We determined their length and height, located their inner and outer foot points, and derived flow velocities along them. The magnetic field lines related to the IEF reach on average a height of 3 Mm over a length of 13 Mm. The inner (outer) foot points are located at 1.2 (1.9) sunspot radii. The average field strength difference $ΔB$ between inner and outer foot points is +400 G. The temperature difference $ΔT$ is anti-correlated with $ΔB$ with an average value of -100 K. The pressure difference $Δp$ is dominated by $ΔB$ and is primarily positive with a driving force towards the inner foot points of 1.7 kPa on average. The velocities predicted from $Δp$ reproduce the LOS velocities of 2-10 km s$^{-1}$ with a square-root dependence. We find that the IEF is driven along magnetic field lines connecting network elements with the outer penumbra by a gas pressure difference that results from a difference in field strength as predicted by the classical siphon flow scenario. △ Less

Submitted 5 March, 2022; originally announced March 2022.

Comments: 23 pages, 27 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 662, A25 (2022)

Quantitative Characterization of Magnetic Flux Rope Properties for Two Solar Eruption Events

Authors: Wen He, Qiang Hu, Chaowei Jiang, Jiong Qiu, Avijeet Prasad

Abstract: In order to bridge the gap between heliospheric and solar observations of coronal mass ejections (CMEs), one of the key steps is to improve the understanding of their corresponding magnetic structures like the magnetic flux ropes (MFRs). But it remains a challenge to confirm the existence of a coherent MFR before or upon the CME eruption on the Sun and to quantitatively characterize the CME-MFR du… ▽ More In order to bridge the gap between heliospheric and solar observations of coronal mass ejections (CMEs), one of the key steps is to improve the understanding of their corresponding magnetic structures like the magnetic flux ropes (MFRs). But it remains a challenge to confirm the existence of a coherent MFR before or upon the CME eruption on the Sun and to quantitatively characterize the CME-MFR due to the lack of direct magnetic field measurement in the corona. In this study, we investigate the MFR structures, originating from two active regions (ARs), AR 11719 and AR 12158, and estimate their magnetic properties quantitatively. We perform the nonlinear force-free field extrapolations with preprocessed photospheric vector magnetograms. In addition, remote-sensing observations are employed to find indirect evidence of MFRs on the Sun and to analyze the time evolution of magnetic reconnection flux associated with the flare ribbons during the eruption. A coherent "pre-existing" MFR structure prior to the flare eruption is identified quantitatively for one event from the combined analysis of the extrapolation and observation. Then the characteristics of MFRs for two events on the Sun before and during the eruption, forming the CME-MFR, including the axial magnetic flux, field-line twist, and reconnection flux, are estimated and compared with the corresponding in situ modeling results. We find that the magnetic reconnection associated with the accompanying flares for both events injects significant amount of flux into the erupted CME-MFRs. △ Less

Submitted 9 January, 2022; originally announced January 2022.

Comments: 28 pages, 10 figures, 3 tables, submitted to ApJ

Comparison of the Hall Magnetohydrodynamics and Magnetohydrodynamics evolution of a flaring solar active region

Authors: Kamlesh Bora, Ramit Bhattacharyya, Avijeet Prasad, Bhuwan Joshi, Qiang Hu

Abstract: This work analyzes the Hall magnetohydrodynamics (HMHD) and magnetohydrodynamics (MHD) numerical simulations of a flaring solar active region as a testbed while idealizing the coronal Alfvén speed to be of two orders of magnitude lesser. HMHD supports faster magnetic reconnection and shows richer complexity in magnetic field line evolution compared to the MHD. The magnetic reconnections triggering… ▽ More This work analyzes the Hall magnetohydrodynamics (HMHD) and magnetohydrodynamics (MHD) numerical simulations of a flaring solar active region as a testbed while idealizing the coronal Alfvén speed to be of two orders of magnitude lesser. HMHD supports faster magnetic reconnection and shows richer complexity in magnetic field line evolution compared to the MHD. The magnetic reconnections triggering the flare are explored by numerical simulations augmented with relevant multi-wavelength observations. The initial coronal magnetic field is constructed by non-force-free extrapolation of photospheric vector magnetic field. Magnetic structure involved in the flare is identified to be a flux rope, with its overlying magnetic field lines constituting the quasi-separatrix layers (QSLs) along with a three-dimensional null point and a null line. Compared to the MHD simulation, the HMHD simulation shows a higher and faster ascend of the rope together with the overlying field lines, which further reconnect at the QSL located higher up in the corona. The foot points of the field lines match better with the observations for the HMHD case with the central part of the flare ribbon located at the chromosphere. Additionally, field lines are found to rotate in a circular pattern in the HMHD, whereas no such rotation is seen in the MHD results. Interestingly, plasma is also observed to be rotating in a co-spatial chromospheric region, which makes the HMHD simulation more credible. Based on the aforementioned agreements, HMHD simulation is found to agree better with observations and, thus, opens up a novel avenue to explore. △ Less

Submitted 3 December, 2021; originally announced December 2021.

Comments: 31 pages, 17 figures (accepted for publication in The Astrophysical Journal)

Small-scale magnetic flux ropes and their properties based on in-situ measurements from Parker Solar Probe

Authors: Yu Chen, Qiang Hu

Abstract: We report small-scale magnetic flux ropes via the Parker Solar Probe in situ measurements during the first six encounters and present additional analyses to supplement our prior work in Chen et al. 2021. These flux ropes are detected by the Grad-Shafranov-based algorithm with the duration and scale size ranging from 10 seconds to $\lesssim$1 hour and from a few hundred kilometers to 10$^{-3}$ au,… ▽ More We report small-scale magnetic flux ropes via the Parker Solar Probe in situ measurements during the first six encounters and present additional analyses to supplement our prior work in Chen et al. 2021. These flux ropes are detected by the Grad-Shafranov-based algorithm with the duration and scale size ranging from 10 seconds to $\lesssim$1 hour and from a few hundred kilometers to 10$^{-3}$ au, respectively. They include both static structures and those with significant field-aligned plasma flows. Most structures tend to possess large cross helicity, while the residual energy distributes in wide ranges. We find that these dynamic flux ropes mostly propagate anti-sunward, with no preferential sign of magnetic helicity. The magnetic flux function follows a power law and is proportional to scale size. We also present case studies showing reconstructed two-dimensional (2D) configurations, which confirm that the static and dynamic flux ropes have the common configuration of spiral magnetic field lines (also streamlines). Moreover, the existence of such events hints at the interchange reconnection as a possible mechanism to generate flux rope-like structures near the Sun. Lastly, we summarize the major findings and discuss the possible correlation between these flux rope-like structures and turbulence due to the process of local Alfvenic alignment. △ Less

Submitted 17 November, 2021; originally announced November 2021.

Comments: Accepted by ApJ on Oct 26 2021

Semianalytical Approach for Sky Localization of Gravitational Waves

Authors: Qian Hu, Cong Zhou, Jhao-Hong Peng, Linqing Wen, Qi Chu, Manoj Kovalam

Abstract: Rapid sky localization of gravitational wave sources is crucial to enable prompt electromagnetic follow-ups. In this article, we present a novel semianalytical approach for sky localization of gravitational waves from compact binary coalescences. We use the Bayesian framework with an analytical approximation to the prior distributions for a given astrophysical model. We derive a semianalytical sol… ▽ More Rapid sky localization of gravitational wave sources is crucial to enable prompt electromagnetic follow-ups. In this article, we present a novel semianalytical approach for sky localization of gravitational waves from compact binary coalescences. We use the Bayesian framework with an analytical approximation to the prior distributions for a given astrophysical model. We derive a semianalytical solution to the posterior distribution of source directions. This method only requires one-fold numerical integral that marginalizes over the merger time, compared to the five-fold numerical integration otherwise needed in the Bayesian localization method. The performance of the method is demonstrated using a set of binary neutron stars (BNS) injections on Gaussian noise using LIGO-Virgo's design and O2 sensitivity. We find the median of 90% confidence area in O2 sensitivity to be $\mathcal{O}(10^2) ~\mathrm{deg}^2$, comparable to that of the existing LIGO-Virgo online localization method Bayestar and parameter estimation toolkit LALInference. In the end, we apply this method to localize the BNS event GW170817 and find the 50% (90%) confidence region of 11 $\mathrm{deg}^2$ (50 $\mathrm{deg}^2$). The detected optical counterpart of GW170817 resides within our 50% confidence area. △ Less

Submitted 3 November, 2021; v1 submitted 5 October, 2021; originally announced October 2021.

Comments: 15 pages, 8 figures

Journal ref: Phys. Rev. D 104, 104008 (2021)

The $^{59}$Fe(n, γ)$^{60}$Fe Cross Section from the Surrogate Ratio Method and Its Effect on the $^{60}$Fe Nucleosynthesis

Authors: S. Q. Yan, X. Y. Li, K. Nishio, M. Lugaro, Z. H. Li, H. Makii, M. Pignatari, Y. B. Wang, R. Orlandi, K. Hirose, K. Tsukada, P. Mohr, G. S. Li, J. G. Wang, B. S. Gao, Y. L. Han, B. Guo, Y. J. Li, Y. P. Shen, T. K. Sato, Y. Ito, F. Suzaki, J. Su, Y. Y. Yang, J. S. Wang , et al. (17 additional authors not shown)

Abstract: The long-lived $^{60}$Fe (with a half-life of 2.62 Myr) is a crucial diagnostic of active nucleosynthesis in the Milky Way galaxy and in supernovae near the solar system. The neutron-capture reaction $^{59}$Fe(n,$γ$)$^{60}$Fe on $^{59}$Fe (half-life = 44.5 days) is the key reaction for the production of $^{60}$Fe in massive stars. This reaction cross section has been previously constrained by the… ▽ More The long-lived $^{60}$Fe (with a half-life of 2.62 Myr) is a crucial diagnostic of active nucleosynthesis in the Milky Way galaxy and in supernovae near the solar system. The neutron-capture reaction $^{59}$Fe(n,$γ$)$^{60}$Fe on $^{59}$Fe (half-life = 44.5 days) is the key reaction for the production of $^{60}$Fe in massive stars. This reaction cross section has been previously constrained by the Coulomb dissociation experiment, which offered partial constraint on the $E$1 $γ$-ray strength function but a negligible constraint on the $M$1 and $E$2 components. In this work, for the first time, we use the surrogate ratio method to experimentally determine the $^{59}$Fe(n,$γ$)$^{60}$Fe cross sections in which all the components are included. We derived a Maxwellian-averaged cross section of 27.5 $\pm$ 3.5 mb at $kT$= 30 keV and 13.4 $\pm$ 1.7 mb at $kT$= 90 keV, roughly 10 - 20% higher than previous estimates. We analyzed the impact of our new reaction rates in nucleosynthesis models of massive stars and found that uncertainties in the production of $^{60}$Fe from the $^{59}$Fe(n,$γ$)$^{60}$Fe rate are at most of 25%. We conclude that stellar physics uncertainties now play a major role in the accurate evaluation of the stellar production of $^{60}$Fe. △ Less

Submitted 26 September, 2021; originally announced September 2021.

Comments: 9 pages with 6 figures

Decoding the Morphological Evolution of Open Clusters

Authors: Qingshun Hu, Yu Zhang, Ali Esamdin

Abstract: Base on Gaia Second Data Release and the combination of nonparametric bivariate density estimation with the least square ellipse fitting, we derive the shape parameters of the sample clusters. By analyzing the dislocation of the sample clusters, the dislocation $d$ is related to the X-axis pointing toward the Galactic center, Y-axis pointing in the direction of Galactic rotation, and the Z-axis (l… ▽ More Base on Gaia Second Data Release and the combination of nonparametric bivariate density estimation with the least square ellipse fitting, we derive the shape parameters of the sample clusters. By analyzing the dislocation of the sample clusters, the dislocation $d$ is related to the X-axis pointing toward the Galactic center, Y-axis pointing in the direction of Galactic rotation, and the Z-axis (log(|H|/pc)) that is positive toward the Galactic north pole. This finding underlines the important role of the dislocation of clusters in tracking the external environment of the Milky Way. The orientation ($q_{pm}$) of the clusters with $e_{pm}$ $\geq$ 0.4 presents an aggregate distribution in the range of -45$\degr$ to 45$\degr$, about 74\% of them. This probably suggests that these clusters tend to deform heavily in the direction of the Galactic plane. NGC 752 is in a slight stage of expansion in the two-dimensional space and will deform itself morphology along the direction perpendicular to the original stretching direction in the future if no other events occur. The relative degree of deformation of the sample clusters in the short-axis direction decreases as their ages increase. On average, the severely distorted sample clusters in each group account for about 26\% $\pm$ 9\%. This possibly implies a uniform external environment in the range of $|$H$|$ $\leq$ 300 pc if the sample completeness of each group is not taken into account. △ Less

Submitted 31 October, 2021; v1 submitted 10 September, 2021; originally announced September 2021.

Comments: 11pages,11 figures, accepted by Astronomy & Astrophysics

Journal ref: A&A 656, A49 (2021)

A fundamental mechanism of solar eruption initiation

Authors: Chaowei Jiang, Xueshang Feng, Rui Liu, XiaoLi Yan, Qiang Hu, Ronald L. Moore, Aiying Duan, Jun Cui, Pingbing Zuo, Yi Wang, Fengsi Wei

Abstract: Solar eruptions are spectacular magnetic explosions in the Sun's corona, and how they are initiated remains unclear. Prevailing theories often rely on special magnetic topologies that may not generally exist in the pre-eruption source region of corona. Here, using fully three-dimensional magnetohydrodynamic simulations with high accuracy, we show that solar eruptions can be initiated in a single b… ▽ More Solar eruptions are spectacular magnetic explosions in the Sun's corona, and how they are initiated remains unclear. Prevailing theories often rely on special magnetic topologies that may not generally exist in the pre-eruption source region of corona. Here, using fully three-dimensional magnetohydrodynamic simulations with high accuracy, we show that solar eruptions can be initiated in a single bipolar configuration with no additional special topology. Through photospheric shearing motion alone, an electric current sheet forms in the highly sheared core field of the magnetic arcade during its quasi-static evolution. Once magnetic reconnection sets in, the whole arcade is expelled impulsively, forming a fast-expanding twisted flux rope with a highly turbulent reconnecting region underneath. The simplicity and efficacy of this scenario argue strongly for its fundamental importance in the initiation of solar eruptions. △ Less

Submitted 17 July, 2021; originally announced July 2021.

Comments: Published in Nature Astronomy (2021)

Configuration of a Magnetic Cloud from Solar Orbiter and Wind Spacecraft In-situ Measurements

Authors: Qiang Hu, Wen He, Lingling Zhao, Edward Lu

Abstract: Coronal mass ejections (CMEs) represent one type of the major eruption from the Sun. Their interplanetary counterparts, the interplanetary CMEs (ICMEs), are the direct manifestations of these structures when they propagate into the heliosphere and encounter one or more observing spacecraft. The ICMEs generally exhibit a set of distinctive signatures from the in-situ spacecraft measurements. A part… ▽ More Coronal mass ejections (CMEs) represent one type of the major eruption from the Sun. Their interplanetary counterparts, the interplanetary CMEs (ICMEs), are the direct manifestations of these structures when they propagate into the heliosphere and encounter one or more observing spacecraft. The ICMEs generally exhibit a set of distinctive signatures from the in-situ spacecraft measurements. A particular subset of ICMEs, the so-called Magnetic Clouds (MCs), is more uniquely defined and has been studied for decades, based on in-situ magnetic field and plasma measurements. By utilizing the latest multiple spacecraft measurements and analysis tools, we report a detailed study of the internal magnetic field configuration of an MC event observed by both the Solar Orbiter (SO) and Wind spacecraft in the solar wind near the Sun-Earth line. Both two-dimensional (2D) and three-dimensional (3D) models are applied to reveal the flux rope configurations of the MC. Various geometrical as well as physical parameters are derived and found to be similar within error estimates for the two methods. These results quantitatively characterize the coherent MC flux rope structure crossed by the two spacecraft along different paths. The implication for the radial evolution of this MC event is also discussed. △ Less

Submitted 4 July, 2021; originally announced July 2021.

Comments: Submitted to Frontiers in Physics, Research Topic: The Magnetic Structures and Their Role in The Evolution of Coronal Mass Ejections

The Inhomogeneity of Composition along the Magnetic Cloud Axis

Authors: Hongqiang Song, Qiang Hu, Xin Cheng, Jie Zhang, Leping Li, Ake Zhao, Bing Wang, Ruisheng Zheng, Yao Chen

Abstract: It is generally accepted that CMEs result from eruptions of magnetic flux ropes, which are dubbed as magnetic clouds in interplanetary space. The composition (including the ionic charge states and elemental abundances) is determined prior to and/or during CME eruptions in the solar atmosphere, and does not alter during magnetic cloud propagation to 1 AU and beyond. It has been known that the compo… ▽ More It is generally accepted that CMEs result from eruptions of magnetic flux ropes, which are dubbed as magnetic clouds in interplanetary space. The composition (including the ionic charge states and elemental abundances) is determined prior to and/or during CME eruptions in the solar atmosphere, and does not alter during magnetic cloud propagation to 1 AU and beyond. It has been known that the composition is not uniform within a cross section perpendicular to magnetic cloud axis, and the distribution of ionic charge states within a cross section provides us an important clue to investigate the formation and eruption processes of flux ropes due to the freeze-in effect. The flux rope is a three dimensional magnetic structure intrinsically, and it remains unclear whether the composition is uniform along the flux rope axis as most magnetic clouds are only detected by one spacecraft. In this paper we report a magnetic cloud that was observed by ACE near 1 AU on 1998 March 4--6 and Ulysses near 5.4 AU on March 24--28 sequentially. At these times, both spacecraft were located around the ecliptic plane, and the latitudinal and longitudinal separations between them were $\sim$2.2$^{\circ}$ and $\sim$5.5$^{\circ}$, respectively. It provides us an excellent opportunity to explore the axial inhomogeneity of flux rope composition, as both spacecraft almost intersected the cloud center at different sites along its axis. Our study shows that the average values of ionic charge states exhibit significant difference along the axis for carbon, and the differences are relatively slight but still obvious for charge states of oxygen and iron, as well as the elemental abundances of iron and helium. Besides the means, the composition profiles within the cloud measured by both spacecraft also exhibit some discrepancies. We conclude that the inhomogeneity of composition exists along the cloud axis. △ Less

Submitted 30 June, 2021; originally announced June 2021.

Comments: 15 pages, 4 figures, accepted by Frontiers in Physics

Small-scale Magnetic Flux Ropes with Field-aligned Flows via the PSP In-situ Observations

Authors: Yu Chen, Qiang Hu, Lingling Zhao, Justin C. Kasper, Jia Huang

Abstract: Magnetic flux rope, formed by the helical magnetic field lines, can sometimes remain its shape while carrying significant plasma flow that is aligned with the local magnetic field. We report the existence of such structures and static flux ropes by applying the Grad-Shafranov-based algorithm to the Parker Solar Probe (PSP) in-situ measurements in the first five encounters. These structures are det… ▽ More Magnetic flux rope, formed by the helical magnetic field lines, can sometimes remain its shape while carrying significant plasma flow that is aligned with the local magnetic field. We report the existence of such structures and static flux ropes by applying the Grad-Shafranov-based algorithm to the Parker Solar Probe (PSP) in-situ measurements in the first five encounters. These structures are detected at heliocentric distances, ranging from 0.13 to 0.66 au, in a total of 4-month time period. We find that flux ropes with field-aligned flows, although occur more frequently, have certain properties similar to those of static flux ropes, such as the decaying relations of the magnetic fields within structures with respect to heliocentric distances. Moreover, these events are more likely with magnetic pressure dominating over the thermal pressure. About one-third of events are detected in the relatively fast-speed solar wind. Taking into account the high Alfvenicity, we also compare with switchback spikes identified during three encounters and interpret their inter-relations. We find that some switchbacks can be detected when the spacecraft traverses flux rope-like structures. The cross-section maps for selected events are presented via the new Grad-Shafranov type reconstruction. Finally, the possible evolution of the magnetic flux rope structures in the inner heliosphere is discussed. △ Less

Submitted 17 May, 2021; v1 submitted 13 May, 2021; originally announced May 2021.

Comments: Accepted by ApJ on 2021 April 28

Deciphering Star Cluster Evolution by Shape Morphology

Authors: Qingshun Hu, Yu Zhang, Ali Esamdin, Jinzhong Liu, Xiangyun Zeng

Abstract: We analyze the morphological evolution of open clusters and provide shape parameters for 265 open clusters. The results show that the overall shape of sample clusters becomes more elliptical as they grow older, while their core remains circular or slightly trend to circularize. There is a negative correlation of the ellipticities with the number of members of the sample clusters. A significant neg… ▽ More We analyze the morphological evolution of open clusters and provide shape parameters for 265 open clusters. The results show that the overall shape of sample clusters becomes more elliptical as they grow older, while their core remains circular or slightly trend to circularize. There is a negative correlation of the ellipticities with the number of members of the sample clusters. A significant negative correlation between the overall ellipticities and masses is also detected for the sample clusters with log(age/year) $\geq$ 8, suggesting that the overall shapes of the clusters are possibly influenced by the number of members and masses, in addition to the external forces and the surrounding environment. For most young sample clusters, the radial stratification degree of the short axis direction is greater than that of the long, implying that the radial stratification degree in the two directions within the young sample cluster may be unevenly affected by an internal evolutionary process. Older sample clusters exhibit lower stratification in the tangential direction, which possibly means those clusters may continue to survive for a long time at a low level of stratification. Our analysis shows that the overall shape of the sample clusters may be more susceptible to the influence of Galactic tides toward the Galactic center than the shear forces embedded in Galactic differential rotation. By analyzing the distribution of the ages and number of members of star clusters, we suggest that NGC 6791 may originate from superclusters. △ Less

Submitted 5 April, 2021; v1 submitted 4 March, 2021; originally announced March 2021.