Showing 1–50 of 409 results for author: Liang, E

The white-light superflares from cool stars in GWAC triggers

Authors: Guang-Wei Li, Liang Wang, Hai-Long Yuan, Li-Ping Xin, Jing Wang, Chao Wu, Hua-Li Li, Hasitieer Haerken, Wei-Hua Wang, Hong-Bo Cai, Xu-Hui Han, Yang Xu, Lei Huang, Xiao-Meng Lu, Jian-Ying Bai, Xiang-Yu Wang, Zi-Gao Dai, En-Wei Liang, Jian-Yan Wei

Abstract: M-type stars are the ones that flare most frequently, but how big their maximum flare energy can reach is still unknown. We present 163 flares from 162 individual M2 through L1-type stars that triggered the GWAC, with flare energies ranging from $10^{32.2}$ to $10^{36.4}$ erg . The flare amplitudes range from $\triangle G = 0.84$ to $\sim 10$ mag. Flare energy increases with stellar surface temper… ▽ More M-type stars are the ones that flare most frequently, but how big their maximum flare energy can reach is still unknown. We present 163 flares from 162 individual M2 through L1-type stars that triggered the GWAC, with flare energies ranging from $10^{32.2}$ to $10^{36.4}$ erg . The flare amplitudes range from $\triangle G = 0.84$ to $\sim 10$ mag. Flare energy increases with stellar surface temperature ($T_{\rm eff}$) but both $\triangle G$ and equivalent duration $\log_{10}(ED)$ seem to be independent of $T_{\rm eff}$. Combining periods detected from light curves of TESS and K2, spectra from LAMOST, SDSS and the 2.16 m Telescope, and the Gaia DR3 data, we found that these GWAC flare stars are young. For the stars that have spectra, we found that these stars are in or very near to the saturation region, and $\log_{10}(L_{\rm Hα}/L_{\rm bol})$ is lower for M7-L1 stars than for M2-M6 stars. We also studied the relation between GWAC flare bolometric energy $E_{\rm bol}$ and stellar hemispherical area $S$, and found that $\log_{10}E_{\rm bol}$ (in erg) increases with increasing $S$ (in cm$^2$), and the maximum flare energy $\log_{10}E_{\rm bol, max} \geqslant \log_{10}S + 14.25$. For M7-L1 stars, there seem to be other factors limiting their maximum flare energies in addition to stellar hemispherical area. △ Less

Submitted 11 July, 2024; originally announced July 2024.

Comments: 18 pages, 11 figures, 4 tables

GRB 221009A/SN 2022xiw: A Supernova Obscured by a Gamma-Ray Burst Afterglow?

Authors: De-Feng Kong, Xiang-Gao Wang, WeiKang Zheng, Hou-Jun Lü, L. P. Xin, Da-Bin Lin, Jia-Xin Cao, Ming-Xuan Lu, B. Ren, Edgar P. Vidal, J. Y. Wei, En-Wei Liang, Alexei V. Filippenko

Abstract: We present optical photometry for the afterglow of GRB 221009A, in some respects the most extraordinary gamma-ray burst (GRB) ever observed. Good quality in the R-band light curve is obtained, covering 0.32-19.57 days since the Fermi-GBM trigger. We find that a weak bump emerges fromthe declining afterglow at $t \approx 11$ days; a supernova (SN) may be responsible. We use a smooth broken power-la… ▽ More We present optical photometry for the afterglow of GRB 221009A, in some respects the most extraordinary gamma-ray burst (GRB) ever observed. Good quality in the R-band light curve is obtained, covering 0.32-19.57 days since the Fermi-GBM trigger. We find that a weak bump emerges fromthe declining afterglow at $t \approx 11$ days; a supernova (SN) may be responsible. We use a smooth broken power-law and $^{56}\mathrm{Ni}$ model to fit the light curve. The best-fitting results reveal that the SN ejected a total mass of $M_\mathrm{ej} = 3.70 M_\odot$, a $^{56}\mathrm{Ni}$ mass of $M_\mathrm{Ni} = 0.23 M_\odot$, and a kinetic energy of $E_\mathrm{SN,K} = 2.35 \times 10^{52} \mathrm{erg}$. We also compare GRB 221009A with other GRB-SN events based on a GRB-associated SN sample, and find that only SN 2003lw and SN 2011kl can be obviously revealed in the afterglow of GRB 221009A by setting these objects at its distance. This suggests that a supernova (SN 2022xiw) is possibly obscured by the brighter afterglow emission from GRB 221009A. △ Less

Submitted 30 June, 2024; originally announced July 2024.

Extended GeV $γ$-ray emission around the star forming region of the W3 complex

Authors: Qihang Wu, Xiaona Sun, Ruizhi Yang, Tingting Ge, Yunfeng Liang, Enwei Liang

Abstract: We analyze the GeV $γ$-ray emission from the W3 complex using about 14 years of Pass 8 data recorded by the $\it Fermi$ Large Area Telescope (\textit{Fermi}-LAT). We resolve the $γ$-ray emissions around W3 into two components: an elliptical Gaussian overlapping with the molecular gas and a point-like source near the cluster W3 Main. The pion-bump feature of SED for the elliptical Gaussian together… ▽ More We analyze the GeV $γ$-ray emission from the W3 complex using about 14 years of Pass 8 data recorded by the $\it Fermi$ Large Area Telescope (\textit{Fermi}-LAT). We resolve the $γ$-ray emissions around W3 into two components: an elliptical Gaussian overlapping with the molecular gas and a point-like source near the cluster W3 Main. The pion-bump feature of SED for the elliptical Gaussian together with the better fitting result of pion decay model favor the hadronic origin. We further argue that the cosmic rays (CRs) could originate from the interactions between cluster winds and the shock produced by the SNR HB3. The point-like source positionally coincident with the star cluster W3 Main indicates it may be directly powered by near clusters, while its fainter $γ$-ray emissions below 10 GeV is possibly due to the shelter from dense gas making the low-energy CRs incapable of penetrating the dense materials. Meanwhile, we cannot rule out that the $γ$-ray emissions originate from the interaction of accelerated protons in SNR with the ambient gas. △ Less

Submitted 27 June, 2024; originally announced June 2024.

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

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

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

Submitted 12 June, 2024; originally announced June 2024.

Comments: 17 pages, 12 figures, accepted by PRL

Embedding-Aligned Language Models

Authors: Guy Tennenholtz, Yinlam Chow, Chih-Wei Hsu, Lior Shani, Ethan Liang, Craig Boutilier

Abstract: We propose a novel approach for training large language models (LLMs) to adhere to objectives defined within a latent embedding space. Our method leverages reinforcement learning (RL), treating a pre-trained LLM as an environment. Our embedding-aligned guided language (EAGLE) agent is trained to iteratively steer the LLM's generation towards optimal regions of the latent embedding space, w.r.t. so… ▽ More We propose a novel approach for training large language models (LLMs) to adhere to objectives defined within a latent embedding space. Our method leverages reinforcement learning (RL), treating a pre-trained LLM as an environment. Our embedding-aligned guided language (EAGLE) agent is trained to iteratively steer the LLM's generation towards optimal regions of the latent embedding space, w.r.t. some predefined criterion. We demonstrate the effectiveness of the EAGLE agent using the MovieLens 25M dataset to surface content gaps that satisfy latent user demand. We also demonstrate the benefit of using an optimal design of a state-dependent action set to improve EAGLE's efficiency. Our work paves the way for controlled and grounded text generation using LLMs, ensuring consistency with domain-specific knowledge and data representations. △ Less

Submitted 24 May, 2024; originally announced June 2024.

Spin, inclination, and magnetic field evolution of magnetar population in vacuum and plasma-filled magnetospheres

Authors: Jun-Xiang Huang, Hou-Jun Lü, Jared Rice, En-Wei Liang

Abstract: Magnetars are potential energy sources or central engines for numerous transient phenomena in the Universe. How newborn magnetars evolve in different environments remains an open question. Based on both observed and candidate magnetars, it is found that the periods of all magnetars or candidates appear as a bimodal distribution, and are defined as the ``long-P'' and ``short-P'' magnetar subclasses… ▽ More Magnetars are potential energy sources or central engines for numerous transient phenomena in the Universe. How newborn magnetars evolve in different environments remains an open question. Based on both observed and candidate magnetars, it is found that the periods of all magnetars or candidates appear as a bimodal distribution, and are defined as the ``long-P'' and ``short-P'' magnetar subclasses, respectively. We find that for the ``short-P'' subclass of magnetars, the $\dot{P}$ values also appear as a bimodal distribution, and therefore can be classified as ``high-$\dot{P}$ short-P'' and ``low-$\dot{P}$ short-P'' magnetar subclasses. In this paper, we use Monte Carlo simulations to generate synthetic magnetar populations and investigate the evolution of the ``high-$\dot{P}$ short-P'' and ``low-$\dot{P}$ short-P'' magnetar subclasses by considering both the magnetar spin and inclination, as well as the decay of their magnetic field within their evolution in both vacuum and plasma-filled magnetospheres. We find that the magnetar evolution is dependent on both spin and magnetic field, but seems to be insensitive to inclination evolution and magnetospheric environment for the ``high-$\dot{P}$ short-P'' sub-class. In comparison for the case of ``high-$\dot{P}$ short-P'', the magnetar evolution is dependent on spin, magnetic field, and inclination evolution, as well as the magnetospheric environment. The best evolution model should be the case of inclination evolution in vacuum with a small value of $\overline{\mathrm{FOM}}$. The differences in the best-fit parameters also suggest that the ``high-$\dot{P}$ short-P'' and ``low-$\dot{P}$ short-P'' magnetar subclasses may be tracking with different evolution channels. △ Less

Submitted 18 June, 2024; v1 submitted 24 May, 2024; originally announced May 2024.

Comments: 18 pages, 2 Tables, and 9 Figures. PRD in press, and matched with the published verison

Influence of quantum correction on the Schwarzschild black hole polarized image

Authors: Sen Guo, Yu-Xiang Huang, Kuan Liu, En-Wei Liang, Kai Lin

