Showing 1–50 of 111 results for author: Xia, C

A timing view of the additional high-energy spectral component discovered in the black hole candidate Swift J1727.8-1613

Authors: Zi-Xu Yang, Liang Zhang, Shuang-Nan Zhang, L. Tao, Shu Zhang, Ruican Ma, Qingcui Bu, Yue Huang, He-Xin Liu, Wei Yu, Guang C. Xiao, Peng-Ju Wang, Hua Feng, Li-Ming Song, Xiang Ma, Mingyu Ge, QingChang Zhao, J. L. Qu

Abstract: We present an energy-dependent analysis for the type-C quasi-periodic oscillations (QPOs) observed in the black hole X-ray binary Swift J1727.8-1613 using Insight-HXMT observations. We find that the QPO fractional rms at energies above 40 keV is significantly higher than that below 20 keV. This is the first report of a high energy (HE)-rms excess in the rms spectrum of a black hole X-ray binary. I… ▽ More We present an energy-dependent analysis for the type-C quasi-periodic oscillations (QPOs) observed in the black hole X-ray binary Swift J1727.8-1613 using Insight-HXMT observations. We find that the QPO fractional rms at energies above 40 keV is significantly higher than that below 20 keV. This is the first report of a high energy (HE)-rms excess in the rms spectrum of a black hole X-ray binary. In the high energy band, an extra hard component is observed in additional to the standard thermal Comptonization component at similar energy band. The value of the QPO HE-rms excess is not only correlated with the disk parameters and the photon index of the standard Comptonization component, but also exhibits a moderate positive correlation with the flux of the additional hard spectral component. No features in the QPO phase-lag spectra are seen corresponding to the additional hard component. We propose that the additional hard component in the spectrum may originate from jet emission and the associated QPO HE-rms excess can be explained by the precession of the jet base. △ Less

Submitted 6 July, 2024; originally announced July 2024.

Estimation of radial velocities of BHB stars

Authors: Tahereh Ramezani, Ernst Paunzen, Caiyun Xia, Katerina Pivonkova, Prapti Mondal

Abstract: We studied blue horizontal branch stars (BHBs), and calculated their radial velocities. Spectra of these stars have been obtained with moderate signal-to-noise ratio for five blue horizontal-branch stars using the 2 meter telescope and Echelle Spectrograph in Ondrejov observatory, Czech republic. We studied blue horizontal branch stars (BHBs), and calculated their radial velocities. Spectra of these stars have been obtained with moderate signal-to-noise ratio for five blue horizontal-branch stars using the 2 meter telescope and Echelle Spectrograph in Ondrejov observatory, Czech republic. △ Less

Submitted 11 June, 2024; v1 submitted 22 May, 2024; originally announced May 2024.

Comments: 8 pages, 7 figures, 3 tables

Chemically peculiar stars on the pre-main sequence

Authors: L. Kueß, E. Paunzen, N. Faltová, D. Jadlovský, M. Labaj, M. Mesarč, P. Mondal, M. Prišegen, T. Ramezani, J. Supíková, K. Svačinková, M. Vítková, C. Xia, K. Bernhard, S. Hümmerich

Abstract: Context. The chemically peculiar (CP) stars of the upper main sequence are defined by spectral peculiarities that indicate unusual elemental abundance patterns in the presence of diffusion in the calm, stellar atmospheres. Some of them have a stable local magnetic field of up to several kiloGauss. The pre-main-sequence evolution of these objects is still a mystery and contains many open questions.… ▽ More Context. The chemically peculiar (CP) stars of the upper main sequence are defined by spectral peculiarities that indicate unusual elemental abundance patterns in the presence of diffusion in the calm, stellar atmospheres. Some of them have a stable local magnetic field of up to several kiloGauss. The pre-main-sequence evolution of these objects is still a mystery and contains many open questions. Aims. We identify CP stars on the pre-main sequence to determine possible mechanisms that lead to the occurrence of chemical peculiarities in the (very) early stages of stellar evolution. Methods. We identified likely pre-main-sequence stars by fitting the spectral energy distributions. The subsequent analysis using stellar spectra and photometric time series helped us to distinguish between CP and non-CP stars. Additionally, we compared our results to the literature to provide the best possible quality assessment. Results. Out of 45 candidates, about 70 % seem to be true CP stars or CP candidates. Furthermore, 9 sources appear to be CP stars on the pre-main sequence, and all are magnetic. We finally report a possible CP2 star that is also a pre-main-sequence star and was not previously in the literature. Conclusions. The evolution of the peculiarities seems to be related to the (strong) magnetic fields in these CP2 stars. △ Less

Submitted 14 May, 2024; originally announced May 2024.

Formation of Polar Crown Filaments Magnetic Fields by Supergranular Helicity Injection

Authors: Huanxin Chen, Chun Xia, Hechao Chen

Abstract: To understand the magnetic fields of the polar crown filaments (PCFs) at high latitudes near polar regions of the Sun, we perform magnetofrictional numerical simulations on the long-term magnetic evolution of bipolar fields with roughly east-west polarity inversion lines (PILs) in a three-dimensional (3D) spherical wedge domain near polar regions. The Coriolis effect induced vortical motions at th… ▽ More To understand the magnetic fields of the polar crown filaments (PCFs) at high latitudes near polar regions of the Sun, we perform magnetofrictional numerical simulations on the long-term magnetic evolution of bipolar fields with roughly east-west polarity inversion lines (PILs) in a three-dimensional (3D) spherical wedge domain near polar regions. The Coriolis effect induced vortical motions at the boundaries of several supergranular cells inject magnetic helicity from the photospheric boundary into the solar atmosphere. Supergranular-scale helicity injection, transfer, and condensation produce strongly sheared magnetic fields. Magnetic reconnections at footpoints of the sheared fields produce magnetic flux ropes (MFRs) with helicity signs consistent with the observed Hemispheric Helicity Rule (HHR). The cross-sectional area of MFRs exhibits an uneven distribution, resembling a "foot-node-foot" periodic configuration. Experiments with different tilt directions of PILs indicate that the PCFs preferably form along PILs with the western end close to the polar region. The bending of PILs caused by supergranular flows, forming S-shape (Z-shape) PIL segments, promotes the formation of dextral (sinistral) MFRs. The realistic magnetic models we got can serve as starting points for the study of the plasma formation and eruption of PCFs. △ Less

Submitted 19 March, 2024; originally announced March 2024.

Comments: 16 pages, 12 figures

Unveiling the Initiation Route of Coronal Mass Ejections through their Slow Rise Phase

Authors: Chen Xing, Guillaume Aulanier, Xin Cheng, Chun Xia, Mingde Ding

Abstract: Understanding the early evolution of coronal mass ejections (CMEs), in particular their initiation, is the key to forecasting solar eruptions and induced disastrous space weather. Although many initiation mechanisms have been proposed, a full understanding of CME initiation, which is identified as a slow rise of CME progenitors in kinematics before the impulsive acceleration, remains elusive. Here… ▽ More Understanding the early evolution of coronal mass ejections (CMEs), in particular their initiation, is the key to forecasting solar eruptions and induced disastrous space weather. Although many initiation mechanisms have been proposed, a full understanding of CME initiation, which is identified as a slow rise of CME progenitors in kinematics before the impulsive acceleration, remains elusive. Here, with a state-of-the-art thermal-magnetohydrodynamics simulation, we determine a complete CME initiation route in which multiple mainstream mechanisms occur in sequence yet are tightly coupled. The slow rise is first triggered and driven by the developing hyperbolic flux tube (HFT) reconnection. Subsequently, the slow rise continues as driven by the coupling of the HFT reconnection and the early development of torus instability. The end of the slow rise, i.e., the onset of the impulsive acceleration, is induced by the start of the fast magnetic reconnection coupled with the torus instability. These results unveil that the CME initiation is a complicated process involving multiple physical mechanisms, thus being hardly resolved by a single initiation mechanism. △ Less

Submitted 26 February, 2024; originally announced February 2024.

Comments: 35 pages, 15 figures, accepted for publication in ApJ

Boosted Dark Matter From Centaurus A and Its Detection

Authors: Chen Xia, Chuan-Yang Xing, Yan-Hao Xu

Abstract: Dark matter can be boosted by high energy particles in astrophysical environments through elastic scattering. We study the production of boosted dark matter via scattering with electrons in the relativistic jet of the closest active galactic nucleus, Centaurus A, and its detection in the Super-Kamiokande experiment. Since there are a huge number of electrons in the jet and dark matter is extremely… ▽ More Dark matter can be boosted by high energy particles in astrophysical environments through elastic scattering. We study the production of boosted dark matter via scattering with electrons in the relativistic jet of the closest active galactic nucleus, Centaurus A, and its detection in the Super-Kamiokande experiment. Since there are a huge number of electrons in the jet and dark matter is extremely dense around the supermassive black hole that powers the jet, the number of boosted dark matter is tremendously large. Compared to boosted dark matter from blazars, the dark matter flux from Centaurus A is enhanced due to the proximity of Centaurus A. The constraint on dark matter-electron scattering cross section set by Super-Kamiokande is more stringent, down to $\sim 10^{-36} \, \mathrm{cm}^2$ for $\mathrm{MeV}$ dark matter. △ Less

Submitted 14 March, 2024; v1 submitted 8 January, 2024; originally announced January 2024.

Comments: 18 pages, 5 figures, version appeared in JHEP

Journal ref: JHEP03(2024)076

Data-driven Modeling of a Coronal Magnetic Flux Rope: from Birth to Death

Authors: J. H. Guo, Y. W. Ni, Y. Guo, C. Xia, B. Schmieder, S. Poedts, Z. Zhong, Y. H. Zhou, F. Yu, P. F. Chen

Abstract: Magnetic flux ropes are a bundle of twisted magnetic field lines produced by internal electric currents, which are responsible for solar eruptions and are the major drivers of geomagnetic storms. As such, it is crucial to develop a numerical model that can capture the entire evolution of a flux rope, from its birth to death, in order to predict whether adverse space weather events might occur or n… ▽ More Magnetic flux ropes are a bundle of twisted magnetic field lines produced by internal electric currents, which are responsible for solar eruptions and are the major drivers of geomagnetic storms. As such, it is crucial to develop a numerical model that can capture the entire evolution of a flux rope, from its birth to death, in order to predict whether adverse space weather events might occur or not. In this paper, we develop a data-driven modeling that combines a time-dependent magneto-frictional approach with a thermodynamic magnetohydrodynamic model. Our numerical modeling successfully reproduces the formation and confined eruption of an observed flux rope, and unveils the physical details behind the observations. Regarding the long-term evolution of the active region, our simulation results indicate that the flux cancellation due to collisional shearing plays a critical role in the formation of the flux rope, corresponding to a substantial increase in magnetic free energy and helicity. Regarding the eruption stage, the deformation of the flux rope during its eruption can cause an increase in the downward tension force, which suppresses it from further rising. This finding may shed light on why some torus-unstable flux ropes lead to failed eruptions after large-angle rotations. Moreover, we find that twisted fluxes can accumulate during the confined eruptions, which would breed the subsequent eruptive flares. △ Less

Submitted 30 October, 2023; originally announced October 2023.

Comments: 30 pages, 10 figures, Accepted for ApJ

Hybrid Strangeon Stars

Authors: Chen Zhang, Yong Gao, Cheng-Jun Xia, Renxin Xu

Abstract: It was conjectured that the basic units of the ground state of bulk strong matter may be strange-clusters called strangeons, and they can form self-bound strangeon stars that are highly compact. Strangeon stars can develop a strange quark matter (SQM) core at high densities, particularly in the color-flavor-locking phase, yielding a branch of hybrid strangeon stars. We explore the stellar structur… ▽ More It was conjectured that the basic units of the ground state of bulk strong matter may be strange-clusters called strangeons, and they can form self-bound strangeon stars that are highly compact. Strangeon stars can develop a strange quark matter (SQM) core at high densities, particularly in the color-flavor-locking phase, yielding a branch of hybrid strangeon stars. We explore the stellar structure and astrophysical implications of hybrid strangeon stars. We find that hybrid strangeon stars can meet various astrophysical constraints on pulsar masses, radii, and tidal deformabilities. Finally, we show that the strangeon-SQM mixed phase is not preferred if the charge-neutrality condition is imposed at the strangeon-SQM transition region. △ Less

