-
The Giant Radio Array for Neutrino Detection (GRAND) Collaboration -- Contributions to the 10th International Workshop on Acoustic and Radio EeV Neutrino Detection Activities (ARENA 2024)
Authors:
Rafael Alves Batista,
Aurélien Benoit-Lévy,
Teresa Bister,
Martina Bohacova,
Mauricio Bustamante,
Washington Carvalho,
Yiren Chen,
LingMei Cheng,
Simon Chiche,
Jean-Marc Colley,
Pablo Correa,
Nicoleta Cucu Laurenciu,
Zigao Dai,
Rogerio M. de Almeida,
Beatriz de Errico,
Sijbrand de Jong,
João R. T. de Mello Neto,
Krijn D de Vries,
Valentin Decoene,
Peter B. Denton,
Bohao Duan,
Kaikai Duan,
Ralph Engel,
William Erba,
Yizhong Fan
, et al. (100 additional authors not shown)
Abstract:
This is an index of the contributions by the Giant Radio Array for Neutrino Detection (GRAND) Collaboration to the 10th International Workshop on Acoustic and Radio EeV Neutrino Detection Activities (ARENA 2024, University of Chicago, June 11-14, 2024). The contributions include an overview of GRAND in its present and future incarnations, methods of radio-detection that are being developed for the…
▽ More
This is an index of the contributions by the Giant Radio Array for Neutrino Detection (GRAND) Collaboration to the 10th International Workshop on Acoustic and Radio EeV Neutrino Detection Activities (ARENA 2024, University of Chicago, June 11-14, 2024). The contributions include an overview of GRAND in its present and future incarnations, methods of radio-detection that are being developed for them, and ongoing joint work between the GRAND and BEACON experiments.
△ Less
Submitted 5 September, 2024;
originally announced September 2024.
-
Measuring the Time Variation of the Fine-structure Constant with Quasars Detected by LAMOST
Authors:
Jin-Nan Wei,
Rui-Jie Chen,
Jun-Jie Wei,
Martin Lopez-Corredoira,
Xue-Feng Wu
Abstract:
The [O III] $λλ 4960,5008$ emission lines in the optical spectra of galaxies and quasars have been widely used to investigate the possible variation of the fine-structure constant $α$ over cosmic time. In this work, we utilize the Large Sky Area Multi-object Fiber Spectroscopic Telescope (LAMOST) quasar survey, for the first time, to measure the relative $α$ variation $Δα/α$ in time through the [O…
▽ More
The [O III] $λλ 4960,5008$ emission lines in the optical spectra of galaxies and quasars have been widely used to investigate the possible variation of the fine-structure constant $α$ over cosmic time. In this work, we utilize the Large Sky Area Multi-object Fiber Spectroscopic Telescope (LAMOST) quasar survey, for the first time, to measure the relative $α$ variation $Δα/α$ in time through the [O III] doublet method. From the LAMOST Data Release 9 quasar catalog, we refine a sample of 209 quasar spectra with strong and narrow [O III] emission lines over a redshift range of $0<z<0.8$. Analysis on all of the 209 spectra obtains $Δα/α= (0.5 \pm 3.7) \times 10^{-4}$, which suggests that there is no evidence of varying $α$ on the explored cosmological timescales. Assuming a linear variation, the mean rate of change in $Δα/α$ is limited to be $(-3.4 \pm 2.4)\times 10^{-13}$ $\mathrm{yr^{-1}}$ in the last 7.0 Gyr. While our LAMOST-based constraint on $Δα/α$ is not competitive with those of the Sloan Digital Sky Survey (SDSS) quasar observations, our analysis serves to corroborate the results of SDSS with another independent survey.
△ Less
Submitted 2 September, 2024;
originally announced September 2024.
-
Energy evolution in the progenitor of galaxy shells: a semi-analytical model
Authors:
Beibei Guo,
Xufen Wu,
HongSheng Zhao,
Lulu Fan
Abstract:
The stellar shells surrounding an elliptical galaxy, as remnants of a dwarf galaxy disrupted during merging, reveal the distribution of energy and angular momentum of the progenitor dwarf galaxy. We develop a semi-analytical model to describe the changes of energy $ΔE_i$ and angular momentum $ΔLz_i$ for particles during the first infall. We show that these changes, induced by the self-gravity of t…
▽ More
The stellar shells surrounding an elliptical galaxy, as remnants of a dwarf galaxy disrupted during merging, reveal the distribution of energy and angular momentum of the progenitor dwarf galaxy. We develop a semi-analytical model to describe the changes of energy $ΔE_i$ and angular momentum $ΔLz_i$ for particles during the first infall. We show that these changes, induced by the self-gravity of the progenitor, are important in broadening the initial energy distribution of the Plummer or Hernquist progenitor model. Consequently, these changes are crucial in shaping the shells. In the free fall stage following the disintegration of the progenitor potential, particles are no longer bound by self-gravity but move within the gravitational potential of the target galaxy. We investigate the relationship between the radial period and the energy of particles undergoing radial motion. We show that an accurate model of the energy range of the dwarf galaxy at disruption is essential to predict the number of observable shells.
△ Less
Submitted 2 September, 2024;
originally announced September 2024.
-
Galaxies Lighting Up: Discovery of Seventy New Turn-on Changing-look Quasars
Authors:
Qian Yang,
Paul J. Green,
Xue-Bing Wu,
Michael Eracleous,
Linhua Jiang,
Yuming Fu
Abstract:
"Changing-look quasars" (CLQs), discovered less than a decade ago, show dramatic, rapid changes in optical/UV continuum and broad line emission. The majority of CLQs have been found dimming as "turn-off" CLQs because most selection methods start from samples of spectroscopically-confirmed quasars. We present here a sample of 82 spectroscopically confirmed "turn-on" CLQs, 70 of which are newly iden…
▽ More
"Changing-look quasars" (CLQs), discovered less than a decade ago, show dramatic, rapid changes in optical/UV continuum and broad line emission. The majority of CLQs have been found dimming as "turn-off" CLQs because most selection methods start from samples of spectroscopically-confirmed quasars. We present here a sample of 82 spectroscopically confirmed "turn-on" CLQs, 70 of which are newly identified. The turn-on CLQs are selected from spectroscopically classified galaxies with subsequent significant and dramatic variability in both the optical and mid-infrared bands, indicating a mechanism of changing accretion rate of the supermassive black holes rather than variable obscuration. Based on their bright state Eddington ratios, turn-on CLQs are associated with lower accretion rates compared to turn-off CLQs or typical SDSS quasars with similar redshift and magnitude distributions, even though turn-on CLQs have lower black hole masses. Most turn-on CLQs reside in host galaxies that follow local relations between the central black hole mass and host galaxy properties, such as stellar mass and velocity dispersion. However, their host galaxies have higher mass than normal inactive galaxies, with star formation rates more similar to hosts of Type 2 AGN than to the overall galaxy population.
△ Less
Submitted 28 August, 2024;
originally announced August 2024.
-
GRANDlib: A simulation pipeline for the Giant Radio Array for Neutrino Detection (GRAND)
Authors:
GRAND Collaboration,
Rafael Alves Batista,
Aurélien Benoit-Lévy,
Teresa Bister,
Martina Bohacova,
Mauricio Bustamante,
Washington Carvalho,
Yiren Chen,
LingMei Cheng,
Simon Chiche,
Jean-Marc Colley,
Pablo Correa,
Nicoleta Cucu Laurenciu,
Zigao Dai,
Rogerio M. de Almeida,
Beatriz de Errico,
Sijbrand de Jong,
João R. T. de Mello Neto,
Krijn D. de Vries,
Valentin Decoene,
Peter B. Denton,
Bohao Duan,
Kaikai Duan,
Ralph Engel,
William Erba
, et al. (90 additional authors not shown)
Abstract:
The operation of upcoming ultra-high-energy cosmic-ray, gamma-ray, and neutrino radio-detection experiments, like the Giant Radio Array for Neutrino Detection (GRAND), poses significant computational challenges involving the production of numerous simulations of particle showers and their detection, and a high data throughput. GRANDlib is an open-source software tool designed to meet these challen…
▽ More
The operation of upcoming ultra-high-energy cosmic-ray, gamma-ray, and neutrino radio-detection experiments, like the Giant Radio Array for Neutrino Detection (GRAND), poses significant computational challenges involving the production of numerous simulations of particle showers and their detection, and a high data throughput. GRANDlib is an open-source software tool designed to meet these challenges. Its primary goal is to perform end-to-end simulations of the detector operation, from the interaction of ultra-high-energy particles, through -- by interfacing with external air-shower simulations -- the ensuing particle shower development and its radio emission, to its detection by antenna arrays and its processing by data-acquisition systems. Additionally, GRANDlib manages the visualization, storage, and retrieval of experimental and simulated data. We present an overview of GRANDlib to serve as the basis of future GRAND analyses.
△ Less
Submitted 20 August, 2024;
originally announced August 2024.
-
Superfluid quantum criticality and the thermal evolution of neutron stars
Authors:
Hao-Fu Zhu,
Guo-Zhu Liu,
Jing-Rong Wang,
Xufen Wu
Abstract:
The neutron star starts to cool down shortly after its birth by emitting neutrinos. When it is cold enough, the Cooper pairs of neutrons are formed, which triggers a superfluid transition. Previous works on neutron superfluidity are focused on the finite-temperature transition. Little attention is paid to the potentially important quantum critical phenomena associated with superfluidity. Here, we…
▽ More
The neutron star starts to cool down shortly after its birth by emitting neutrinos. When it is cold enough, the Cooper pairs of neutrons are formed, which triggers a superfluid transition. Previous works on neutron superfluidity are focused on the finite-temperature transition. Little attention is paid to the potentially important quantum critical phenomena associated with superfluidity. Here, we provide the first theoretical analysis of superfluid quantum criticality, concentrating on its impact on neutron star cooling. Extensive calculations found that superfluidity occurs within a finite range of nuclear density $ρ$. The density serves as a non-thermal parameter for superfluid quantum phase transition. In a broad quantum critical region, there is a strong coupling between the gapless neutrons and the quantum critical fluctuations of superfluid order parameter. We handle this coupling by using perturbation theory and renormalization group method, and show that it leads to non-Fermi liquid behavior, which yields a logarithmic $T\ln(1/T)$ correction to the Fermi-liquid specific heat $c_{n}\propto T$ and also markedly alters the neutrino emissivity. Quantum critical phenomena emerge at a time much earlier than the onset of superfluidity and persist throughout almost the entire life of a neutron star. At a low $T$, such phenomena coexist with superfluidity in the neutron star interior, yet occupying different layers. We include superfluid quantum criticality into the theoretical description of neutron star cooling and find that it obviously prolongs the thermal relaxation time. By varying the strength of superfluid fluctuations and other quantities, we obtain an excellent fit to the observed cooling data of a variety of neutron stars. Our results indicate the existence of an intriguing correlation between superfluid quantum criticality and the thermal evolution of neutron stars.
△ Less
Submitted 7 August, 2024;
originally announced August 2024.
-
Long-term radio monitoring of the fast X-ray transient EP240315a: evidence for a relativistic jet
Authors:
R. Ricci,
E. Troja,
Y. Yang,
M. Yadav,
Y. Liu,
H. Sun,
X. Wu,
H. Gao,
B. Zhang,
W. Yuan
Abstract:
The recent launch of Einstein Probe (EP) in early 2024 opened up a new window onto the transient X-ray sky, allowing for real-time discovery and follow-up of fast X-ray transients (FXRTs). Multi-wavelength observations of FXRTs and their counterparts are key to characterize the properties of their outflows and, ultimately, identify their progenitors. Here, we report our long-term radio monitoring…
▽ More
The recent launch of Einstein Probe (EP) in early 2024 opened up a new window onto the transient X-ray sky, allowing for real-time discovery and follow-up of fast X-ray transients (FXRTs). Multi-wavelength observations of FXRTs and their counterparts are key to characterize the properties of their outflows and, ultimately, identify their progenitors. Here, we report our long-term radio monitoring of EP240315A, a long-lasting ($\sim 1000$ s) high redshift ($z=4.9$) FXRT associated to GRB~240315C. Our campaign, carried out with the Australian Telescope Compact Array (ATCA), followed the transient's evolution at two different frequencies (5.5 GHz and 9~GHz) for three months. In the radio lightcurves we identify an unusual steep rise at 9 GHz, possibly due to a refreshed reverse shock, and a late-time rapid decay of the radio flux, which we interpret as a jet break due to the outflow collimation. We find that the multi-wavelength counterpart of EP240315A is well described by a model of relativistic jet seen close to its axis, with jet half-opening angle $θ_j \approx 3 ^{\circ}$ and beaming-corrected total energy $E \simeq 4\times 10^{51}$~erg, typical of GRBs. These results show that a substantial fraction of FXRTs may be associated to standard GRBs and that sensitive X-ray monitors, such as Einstein Probe and the proposed HiZ-GUNDAM and Theseus missions, can successfully pinpoint their relativistic outflows up to high-redshifts.
△ Less
Submitted 25 July, 2024;
originally announced July 2024.
-
Tracing the evolution of the cool gas in CGM and IGM environments through Mg II absorption from redshift z=0.75 to z=1.65 using DESI-Y1 data
Authors:
X. Wu,
Z. Cai,
T. -W. Lan,
S. Zou,
A. Anand,
Biprateep Dey,
Z. Li,
J. Aguilar,
S. Ahlen,
D. Brooks,
T. Claybaugh,
A. de la Macorra,
P. Doel,
S. Ferraro,
J. E. Forero-Romero,
S. Gontcho A Gontcho,
K. Honscheid,
S. Juneau,
R. Kehoe,
T. Kisner,
A. Lambert,
M. Landriau,
L. Le Guillou,
M. Manera,
A. Meisner
, et al. (13 additional authors not shown)
Abstract:
We present a measurement of the mean absorption of cool gas traced by Mg II (${λλ2796, 2803}$) around emission line galaxies (ELGs), spanning spatial scales from 20 kpc to 10 Mpc. The measurement is based on cross-matching the positions of about 2.5 million ELGs at $z = 0.75-1.65$ and the metal absorption in the spectra of 1.4 million background quasars with data provided by the Year 1 sample of t…
▽ More
We present a measurement of the mean absorption of cool gas traced by Mg II (${λλ2796, 2803}$) around emission line galaxies (ELGs), spanning spatial scales from 20 kpc to 10 Mpc. The measurement is based on cross-matching the positions of about 2.5 million ELGs at $z = 0.75-1.65$ and the metal absorption in the spectra of 1.4 million background quasars with data provided by the Year 1 sample of the Dark Energy Spectroscopic Instrument (DESI). The ELGs are divided into two redshift intervals: $0.75 < z < 1.0$ and $1.0 < z < 1.65$. We find that the composite spectra constructed by stacking the ELG-QSO pairs show evolution with redshift, with $z>1$ having a systematically higher signal of Mg II absorption. Within 1 Mpc, the covering fraction of the cool gas at $z > 1$ is higher than that of $z < 1$. The enhancement becomes less apparent especially if the projected distance $r_{p}>$1 Mpc. Also, ELGs with higher stellar mass and star formation rate (SFR) yield higher clustering of Mg II absorbers at $z<1$. For $z>1$, the covering fractions with different SFRs show little difference. The higher Mg II absorption at higher redshift also supports the observations of higher star formation at cosmic noon. Besides, the profile of Mg II absorption reveals a change of slope on scales of about 1 Mpc, consistent with the expected transition from a dark matter halo-dominated environment to a regime where clustering is dominated by halo-halo correlations. We estimate the cool gas density profile and derive the metal abundance at different redshifts. The growth of metal abundance suggests an increased presence of cool gas in the intergalactic medium (IGM) towards higher redshifts.
△ Less
Submitted 25 July, 2024;
originally announced July 2024.
-
Triggering the Untriggered: The First Einstein Probe-Detected Gamma-Ray Burst 240219A and Its Implications
Authors:
Yi-Han Iris Yin,
Bin-Bin Zhang,
Jun Yang,
Hui Sun,
Chen Zhang,
Yi-Xuan Shao,
You-Dong Hu,
Zi-Pei Zhu,
Dong Xu,
Li An,
He Gao,
Xue-Feng Wu,
Bing Zhang,
Alberto Javier Castro-Tirado,
Shashi B. Pandey,
Arne Rau,
Weihua Lei,
Wei Xie,
Giancarlo Ghirlanda,
Luigi Piro,
Paul O'Brien,
Eleonora Troja,
Peter Jonker,
Yun-Wei Yu,
Jie An
, et al. (26 additional authors not shown)
Abstract:
The Einstein Probe (EP) achieved its first detection and localization of a bright X-ray flare, EP240219a, on February 19, 2024, during its commissioning phase. Subsequent targeted searches triggered by the EP240219a alert identified a faint, untriggered gamma-ray burst (GRB) in the archived data of Fermi/GBM, Swift/BAT, Insight-HXMT/HE and INTEGRAL/SPI-ACS. The EP/WXT light curve reveals a long du…
▽ More
The Einstein Probe (EP) achieved its first detection and localization of a bright X-ray flare, EP240219a, on February 19, 2024, during its commissioning phase. Subsequent targeted searches triggered by the EP240219a alert identified a faint, untriggered gamma-ray burst (GRB) in the archived data of Fermi/GBM, Swift/BAT, Insight-HXMT/HE and INTEGRAL/SPI-ACS. The EP/WXT light curve reveals a long duration of approximately 160 seconds with a slow decay, whereas the Fermi/GBM light curve shows a total duration of approximately 70 seconds. The peak in the Fermi/GBM light curve occurs slightly later with respect to the peak seen in the EP/WXT light curve. Our spectral analysis shows that a single cutoff power-law model effectively describes the joint EP/WXT-Fermi/GBM spectra in general, indicating coherent broad emission typical of GRBs. The model yielded a photon index of $\sim -1.70 \pm 0.05$ and a peak energy of $\sim 257 \pm 134$ keV. After detection of GRB 240219A, long-term observations identified several candidates in optical and radio wavelengths, none of which was confirmed as the afterglow counterpart during subsequent optical and near-infrared follow-ups. The analysis of GRB 240219A classifies it as an X-ray rich GRB with a high peak energy, presenting both challenges and opportunities for studying the physical origins of X-ray flashes (XRFs), X-ray rich GRBs (XRRs), and classical GRBs (C-GRBs). Furthermore, linking the cutoff power-law component to non-thermal synchrotron radiation suggests that the burst is driven by a Poynting flux-dominated outflow.
△ Less
Submitted 14 July, 2024;
originally announced July 2024.
-
Detection of hidden emissions in two rotating radio transients with high surface magnetic fields
Authors:
S. B. Zhang,
X. Yang,
J. J. Geng,
Y. P. Yang,
X. F. Wu
Abstract:
Rotating Radio Transients (RRATs) are neutron stars emitting sporadic radio pulses. The unique emission of RRATs has been proposed to resemble those of known pulsar types, such as extreme nulling pulsars or pulsars with giant pulses. However, the presence of additional radiation beyond these sporadic pulses remains unclear. Through high-sensitivity observations and extended tracking, we detected t…
▽ More
Rotating Radio Transients (RRATs) are neutron stars emitting sporadic radio pulses. The unique emission of RRATs has been proposed to resemble those of known pulsar types, such as extreme nulling pulsars or pulsars with giant pulses. However, the presence of additional radiation beyond these sporadic pulses remains unclear. Through high-sensitivity observations and extended tracking, we detected the sequential weak emissions in two RRATs with relatively high surface magnetic fields (Bs > 10^13 G): J1846-0257 and J1854+0306. These emissions show peak flux densities of 0.15 and 0.41 mJy, up to 687 and 512 times weaker than our detected RRAT single pulses, respectively. The weak emissions contribute small fractions (~ 16% and 5%) to the total radio pulse energy releases, contrasting significantly with giant-pulse pulsars where normal pulses dominate. Polarization analysis of J1854+0306 suggests that its sporadic RRAT pulses may originate from intermittent enhanced sparking processes due to magnetospheric evolution. Our findings indicate that some RRATs may represent a novel class of pulsars, distinct from any previously known subclass. Further observations of sources with similar rotational properties using high-sensitivity instruments could validate the generality of these hidden emissions.
△ Less
Submitted 13 July, 2024;
originally announced July 2024.
-
A lensed FRB candidate in the first CHIME/FRB Catalogue and its potential implications
Authors:
Chenming Chang,
Songbo Zhang,
Di Xiao,
Zhenfan Tang,
Ye Li,
Junjie Wei,
Xuefeng Wu
Abstract:
Fast radio bursts (FRBs) are immensely energetic radio pulses with durations of milliseconds. Given their high all-sky rate, the probability of an FRB being lensed by an intervening massive object is non-negligible. In this study, we search for possible lensing candidates within the first Canadian Hydrogen Intensity Mapping Experiment FRB catalogue using an autocorrelation algorithm and verificati…
▽ More
Fast radio bursts (FRBs) are immensely energetic radio pulses with durations of milliseconds. Given their high all-sky rate, the probability of an FRB being lensed by an intervening massive object is non-negligible. In this study, we search for possible lensing candidates within the first Canadian Hydrogen Intensity Mapping Experiment FRB catalogue using an autocorrelation algorithm and verification through signal simulations. We identify FRB 20190308C as a lensed candidate with a significance of 3.4 sigma. Furthermore, we constrain the mass of the lensing object using the Chang-Refsdal lens model, based on the flux ratio and time delay between the substructures of FRB 20190308C. Future long-term and high-precision observations are expected to reveal more lensed FRBs.
△ Less
Submitted 28 June, 2024;
originally announced June 2024.
-
Effects of Two Quantum Correction Parameters on Chaotic Dynamics of Particles Near Renormalized Group Improved Schwarzschild Black Holes
Authors:
Junjie Lu,
Xin Wu
Abstract:
A renormalized group improved Schwarzschild black hole spacetime contains two quantum correction parameters. One parameter $γ$ represents the identification of cutoff of the distance scale, and another parameter $Ω$ stems from nonperturbative renormalization group theory. The two parameters are constrained by the data from the shadow of M87* central black hole. The dynamics of electrically charged…
▽ More
A renormalized group improved Schwarzschild black hole spacetime contains two quantum correction parameters. One parameter $γ$ represents the identification of cutoff of the distance scale, and another parameter $Ω$ stems from nonperturbative renormalization group theory. The two parameters are constrained by the data from the shadow of M87* central black hole. The dynamics of electrically charged test particles around the black hole are integrable. However, when the black hole is immersed in an external asymptotically uniform magnetic field, the dynamics are not integrable and may allow for the occurrence of chaos. % English Editor: %Please verify that the intended meaning has been retained. Employing an explicit symplectic integrator, we survey the contributions of the two parameters to the chaotic dynamical behavior. It is found that a small change of the parameter $γ$ constrained by the shadow of M87* black hole has an almost negligible effect on the dynamical transition of particles from order to chaos. However, a small decrease in the parameter $Ω$ leads to an enhancement in the strength of chaos from the global phase space structure. A theoretical interpretation is given to the different contributions. The term with the parameter $Ω$ dominates the term with the parameter $γ$, even if the two parameters have same values. % English Editor: Please verify %that the intended meaning has been retained. In particular, the parameter $Ω$ acts as a repulsive force, and its decrease means a weakening of the repulsive force or equivalently enhancing the attractive force from the black hole. On the other hand, there is a positive Lyapunov exponent that is universally given by the surface gravity of the black hole when $Ω\geq 0$ is small and the external magnetic field vanishes.
△ Less
Submitted 27 June, 2024;
originally announced June 2024.
-
Searching anomalies using nonlinear dimensionality reduction techniques
Authors:
X. Yang,
G. Hobbs,
S. -B. Zhang,
A. Zic,
Lawrence Toomey,
Y. Li,
J. -S. Wang,
S. Dai,
X. -F. Wu
Abstract:
We have searched for anomalous events using 2,520 hours of archival observations from Murriyang, CSIRO's Parkes radio telescope. These observations were originally undertaken to search for pulsars. We used a machine-learning algorithm based on ResNet and Uniform Manifold Approximation and Projection (UMAP) in order to identify parts of the data stream that potentially contain anomalous signals. Ma…
▽ More
We have searched for anomalous events using 2,520 hours of archival observations from Murriyang, CSIRO's Parkes radio telescope. These observations were originally undertaken to search for pulsars. We used a machine-learning algorithm based on ResNet and Uniform Manifold Approximation and Projection (UMAP) in order to identify parts of the data stream that potentially contain anomalous signals. Many of these anomalous events are radio frequency interference, which were subsequently filtered using multibeam information. We detected 202 anomalous events and provide their positions and event times. We discuss the possibility that one of the events comes from radio emission from a white dwarf star.The other events are currently of unknown type.
△ Less
Submitted 1 July, 2024; v1 submitted 26 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.
-
Response of the first POLAR-2 Prototype to Polarized Beams
Authors:
Merlin Kole,
Nicolas de Angelis,
Ana Bacelj,
Franck Cadoux,
Agnieszka Elwertowska,
Johannes Hulsman,
Hancheng Li,
Grzegorz Łubian,
Tomasz Kowalski,
Gilles Koziol,
Agnieszka Pollo,
Nicolas Produit,
Dominik Rybka,
Adrien Stil,
Jianchao Sun,
Xin Wu,
Kacper Zezuliński,
Shuang-Nan Zhang
Abstract:
POLAR-2 is a dedicated gamma-ray polarimeter currently foreseen to be launched towards the China Space Station around 2027. The design of the detector is based on the legacy of its predecessor mission POLAR which was launched in 2016. POLAR-2 aims to measure the polarization of the Gamma-ray Burst prompt emission within the 30-800 keV energy range. Thanks to its high sensitivity to gamma-ray polar…
▽ More
POLAR-2 is a dedicated gamma-ray polarimeter currently foreseen to be launched towards the China Space Station around 2027. The design of the detector is based on the legacy of its predecessor mission POLAR which was launched in 2016. POLAR-2 aims to measure the polarization of the Gamma-ray Burst prompt emission within the 30-800 keV energy range. Thanks to its high sensitivity to gamma-ray polarization, as well as its large effective area, POLAR-2 will provide the most precise measurements of this type to date. Such measurements are key to improve our understanding of the astrophysical processes responsible for Gamma-Ray Bursts. The detector consists of a segmented array of plastic scintillator bars, each one of which is read out by a Silicon PhotoMultiplier channel. The flight model of POLAR-2 will contain a total of 6400 scintillators. These are divided into 100 groups of 64 bars each, in so-called polarimeter modules. In recent years, the collaboration has designed and produced the first prototypes of these polarimeter modules and subjected these to space qualification tests. In addition, in April 2023, the first of these modules were calibrated using fully polarized gamma-ray beams at the European Synchrotron Radiation Facility (ESRF) in France. In this work, we will present the results of this calibration campaign and compare these to the simulated performance of the POLAR-2 modules. Potential improvements to the design are also discussed. Finally, the measurements are used, in combination with the verified simulation framework, to estimate the scientific performance of the full POLAR-2 detector and compare it to its predecessor.
△ Less
Submitted 9 June, 2024;
originally announced June 2024.
-
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.
-
The Intermediate-Mass Black Hole Reverberation Mapping Project: Initial Results for a candidate IMBH in a nearby Seyfert 1 Galaxy
Authors:
Wenwen Zuo,
Hengxiao Guo,
Jingbo Sun,
Qi Yuan,
Paulina Lira,
Minfeng Gu,
Philip G. Edwards,
Alok C. Gupta,
Shubham Kishore,
Jamie Stevens,
Tao An,
Zhen-Yi Cai,
Haicheng Feng,
Luis C. Ho,
Dragana Ilić,
Andjelka B. Kovačević,
ShaSha Li,
Mar Mezcua,
Luka Č. Popović,
Mouyuan Sun,
Tushar Tripathi,
Vivian U.,
Oliver Vince,
Jianguo Wang,
Junxian Wang
, et al. (3 additional authors not shown)
Abstract:
To investigate the short-term variability and determine the size of the optical continuum emitting size of intermediate-mass black holes (IMBHs), we carried out high-cadence, multi-band photometric monitoring of a Seyfert 1 galaxy J0249-0815 across two nights, together with a one-night single-band preliminary test. The presence of the broad Ha component in our target was confirmed by recent Paloma…
▽ More
To investigate the short-term variability and determine the size of the optical continuum emitting size of intermediate-mass black holes (IMBHs), we carried out high-cadence, multi-band photometric monitoring of a Seyfert 1 galaxy J0249-0815 across two nights, together with a one-night single-band preliminary test. The presence of the broad Ha component in our target was confirmed by recent Palomar/P200 spectroscopic observations, 23 years after Sloan Digital Sky Survey, ruling out the supernovae origin of the broad Ha line. The photometric experiment was primarily conducted utilizing four-channel imagers MuSCAT 3 & 4 mounted on 2-meter telescopes within the Las Cumbres Observatory Global Telescope Network. Despite the expectation of variability, we observed no significant variation (<1.4%) on timescales of 6-10 hours. This non-detection is likely due to substantial host galaxy light diluting the subtle AGN variability. Dual-band preliminary tests and tailored simulations may enhance the possibility of detecting variability and lag in future IMBH reverberation campaigns.
△ Less
Submitted 19 May, 2024;
originally announced May 2024.
-
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.
-
Large Scale Overdensity of Lyman Break Galaxies Around the z=6.3 Ultraluminous Quasar J0100+2802
Authors:
Maria Pudoka,
Feige Wang,
Xiaohui Fan,
Jinyi Yang,
Jaclyn Champagne,
Victoria Jones,
Fuyan Bian,
Zheng Cai,
Linhua Jiang,
Dezi Liu,
Xue-Bing Wu
Abstract:
We study the environment of the z=6.33 ultraluminous quasar SDSS J010013.02+280225.8 (J0100) to understand its association with large-scale structure. Theoretical models propose high-redshift quasars as markers of galaxy overdensities residing in the most massive dark matter halos (DMHs) in the early universe. J0100 is an ultraluminous quasar with the most massive black hole known at z>6, suggesti…
▽ More
We study the environment of the z=6.33 ultraluminous quasar SDSS J010013.02+280225.8 (J0100) to understand its association with large-scale structure. Theoretical models propose high-redshift quasars as markers of galaxy overdensities residing in the most massive dark matter halos (DMHs) in the early universe. J0100 is an ultraluminous quasar with the most massive black hole known at z>6, suggesting a high likelihood of residing in a massive DMH. We present wide-field ($\sim$522 square arcminute) imaging in the r-, i-, and z-bands from the Large Binocular Camera on the Large Binocular Telescope, with Y- and J-band imaging from the Wide-field Infrared Camera on the Canada-France-Hawaii Telescope, centered on J0100. Applying color selections, we identify 23 objects as i-droput Lyman Break Galaxy (LBG) candidates in the J0100 field. We use the deep photometric catalog in the 1.27 square degree COSMOS field to calculate the density of LBGs in a blank field, and to estimate the selection completeness and purity. The observed surface density of LBG candidates in the J0100 field corresponds to a galaxy overdensity of $δ$=4 (at 8.4$σ$). This large-scale overdensity suggests that the $\sim$ 22 square arcminute overdensity found by Kashino et al. using JWST data extends out to much larger scales. We calculate the angular auto-correlation function of the candidates and find a positive correlation on $\lesssim$ 10 arcminute scales as well as evidence of asymmetries in their spatial distribution, further suggesting a direct detection of large-scale structure around the ultra-luminous quasar J0100.
△ Less
Submitted 6 May, 2024;
originally announced May 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.
-
A Spatially Resolved [CII] Survey of 31 $z\sim7$ Massive Galaxies Hosting Luminous Quasars
Authors:
Feige Wang,
Jinyi Yang,
Xiaohui Fan,
Bram Venemans,
Roberto Decarli,
Eduardo Bañados,
Fabian Walter,
Aaron J. Barth,
Fuyan Bian,
Frederick B. Davies,
Anna-Christina Eilers,
Emanuele Paolo Farina,
Joseph F. Hennawi,
Jiang-Tao Li,
Chiara Mazzucchelli,
Ran Wang,
Xue-Bing Wu,
Minghao Yue
Abstract:
The [CII] 158 $μ$m emission line and the underlying far-infrared (FIR) dust continuum are important tracers for studying star formation and kinematic properties of early galaxies. We present a survey of the [CII] emission lines and FIR continua of 31 luminous quasars at $z>6.5$ using the Atacama Large Millimeter Array (ALMA) and the NOrthern Extended Millimeter Array (NOEMA) at sub-arcsec resoluti…
▽ More
The [CII] 158 $μ$m emission line and the underlying far-infrared (FIR) dust continuum are important tracers for studying star formation and kinematic properties of early galaxies. We present a survey of the [CII] emission lines and FIR continua of 31 luminous quasars at $z>6.5$ using the Atacama Large Millimeter Array (ALMA) and the NOrthern Extended Millimeter Array (NOEMA) at sub-arcsec resolution. This survey more than doubles the number of quasars with [CII] and FIR observations at these redshifts and enables statistical studies of quasar host galaxies deep into the epoch of reionization. We detect [CII] emission in 27 quasar hosts with a luminosity range of $L_{\rm [CII]}=(0.3-5.5)\times10^9~L_\odot$ and detect the FIR continuum of 28 quasar hosts with a luminosity range of $L_{\rm FIR}=(0.5-13.0)\times10^{12}~L_\odot$. Both $L_{\rm [CII]}$ and $L_{\rm FIR}$ are correlated ($ρ\simeq0.4$) with the quasar bolometric luminosity, albeit with substantial scatter. The quasar hosts detected by ALMA are clearly resolved with a median diameter of $\sim$5 kpc. About 40% of the quasar host galaxies show a velocity gradient in [CII] emission, while the rest show either dispersion-dominated or disturbed kinematics. Basic estimates of the dynamical masses of the rotation-dominated host galaxies yield $M_{\rm dyn}=(0.1-7.5)\times10^{11}~M_\odot$. Considering our findings alongside those of literature studies, we found that the ratio between $M_{\rm BH}$ and $M_{\rm dyn}$ is about ten times higher than that of local $M_{\rm BH}-M_{\rm dyn}$ relation on average but with substantial scatter (the ratio difference ranging from $\sim$0.6 to 60) and large uncertainties.
△ Less
Submitted 23 April, 2024;
originally announced April 2024.
-
Efficient Identification of Broad Absorption Line Quasars using Dimensionality Reduction and Machine Learning
Authors:
Wei-Bo Kao,
Yanxia Zhang,
Xue-Bing Wu
Abstract:
Broad Absorption Line Quasars (BALQSOs) displaying distinct blue-shifted broad absorption lines. These serve as invaluable probes for unraveling the intricate structure and evolution of quasars, shedding light on the profound influence exerted by supermassive black holes on galaxy formation. The proliferation of large-scale spectroscopic surveys such as LAMOST, SDSS, and DESI has exponentially exp…
▽ More
Broad Absorption Line Quasars (BALQSOs) displaying distinct blue-shifted broad absorption lines. These serve as invaluable probes for unraveling the intricate structure and evolution of quasars, shedding light on the profound influence exerted by supermassive black holes on galaxy formation. The proliferation of large-scale spectroscopic surveys such as LAMOST, SDSS, and DESI has exponentially expanded the repository of quasar spectra at our disposal. In this study, we present an innovative approach to streamline the identification of BALQSOs, leveraging the power of dimensionality reduction and machine learning algorithms. Our dataset is curated from the SDSS DR16, amalgamating quasar spectra with classification labels sourced from the DR16Q quasar catalog. We employ a diverse array of dimensionality reduction techniques, including Principal Component Analysis (PCA), t-Distributed Stochastic Neighbor Embedding (t-SNE), Locally Linear Embedding (LLE), and Isometric Mapping (ISOMAP), to distill the essence of the original spectral data. The resultant low-dimensional representations serve as inputs for a suite of machine learning classifiers, including XGBoost and Random Forest models. Through experimentation, we unveil PCA as the most effective dimensionality reduction methodology, adeptly navigating the intricate balance between dimensionality reduction and preservation of vital spectral information. Notably, the synergistic fusion of PCA with the XGBoost classifier emerges as the pinnacle of efficacy in the BALQSO classification endeavor, boasting impressive accuracy rates of 97.60% by 10-cross validation and 96.92% on the outer test sample. This study not only introduces a novel machine learning-based paradigm for quasar classification but also offers invaluable insights transferrable to a myriad of spectral classification challenges pervasive in the realm of astronomy.
△ Less
Submitted 18 April, 2024;
originally announced April 2024.
-
H$α$ Reverberation Mapping from Broad-band Photometry of Dwarf Seyfert 1 Galaxy NGC 4395
Authors:
Huapeng Gu,
Xue-Bing Wu,
Yuhan Wen,
Qinchun Ma,
Hengxiao Guo
Abstract:
NGC 4395 is a dwarf Seyfert 1 galaxy with a possible intermediate-mass black hole of several $\rm{10^4}$ solar masses in its center. As a well-studied object, its broad line region size has been measured via H$\rmα$ time lag in numerous spectroscopic reverberation mapping (SRM) and narrow-band photometric reverberation mapping (PRM) campaigns. Here we present its H$\rmα$ time lag measurement using…
▽ More
NGC 4395 is a dwarf Seyfert 1 galaxy with a possible intermediate-mass black hole of several $\rm{10^4}$ solar masses in its center. As a well-studied object, its broad line region size has been measured via H$\rmα$ time lag in numerous spectroscopic reverberation mapping (SRM) and narrow-band photometric reverberation mapping (PRM) campaigns. Here we present its H$\rmα$ time lag measurement using broad-band photometric data, with the application of our newly-developed ICCF-Cut method as well as the JAVELIN and $χ^2$ methods. utilizing the minute-cadence multi-band light curves obtained from the $\rm{2}$m FTN and $\rm{10.4}$m GTC telescopes in recent works, we measured its H$\rmα$ lag as approximately $40 \sim 90$ minutes from broad-band PRM. With the H$\rmα$ emission line velocity dispersion, we calculated its central black hole mass as $\rm M_{\rm BH} = (8\pm4) \times 10^3\, M_{\rm \odot}$. These results are comparable with previous results obtained by narrow-band PRM and SRM, providing further support to an intermediate-mass black hole in NGC 4395. In addition, our study also validates the ICCF-Cut as an effective method for broad-band PRM, which holds the potential for widespread application in the era of large multi-epoch, high-cadence photometric surveys.
△ Less
Submitted 17 April, 2024; v1 submitted 17 April, 2024;
originally announced April 2024.
-
LHAASO-KM2A detector simulation using Geant4
Authors:
Zhen Cao,
F. Aharonian,
Q. An,
Axikegu,
Y. X. Bai,
Y. W. Bao,
D. Bastieri,
X. J. Bi,
Y. J. Bi,
J. T. Cai,
Q. Cao,
W. Y. Cao,
Zhe Cao,
J. Chang,
J. F. Chang,
A. M. Chen,
E. S. Chen,
Liang Chen,
Lin Chen,
Long Chen,
M. J. Chen,
M. L. Chen,
Q. H. Chen,
S. H. Chen,
S. Z. Chen
, et al. (254 additional authors not shown)
Abstract:
KM2A is one of the main sub-arrays of LHAASO, working on gamma ray astronomy and cosmic ray physics at energies above 10 TeV. Detector simulation is the important foundation for estimating detector performance and data analysis. It is a big challenge to simulate the KM2A detector in the framework of Geant4 due to the need to track numerous photons from a large number of detector units (>6000) with…
▽ More
KM2A is one of the main sub-arrays of LHAASO, working on gamma ray astronomy and cosmic ray physics at energies above 10 TeV. Detector simulation is the important foundation for estimating detector performance and data analysis. It is a big challenge to simulate the KM2A detector in the framework of Geant4 due to the need to track numerous photons from a large number of detector units (>6000) with large altitude difference (30 m) and huge coverage (1.3 km^2). In this paper, the design of the KM2A simulation code G4KM2A based on Geant4 is introduced. The process of G4KM2A is optimized mainly in memory consumption to avoid memory overffow. Some simpliffcations are used to signiffcantly speed up the execution of G4KM2A. The running time is reduced by at least 30 times compared to full detector simulation. The particle distributions and the core/angle resolution comparison between simulation and experimental data of the full KM2A array are also presented, which show good agreement.
△ Less
Submitted 7 April, 2024;
originally announced April 2024.
-
Evidence for an Outer Component in the Continuum Reverberation Mapping of Active Galactic Nuclei
Authors:
Yuanzhe Jiang,
Xue-Bing Wu,
Qinchun Ma,
Huapeng Gu,
Yuhan Wen
Abstract:
The continuum reverberation mapping is widely used in studying accretion disk of active galactic nuclei (AGN). While some indirect evidence and simulations indicated that the diffuse continuum, especially the strong Balmer continuum from the broad line region (BLR), may contribute to the continuum in the u/U band. Here, we present direct evidence for this contribution. In this work, we apply the I…
▽ More
The continuum reverberation mapping is widely used in studying accretion disk of active galactic nuclei (AGN). While some indirect evidence and simulations indicated that the diffuse continuum, especially the strong Balmer continuum from the broad line region (BLR), may contribute to the continuum in the u/U band. Here, we present direct evidence for this contribution. In this work, we apply the ICCF-Cut method to continuum reverberation mapping to extract the possible diffuse continuum light curves of 6 AGNs with high cadence, high quality and multi-band observations. We find the existence of an outer component out of the accretion disk for each of 6 AGNs in the Swift U band. Meanwhile, similar results can be derived by JAVELIN Photometric Reverberation Mapping Model for 4 of them. The lags of the outer components are consistent with the predicted Balmer continuum lags, which are about half of the H$β$ lag values. Our result directly reinforces that an outer component, especially the Balmer continuum in the rest-frame u/U band, can contribute significantly to the continuum reverberation lags of AGNs.
△ Less
Submitted 20 March, 2024;
originally announced March 2024.
-
H$α$ Time Delays of AGNs from the Zwicky Transcient Facility Broadband Photometry
Authors:
Qinchun Ma,
Yuhan Wen,
Xue-Bing Wu,
Huapeng Gu,
Yuming Fu
Abstract:
In our previous work on broadband photometric reverberation mapping (PRM), we proposed the ICCF-Cut process to obtain the time lags of H$α$ emission line from two broadband lightcurves via subtracting the continuum emission from the line band. Extending the work, we enlarge our sample to the Zwicky Transient Facility (ZTF) database. We adopt two criteria to select 123 type 1 AGNs with sufficient v…
▽ More
In our previous work on broadband photometric reverberation mapping (PRM), we proposed the ICCF-Cut process to obtain the time lags of H$α$ emission line from two broadband lightcurves via subtracting the continuum emission from the line band. Extending the work, we enlarge our sample to the Zwicky Transient Facility (ZTF) database. We adopt two criteria to select 123 type 1 AGNs with sufficient variability and smooth lightcurves from 3537 AGNs at $z<0.09$ with more than 100 epoch observations in the $g$ and $r$ bands from the ZTF database. We calculate the H$α$ time lags for 23 of them which have previous spectroscopic reverberation mapping (SRM) results using ICCF-Cut, JAVELIN and $χ^2$ methods. Our obtained H$α$ time lags are slightly larger than the H$β$ time lags, which is consistent with the previous SRM results and the theoretical model of the AGN broad line region. The comparisons between SRM and PRM lag distributions and between the subtracted emission line lightcurves indicate that after selecting AGNs with the two criteria, combining the ICCF-Cut, JAVELIN and $χ^2$ methods provides an efficient way to get the reliable H$α$ lags from the broadband PRM. Such techniques can be used to estimate the black hole masses of a large sample of AGNs in the large multi-epoch photometric sky surveys such as the Legacy Survey of Space and Time (LSST) and the survey from the Wide Field Survey Telescope (WFST) in the near future.
△ Less
Submitted 15 March, 2024;
originally announced March 2024.
-
Measurements of All-Particle Energy Spectrum and Mean Logarithmic Mass of Cosmic Rays from 0.3 to 30 PeV with LHAASO-KM2A
Authors:
The LHAASO Collaboration,
Zhen Cao,
F. Aharonian,
Q. An,
A. Axikegu,
Y. X. Bai,
Y. W. Bao,
D. Bastieri,
X. J. Bi,
Y. J. Bi,
J. T. Cai,
Q. Cao,
W. Y. Cao,
Zhe Cao,
J. Chang,
J. F. Chang,
A. M. Chen,
E. S. Chen,
Liang Chen,
Lin Chen,
Long Chen,
M. J. Chen,
M. L. Chen,
Q. H. Chen,
S. H. Chen
, et al. (256 additional authors not shown)
Abstract:
We present the measurements of all-particle energy spectrum and mean logarithmic mass of cosmic rays in the energy range of 0.3-30 PeV using data collected from LHAASO-KM2A between September 2021 and December 2022, which is based on a nearly composition-independent energy reconstruction method, achieving unprecedented accuracy. Our analysis reveals the position of the knee at…
▽ More
We present the measurements of all-particle energy spectrum and mean logarithmic mass of cosmic rays in the energy range of 0.3-30 PeV using data collected from LHAASO-KM2A between September 2021 and December 2022, which is based on a nearly composition-independent energy reconstruction method, achieving unprecedented accuracy. Our analysis reveals the position of the knee at $3.67 \pm 0.05 \pm 0.15$ PeV. Below the knee, the spectral index is found to be -$2.7413 \pm 0.0004 \pm 0.0050$, while above the knee, it is -$3.128 \pm 0.005 \pm 0.027$, with the sharpness of the transition measured with a statistical error of 2%. The mean logarithmic mass of cosmic rays is almost heavier than helium in the whole measured energy range. It decreases from 1.7 at 0.3 PeV to 1.3 at 3 PeV, representing a 24% decline following a power law with an index of -$0.1200 \pm 0.0003 \pm 0.0341$. This is equivalent to an increase in abundance of light components. Above the knee, the mean logarithmic mass exhibits a power law trend towards heavier components, which is reversal to the behavior observed in the all-particle energy spectrum. Additionally, the knee position and the change in power-law index are approximately the same. These findings suggest that the knee observed in the all-particle spectrum corresponds to the knee of the light component, rather than the medium-heavy components.
△ Less
Submitted 26 March, 2024; v1 submitted 15 March, 2024;
originally announced March 2024.
-
Partial Tidal Disruption Events by Intermediate-mass Black Holes in Supermassive and Intermediate-mass Black Hole Binaries
Authors:
Xiao-Jun Wu,
Ye-Fei Yuan,
Yan Luo,
Wenbin Lin
Abstract:
In the centers of galaxies, stars that orbit supermassive black hole binaries (SMBHBs) can undergo tidal disruptions due to the Lidov-Kozai mechanism. Nevertheless, most previous researches have predominantly focused on full tidal disruption events (FTDEs). In this study, we employ N-body simulations to investigate partial tidal disruption events (PTDEs) induced by intermediate-mass black holes (I…
▽ More
In the centers of galaxies, stars that orbit supermassive black hole binaries (SMBHBs) can undergo tidal disruptions due to the Lidov-Kozai mechanism. Nevertheless, most previous researches have predominantly focused on full tidal disruption events (FTDEs). In this study, we employ N-body simulations to investigate partial tidal disruption events (PTDEs) induced by intermediate-mass black holes (IMBHs) in SMBH-IMBH binaries, taking into account consideration the IMBH's mass, semi-major axis, and eccentricity of the outer orbit. Our findings indicate that, in comparison to FTDEs, the majority of tidal disruption events are actually PTDEs. Furthermore, we find that a significant number of stars experiencing partial disruption ultimately get captured by the IMBH, potentially leading to repeating flares. By comparing the period of the periodic eruptions observed in ASASSN-14ko, we find that PTDEs in a specific SMBH-IMBH binary system can align with the observed period if the SMBH has a mass of $10^7\rm{\ M_\odot}$, the IMBH has a mass smaller than approximately $10^5\rm{\ M_\odot}$, the eccentricity of the SMBH-IMBH binary exceeds approximately $0.5$, and the semi-major axis of the SMBH-IMBH binary is larger than approximately $0.001\rm{\ pc}$. Moreover, our model effectively accounts for the observed period derivative for ASASSN-14ko ($\dot{P}=-0.0026\pm 0.0006$), and our results also imply that some quasi-periodic eruptions may be attributed to PTDEs occurring around SMBH-IMBH binaries.
△ Less
Submitted 28 February, 2024;
originally announced February 2024.
-
Quasar winds caught on acceleration and deceleration
Authors:
Weimin Yi,
P. B. Hall,
Zunlin Yuan,
W. N. Brandt,
D. P. Schneider,
Zhicheng He,
Jin-Ming Bai,
Xue-Bing Wu
Abstract:
We present an observational study of wind acceleration based on four low-ionization broad absorption line (LoBAL) quasars (J0136, J1238, J1259, J1344). J0136 and J1344 (group-1) are radio quiet and show large BAL-velocity shifts as opposed to stable line-locking associated absorption lines (AALs). Notably, J1344 displays a linear relation between BAL-velocity shift and time interval over three con…
▽ More
We present an observational study of wind acceleration based on four low-ionization broad absorption line (LoBAL) quasars (J0136, J1238, J1259, J1344). J0136 and J1344 (group-1) are radio quiet and show large BAL-velocity shifts as opposed to stable line-locking associated absorption lines (AALs). Notably, J1344 displays a linear relation between BAL-velocity shift and time interval over three consecutive epochs, characteristic of compelling evidence for BAL acceleration. J1238 and J1259 (group 2) exhibit small BAL-velocity shifts along with steep-spectrum, weak radio emission at 3.0 and 1.4 GHz. All four quasars have spectral energy distributions (SEDs) with a peak at $λ_{\rm rest }\sim10~μ$m, suggesting a link between the BAL acceleration and hot dust emission. The group-2 quasars are redder than group-1 quasars and have a steeper rise at $1<λ_{\rm rest }<3~μ$m in their SEDs. All but J1238 exhibit a steep rise followed by a plateau-like time evolution in BAL-velocity shift. Our investigations, combined with previous studies of BAL acceleration, indicate that (1) the BAL-ISM coupling process is one of the major avenues for the origin of quasar reddening and patchy obscuration, (2) AAL outflows are ubiquitous and likely signify large-scale remnants of BAL winds coupled to interstellar medium (ISM), and (3) wind deceleration that is closely linked to the BAL-ISM coupling process may produce weak radio emission in otherwise radio-quiet quasars.
△ Less
Submitted 15 February, 2024;
originally announced February 2024.
-
Dust and Cold Gas Properties of Starburst HyLIRG-Quasars at $z \sim 2.5$
Authors:
Feng-Yuan Liu,
Y. Sophia Dai,
Alain Omont,
Daizhong Liu,
Pierre Cox,
Roberto Neri,
Melanie Krips,
Chentao Yang,
Xue-Bing Wu,
Jia-Sheng Huang
Abstract:
Some high-z active galactic nuclei (AGNs) are found to reside in extreme star-forming galaxies, such as hyper-luminous infrared galaxies (HyLIRGs), with AGN-removed $L_{\rm{IR}}$ of $>10^{13} L_{\rm{\odot}}$. In this paper, we report NOEMA observations of six apparent starburst HyLIRGs associated with optical quasars at $z\sim2-3$ in the Stripe 82 field, to study their dust and molecular CO proper…
▽ More
Some high-z active galactic nuclei (AGNs) are found to reside in extreme star-forming galaxies, such as hyper-luminous infrared galaxies (HyLIRGs), with AGN-removed $L_{\rm{IR}}$ of $>10^{13} L_{\rm{\odot}}$. In this paper, we report NOEMA observations of six apparent starburst HyLIRGs associated with optical quasars at $z\sim2-3$ in the Stripe 82 field, to study their dust and molecular CO properties. Five out of the six candidates are detected with CO(4-3) or CO(5-4) emission, and four in 2mm dust continuum. Based on the linewidth-$L'_{\rm{CO(1-0)}}$ diagnostics, we find that four galaxies are likely unlensed or weakly lensed sources. The molecular gas mass is in the range of $μM_{\rm{H_2}} \sim0.8-9.7\times10^{10} M_{\odot}$ (with $α= 0.8 M_{\odot} (\rm{K km s^{-1} pc^2})^{-1}$ and $μ$ is the unknown possible gravitational magnification factor). We fit their SEDs, after including the observed 2mm fluxes and upper limits, and estimate their apparent (uncorrected for possible lensing effect) star formation rates ($μ$SFRs) to be $\sim400-2500$ $M_{\rm{\odot}} \rm{yr^{-1}}$ with depletion time of $\sim20-110$ Myr. We notice interesting offsets, of $\sim10-40$ kpc spatially or $\sim1000-2000$ km s$^{-1}$ spectroscopically, between the optical quasar and the mm continuum or CO emissions. The observed velocity shift is likely related to the blueshifted broad-emission-line region of quasars, though mergers or recoiling black holes are also possible causes, which can explain the spatial offset and the high intrinsic SFRs in the HyLIRG-quasar systems.
△ Less
Submitted 6 February, 2024;
originally announced February 2024.
-
Nuclear mass table in deformed relativistic Hartree-Bogoliubov theory in continuum, II: Even-$Z$ nuclei
Authors:
DRHBc Mass Table Collaboration,
Peng Guo,
Xiaojie Cao,
Kangmin Chen,
Zhihui Chen,
Myung-Ki Cheoun,
Yong-Beom Choi,
Pak Chung Lam,
Wenmin Deng,
Jianmin Dong,
Pengxiang Du,
Xiaokai Du,
Kangda Duan,
Xiaohua Fan,
Wei Gao,
Lisheng Geng,
Eunja Ha,
Xiao-Tao He,
Jinniu Hu,
Jingke Huang,
Kun Huang,
Yanan Huang,
Zidan Huang,
Kim Da Hyung,
Hoi Yat Chan
, et al. (58 additional authors not shown)
Abstract:
The mass table in the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc) with the PC-PK1 density functional has been established for even-$Z$ nuclei with $8\le Z\le120$, extended from the previous work for even-even nuclei [Zhang $\it{et.~al.}$ (DRHBc Mass Table Collaboration), At. Data Nucl. Data Tables 144, 101488 (2022)]. The calculated binding energies, two-nucleon and one-ne…
▽ More
The mass table in the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc) with the PC-PK1 density functional has been established for even-$Z$ nuclei with $8\le Z\le120$, extended from the previous work for even-even nuclei [Zhang $\it{et.~al.}$ (DRHBc Mass Table Collaboration), At. Data Nucl. Data Tables 144, 101488 (2022)]. The calculated binding energies, two-nucleon and one-neutron separation energies, root-mean-square (rms) radii of neutron, proton, matter, and charge distributions, quadrupole deformations, and neutron and proton Fermi surfaces are tabulated and compared with available experimental data. A total of 4829 even-$Z$ nuclei are predicted to be bound, with an rms deviation of 1.477 MeV from the 1244 mass data. Good agreement with the available experimental odd-even mass differences, $α$ decay energies, and charge radii is also achieved. The description accuracy for nuclear masses and nucleon separation energies as well as the prediction for drip lines is compared with the results obtained from other relativistic and nonrelativistic density functional. The comparison shows that the DRHBc theory with PC-PK1 provides an excellent microscopic description for the masses of even-$Z$ nuclei. The systematics of the nucleon separation energies, odd-even mass differences, pairing energies, two-nucleon gaps, $α$ decay energies, rms radii, quadrupole deformations, potential energy curves, neutron density distributions, and neutron mean-field potentials are discussed.
△ Less
Submitted 10 June, 2024; v1 submitted 5 February, 2024;
originally announced February 2024.
-
On the Broadening of the Pulse Width of FRB 20121102A due to Propagation and Instrumental Effects
Authors:
Jia-Peng Wei,
Yong-Feng Huang,
Lang Cui,
Xiang Liu,
Jin-Jun Geng,
Xue-Feng Wu
Abstract:
The pulse widths of fast radio bursts are always broadened due to the scattering of the plasma medium through which the electromagnetic wave passes. The recorded pulse width will be further affected by the radio telescopes since the sampling time and the bandwidth cannot be infinitely small. In this study, we focus on the pulse widths of the 3287 bursts detected from FRB 20121102A as of October 20…
▽ More
The pulse widths of fast radio bursts are always broadened due to the scattering of the plasma medium through which the electromagnetic wave passes. The recorded pulse width will be further affected by the radio telescopes since the sampling time and the bandwidth cannot be infinitely small. In this study, we focus on the pulse widths of the 3287 bursts detected from FRB 20121102A as of October 2023. Various effects such as the scattering broadening, the redshift broadening and the instrumental broadening are examined. It is found that the instrumental broadening only contributes a fraction of $10^{-3}$--$10^{-1}$ to the observed pulse width. The scattering broadening is even smaller, which constitutes a tiny fraction of $10^{-5}$--$10^{-2}$ in the observed pulse width. After correcting for these broadenings, the intrinsic pulse width is derived for each burst. The maximum and minimum pulse widths at different frequencies are highlighted. Interestingly, both the mean value and the dispersion range of intrinsic pulse width are found to be inversely proportional to the square of the central frequency. The intrinsic widths of most bursts are in a narrow range of 1--10 ms, which leads to a quasi-linear correlation between the fluence and the peak flux.
△ Less
Submitted 4 February, 2024;
originally announced February 2024.
-
Detecting Extragalactic Axion-like Dark Matter with Polarization Measurements of Fast Radio Bursts
Authors:
Bao Wang,
Xuan Yang,
Jun-Jie Wei,
Song-Bo Zhang,
Xue-Feng Wu
Abstract:
Axions or axion-like particles (ALPs) are one of the most promising candidates of dark matter (DM). A prevalent method to detect axion-like DM is to seek the periodic oscillation feature in the polarization angles of linearly polarized light emitted from astrophysical sources. In this work, we use the time-resolved polarization measurements of the hyperactive repeating fast radio burst, FRB 202209…
▽ More
Axions or axion-like particles (ALPs) are one of the most promising candidates of dark matter (DM). A prevalent method to detect axion-like DM is to seek the periodic oscillation feature in the polarization angles of linearly polarized light emitted from astrophysical sources. In this work, we use the time-resolved polarization measurements of the hyperactive repeating fast radio burst, FRB 20220912A, detected by the Five-hundred-meter Aperture Spherical radio Telescope (FAST) to search for extragalactic axion-like DM for the first time. Given a DM density profile of FRB 20220912A's host, we obtain upper limits on the ALP-photon coupling constant of $g_{a γ}<(2.9 \times 10^{-11}-1.1\times 10^{-9})\,\mathrm{GeV}^{-1}$ for the ALP masses $m_a \sim (1.4\times10^{-21}-5.2\times10^{-20})$ eV. Persistent polarimetric observations with FAST would further improve the constraints. We prove that FRBs offer an alternative way to detect axion-like DM on extragalactic distance scales, complementary to other galactic DM probes.
△ Less
Submitted 1 February, 2024;
originally announced February 2024.
-
A spectral data release for 104 Type II Supernovae from the Tsinghua Supernova Group
Authors:
Han Lin,
Xiaofeng Wang,
Jujia Zhang,
Danfeng Xiang,
Tianmeng Zhang,
Xulin Zhao,
Xinghan Zhang,
Hanna Sai,
Liming Rui,
Jun Mo,
Gaobo Xi,
Fang Huang,
Xue Li,
Yongzhi Cai,
Weili Lin,
Jie Lin,
Chengyuan Wu,
Jicheng Zhang,
Zhihao Chen,
Zhitong Li,
Wenxiong Li,
Linyi Li,
Kaicheng Zhang,
Cheng Miao,
Juncheng Chen
, et al. (11 additional authors not shown)
Abstract:
We present 206 unpublished optical spectra of 104 type II supernovae obtained by the Xinglong 2.16m telescope and Lijiang 2.4m telescope during the period from 2011 to 2018, spanning the phases from about 1 to 200 days after the SN explosion. The spectral line identifications, evolution of line velocities and pseudo equivalent widths, as well as correlations between some important spectral paramet…
▽ More
We present 206 unpublished optical spectra of 104 type II supernovae obtained by the Xinglong 2.16m telescope and Lijiang 2.4m telescope during the period from 2011 to 2018, spanning the phases from about 1 to 200 days after the SN explosion. The spectral line identifications, evolution of line velocities and pseudo equivalent widths, as well as correlations between some important spectral parameters are presented. Our sample displays a large range in expansion velocities. For instance, the Fe~{\sc ii} $5169$ velocities measured from spectra at $t\sim 50$ days after the explosion vary from ${\rm 2000\ km\ s^{-1}}$ to ${\rm 5500\ km\ s^{-1}}$, with an average value of ${\rm 3872 \pm 949\ km\ s^{-1}}$. Power-law functions can be used to fit the velocity evolution, with the power-law exponent quantifying the velocity decline rate. We found an anticorrelation existing between H$β$ velocity at mid-plateau phase and its velocity decay exponent, SNe II with higher velocities tending to have smaller velocity decay rate. Moreover, we noticed that the velocity decay rate inferred from the Balmer lines (i.e., H$α$ and H$β$) have moderate correlations with the ratio of absorption to emission for H$α$ (a/e). In our sample, two objects show possibly flash-ionized features at early phases. Besides, we noticed that multiple high-velocity components may exist on the blue side of hydrogen lines of SN 2013ab, possibly suggesting that these features arise from complex line forming region. All our spectra can be found in WISeREP and Zenodo.
△ Less
Submitted 11 January, 2024;
originally announced January 2024.
-
Formation of collisional ring galaxies in Milgromian dynamics
Authors:
Li Ma,
Xufen Wu
Abstract:
Ring galaxies are rare in the Universe. A head-on or off-centre collision between an intruder galaxy and a disc galaxy leads to a collisional ring galaxy (CRG) when the intruder-to-target mass ratio (ITMR) is over 0.1 in Newtonian dynamics. Within the framework of Milgromian dynamics, the strong equivalence principle is violated due to the presence of an external field. When two galaxies collide t…
▽ More
Ring galaxies are rare in the Universe. A head-on or off-centre collision between an intruder galaxy and a disc galaxy leads to a collisional ring galaxy (CRG) when the intruder-to-target mass ratio (ITMR) is over 0.1 in Newtonian dynamics. Within the framework of Milgromian dynamics, the strong equivalence principle is violated due to the presence of an external field. When two galaxies collide towards each other, the dynamical mass of the phantom dark halo in a galaxy is suppressed by the external field induced by the other galaxy. As a consequence of such suppression, the gravitational perturbation for the target galaxy introduced by a collision is weakened. In this case, a minor collision may not be capable of generating a CRG. In this work, we address this issue and perform a series of numerical simulations of collisions by tuning the values of ITMR. We find that the critical ITMR is 0.5 in MOND, which is much larger than that in Newtonian dynamics. The observed massive ring galaxies, such as Arp 147, can be effectively interpreted by CRGs in MOND. This interpretation does not necessitate the presence of dark matter halos for either the target or intruder galaxies. Moreover, for a larger inclination angle or a larger impact parameter, the off-centred ring structure is fainter. The larger critical ITMR indicates that it is harder to form a CRG in Milgrom's Modified Newtonian Dynamics (MOND). To account for the observed ring structures of the NGC 922-like galaxies in MOND, it is necessary to invoke other scenarios than a dry minor collision.
△ Less
Submitted 3 January, 2024;
originally announced January 2024.
-
The Propagation of Fast Radio Bursts in the Magnetosphere Shapes Their Waiting-time and Flux Distributions
Authors:
Di Xiao,
Zi-Gao Dai,
Xue-Feng Wu
Abstract:
The field of fast radio bursts (FRBs) has entered the age of fine characterization as observational results from different radio telescopes become more and more abundant. The large FRB sample is suitable for a statistical study. There is an interesting finding that the waiting-time distributions of very active repeating FRBs show a universal double-peaked feature, with left peaks lower than right…
▽ More
The field of fast radio bursts (FRBs) has entered the age of fine characterization as observational results from different radio telescopes become more and more abundant. The large FRB sample is suitable for a statistical study. There is an interesting finding that the waiting-time distributions of very active repeating FRBs show a universal double-peaked feature, with left peaks lower than right ones. Assuming these two peaks are independent and initially comparable, we show that the observed asymmetric shape can be ascribed to the propagational effect in the magnetosphere. An FRB passing through the magnetized plasma will induce the circular motion of charged particles to form a current loop. This further leads to an induced magnetic field with opposite direction respect to the background field. As the effective field strength changes, the scattering absorption probability of the following FRB will be influenced. The absorption can be important under certain physical conditions and bursts with smaller time-lags are easier to be absorbed. Also there will be an imprint on the flux distribution as the scattering optical depth depends on burst luminosity as well.
△ Less
Submitted 21 December, 2023;
originally announced December 2023.
-
Non-perturbative study of quantum many-body correlation effects in neutron stars: Equation of state
Authors:
Hao-Fu Zhu,
Xufen Wu,
Guo-Zhu Liu
Abstract:
Although neutron stars have been studied for decades, their internal structure remains enigmatic, mainly due to large uncertainties in the equation of state. In neutron stars, nucleons are strongly interacting by exchanging mesons, which leads to significant quantum many-body correlation effects. Previous mean-field calculations failed to capture these effects. Here, we develop a non-perturbative…
▽ More
Although neutron stars have been studied for decades, their internal structure remains enigmatic, mainly due to large uncertainties in the equation of state. In neutron stars, nucleons are strongly interacting by exchanging mesons, which leads to significant quantum many-body correlation effects. Previous mean-field calculations failed to capture these effects. Here, we develop a non-perturbative quantum field-theoretic approach to handle strongly correlated nuclear matters. Based on a linear $σ$-$ω$-$ρ$ model with a density-dependent isovector coupling $Γ_ρ$, we reproduce all the six frequently-used nuclear matter observables by tuning six parameters. After including quantum many-body effects into the equation of state of realistic neutron-star matters, we obtain a mass-radius relation that is in good agreement with recent astrophysical observations. Our results show that the maximal neutron-star mass supported by neutron-dominated matters is $\sim 2.67M_{\odot}$, where $M_{\odot}$ is solar mass.
△ Less
Submitted 15 April, 2024; v1 submitted 19 December, 2023;
originally announced December 2023.
-
Probing Thermal Electrons in GRB Afterglows
Authors:
Hao-Xuan Gao,
Jin-Jun Geng,
Tian-Rui Sun,
Liang Li,
Yong-Feng Huang,
Xue-Feng Wu
Abstract:
Particle-in-cell simulations have unveiled that shock-accelerated electrons do not follow a pure power-law distribution, but have an additional low-energy "thermal" part, which owns a considerable portion of the total energy of electrons. Investigating the effects of these thermal electrons on gamma-ray burst (GRB) afterglows may provide valuable insights into the particle acceleration mechanisms.…
▽ More
Particle-in-cell simulations have unveiled that shock-accelerated electrons do not follow a pure power-law distribution, but have an additional low-energy "thermal" part, which owns a considerable portion of the total energy of electrons. Investigating the effects of these thermal electrons on gamma-ray burst (GRB) afterglows may provide valuable insights into the particle acceleration mechanisms. We solve the continuity equation of electrons in the energy space, from which multi-wavelength afterglows are derived by incorporating processes including synchrotron radiation, synchrotron self-absorption, synchrotron self-Compton scattering, and gamma-gamma annihilation. First, there is an underlying positive correlation between temporal and spectral indices due to the cooling of electrons. Moreover, thermal electrons would result in the simultaneous non-monotonic variation in both spectral and temporal indices at multi-wavelength, which could be individually recorded by the 2.5-meter Wide Field Survey Telescope and Vera Rubin Observatory Legacy Survey of Space and Time (LSST). The thermal electrons could also be diagnosed from afterglow spectra by synergy observation in the optical (with LSST) and X-ray bands (with the Microchannel X-ray Telescope on board the Space Variable Objects Monitor). Finally, we use Monte Carlo simulations to obtain the distribution of peak flux ratio ($R_{\rm X}$) between soft and hard X-rays, and of the time delay ($Δt$) between peak times of soft X-ray and optical light curves. The thermal electrons significantly raise the upper limits of both $R_{\rm X}$ and $Δt$. Thus the distribution of GRB afterglows with thermal electrons is more dispersive in the $R_{\rm X} - Δt$ plane.
△ Less
Submitted 10 June, 2024; v1 submitted 10 November, 2023;
originally announced November 2023.
-
Does or did the supernova remnant Cassiopeia A operate as a PeVatron?
Authors:
Zhen Cao,
F. Aharonian,
Q. An,
Axikegu,
Y. X. Bai,
Y. W. Bao,
D. Bastieri,
X. J. Bi,
Y. J. Bi,
J. T. Cai,
Q. Cao,
W. Y. Cao,
Zhe Cao,
J. Chang,
J. F. Chang,
A. M. Chen,
E. S. Chen,
Liang Chen,
Lin Chen,
Long Chen,
M. J. Chen,
M. L. Chen,
Q. H. Chen,
S. H. Chen,
S. Z. Chen
, et al. (255 additional authors not shown)
Abstract:
For decades, supernova remnants (SNRs) have been considered the prime sources of Galactic Cosmic rays (CRs). But whether SNRs can accelerate CR protons to PeV energies and thus dominate CR flux up to the knee is currently under intensive theoretical and phenomenological debate. The direct test of the ability of SNRs to operate as CR PeVatrons can be provided by ultrahigh-energy (UHE;…
▽ More
For decades, supernova remnants (SNRs) have been considered the prime sources of Galactic Cosmic rays (CRs). But whether SNRs can accelerate CR protons to PeV energies and thus dominate CR flux up to the knee is currently under intensive theoretical and phenomenological debate. The direct test of the ability of SNRs to operate as CR PeVatrons can be provided by ultrahigh-energy (UHE; $E_γ\geq 100$~TeV) $γ$-rays. In this context, the historical SNR Cassiopeia A (Cas A) is considered one of the most promising target for UHE observations. This paper presents the observation of Cas A and its vicinity by the LHAASO KM2A detector. The exceptional sensitivity of LHAASO KM2A in the UHE band, combined with the young age of Cas A, enabled us to derive stringent model-independent limits on the energy budget of UHE protons and nuclei accelerated by Cas A at any epoch after the explosion. The results challenge the prevailing paradigm that Cas A-type SNRs are major suppliers of PeV CRs in the Milky Way.
△ Less
Submitted 25 October, 2023;
originally announced October 2023.
-
CatNorth: An Improved Gaia DR3 Quasar Candidate Catalog with Pan-STARRS1 and CatWISE
Authors:
Yuming Fu,
Xue-Bing Wu,
Yifan Li,
Yuxuan Pang,
Ravi Joshi,
Shuo Zhang,
Qiyue Wang,
Jing Yang,
FanLam Ng,
Xingjian Liu,
Yu Qiu,
Rui Zhu,
Huimei Wang,
Christian Wolf,
Yanxia Zhang,
Zhi-Ying Huo,
Y. L. Ai,
Qinchun Ma,
Xiaotong Feng,
R. J. Bouwens
Abstract:
A complete and pure sample of quasars with accurate redshifts is crucial for quasar studies and cosmology. In this paper, we present CatNorth, an improved Gaia DR3 quasar candidate catalog with more than 1.5 million sources in the 3$π$ sky built with data from Gaia, Pan-STARRS1, and CatWISE2020. The XGBoost algorithm is used to reclassify the original Gaia DR3 quasar candidates as stars, galaxies,…
▽ More
A complete and pure sample of quasars with accurate redshifts is crucial for quasar studies and cosmology. In this paper, we present CatNorth, an improved Gaia DR3 quasar candidate catalog with more than 1.5 million sources in the 3$π$ sky built with data from Gaia, Pan-STARRS1, and CatWISE2020. The XGBoost algorithm is used to reclassify the original Gaia DR3 quasar candidates as stars, galaxies, and quasars. To construct training/validation datasets for the classification, we carefully built two different master stellar samples in addition to the spectroscopic galaxy and quasar samples. An ensemble classification model is obtained by averaging two XGBoost classifiers trained with different master stellar samples. Using a probability threshold of $p_{\mathrm{QSO\_mean}}>0.95$ in our ensemble classification model and an additional cut on the logarithmic probability density of zero proper motion, we retrieved 1,545,514 reliable quasar candidates from the parent Gaia DR3 quasar candidate catalog. We provide photometric redshifts for all candidates with an ensemble regression model. For a subset of 89,100 candidates, accurate spectroscopic redshifts are estimated with the Convolutional Neural Network from the Gaia BP/RP spectra. The CatNorth catalog has a high purity of ~ 90% while maintaining high completeness, which is an ideal sample to understand the quasar population and its statistical properties. The CatNorth catalog is used as the main source of input catalog for the LAMOST phase III quasar survey, which is expected to build a highly complete sample of bright quasars with $i < 19.5$.
△ Less
Submitted 13 February, 2024; v1 submitted 19 October, 2023;
originally announced October 2023.
-
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.
-
A bright burst from FRB 20200120E in a globular cluster of the nearby galaxy M81
Authors:
S. B. Zhang,
J. S. Wang,
X. Yang,
Y. Li,
J. J. Geng,
Z. F. Tang,
C. M. Chang,
J. T. Luo,
X. C. Wang,
X. F. Wu,
Z. G. Dai,
B. Zhang
Abstract:
Fast radio bursts (FRBs) are immensely energetic millisecond-duration radio pulses. Observations indicate that nearby FRBs can be produced by old stellar populations, as suggested by the localization of the repeating source FRB 20200120E in a globular cluster of M81. Nevertheless, the burst energies of FRB 20200120E are significantly smaller than those of other cosmological FRBs, even falling belo…
▽ More
Fast radio bursts (FRBs) are immensely energetic millisecond-duration radio pulses. Observations indicate that nearby FRBs can be produced by old stellar populations, as suggested by the localization of the repeating source FRB 20200120E in a globular cluster of M81. Nevertheless, the burst energies of FRB 20200120E are significantly smaller than those of other cosmological FRBs, even falling below the energy of the Galactic event FRB 20200428. Here, we report the detection of a bright burst from FRB 20200120E in 1.1 -- 1.7 GHz, with a fluence of about 30 Jy ms, which is more than 42 times larger than the previously detected bursts near 1.4 GHz frequency. It reaches one-third of the energy of the weakest burst from FRB 20121102A and is detectable at a distance exceeding 200 Mpc. Our finding bridges the gap between nearby and cosmological FRBs and indicates that FRBs hosted in globular clusters can be bright enough to be observable at cosmological distances.
△ Less
Submitted 31 July, 2024; v1 submitted 2 October, 2023;
originally announced October 2023.
-
Periodic Variable Star Classification with Deep Learning: Handling Data Imbalance in an Ensemble Augmentation Way
Authors:
Zihan Kang,
Yanxia Zhang,
Jingyi Zhang,
Changhua Li,
Minzhi Kong,
Yongheng Zhao,
Xue-Bing Wu
Abstract:
Time-domain astronomy is progressing rapidly with the ongoing and upcoming large-scale photometric sky surveys led by the Vera C. Rubin Observatory project (LSST). Billions of variable sources call for better automatic classification algorithms for light curves. Among them, periodic variable stars are frequently studied. Different categories of periodic variable stars have a high degree of class i…
▽ More
Time-domain astronomy is progressing rapidly with the ongoing and upcoming large-scale photometric sky surveys led by the Vera C. Rubin Observatory project (LSST). Billions of variable sources call for better automatic classification algorithms for light curves. Among them, periodic variable stars are frequently studied. Different categories of periodic variable stars have a high degree of class imbalance and pose a challenge to algorithms including deep learning methods. We design two kinds of architectures of neural networks for the classification of periodic variable stars in the Catalina Survey's Data Release 2: a multi-input recurrent neural network (RNN) and a compound network combing the RNN and the convolutional neural network (CNN). To deal with class imbalance, we apply Gaussian Process to generate synthetic light curves with artificial uncertainties for data augmentation. For better performance, we organize the augmentation and training process in a "bagging-like" ensemble learning scheme. The experimental results show that the better approach is the compound network combing RNN and CNN, which reaches the best result of 86.2% on the overall balanced accuracy and 0.75 on the macro F1 score. We develop the ensemble augmentation method to solve the data imbalance when classifying variable stars and prove the effectiveness of combining different representations of light curves in a single model. The proposed methods would help build better classification algorithms of periodic time series data for future sky surveys (e.g., LSST).
△ Less
Submitted 24 September, 2023;
originally announced September 2023.
-
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.
-
Refractive lensing of scintillating FRBs by sub-parsec cloudlets in the multi-phase CGM
Authors:
Dylan L. Jow,
Xiaohan Wu,
Ue-Li Pen
Abstract:
We consider the refractive lensing effects of ionized cool ($T \sim 10^4\,{\rm K}$) gas cloudlets in the circumgalactic medium (CGM) of galaxies. In particular, we discuss the combined effects of lensing from these cloudlets and scintillation from plasma screens in the Milky Way interstellar medium (ISM). We show that, if the CGM comprises a mist of sub-parsec cloudlets with column densities of or…
▽ More
We consider the refractive lensing effects of ionized cool ($T \sim 10^4\,{\rm K}$) gas cloudlets in the circumgalactic medium (CGM) of galaxies. In particular, we discuss the combined effects of lensing from these cloudlets and scintillation from plasma screens in the Milky Way interstellar medium (ISM). We show that, if the CGM comprises a mist of sub-parsec cloudlets with column densities of order $10^{17}\,{\rm cm}^{-2}$ (as predicted by McCourt et al. (2018)), then FRBs whose sightlines pass within a virial radius of a CGM halo will generically be lensed into tens of refractive images with a $\sim 10\,{\rm ms}$ scattering timescale. These images will be resolved by scintillating screens in the Milky Way ISM, and therefore are expected to suppress scintillation. From this, we argue that positive detections of FRB scintillation can constrain the properties of these cool-gas cloudlets, with current scintillation observation disfavouring the cloudlet model. We propose that sheet-like geometries for the cool gas in the CGM can reconcile quasar absorption measurements (from which we infer the presence of the cool gas with structure on sub-parsec scales) and the unexpected lack of lensing signals from this gas thus far observed.
△ Less
Submitted 13 September, 2023;
originally announced September 2023.
-
POLAR-2, the next generation of GRB polarization detector
Authors:
Nicolas Produit,
Merlin Kole,
Xin Wu,
Nicolas De Angelis,
Hancheng Li,
Dominik Rybka,
Agnieszka Pollo,
Slawomir Mianowski,
Jochen Greiner,
J. Michael Burgess,
Jianchao Sun,
Shuang-Nan Zhang
Abstract:
The POLAR-2 Gamma-Ray Burst (GRB) Polarimetry mission is a follow-up to the successful POLAR mission. POLAR collected six months of data in 2016-2017 on board the Tiangong-2 Chinese Space laboratory. From a polarization study on 14 GRBs, POLAR measured an overall low polarization and a hint for an unexpected complexity in the time evolution of polarization during GRBs. Energy-dependent measurement…
▽ More
The POLAR-2 Gamma-Ray Burst (GRB) Polarimetry mission is a follow-up to the successful POLAR mission. POLAR collected six months of data in 2016-2017 on board the Tiangong-2 Chinese Space laboratory. From a polarization study on 14 GRBs, POLAR measured an overall low polarization and a hint for an unexpected complexity in the time evolution of polarization during GRBs. Energy-dependent measurements of the GRB polarization will be presented by N. de Angelis in GA21-09 (August 2nd). These results demonstrate the need for measurements with significantly improved accuracy. Moreover, the recent discovery of gravitational waves and their connection to GRBs justifies a high-precision GRB polarimeter that can provide both high-precision polarimetry and detection of very faint GRBs. The POLAR-2 polarimeter is based on the same Compton scattering measurement principle as POLAR, but with an extended energy range and an order of magnitude increase in total effective area for polarized events. Proposed and developed by a joint effort of Switzerland, China, Poland and Germany, the device was selected for installation on the China Space Station and is scheduled to start operation for at least 2 years in 2025.
△ Less
Submitted 1 September, 2023;
originally announced September 2023.
-
Energy-dependent polarization of Gamma-Ray Bursts' prompt emission with the POLAR and POLAR-2 instruments
Authors:
Nicolas De Angelis,
J. Michael Burgess,
Franck Cadoux,
Jochen Greiner,
Merlin Kole,
Hancheng Li,
Slawomir Mianowski,
Agnieszka Pollo,
Nicolas Produit,
Dominik Rybka,
Jianchao Sun,
Xin Wu,
Shuang-Nan Zhang
Abstract:
Gamma-Ray Bursts are among the most powerful events in the Universe. Despite half a century of observations of these transient sources, many open questions remain about their nature. Polarization measurements of the GRB prompt emission have long been theorized to be able to answer most of these questions. With the aim of characterizing the polarization of these prompt emissions, a compact Compton…
▽ More
Gamma-Ray Bursts are among the most powerful events in the Universe. Despite half a century of observations of these transient sources, many open questions remain about their nature. Polarization measurements of the GRB prompt emission have long been theorized to be able to answer most of these questions. With the aim of characterizing the polarization of these prompt emissions, a compact Compton polarimeter, called POLAR, has been launched to space in September 2016. Time integrated polarization analysis of the POLAR GRB catalog have shown that the prompt emission is lowly polarized or fully unpolarized. However, time resolved analysis depicted strong hints of an evolving polarization angle within single pulses, washing out the polarization degree in time integrated analyses. Here we will for the first time present energy resolved polarization measurements with the POLAR data. The novel analysis, performed on several GRBs, will provide new insights and alter our understanding of GRB polarization. The analysis was performed using the 3ML framework to fit polarization parameters versus energy in parallel to the spectral parameters. Although limited by statistics, the results could provide a very relevant input to disentangle between existing theoretical models. In order to gather more statistics per GRB and perform joint time and energy resolved analysis, a successor instrument, called POLAR-2, is under development with a launch window early 2025 to the CSS. After presenting the first energy resolved polarization results of the POLAR mission, we will present the prospects for such measurements with the upcoming POLAR-2 mission.
△ Less
Submitted 1 September, 2023;
originally announced September 2023.
-
Electromagnetic signatures of white dwarf collisions in AGN discs
Authors:
Shu-Rui Zhang,
Yan Luo,
Xiao-Jun Wu,
Jian-Min Wang,
Luis C. Ho,
Ye-Fei Yuan
Abstract:
In the inner region of the disc of an active galactic nucleus (AGN), the collision of two white dwarfs (WDs) through Jacobi capture might be inevitable, leading to a Type Ia supernova (SN Ia) explosion. This transient event, influenced by the disc gas and the gravity of the supermassive black hole (SMBH), exhibits distinct characteristics compared to normal SNe Ia. The energy of the explosion is m…
▽ More
In the inner region of the disc of an active galactic nucleus (AGN), the collision of two white dwarfs (WDs) through Jacobi capture might be inevitable, leading to a Type Ia supernova (SN Ia) explosion. This transient event, influenced by the disc gas and the gravity of the supermassive black hole (SMBH), exhibits distinct characteristics compared to normal SNe Ia. The energy of the explosion is mainly stored in the ejecta in the form of kinetic energy. Typically, the ejecta is not effectively decelerated by the AGN disc and rushes rapidly out of the AGN disc. However, under the influence of the SMBH, most of the ejecta falls back toward the AGN disc. As the fallback ejecta becomes more dispersed, it interacts with the disc gas, converting its kinetic energy into thermal energy. This results in a high-energy transient characterized by a rapid initial rise followed by a decay with $L\propto t^{-2.8}$. The time-scale of the transient ranges from hours to weeks, depending on the mass of the SMBH. This process generates high-energy radiation spanning from hard X-rays to the soft $γ$ range. Additionally, the subsequent damage to the disc may result in changing-look AGNs. Moreover, the falling back of SNe Ia ejecta onto the AGN disc significantly increases the metallicity of the AGN and can even generate heavy elements within the AGN discs.
△ Less
Submitted 28 August, 2023;
originally announced August 2023.
-
White Dwarf--White Dwarf collisions in AGN discs via close encounters
Authors:
Yan Luo,
Xiao-Jun Wu,
Shu-Rui Zhang,
Jian-Min Wang,
Luis C. Ho,
Ye-Fei Yuan
Abstract:
White dwarfs (WDs) in active galactic nucleus (AGNs) discs might migrate to the inner radii of the discs and form restricted three-body systems with two WDs moving around the central supermassive black hole (SMBH) in close orbits. These systems could be dynamical unstable, which can lead to very close encounters or direct collisions. In this work, we use N-body simulations to study the evolution o…
▽ More
White dwarfs (WDs) in active galactic nucleus (AGNs) discs might migrate to the inner radii of the discs and form restricted three-body systems with two WDs moving around the central supermassive black hole (SMBH) in close orbits. These systems could be dynamical unstable, which can lead to very close encounters or direct collisions. In this work, we use N-body simulations to study the evolution of such systems with the different initial orbital separation $p$, relative orbital inclination $Δ{i}$ and SMBH mass $M$. It is found that the close encounters of WDs mainly occur at $1.1R_{\rm H} \lesssim p \lesssim 2\sqrt{3}R_{\rm H}$, where $R_{\rm H}$ is the mutual Hill radius. For $p<1.1R_{\rm H}$, the majority of WDs move in horseshoe or tadpole orbits, and only few of them with small initial orbital phase difference undergo close encounters. For $p=3.0R_{\rm H}$, WD-WD collisions occur in most of the samples within a time of $10^5P_1$, and considerable collisions occur within a time of $t<62P_1$ for small orbital radii, where $P_1$ is the orbital period. The peak of the closest separation distribution increase and the WD-WD collision fraction decreases with an increase of the relative inclination. The closest separation distribution is similar in cases with the different SMBH mass, but the WD-WD collision fraction decreases as the mass of SMBHs increases. According to our estimation, the event rate of the cosmic WD-WD collision in AGN discs is about $300{\rm Gpc^{-3}yr^{-1}}$, roughly $1\%$ of the one of the observed type Ia supernova. The corresponding electromagnetic emission signals can be observed by large surveys of AGNs.
△ Less
Submitted 28 August, 2023;
originally announced August 2023.
-
Effects of coupling constants on chaos of charged particles in the $Einstein-Æ$ ther theory
Authors:
Caiyu Liu,
Xin Wu
Abstract:
There are two free coupling parameters $c_{13}$ and $c_{14}$ in the Einstein-Æther metric describing a non-rotating black hole. This metric is the Reissner-Nordström black hole solution when $0\leq 2c_{13}<c_{14}<2$, but it is not for $0\leq c_{14}<2c_{13}<2$. When the black hole is immersed in an external asymptotically uniform magnetic field, the Hamiltonian system describing the motion of charg…
▽ More
There are two free coupling parameters $c_{13}$ and $c_{14}$ in the Einstein-Æther metric describing a non-rotating black hole. This metric is the Reissner-Nordström black hole solution when $0\leq 2c_{13}<c_{14}<2$, but it is not for $0\leq c_{14}<2c_{13}<2$. When the black hole is immersed in an external asymptotically uniform magnetic field, the Hamiltonian system describing the motion of charged particles around the black hole is not integrable. However, the Hamiltonian allows for the construction of explicit symplectic integrators. The proposed fourth-order explicit symplectic scheme is used to investigate the dynamics of charged particles because it exhibits excellent long-term performance in conserving the Hamiltonian. No universal rule can be given to the dependence of regular and chaotic dynamics on varying one or two parameters $c_{13}$ and $c_{14}$ in the two cases of $0\leq 2c_{13}<c_{14}<2$ and $0\leq c_{14}<2c_{13}<2$. The distributions of order and chaos in the binary parameter space $(c_{13},c_{14})$ rely on different combinations of the other parameters and the initial conditions.
△ Less
Submitted 9 August, 2023;
originally announced August 2023.