Abstract: Using a model of an accretion disk around a Schwarzschild black hole, the analytic estimates for image polarization were derived by Narayan $et~al.$. [Astrophys. J, 102, 912 (2021)]. Recently, the EHT team also obtained polarization images of the Sgr A$^{*}$ and measured both linear and circular polarization [Astrophys. J. Lett, 964, L25 (2024)]. We find that quantum correction effects can also in… ▽ More Using a model of an accretion disk around a Schwarzschild black hole, the analytic estimates for image polarization were derived by Narayan $et~al.$. [Astrophys. J, 102, 912 (2021)]. Recently, the EHT team also obtained polarization images of the Sgr A$^{*}$ and measured both linear and circular polarization [Astrophys. J. Lett, 964, L25 (2024)]. We find that quantum correction effects can also influence polarization information. Considering the quantum corrected Schwarzschild black hole (Kazakov-Solodukhin black hole), we derive the polarization intensity of the target black hole and investigate polarization images under different parameters. It is found that a larger quantum deformation leads to an expansion of the polarization region, while the polarization intensity value decrease. Under different observation angles, magnetic fields, fluid direction angles, and fluid velocity conditions, we also derive polarization images of corrected black holes. These key indicators not only affect the intensity of polarization but also the direction of polarization. We establish the relationship between polarization intensity and quantum correction deformation parameters, revealing a gradual decline in polarization intensity with reduced radius and an anti-polarization behavior induced by the progressive increase in deformation parameters at a constant radius. Our analysis may provide observational evidence for quantum effect of general relativity. △ Less

Submitted 21 May, 2024; originally announced May 2024.

Comments: 20 pages, 8 figures

Report number: Accepted European Physical Journal C (EPJC) 2024

SN 2019tua : A Type IIb Supernova with Multiple Bumps in the Light Curves

Authors: Xin-Bo Huang, Xiang-Gao Wang, Long Li, Li-Ping Xin, Jing Wang, Tian-Ci Zheng, Qi Wang, Hui-Ya Liu, Zi-Min Zhou, Xiao-meng Lu, jian-yan Wei, En-Wei Liang

Abstract: We present photometric and spectroscopic observations and analysis of the type IIb supernova (SN) SN 2019tua, which exhibits multiple bumps in its declining light curves between 40 and 65 days after discovery. SN 2019tua shows a time to peak of about 25 days similar to other type IIb SNe. Our observations indicate a decrease in its brightness of about 1 magnitude in the 60 days after the peak. At… ▽ More We present photometric and spectroscopic observations and analysis of the type IIb supernova (SN) SN 2019tua, which exhibits multiple bumps in its declining light curves between 40 and 65 days after discovery. SN 2019tua shows a time to peak of about 25 days similar to other type IIb SNe. Our observations indicate a decrease in its brightness of about 1 magnitude in the 60 days after the peak. At about days 50, and 60, its multiband light curves exhibit bumpy behavior. The complex luminosity evolution of SN 2019tua could not be well modeled with a single currently popular energy source model, e.g., radioactive decay of $^{56}$Ni, magnetar, interaction between the ejecta and a circumstellar shell. Even though the magnetar model has a smaller $χ^2 / \text{dof}$ value, the complex changes in SN 2019tua's brightness suggest that more than one physical process might be involved. We propose a hybrid CSM interaction plus $^{56}$Ni model to explain the bolometric light curve (LC) of SN 2019tua. The fitting results show that the ejecta mass $M_{\rm ej} \approx 2.4~M_\odot$, the total CSM mass $M_{\rm CSM} \approx 1.0~M_\odot$, and the $^{56}$Ni mass $M_{\rm Ni} \approx 0.4~M_\odot$. The total kinetic energy of the ejecta is $E_k\approx 0.5 \times 10^{51}\rm~erg$. Pre-existing multiple shells suggest that the progenitor of SN 2019tua experienced mass ejections within approximately $\sim6 - 44$ years prior to the explosion. △ Less

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

Comments: Accepted for publication in ApJ. 24 pages, 17 figures, 6 tables

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

Quantum criticality and Kibble-Zurek scaling in the Aubry-André-Stark model

Authors: En-Wen Liang, Ling-Zhi Tang, Dan-Wei Zhang

Abstract: We explore quantum criticality and Kibble-Zurek scaling (KZS) in the Aubry-Andre-Stark (AAS) model, where the Stark field of strength $\varepsilon$ is added onto the one-dimensional quasiperiodic lattice. We perform scaling analysis and numerical calculations of the localization length, inverse participation ratio (IPR), and energy gap between the ground and first excited states to characterize cr… ▽ More We explore quantum criticality and Kibble-Zurek scaling (KZS) in the Aubry-Andre-Stark (AAS) model, where the Stark field of strength $\varepsilon$ is added onto the one-dimensional quasiperiodic lattice. We perform scaling analysis and numerical calculations of the localization length, inverse participation ratio (IPR), and energy gap between the ground and first excited states to characterize critical properties of the delocalization-localization transition. Remarkably, our scaling analysis shows that, near the critical point, the localization length $ξ$ scales with $\varepsilon$ as $ξ\propto\varepsilon^{-ν}$ with $ν\approx0.3$ a new critical exponent for the AAS model, which is different from the counterparts for both the pure Aubry-Andre (AA) model and Stark model. The IPR $\mathcal{I}$ scales as $\mathcal{I}\propto\varepsilon^{s}$ with the critical exponent $s\approx0.098$, which is also different from both two pure models. The energy gap $ΔE$ scales as $ΔE\propto \varepsilon^{νz}$ with the same critical exponent $z\approx2.374$ as that for the pure AA model. We further reveal hybrid scaling functions in the overlap between the critical regions of the Anderson and Stark localizations. Moreover, we investigate the driven dynamics of the localization transitions in the AAS model. By linearly changing the Stark (quasiperiodic) potential, we calculate the evolution of the localization length and the IPR, and study their dependence on the driving rate. We find that the driven dynamics from the ground state is well described by the KZS with the critical exponents obtained from the static scaling analysis. When both the Stark and quasiperiodic potentials are relevant, the KZS form includes the two scaling variables. This work extends our understanding of critical phenomena on localization transitions and generalizes the application of the KZS to hybrid models. △ Less

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

Comments: 10 pages, 11 figures; accepted by PRB

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

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

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

Submitted 13 May, 2024; originally announced May 2024.

Comments: 11 pages, 5 figures

The Energy Sources, the Physical Properties, and the Mass-loss History of SN 2017dio

Authors: Deng-Wang Shi, Shan-Qin Wang, Wen-Pei Gan, En-Wei Liang

Abstract: We study the energy sources, the physical properties of the ejecta and the circumstellar medium (CSM), as well as the mass-loss history of the progenitor of SN 2017dio which is a broad-lined Ic (Ic-BL) supernova (SN) having unusual light curves (LCs) and signatures of hydrogen-rich CSM in its early spectrum. We find that the temperature of SN 2017dio began to increase linearly about 20 days after… ▽ More We study the energy sources, the physical properties of the ejecta and the circumstellar medium (CSM), as well as the mass-loss history of the progenitor of SN 2017dio which is a broad-lined Ic (Ic-BL) supernova (SN) having unusual light curves (LCs) and signatures of hydrogen-rich CSM in its early spectrum. We find that the temperature of SN 2017dio began to increase linearly about 20 days after the explosion. We use the $^{56}$Ni plus the ejecta-CSM interaction (CSI) model to fit the LCs of SN 2017dio, finding that the masses of the ejecta, the $^{56}$Ni, and the CSM are $\sim$ 12.41 M$_\odot$, $\sim$ 0.17 M$_\odot$, and $\sim$ 5.82 M$_\odot$, respectively. The early-time photosphere velocity and the kinetic energy of the SN are respectively {$\sim$ 1.89 $\times 10^4$ km s$^{-1}$} and $\sim$ 2.66 $\times 10^{52}$ erg, which are respectively comparable to those of SNe Ic-BL and hypernovae (HNe). We suggest that the CSM of SN 2017dio might be {from an luminous-blue-variable-like outburst or} pulsational pair instability $\sim$ 1.2$-$11.4 yr prior to the SN explosion{, or binary mass transfer}. {Moreover,} we find that its ejecta mass is larger than those of many SNe Ic-BL, and that its $^{56}$Ni mass ($M_{\rm Ni}$) is approximately equal to the mean (or median) value of $M_{\rm Ni}$ of SNe Ic-BL in the literature, but lower than $M_{\rm Ni}$ of prototype HNe (e.g., SN 1998bw and SN 2003dh). △ Less

Submitted 13 May, 2024; originally announced May 2024.

Comments: Accepted for publication in ApJ, 17 pages, 4 figures, 3 tables

Radar Fields: Frequency-Space Neural Scene Representations for FMCW Radar

Authors: David Borts, Erich Liang, Tim Brödermann, Andrea Ramazzina, Stefanie Walz, Edoardo Palladin, Jipeng Sun, David Bruggemann, Christos Sakaridis, Luc Van Gool, Mario Bijelic, Felix Heide

Abstract: Neural fields have been broadly investigated as scene representations for the reproduction and novel generation of diverse outdoor scenes, including those autonomous vehicles and robots must handle. While successful approaches for RGB and LiDAR data exist, neural reconstruction methods for radar as a sensing modality have been largely unexplored. Operating at millimeter wavelengths, radar sensors… ▽ More Neural fields have been broadly investigated as scene representations for the reproduction and novel generation of diverse outdoor scenes, including those autonomous vehicles and robots must handle. While successful approaches for RGB and LiDAR data exist, neural reconstruction methods for radar as a sensing modality have been largely unexplored. Operating at millimeter wavelengths, radar sensors are robust to scattering in fog and rain, and, as such, offer a complementary modality to active and passive optical sensing techniques. Moreover, existing radar sensors are highly cost-effective and deployed broadly in robots and vehicles that operate outdoors. We introduce Radar Fields - a neural scene reconstruction method designed for active radar imagers. Our approach unites an explicit, physics-informed sensor model with an implicit neural geometry and reflectance model to directly synthesize raw radar measurements and extract scene occupancy. The proposed method does not rely on volume rendering. Instead, we learn fields in Fourier frequency space, supervised with raw radar data. We validate the effectiveness of the method across diverse outdoor scenarios, including urban scenes with dense vehicles and infrastructure, and in harsh weather scenarios, where mm-wavelength sensing is especially favorable. △ Less

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

Comments: 8 pages, 6 figures, to be published in SIGGRAPH 2024

Revisiting the Constraint on Equation of State of Neutron Star based on the Binary Neutron Star Mergers