Submitted 8 January, 2024; v1 submitted 25 September, 2023; originally announced September 2023.

Comments: 11 pages, 5 figures. Published version

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

Quasi-periodic oscillations during magnetar giant flares in the strangeon star model

Authors: Hong-Bo Li, Yacheng Kang, Zexin Hu, Lijing Shao, Cheng-Jun Xia, Ren-Xin Xu

Abstract: Soft gamma-ray repeaters (SGRs) are widely understood as slowly rotating isolated neutron stars. Their generally large spin-down rates, high magnetic fields, and strong outburst energies render them different from ordinary pulsars. In a few giant flares (GFs) and short bursts of SGRs, high-confidence quasi-periodic oscillations (QPOs) were observed. Although remaining an open question, many theore… ▽ More Soft gamma-ray repeaters (SGRs) are widely understood as slowly rotating isolated neutron stars. Their generally large spin-down rates, high magnetic fields, and strong outburst energies render them different from ordinary pulsars. In a few giant flares (GFs) and short bursts of SGRs, high-confidence quasi-periodic oscillations (QPOs) were observed. Although remaining an open question, many theoretical studies suggest that the torsional oscillations caused by starquakes could explain QPOs. Motivated by this scenario, we systematically investigate torsional oscillation frequencies based on the strangeon-star (SS) model with various values of harmonic indices and overtones. To characterize the strong-repulsive interaction at short distances and the non-relativistic nature of strangeons, a phenomenological Lennard-Jones model is adopted. We show that, attributing to the large shear modulus of SSs, our results explain well the high-frequency QPOs ($\gtrsim 150\,\mathrm{Hz}$) during the GFs. The low-frequency QPOs ($\lesssim 150\,\mathrm{Hz}$) can also be interpreted when the ocean-crust interface modes are included. We also discuss possible effects of the magnetic field on the torsional mode frequencies. Considering realistic models with general-relativistic corrections and magnetic fields, we further calculate torsional oscillation frequencies for quark stars. We show that it would be difficult for quark stars to explain all QPOs in GFs. Our work advances the understanding of the nature of QPOs and magnetar asteroseismology. △ Less

Submitted 28 October, 2023; v1 submitted 18 September, 2023; originally announced September 2023.

Comments: 9 pages, 7 figures; accepted by MNRAS

Journal ref: MNRAS 527 (2024) 855

Data-constrained Magnetohydrodynamic Simulation of an Intermediate Solar Filament Eruption

Authors: Yang Guo, Jinhan Guo, Yiwei Ni, M. D. Ding, P. F. Chen, Chun Xia, Rony Keppens, Kai E. Yang

Abstract: Solar eruptive activities could occur in weak magnetic field environments and over large spatial scales, especially relevant to eruptions involving intermediate or quiescent solar filaments. To handle the large scales, we implement and apply a flux rope embedding method using regularized Biot-Savart laws in the spherical coordinate system. Combined with a potential field source surface model and a… ▽ More Solar eruptive activities could occur in weak magnetic field environments and over large spatial scales, especially relevant to eruptions involving intermediate or quiescent solar filaments. To handle the large scales, we implement and apply a flux rope embedding method using regularized Biot-Savart laws in the spherical coordinate system. Combined with a potential field source surface model and a magneto-frictional method, a nonlinear force-free field comprising a flux rope embedded in a potential field is constructed. Using the combined nonlinear force-free field as the initial condition, we then perform a zero-$β$ data-constrained magnetohydrodynamic (MHD) simulation for an M8.7 flare at 03:38 UT on 2012 January 23. The MHD model reproduces the eruption process, flare ribbon evolution (represented by the quasi-separatrix layer evolution) and kinematics of the flux rope. This approach could potentially model global-scale eruptions from weak field regions. △ Less

Submitted 3 September, 2023; originally announced September 2023.

Comments: 23 pages, 7 figures, accepted for publicaiton in ApJ

Diurnal modulation of electron recoils from DM-nucleon scattering through the Migdal effect

Authors: Mai Qiao, Chen Xia, Yu-Feng Zhou

Abstract: Halo dark matter (DM) particles could lose energy due to the scattering off nuclei within the Earth before reaching the underground detectors of DM direct detection experiments. This Earth shielding effect can result in diurnal modulation of the DM-induced recoil event rates observed underground due to the self-rotation of the Earth. For electron recoil signals from DM-electron scatterings, the cu… ▽ More Halo dark matter (DM) particles could lose energy due to the scattering off nuclei within the Earth before reaching the underground detectors of DM direct detection experiments. This Earth shielding effect can result in diurnal modulation of the DM-induced recoil event rates observed underground due to the self-rotation of the Earth. For electron recoil signals from DM-electron scatterings, the current experimental constraints are very stringent such that the diurnal modulation cannot be observed for halo DM. We propose a novel type of diurnal modulation effect: diurnal modulation in electron recoil signals induced by DM-nucleon scattering via the Migdal effect. We set so far the most stringent constraints on DM-nucleon scattering cross section via the Migdal effect for sub-GeV DM using the S2-only data of PandaX-II and PandaX-4T with improved simulations of the Earth shielding effect. Based on the updated constraints, we show that the Migdal effect induced diurnal modulation of electron events can still be significant in the low energy region, and can be probed by experiments such as PandaX-4T in the near future. △ Less

Submitted 1 November, 2023; v1 submitted 24 July, 2023; originally announced July 2023.

Comments: Comments on the Migdal effects added, figures and text improved, version accepted by JCAP

Rescaling strange-cluster stars and its implications on gravitational-wave echoes

Authors: Chen Zhang, Yong Gao, Cheng-Jun Xia, Renxin Xu

Abstract: Solid states of strange-cluster matter called strangeon matter can form strangeon stars that are highly compact. We show that strangeon matter and strangeon stars can be recast into dimensionless forms by a simple reparametrization and rescaling, through which we manage to maximally reduce the number of degrees of freedom. With this dimensionless scheme, we find that strangeon stars are generally… ▽ More Solid states of strange-cluster matter called strangeon matter can form strangeon stars that are highly compact. We show that strangeon matter and strangeon stars can be recast into dimensionless forms by a simple reparametrization and rescaling, through which we manage to maximally reduce the number of degrees of freedom. With this dimensionless scheme, we find that strangeon stars are generally compact enough to feature a photon sphere that is essential to foster gravitational-wave (GW) echoes. Rescaling the dimension back, we illustrate its implications on the expanded dimensional parameter space, and calculate the GW echo frequencies associated with strangeon stars, showing that the minimum echo frequency is $\sim 8$ kHz for empirical parameter space that satisfies the GW170817 constraint, and can reduce to $\mathcal O(100)$ Hertz at the extended limit. △ Less

Submitted 2 October, 2023; v1 submitted 16 May, 2023; originally announced May 2023.

Comments: 8 pages, 6 figures. Published version. Typo fixed

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

Ultra low-mass and small-radius white dwarfs made of heavy elements

Authors: Cheng-Jun Xia, Yong-Feng Huang, Hong-Bo Li, Lijing Shao, Ren-Xin Xu

Abstract: Seven ultra low-mass and small-radius white dwarfs (LSPM J0815+1633, LP 240-30, BD+20 5125B, LP 462-12, WD J1257+5428, 2MASS J13453297+4200437, and SDSS J085557.46+053524.5) have been recently identified with masses ranging from $\sim$0.02 $M_\odot$ to $\sim$0.08 $M_\odot$ and radii from $\sim$ 4270 km to 10670 km. The mass-radius measurements of these white dwarfs pose challenges to traditional w… ▽ More Seven ultra low-mass and small-radius white dwarfs (LSPM J0815+1633, LP 240-30, BD+20 5125B, LP 462-12, WD J1257+5428, 2MASS J13453297+4200437, and SDSS J085557.46+053524.5) have been recently identified with masses ranging from $\sim$0.02 $M_\odot$ to $\sim$0.08 $M_\odot$ and radii from $\sim$ 4270 km to 10670 km. The mass-radius measurements of these white dwarfs pose challenges to traditional white dwarf models assuming they are mostly made of nuclei lighter than $^{56}$Fe. In this work we consider the possibility that those white dwarfs are made of heavier elements. Due to the small charge-to-mass ratios in heavy elements, the electron number density in white dwarf matter is effectively reduced, which reduces the pressure with additional contributions of lattice energy and electron polarization corrections. This consequently leads to white dwarfs with much smaller masses and radii, which coincide with the seven ultra low-mass and small-radius white dwarfs. The corresponding equation of state and matter contents of dense stellar matter with and without reaching the cold-catalyzed ground state are presented, which are obtained using the latest Atomic Mass Evaluation (AME 2020). Further observations are necessary to unveil the actual matter contents in those white dwarfs via, e.g., spectroscopy, asteroseismology, and discoveries of other ultra low-mass and small-radius white dwarfs. △ Less

Submitted 9 May, 2023; originally announced May 2023.

Journal ref: Front. Astron. Space Sci. 10 (2023) 1334642

Coronal Heating as Determined by the Solar Flare Frequency Distribution Obtained by Aggregating Case Studies

Authors: James Paul Mason, Alexandra Werth, Colin G. West, Allison A. Youngblood, Donald L. Woodraska, Courtney Peck, Kevin Lacjak, Florian G. Frick, Moutamen Gabir, Reema A. Alsinan, Thomas Jacobsen, Mohammad Alrubaie, Kayla M. Chizmar, Benjamin P. Lau, Lizbeth Montoya Dominguez, David Price, Dylan R. Butler, Connor J. Biron, Nikita Feoktistov, Kai Dewey, N. E. Loomis, Michal Bodzianowski, Connor Kuybus, Henry Dietrick, Aubrey M. Wolfe , et al. (977 additional authors not shown)

Abstract: Flare frequency distributions represent a key approach to addressing one of the largest problems in solar and stellar physics: determining the mechanism that counter-intuitively heats coronae to temperatures that are orders of magnitude hotter than the corresponding photospheres. It is widely accepted that the magnetic field is responsible for the heating, but there are two competing mechanisms th… ▽ More Flare frequency distributions represent a key approach to addressing one of the largest problems in solar and stellar physics: determining the mechanism that counter-intuitively heats coronae to temperatures that are orders of magnitude hotter than the corresponding photospheres. It is widely accepted that the magnetic field is responsible for the heating, but there are two competing mechanisms that could explain it: nanoflares or Alfvén waves. To date, neither can be directly observed. Nanoflares are, by definition, extremely small, but their aggregate energy release could represent a substantial heating mechanism, presuming they are sufficiently abundant. One way to test this presumption is via the flare frequency distribution, which describes how often flares of various energies occur. If the slope of the power law fitting the flare frequency distribution is above a critical threshold, $α=2$ as established in prior literature, then there should be a sufficient abundance of nanoflares to explain coronal heating. We performed $>$600 case studies of solar flares, made possible by an unprecedented number of data analysts via three semesters of an undergraduate physics laboratory course. This allowed us to include two crucial, but nontrivial, analysis methods: pre-flare baseline subtraction and computation of the flare energy, which requires determining flare start and stop times. We aggregated the results of these analyses into a statistical study to determine that $α= 1.63 \pm 0.03$. This is below the critical threshold, suggesting that Alfvén waves are an important driver of coronal heating. △ Less

Submitted 9 May, 2023; originally announced May 2023.

Comments: 1,002 authors, 14 pages, 4 figures, 3 tables, published by The Astrophysical Journal on 2023-05-09, volume 948, page 71

Thermodynamic and Magnetic Topology Evolution of the X1.0 Flare on 2021 October 28 Simulated by a Data-driven Radiative Magnetohydrodynamic Model

Authors: Jin-han Guo, Yi-wei Ni, Ze Zhong, Yang Guo, Chun Xia, Hai-tang Li, Stefaan Poedts, Brigitte Schmieder, Peng-fei Chen

