Electromagnetic Counterparts Research Articles

Origin of the black hole spin in lower-mass-gap black hole-neutron star binaries

During the fourth observing run, the LIGO-Virgo-KAGRA Collaboration reported the detection of a coalescing compact binary (GW230529$_ $181500) with component masses estimated at $2.5-4.5\, M_ and $1.2-2.0\, M_ with 90<!PCT!> credibility. Given the current constraints on the maximum neutron star (NS) mass, this event is most likely a lower-mass-gap (LMG) black hole-neutron star (BHNS) binary. The spin magnitude of the BH, especially when aligned with the orbital angular momentum, is critical in determining whether the NS is tidally disrupted. An LMG BHNS merger with a rapidly spinning BH is an ideal candidate for producing electromagnetic counterparts. However, no such signals have been detected. In this study, we employ a detailed binary evolution model that incorporates new dynamical tide implementations to explore the origin of BH spin in an LMG BHNS binary. If the NS forms first, the BH progenitor (He-rich star) must begin in orbit shorter than 0.35 days to spin up efficiently, potentially achieving a spin magnitude of $ BH > 0.3$. Alternatively, if a nonspinning BH (e.g., $M_ BH = 3.6\, M_ forms first, it can accrete up to $ 0.2\, M_ via case BA mass transfer (MT), reaching a spin magnitude of $ BH 0.18$ under Eddington-limited accretion. With a higher Eddington accretion limit (i.e., 10.0 $ M Edd $), the BH can attain a significantly higher spin magnitude of $ BH by accreting approximately $1.0\, M_ during case BA MT phase.

Host Galaxy Demographics Of Individually Detectable Supermassive Black-hole Binaries with Pulsar Timing Arrays

Abstract Massive black hole binaries (MBHBs) produce gravitational waves (GWs) that are detectable with pulsar timing arrays. We determine the properties of the host galaxies of simulated MBHBs at the time they are producing detectable GW signals. The population of MBHB systems we evaluate is from the \textit{Illustris} cosmological simulations taken in tandem with post processing semi-analytic models of environmental factors in the evolution of binaries. Upon evolving to the GW frequency regime accessible by pulsar timing arrays, we calculate the detection probability of each system using a variety of different values for pulsar noise characteristics in a plausible near-future International Pulsar Timing Array dataset. We find that detectable systems have host galaxies that are clearly distinct from the overall binary population and from most galaxies in general. With conservative noise factors, we find that host stellar metallicity, for example, peaks at $\sim2Z_\odot$ as opposed to the total population of galaxies which peaks at $\sim0.6Z_\odot$. Additionally, the most detectable systems are much brighter in magnitude and more red in color than the overall population, indicating their likely identity as large ellipticals with diminished star formation. These results can be used to develop effective search strategies for identifying host galaxies and electromagnetic counterparts following GW detection by pulsar timing arrays.

Characterizing gravitational wave detector networks: from A♯ to cosmic explorer

Abstract Gravitational-wave observations by the Laser Interferometer Gravitational-Wave Observatory (LIGO) and Virgo have provided us a new tool to explore the Universe on all scales from nuclear physics to the cosmos and have the massive potential to further impact fundamental physics, astrophysics, and cosmology for decades to come. In this paper we have studied the science capabilities of a network of LIGO detectors when they reach their best possible sensitivity, called A#, given the infrastructure in which they exist and a new generation of observatories that are factor of 10 to 100 times more sensitive (depending on the frequency), in particular a pair of L-shaped Cosmic Explorer observatories (one 40 km and one 20 km arm length) in the US and the triangular Einstein Telescope with 10 km arms in Europe. We use science metrics derived from the top priorities of several funding agencies to characterize the science capabilities of different networks. The presence of one or two A# observatories in a network containing two or one next generation observatories, respectively, will provide good localization capabilities for facilitating multimessenger astronomy and precision measurement of the Hubble parameter. Two Cosmic Explorer observatories are indispensable for achieving precise localization of binary neutron stars, facilitating detection of electromagnetic counterparts and transforming multimessenger astronomy. Their combined operation is even more important in the detection and localization of high-redshift sources, such as binary neutron stars, beyond the star-formation peak, and primordial black hole mergers, which may occur roughly 100 million years after the Big Bang. The addition of the Einstein Telescope to a network of two Cosmic Explorer observatories is critical for accomplishing all the identified science metrics. For most metrics the triple network of next generation terrestrial observatories are a factor 100 better than what can be accomplished by a network of three A# observatories.

Cross-correlating Dark Sirens and Galaxies: Constraints on H 0 from GWTC-3 of LIGO–Virgo–KAGRA

We apply the cross-correlation technique to infer the Hubble constant (H 0) of the Universe using gravitational-wave (GW) sources without electromagnetic counterparts (dark sirens) from the third GW Transient Catalog (GWTC-3) and the photometric galaxy surveys 2MPZ and WISE-SuperCOSMOS, and combine these with the bright siren measurement of H 0 from GW170817. The posterior on H 0 with only dark sirens is uninformative due to the small number of well-localized GW sources. Using the eight well-localized dark sirens and the binary neutron star GW170817 with electromagnetic counterpart, we obtain a value of the Hubble constant H0=75.4−6+11 km s−1 Mpc−1 (median and 68.3% equal-tailed interval) after marginalizing over the matter density and the GW bias parameters. This measurement is mostly driven by the bright siren measurement, and any constraint from dark sirens is not statistically significant. In the future, with more well-localized GW events, the constraints on expansion history will improve.

The glow of axion quark nugget dark matter

Context. The existence of axion quark nuggets is a potential consequence of the axion field, which provides a possible solution to the charge-conjugation parity violation in quantum chromodynamics. In addition to explaining the cosmological discrepancy of matter-antimatter asymmetry and a visible-to-dark-matter ratio of Ωdark/Ωvisible ≃ 5, these composite compact objects are expected to represent a potentially ubiquitous electromagnetic background radiation by interacting with ordinary baryonic matter. We conducted an in-depth analysis of axion quark nugget-baryonic matter interactions in the environment of the intracluster medium in the constrained cosmological Simulation of the LOcal Web (SLOW). Aims. Here, we aim to provide upper limit predictions on electromagnetic counterparts of axion quark nuggets in the environment of galaxy clusters by inferring their thermal and non-thermal emission spectrum originating from axion quark nugget-cluster gas interactions. Methods. We analyzed the emission of axion quark nuggets in a large sample of 161 simulated galaxy clusters using the SLOW simulation. These clusters are divided into a sub-sample of 150 galaxy clusters, ordered in five mass bins ranging from 0.8 to 31.7 × 1014 M⊙, along with 11 cross-identified galaxy clusters from observations. We investigated dark matter-baryonic matter interactions in galaxy clusters in their present stage at the redshift of z = 0 by assuming all dark matter consists of axion quark nuggets. The resulting electromagnetic signatures were compared to thermal Bremsstrahlung and non-thermal cosmic ray (CR) synchrotron emission in each galaxy cluster. We further investigated individual frequency bands imitating the observable range of the WMAP, Planck, Euclid, and XRISM telescopes for the most promising cross-identified galaxy clusters hosting detectable signatures of axion quark nugget emission. Results. We observed a positive excess in the low- and high-energy frequency windows, where thermal and non-thermal axion quark nugget emission can significantly contribute to (or even outshine) the emission of the intracluster medium (ICM) in frequencies up to νT ≲ 3842.19 GHz and νT ϵ [3.97, 10.99] × 1010GHz, respectively. Emission signatures of axion quark nuggets are found to be observable if CR synchrotron emission of individual clusters is sufficiently low. The degeneracy in the parameters contributing to an emission excess makes it challenging to offer predictions with respect to pinpointing specific regions of a positive axion quark nugget excess; however, a general increase in the total galaxy cluster emission is expected based on this dark matter model. Axion quark nuggets constitute an increment of 4.80% of the total galaxy cluster emission in the low-energy regime of νT ≲ 3842.19 GHz for a selection of cross-identified galaxy clusters. We propose that the Fornax and Virgo clusters represent the most promising candidates in the search for axion quark nugget emission signatures. Conclusions. The results from our simulations point towards the possibility of detecting an axion quark nugget excess in galaxy clusters in observations if their signatures can be sufficiently disentangled from the ICM radiation. While this model proposes a promising explanation for the composition of dark matter, with the potential to have this outcome verified by observations, we propose further changes that are aimed at refining our methods. Our ultimate goal is to identify the extracted electromagnetic counterparts of axion quark nuggets with even greater precision in the near future.

Off-Axis Color Characteristics of Binary Neutron Star Merger Events: Applications for Space Multi-Band Variable Object Monitor and James Webb Space Telescope

With advancements in gravitational wave detection technology, an increasing number of binary neutron star (BNS) merger events are expected to be detected. Due to the narrow opening angle of jet cores, many BNS merger events occur off-axis, resulting in numerous gamma-ray bursts (GRBs) going undetected. Models suggest that kilonovae, which can be observed off-axis, offer more opportunities to be detected in the optical/near-infrared band as electromagnetic counterparts of BNS merger events. In this study, we calculate kilonova emission using a three-dimensional semi-analytical code and model the GRB afterglow emission with the open-source Python package afterglowpy at various inclination angles. Our results show that it is possible to identify the kilonova signal from the observed color evolution of BNS merger events. We also deduce the optimal observing window for SVOM/VT and JWST/NIRCam, which depends on the viewing angle, jet opening angle, and circumburst density. These parameters can be cross-checked with the multi-band afterglow fitting. We suggest that kilonovae are more likely to be identified at larger inclination angles, which can also help determine whether the observed signals without accompanying GRBs originate from BNS mergers.

A simulation study on the sub-threshold joint gravitational wave-electromagnetic wave observation on binary neutron star mergers

Abstract The coalescence of binary neutron stars (BNS) is a prolific source of gravitational waves (GWs) and electromagnetic (EM) radiation, offering a dual observational window into the Universe. Lowering the signal-to-noise ratio (S/N) threshold is a simple and cost-effective way to enhance the detection probability of GWs from BNS mergers. In this study, we introduce a metric of the purity of joint GW and EM detections Pjoint, which is in analogue to Pastro in GW only observations. By simulating BNS merger GWs jointly detected by the HLV network and EM counterparts (kilonovae and short Gamma-ray bursts, sGRBs) with an assumed merger rate density of BNS, we generate catalogs of GW events and EM counterparts. Through this simulation, we analyze joint detection pairs, both correct and misidentified. We find the following: 1. For kilonovae, requiring Pjoint > 95% instead of $P_{\rm astro}>95~{{\%}}$ reduces the S/N from 9.2 to 8.5-8.8, allowing 5-13 additional joint detections per year and increasing the GW detection volume by 9-17%; 2. For sGRBs, requiring Pjoint > 95% instead of Pastro reduces the S/N from 9.2 to 8.1-8.5; 3. Increasing kilonova or sGRB detection capability does not improve Pjoint due to a higher rate of misidentifications. We also show that sub-threshold GW and kilonova detections can reduce the uncertainty in measuring the Hubble constant to 89-92% of its original value, and sub-threshold GW and sGRB observations can enhance the precision of constraining the speed of GWs to 88% of previously established values.

An acceleration-radiation model for nonthermal flares from Sgr A*

is the electromagnetic counterpart of the accreting supermassive black hole in the Galactic center. Its emission is variable in the near-infrared (NIR) and X-ray wavelengths on short timescales (several minutes to a few hours). The NIR light curve displays red-noise variability, while the X-ray light curve exhibits bright flares that rise by many orders of magnitude upon the stable X-ray quiescent emission. Every X-ray flare is associated with a bright NIR flux change, but the opposite is not always true. The physical origin of NIR and X-ray flares is still under debate. We introduce a model for the production of NIR and X-ray flares flares from an active region in where particle acceleration takes place intermittently. A fraction of electrons from their thermal pool is accelerated to higher energies while they radiate via synchrotron and synchrotron self-Compton (SSC) processes. in contrast to other radiation models for flares, the particle acceleration is not assumed to be instantaneous. We studied the evolution of the particle distribution and the emitted electromagnetic radiation from the flaring region by numerically solving the kinetic equations for electrons and photons. Our calculations took the finite duration of particle acceleration, radiative energy losses, and physical escape from the flaring region into account. To gain better insight into the relation of the model parameters, we complemented our numerical study with analytical calculations. Flares are produced when the acceleration episode has a finite duration. The rising part in the light curve of a flare is related to the particle acceleration timescale, while the decay is controlled by the cooling or escape timescale of particles. The emitted synchrotron spectra are power laws whose photon index is determined by the ratio of the acceleration and escape timescales, followed by an exponential cutoff. This occurs at the characteristic synchrotron photon energy emitted by particles with the maximum Lorentz factor (where energy loss and gain rates become equal). The NIR flux increases before the onset of the X-ray flare, and the time lag is linked to the particle acceleration timescale. Bright X-ray flares, such as the one observed in 2014, have gamma -ray counterparts that might be detected by the Cherenkov Telescope Array Observatory. Our generic model for NIR and X-ray flares favors an interpretation of diffusive nonresonant particle acceleration in magnetized turbulence. If direct acceleration by the reconnection electric field in macroscopic current sheets causes the energization of particles during flares in then models considering the injection of preaccelerated particles into a blob where particles cool and/or escape would be appropriate to describe the flare.

Formation of GW230529 from Isolated Binary Evolution

In this paper, we explore the formation of the mass-gap black hole-neutron star (mgBHNS) merger detected in gravitational wave (GW) event, i.e., GW230529, from the isolated binary evolution channel, and study potential signatures of its electromagnetic counterparts. By adopting the “delayed” supernova prescription and reasonable model realizations, our population synthesis simulation results can simultaneously match the rate densities of mgBHNS and total BHNS mergers inferred from the population analyses, along with the population distribution of the BH mass in BHNS mergers reported by the LIGO–Virgo–KAGRA Collaboration. Because GW230529 contributes significantly to the inferred mgBHNS rate densities, we suggest that GW230529 can be explained through the isolated binary evolution channel. Considering the AP4 (DD2) equation of state, the probability that GW230529 can make tidal disruption is 12.8% (63.2%). If GW230529 is a disrupted event, its kilonova peak apparent magnitude is predicted ∼23–24 mag, and hence, can be detected by the present survey projects and Large Synoptic Survey Telescope. Since GW230529 could be an off-axis event inferred from the GW observation, its associated gamma-ray burst (GRB) might be too dim to be observed by γ-ray detectors, interpreting the lack of GRB observations. Our study suggests the existence of mgBHNS mergers formed through the isolated binary evolution channel due to the discovery of GW230529, indicating that BHNS mergers are still likely to be multimessenger sources that emit GWs, GRBs, and kilonovae. Although mgBHNS mergers account for ∼50% of the cosmological BHNS population, we find that ≳90% of disrupted BHNS mergers are expected to originate from mgBHNS mergers.

Signatures of low-mass black hole–neutron star mergers

ABSTRACT The recent observation of the GW230529 event indicates that black hole–neutron star binaries can contain low-mass black holes. Since lower mass systems are more favourable for tidal disruption, such events are promising candidates for multimessenger observations. In this study, we employ five finite-temperature, composition-dependent matter equations of state and present results from ten 3D general relativistic hydrodynamic simulations for the mass ratios $q = 2.6$ and 5. Two of these simulations target the chirp mass and effective spin parameter of the GW230529 event, while the remaining eight contain slightly higher mass black holes, including both spinning ($a_{\mathrm{ BH}} = 0.7$) and non-spinning ($a_{\mathrm{ BH}} = 0$) models. We discuss the impact of the equation of state, spin, and mass ratio on black hole–neutron star mergers by examining both gravitational-wave and ejected matter properties. For the low-mass ratio model we do not see fast-moving ejecta for the softest equation of state model, but the stiffer model produces on the order of $10^{-6}\,\mathrm{ M}_\odot$ of fast-moving ejecta, expected to contribute to an electromagnetic counterpart. Notably, the high-mass ratio model produces nearly the same amount of total dynamical ejecta, but yields 52 times more fast-moving ejecta than the low-mass ratio system. In addition, we observe that the black-hole spin tends to decrease the amount of fast-moving ejecta while increasing significantly the total ejected mass. Finally, we note that the disc mass tends to increase as the neutron star compactness decreases.

Electrohydraulic musculoskeletal robotic leg for agile, adaptive, yet energy-efficient locomotion

Robotic locomotion in unstructured terrain demands an agile, adaptive, and energy-efficient architecture. To traverse such terrains, legged robots use rigid electromagnetic motors and sensorized drivetrains to adapt to the environment actively. These systems struggle to compete with animals that excel through their agile and effortless motion in natural environments. We propose a bio-inspired musculoskeletal leg architecture driven by antagonistic pairs of electrohydraulic artificial muscles. Our leg is mounted on a boom arm and can adaptively hop on varying terrain in an energy-efficient yet agile manner. It can also detect obstacles through capacitive self-sensing. The leg performs powerful and agile gait motions beyond 5 Hz and high jumps up to 40 % of the leg height. Our leg’s tunable stiffness and inherent adaptability allow it to hop over grass, sand, gravel, pebbles, and large rocks using only open-loop force control. The electrohydraulic leg features a low cost of transport (0.73), and while squatting, it consumes only a fraction of the energy (1.2 %) compared to its conventional electromagnetic counterpart. Its agile, adaptive, and energy-efficient properties would open a roadmap toward a new class of musculoskeletal robots for versatile locomotion and operation in unstructured natural environments.

The Design of GECAM Scientific Ground Segment

The Gravitational wave burst high-energy Electromagnetic Counterpart All-sky Monitor (GECAM) is a dedicated mission for monitoring high-energy transients. Here we report the design of the GECAM Scientific Ground Segment (GSGS) in terms of the scientific requirements, including the architecture, the external interfaces, the main function, and workflow. Judging from the analysis and verification results during the commissioning phase, the GSGS functions well and is able to monitor the status of the payloads, adjust the parameters, develop the scientific observation plans, generate the scientific data products, analyze the data, etc. Thus, the on-orbit operation and scientific researches of GECAM are guaranteed.

In-flight Energy Calibration of the GECAM Gamma-ray Detectors

The Gravitational wave high-energy Electromagnetic Counterpart All-sky Monitor (GECAM) mission is designed to monitor the Gamma-Ray Bursts (GRBs) associated with gravitational waves and other high-energy transient sources. The mission consists of two microsatellites which are planned to operate at the opposite sides of the Earth. Each GECAM satellite could detect and localize GRBs in about 8 keV–5 MeV with its 25 Gamma-Ray Detectors (GRDs). In this work, we report the in-flight energy calibration of GRDs using the characteristic gamma-ray lines in the background spectra, and show their performance evolution during the commissioning phase. Besides, a preliminary cross-calibration of energy response with Fermi GBM data is also presented, validating the energy response of GRDs.

SiPM-based Gamma-ray Detectors of GECAM

The Gravitational Wave High-Energy Electromagnetic Counterparts All-sky Monitor (GECAM) is a space mission dedicated to detecting gamma-ray bursts associated with gravitational wave events and various cosmic phenomena. GECAM consists of several satellites, with three currently in orbit and a fourth in the orbital insertion stage. GECAM has yielded numerous scientific discoveries, including observing the brightest gamma-ray burst recorded to date. The compact Gamma-ray Detectors (GRDs) based on silicon photomultipliers (SiPMs) are a crucial component of GECAM, representing the first extensive application of SiPM technology in a spaceborne gamma-ray scientific satellite. This review article first introduces the current status of GECAM and its observational achievements, then focuses on the design and in-flight performance of the gamma-ray detectors, irradiation damage to SiPMs and mitigation strategies, and outlines future plans.

Afterglows from Binary Neutron Star Postmerger Systems Embedded in Active Galactic Nuclei Disks

The observability of afterglows from binary neutron star mergers occurring within active galactic nuclei (AGN) disks is investigated. We perform 3D GRMHD simulations of a postmerger system and follow the jet launched from the compact object. We use semianalytic techniques to study the propagation of the blast wave powered by the jet through an AGN disk-like external environment, extending to distances beyond the disk scale height. The synchrotron emission produced by the jet-driven forward shock is calculated to obtain the afterglow emission. The observability of this emission at different frequencies is assessed by comparing it to the quiescent AGN emission. In the scenarios where the afterglow could temporarily outshine the AGN, we find that detection will be more feasible at higher frequencies (≳1014 Hz) and the electromagnetic counterpart could manifest as a fast variability in the AGN emission, on timescales less than a day.

Identifying the electromagnetic counterparts of LISA massive black hole binaries in archival LSST data

ABSTRACT LSST will catalogue the light curves of up to 100 million quasars. Among these there can be $\sim$100 ultra-compact massive black hole (MBH) binaries, whose gravitational waves (GWs) can be detected 5–15 yr later by LISA. Here, we assume such a LISA detection occurred, and assess whether or not its electromagnetic (EM) counterpart can be identified as a periodic quasar in archival LSST data. We use the binary’s properties derived from the LISA waveform, including the evolution of its orbital frequency, its total mass, distance, and sky localization, to predict the redshift, magnitude, and historical periodicity of the quasar expected in the LSST data. We then use Monte Carlo simulations to compute the false alarm probability (FAP), i.e. the number of quasars in the LSST catalogue matching these properties by chance, based on the (extrapolated) quasar luminosity function, the cadence of LSST, and intrinsic ‘damped random walk’ quasar variability. We analyse four fiducial LISA binaries, with masses and redshifts of $(M_{\rm bin}/{\rm M_{\odot }},z) = (3\times 10^5,0.3)$, $(3\times 10^6,0.3)$, $(10^7,0.3)$, and $(10^7,1)$. While noise and aliasing due to LSST’s cadence produces false periodicities by chance, we find that the frequency chirp of the LISA source during the LSST observations washes out these noise peaks and allows the genuine source to stand out in appropriately scaled Lomb–Scargle periodograms. We find that all four fiducial binaries can be uniquely identified, with ${\rm FAP}\lt 10^{-5}$, a week or more before merger. This should enable follow-up EM observations targeting individual EM counterparts during their inspiral stage.

Enhancing GWOPS Capabilities for Coordinated Multi-Telescope Detection of Gravitational Wave Electromagnetic Counterparts

The groundbreaking detection of gravitational waves (GWs) has ushered in a new era of astronomical observation, granting us access to cosmic phenomena that are imperceptible to electromagnetic waves. The inherently weak GW signals coupled with the substantial uncertainties in source localization pose significant challenges to the field of astronomy. In this paper, we introduce innovative strategies to enhance the efficiency of observing electromagnetic counterparts to GW events, thereby unlocking further secrets of the cosmos. We present a novel technique for designing observation targets and establishing priorities, progressing from the epicenter to the periphery within the boundaries of the GW error sky region. This method has significantly reduced the average slewing distance of telescopes by 41% compared to traditional methods, thus enhancing observational efficiency. Additionally, we have developed a collaborative observation strategy for telescope networks, allocating observation targets based on the field-of-view (FOV) sizes of individual telescopes. This ensures comprehensive coverage without redundancy, allowing a network of four telescopes to cover a sky area and accumulate observation probability more than four times that of a single telescope operating independently over an equivalent period. Building upon these strategies, we have significantly upgraded GWOPS, the GW Follow-up Observation Planning System developed by the China-VO team, to provide precise observational planning for large FOV (greater than 1 square degree) telescope networks. The system also features a web-based user interface that presents the GW error sky area and observation planning results in a graphical format, significantly improving user interaction and experience. The research presented herein equips astronomers with a robust toolkit, advancing the efficiency of searching for and studying electromagnetic counterparts to GW events, and heralding new frontiers in the research of astrophysics and cosmology.

Direct In‐Situ Estimates of Energy and Force Balance Associated With Magnetopause Reconnection

AbstractFundamental processes in plasmas act to convert energies into different forms, for example, electromagnetic, kinetic and thermal. Direct derivation from the Vlasov‐Maxwell equation yields sets of equations that describe the temporal evolution of magnetic, kinetic and internal energies in either the monofluid or multifluid frameworks. In this work, we focus on the main terms affecting the changes in kinetic energy. These are pressure‐gradient‐related terms and electromagnetic terms. The former account for plasma acceleration/deceleration from a pressure gradient, while the latter from an electric field. Although limited spatial and temporal deviations are expected, a statistical balance between these terms is fundamental to ensure the overall conservation of energy and momentum. We use in‐situ observations from the Magnetospheric MultiScale (MMS) mission to study the relationship between these terms. We perform a statistical analysis of those parameters in the context of magnetic reconnection by focusing on small‐scale Electron Diffusion Regions and large‐scale Flux Transfer Events. The analysis reveals a correlation between the two terms in the monofluid force balance, and in the ion force and energy balance. However, the expected relationship cannot be verified from electron measurements. Generally, the pressure‐gradient‐related terms are smaller than their electromagnetic counterparts. We perform an error analysis to quantify the expected underestimation of gradient values as a function of the spacecraft separation compared to the gradient scale. Our findings highlight that MMS is capable of capturing energy and force balance for the ion fluid, but that care should be taken for energy conversion terms based on electron pressure gradients.

Kilonova Spectral Inverse Modelling with Simulation-based Inference: An Amortized Neural Posterior Estimation Analysis

Kilonovae represent a category of astrophysical transients, identifiable as the electromagnetic (EM) counterparts associated with the coalescence events of binary systems comprising neutron stars and neutron star–black hole pairs. They act as probes for heavy-element nucleosynthesis in astrophysical environments. These studies rely on an inference of the physical parameters (e.g., ejecta mass, velocity, composition) that describe kilonovae-based on EM observations. This is a complex inverse problem typically addressed with sampling-based methods such as Markov Chain Monte Carlo or nested sampling algorithms. However, repeated inferences can be computationally expensive, due to the sequential nature of these methods. This poses a significant challenge to ensuring the reliability and statistical validity of the posterior approximations and, thus, the inferred kilonova parameters themselves. We present a novel approach: simulation-based inference using simulations produced by KilonovaNet. Our method employs an ensemble of amortized neural posterior estimation (ANPE) with an embedding network to directly predict posterior distributions from simulated spectral energy distributions. We take advantage of the quasi-instantaneous inference time of ANPE to demonstrate the reliability of our posterior approximations using diagnostics tools, including coverage diagnostic and posterior predictive checks. We further test our model with real observations from AT 2017gfo, the only kilonova with multimessenger data, demonstrating agreement with previous likelihood-based methods while reducing inference time down to a few seconds. The inference results produced by ANPE appear to be conservative and reliable, paving the way for testable and more efficient kilonova parameter inference.

Insights from the Gaussian Process Method for the Fast Radio Burst–associated X-Ray Burst of SGR 1935+2154

The Gaussian process method is employed to analyze the light curves of bursts detected by Insight-HXMT, Neutron Star Interior Composition Explorer (NICER), and Gravitational Wave High-energy Electromagnetic Counterpart All-sky Monitor from SGR 1935+2154 between 2020 and 2022. It is found that a stochastically driven damped simple harmonic oscillator (SHO) is necessary to capture the characteristics of the X-ray bursts (XRBs). A variability timescale of the XRBs, corresponding to the broken frequencies in the SHO power spectral densities (PSDs), is extracted. In particular, a high broken frequency of 35 Hz where the index of the SHO PSD changes from −4 to −2 is constrained by the HXMT-HE burst associated with fast radio burst (FRB) 200428. It is suggested that the corresponding timescale of 0.03 s could be the retarding timescale of the system driven by some energy release, and the production of the HE photon should be quasi-simultaneous with the response. The other special event is a NICER burst with a retarding timescale of 1/(39 Hz) ≈ 0.02 s. In the normal XRBs, no retarding timescale is constrained; a long relax/equilibrium timescale (corresponding to a broken frequency of 1–10 Hz, where the index of the SHO PSD changes from −4/−2 to 0 in the SHO PSD) is obtained. The results indicate that the FRB-associated HXMT-HE XRB could be produced immediately when the system is responding to the energy disturbance, far before the equilibrium state.