Authors: Yun-Peng Li, Zhi-Lin Chen, Da-Bin Lin, En-Wei Liang

Abstract: The merger of neutron star (NS)-NS binary can form different production of the compact remnant, among which the supramassive NS (SMNS) could create an internal plateau and the followed steep decay marks the collapse of the SMNS. The proportion of SMNS and the corresponding collapse-time are often used to constrain the NS equation of state (EoS). This paper revisits this topic by considering the ef… ▽ More The merger of neutron star (NS)-NS binary can form different production of the compact remnant, among which the supramassive NS (SMNS) could create an internal plateau and the followed steep decay marks the collapse of the SMNS. The proportion of SMNS and the corresponding collapse-time are often used to constrain the NS equation of state (EoS). This paper revisits this topic by considering the effect of an accretion disk on the compact remnant, which is not considered in previous works. Compared with previous works, the collapse-time distribution (peaks $\sim$100 s) of the SMNSs formed from NS-NS merger is almost unaffected by the initial surface magnetic ($B_{{\rm s},i}$) of NS, but the total energy output of the magnetic dipole radiation from the SMNSs depends on $B_{{\rm s},i}$ significantly. Coupling the constraints from the SMNS fraction, we exclude some EoSs and obtain three candidate EoSs, i.e., DD2, ENG, and MPA1. By comparing the distributions of the collapse-time and the luminosity of the internal plateau (in the short gamma-ray bursts) for those from observations with those obtained based on the three candidate EoSs, it is shown that only the EoS of ENG is favored. Our sample based on the ENG EOS and a mass distribution motivated by Galactic systems suggests that approximately $99\%$ of NS-NS mergers collapse to form a black hole within $10^7$s. This includes scenarios forming a BH promptly ($36.5\%$), a SMNS ($60.7\%$), or a stable NS that transitions into a BH or a SMNS following accretion ($2.1\%$). It also indicates that the remnants for GW170817 and GW190425, and the second object of GW190814 are more likely to be BHs. △ Less

Submitted 30 April, 2024; originally announced April 2024.

Soft X-ray prompt emission from a high-redshift gamma-ray burst EP240315a

Authors: Y. Liu, H. Sun, D. Xu, D. S. Svinkin, J. Delaunay, N. R. Tanvir, H. Gao, C. Zhang, Y. Chen, X. -F. Wu, B. Zhang, W. Yuan, J. An, G. Bruni, D. D. Frederiks, G. Ghirlanda, J. -W. Hu, A. Li, C. -K. Li, J. -D. Li, D. B. Malesani, L. Piro, G. Raman, R. Ricci, E. Troja , et al. (170 additional authors not shown)

Abstract: Long gamma-ray bursts (GRBs) are believed to originate from core collapse of massive stars. High-redshift GRBs can probe the star formation and reionization history of the early universe, but their detection remains rare. Here we report the detection of a GRB triggered in the 0.5--4 keV band by the Wide-field X-ray Telescope (WXT) on board the Einstein Probe (EP) mission, designated as EP240315a,… ▽ More Long gamma-ray bursts (GRBs) are believed to originate from core collapse of massive stars. High-redshift GRBs can probe the star formation and reionization history of the early universe, but their detection remains rare. Here we report the detection of a GRB triggered in the 0.5--4 keV band by the Wide-field X-ray Telescope (WXT) on board the Einstein Probe (EP) mission, designated as EP240315a, whose bright peak was also detected by the Swift Burst Alert Telescope and Konus-Wind through off-line analyses. At a redshift of $z=4.859$, EP240315a showed a much longer and more complicated light curve in the soft X-ray band than in gamma-rays. Benefiting from a large field-of-view ($\sim$3600 deg$^2$) and a high sensitivity, EP-WXT captured the earlier engine activation and extended late engine activity through a continuous detection. With a peak X-ray flux at the faint end of previously known high-$z$ GRBs, the detection of EP240315a demonstrates the great potential for EP to study the early universe via GRBs. △ Less

Submitted 25 April, 2024; originally announced April 2024.

Comments: 41 pages, 8 figures, 7 tables

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.

GeV gamma-ray emission in the field of young massive star cluster RCW 38

Authors: Ting-Ting Ge, Xiao-Na Sun, Rui-Zhi Yang, Pak-Hin Thomas Tam, Ming-Xuan Lu, En-Wei Liang

Abstract: We report the detection of gamma-ray emission by the Fermi Large Area Telescope (Fermi-LAT) towards the young massive star cluster RCW 38 in the 1-500 GeV photon energy range. We found spatially extended GeV emission towards the direction of RCW 38, which is best modelled by a Gaussian disc of 0.23$°$ radius with a significance of the extension is $\sim 11.4 σ$. Furthermore, the spatial correlatio… ▽ More We report the detection of gamma-ray emission by the Fermi Large Area Telescope (Fermi-LAT) towards the young massive star cluster RCW 38 in the 1-500 GeV photon energy range. We found spatially extended GeV emission towards the direction of RCW 38, which is best modelled by a Gaussian disc of 0.23$°$ radius with a significance of the extension is $\sim 11.4 σ$. Furthermore, the spatial correlation with the ionized and molecular gas content favors the hadronic origin of the gamma-ray emission. The gamma-ray spectrum of RCW 38 has a relatively hard photon index of $2.44 \pm 0.03$, which is similar to other young massive star clusters. We argue that the diffuse GeV gamma-ray emission in this region likely originates from the interaction of accelerated protons in the stellar cluster with the ambient gas. △ Less

Submitted 28 March, 2024; originally announced March 2024.

Comments: 8 pages, 5 figures, accepted for publication in MNRAS. arXiv admin note: text overlap with arXiv:2210.01352

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)

Narrowly-Banded Spectra with Peak Frequency Around 1 GHz of FRB 20201124A: Implications for Energy Function and Radiation Physics

Authors: Fen Lyu, En-Wei Liang, D. Li

Abstract: The radiation physics of fast radio bursts (FRBs) remains an open question. Current observations have discovered that narrowly-banded bursts of FRB 20201124A are active in 0.4-2 GHz and their spectral peak frequency ($ν^{\rm obs}_{p}$) are mostly toward $\sim 1$ GHz. Utilizing a sample of 1268 bursts of FRB 20201124A detected with the FAST telescope, we show that the $1σ$ spectral regime of 71.4\%… ▽ More The radiation physics of fast radio bursts (FRBs) remains an open question. Current observations have discovered that narrowly-banded bursts of FRB 20201124A are active in 0.4-2 GHz and their spectral peak frequency ($ν^{\rm obs}_{p}$) are mostly toward $\sim 1$ GHz. Utilizing a sample of 1268 bursts of FRB 20201124A detected with the FAST telescope, we show that the $1σ$ spectral regime of 71.4\% events (in-band bursts) is within the FAST bandpass. Their intrinsic burst energies ($E^{\rm obs}_{\rm BWe}$) and spectral widths ($σ_s^{\rm obs}$) are well measured by fitting the spectral profile with a Gaussian function. The derived $E^{\rm obs}_{\rm BWe}$ and $σ_s^{\rm obs}$ distributions are log-normal and centering at $\log E^{\rm obs}_{\rm BWe}/{\rm erg}=37.2~ (σ=0.76)$ and $\log σ_s^{\rm obs}/{\rm GHz}=-1.16~ (σ=0.17)$. Our Monte Carlo simulation analysis infers its intrinsic $ν_p$ distribution as a normal function centered at $ν_{p,c}=1.16$ GHz ($σ=0.22$) and its intrinsic energy function as $Φ(E)\propto E^{-0.60}e^{-E/E_c}$ with $E_c=9.49 \times 10^{37}$ erg. We compare these results with that of typical repeating FRBs 20121102A and 20190520B that are active over a broad frequency range at several specific frequencies and discuss possible observational biases on the estimation of the event rate and energy function. Based on these results, we argue that FRB 20201124A likely occurs in a fine-tuned plasma for maser radiations at a narrow frequency range, while FRB 20121102A and FRB 20190520B could involve clumpy plasma conditions that make maser emission around several specific frequencies in a broad range. △ Less

Submitted 13 March, 2024; originally announced March 2024.

Comments: 21 pages, 9 figures. Accepted for publication in The Astrophysical Journal

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

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

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

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

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

Updated kinematics of the Radcliffe Wave: non-synchronous, dipole-like vertical oscillations

Authors: Zhi-Kai Zhu, Min Fang, Zu-Jia Lu, Junzhi Wang, Guang-Xing Li, Shiyu Zhang, Veli-Matti Pelkonen, Paolo Padoan, En-Wei Liang

Abstract: The kinematic structure of the Radcliffe Wave (RW) is crucial for understanding its origin and evolution. In this work, we present an accurate measurement of the vertical velocity $V_Z$ by where the radial velocity (RV) measures are taken into consideration. This is achieved in two ways. First, the velocities are measured towards Young Stellar Objects (YSOs), using their RV and proper motion measu… ▽ More The kinematic structure of the Radcliffe Wave (RW) is crucial for understanding its origin and evolution. In this work, we present an accurate measurement of the vertical velocity $V_Z$ by where the radial velocity (RV) measures are taken into consideration. This is achieved in two ways. First, the velocities are measured towards Young Stellar Objects (YSOs), using their RV and proper motion measurements from APOGEE-2 and Gaia DR3. Second, we combine RV measurements toward clouds with proper motion measurements of associated YSOs to determine the vertical velocities of the clouds. The results reveal that the oscillations in $V_Z$ are not synchronous with the vertical coordinate. The difference is caused by a combination of the effect of the radial velocity which we include in this paper, and the difference in models. By supplementing our analysis with additional young star samples, we find a consistent dipole pattern in $V_Z$. The fact that no significant amplitude differences are found among the analyzed samples indicates that there is no apparent age gradient within the dipole. We propose that RW evolves at a relatively slow rate. The fact that it will take a much longer time for RW to complete a full period compared to the cloud lifetimes challenges its classification as a traditional "wave". This age discrepancy should explain the phase difference, and non-synchronous oscillation found in kinematic studies. △ Less

Submitted 23 February, 2024; originally announced February 2024.

Comments: 13 pages, 9 figures, submitted on 6 Feb 2024

Multiwavelength Polarization Observations of Mrk 501