Abstract: Solar filament eruptions, flares and coronal mass ejections (CMEs) are manifestations of drastic release of energy in the magnetic field, which are related to many eruptive phenomena from the Earth magnetosphere to black hole accretion disks. With the availability of high-resolution magnetograms on the solar surface, observational data-based modelling is a promising way to quantitatively study the… ▽ More Solar filament eruptions, flares and coronal mass ejections (CMEs) are manifestations of drastic release of energy in the magnetic field, which are related to many eruptive phenomena from the Earth magnetosphere to black hole accretion disks. With the availability of high-resolution magnetograms on the solar surface, observational data-based modelling is a promising way to quantitatively study the underlying physical mechanisms behind observations. By incorporating thermal conduction and radiation losses in the energy equation, we develop a new data-driven radiative magnetohydrodynamic (MHD) model, which has the capability to capture the thermodynamic evolution compared to our previous zero-\b{eta} model. Our numerical results reproduce major observational characteristics of the X1.0 flare on 2021 October 28 in NOAA active region (AR) 12887, including the morphology of the eruption, kinematic of flare ribbons, extreme-ultraviolet (EUV) radiations, and two components of the EUV waves predicted by the magnetic stretching model, i.e., a fast-mode shock wave and a slower apparent wave due to successive stretching of magnetic field lines. Moreover, some intriguing phenomena are revealed in the simulation. We find that flare ribbons separate initially and ultimately stop at the outer stationary quasi-separatrix layers (QSLs). Such outer QSLs correspond to the border of the filament channel and determine the final positions of flare ribbons, which can be used to predict the size and the lifetime of a flare before it occurs. In addition, the side view of the synthesized EUV and white-light images exhibit typical three-part structures of CMEs, where the bright leading front is roughly cospatial with the non-wave component of the EUV wave, reinforcing the magnetic stretching model for the slow component of EUV waves. △ Less

Submitted 24 March, 2023; originally announced March 2023.

Comments: 46 pages, 15 figures, accepted for publication in The Astrophysical Journal Supplement

MPI-AMRVAC 3.0: updates to an open-source simulation framework

Authors: Rony Keppens, Beatrice Popescu Braileanu, Yuhao Zhou, Wenzhi Ruan, Chun Xia, Yang Guo, Niels Claes, Fabio Bacchini

Abstract: Computational astrophysics routinely combines grid-adaptive capabilities with modern shock-capturing, high resolution spatio-temporal schemes on multi-dimensional hydro- and magnetohydrodynamics. We provide an update on developments within the open-source MPI-AMRVAC code. With online documentation, the MPI-AMRVAC 3.0 release includes several added equation sets, and many options to explore and qua… ▽ More Computational astrophysics routinely combines grid-adaptive capabilities with modern shock-capturing, high resolution spatio-temporal schemes on multi-dimensional hydro- and magnetohydrodynamics. We provide an update on developments within the open-source MPI-AMRVAC code. With online documentation, the MPI-AMRVAC 3.0 release includes several added equation sets, and many options to explore and quantify the influence of implementation details. Showcasing this on a variety of hydro and MHD tests, we document new modules of interest for state-of-the-art solar applications. Test cases address how higher order reconstructions impact long term simulations of shear layers, with and without gas-dust coupling, how runaway radiative losses transit to intricate multi-temperature, multi-phase dynamics, and how different flavors of spatio-temporal schemes and magnetic monopole control produce consistent MHD results in combination with adaptive meshes. We demonstrate Super-Time-Stepping strategies for specific parabolic terms and give details on all implemented Implicit-Explicit integrators. A new magnetofrictional module can be used for computing force-free magnetic fields or for data-driven time-dependent evolutions, while the Regularized-Biot-Savart-Law approach can insert fluxropes in 3D domains. Synthetic observations of 3D MHD simulations can be rendered on-the-fly, or in post-processing, in many spectral wavebands. A particle module and a generic fieldline tracing, compatible with the hierarchical meshes, can be used to sample information at prescribed locations, to follow dynamics of charged particles, or realize two-way coupled simulations between MHD setups and field-aligned non-thermal processes. Highlighting the latest additions and various technical aspects, our open-source strategy welcomes any further code usage, contribution, or spin-off development. △ Less

Submitted 6 March, 2023; originally announced March 2023.

Comments: Accepted for A&A

Journal ref: A&A 673, A66 (2023)

Probing phase transition in neutron stars via the crust-core interfacial mode

Authors: Jiaxiang Zhu, Chuming Wang, Chengjun Xia, Enping Zhou, Yiqiu Ma

Abstract: Gravitational waves emitted from the binary neutron star (BNS) systems can carry information about the dense matter phase in these compact stars. The crust-core interfacial mode is an oscillation mode in a neutron star and it depends mostly on the equation of the state of the matter in the crust-core transition region. This mode can be resonantly excited by the tidal field of an inspiraling-in BNS… ▽ More Gravitational waves emitted from the binary neutron star (BNS) systems can carry information about the dense matter phase in these compact stars. The crust-core interfacial mode is an oscillation mode in a neutron star and it depends mostly on the equation of the state of the matter in the crust-core transition region. This mode can be resonantly excited by the tidal field of an inspiraling-in BNS system, thereby affecting the emitted gravitational waves, and hence could be used to probe the equation of state in the crust-core transition region. In this work, we investigate in detail how the first-order phase transition inside the neutron star affects the properties of the crust-core interfacial mode, using a Newtonian fluid perturbation theory on a general relativistic background solution of the stellar structure. Two possible types of phase transitions are considered: (1) the phase transitions happen in the fluid core but near the crust-core interface, which results in density discontinuities; and (2) the strong interaction phase transitions in the dense core (as in the conventional hybrid star case). These phase transitions' impacts on interfacial mode properties are discussed. In particular, the former phase transition has a minor effect on the M-R relation and the adiabatic tidal deformability, but can significantly affect the interfacial mode frequency and thereby could be probed using gravitational waves. For the BNS systems, we discuss the possible observational signatures of these phase transitions in the gravitational waveforms and their detectability. Our work enriches the exploration of the physical properties of the crust-core interfacial mode and provides a promising method for probing the phase transition using the seismology of a compact star. △ Less

Submitted 21 November, 2022; originally announced November 2022.

Comments: 18 pages, 14 figures

Bulk Strong Matter: the Trinity

Authors: Xiaoyu Lai, Chengjun Xia, Renxin Xu