Authors: Xin-Ke Hu, Yu-Wei Yu, Jin Zhang, Xiang-Gao Wang, Kishore C. Patra, Thomas G. Brink, Wei-Kang Zheng, Qi Wang, De-Feng Kong, Liang-Jun Chen, Ji-Wang Zhou, Jia-Xin Cao, Ming-Xuan Lu, Zi-Min Zhou, Yi-Ning Wei, Xin-Bo Huang, Xing-Lin Li, Hao Lou, Ji-Rong Mao, En-Wei Liang, Alexei V. Filippenko

Abstract: Mrk 501 is a prototypical high-synchrotron-peaked blazar (HBL) and serves as one of the primary targets for the {\it Imaging X-ray Polarimetry Explorer} ({\it IXPE}). In this study, we report X-ray polarization measurements of Mrk 501 based on six {\it IXPE} observations. The detection of X-ray polarization at a confidence level exceeding 99\% is achieved in four out of the six observations conduc… ▽ More Mrk 501 is a prototypical high-synchrotron-peaked blazar (HBL) and serves as one of the primary targets for the {\it Imaging X-ray Polarimetry Explorer} ({\it IXPE}). In this study, we report X-ray polarization measurements of Mrk 501 based on six {\it IXPE} observations. The detection of X-ray polarization at a confidence level exceeding 99\% is achieved in four out of the six observations conducted across the entire energy range (2--8 keV) of {\it IXPE}. The maximum polarization degree ($Π_{\rm X}$) is measured to be $15.8\%\pm2.8\%$, accompanied by a polarization angle ($ψ_{\rm X}$) of $98.0°\pm5.1°$ at a confidence level of $5.6 σ$. During the remaining two observations, only an upper limit of $Π_{\rm X}<$12\% could be derived at the 99\% confidence level. No temporal variability in polarization is observed throughout all six {\it IXPE} observations for Mrk 501. A discernible trend of energy-dependent variation in the polarization degree is detected in optical spectropolarimetry; however, no analogous indication is observed in $Π_{\rm X}$. The chromatic behavior of $Π$ and the consistent values of $ψ$ across different frequencies from X-rays to radio waves, along with the agreement between $ψ$ and jet position angle, strongly support the interpretation of the energy-stratified model with shock-accelerated particles in the jet of Mrk 501. Additionally, the possibility of the presence of a global helical magnetic field in the jet of Mrk 501 is discussed. △ Less

Submitted 3 July, 2024; v1 submitted 19 February, 2024; originally announced February 2024.

Comments: 20 pages, 8 figures, 4 tables, accepted for publication in ApJL

The progenitor and central engine of a peculiar GRB 230307A

Authors: ZhaoWei Du, HouJun Lü, Yong Yuan, Xing Yang, EnWei Liang

Abstract: Recently, a lack of supernova-associated long-duration gamma-ray burst (GRB 230307A) at such a low redshift $z=0.065$, but associated with a possible kilonova emission, has attracted great attention. Its heavy element nucleosynthesis and the characteristic of soft X-ray emission suggests that the central engine of GRB 230307A is magnetar which is originated from a binary compact star merger. The c… ▽ More Recently, a lack of supernova-associated long-duration gamma-ray burst (GRB 230307A) at such a low redshift $z=0.065$, but associated with a possible kilonova emission, has attracted great attention. Its heavy element nucleosynthesis and the characteristic of soft X-ray emission suggests that the central engine of GRB 230307A is magnetar which is originated from a binary compact star merger. The calculated lower value of $\varepsilon \sim 0.05$ suggests that the GRB 230307A seems to be with ambiguous progenitor. The lower value of $f_{\rm eff}=1.23$ implies that the GRB 230307A is not likely to be from the effect of "tip of iceberg". We adopt the magnetar central engine model to fit the observed soft X-ray emission with a varying efficiency and find that the parameters constraints of magnetar fall into a reasonable range, i.e., $B<9.4\times10^{15}$ G and $P<2.5$ ms for $Γ_{\rm sat} = 10^3$, and $B<3.6\times10^{15}$ G and $P<1.05$ ms for $Γ_{\rm sat} = 10^4$. Whether the progenitor of GBR 230307A is from the mergers of neutron star - white dwarf (NS - WD) or neutron star - neutron star (NS - NS) remains unknown. △ Less

Submitted 16 February, 2024; originally announced February 2024.

Comments: 7 pages, 1 table, 4 figures. Accepted for publication in ApJ Letters, 962:L27, 2024

Neutron Star Mergers as the Dominant Contributor to the Production of Heavy $r$-Process Elements

Authors: Meng-Hua Chen, Li-Xin Li, Qiu-Hong Chen, Rui-Chong Hu, En-Wei Liang

Abstract: The discovery of the radioactively powered kilonova AT2017gfo, associated with the short-duration gamma-ray burst GRB 170817A and the gravitational wave source GW170817, has provided the first direct evidence supporting binary neutron star mergers as crucial astrophysical sites for the synthesis of heavy elements beyond iron through $r$-process nucleosysthesis in the universe. However, recent iden… ▽ More The discovery of the radioactively powered kilonova AT2017gfo, associated with the short-duration gamma-ray burst GRB 170817A and the gravitational wave source GW170817, has provided the first direct evidence supporting binary neutron star mergers as crucial astrophysical sites for the synthesis of heavy elements beyond iron through $r$-process nucleosysthesis in the universe. However, recent identifications of kilonovae following long-duration gamma-ray bursts, such as GRB 211211A and GRB 230307A, has sparked discussions about the potential of neutron star-white dwarf mergers to also produce neutron-rich ejecta and contribute to the production of heavy $r$-process elements. In this work, we estimate the contribution of binary neutron star mergers to the total mass of $r$-process elements in the Milky Way and investigate the possibility of neutron star-white dwarf mergers as alternative astrophysical sites for $r$-process nucleosynthesis through an analysis of the total mass of the $r$-process elements in the Milky Way. Our results reveal that binary neutron star mergers can sufficiently account for the Galactic heavy $r$-process elements, suggesting that these events are the dominant contributor to the production of heavy $r$-process elements in the Milky Way. Considering the total mass of $r$-process elements in the Milky Way and the higher occurrence rate of neutron star-white dwarf mergers, it is unlikely that such mergers can produce a significant amount of neutron-rich ejecta, with the generated mass of $r$-process elements being lower than $0.005M_{\odot}$. △ Less

Submitted 12 February, 2024; originally announced February 2024.

Comments: 7 pages, 4 figures, Accepted for publication in MNRAS

Revisit the Emission Polarization of the Internal-shock for the blazars' Jet

Authors: Hao-Qiang Zhang, Da-Bin Lin, Kuan Liu, En-Wei Liang

Abstract: Recent Imaging X-ray Polarimetry Explorer (IXPE) observations of blazars tend to support the shock model for the X-ray emission, but report a low polarization degree ($Π\sim 10\%$) in X-rays compared with the previous theoretical expectations in the shock model. In order to reconcile the theoretical expectations with observations, we revisit the polarization of the shock emission by considering di… ▽ More Recent Imaging X-ray Polarimetry Explorer (IXPE) observations of blazars tend to support the shock model for the X-ray emission, but report a low polarization degree ($Π\sim 10\%$) in X-rays compared with the previous theoretical expectations in the shock model. In order to reconcile the theoretical expectations with observations, we revisit the polarization of the shock emission by considering different kind of direction distribution for the shock-generated magnetic fields (sgMFs). Here, $w'_{\rm sg}\propto(\sinθ')^{ζ_{\rm sg}}$ with $θ'=0$ along the shock normal direction is used to describe the direction distribution of the sgMFs in the shock co-moving frame. It is found that the polarization in the X-ray and radio emission for a general jet in blazars can be described as $Π\sim 44.5[1-\exp(-ζ_{\rm sg}/2.6)]\%$ and $Π\sim 20[1-\exp(-ζ_{\rm sg}/2.4)]\%$, respectively. Correspondingly, one can have $ζ_{\rm sg}\sim 1-1.5$ according to the IXPE observations. Besides the sgMFs, the magnetic fields generated by the Richmyer-Meshkov instability (rmMFs) is supposed to present in the jets. The direction of the rmMFs is mainly distributed along the shock normal in the simulations and thus $w'_{\rm rm}\propto(\cos θ')^{ζ_{\rm rm}}$ is adopted to describe the direction distribution of rmMFs. We find that the rmMFs is likely to significantly affect the polarization properties at the low-frequency emission, especially when the sgMFs decay rapidly. Based on the contemporaneous radio and X-ray observations, we find the the emission of the electrons in the rmMFs make a significant contribution in the low-frequency emission and the the ordered background magnetic fields (obMFs) can be neglected. △ Less

Submitted 10 February, 2024; originally announced February 2024.

Comprehensive Forecasting of California's Energy Consumption: A Multi-Source and Sectoral Analysis Using ARIMA and ARIMAX Models

Authors: Zahra Moslemi, Logan Clark, Sarah Kernal, Samantha Rehome, Scott Sprengel, Ahoora Tamizifar, Shawna Tuli, Vish Chokshi, Mo Nomeli, Ella Liang, Moury Bidgoli, Jeff Lu, Manish Dasaur, Marty Hodgett

Abstract: California's significant role as the second-largest consumer of energy in the United States underscores the importance of accurate energy consumption predictions. With a thriving industrial sector, a burgeoning population, and ambitious environmental goals, the state's energy landscape is dynamic and complex. This paper presents a comprehensive analysis of California's energy consumption trends an… ▽ More California's significant role as the second-largest consumer of energy in the United States underscores the importance of accurate energy consumption predictions. With a thriving industrial sector, a burgeoning population, and ambitious environmental goals, the state's energy landscape is dynamic and complex. This paper presents a comprehensive analysis of California's energy consumption trends and provides detailed forecasting models for different energy sources and sectors. The study leverages ARIMA and ARIMAX models, considering both historical consumption data and exogenous variables. We address the unique challenges posed by the COVID-19 pandemic and the limited data for 2022, highlighting the resilience of these models in the face of uncertainty. Our analysis reveals that while fossil fuels continue to dominate California's energy landscape, renewable energy sources, particularly solar and biomass, are experiencing substantial growth. Hydroelectric power, while sensitive to precipitation, remains a significant contributor to renewable energy consumption. Furthermore, we anticipate ongoing efforts to reduce fossil fuel consumption. The forecasts for energy consumption by sector suggest continued growth in the commercial and residential sectors, reflecting California's expanding economy and population. In contrast, the industrial sector is expected to experience more moderate changes, while the transportation sector remains the largest energy consumer. △ Less

Submitted 6 February, 2024; originally announced February 2024.

A Channel to Form Fast-spinning Black Hole-Neutron Star Binary Mergers as Multimessenger Sources. II. Accretion-induced Spin-up

Authors: Zhen-Han-Tao Wang, Rui-Chong Hu, Ying Qin, Jin-Ping Zhu, Bing Zhang, Shuang-Xi Yi, Qin-Wen Tang, Xin-Wen Shu, Fen Lyu, En-Wei Liang

Abstract: In this work, we investigate an alternative channel for the formation of fast-spinning black hole-neutron star (BHNS) binaries, in which super-Eddington accretion is expected to occur in accreting BHs during the stable mass transfer phase within BH-stripped helium (BH--He-rich) star binary systems. We evolve intensive \texttt{MESA} grids of close-orbit BH--He-rich star systems to systematically ex… ▽ More In this work, we investigate an alternative channel for the formation of fast-spinning black hole-neutron star (BHNS) binaries, in which super-Eddington accretion is expected to occur in accreting BHs during the stable mass transfer phase within BH-stripped helium (BH--He-rich) star binary systems. We evolve intensive \texttt{MESA} grids of close-orbit BH--He-rich star systems to systematically explore the projected aligned spins of BHs in BHNS binaries, as well as the impact of different accretion limits on the tidal disruption probability and electromagnetic (EM) signature of BHNS mergers. Most of the BHs in BHNS mergers cannot be effectively spun up through accretion, if the accretion rate is limited to $\lesssim10\,\dot{M}_{\rm Edd}$, where $\dot{M}_{\rm Edd}$ is the standard Eddington accretion limit. In order to reach high spins (e.g., $χ_{\rm BH} \gtrsim 0.5$), the BHs are required to be born less massive (e.g., $\lesssim3.0\,M_\odot$) in binary systems with initial periods of $\lesssim0.2-0.3\,{\rm days}$ and accrete material at $\sim100\,\dot{M}_{\rm Edd}$. However, even under this high accretion limit, $\gtrsim6\,M_\odot$ BHs are typically challenging to significantly spin up and generate detectable associated EM signals. Our population simulations suggest that different accretion limits have a slight impact on the ratio of tidal disruption events. However, as the accretion limit increases, the EM counterparts from the cosmological BHNS population can become bright overall. △ Less

Submitted 23 February, 2024; v1 submitted 30 January, 2024; originally announced January 2024.

Comments: 18 pages, 9 figures, 1 table, accepted for publication in ApJ

Ejecta-circumstellar medium interaction in high-density environment contribution to kilonova emission: Application to GRB 191019A

Authors: Suo-Ning Wang, Hou-Jun Lü, Yong Yuan, Hao-Yu Yuan, Jared Rice, Meng-Hua Chen, En-Wei Liang

Abstract: The nearby long-duration GRB 191019A recently detected by Swift lacks an associated supernova and belongs to a host galaxy with little star formation activity, suggesting that the origin of this burst is the result of a merger of two compact objects with dynamical interactions in a high-density medium of an active galactic nucleus. Given the potential motivation of this event, and given that it oc… ▽ More The nearby long-duration GRB 191019A recently detected by Swift lacks an associated supernova and belongs to a host galaxy with little star formation activity, suggesting that the origin of this burst is the result of a merger of two compact objects with dynamical interactions in a high-density medium of an active galactic nucleus. Given the potential motivation of this event, and given that it occurs in such a high-density environment, the ejecta-circumstellar medium (CSM) interaction cannot be ignored as possibly contributing to the kilonova emission. Here, we theoretically calculate the kilonova emission by considering the contribution of the ejecta-CSM interaction in a high-density environment. We find that the contribution to the kilonova emission from the ejecta-CSM interaction will dominate at a later time, and a smaller ejecta mass will have a stronger kilonova emission from the ejecta-CSM interaction. Moreover, we try to apply it to GRB 191019A, but we find that it is difficult to identify the possible kilonova emission from the observations, due to the contribution of the bright host galaxy. On the other hand, less injected mass (less than $M_{\rm ej}=2\times10^{-5}M_{\odot}$) will be required if one can detect the kilonova emission associated with a GRB 191019A-like event in the future. The {\em r}-process-powered and spin energy contributions from the magnetar are also discussed. △ Less

Submitted 20 February, 2024; v1 submitted 22 January, 2024; originally announced January 2024.

Comments: 10 pages, 6 figures, updated the rederence list. ApJ in press, and matched with the published verison

Non-symmetrical sparking may hint "zits'' on a pulsar surface

Authors: Zhengli Wang, Jiguang Lu, Jingchen Jiang, Shunshun Cao, Weiyang Wang, Enwei Liang, Renxin Xu

Abstract: Pulsar electrodynamics could be relevant to the physics of stellar surface, which remains poorly understood for more than half a centenary and is difficult to probe due to the absence of direct and clear observational evidence. Nevertheless, highly-sensitive telescopes (e.g., China's Five-hundred-meter Aperture Spherical radio Telescope, FAST) may play an essential role in solving the problem sinc… ▽ More Pulsar electrodynamics could be relevant to the physics of stellar surface, which remains poorly understood for more than half a centenary and is difficult to probe due to the absence of direct and clear observational evidence. Nevertheless, highly-sensitive telescopes (e.g., China's Five-hundred-meter Aperture Spherical radio Telescope, FAST) may play an essential role in solving the problem since the predicted surface condition would have quite different characteristics in some models of pulsar structure, especially after the establishment of the standard model of particle physics. For instance, small hills (or ``zit'') may exist on solid strangeon star surface with rigidity, preferential discharge, i.e., gap sparking, may occur around the hills in the polar cap region. In this work, with the 110-min polarization observation of PSR B0950+08 targeted by FAST, we report that the gap sparking is significantly non-symmetrical to the meridian plane on which the rotational and magnetic axes lie. It is then speculated that this asymmetry could be the result of preferential sparking around zits which might rise randomly on pulsar surface. Some polarization features of both single pulses and the mean pulse, as well as the cross-correlation function of different emission regions, have also been presented. △ Less

Submitted 25 January, 2024; originally announced January 2024.

Comments: Accepted for publication in Astronomische Nachrichten

The Jet Composition of GRB 230307A: Poynting-Flux-Dominated Outflow?

Authors: Zhao-Wei Du, HouJun Lü, Xiaoxuan Liu, EnWei Liang

Abstract: The jet composition of GRB plays an important role in understanding the energy dissipation and radiation mechanisms in GRB physics, but it is poorly constrained from the observational data. Recently, an interesting long-duration GRB 230307A with redshift $z=$0.065 has attracted great attention. The lack of detected thermal emission and mini-structure of prompt emission lightcurve of this burst sug… ▽ More The jet composition of GRB plays an important role in understanding the energy dissipation and radiation mechanisms in GRB physics, but it is poorly constrained from the observational data. Recently, an interesting long-duration GRB 230307A with redshift $z=$0.065 has attracted great attention. The lack of detected thermal emission and mini-structure of prompt emission lightcurve of this burst suggest that the outflow is Poynting-flux-dominated and point towards the ICMART model. In this paper, we invoke two independent methods to investigate the jet composition of GRB 230307A. The high magnetization parameter ($σ>7$ or ever large) for$R_0=10^{10}$ cm that is used to suppress thermal component, strongly suggests that a significant fraction of the outflow energy is likely in a Poynting flux entrained with the baryonic matter. Moreover, it is found that the radiation efficiency of this burst for typical values $ε_e=0.1$ and $ε_B=0.01$ can reach as high as $~50\%$ which disfavors the internal shock model, but is consistent with ICMART model. Finally, a possible unified picture to produce GRB 230307A originated from a compact star merger is also discussed. △ Less

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

Comments: 6 pages, 2 figures, 1 table, updated references, and matched with the published verison

X-Ray Polarization Variability of High Spectral Peak BL Lacertaes: Cases of 1ES 1959+650 and PKS 2155-304

Authors: Xin-Ke Hu, Yu-Wei Yu, Jin Zhang, Tan-Zheng Wu, Ji-Shun Lian, Xiang-Gao Wang, Hai-Ming Zhang, En-Wei Liang

Abstract: The high-energy-peaked BL Lacertae objects (HBLs) are the main targets of the Imaging X-ray Polarimetry Explorer (IXPE) for investigating the mechanisms of radiation and particle acceleration in jets. In this paper, we report the first IXPE observations of two HBLs, 1ES 1959+650 and PKS 2155--304. Both sources exhibit X-ray polarization with a confidence level exceeding 99\%, as well as significan… ▽ More The high-energy-peaked BL Lacertae objects (HBLs) are the main targets of the Imaging X-ray Polarimetry Explorer (IXPE) for investigating the mechanisms of radiation and particle acceleration in jets. In this paper, we report the first IXPE observations of two HBLs, 1ES 1959+650 and PKS 2155--304. Both sources exhibit X-ray polarization with a confidence level exceeding 99\%, as well as significant variability in polarization across different time intervals and energy ranges. Notably, PKS 2155--304 demonstrates the highest X-ray polarization among all blazars detected by IXPE within its entire energy band (2--8 keV), with a polarization degree of $Π_{\rm X}=21.9\%\pm1.9\%$ (MDP$_{99}\sim$6.0\%). An even higher polarization is observed in the 3--4 keV band, reaching $Π_{\rm X}=28.6\%\pm2.7\%$ (MDP$_{99}\sim$8.1\%) with a confidence level of 10.8$σ$. Furthermore, no polarization is detected above 5 keV energy band. For 1ES 1959+650, the highest detected polarization degree in the 2--8 keV band is $Π_{\rm X}=12.4\%\pm0.7\%$ (MDP$_{99}\sim$2.2\%), with an electric vector position angle (EVPA) of $ψ_{\rm X}=19.7^{\circ}\pm1.6^{\circ}$. The X-ray polarization of 1ES 1959+650 exhibits evident variability, accompanied by the variations of $ψ_{\rm X}$, flux, spectrum, and energy bin. We discuss possible implications of these observational findings, including the variability in polarization, rotation of EVPA, and transition between synchrotron and synchrotron-self-Compton. We speculate that the X-rays observed during different IXPE observations originate from distinct regions in the jet and may involve diverse mechanisms for particle acceleration. △ Less

Submitted 28 February, 2024; v1 submitted 4 December, 2023; originally announced December 2023.

Comments: The Astrophysical Journal Letters, In Press

A Spatially resolved X-ray Polarization map of the Vela Pulsar Wind Nebula

Authors: Kuan Liu, Fei Xie, Yi-Han Liu, Chi-Yung Ng, Niccolo' Bucciantini, Roger W. Romani, Martin C. Weisskopf, Enrico Costa, Alessandro Di Marco, Fabio La Monaca, Fabio Muleri, Paolo Soffitta, Wei Deng, Yu Meng, En-wei Liang

Abstract: In this paper, we present a full spatially resolved polarization map for the Vela Pulsar Wind Nebula (PWN) observed by IXPE. By employing effective background discrimination techniques, our results show a remarkably high degree of local polarization in the outskirt region, exceeding 60% (55%) with a probability of 95% (99%), which approaches the upper limit predicted by the synchrotron emission me… ▽ More In this paper, we present a full spatially resolved polarization map for the Vela Pulsar Wind Nebula (PWN) observed by IXPE. By employing effective background discrimination techniques, our results show a remarkably high degree of local polarization in the outskirt region, exceeding 60% (55%) with a probability of 95% (99%), which approaches the upper limit predicted by the synchrotron emission mechanism. The high degree of polarization suggests that the turbulent magnetic energy is at most 33% of the ordered one. In addition, the X-ray polarization map exhibits a toroidal magnetic field pattern that is consistent with the field revealed by radio observations across the entire nebula. This consistency reveals that the observed X-ray and radio emissions are radiated by electrons from the same magnetic field. Different from the Crab PWN, the consistency observed in the Vela PWN may be attributed to the interaction between the reverse shock of supernova blast wave and the PWN, which leads to a displacement between the synchrotron-cooled nebula and the fresh nebula close to the pulsar. These findings deepen our understanding of the structure and evolution of the Vela PWN, and the magnetohydrodynamic interaction in PWNe. △ Less

Submitted 21 November, 2023; originally announced November 2023.

Comments: Accepted by ApJL

Plasma Radiation Model of Fast Radio Bursts from Magnetars

Authors: Edison Liang

Abstract: We propose a novel idea for the coherent intense millisecond radio emission of cosmic fast radio bursts (FRBs), which have recently been identified with flares from a magnetar. Motivated by the conventional paradigm of Type III solar radio bursts, we will explore the emission of coherent plasma line radiation at the electron plasma frequency and its harmonic as potential candidates of the coherent… ▽ More We propose a novel idea for the coherent intense millisecond radio emission of cosmic fast radio bursts (FRBs), which have recently been identified with flares from a magnetar. Motivated by the conventional paradigm of Type III solar radio bursts, we will explore the emission of coherent plasma line radiation at the electron plasma frequency and its harmonic as potential candidates of the coherent FRB emissions associated with magnetar flares. We discuss the emissions region parameters in relativistic strongly magnetized plasmas consisting of electrons, positrons and protons. The goal is to make observable predictions of this model to confront the multi-wavelength observations of FRBs from magnetars. These results will impact both observational radio astronomy and space-based astrophysics △ Less

Submitted 15 November, 2023; originally announced November 2023.

Comments: 13 pages 7 figures

Radioactive Decay of Specific Heavy Elements as an Energy Source for Late-Time Kilonovae and Potential JWST Observations

Authors: Meng-Hua Chen, En-Wei Liang

Abstract: Revealing the temporal evolution of individual heavy elements synthesized in the merger ejecta from binary neutron star mergers not only improves our understanding of the origin of heavy elements beyond iron but also clarifies the energy sources of kilonovae. In this work, we present a comprehensive analysis of the temporal evolution of the energy fraction of each nuclide based on the $r$-process… ▽ More Revealing the temporal evolution of individual heavy elements synthesized in the merger ejecta from binary neutron star mergers not only improves our understanding of the origin of heavy elements beyond iron but also clarifies the energy sources of kilonovae. In this work, we present a comprehensive analysis of the temporal evolution of the energy fraction of each nuclide based on the $r$-process nucleosynthesis simulations. The heavy elements dominating the kilonova emission within $\sim100$~days are identified, including $^{127}$Sb, $^{128}$Sb, $^{129}$Sb, $^{130}$Sb, $^{129}$Te, $^{132}$I, $^{222}$Rn, $^{223}$Ra, $^{224}$Ra, and $^{225}$Ac. It is found that the late-time kilonova light curve ($t\gtrsim20$~days) is highly sensitive to the presence of the heavy element $^{225}$Ac (with a half-life of 10.0~days). Our analysis shows that the James Webb Space Telescope (JWST), with its high sensitivity in the near-infrared band, is a powerful instrument for the identification of these specific heavy elements. △ Less

Submitted 14 November, 2023; originally announced November 2023.

Comments: 7 pages, 8 figures, accepted by MNRAS

Multipole expansion of the gravitational field in a general class of fourth-order theories of gravity and the application in gyroscopic precession

Authors: Bofeng Wu, En-Wei Liang

Abstract: A viable weak-field and slow-motion approximation method is constructed in $F(R,R_{μν}R^{μν}, R_{μνρσ}R^{μνρσ})$ gravity, a general class of fourth-order theories of gravity. By applying this method, the metric, presented in the form of the multipole expansion, outside a spatially compact source up to $1/c^3$ order is provided, and the closed-form expressions for the source multipole moments are a… ▽ More A viable weak-field and slow-motion approximation method is constructed in $F(R,R_{μν}R^{μν}, R_{μνρσ}R^{μνρσ})$ gravity, a general class of fourth-order theories of gravity. By applying this method, the metric, presented in the form of the multipole expansion, outside a spatially compact source up to $1/c^3$ order is provided, and the closed-form expressions for the source multipole moments are all presented explicitly. The metric consists of the massless tensor part, the massive scalar part, and the massive tensor part, where the former is exactly the metric in General Relativity, and the latter two are the corrections to it. It is shown that the corrections bear the Yukawa-like dependence on the two massive parameters and predict the appearance of six additional sets of source multipole moments, which indicates that up to $1/c^3$ order, there exist six degrees of freedom beyond General Relativity within $F(R,R_{μν}R^{μν}, R_{μνρσ}R^{μνρσ})$ gravity. By means of the metric, for a gyroscope moving around the source without experiencing any torque, the multipole expansions of its spin's angular velocities of the Thomas precession, the geodetic precession, and the Lense-Thirring precession are derived, and from them, the corrections to the angular velocities of the three types of precession in General Relativity can be read off. These results indicate that differently from $f(R)$ or $f(R,\mathcal{G})$ gravity, the most salient feature of the general $F(R,R_{μν}R^{μν}, R_{μνρσ}R^{μνρσ})$ gravity is that it gives the nonvanishing correction to the gyroscopic spin's angular velocity of the Lense-Thirring precession in General Relativity. △ Less

Submitted 29 April, 2024; v1 submitted 13 November, 2023; originally announced November 2023.

Comments: 46 pages,0 figure,accepted by JCAP

Unveiling the unconventional optical signatures of regular black holes within accretion disk

Authors: Sen Guo, Yu-Xiang Huang, Yu-Hao Cui, Yan Han, Qing-Quan Jiang, En-Wei Liang, Kai Lin

Abstract: The optical characteristics of three types of black holes (BHs) surrounded by a thin accretion disk are discussed, namely the Schwarzschild BH, Bardeen BH, and Hayward BH. We calculate the deflection angle of light as it traverses the vicinity of each BH using numerical integration and semi-analytical methods, revealing that both approaches can effectively elucidate the deflection of light around… ▽ More The optical characteristics of three types of black holes (BHs) surrounded by a thin accretion disk are discussed, namely the Schwarzschild BH, Bardeen BH, and Hayward BH. We calculate the deflection angle of light as it traverses the vicinity of each BH using numerical integration and semi-analytical methods, revealing that both approaches can effectively elucidate the deflection of light around the BH. We investigate the optical appearance of the accretion disk and its corresponding observational images at various viewing angles, discovering that the luminosity in the region near the BH on the inner side of the accretion disk is higher than that on the outer side owing to higher material density in closer proximity to the BH. We observe a significant accumulation of brightness on the left side of the accretion disk, attributed to the motion of matter and geometric effects. Our findings emphasize the significant influence of the observation inclination angle on the observed outcomes. An increase in the observation inclination angle results in the separation of higher-order images. With the improvement in EHT observation accuracy, we believe that the feature of a minimal distance between the innermost region of the direct image of the Hayward BH and the outermost region of the secondary image can be used as an indicator for identifying Hayward BHs. △ Less

Submitted 31 October, 2023; originally announced October 2023.

Comments: 27 pages, 10 figures

Journal ref: Published European Physical Journal C (2023)

Searching for Associations Between Short Gamma-ray Bursts and Fast Radio Burst

Authors: Ming-Xuan Lu, Long Li, Xiang-Gao Wang, Can-Min Deng, Yun-Feng Liang, Da-Bin Lin, En-Wei Liang

Abstract: The physical origin of fast radio bursts (FRBs) is still unclear. However, young magnetars associated with short-duration gamma-ray bursts (SGRBs) have been thought to be possible central engines for some FRBs. In this paper, we perform a systematic search for SGRBs that are associated with FRBs in a sample including 623 FRBs (601 one-off bursts and 22 repeaters) and 168 SGRBs with precise localiz… ▽ More The physical origin of fast radio bursts (FRBs) is still unclear. However, young magnetars associated with short-duration gamma-ray bursts (SGRBs) have been thought to be possible central engines for some FRBs. In this paper, we perform a systematic search for SGRBs that are associated with FRBs in a sample including 623 FRBs (601 one-off bursts and 22 repeaters) and 168 SGRBs with precise localizations. We find that FRB 190309A is spatially associated with GRB 060502B, with a chance probability of 0.05 when temporal and redshift information is taken into account. Considering the high chance probability (the statistical significance is < 3σ), we examine other observational properties such as the host galaxy, the dispersion measure, and the energy budget of the central engine to check the possibility of their association. Although the available observational information is insufficient to determine whether they are physically associated, it does not rule out such a possibility. As the only pair of FRB and GRB that are spatially associated, it remains an interesting case worthy of further attention △ Less

Submitted 30 October, 2023; originally announced October 2023.

Comments: 19 pages, 4 figures, accepted to ApJ

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

Density Functional Theory Calculations on the Interstellar Formation of Biomolecules

Authors: Qingli Liao, Junzhi Wang, Peng Xie, Enwei Liang, Zhao Wang

Abstract: The density functional theory (DFT) is the most versatile electronic structure method used in quantum chemical calculations, and is increasingly applied in astrochemical research. This mini-review provides an overview of the applications of DFT calculations in understanding the chemistry that occurs in star-forming regions. We survey investigations into the formation of biologically-relevant compo… ▽ More The density functional theory (DFT) is the most versatile electronic structure method used in quantum chemical calculations, and is increasingly applied in astrochemical research. This mini-review provides an overview of the applications of DFT calculations in understanding the chemistry that occurs in star-forming regions. We survey investigations into the formation of biologically-relevant compounds such as nucleobases in the interstellar medium, and also covers the formation of both achiral and chiral amino acids, as well as biologically-relevant molecules such as sugars, and nitrogen-containing polycyclic aromatic hydrocarbons. Additionally, DFT calculations are used to estimate the potential barriers for chemical reactions in astronomical environments. We conclude by noting several areas that require more research, such as the formation pathways of chiral amino acids, complex sugars and other biologically-important molecules, and the role of environmental factors in the formation of interstellar biomolecules. △ Less

Submitted 22 October, 2023; originally announced October 2023.

Maximum Number of Quads

Authors: Nikhil Byrapuram, Hwiseo, Choi, Adam Ge, Selena Ge, Tanya Khovanova, Sylvia Zia Lee, Evin Liang, Rajarshi Mandal, Aika Oki, Daniel Wu, Michael Yang

Abstract: We study the maximum number of quads among $\ell$ cards from an EvenQuads deck of size $2^n$. This corresponds to enumerating quadruples of integers in the range $[0,\ell-1]$ such that their bitwise XOR is zero. In this paper, we conjecture a formula that calculates the maximum number of quads among $\ell$ cards. We study the maximum number of quads among $\ell$ cards from an EvenQuads deck of size $2^n$. This corresponds to enumerating quadruples of integers in the range $[0,\ell-1]$ such that their bitwise XOR is zero. In this paper, we conjecture a formula that calculates the maximum number of quads among $\ell$ cards. △ Less

Submitted 14 October, 2023; originally announced October 2023.

Comments: 13 pages. 2 tables

MSC Class: 05A99; 00A08

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

EvenQuads Game and Error-Correcting Codes

Authors: Nikhil Byrapuram, Hwiseo, Choi, Adam Ge, Selena Ge, Tanya Khovanova, Sylvia Zia Lee, Evin Liang, Rajarshi Mandal, Aika Oki, Daniel Wu, Michael Yang

Abstract: EvenQuads is a new card game that is a generalization of the SET game, where each card is characterized by three attributes, each taking four possible values. Four cards form a quad when, for each attribute, the values are the same, all different, or half and half. Given $\ell$ cards from the deck of EvenQuads, we can build an error-correcting linear binary code of length $\ell$ and Hamming distan… ▽ More EvenQuads is a new card game that is a generalization of the SET game, where each card is characterized by three attributes, each taking four possible values. Four cards form a quad when, for each attribute, the values are the same, all different, or half and half. Given $\ell$ cards from the deck of EvenQuads, we can build an error-correcting linear binary code of length $\ell$ and Hamming distance 4. The quads correspond to codewords of weight 4. Error-correcting codes help us calculate the possible number of quads when given up to 8 cards. We also estimate the number of cards that do not contain quads for decks of different sizes. In addition, we discuss properties of error-correcting codes built on semimagic, magic, and strongly magic quad squares. △ Less

Submitted 2 October, 2023; originally announced October 2023.

Comments: 20 pages, 8 tables

MSC Class: 05A99; 94B05; 00A08

In-orbit background and sky survey simulation study of POLAR-2/LPD

Authors: Zu-Ke Feng, Hong-Bang Liu, Fei Xie, Huan-Bo Feng, Qian-Nan Mai, Jiang-Chuan Tuo, Qian Zhong, Jian-Chao Sun, Jiang He, Yuan-Hao Wang, Qian Liu, Di-Fan Yi, Rui-Ting Ma, Bin-Long Wang, Zhen-Yu Tang, Shuang-Nan Zhang, En-Wei Liang

Abstract: The Low-Energy X-ray Polarization Detector (LPD) is one of the payloads in the POLAR-2 experiment, designed as an external payload for the China Space Station (CSS) deployment in early 2024. LPD is specifically designed to observe the polarization of Gamma-Ray Bursts (GRBs) prompt emission in the energy range of 2-10 keV, with a wide field of view (FoV) of 90 degrees in preliminary design. This ob… ▽ More The Low-Energy X-ray Polarization Detector (LPD) is one of the payloads in the POLAR-2 experiment, designed as an external payload for the China Space Station (CSS) deployment in early 2024. LPD is specifically designed to observe the polarization of Gamma-Ray Bursts (GRBs) prompt emission in the energy range of 2-10 keV, with a wide field of view (FoV) of 90 degrees in preliminary design. This observation is achieved using an array of X-ray photoelectric polarimeters based on gas pixel detectors. Due to the wide FoV configuration, the in-orbit background count rate in the soft X-ray range is high, while GRBs themselves also exhibit a high flux in this energy band. In order to assess the contribution of various background components to the total count rate, we conducted detailed simulations using the GEANT4 C++ package. Our simulations encompassed the main interactions within the instrument materials and provided insights into various background components within the wide FoV scheme. The simulation results reveal that among the background components, the primary contributors are the cosmic X-ray background (CXB) and bright X-ray sources. The total background count rate of LPD, after applying the charged particle background rejection algorithm, is approximately 0.55 counts/cm^2/s on average, and it varies with the detector's orbit and pointing direction. Furthermore, we performed comprehensive simulations and comparative analyses of the CXB and X-ray bright sources under different FoVs and detector pointings. These analyses provide valuable insights into the background characteristic for soft X-ray polarimeter with wide FoV. △ Less

Submitted 1 October, 2023; originally announced October 2023.

Constraining the Charge of a Black Hole with Electromagnetic Radiation from a Black Hole-Neutron Star System

Authors: Hao-Yu Yuan, Hou-Jun Lü, Jared Rice, En-Wei Liang

Abstract: Black hole-neutron star (BH-NS) mergers are expected to emit gravitational-wave (GW) and electromagnetic (EM) counterparts when the NS is tidally disrupted or plunges into the BH. Recently, GW 200105 and GW200115 were claimed as originating in BH-NS mergers, even GW 200105 remains in debate. Several optical source candidates are reported to possible associate with the two GW events, but not confir… ▽ More Black hole-neutron star (BH-NS) mergers are expected to emit gravitational-wave (GW) and electromagnetic (EM) counterparts when the NS is tidally disrupted or plunges into the BH. Recently, GW 200105 and GW200115 were claimed as originating in BH-NS mergers, even GW 200105 remains in debate. Several optical source candidates are reported to possible associate with the two GW events, but not confirmed yet. In this work, we assume that the BH is charged (the NS is naturally charged) and try to constrain the charge of the BH by using the possible associated EM emission from the charged BH and NS system working in the inspiral regime. We adopt electric and magnetic dipole radiations for the binaries which power a Poynting-flux-dominated outflow to accelerate electrons. Then, it produces the observed EM radiation via synchrotron radiation. We find that the conversion efficiency in the X-ray band is much higher than that of the ultraviolet (UV), near-infrared, and radio bands. The estimated maximum charge-to-mass ratio (the charge for unit mass) of the BH is $1.12\times 10^{-6}$ and $1.53\times 10^{-6}$ esu for the binary systems of GW200105 and GW200115, respectively, if magnetic field strength $B_{p}\lesssim ~10^{16}$ G and period $P>~1$ ms for the NS spin. △ Less

Submitted 24 September, 2023; originally announced September 2023.

Comments: 13 pages, 4 Tables, and 2 Figures. Accepted for publication in PRD

Influences of tilted thin accretion disks on the observational appearance of hairy black holes in Horndeski gravity

Authors: Shiyang Hu, Dan Li, Chen Deng, Xin Wu, Enwei Liang

Abstract: Research on the observational appearance of black holes, both in general relativity and modified gravity, has been in full swing since the Event Horizon Telescope Collaboration announced photos of M87$^{*}$ and Sagittarius A$^{*}$. Nevertheless, limited attention has been given to the impact of tilted accretion disks on black hole images. This paper investigates the $230$ GHz images of non-rotatin… ▽ More Research on the observational appearance of black holes, both in general relativity and modified gravity, has been in full swing since the Event Horizon Telescope Collaboration announced photos of M87$^{*}$ and Sagittarius A$^{*}$. Nevertheless, limited attention has been given to the impact of tilted accretion disks on black hole images. This paper investigates the $230$ GHz images of non-rotating hairy black holes illuminated by tilted, thin accretion disks in Horndeski gravity with the aid of a ray tracing method. The results indicate that reducing the scalar hair parameter effectively diminishes image luminosity and extends both the critical curve and the inner shadow. This trend facilitates the differentiation between hairy black holes and Schwarzschild black holes, especially in certain parameter spaces where the current Event Horizon Telescope array is capable of capturing such variations. Furthermore, we observe that the inclination of the tilted accretion disk can mimic the observation angle, consequently affecting image brightness and the morphology of the inner shadow. In specific parameter spaces, alterations in the tilt or position of the accretion disk can lead to a drift in the light spot within the images of hairy black holes. This finding may establish a potential correlation between the precession of the tilted accretion disk and image features. Additionally, through an examination of images depicting hairy black holes surrounded by two thin accretion disks, we report the obscuring effect of the accretion environment on the inner shadow of the black hole. △ Less

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

Comments: 32 pages, 17 figures, accepted for publication in JCAP

Super-Eddington Accretion as a Possible Scenario to Form GW190425

Authors: W. T. Zhang, Z. H. T. Wang, J. -P. Zhu, R. -C. Hu, X. W. Shu, Q. W. Tang, S. X. Yi, F. Lyu, E. W. Liang, Y. Qin

Abstract: On 2019 April 25, the LIGO/Virgo Scientific Collaboration detected a compact binary coalescence, GW190425. Under the assumption of the binary neutron star (BNS), the total mass of $3.4^{+0.3}_{-0.1}\, M_\odot$ lies five standard deviations away from the known Galactic population mean. In the standard common envelope scenario, the immediate progenitor of GW190425 is a close binary system composed o… ▽ More On 2019 April 25, the LIGO/Virgo Scientific Collaboration detected a compact binary coalescence, GW190425. Under the assumption of the binary neutron star (BNS), the total mass of $3.4^{+0.3}_{-0.1}\, M_\odot$ lies five standard deviations away from the known Galactic population mean. In the standard common envelope scenario, the immediate progenitor of GW190425 is a close binary system composed of an NS and a He-rich star. With the detailed binary evolutionary modeling, we find that in order to reproduce GW190425-like events, super-Eddington accretion (e.g., $1,000\,\dot{M}_{\rm Edd}$) from a He-rich star onto the first-born NS with a typical mass of 1.33 $M_\odot$ via stable Case BB mass transfer (MT) is necessarily required. Furthermore, the immediate progenitors should potentially have an initial mass of $M_{\rm ZamsHe}$ in a range of $3.0-3.5$ $M_\odot$ and an initial orbital period of $P_{\rm init}$ from 0.08 days to 0.12 days, respectively. The corresponding mass accreted onto NSs via stable Case BB MT phase varies from $0.70\, M_\odot$ to $0.77\, M_\odot$. After the formation of the second-born NS, the BNSs are expected to be merged due to gravitational wave emission from $\sim$ 11 Myr to $\sim$ 190 Myr. △ Less

Submitted 28 September, 2023; v1 submitted 10 September, 2023; originally announced September 2023.

Comments: 8 pages, 10 figures, updated to add one reference

Search for the gamma-ray spectral lines with the DAMPE and the Fermi-LAT observations

Authors: Ji-Gui Cheng, Yun-Feng Liang, En-Wei Liang

Abstract: Weakly interacting massive particles, as a major candidate of dark matter (DM), may directly annihilate or decay into high-energy photons, producing monochromatic spectral lines in the gamma-ray band. These spectral lines, if detected, are smoking-gun signatures for the existence of new physics. Using the 5 years of DAMPE and 13 years of Fermi-LAT data, we search for line-like signals in the energ… ▽ More Weakly interacting massive particles, as a major candidate of dark matter (DM), may directly annihilate or decay into high-energy photons, producing monochromatic spectral lines in the gamma-ray band. These spectral lines, if detected, are smoking-gun signatures for the existence of new physics. Using the 5 years of DAMPE and 13 years of Fermi-LAT data, we search for line-like signals in the energy range of 3 GeV to 1 TeV from the Galactic halo. Different regions of interest are considered to accommodate different DM density profiles. We do not find any significant line structure, and the previously reported line-like feature at $\sim$133 GeV is also not detected in our analysis. Adopting a local DM density of $ρ_{\rm local}=0.4\,{\rm GeV\,cm^{-3}}$, we derive 95% confidence level constraints on the velocity-averaged cross-section of $\langle{σv}\rangle_{γγ} \lesssim 4 \times 10^{-28}\,{\rm cm^{3}\,s^{-1}}$ and the decay lifetime of $τ_{γν} \gtrsim 5 \times 10^{29}\,{\rm s}$ at 100 GeV, achieving the strongest constraints to date for the line energies of 6-660 GeV. The improvement stems from the longer Fermi-LAT data set used and the inclusion of DAMPE data in the analysis. The simultaneous use of two independent data sets could also reduce the systematic uncertainty of the search. △ Less

Submitted 31 August, 2023; originally announced August 2023.

Comments: 14 pages, 10 figures. Accepted for publication in PRD

Journal ref: Phys. Rev. D 108, 063015 (2023)

Studying X-ray spectra from large-scale jets of FR II radio galaxies: application of shear particle acceleration

Authors: Jia-Chun He, Xiao-Na Sun, Jie-Shuang Wang, Frank M. Rieger, Ruo-Yu Liu, En-Wei Liang

Abstract: Shear particle acceleration is a promising candidate for the origin of extended high-energy emission in extra-galactic jets. In this paper, we explore the applicability of a shear model to 24 X-ray knots in the large-scale jets of FR II radio galaxies, and study the jet properties by modeling the multi-wavelength spectral energy distributions (SEDs) in a leptonic framework including synchrotron an… ▽ More Shear particle acceleration is a promising candidate for the origin of extended high-energy emission in extra-galactic jets. In this paper, we explore the applicability of a shear model to 24 X-ray knots in the large-scale jets of FR II radio galaxies, and study the jet properties by modeling the multi-wavelength spectral energy distributions (SEDs) in a leptonic framework including synchrotron and inverse Compton - CMB processes. In order to improve spectral modelling, we analyze Fermi-LAT data for five sources and reanalyzed archival data of Chandra on 15 knots, exploring the radio to X-ray connection. We show that the X-ray SEDs of these knots can be satisfactorily modelled by synchrotron radiation from a second, shear-accelerated electron population reaching multi-TeV energies. The inferred flow speeds are compatible with large-scale jets being mildly relativistic. We explore two different shear flow profiles (i.e., linearly decreasing and power-law) and find that the required spine speeds differ only slightly, supporting the notion that for higher flow speeds the variations in particle spectral indices are less dependent on the presumed velocity profile. The derived magnetic field strengths are in the range of a few to ten microGauss, and the required power in non-thermal particles typically well below the Eddington constraint. Finally, the inferred parameters are used to constrain the potential of FR II jets as possible UHECR accelerators. △ Less

Submitted 22 August, 2023; originally announced August 2023.

Comments: 14 pages, 4 figures, 6 table, accepted for publication in MNRAS

Revisiting the Properties of GW190814 and Its Formation History

Authors: F. Lyu, L. Yuan, D. H. Wu, W. H. Guo, Y. Z. Wang, S. X. Yi, Q. W. Tang, R. -C. Hu, J. -P. Zhu, X. W. Shu, Y. Qin, E. W. Liang

Abstract: GW190814 was reported during LIGO's and Virgo's third observing run with the most asymmetric component masses (a $\sim 23$ $M_{\odot}$ black hole and a $\sim2.6$ $M_{\odot}$ compact object). Under the assumption that this event is a binary black hole (BBH) merger formed through the isolated binary evolution channel, we reanalyze the publicly released data of GW190814 with the modified astrophysica… ▽ More GW190814 was reported during LIGO's and Virgo's third observing run with the most asymmetric component masses (a $\sim 23$ $M_{\odot}$ black hole and a $\sim2.6$ $M_{\odot}$ compact object). Under the assumption that this event is a binary black hole (BBH) merger formed through the isolated binary evolution channel, we reanalyze the publicly released data of GW190814 with the modified astrophysical priors on the effective spin $χ_{\rm eff}$, and further explore its formation history using detailed binary modeling. We show that GW190814 is likely to have been formed through the classical common envelope channel. Our findings show that the properties inferred using the modified astrophysical priors are consistent with those inferred by the uniform priors. With the newly-inferred properties of GW190814, we perform detailed binary evolution of the immediate progenitor of the BBH (namely a close binary system composed of a BH and a helium star) in a large parameter space, taking into account mass-loss, internal differential rotation, supernova kicks, and tidal interactions between the helium star and the BH companion. Our findings show that GW190814-like events could be formed in limited initial conditions just after the common envelope phase: a $\sim 23$ $M_{\odot}$ BH and a helium star of $M_{\rm ZamsHe}$ $\sim$ 8.5 $M_{\odot}$ at solar metallicity ($\sim$ 7.5 $M_{\odot}$ at 10\% solar metallicity) with an initial orbital period at around 1.0 day. Additionally, the inferred low spin of the secondary indicates that the required metallicity for reproducing GW190814-like events should not be too low (e.g., Z $\gtrsim$ 0.1 $Z_{\odot}$). △ Less

Submitted 3 September, 2023; v1 submitted 18 August, 2023; originally announced August 2023.

Comments: 8 pages, 8 figures, 1 table, final version published in MNRANS

Atypical radio pulsations from magnetar SGR 1935+2154

Authors: Pei Wang, Jian Li, Long Ji, Xian Hou, Erbil Gugercinoglu, Di Li, Diego F. Torres, Yutong Chen, Jiarui Niu, Weiwei Zhu, Bing Zhang, En-wei Liang, Li Zhang, Mingyu Ge, Zigao Dai, Lin Lin, Jinlin Han, Yi Feng, Chenhui Niu, Yongkun Zhang, Dengjiang Zhou, Heng Xu, Chunfeng Zhang, Jinchen Jiang, Chenchen Miao , et al. (10 additional authors not shown)

Abstract: Magnetars are neutron stars with extremely strong magnetic fields, frequently powering high-energy activity in X-rays. Pulsed radio emission following some X-ray outbursts have been detected, albeit its physical origin is unclear. It has long been speculated that the origin of magnetars' radio signals is different from those from canonical pulsars, although convincing evidence is still lacking. Fi… ▽ More Magnetars are neutron stars with extremely strong magnetic fields, frequently powering high-energy activity in X-rays. Pulsed radio emission following some X-ray outbursts have been detected, albeit its physical origin is unclear. It has long been speculated that the origin of magnetars' radio signals is different from those from canonical pulsars, although convincing evidence is still lacking. Five months after magnetar SGR 1935+2154's X-ray outburst and its associated Fast Radio Burst (FRB) 20200428, a radio pulsar phase was discovered. Here we report the discovery of X-ray spectral hardening associated with the emergence of periodic radio pulsations from SGR 1935+2154 and a detailed analysis of the properties of the radio pulses. The complex radio pulse morphology, which contains both narrow-band emission and frequency drifts, has not been seen before in other magnetars, but is similar to those of repeating FRBs - even though the luminosities are many orders of magnitude different. The observations suggest that radio emission originates from the outer magnetosphere of the magnetar, and the surface heating due to the bombardment of inward-going particles from the radio emission region is responsible for the observed X-ray spectral hardening. △ Less

Submitted 17 August, 2023; originally announced August 2023.

Comments: 47 pages, 11 figures