Abstract: Our world is wonderful because of the normal but negligibly small baryonic part (i.e., atoms) although unknown dark matter and dark energy dominate the Universe. A stable atomic nucleus could be simply termed as ``strong matter'' since its nature is dominated by the fundamental strong interaction. Is there any other form of strong matter? Although nuclei are composed of 2-flavoured (i.e., up and d… ▽ More Our world is wonderful because of the normal but negligibly small baryonic part (i.e., atoms) although unknown dark matter and dark energy dominate the Universe. A stable atomic nucleus could be simply termed as ``strong matter'' since its nature is dominated by the fundamental strong interaction. Is there any other form of strong matter? Although nuclei are composed of 2-flavoured (i.e., up and down flavours of valence quarks) nucleons, it is conjectured that bulk strong matter could be 3-flavoured (with additional strange quarks) if the baryon number exceeds the critical value, $A_{\rm c}$, in which case quarks could be either free (so-called strange quark matter) or localized (in strangeons, coined by combining ``strange nucleon''). Bulk strong matter could be manifested in the form of compact stars, cosmic rays, and even dark matter. This trinity will be explained in this brief review, that may impact dramatically on today's physics, particularly in the era of multi-messenger astronomy after the discovery of gravitational wave. △ Less

Submitted 6 November, 2022; v1 submitted 4 October, 2022; originally announced October 2022.

Comments: 12 pages, 2 figures. Accepted by Advances in Physics: X

Journal ref: ADVANCES IN PHYSICS: X (2023) VOL. 8, NO. 1, 2137433

Elastic properties of nuclear pasta in a fully three-dimensional geometry

Authors: Cheng-Jun Xia, Toshiki Maruyama, Nobutoshi Yasutake, Toshitaka Tatsumi, Ying-Xun Zhang

Abstract: Realistic estimations on the elastic properties of neutron star matter are carried out with a large strain ($\varepsilon \lesssim 0.5$) in the framework of relativistic-mean-field model with Thomas-Fermi approximation, where various crystalline configurations are considered in a fully three-dimensional geometry with reflection symmetry. Our calculation confirms the validity of assuming Coulomb cry… ▽ More Realistic estimations on the elastic properties of neutron star matter are carried out with a large strain ($\varepsilon \lesssim 0.5$) in the framework of relativistic-mean-field model with Thomas-Fermi approximation, where various crystalline configurations are considered in a fully three-dimensional geometry with reflection symmetry. Our calculation confirms the validity of assuming Coulomb crystals for the droplet phase above neutron drip density, which nonetheless does not work at large densities since the elastic constants are found to be decreasing after reaching their peaks. Similarly, the analytic formulae derived in the incompressible liquid-drop model give excellent description for the rod phase at small densities, which overestimates the elastic constants at larger densities. For slabs, due to the negligence on the variations of their thicknesses, the analytic formulae from liquid-drop model agree qualitatively but not quantitatively with our numerical estimations. By fitting to the numerical results, these analytic formulae are improved by introducing dampening factors. The impacts of nuclear symmetry energy are examined adopting two parameter sets, corresponding to the slope of symmetry energy $L = 41.34$ and 89.39 MeV. Even with the uncertainties caused by the anisotropy in polycrystallines, the elastic properties of neutron star matter obtained with $L = 41.34$ and 89.39 MeV are distinctively different, results in detectable differences in various neutron star activities. △ Less

Submitted 14 February, 2023; v1 submitted 27 September, 2022; originally announced September 2022.

Unified neutron star EOSs and neutron star structures in RMF models

Authors: Cheng-Jun Xia, Toshiki Maruyama, Ang Li, Bao Yuan Sun, Wen-Hui Long, Ying-Xun Zhang

Abstract: In the framework of Thomas-Fermi approximation, we study systematically the EOSs and microscopic structures of neutron star matter in a vast density range with $n_\mathrm{b}\approx 10^{-10}$-2 $\mathrm{fm}^{-3}$, where various covariant density functionals are adopted, i.e., those with nonlinear self couplings (NL3, PK1, TM1, GM1, MTVTC) and density-dependent couplings (DD-LZ1, DDME-X, PKDD, DD-ME… ▽ More In the framework of Thomas-Fermi approximation, we study systematically the EOSs and microscopic structures of neutron star matter in a vast density range with $n_\mathrm{b}\approx 10^{-10}$-2 $\mathrm{fm}^{-3}$, where various covariant density functionals are adopted, i.e., those with nonlinear self couplings (NL3, PK1, TM1, GM1, MTVTC) and density-dependent couplings (DD-LZ1, DDME-X, PKDD, DD-ME2, DD2, TW99). It is found that the EOSs generally coincide with each other at $n_\mathrm{b}\lesssim 10^{-4}$ fm${}^{-3}$ and 0.1 fm${}^{-3}\lesssim n_\mathrm{b} \lesssim 0.3$ fm${}^{-3}$, while in other density regions they are sensitive to the effective interactions between nucleons. By adopting functionals with larger slope of symmetry energy $L$, the curvature parameter $K_\mathrm{sym}$ and neutron drip density generally increase, while the droplet size, proton number of nucleus, core-crust transition density, and onset density of non-spherical nuclei decrease. All functionals predict neutron stars with maximum masses exceeding the two-solar-mass limit, while those of DD2, DD-LZ1, DD-ME2, and DDME-X predict optimum neutron star radii according to the observational constraints. Nevertheless, the corresponding skewness coefficients $J$ are much lager than expected, while only the functionals MTVTC and TW99 meet the start-of-art constraints on $J$. More accurate measurements on the radius of PSR J0740+6620 and the maximum mass of neutron stars are thus essential to identify the functional that satisfies all constraints from nuclear physics and astrophysical observations. Approximate linear correlations between neutron stars' radii at $M=1.4 M_{\odot}$ and $2 M_{\odot}$, the slope $L$ and curvature parameter $K_\mathrm{sym}$ of symmetry energy are observed as well, which is mainly attributed to the curvature-slope correlations in the functionals adopted here. △ Less

Submitted 26 August, 2022; originally announced August 2022.

Comments: The EOS tables can be found at: teacher.yzu.edu.cn/CJXia/en/lwcg/492393/content/79242.htm

Journal ref: Commun. Theor. Phys. 74 (2022) 095303

Transitions and Origin of the Type-B Quasi-Periodic Oscillation in the Black Hole X-ray Binary MAXI~ J1348--630

Authors: H. X. Liu, Y. Huang, Q. C. Bu, W. Yu, Z. X. Yang, L. Zhang, L. D. Kong, G. C. Xiao, J. L. Qu, S. N. Zhang, S. Zhang, L. M. Song, S. M. Jia, X. Ma, L. Tao, M. Y. Ge, Q. Z. Liu, J. Z. Yan, R. C. Ma, X. Q. Ren, D. K. Zhou, T. M. Li, B. Y. Wu, Y. C. Xu, Y. F. Du , et al. (4 additional authors not shown)

Abstract: The fast transitions between different types of quasi-periodic oscillations (QPOs) are generally observed in black hole transient sources (BHTs). We present a detailed study on the timing and spectral properties of the transitions of type-B QPOs in MAXI~J1348--630, observed by \emph{Insight}-HXMT. The fractional rms variability--energy relationship and energy spectra reveal that type-B QPOs probab… ▽ More The fast transitions between different types of quasi-periodic oscillations (QPOs) are generally observed in black hole transient sources (BHTs). We present a detailed study on the timing and spectral properties of the transitions of type-B QPOs in MAXI~J1348--630, observed by \emph{Insight}-HXMT. The fractional rms variability--energy relationship and energy spectra reveal that type-B QPOs probably originate from jet precession. Compared to weak power-law dominated power spectrum, when type-B QPO is present, the corresponding energy spectrum shows an increase in Comptonization component and the need for {\tt\string xillverCp} component, and a slight increase of height of the corona when using {\tt\string relxilllp} model. Therefore, we suggest that a coupled inner disk-jet region is responsible for the observed type-B QPOs transitions. The time scale for the appearance/disappearance of type-B QPOs is either long or short (seconds), which may indicate an instability of disk-jet structure. For these phenomena, we give the hypothesis that the Bardeen-Petterson effect causes disk-jet structure to align with BH spin axis, or that the disappearance of small-scale jets bound by the magnetic flux tubes lead to the disappearance of type-B QPOs. We observed three events regarding the B/C transitions, one of which occurred in a short time from $\sim 9.2$ Hz (C) to $\sim 4.8$ Hz (B). The energy spectral analysis for the other two transitions shows that when type-C QPO is present, the Comptonization flux is higher, the spectrum is harder and the inner radius of disk changes insignificantly. We suggest that type-C QPOs probably originate from relatively stronger jets or corona. △ Less

Submitted 15 August, 2022; originally announced August 2022.

Comments: 26 pages, 9 figures

Formation of Quiescent Prominence Magnetic Field by Supergranulations

Authors: Qingjun Liu, Chun Xia

Abstract: To understand the formation of quiescent solar prominences, the origin of their magnetic field structures, i.e., magnetic flux ropes (MFRs), must be revealed. We use three-dimensional magnetofriction simulations in a spherical subdomain to investigate the role of typical supergranular motions in the long-term formation of a prominence magnetic field. Time-dependant horizontal supergranular motions… ▽ More To understand the formation of quiescent solar prominences, the origin of their magnetic field structures, i.e., magnetic flux ropes (MFRs), must be revealed. We use three-dimensional magnetofriction simulations in a spherical subdomain to investigate the role of typical supergranular motions in the long-term formation of a prominence magnetic field. Time-dependant horizontal supergranular motions with and without the effect of Coriolis force are simulated on the solar surface via Voronoi tessellation. The vortical motions by the Coriolis effect at boundaries of supergranules inject magnetic helicity into the corona. The helicity is transferred and accumulated along the polarity inversion line (PIL) as strongly sheared magnetic field via helicity condensation. The diverging motions of supergranules converge opposite magnetic polarities at the PIL and drive the magnetic reconnection between footpoints of the sheared magnetic arcades to form an MFR. The magnetic network, negative-helicity MFR in the northern hemisphere, and fragmented-to-continuous formation process of magnetic dip regions are in agreement with observations. Although diverging supergranulations, differential rotation, and meridional flows are included, the simulation without the Coriolis effect can not produce an MFR or sheared arcades to host a prominence. Therefore Coriolis force is a key factor for helicity injection and the formation of magnetic structures of quiescent solar prominences. △ Less

Submitted 13 July, 2022; originally announced July 2022.

Comments: 5 figures

Nuclear pasta structures at high temperatures

Authors: Cheng-Jun Xia, Toshiki Maruyama, Nobutoshi Yasutake, Toshitaka Tatsumi

Abstract: We investigate nuclear pasta structures at high temperatures in the framework of relativistic mean field model with Thomas-Fermi approximation. Typical pasta structures (droplet, rod, slab, tube, and bubble) are obtained, which form various crystalline configurations. The properties of those nuclear pastas are examined in a three-dimensional geometry with reflection symmetry, where the optimum lat… ▽ More We investigate nuclear pasta structures at high temperatures in the framework of relativistic mean field model with Thomas-Fermi approximation. Typical pasta structures (droplet, rod, slab, tube, and bubble) are obtained, which form various crystalline configurations. The properties of those nuclear pastas are examined in a three-dimensional geometry with reflection symmetry, where the optimum lattice constants are fixed by reproducing the droplet/bubble density that minimizes the free energy adopting spherical or cylindrical approximations for Wigner-Seitz cells. It is found that different crystalline structures can evolve into each other via volume conserving deformations. For fixed densities and temperatures, the differences of the free energies per baryon of nuclear pasta in various shapes and lattice structures are typically on the order of tens of keV, suggesting the possible coexistence of those structures. As temperature increases, the thermodynamic fluctuations are expected to disrupt the long-range ordering in nuclear pasta structures. We then estimate the critical conditions for nuclear pasta to become disordered and behave like liquid, which are found to be sensitive to the densities, temperatures, proton fractions, and nuclear shapes. If we further increase temperature, eventually the nonuniform structures of nuclear pasta become unstable and are converted into uniform nuclear matter. The phase diagrams of nuclear matter are then estimated, which should be useful for understanding the evolutions of neutron stars, supernova dynamics, and binary neutron star mergers. △ Less

Submitted 27 September, 2022; v1 submitted 19 June, 2022; originally announced June 2022.

Journal ref: Phys. Rev. D 106 (2022) 063020

Interface effects of quark matter: Light-quark nuggets and compact stars

Authors: Cheng-Jun Xia, Jian-Feng Xu, Guang-Xiong Peng, Ren-Xin Xu

Abstract: The interface effects of quark matter play important roles in the properties of compact stars and small nuggets such as strangelets and $ud$QM nuggets. By introducing a density derivative term to the Lagrangian density and adopting Thomas-Fermi approximation, we find it is possible to reproduce the results obtained by solving Dirac equations. Adopting certain parameter sets, the energy per baryon… ▽ More The interface effects of quark matter play important roles in the properties of compact stars and small nuggets such as strangelets and $ud$QM nuggets. By introducing a density derivative term to the Lagrangian density and adopting Thomas-Fermi approximation, we find it is possible to reproduce the results obtained by solving Dirac equations. Adopting certain parameter sets, the energy per baryon of $ud$QM nuggets decreases with baryon number $A$ and become more stable than nuclei at $A\gtrsim 300$. The effects of quark matter symmetry energy are examined, where $ud$QM nuggets at $A\approx 1000$ can be more stable than others if large symmetry energy is adopted. In such cases, larger $ud$QM nuggets will decay via fission and the surface of an $ud$QM star will fragment into a crust made of $ud$QM nuggets and electrons, which resembles the cases of a strange star's crust. The corresponding microscopic structures are then investigated adopting spherical and cylindrical approximations for the Wigner-Seitz cells, where the droplet phase is found to be the most stable configuration with $ud$QM stars' crusts and $ud$QM dwarfs made of $ud$QM nuggets ($A\approx 1000$) and electrons. For the cases considered here, the crust thickness of $ud$QM stars is typically $\sim$200 m, which reaches a few kilometers if we neglect the interface effects and adopt Gibbs construction. The masses and radii of $ud$QM dwarfs are smaller than typical white dwarfs, which would increase if the interface effects are neglected. △ Less

Submitted 26 August, 2022; v1 submitted 21 May, 2022; originally announced May 2022.

Journal ref: Phys. Rev. D 106 (2022) 034016

A formation mechanism for the large plumes in the prominence

Authors: Jincheng Wang, Xiaoli Yan, Zhike Xue, Liheng Yang, Qiaoling Li, Hechao Chen, Chun Xia, Zhong Liu

Abstract: To understand the formation mechanism of large plumes in solar prominences, we investigate the formation process of two such phenomena. We studied the dynamic and thermal properties of two large plumes using observations from New Vacuum Solar Telescope, the Solar Dynamic Observatory, and the Solar Terrestrial Relations Observatory-Ahead. We find that two large plumes observed with high-resolution… ▽ More To understand the formation mechanism of large plumes in solar prominences, we investigate the formation process of two such phenomena. We studied the dynamic and thermal properties of two large plumes using observations from New Vacuum Solar Telescope, the Solar Dynamic Observatory, and the Solar Terrestrial Relations Observatory-Ahead. We find that two large plumes observed with high-resolution data are quite different from previously studied small-scale plumes. They are born at the top of a prominence bubble with a large projected area of 10-20 Mm^2 . Before the occurrence of each large plume, the bubble expands and takes on a quasi-semicircular appearance. Meanwhile, the emission intensity of extreme-ultra-violet (EUV) bands increases in the bubble. A small-scale filament is found to erupt in the bubble during the second large plume. At the point at which the height of the bubble is comparable with half the width of the bubble, the bubble becomes unstable and generates the plumes. During the formation of plumes, two side edges of the top of the bubble, which are dominated by opposite Doppler signals, approach each other. The large plume then emerges and keeps rising up with a constant speed of about 13-15 km/s. These two large plumes have temperatures of 1.3 x 10^6 Kelvin and densities of 2.0 x 10^9 cm^-3, two orders hotter and one order less dense than the typical prominence. We also find that the bubble is a hot, low-density volume instead of a void region beneath the cold and dense prominence. Therefore, we conclude that these two large plumes are the result of the breakup of the prominence bubble triggered by an enhancement of thermal pressure; they separate from the bubble, most likely by magnetic reconnection. △ Less

Submitted 17 February, 2022; originally announced February 2022.

Comments: 14 pages, 7 figures, accepted for publication in A&A

Peculiar disk behaviors of the black hole candidate MAXI J1348-630 in the hard state observed by Insight-HXMT and Swift

Authors: W. Zhang, L. Tao, R. Soria, J. L. Qu, S. N. Zhang, S. S. Weng, L. zhang, Y. N. Wang, Y. Huang, R. C. Ma, S. Zhang, M. Y. Ge, L. M. Song, X. Ma, Q. C. Bu, C. Cai, X. L. Cao, Z. Chang, L. Chen, T. X. Chen, Y. B. Chen, Y. Chen, Y. P. Chen, W. W. Cui, Y. Y. Du , et al. (72 additional authors not shown)

Abstract: We present a spectral study of the black hole candidate MAXI J1348-630 during its 2019 outburst, based on monitoring observations with Insight-HXMT and Swift. Throughout the outburst, the spectra are well fitted with power-law plus disk-blackbody components. In the soft-intermediate and soft states, we observed the canonical relation L ~ T_in^4 between disk luminosity L and peak colour temperature… ▽ More We present a spectral study of the black hole candidate MAXI J1348-630 during its 2019 outburst, based on monitoring observations with Insight-HXMT and Swift. Throughout the outburst, the spectra are well fitted with power-law plus disk-blackbody components. In the soft-intermediate and soft states, we observed the canonical relation L ~ T_in^4 between disk luminosity L and peak colour temperature T_in, with a constant inner radius R_in (traditionally identified with the innermost stable circular orbit). At other stages of the outburst cycle, the behaviour is more unusual, inconsistent with the canonical outburst evolution of black hole transients. In particular, during the hard rise, the apparent inner radius is smaller than in the soft state (and increasing), and the peak colour temperature is higher (and decreasing). This anomalous behaviour is found even when we model the spectra with self-consistent Comptonization models, which take into account the up-scattering of photons from the disk component into the power-law component. To explain both those anomalous trends at the same time, we suggest that the hardening factor for the inner disk emission was larger than the canonical value of ~1.7 at the beginning of the outburst. A more physical trend of radii and temperature evolution requires a hardening factor evolving from ~3.5 at the beginning of the hard state to ~1.7 in the hard intermediate state. This could be evidence that the inner disk was in the process of condensing from the hot, optically thin medium and had not yet reached a sufficiently high optical depth for its emission spectrum to be described by the standard optically-thick disk solution. △ Less

Submitted 27 January, 2022; originally announced January 2022.

Comments: 20 pages, 10 figures, 4 tables, accepted by ApJ

Where have all the interstellar silicon carbides gone?

Authors: Tao Chen, C. Y. Xiao, Aigen Li, C. T. Zhou

Abstract: The detection of the 11.3-micron emission feature characteristic of the Si--C stretch in carbon-rich evolved stars reveals that silicon carbide (SiC) dust grains are condensed in the outflows of carbon stars. SiC dust could be a significant constituent of interstellar dust since it is generally believed that carbon stars inject a considerable amount of dust into the interstellar medium (ISM). The… ▽ More The detection of the 11.3-micron emission feature characteristic of the Si--C stretch in carbon-rich evolved stars reveals that silicon carbide (SiC) dust grains are condensed in the outflows of carbon stars. SiC dust could be a significant constituent of interstellar dust since it is generally believed that carbon stars inject a considerable amount of dust into the interstellar medium (ISM). The presence of SiC dust in the ISM is also supported by the identification of presolar SiC grains of stellar origin in primitive meteorites. However, the 11.3-micron absorption feature of SiC has never been seen in the ISM and oxidative destruction of SiC is often invoked. In this work we quantitatively explore the destruction of interstellar SiC dust through oxidation based on molecular dynamics simulations and density functional theory calculations. We find that the reaction of an oxygen atom with SiC molecules and clusters is exothermic and could cause CO-loss. Nevertheless, even if this is extrapolable to bulk SiC dust, the destruction rate of SiC dust through oxidation could still be considerably smaller than the (currently believed) injection rate from carbon stars. Therefore, the lack of the 11.3-micron absorption feature of SiC dust in the ISM remains a mystery. A possible solution may lie in the currently believed stellar injection rate of SiC (which may have been overestimated) and/or the size of SiC dust (which may actually be considerably smaller than submicron in size). △ Less

Submitted 6 November, 2021; originally announced November 2021.

Comments: 7 pages, 5 figures; accepted for publication in MNRAS

Search for Gamma-Ray Bursts and Gravitational Wave Electromagnetic Counterparts with High Energy X-ray Telescope of \textit{Insight}-HXMT

Authors: C. Cai, S. L. Xiong, C. K. Li, C. Z. Liu, S. N. Zhang, X. B. Li, L. M. Song, B. Li, S. Xiao, Q. B. Yi, Y. Zhu, Y. G. Zheng, W. Chen, Q. Luo, Y. Huang, X. Y. Song, H. S. Zhao, Y. Zhao, Z. Zhang, Q. C. Bu, X. L. Cao, Z. Chang, L. Chen, T. X. Chen, Y. B. Chen , et al. (74 additional authors not shown)

Abstract: The High Energy X-ray telescope (HE) on-board the Hard X-ray Modulation Telescope (\textit{Insight}-HXMT) can serve as a wide Field of View (FOV) gamma-ray monitor with high time resolution ($μ$s) and large effective area (up to thousands cm$^2$). We developed a pipeline to search for Gamma-Ray Bursts (GRBs), using the traditional signal-to-noise ratio (SNR) method for blind search and the coheren… ▽ More The High Energy X-ray telescope (HE) on-board the Hard X-ray Modulation Telescope (\textit{Insight}-HXMT) can serve as a wide Field of View (FOV) gamma-ray monitor with high time resolution ($μ$s) and large effective area (up to thousands cm$^2$). We developed a pipeline to search for Gamma-Ray Bursts (GRBs), using the traditional signal-to-noise ratio (SNR) method for blind search and the coherent search method for targeted search. By taking into account the location and spectrum of the burst and the detector response, the targeted coherent search is more powerful to unveil weak and sub-threshold bursts, especially those in temporal coincidence with Gravitational Wave (GW) events. Based on the original method in literature, we further improved the coherent search to filter out false triggers caused by spikes in light curves, which are commonly seen in gamma-ray instruments (e.g. \textit{Fermi}/GBM, \textit{POLAR}). We show that our improved targeted coherent search method could eliminate almost all false triggers caused by spikes. Based on the first two years of \textit{Insight}-HXMT/HE data, our targeted search recovered 40 GRBs, which were detected by either \textit{Swift}/BAT or \textit{Fermi}/GBM but too weak to be found in our blind search. With this coherent search pipeline, the GRB detection sensitivity of \textit{Insight}-HXMT/HE is increased to about 1.5E-08 erg/cm$^2$ (200 keV$-$3 MeV). We also used this targeted coherent method to search \textit{Insight}-HXMT/HE data for electromagnetic (EM) counterparts of LIGO-Virgo GW events (including O2 and O3a runs). However, we did not find any significant burst associated with GW events. △ Less

Submitted 25 September, 2021; originally announced September 2021.

Comments: 12 pages, 14 figures, 5 tables; accepted for publication in MNRAS

Observations of an Electron-cold Ion Component Reconnection at the Edge of an Ion-scale Antiparallel Reconnection at the Dayside Magnetopause

Authors: S. Q. Zhao, H. Zhang, Terry Z. Liu, Huirong Yan, C. J. Xiao, Mingzhe Liu, Q. -G. Zong, Xiaogang Wang, Mijie Shi, Shangchun Teng, Huizi Wang, R. Rankin, C. Pollock, G. Le

Abstract: Solar wind parameters play a dominant role in reconnection rate, which controls the solar wind-magnetosphere coupling efficiency at Earth's magnetopause. Besides, low-energy ions from the ionosphere, frequently detected on the magnetospheric side of the magnetopause, also affect magnetic reconnection. However, the specific role of low-energy ions in reconnection is still an open question under act… ▽ More Solar wind parameters play a dominant role in reconnection rate, which controls the solar wind-magnetosphere coupling efficiency at Earth's magnetopause. Besides, low-energy ions from the ionosphere, frequently detected on the magnetospheric side of the magnetopause, also affect magnetic reconnection. However, the specific role of low-energy ions in reconnection is still an open question under active discussion. In the present work, we report in situ observations of a multiscale, multi-type magnetopause reconnection in the presence of low-energy ions using NASA's Magnetospheric Multiscale data on 11 September 2015. This study divides ions into cold and hot populations. The observations can be interpreted as a secondary reconnection dominated by electrons and cold ions located at the edge of an ion-scale reconnection. This analysis demonstrates a dominant role of cold ions in the secondary reconnection without hot ions' response. Cold ions and electrons are accelerated and heated by the secondary process. The case study provides observational evidence for the simultaneous operation of antiparallel and component reconnection. Our results imply that the pre-accelerated and heated cold ions and electrons in the secondary reconnection may participate in the primary ion-scale reconnection affecting the solar wind-magnetopause coupling and the complicated magnetic field topology affect the reconnection rate. △ Less

Submitted 22 September, 2021; originally announced September 2021.

Optimizing the hybrid parallelization of BHAC

Authors: Salvatore Cielo, Oliver Porth, Luigi Iapichino, Anupam Karmakar, Hector Olivares, Chun Xia

Abstract: We present our experience with the modernization on the GR-MHD code BHAC, aimed at improving its novel hybrid (MPI+OpenMP) parallelization scheme. In doing so, we showcase the use of performance profiling tools usable on x86 (Intel-based) architectures. Our performance characterization and threading analysis provided guidance in improving the concurrency and thus the efficiency of the OpenMP paral… ▽ More We present our experience with the modernization on the GR-MHD code BHAC, aimed at improving its novel hybrid (MPI+OpenMP) parallelization scheme. In doing so, we showcase the use of performance profiling tools usable on x86 (Intel-based) architectures. Our performance characterization and threading analysis provided guidance in improving the concurrency and thus the efficiency of the OpenMP parallel regions. We assess scaling and communication patterns in order to identify and alleviate MPI bottlenecks, with both runtime switches and precise code interventions. The performance of optimized version of BHAC improved by $\sim28\%$, making it viable for scaling on several hundreds of supercomputer nodes. We finally test whether porting such optimizations to different hardware is likewise beneficial on the new architecture by running on ARM A64FX vector nodes. △ Less

Submitted 27 August, 2021; originally announced August 2021.

Comments: 10 pages, 9 figures, 1 table; in review

A case study of ACV variables discovered in the Zwicky Transient Facility survey

Authors: N. Faltová, K. Kallová, M. Prišegen, P. Staněk, J. Supíková, C. Xia, K. Bernhard, S. Hümmerich, E. Paunzen

Abstract: Magnetic chemically peculiar (mCP) stars exhibit complex atmospheres that allow the investigation of the interplay of atomic diffusion, magnetic fields, and stellar rotation. A non-uniform surface distribution of chemical elements and the non-alignment of the rotational and magnetic axes result in the variability of several observables. Photometrically variable mCP stars are referred to as alpha2… ▽ More Magnetic chemically peculiar (mCP) stars exhibit complex atmospheres that allow the investigation of the interplay of atomic diffusion, magnetic fields, and stellar rotation. A non-uniform surface distribution of chemical elements and the non-alignment of the rotational and magnetic axes result in the variability of several observables. Photometrically variable mCP stars are referred to as alpha2 Canum Venaticorum (ACV) variables. The present work presents a case study of known variables from the Zwicky Transient Facility (ZTF) survey, with the aim of investigating the survey's suitability for the detection and study of new ACV variables. Using suitable selection criteria based on the known characteristics of ACV variables, candidate ACV stars were selected from the ZTF Catalog of Periodic Variable Stars. All light curves were inspected in detail to select the most promising candidates. Where available, low-resolution spectra from the LAMOST were employed to classify the stars on the MK system and confirm their status as mCP stars. We have identified 86 new promising ACV star candidates. 15 of these stars have LAMOST spectra available, which, in all cases, confirm them as classical mCP stars, which highlights the viability of our approach. The sample stars can be sorted into four subgroups characterized by distinct light curve shapes. Anti-phase variations in different photometric passbands, in particular, is a unique characteristic of a subset of ACV stars readily usable for their identification. The availability of data in three different passbands (g, r, and i) is a major advantage of the ZTF survey. On the basis of our experience with other photometric surveys and the analysis of light curves, we conclude that the ZTF is well suited for the search for, and the analysis of, ACV variables, which, however, are not considered in the available ZTF variable star catalogues. △ Less

Submitted 25 August, 2021; originally announced August 2021.

Comments: 17 pages, 8 figures, 2 tables, accepted for Astronomy & Astrophysics

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

Data-constrained Magnetohydrodynamic Simulation of a Long Duration Eruptive Flare

Authors: Yang Guo, Ze Zhong, M. D. Ding, P. F. Chen, Chun Xia, Rony Keppens

Abstract: We perform a zero-$β$ magnetohydrodynamic simulation for the C7.7 class flare initiated at 01:18 UT on 2011 June 21 using the Message Passing Interface Adaptive Mesh Refinement Versatile Advection Code (MPI-AMRVAC). The initial condition for the simulation involves a flux rope which we realize through the regularized Biot-Savart laws, whose parameters are constrained by observations from the Atmos… ▽ More We perform a zero-$β$ magnetohydrodynamic simulation for the C7.7 class flare initiated at 01:18 UT on 2011 June 21 using the Message Passing Interface Adaptive Mesh Refinement Versatile Advection Code (MPI-AMRVAC). The initial condition for the simulation involves a flux rope which we realize through the regularized Biot-Savart laws, whose parameters are constrained by observations from the Atmospheric Imaging Assembly (AIA) on the Solar Dynamics Observatory (SDO) and the Extreme Ultraviolet Imager (EUVI) on the twin Solar Terrestrial Relations Observatory (STEREO). This data-constrained initial state is then relaxed to a force-free state by the magneto-frictional module in MPI-AMRVAC. The further time-evolving simulation results reproduce the eruption characteristics obtained by SDO/AIA 94 A, 304 A, and STEREO/EUVI 304 A observations fairly well. The simulated flux rope possesses similar eruption direction, height range, and velocities to the observations. Especially, the two phases of slow evolution and fast eruption are reproduced by varying the density distribution in light of the filament material draining process. Our data-constrained simulations also show other advantages, such as a large field of view (about 0.76 solar radii). We study the twist of the magnetic flux rope and the decay index of the overlying field, and find that in this event, both the magnetic strapping force and the magnetic tension force are sufficiently weaker than the magnetic hoop force, thus allowing the successful eruption of the flux rope. We also find that the anomalous resistivity is necessary in keeping the correct morphology of the erupting flux rope. △ Less

Submitted 29 June, 2021; originally announced June 2021.

Comments: 25 pages, 7 figures, 1 table, accepted by ApJ

Study on the energy limits of kHz QPOs in Sco X-1 with $RXTE$ and $Insight$-HXMT observations

Authors: S. M. Jia, J. L. Qu, F. J. Lu, S. N. Zhang, S. Zhang, Y. Huang, D. H. Wang, D. K. Zhou, G. C. Xiao, Q. C. Bu, L. Chen, X. Ma, L. M. Song, L. Tao, X. L. Cao, Y. Chen, C. Z. Liu, Y. P. Xu

Abstract: We present a detailed spectral-timing analysis of the Kilohertz quasi-periodic oscillations (kHz QPOs) in Sco X-1 using the data of Rossi X-ray Timing Explorer ($RXTE$) and the Hard X-ray Modulation Telescope ($Insight$-HXMT). The energy band with detectable kHz QPOs is studied for the first time: on the horizontal branch, it is $\sim$ 6.89--24.01 keV and $\sim$ 8.68--21.78 keV for the upper and l… ▽ More We present a detailed spectral-timing analysis of the Kilohertz quasi-periodic oscillations (kHz QPOs) in Sco X-1 using the data of Rossi X-ray Timing Explorer ($RXTE$) and the Hard X-ray Modulation Telescope ($Insight$-HXMT). The energy band with detectable kHz QPOs is studied for the first time: on the horizontal branch, it is $\sim$ 6.89--24.01 keV and $\sim$ 8.68--21.78 keV for the upper and lower kHz QPOs detected by $RXTE$, and $\sim$ 9--27.5 keV for the upper kHz QPOs by $Insight$-HXMT; on the lower normal branch, the energy band is narrower. The fractional root mean square (rms) of the kHz QPOs increases with energy at lower energy, reaches a plateau at about 16 keV and 20 keV for the lower and upper peaks, and then levels off though with a large uncertainty. The simulation of the deadtime effect of $RXTE$/PCA shows that the deadtime does not affect much the search of the kHz QPOs but makes the rms amplitude underestimated. No significant QPO is detected below $\sim$ 6 keV as shown by the $RXTE$ data, implying that the kHz QPOs do not originate from the black body emission of the accretion disk and neutron star surface. In addition, with the combined analysis of the energy spectra and the absolute rms spectra of kHz QPOs, we suggest that the kHz QPOs in Sco X-1 originate from the Comptonization of the inner part of the transition layer, where the rotation sets the frequency and the inward bulk motion makes the spectrum harder. △ Less

Submitted 17 April, 2021; originally announced April 2021.

Measuring three-dimensional shapes of stable solar prominences using stereoscopic observations from SDO and STEREO

Authors: Chengrui Zhou, Chun Xia, Yuandeng Shen

Abstract: Although the real shapes and trajectories of erupting solar prominences in three dimensions have been intensively studied, the three-dimensional (3D) shapes of stable prominences before eruptions have not been measured accurately so far. We intend to make such a measurement to constrain 3D prominence models and to extend our knowledge of prominences. Using multiperspective observations from the At… ▽ More Although the real shapes and trajectories of erupting solar prominences in three dimensions have been intensively studied, the three-dimensional (3D) shapes of stable prominences before eruptions have not been measured accurately so far. We intend to make such a measurement to constrain 3D prominence models and to extend our knowledge of prominences. Using multiperspective observations from the Atmospheric Imaging Assembly on board SDO and the Extreme Ultraviolet Imager on board STEREO, we reconstructed 3D coordinates of three stable prominences: a quiescent, an intermediate, and a mixed type. Based on the 3D coordinates, we measured the height, length, and inclination angle of the legs of these prominences. To study the spatial relationship between the footpoints of prominences and photospheric magnetic structures, we also used the Global Oscillation Network Group H alpha images and magnetograms from the HMI on board the SDO. In three stable prominences, we find that the axes of the prominence legs are inclined by 68 degrees on average to the solar surface. Legs at different locations along a prominence axis have different heights with a two- to threefold difference. Our investigation suggests that over 96% of prominence footpoints in a sample of 70 footpoints are located at supergranular boundaries. The widths of two legs have similar values measured in two orthogonal lines of sight. We also find that a prominence leg above the solar limb showed horizontal oscillations with larger amplitudes at higher locations. With a limited image resolution and number of cases, our measurement suggests that the legs of prominences may have various orientations and do not always stand vertically on the surface of the sun. Moreover, the locations of prominence legs are closely related to supergranules. △ Less

Submitted 12 March, 2021; originally announced March 2021.

Comments: 11 figures

Journal ref: A&A 647, A112 (2021)

Optical properties of elongated conducting grains

Authors: X. M. Huang, Qi Li, Aigen Li, J. H. Chen, F. Z. Liu, C. Y. Xiao

Abstract: Extremely elongated, conducting dust particles (also known as metallic "needles" or "whiskers") are seen in carbonaceous chondrites and in samples brought back from the Itokawa asteroid. Their formation in protostellar nebulae and subsequent injection into the interstellar medium have been demonstrated, both experimentally and theoretically. Metallic needles have been suggested to explain a wide v… ▽ More Extremely elongated, conducting dust particles (also known as metallic "needles" or "whiskers") are seen in carbonaceous chondrites and in samples brought back from the Itokawa asteroid. Their formation in protostellar nebulae and subsequent injection into the interstellar medium have been demonstrated, both experimentally and theoretically. Metallic needles have been suggested to explain a wide variety of astrophysical phenomena, ranging from the mid-infrared interstellar extinction at ~3--8 micron to the thermalization of starlight to generate the cosmic microwave background. To validate (or invalidate) these suggestions, an accurate knowledge of the optics (e.g., the amplitude and the wavelength dependence of the absorption cross sections) of metallic needles is crucial. Here we calculate the absorption cross sections of iron needles of various aspect ratios over a wide wavelength range, by exploiting the discrete dipole approximation, the most powerful technique for rigorously calculating the optics of irregular or nonspherical grains. Our calculations support the earlier findings that the antenna theory and the Rayleigh approximation, which are often taken to approximate the optical properties of metallic needles are indeed inapplicable. △ Less

Submitted 11 March, 2021; originally announced March 2021.

Comments: 7 pages, 5 figures; accepted for publication in MNRAS

Accretion Torque Reversals in GRO J1008-57 Revealed by Insight-HXMT

Authors: W. Wang, Y. M. Tang, Y. L. Tuo, P. R. Epili, S. N. Zhang, L. M. Song, F. J. Lu, J. L. Qu, S. Zhang, M. Y. Ge, Y. Huang, B. Li, Q. C. Bu, C. Cai, X. L. Cao, Z. Chang, L. Chen, T. X. Chen, Y. B. Chen, Y. Chen, Y. P. Chen, W. W. Cui, Y. Y. Du, G. H. Gao, H. Gao , et al. (70 additional authors not shown)

Abstract: GRO J1008-57, as a Be/X-ray transient pulsar, is considered to have the highest magnetic field in known neutron star X-ray binary systems. Observational data of the X-ray outbursts in GRO J1008-57 from 2017 to 2020 were collected by the Insight-HXMT satellite. In this work, the spin period of the neutron star in GRO J1008-57 was determined to be about 93.28 seconds in August 2017, 93.22 seconds in… ▽ More GRO J1008-57, as a Be/X-ray transient pulsar, is considered to have the highest magnetic field in known neutron star X-ray binary systems. Observational data of the X-ray outbursts in GRO J1008-57 from 2017 to 2020 were collected by the Insight-HXMT satellite. In this work, the spin period of the neutron star in GRO J1008-57 was determined to be about 93.28 seconds in August 2017, 93.22 seconds in February 2018, 93.25 seconds in June 2019 and 93.14 seconds in June 2020. GRO J1008-57 evolved in the spin-up process with a mean rate of $-(2.10\pm 0.05)\times$10$^{-4}$ s/d from 2009 -- 2018, and turned into a spin down process with a rate of $(6.7\pm 0.6)\times$10$^{-5}$ s/d from Feb 2018 to June 2019. During the type II outburst of 2020, GRO J1008-57 had the spin-up torque again. During the torque reversals, the pulse profiles and continuum X-ray spectra did not change significantly, and the cyclotron resonant scattering feature around 80 keV was only detected during the outbursts in 2017 and 2020. Based on the observed mean spin-up rate, we estimated the inner accretion disk radius in GRO J1008-57 (about 1 - 2 times of the Alfvén radius) by comparing different accretion torque models of magnetic neutron stars. During the spin-down process, the magnetic torque should dominate over the matter accreting inflow torque, and we constrained the surface dipole magnetic field $B\geq 6\times 10^{12}$ G for the neutron star in GRO J1008-57, which is consistent with the magnetic field strength obtained by cyclotron line centroid energy. △ Less

Submitted 24 February, 2021; originally announced February 2021.

Comments: 10 pages, 5 figures, and 3 tables, the Journal of High Energy Astrophysics in press

QPOs and Orbital elements of X-ray binary 4U 0115+63 during the 2017 outburst observed by Insight-HXMT

Authors: Y. Z. Ding, W. Wang, P. Zhang, Q. C. Bu, C. Cai, X. L. Cao, C. Zhi, L. Chen, T. X. Chen, Y. B. Chen, Y. Chen, Y. P. Chen, W. W. Cui, Y. Y. Du, G. H. Gao, H. Gao, M. Y. Ge, Y. D. Gu, J. Guan, C. C. Guo, D. W. Han, Y. Huang, J. Huo, S. M. Jia, W. C. Jiang , et al. (69 additional authors not shown)

Abstract: In this paper, we presented a detailed timing analysis of a prominent outburst of 4U 0115+63 detected by \textit{Insight}-HXMT in 2017 August. The spin period of the neutron star was determined to be $3.61398\pm 0.00002$ s at MJD 57978. We measured the period variability and extract the orbital elements of the binary system. The angle of periastron evolved with a rate of $0.048\pm0.003$ $yr^{-1}$.… ▽ More In this paper, we presented a detailed timing analysis of a prominent outburst of 4U 0115+63 detected by \textit{Insight}-HXMT in 2017 August. The spin period of the neutron star was determined to be $3.61398\pm 0.00002$ s at MJD 57978. We measured the period variability and extract the orbital elements of the binary system. The angle of periastron evolved with a rate of $0.048\pm0.003$ $yr^{-1}$. The light curves are folded to sketch the pulse profiles in different energy ranges. A multi-peak structure in 1-10 keV is clearly illustrated. We introduced wavelet analysis into our data analysis procedures to study QPO signals and perform a detailed wavelet analysis in many different energy ranges. Through the wavelet spectra, we report the discovery of a QPO at the frequency $\sim 10$ mHz. In addition, the X-ray light curves showed multiple QPOs in the period of $\sim 16-32 $ s and $\sim 67- 200 $ s. We found that the $\sim100$ s QPO was significant in most of the observations and energies. There exist positive relations between X-ray luminosity and their Q-factors and S-factors, while the QPO periods have no correlation with X-ray luminosity. In wavelet phase maps, we found that the pulse phase of $\sim 67- 200 $ s QPO drifting frequently while the $\sim 16-32 $ s QPO scarcely drifting. The dissipation of oscillations from high energy to low energy was also observed. These features of QPOs in 4U 0115+63 provide new challenge to our understanding of their physical origins. △ Less

Submitted 18 February, 2021; originally announced February 2021.

Comments: 14 pages, 9 figures, and 6 tables. This work has been submitted to MNRAS after the referee's report

Transition region adaptive conduction (TRAC) in multidimensional magnetohydrodynamic simulations

Authors: Yu-Hao Zhou, Wen-Zhi Ruan, Chun Xia, Rony Keppens

Abstract: In solar physics, a severe numerical challenge for modern simulations is properly representing a transition region between the million-degree hot corona and a much cooler plasma of about 10000 K (e.g., the upper chromosphere or a prominence). In previous 1D hydrodynamic simulations, the transition region adaptive conduction (TRAC) method has been proven to capture aspects better that are related t… ▽ More In solar physics, a severe numerical challenge for modern simulations is properly representing a transition region between the million-degree hot corona and a much cooler plasma of about 10000 K (e.g., the upper chromosphere or a prominence). In previous 1D hydrodynamic simulations, the transition region adaptive conduction (TRAC) method has been proven to capture aspects better that are related to mass evaporation and energy exchange. We aim to extend this method to fully multidimensional magnetohydrodynamic (MHD) settings, as required for any realistic application in the solar atmosphere. Because modern MHD simulation tools efficiently exploit parallel supercomputers and can handle automated grid refinement, we design strategies for any-dimensional block grid-adaptive MHD simulations. We propose two different strategies and demonstrate their working with our open-source MPI-AMRVAC code. We benchmark both strategies on 2D prominence formation based on the evaporation--condensation scenario, where chromospheric plasma is evaporated through the transition region and then is collected and ultimately condenses in the corona. A field-line-based TRACL method and a block-based TRACB method are introduced and compared in block grid-adaptive 2D MHD simulations. Both methods yield similar results and are shown to satisfactorily correct the underestimated chromospheric evaporation, which comes from a poor spatial resolution in the transition region. Because fully resolving the transition region in multidimensional MHD settings is virtually impossible, TRACB or TRACL methods will be needed in any 2D or 3D simulations involving transition region physics. △ Less

Submitted 15 February, 2021; originally announced February 2021.

Comments: Accepted by A&A

Journal ref: A&A 648, A29 (2021)

Insight-HXMT observations of Swift J0243.6+6124: the evolution of RMS pulse fractions at super-Eddington luminosity

Authors: P. J. Wang, L. D. Kong, S. Zhang, Y. P. Chen, S. N. Zhang, J. L. Qu, L. Ji, L. Tao, M. Y. Ge, F. J. Lu, L. Chen, L. M. Song, T. P. Li, Y. P. Xu, X. L. Cao, Y. Chen, C. Z. Liu, Q. C. Bu, C. Cai, Z. Chang, G. Chen, T. X. Chen, Y. B. Chen, W. Cui, W. W. Cui , et al. (95 additional authors not shown)

Abstract: Based on Insight-HXMT data, we report on the pulse fraction evolution during the 2017-2018 outburst of the newly discovered first Galactic ultraluminous X-ray source (ULX) Swift J0243.6+6124. The pulse fractions of 19 observation pairs selected in the rising and fading phases with similar luminosity are investigated. The results show a general trend of the pulse fraction increasing with luminosity… ▽ More Based on Insight-HXMT data, we report on the pulse fraction evolution during the 2017-2018 outburst of the newly discovered first Galactic ultraluminous X-ray source (ULX) Swift J0243.6+6124. The pulse fractions of 19 observation pairs selected in the rising and fading phases with similar luminosity are investigated. The results show a general trend of the pulse fraction increasing with luminosity and energy at super-critical luminosity. However, the relative strength of the pulsation between each pair evolves strongly with luminosity. The pulse fraction in the rising phase is larger at luminosity below $7.71\times10^{38}$~erg~s$^{-1}$, but smaller at above. A transition luminosity is found to be energy independent. Such a phenomena is firstly confirmed by Insight-HXMT observations and we speculate it may have relation with the radiation pressure dominated accretion disk. △ Less

Submitted 24 December, 2020; originally announced December 2020.

Physical origin of the nonphysical spin evolution of MAXI J1820+070

Authors: J. Guan, L. Tao, J. L. Qu, S. N. Zhang, W. Zhang, S. Zhang, R. C. Ma, M. Y. Ge, L. M. Song, F. J. Lu, T. P. Li, Y. P. Xu, Y. Chen, X. L. Cao, C. Z. Liu, L. Zhang, Y. N. Wang, Y. P. Chen, Q. C. Bu, C. Cai, Z. Chang, L. Chen, T. X. Chen, Y. B. Chen, W. W. Cui , et al. (70 additional authors not shown)

Abstract: We report on the Insight-HXMT observations of the new black hole X-ray binary MAXI J1820+070 during its 2018 outburst. Detailed spectral analysis via the continuum fitting method shows an evolution of the inferred spin during its high soft sate. Moreover, the hardness ratio, the non-thermal luminosity and the reflection fraction also undergo an evolution, exactly coincident to the period when the… ▽ More We report on the Insight-HXMT observations of the new black hole X-ray binary MAXI J1820+070 during its 2018 outburst. Detailed spectral analysis via the continuum fitting method shows an evolution of the inferred spin during its high soft sate. Moreover, the hardness ratio, the non-thermal luminosity and the reflection fraction also undergo an evolution, exactly coincident to the period when the inferred spin transition takes place. The unphysical evolution of the spin is attributed to the evolution of the inner disc, which is caused by the collapse of a hot corona due to condensation mechanism or may be related to the deceleration of a jet-like corona. The studies of the inner disc radius and the relation between the disc luminosity and the inner disc radius suggest that, only at a particular epoch, did the inner edge of the disc reach the innermost stable circular orbit and the spin measurement is reliable. We then constrain the spin of MAXI J1820+070 to be a*=0.2^{+0.2}_{-0.3}. Such a slowly spinning black hole possessing a strong jet suggests that its jet activity is driven mainly by the accretion disc rather than by the black hole spin. △ Less

Submitted 31 March, 2021; v1 submitted 22 December, 2020; originally announced December 2020.

Comments: 14 pages, 13 figures, 5 tables, accepted for publication in MNRAS

MSC Class: 83C57

Flare stars in nearby Galactic open clusters based on {\it TESS} data

Authors: O. Maryeva, K. Bicz, C. Xia, M. Baratella, P. Čechvala, K. Vida

Abstract: The study is devoted to search for flare stars among confirmed members of Galactic open clusters using high-cadence photometry from {\it TESS} mission. We analyzed 957 high-cadence light curves of members from 136 open clusters. As a result, 56 flare stars were found, among them 8 hot B-A type objects. Of all flares, 63\% were detected in a sample of cool stars ($T_{\rm eff}<5000$~K), and 29\% --… ▽ More The study is devoted to search for flare stars among confirmed members of Galactic open clusters using high-cadence photometry from {\it TESS} mission. We analyzed 957 high-cadence light curves of members from 136 open clusters. As a result, 56 flare stars were found, among them 8 hot B-A type objects. Of all flares, 63\% were detected in a sample of cool stars ($T_{\rm eff}<5000$~K), and 29\% -- in stars of spectral type G, while 23\% in K-type stars and approximately 34\% of all detected flares are in M-type stars. Using the FLATW'RM (FLAre deTection With Ransac Method) flare finding algorithm, we estimated parameters of flares and rotation period of detected flare stars. The flare with the largest amplitude appears on the M3 type EQ\,Cha star. Statistical analysis did not reveal any direct correlation between ages, rotation periods and flaring activity. △ Less

Submitted 17 December, 2020; originally announced December 2020.

Comments: 20 pages, 9 figures, article based on the project work at ERASMUS+ GATE 2020 Summer school organised by Masaryk university Brno, Czech Republic. Accepted for publication to Contrib. Astron. Obs. Skalnate Pleso

X-ray reprocessing in accreting pulsar GX 301-2 observed with Insight-HXMT

Authors: L. Ji, V. Doroshenko, V. Suleimanov, A. Santangelo, M. Orlandini, J. Liu, L. Ducci, S. N. Zhang, A. Nabizadeh, D. Gavran, S. Zhang, M. Y. Ge, X. B. Li, L. Tao, Q. C. Bu, J. L. Qu, F. J. Lu, L. Chen, L. M. Song, T. P. Li, Y. P. Xu, X. L. Cao, Y. Chen, C. Z. Liu, C. Cai , et al. (78 additional authors not shown)

Abstract: We investigate the absorption and emission features in observations of GX 301-2 detected with Insight-HXMT/LE in 2017-2019. At different orbital phases, we found prominent Fe Kalpha, Kbeta and Ni Kalpha lines, as well as Compton shoulders and Fe K-shell absorption edges. These features are due to the X-ray reprocessing caused by the interaction between the radiation from the source and surrounding… ▽ More We investigate the absorption and emission features in observations of GX 301-2 detected with Insight-HXMT/LE in 2017-2019. At different orbital phases, we found prominent Fe Kalpha, Kbeta and Ni Kalpha lines, as well as Compton shoulders and Fe K-shell absorption edges. These features are due to the X-ray reprocessing caused by the interaction between the radiation from the source and surrounding accretion material. According to the ratio of iron lines Kalpha and Kbeta, we infer the accretion material is in a low ionisation state. We find an orbital-dependent local absorption column density, which has a large value and strong variability around the periastron. We explain its variability as a result of inhomogeneities of the accretion environment and/or instabilities of accretion processes. In addition, the variable local column density is correlated with the equivalent width of the iron Kalpha lines throughout the orbit, which suggests that the accretion material near the neutron star is spherically distributed. △ Less

Submitted 4 December, 2020; originally announced December 2020.

Comments: 10 pages, 5 figures, 2 tables, accepted for publication in MNRAS

Insight-HXMT observations of a possible fast transition from jet to wind dominated state during a huge flare of GRS~1915+105

Authors: L. D. Kong, S. Zhang, Y. P. Chen, S. N. Zhang, L. Ji, P. J. Wang, L. Tao, M. Y. Ge, C. Z. Liu, L. M. Song, F. J. Lu, J. L. Qu, T. P. Li, Y. P. Xu, X. L. Cao, Y. Chen, Q. C. Bu, C. Cai, Z. Chang, G. Chen, L. Chen, T. X. Chen, W. W. Cui, Y. Y. Du, G. H. Gao , et al. (71 additional authors not shown)

Abstract: We present the analysis of the brightest flare that was recorded in the \emph{Insight}-HMXT data set, in a broad energy range (2$-$200 keV) from the microquasar GRS~1915+105 during an unusual low-luminosity state. This flare was detected by \emph{Insight}-HXMT among a series of flares during 2 June 2019 UTC 16:37:06 to 20:11:36, with a 2-200 keV luminosity of 3.4$-$7.27$\times10^{38}$ erg s… ▽ More We present the analysis of the brightest flare that was recorded in the \emph{Insight}-HMXT data set, in a broad energy range (2$-$200 keV) from the microquasar GRS~1915+105 during an unusual low-luminosity state. This flare was detected by \emph{Insight}-HXMT among a series of flares during 2 June 2019 UTC 16:37:06 to 20:11:36, with a 2-200 keV luminosity of 3.4$-$7.27$\times10^{38}$ erg s$^{-1}$. Basing on the broad-band spectral analysis, we find that the flare spectrum shows different behaviors during bright and faint epochs. The spectrum of the flare can be fitted with a model dominated by a power-law component. Additional components show up in the bright epoch with a hard tail and in the faint epoch with an absorption line $\sim$ 6.78 keV. The reflection component of the latter is consistent with an inner disk radius $\sim$ 5 times larger than that of the former. These results on the giant flare during the "unusual" low-luminosity state of GRS~1915+105 may suggest that the source experiences a possible fast transition from a jet-dominated state to a wind-dominated state. We speculate that the evolving accretion disk and the large-scale magnetic field may play important roles in this peculiar huge flare. △ Less

Submitted 4 December, 2020; originally announced December 2020.

Nuclear pasta structures and symmetry energy

Authors: Cheng-Jun Xia, Toshiki Maruyama, Nobutoshi Yasutake, Toshitaka Tatsumi, Ying-Xun Zhang

Abstract: In the framework of the relativistic mean field model with Thomas-Fermi approximation, we study the structures of low density nuclear matter in a three-dimensional geometry with reflection symmetry. The numerical accuracy and efficiency are improved by expanding the mean fields according to fast cosine transformation and considering only one octant of the unit cell. The effect of finite cell size… ▽ More In the framework of the relativistic mean field model with Thomas-Fermi approximation, we study the structures of low density nuclear matter in a three-dimensional geometry with reflection symmetry. The numerical accuracy and efficiency are improved by expanding the mean fields according to fast cosine transformation and considering only one octant of the unit cell. The effect of finite cell size is treated carefully by searching for the optimum cell size. Typical pasta structures (droplet, rod, slab, tube, and bubble) arranged in various crystalline configurations are obtained for both fixed proton fractions and $β$-equilibration. It is found that the properties of droplets/bubbles are similar in body-centered cubic (BCC) and face-centered cubic (FCC) lattices, where the FCC lattice generally becomes more stable than BCC lattice as density increases. For the rod/tube phases, the honeycomb lattice is always more stable than the simple one. By introducing an $ω$-$ρ$ cross coupling term, we further examine the pasta structures with a smaller slope of symmetry energy $L = 41.34$ MeV, which predicts larger onset densities for core-crust transition and non-spherical nuclei. Such a variation due to the reduction of $L$ is expected to have impacts on various properties in neutron stars, supernova dynamics, and binary neutron star mergers. △ Less

Submitted 19 May, 2021; v1 submitted 2 December, 2020; originally announced December 2020.

Journal ref: Phys. Rev. C 103, 055812 (2021)

A roadmap to strange star

Authors: Renxin Xu, Xiaoyu Lai, Chengjun Xia

Abstract: What if normal baryonic matter is compressed so tightly that atomic nuclei come into close contact? This question has been asked since 1930s. The fist answer was presented by Lev Landau whose speculation has been developed, and the concept of neutron star is then popularized. However, another answer is related to strange star, which becomes worthy of attention especially after the establishment of… ▽ More What if normal baryonic matter is compressed so tightly that atomic nuclei come into close contact? This question has been asked since 1930s. The fist answer was presented by Lev Landau whose speculation has been developed, and the concept of neutron star is then popularized. However, another answer is related to strange star, which becomes worthy of attention especially after the establishment of the standard model of particle physics in 1960s. The basic ideas of this study are introduced pedagogically. We must point out emphatically that flavour symmetry of and strong coupling between quarks would be essential in seeking true answer to the question. The final answer is expected to appear in the era of multimessenger astronomy. It is emphasized too that, besides the differences of global properties (e.g., mass-radius relation, maximum mass, tidal deformability), the strong-bound surface of strange star (rather than the gravity-bound one for conventional neutron star) could play an important role in identifying a strange star by astronomical observations. △ Less

Submitted 14 November, 2020; originally announced November 2020.

Comments: to appear in Astronomische Nachrichten

A variable ionized disk wind in the black-hole candidate EXO 1846-031

Authors: Yanan Wang, Long Ji, Javier A. Garcia, Thomas Dauser, Mariano Mendez, Junjie Mao, L. Tao, Diego Altamirano, Pierre Maggi, S. N. Zhang, M. Y. Ge, L. Zhang, J. L. Qu, S. Zhang, X. Ma, F. J. Lu, T. P. Li, Y. Huang, S. J. Zheng, Z. Chang, Y. L. Tuo, L. M. Song, Y. P. Xu, Y. Chen, C. Z. Liu , et al. (66 additional authors not shown)

Abstract: After 34 years, the black-hole candidate EXO 1846-031 went into outburst again in 2019. We investigate its spectral properties in the hard intermediate and the soft states with NuSTAR and Insight-HXMT. A reflection component has been detected in the two spectral states but possibly originating from different illumination spectra: in the intermediate state, the illuminating source is attributed to… ▽ More After 34 years, the black-hole candidate EXO 1846-031 went into outburst again in 2019. We investigate its spectral properties in the hard intermediate and the soft states with NuSTAR and Insight-HXMT. A reflection component has been detected in the two spectral states but possibly originating from different illumination spectra: in the intermediate state, the illuminating source is attributed to a hard coronal component, which has been commonly observed in other X-ray binaries, whereas in the soft state the reflection is probably produced by the disk self-irradiation. Both cases support EXO 1846-031 as a low inclination system of ~40 degrees. An absorption line is clearly detected at ~7.2 keV in the hard intermediate state, corresponding to a highly ionized disk wind (log ξ > 6.1) with a velocity up to 0.06c. Meanwhile, quasi-simultaneous radio emissions have been detected before and after the X-rays, implying the co-existence of disk winds and jets in this system. Additionally, the observed wind in this source is potentially driven by magnetic forces. The absorption line disappeared in the soft state and a narrow emission line appeared at ~6.7 keV on top of the reflection component, which may be evidence for disk winds, but data with the higher spectral resolution are required to examine this. △ Less

Submitted 28 December, 2020; v1 submitted 27 October, 2020; originally announced October 2020.

Comments: Accepted for publication in ApJ

Major Scientific Challenges and Opportunities in Understanding Magnetic Reconnection and Related Explosive Phenomena in Solar and Heliospheric Plasmas

Authors: H. Ji, J. Karpen, A. Alt, S. Antiochos, S. Baalrud, S. Bale, P. M. Bellan, M. Begelman, A. Beresnyak, A. Bhattacharjee, E. G. Blackman, D. Brennan, M. Brown, J. Buechner, J. Burch, P. Cassak, B. Chen, L. -J. Chen, Y. Chen, A. Chien, L. Comisso, D. Craig, J. Dahlin, W. Daughton, E. DeLuca , et al. (83 additional authors not shown)

Abstract: Magnetic reconnection underlies many explosive phenomena in the heliosphere and in laboratory plasmas. The new research capabilities in theory/simulations, observations, and laboratory experiments provide the opportunity to solve the grand scientific challenges summarized in this whitepaper. Success will require enhanced and sustained investments from relevant funding agencies, increased interagen… ▽ More Magnetic reconnection underlies many explosive phenomena in the heliosphere and in laboratory plasmas. The new research capabilities in theory/simulations, observations, and laboratory experiments provide the opportunity to solve the grand scientific challenges summarized in this whitepaper. Success will require enhanced and sustained investments from relevant funding agencies, increased interagency/international partnerships, and close collaborations of the solar, heliospheric, and laboratory plasma communities. These investments will deliver transformative progress in understanding magnetic reconnection and related explosive phenomena including space weather events. △ Less

Submitted 16 September, 2020; originally announced September 2020.

Comments: 4 pages (including a title page), white paper submitted to Helio2050 workshop at https://www.hou.usra.edu/meetings/helio2050/. arXiv admin note: substantial text overlap with arXiv:2004.00079

On the Absorption Properties of Metallic Needles

Authors: C. Y. Xiao, Qi Li, Aigen Li, J. H. Chen

Abstract: Needle-like metallic particles have been suggested to explain a wide variety of astrophysical phenomena, ranging from the mid-infrared interstellar extinction to the thermalization of starlight to generate the cosmic microwave background. These suggestions rely on the amplitude and the wavelength dependence of the absorption cross sections of metallic needles. On the absence of an exact solution t… ▽ More Needle-like metallic particles have been suggested to explain a wide variety of astrophysical phenomena, ranging from the mid-infrared interstellar extinction to the thermalization of starlight to generate the cosmic microwave background. These suggestions rely on the amplitude and the wavelength dependence of the absorption cross sections of metallic needles. On the absence of an exact solution to the absorption properties of metallic needles, their absorption cross sections are often derived from the antenna approximation. However, it is shown here that the antenna approximation is not an appropriate representation since it violates the Kramers-Kronig relation. Stimulated by the recent discovery of iron whiskers in asteroid Itokawa and graphite whiskers in carbonaceous chondrites, we call for rigorous calculations of the absorption cross sections of metallic needle-like particles, presumably with the discrete dipole approximation. We also call for experimental studies of the formation and growth mechanisms of metallic needle-like particles as well as experimental measurements of the absorption cross sections of metallic needles of various aspect ratios over a wide wavelength range to bound theoretical calculations. △ Less

Submitted 17 September, 2020; originally announced September 2020.

Comments: 5 pages, 1 figure; accepted for publication in MNRAS

Merging strangeon stars II: the ejecta and light curves

Authors: X. Y. Lai, C. J. Xia, Y. W. Yu, R. X. Xu

Abstract: The state of supranuclear matter in compact stars remains puzzling, and it is argued that pulsars could be strangeon stars. The consequences of merging double strangeon stars are worth exploring, especially in the new era of multi-messenger astronomy. To develop the "strangeon kilonova" scenario proposed in Paper I, we make a qualitative description about the evolution of ejecta and light curves f… ▽ More The state of supranuclear matter in compact stars remains puzzling, and it is argued that pulsars could be strangeon stars. The consequences of merging double strangeon stars are worth exploring, especially in the new era of multi-messenger astronomy. To develop the "strangeon kilonova" scenario proposed in Paper I, we make a qualitative description about the evolution of ejecta and light curves for merging double strangeon stars. In the hot environment of the merger, the strangeon nuggets ejected by tidal disruption and hydrodynamical squeezing would suffer from evaporation, in which process particles, such as strangeons, neutrons and protons, are emitted. Taking into account both the evaporation of strangeon nuggets and the decay of strangeons, most of the strangeon nuggets would turn into neutrons and protons, within dozens of milliseconds after being ejected. The evaporation rates of different particles depend on temperature, and we find that the ejecta could end up with two components, with high and low opacity respectively. The high opacity component would be in the directions around the equatorial plane, and the low opacity component would be in a broad range of angular directions. The bolometric light curves show that even if the total ejected mass would be as low as $\sim 10^{-4} M_\odot$, the spin-down power of the long-lived remnant would account for the whole emission of kilonova AT2017gfo associated with GW 170817. The detailed picture of merging double strangeon stars is expected to be tested by future numerical simulations. △ Less

Submitted 26 June, 2021; v1 submitted 13 September, 2020; originally announced September 2020.

Constraining light dark matter upscattered by ultrahigh-energy cosmic rays

Authors: Chen Xia, Yan-Hao Xu, Yu-Feng Zhou

Abstract: Light halo dark matter (DM) particles upscattered by high-energy cosmic rays (CRs) can be energetic, and become detectable by conventional direct detection experiments. The current constraints derived from space-based direct CR measurements can reach $\mathcal{O}(10^{-31})\text{ cm}^{2}$ for a constant DM-nucleon scattering cross section. We show that if the CR energy spectrum follows a power law… ▽ More Light halo dark matter (DM) particles upscattered by high-energy cosmic rays (CRs) can be energetic, and become detectable by conventional direct detection experiments. The current constraints derived from space-based direct CR measurements can reach $\mathcal{O}(10^{-31})\text{ cm}^{2}$ for a constant DM-nucleon scattering cross section. We show that if the CR energy spectrum follows a power law of type $\sim E^{-3}$, the derived constraints on the scattering cross section will be highly insensitive to DM particle mass. This suggests that ultrahigh-energy CRs (UHECRs) indirectly measured by ground-based detectors can be used to place constraints on ultralight DM particles, as $E^{-3}$ is a very good approximation of the UHECR energy spectrum up to energy $\sim10^{20}\text{ eV}$. Using the recent UHECR flux data, we show that the current constraints derived from space-based CR measurements can in principle be extended to ultralight DM particles far below eV scale. △ Less

Submitted 11 July, 2021; v1 submitted 1 September, 2020; originally announced September 2020.

Comments: 25 pages, 6 figures, match the publsied version

Journal ref: Nucl.Phys.B 969 (2021) 115470