Peer Revenue Fully funded

Papers

Fresh from arXiv, typeset properly. Highlighted words trade on the market.

Physics Astrophysics

astro-ph Watch Astrophysics (sign in to use)

Tuesday 8 September

Cosmology and Nongalactic Astrophysics · Cross-listed from hep-ph

A Vector-Like Lepton Interpretation of the High-Energy Nuclear Recoil Candidate in LUX-ZEPLIN

The LUX-ZEPLIN (LZ) experiment has recently reported a nuclear-recoil candidate at $E_R=248\pm23_{\rm stat}\pm23_{\rm sys}\,{\rm keV}$ in an exposure of $2.84\,{\rm tonne\,yr}$. At this unusually high recoil energy, the spin-independent (SI) xenon response is strongly suppressed, motivating dark matter interactions with a harder recoil spectrum. We show that elastic spin-dependent (SD) scattering through the nonrelativistic operator $\mathcal O_4=\mathbf S_{\rm DM}\cdot\mathbf S_N$ provides a simple realization of such a spectrum. We consider singlet-doublet Majorana dark matter, for which a diagonal $Z$ coupling generates the required SD interaction. The accompanying Higgs-mediated SI interaction would generically produce too many low-energy recoils, but can be suppressed along a Higgs blind spot while retaining a nonzero $Z$ coupling. We find a region of the Higgs blind spot that simultaneously predicts an $\mathcal O(1)$ event rate in the LZ high-recoil window and reproduces the observed thermal relic abundance. The same SD interaction leads to solar capture, allowing this region to be tested independently with solar-neutrino searches. Current IceCube limits already probe part of the LZ-compatible thermal parameter space, while a viable region remains. Additional xenon exposure and improved solar-neutrino searches can therefore provide complementary tests of this interpretation.
Read the full abstract

Instrumentation and Methods for Astrophysics · Cross-listed from gr-qc

TARTARUS: A High-Performance Python Code for Ray Tracing in Curved Spacetimes

The recent observations of the supermassive black holes M87* and Sagittarius A* by the Event Horizon Telescope collaboration (EHT) have shown the need for robust, efficient, and accessible numerical tools to model light propagation in strongly curved spacetimes. Comparing theoretical accretion models with observational data relies heavily on ray-tracing algorithms that calculate the trajectories of photons traveling from a source plasma to a distant observer. In this paper, we present TARTARUS (Tracer for Astrophysical Ray Trajectories Around Relativistic Ultra-compact Sources), a new, highly modular, and extensible Python-based framework for computing null geodesics around compact objects. By leveraging Just-In-Time (JIT) compilation via Numba, TARTARUS bridges the gap between Python's high-level accessibility and the execution speed of compiled languages like C++ or FORTRAN. The code supports multiple native solvers (including adaptive embedded Runge-Kutta pairs, a Bulirsch-Stoer extrapolator, and a Verlet scheme), event-detection mechanics for handling structural intersections (e.g. event horizons and accretion disks), and support for both analytical and numerically generated background metrics. We demonstrate the code's physical accuracy and computational efficiency through some tests, including the evaluation of the Hamiltonian constraint conservation and ray-tracing of the shadow and the Novikov-Thorne thin accretion disks around a Kerr black hole.
Read the full abstract

Astrophysics of Galaxies

Triaxiality of the Milky Way dark matter halo and the microlensing optical depth towards the Large Magellanic Cloud

The microlensing optical depth towards the Large Magellanic Cloud (LMC), $\tau_{\rm LMC}$, bounds the compact-object fraction of the Galactic dark matter halo, but is almost always evaluated for a spherical halo. We compute it for triaxial haloes and show that the shape dependence is carried entirely by the projection of the line of sight onto the halo principal axes: for a scale-free halo $\rho\propto R^{-\alpha}$ at fixed local density it factorises exactly as $\tau/\tau_{\rm sph}=S_\odot^{\alpha}\langle S^{-\alpha}\rangle_w$, with $S(\Omega)$ the ellipsoidal shape function. The LMC lies only $0\fdg78$ from perpendicular to the Sun-Galactic Centre line, making it a near-optimal probe of in-plane elongation and a poor probe of vertical flattening. For the halo shapes favoured by Gaia we find $\tau_{\rm LMC}/\tau_{\rm sph}=0.95$–$1.15$ at fixed local density, $1.57$ for the Law-Majewski logarithmic potential, and $0.67$–$1.54$ when the major-axis azimuth is left free. For the Einasto halo required by the declining Gaia DR3 rotation curve, truncation near $25\kpc$ is offset by the larger local density it implies and $\tau$ falls by only $7$–$15$ per cent; but the median lens distance moves inwards from $13.4$ to $8.5\kpc$, where the logarithmic slope is steep and rising, amplifying the shape sensitivity: the Gaia RR~Lyrae enhancement grows from $15$ to $23$ per cent and the orientation-free range widens to $0.62$–$1.62$. The systematic on the inferred MACHO or primordial-black-hole fraction is two-sided, of order $\pm15$ per cent for the measured shapes and a factor $2.6$ under full orientation freedom, exceeding the statistical errors of the EROS-2 and OGLE bounds.
Read the full abstract

Earth and Planetary Astrophysics

TEMPUS: Relativistic coordinate time scales for any solar-system body from arbitrary ephemerides

Precise time-scale transformations between solar-system bodies are an essential requirement for deep-space science and solar system navigation. The post-Newtonian relativistic expression relating a body-centered coordinate time TCX to the Barycentric Coordinate Time TCB depends on the gravitational environment and velocity of the target body, and must be evaluated consistently with the chosen planetary ephemerides. So far only INPOP has provided such a solution as part of its ephemerides releases, albeit only for Earth with TCG-TCB. We present TEMPUS, a tool that numerically integrates the IAU 2000 time-transformation for any solar-system body using positions and velocities from an arbitrary ephemerides. We provide a validated, ephemerides-agnostic implementation of TCX-TCB and assess its accuracy against published independent solutions. The rate $d\,\mathrm{TCX}/d\,\mathrm{TCB}$ is integrated with an Adams-Moulton order-12 predictor-corrector. The gravitational potential sum tested here includes the Sun, eight planets, Pluto, up to 343 main-belt asteroids, and 30 trans-Neptunian objects. Outputs can be stored in readable format tables of sampled values or as Chebyshev coefficient files. Comparisons are performed against the TU Dresden Klioner time solutions, the TCG-TCB solution provided with INPOP19a and INPOP21a, and the LTE440 lunar solution. For the planetary bodies TEMPUS agrees with Klioner at the $10^{-21} s\,s^{-1}$ level over 200 yr. Against the INPOP19a released TCG-TCB we also agree to the same level. The main-belt asteroids as well as the trans-Neptunian objects shift the planetary time scales at the $10^{-18} s\,s^{-1}$ level, and the choice of planetary and lunar ephemerides contributes differences on the $10^{-18} s\,s^{-1}$ level as well. TEMPUS provides a validated, ephemeris-agnostic post-Newtonian time-transformation for any solar-system body.
Read the full abstract

Astrophysics of Galaxies

Euclid: Quick Data Release (Q1) – Exploring the detailed visual morphology of galaxies in clusters

Galaxy clusters provide unique laboratories for studying environmental effects on galaxy evolution. The morphology–density ($T$–$\Sigma$) and morphology–cluster-centric radius ($T$–$R$) relations trace how galaxy morphology is influenced by the environment, but previous studies at intermediate redshifts have been limited in both sample size and radial coverage. We use the Euclid Quick Data Release 1 (Q1) visual morphology catalogue to measure the $T$–$\Sigma$ and $T$–$R$ relations for smooth, featured-or-disc, and barred galaxies in known clusters at $0.2\leq z\leq 0.5$, extending this analysis to large cluster-centric distances ($3 R_{500c}$) and studying their dependence on stellar mass. Using photometric redshifts and stellar masses provided by Euclid, we identify 1754 cluster members within $1.5 R_{500c}$, distributed across 71 clusters. We classify the identified galaxies as smooth, featured-or-disc, or barred using the predicted vote fractions provided by the Zoobot deep learning foundation model in the Q1 visual morphology catalogue. We confirm the $T$–$\Sigma$ relation in all the stellar mass ranges studied, with a stronger influence of the cluster environment on galaxy morphology in the densest parts of the cluster and closer to the cluster centre. Beyond $1.5 R_{500c}$, the morphological segregation weakens, with featured-or-disc galaxies overtaking smooth galaxies at the lowest densities, consistent with the growing contribution of field interlopers in the cluster outskirts. For barred galaxies, we find a tentative decline of the bar fraction toward lower densities that is most pronounced for the most massive galaxies. In conclusion, our results show that the cluster environment drives a progressive transformation of galaxy morphology, with the loss of disc structure becoming more pronounced towards the cluster core, where environmental processes act most efficiently.
Read the full abstract

All authors P. Mas-Buitrago, S. Kruk, D. O'Ryan, M. T. Nardone, A. La Marca, P. Awad, M. Baes, M. Fossati, T. Géron, Z. Ghaffari, M. Girardi, C. J. Lintott, J. G. Sorce, B. Altieri, F. Durret, V. Pettorino, B. Y. Irureta-Goyena, B. Simmons, J. Reerink, H. Domínguez Sánchez, J. H. Knapen, F. Shankar, E. Duran-Camacho, J. M. Pérez-Martínez, M. Walmsley, N. Auricchio, C. Baccigalupi, M. Baldi, S. Bardelli, P. Battaglia, A. Biviano, M. Bolzonella, A. Bonchi, E. Branchini, M. Brescia, J. Brinchmann, S. Camera, G. Cañas-Herrera, V. Capobianco, C. Carbone, J. Carretero, M. Castellano, G. Castignani, S. Cavuoti, A. Cimatti, C. Colodro-Conde, G. Congedo, C. J. Conselice, L. Conversi, Y. Copin, A. Costille, F. Courbin, H. M. Courtois, M. Cropper, H. Degaudenzi, G. De Lucia, H. Dole, F. Dubath, X. Dupac, S. Dusini, A. Ealet, M. Fabricius, M. Farina, R. Farinelli, S. Ferriol, P. Fosalba, S. Fotopoulou, M. Frailis, E. Franceschi, M. Fumana, S. Galeotta, K. George, B. Gillis, C. Giocoli, P. Gómez-Alvarez, J. Gracia-Carpio, A. Grazian, F. Grupp, S. Gwyn, W. Holmes, I. M. Hook, F. Hormuth, A. Hornstrup, M. Huertas-Company, K. Jahnke, M. Jhabvala, B. Joachimi, S. Kermiche, A. Kiessling, B. Kubik, M. Kümmel, M. Kunz, H. Kurki-Suonio, A. M. C. Le Brun, S. Ligori, P. B. Lilje, V. Lindholm, I. Lloro, G. Mainetti, O. Mansutti, O. Marggraf, M. Martinelli, N. Martinet, F. Marulli, R. J. Massey, S. Maurogordato, E. Medinaceli, M. Melchior, M. Meneghetti, E. Merlin, G. Meylan, A. Mora, M. Moresco, L. Moscardini, E. Munari, C. Neissner, S.-M. Niemi, J. W. Nightingale, C. Padilla, S. Paltani, F. Pasian, K. Pedersen, W. J. Percival, A. Pezzotta, S. Pires, G. Polenta, M. Poncet, L. A. Popa, L. Pozzetti, F. Raison, A. Renzi, J. Rhodes, G. Riccio, E. Romelli, M. Roncarelli, B. Rusholme, R. Saglia, Z. Sakr, A. G. Sánchez, D. Sapone, B. Sartoris, P. Schneider, T. Schrabback, M. Scodeggio, A. Secroun, E. Sihvola, P. Simon, C. Sirignano, G. Sirri, P. Tallada-Crespí, A. N. Taylor, H. I. Teplitz, I. Tereno, N. Tessore, S. Toft, R. Toledo-Moreo, F. Torradeflot, I. Tutusaus, J. Valiviita, T. Vassallo, G. Verdoes Kleijn, A. Veropalumbo, Y. Wang, J. Weller, A. Zacchei, G. Zamorani, F. M. Zerbi, E. Zucca, J. Macias-Perez, M. Sereno

Earth and Planetary Astrophysics

Interplanetary Dust Ablation Seeds Silicate Clouds in Exoplanet Atmospheres

Interplanetary dust particles (IDPs) are present throughout the Solar System, with Earth experiencing a mass flux of roughly 100 tonnes of IDPs per day. These particles ablate in the upper atmosphere, delivering aerosols that form noctilucent and stratospheric clouds. Since IDPs are present in exoplanet systems as exozodiacal dust disks (exozodis), it is likely that exoplanets experience a constant infall of IDPs as well. In this work, we investigate the effect of IDPs on the cloud structure of gas-giant exoplanets and brown dwarfs, and determine their observability. We apply the microphysical cloud model Nimbus to the Sonora Diamondback model grid to cover a wide range of effective temperatures, surface gravities, atmospheric mixing constants, and IDP mass fluxes. Our results show that IDP infall increases cloud particle number densities, decreases cloud particle radii, and enhances silicate features in transmission and thermal emission spectra. Exoplanets with lower atmospheric mixing rates, effective temperatures, and gravity are more sensitive to IDPs. Most notably, they lead to observable silicate absorption features in colder planets where silicate clouds would otherwise not be observable. WASP-107 b is an example of such an exoplanet with observational evidence of silicate particles. Using Nimbus, we found that a 100 times Earth-like IDP mass flux with a Solar-System-like IDP composition of quartz (SiO2), enstatite (MgSiO3), and iron (Fe) can explain the observations of high altitude silicate particles in WASP-107 b.
Read the full abstract

All authors Sven Kiefer, Gregory J. Gosselin, Caroline V. Morley, Sean P. S. Gulick

Solar and Stellar Astrophysics

An analysis of four stellar rings

About 50 years ago, one thousand ring-like structures (called stellar-rings) were discovered by Isserstedt (1968). They were believed to be groups of young stars formed by shell-like triggered star formation, which would make them excellent tracers of spiral arms, for example. Neglected for 40 years, we used highly accurate kinematic, astrometric, and photometric data to investigate the four most prominent stellar rings. The aim is to investigate if those structures are indeed physically related groups of stars. We used proper motions and parallaxes from the Gaia DR2 to calculate distances and to search for common properties. Color-magnitude diagrams using BVJHKs measurements were investigated and isochrones fitted. None of the four stellar rings consists of a physically related group of young stars. The location of stars in the line-of-sight mimics a ring-like structure on the sky. The color-magnitude diagrams are typical for an integrated field population and not for a young star cluster, for example. The currently available data are sufficient to analyze ring-like structures with a high statistical significance. This allows a new search for such structures in the Milky Way.
Read the full abstract

All authors E. Paunzen, J. Florian, A. Gütl-Wallner, A. Herdin, E. Kralofsky, K. Leschinski, M. Mach, H.M. Maitzen, M. Prišegen, M. Rockenbauer, M. Rode-Paunzen, S. Wallner

Cosmology and Nongalactic Astrophysics

Optical Depths from the Thermal Sunyaev-Zel'dovich Effect with ACT DR6 and DESI DR1 Spectroscopic Galaxies and Optically-Selected Clusters

We present stacked thermal Sunyaev-Zel'dovich (tSZ) effect measurements for three samples of galaxy groups and clusters: those traced by the Dark Energy Spectroscopic Instrument Data Release 1 (DESI DR1) luminous red galaxies (LRG) and the DESI DR1 Bright Galaxy Sample (BGS), and an eROMaPPer optically-selected sample from the DESI Legacy Imaging Survey. We use the latest Atacama Cosmology Telescope DR6 (ACT)+Planck component-separated internal linear combination (ILC) Compton-$y$ maps and ACT+Planck coadded 90, 150, and 220 GHz temperature maps to extract the tSZ signal within a $\sim2'$ disk aperture for sources binned by luminosity, richness, or mass. We measure the average tSZ signal with high statistical significance, with signal-to-noise ratios surpassing 38 for LRG, 27 for BGS, and 39 for the eROMaPPer sample using the 90 GHz ACT DR6+Planck map. We conduct a detailed study of systematics and foregrounds such as dust and cosmic infrared background (CIB) contamination, which remain a core challenge for tSZ analysis. For the LRG and BGS samples, we find that dust and radio source emission dominate the tSZ signal at scales near and below the disk aperture radius. Large-scale ($R>4'$) contamination from the CIB is less significant. We mitigate these contaminants to isolate the tSZ signal and use a combination of simulated and real measurements to develop Compton-$y-$optical depth ($\bar y-\bar \tau$) scaling relations to infer optical depths, which are found to be in agreement with values measured using the pairwise kinematic SZ effect for the same tracer samples. The $\bar y-\bar \tau$ scaling relation for the eROMaPPer sample is the first such relationship to be derived directly from SZ measurements.
Read the full abstract

All authors J. E. Moore, C. Popik, Y. Gong, Y-H. Hsu, E. M. Vavagiakis, N. Battaglia, R. Bean, D. Gruen, P. Gallardo, B. Hadzhiyska, G. Ockert, J. Aguilar, S. Ahlen, A. Aviles, F. Beutler, D. Bianchi, J. R. Bond, D. Brooks, E. Bulbul, A. Carnero Rosell, E. Chaussidon, T. Claybaugh, J. Comparat, A. de la Macorra, Biprateep Dey, P. Doel, S. Ferraro, A. Font-Ribera, J. E. Forero-Romero, E. Gaztañaga, Satya Gontcho A Gontcho, G. Gutierrez, K. Honscheid, M. Ishak, S. Juneau, T. Karim, R. Kehoe, M. Kluge, A. Kremin, A. Kusiak, O. Lahav, M. Landriau, L. Le Guillou, M. Manera, A. Meisner, R. Miquel, T. Mroczkowski, S. Nadathur, J. A. Newman, M. D. Niemack, H. E. Noriega, E. Paillas, W. J. Percival, F. Prada, I. Pérez-Ràfols, C. Ravoux, B. Ried Guachalla, G. Rossi, R. Ruggeri, L. Samushia, E. Sanchez, C. Saulder, D. Schlegel, J. Silber, M. Siudek, G. Tarlé, B. A. Weaver, E. J. Wollack, H. Zou

Solar and Stellar Astrophysics

Now you see me, now you don't: Deep, long-lived obscuration events of white dwarfs by dust and debris from disintegrating planetesimals

Transiting dust and debris around white dwarfs provide a direct view of planetary disruption. We report six previously non-variable white dwarfs that underwent single, deep obscuration events, including ZTF~J2327+0019, which fully disappeared for more than a year. These discoveries raise the number of confirmed debris-transit white dwarfs to 20. In addition, the systematic search of negative alerts from the Zwicky Transient Facility (ZTF), followed by photometric and spectroscopic observations, also identified three hot ($T\approx 55000$–$92000$ K) irregular variables as candidate debris hosts. The six confirmed events last weeks to years, reach transit depths of 25–100%, and show no recurrence over the 8-yr ZTF baseline. They are therefore much deeper and longer lived than the recurring, minute-scale transits typical of previously known systems. The observed reddening is consistent with dust extinction, and high-speed photometry reveals minute-scale structure in several systems, indicating rapidly moving material close to the white dwarf. ZTF~J2327+0019 also shows a simultaneous WISE infrared brightening, the first observed coincidence of infrared brightening and optical obscuration in a white dwarf. Two of the six hosts, one DAZ and one DBZ, are strongly metal enriched, with accretion rates among the highest known for their spectral types and minimum accreted masses of $3.5\times 10^{19}$ and $2\times 10^{24}$ g. We interpret the one-off transit events as collisions between planetesimal remnants and pre-existing circumstellar debris that trigger collisional cascades: short-timescale variability traces dust-producing fragments crossing the star, while slower changes trace the build-up and clearing of an optically thick dusty structure. In contrast, ZTF~J0436+5728 shows a coherent 5.910-h period, consistent with a rocky body shedding material near the tidal-disruption radius.
Read the full abstract

All authors J. van Roestel, S. Bhattacharjee, J. A. Guidry, I. Caiazzo, T. Cunningham, J. J. Hermes, S. Torres, Z. Vanderbosch, A. J. Drake, M. Hollands, B. T. Gänsicke, K. De, T. Prince, A. C. Rodriguez, L. B. Ould Rouis, A. A. Mahabal, E. C. Bellm, M. W. Coughlin, M. J. Graham, D. Hale, G. Helou, J. C. Jaimes, M. Kasliwal, R. R. Laher, R. L. Reed

High Energy Astrophysical Phenomena

Resolving structure within the iron line profile of Serpens X-1 with XRISM and NuSTAR

We present the first simultaneous XRISM and NuSTAR observations of the neutron star low mass X-ray binary Serpens X-1. We perform spectral modeling on the joint observations of XRISM/Resolve, XRISM/Xtend, and NuSTAR, testing Comptonization and double thermal continuum model prescriptions. We find that a hybrid double thermal model adequately describes the underlying continuum of the source while the Comptonization models predict largely unphysical parameter values. We perform reflection modeling and confirm with tight constraints an inclination of 5^{\circ} \pm 1^{\circ} and inner disk radius of 6.6 \pm 0.6 Rg for the source. We show that the spectral resolution of XRISM/Resolve allows for the determination of a unique inner disk radius, which has the potential to constrain neutron star spin in other systems. We discuss discrepancies in spectral shape between XRISM/Resolve and NuSTAR above \sim 8 keV. We also present lightcurve and spectral analyses of 8 Type-I X-ray bursts that occurred during the observations.
Read the full abstract

All authors Hayden Hall, Renee M. Ludlam, Edward Cackett, Jon M. Miller, Javier A. García

High Energy Astrophysical Phenomena

Strong Polarization Signatures from Magnetically Stabilized Luminous Thin Accretion Disks

We utilize the Monte Carlo radiation transport capabilities of Cosmos++ to explore the polarization and Faraday rotation of radiation emitted from a set of general relativistic radiation magnetohydrodynamic simulations of magnetically stabilized, thin, black hole accretion disks. The guiding question is whether or not the Faraday rotation depolarizes the radiation to such a degree as to be inconsistent with the relatively high polarization measurements coming from the Imaging X-ray Polarimetry Explorer (IXPE). After first confirming that our code reproduces expected polarization results for electron clouds, planar scattering atmospheres, and AGN disks, we demonstrate the polarization and Faraday rotation capabilities using a Novikov-Thorne disk threaded with a purely toroidal magnetic field. We then analyze temporal snapshots from three different simulations of luminous, thin accretion disks threaded with different magnetic field configurations. We find that the effects of Faraday rotation are generally modest over the energy range of interest, since the strongest fields present in these simulations ($\gtrsim 10^8$ G) are mostly hidden beneath the photosphere. We can easily produce polarization degrees $\ge 4$% for sources seen mostly edge-on. As one example of a relevant target, we compare our results to an IXPE observation of Cyg X-1, finding that we can match the polarization degree quite well, though our polarization angle is rotated approximately $90^\circ$ with respect to the observed data. We speculate on a few possible explanations for this.
Read the full abstract

All authors P. Chris Fragile, Peter Anninos, Madeline Breiling, Derrick Pickrel, Gabrielle Ruminski

Astrophysics of Galaxies

Turbulence Cascade in Cygnus X Revealed by Multi-point VDF Method

Turbulence plays a crucial role in regulating star formation activities within molecular clouds, yet few methods can directly reveal its properties and underlying processes. We use molecular line data from the Nobeyama 45m Cygnus X CO Survey to study the turbulence properties and their relationships with star-forming activities and/or other non-thermal motions. In this work, we apply the multi-point velocity dispersion function (VDF), rather than direct linewidth measurements, to investigate the non-thermal properties of molecular cloud motions. We filter out the large-scale ordered structure and isolate a relatively small-scale turbulence component. Through the Friends In Velocity (FIVe) algorithm, we identify 10 substructures of the clouds and derive the turbulent properties of each cloud using the VDF method. We find that both the cloud-complex regions and the 10 velocity substructures exhibit turbulence correlation lengths of $\sim 2$–5 pc. This plateau scale suggests a parsec-scale turbulence correlation or driving scale in Cygnus X. Below this scale, the rising VDFs trace the velocity scaling of the turbulent cascade, whereas larger-scale VDF variations likely reflect cloud-scale motions. The comparison between cloud complexes and substructures further suggests that, in observational data, the VDF may constrain the turbulence correlation scale more robustly than the turbulence velocity dispersion.
Read the full abstract

All authors Junjie Huang, Yangjun Pu, Keping Qiu, Junhao Liu, Yingxi Li, Mengke Zhao

High Energy Astrophysical Phenomena · Cross-listed from nucl-th

How strange: Phase diagrams with 3 critical points

We show that the Chiral Mean-Field model (CMF) can produce a phase diagram with three critical points: the nuclear liquid-gas transition, quark deconfinement, and a strangeness driven transition. The strangeness driven transition separates a mainly nucleonic phase from one dominated by hyperons and baryon resonances. We discuss the compositional change in these transitions and their possible signatures in heavy-ion collisions and neutron star mergers.
Read the full abstract

All authors Mateus Reinke Pelicer, Nikolas Cruz Camacho, Rajesh Kumar, Veronica Dexheimer, Jacquelyn Noronha-Hostler

Solar and Stellar Astrophysics

The role of specific entropy in the onset of the luminosity bump

During their evolution, stars follow a distinct path in luminosity–temperature space. Low-mass stars that have exhausted hydrogen in their core follow the so-called red-giant branch which is predominantly in the direction of increasing luminosity and decreasing surface temperature. The luminosity bump is a temporary decrease in luminosity against this otherwise increasing trend. The bump is present in both observations and stellar evolution models. However, with canonical physics included in the models, the bump in computed evolutionary tracks appears at higher luminosities than the observed bump. To understand why the bump appears later in the models than in observations, the physics of the bump needs to be unravelled. The end point of the bump is well-understood; however, the onset of the bump is still an enigma. Here, we report on the physical origin of the onset of the luminosity bump. We show that the difference in specific entropy at the mean molecular weight discontinuity decreases due to the discontinuity moving in. This decrease is attributed to the decrease in the ratio of the temperature to the pressure at smaller radii. Upon reaching a critical value the specific entropy difference at the mean molecular weight discontinuity is reduced sufficiently to reduce the specific entropy in the convective envelope. The latter is a key signature of the bump. Hence the evolution of the specific entropy at the mean molecular weight discontinuity provides a viable description for the onset of the bump.
Read the full abstract

High Energy Astrophysical Phenomena · Cross-listed from hep-ph

Axion Hair and Pulsar Electrodynamics: modelling, discharge dynamics, and particle-in-cell simulations

In a companion paper, we demonstrated that static axion field gradients sourced by dense nuclear matter (axion hair) can dominate the near-field electrodynamics of old rotation-powered pulsars, leading to new constraints on light QCD axions and on CP-violating axion-nucleon interactions. This article provides the extended theoretical and numerical framework underlying those results. We begin by providing a detailed description of the sourcing of axion hair from dense nuclear matter, computing self-consistent field profiles for each interaction across the relevant parameter space. We then study the modification to the electrodynamics induced in the polar gap region by these axion gradients; this is done at the analytic level by studying the modification induced by axion field gradients on the effective discharge parameter (computed in the force-free limit of the split monopole magnetic field configuration, and looking at leading deviations from the force-free limit for dipolar field configurations), and numerically by developing dedicated 1D particle-in-cell simulations which capture the leading order dynamical behavior near the star. Our results demonstrate that axion hair serves to either enhance acceleration, or enhance screening, where the relevant effect changes between the northern and southern hemispheres of the star, and between the field lines which support out-flowing and return currents.
Read the full abstract

All authors Samuel J. Witte, Andrea Caputo, Stefan Stelzl, Alexander Chernoglazov, Alexander A. Philippov, Surjeet Rajendran

Cosmology and Nongalactic Astrophysics

Why Azimuthal Averaging Works in Halo Lensing: Symmetry and Power Counting for Nonlinear Shear and Magnification

Cluster weak-lensing analyses often compress two-dimensional lensing fields into azimuthally averaged radial profiles and evaluate nonlinear observables from the averaged convergence and shear. However, the ring average of a nonlinear observable generally differs from the mean-field prediction formed from the averaged fields. This difference, divided by the mean-field prediction, defines the fractional residual. For a complete ring, the difference begins at second order in angular fluctuations. For centered elliptical halos, rotational symmetry makes the shape contribution even in signed ellipticity, while the leading miscentering dipole is orthogonal to the shape quadrupole. We test these predictions using projected triaxial NFW halos at six mass-redshift grid points spanning $3\leq M_{200\mathrm c}/(10^{14}\,h^{-1}M_\odot)\leq20$ and $0.2\leq z_l\leq0.5$, for $z_s=1$. For the reference offset model, centering offsets follow a Rayleigh distribution with scale $0.05r_{200\mathrm c}$. For $0.366\leq R/r_{200\mathrm c}\leq1$, where at least $95\%$ of each population satisfies a conservative subcriticality criterion, we calculate the median of the fractional residuals across the retained halos for each population and radius. The largest population-median magnitudes are $0.90\%$ for reduced shear, $0.0093\%$ for inverse magnification, $0.64\%$ for magnification, and $0.18\%$ across the two magnification-bias cases $\mu^{\alpha-1}$ ($\alpha=0.3, 1.4$). In the paired calculation at $z_s=2$, the maximum magnitude increases for every observable, consistent with weak-lensing power counting. Across both source planes and the common radial domain, every population median remains below $2\%$ in magnitude. This accuracy follows from first-order cancellation, rotational symmetry, and harmonic orthogonality, with further weak-lensing suppression of the remaining nonlinear terms.
Read the full abstract

High Energy Astrophysical Phenomena

Phase-Resolved Ultra-High-Energy Emission Defines an Energetic Boundary for Compact-Object Engines and Dense-Matter Structure

We show that phase-resolved ultra-high-energy emission from LS I \(+61^\circ303\) defines a critical energetic boundary linking the compact-object engine, circumstellar interaction depth, and dense-matter structure. Combining the GeV–TeV–UHE spectral transition with an orbital acceleration–transparency sieve and branch-resolved stellar sequences for 116 certified EoS, we find that the measured \(25\)–\(100\) TeV component remains energetically viable for all 116 canonical-mass EoS and therefore does not discriminate among them. EoS discrimination arises only in the deliberately more demanding pure-hadronic \(>1\) TeV stress-test construction: at 2.0 kpc and \(1.4\,M_\odot\), this boundary gives \(I_{45,\rm crit}=1.5504\), with 82 of 116 EoS satisfying the energetic condition at\(\dot P_{\rm ref}=3\times10^{-15}\,{\rm s\,s^{-1}}\). A separately applied tidal diagnostic, \(\Lambda_{1.4}\leq580\), is satisfied by 43 EoS, with nine satisfying both conditions. These models span several microscopic classes but occupy a narrow structural region, indicating structural overlap rather than selection of a unique EoS family. A separate GeV power-budget test requires \(\dot P_{\rm rot}\gtrsim 5.65\times10^{-14}f_\Omega\,{\rm s\,s^{-1}}\) at 2.0 kpc. For \(f_\Omega\sim1\), the source lies well above the UHE crossing and the EoS discrimination disappears. The UHE observation therefore identifies the engine–environment regime in which differences in neutron-star structure can become observationally discriminating, rather than defining a universal dense-matter selector.
Read the full abstract

High Energy Astrophysical Phenomena

Enabling new sampling strategies for continuous-wave analyses by integrating Bilby into PyFstat

Continuous-wave searches for rotating neutron stars usually cover wide parameter spaces and produce a large number of candidate signals. Stochastic sampling methods therefore play an important role in candidate follow-up and parameter estimation. We present the integration of the Bayesian inference library BILBY into the open-source continuous-wave analysis package PyFstat, enabling access to the broad range of stochastic samplers available through BILBY while preserving the existing PyFstat analysis infrastructure. The implementation is validated through injection studies in Gaussian noise for two representative single-stage follow-up configurations: candidates from a directed search and candidates from an all-sky search. Using the nested sampler DYNESTY, we obtain detection efficiencies consistent with the theoretical sensitivity predictions. In this setup, DYNESTY performs comparably to the existing PTEMCEE-based implementation of PyFstat for the directed case, and also provides an effective single-stage approach for the all-sky candidate follow-up considered here, where we did not identify a PTEMCEE setup with similar recovery performance. Percentile-percentile tests for DYNESTY further show well-calibrated credible intervals, demonstrating reliable parameter estimation. While a full exploration of DYNESTY settings, multi-stage setups, and other samplers is left for future work, these results demonstrate that this flexible and publicly available framework opens useful new possibilities for Bayesian continuous-wave analyses.
Read the full abstract

High Energy Astrophysical Phenomena

Compton-induced HE – VHE $\gamma$-rays from pair cascade emissions in the SEDs of radio galaxies: NGC 1275

In the scheme of active galactic nuclei (AGNi), the blazar subclass is particularly bright in $\gamma$-rays. This is because of their emissions being beamed into our line of sight, as opposed to their misaligned parent population, radio galaxies. This work presents results of a Monte-Carlo code that propagates jet-collimated $\gamma$-rays through the magnetized AGN environment, leading to secondary cascade emissions from relativistic electron-positron pairs. The resulting multi-wavelength emissions are used to fit the broadband spectral energy distribution (SED) of the radio galaxy NGC 1275. Using known physical parameters of NGC 1275 from the literature, we find that cascade emissions from stochastic processes have the potential to reproduce the broadband SEDs of radio galaxies. We also find that the interplay between accretion disk and broad-line region (BLR) photons may provide insights into the production region(s) of $\gamma$-rays in AGNi. While dense accretion disk and BLR radiation fields are vital to the development of cascades, excessive densities lead to the suppression of cascade synchrotron emission and attenuation of \,TeV $\gamma$-rays. Furthermore, the magnetic field is important in shaping the cascade $\gamma$-ray spectrum.
Read the full abstract

Astrophysics of Galaxies

Design and commissioning of a windowless gas-target system for high-current beams at JUNA

Windowless gas targets avoid the beam-energy loss and straggling introduced by entrance foils and are therefore well suited for direct measurements of low-energy nuclear reactions. A windowless gas-target system designed for operation with milliampere beams has been developed for the Jinping Underground Nuclear Astrophysics facility (JUNA). The system combines three-stage differential pumping, closed-loop gas recovery and purification, a constant-temperature power-compensation calorimeter, and a position-resolved target-thickness monitor based on secondary elastic scattering. Stable operation was achieved over a target-pressure range of 1-3 mbar, with pressure fluctuations below 1% during 8 h of continuous circulation, while the accelerator-side pressure was maintained at approximately \(10^{-4}\) Pa. The closed-loop gas-circulation system maintained stable target conditions, while gas-transport calculations indicated that the axial pressure nonuniformity remained within approximately 1.6% under representative operating conditions. Calorimeter measurements were consistent with the thermal calculations, supporting the sensitivity correction used for beam-power determination. Beam commissioning with \(^{14}\mathrm{N}(p,\gamma)^{15}\mathrm{O}\) and \(^{12}\mathrm{C}(p,\gamma)^{13}\mathrm{N}\) at the 600 kV Cockcroft-Walton accelerator of the China Institute of Atomic Energy (CIAE) demonstrated stable operation of the gas-target and \(\gamma\)-ray detection systems and provided information on the influences of reaction position and beam heating. These results demonstrate the operating stability and diagnostic capability of the system for future high-current, low-energy nuclear-reaction measurements at JUNA.
Read the full abstract

All authors Xuyang Wang, Fuqiang Cao, Weike Nan, Gang Lian, Wei Nan, Hongwei Huang, Yuchen Jiang, Jingyu Dong, Minghao Zhu, Yuwen Chen, Yihui Liu, Jinlong Ma, Yueran Wang, Jie Chen, Yunju Li, Shengquan Yan, Youbao Wang, Xiao Fang, Luohuan Wang, Yangping Shen, Bing Guo, Weiping Liu, JUNA collaboration

Astrophysics of Galaxies

Weak cool core, weak turbulence: XRISM, XMM-Newton and Chandra spectroscopy of the core of Abell 1413

We present the first XRISM spectroscopy of the core of Abell 1413 (z = 0.143), combining 104.7 ks of Resolve with archival XMM-Newton and Chandra spectra of a matched sky region. A1413 is hot, massive and relaxed but lacks a strong cool core (t_cool = 6.4 Gyr), a regime in which the low turbulence reported to date has been tested only once, in A3571. The Resolve spectrum requires a single velocity-broadened component with kT = 7.06 (+0.33, -0.43) keV and a line-of-sight velocity dispersion sigma_v = 170 (+25, -24) km/s (M_3D = 0.217). The bulk velocity relative to the brightest cluster galaxy, v_bulk = -24 +/- 26 (z_ICM) +/- 71 (z_BCG) km/s, is limited by the optical redshift rather than by Resolve, and disfavours coherent line-of-sight sloshing. The non-thermal pressure fraction, 2.5 +/- 0.7 per cent, implies a hydrostatic mass bias b < 2.5 per cent, comparable to the strong cool-core cluster A2029 (2.6 +/- 0.3 per cent) and, at matched sampling scale, if anything below it. With A3571 (1.3 +/- 0.3 per cent), to which A1413 rescales within 0.2 sigma, two weak cool cores now show that this quiescence does not require central cooling. The resolved Fe-K He-alpha complex gives Fe = 0.38 +/- 0.03 (stat) +/- 0.04 (sys) solar; Ni/Fe, separated from Fe K-beta for the first time, is solar, as is the pattern as a whole. The Ne, Mg and Si lines that set Ne/Fe, Mg/Fe and Si/Fe lie below 2 keV, where our instruments disagree by 10-20 per cent; we therefore do not regard the SNIa fraction as constrained.
Read the full abstract

Astrophysics of Galaxies

Detection prospects for heavy WIMP dark matter around M31* in microwave band

This work analyzes the detection prospects for weakly interacting massive particles (WIMPs) in dark matter (DM) density spike around the supermassive black hole (SMBH) in Andromeda galaxy M31 in microwave band. WIMP annihilation produces relativistic electrons and positrons, which generate the synchrotron emission in strong magnetic field of the accretion flow medium at microwave frequencies. Properties of this emission, as well as their dependence on DM and medium parameters, were modeled. Atacama Large Millimeter Array (ALMA) sensitivity to this emission and underlying WIMP parameters was estimated. It was obtained that M31* in microwave band represents the unique target, where DM spike signal can be detected in case of realistic spike and halo density profiles. In such case ALMA may reach the thermal s-wave annihilating WIMPs with masses up to around 1 TeV with just 10–100 hours exposure. Optimistic parameter configurations may allow probing all possible WIMP masses up to 100 TeV. Related systematic uncertainties are discussed.
Read the full abstract

High Energy Astrophysical Phenomena

A Two-Parameter Framework for the Diffuse Supernova Neutrino Background

Following important experimental developments, tests of the Diffuse Supernova Neutrino Background (DSNB) are entering a new phase that might lead to a discovery. This motivates the development of simple and broadly applicable tools for comparing theoretical predictions and interpreting experimental results. We show that, over the energy range relevant to current detectors, the DSNB is well described by an exponential spectrum, characterized by the differential flux at 20 MeV, $\phi_{20}$, and a spectral energy scale, $E_0$. We fit the spectra of 24 representative models and demonstrate that the exponential description provides an excellent approximation. We then recast the recent Super-Kamiokande results in the $(\phi_{20},E_0)$ plane, showing that model-specific best fits and upper bounds can be interpreted naturally within this framework. Our results indicate that, above realistic detection thresholds, and until the precision phase is reached, the observable DSNB is effectively a two-parameter phenomenon. We advocate that the $(\phi_{20},E_0)$ parameterization be adopted as a standard, to be used routinely to report future theoretical predictions and experimental results.
Read the full abstract

High Energy Astrophysical Phenomena

A Unified Timescale Relation for Quasi-Periodic Eruptions and Repeated Nuclear Transients

Quasi-periodic eruptions (QPEs) and recurrent nuclear transients (RNTs) exhibit recurrent high-amplitude flares from galactic nuclei, yet their characteristic timescales remain poorly understood. In this work, we compile a sample of these systems and investigate empirical scaling relations between flare timescales and black hole masses. We find that the recurrence timescale exhibits a positive but highly scattered dependence on black hole mass in the combined QPE and RNT sample, approximately following $t_{\rm rec}\propto M_{\rm BH}^{1.19^{+0.52}_{-0.47}}$ with an intrinsic scatter of 0.95 dex. Remarkably, we uncover a tight nearly linear relation between recurrence time and flare timescale for QPEs and RNTs, described by $t_{\rm rec}\propto t_{\rm rise}^{1.00\pm0.07}$ with an intrinsic scatter of 0.29 dex. This relation extends across timescales from hours for QPEs to months and years for nuclear transients. We further find that QPEs and RNTs approximately follow a common empirical relation between recurrence time and flare rise time, although the physical origin of this relation remains uncertain. Our findings reveal a common phenomenological timescale link across RNTs, providing a practical framework for characterizing their temporal behavior.
Read the full abstract

All authors Shifeng Huang, Tinggui Wang, Ning Jiang, Yibo Wang, Zhenfeng Sheng, Tian-Yu Xia, Jiazheng Zhu, Zheyu Lin, Jie Lin, Ji-an Jiang, Kenta Taguchi, Keiichi Maeda

High Energy Astrophysical Phenomena

Stratified Wind Geometries of GRS 1915+105: Density Profiles and Radii

The low-mass X-ray binary GRS 1915+105 exhibits dramatic spectral and flux variability, often accompanied by highly ionized outflows. In this work, we analyze four archival Chandra/HETG observations spanning over a decade, which vary greatly in their ionizing spectrum from soft (low inner disk temperature) to hard (high inner disk temperature). We measure the wind properties, using the insightful tool of the absorption measure distribution (AMD, a distribution of column density with ionization). While the soft epoch reveals a steep, monotonically rising AMD, hard epochs feature a distinct turnover at high ionization parameter ($\log\xi>3.3$). This is the first time we observe a negative AMD slope in an outflow, which we interpret as a localized density inversion, which could point to radiation-pressure compression zones. We utilize the AMD shape and total column density measurements to calculate the density profile and distance of the wind from the central source. For the hard states, we obtain densities of $10^{12}-10^{13}$ cm$^{-3}$ at radii of $10^{10}-10^{11}$ cm. Such measurements are usually not directly possible in X-ray spectroscopy. Here we constrain the inner radius and density of the outflow to a factor of two. The high-$\xi$ values up to $\log\xi\sim6.5$ with no turnover in the soft state imply a smaller radius of a few $10^{8}$cm, which is in tension with the slow velocities observed. We demonstrate that this is a generic tension between compact launch radii and sub-escape velocities of high-column, high-ionization outflows observed in stellar-mass black hole binaries.
Read the full abstract

All authors Noa Keshet, Ehud Behar, Jon M. Miller, Joey Neilsen

Earth and Planetary Astrophysics

Centaur Longevity Revisited: Lifetime Dispersion and the Bailey & Malhotra Classifications under Modern Orbital Solutions

Centaurs occupy dynamically unstable orbits between Jupiter and Neptune and are typically removed by giant-planet encounters within a few Myr. Bailey and Malhotra (hereafter BM09) classified the known Centaurs, from 2007 orbits, into a short-lived diffusing class (D), a resonance-hopping class (R), and a long-lived, quasi-stable class (Q); their sole Q object, 2005 TH173, was fit to a 17\,d arc. We revisit both halves of that analysis on modern Jet Propulsion Laboratory Small-Body Database (SBDB) orbits (retrieved 2026-07-26), integrating clone ensembles of each object with the REBOUND WHFast $N$-body integrator. We ask (i) whether tighter orbital precision narrows the dispersion in clone lifetimes, and (ii) whether the BM09 labels survive on updated elements. Among 65 genuine Centaurs with arc $\ge 30$\,d and no censored clones, $\sigma(\log_{10}\mathrm{lifetime})=0.54$\,dex (a factor of 3.5). Conditioned on present-day $a$, $e$, $i$, and $q$, this dispersion is insensitive to relative semimajor-axis uncertainty (slope $=-0.002$\,dex per decade in $\sigma_a/a$; 95% confidence interval $-0.034$ to $+0.030$) — a chaos floor. On the D/Q axis, 49/50 comparable objects retain their BM09 class; 2005 TH173 is the exception (Q$\to$D): 80/100 clones escape before 10\,Myr (median 2.99\,Myr), and a 40\,Myr integration of the same initial conditions returns 88\,D / 12\,R / 0\,Q. None of the 30 multi-opposition 2007 orbits changes class; two short-arc objects leave the Centaur zone.}
Read the full abstract

High Energy Astrophysical Phenomena

Supermassive black hole binaries in the multi-messenger context of ground and spaceborne VLBI

Evidence of a gravitational wave background suggests the existence of a population of sub-parsec supermassive black hole binaries (SMBHBs), with characteristic angular separations on the order of 1-10 $\mu$as. Spaceborne extensions of Very Long Baseline Interferometry (VLBI) and the next generation ground arrays introduce the possibility of directly imaging these SMBHBs. In this work, a binary spectral energy distribution model is used to predict the detectability of SMBHBs with ground and spaceborne VLBI. We consider the Black Hole Explorer (BHEX), a proposed spaceborne VLBI mission, as our primary case study. We explore the detectable SMBHB parameter space and identify distinguishable binary signatures in the visibility domain. We find that for a flux-density-limited sample, ground array observations are more effective at detecting a wider region of the binary parameter space, with $M_\mathrm{tot} \gtrsim 10^9$ solar mass systems detectable out to redshift $z=0.075$ and beyond, depending on mass ratio. Conversely, inclusion of a spaceborne element such as BHEX, offering finer angular resolution ($\sim6\,\mu$as) and sampling of the (u,v) plane not limited by Earth rotation synthesis, will provide significant benefits in constraining binary properties. An illustrative Fisher analysis shows improvements in characterisation of the separation and position angle of SMBHBs by a factor of $\sim4$. Near-future ground and/or spaceborne VLBI may achieve the first direct observation of a SMBHB, contributing significantly to multi-messenger studies of such systems with pulsar timing arrays and observations across the electromagnetic spectrum.
Read the full abstract

All authors B. Hudson, L. I. Gurvits, D. J. D'Orazio, C. Tiede, T. M. Eubanks, E. Mooij

Instrumentation and Methods for Astrophysics

Full Characterisation of the Polarisation Primary Beam of the GRAO 32-m Telescope

Direction-dependent instrumental polarisation is a major systematic limitation in high-fidelity single-dish radio polarimetry, yet a unified characterisation of beam leakage and the conditioning of polarisation recovery is lacking for the Ghana Radio Astronomy Observatory (GRAO) 32-m telescope. We aim to establish a quantitative, direction-dependent polarimetric beam model at 5.0 and 6.7 GHz. High-resolution GRASP electromagnetic simulations of the nominal telescope configuration are used to derive the full complex Jones response and corresponding Mueller matrices, from which beam shape, beam squint, instrumental Stokes leakage, and intrinsic cross-polarisation ratio (IXR) are evaluated across the primary beam. The Stokes $I$ beams have half-power beam widths of $396.4$ and $295.3$ arcsec at 5.0 and 6.7 GHz, respectively, with main-beam efficiencies of $56.4$ and $55.6$ per cent. At 5.0 GHz, the circular-polarisation beam squint is $9.96$ arcsec ($2.5$ per cent of the HPBW), whereas no statistically significant squint is detected at 6.7 GHz. Leakage from Stokes $I$ into linear polarisation remains below $0.1$ per cent within the half-power beam but increases substantially towards the sidelobes. IXR reaches approximately $80$ dB on axis and decreases with angular offset. These results establish the intrinsic electromagnetic polarimetric response of the GRAO 32-m telescope and provide a quantitative baseline for direction-dependent calibration and subsequent observational validation.
Read the full abstract

All authors Theophilus Ansah-Narh, Nia Imara, Benedicta Woode, Joseph Bremang Tandoh, Emmanuel Proven Adzri, Diana Klutse, Evaristus Uzochukwu Iyida

Instrumentation and Methods for Astrophysics

A Content-Addressed Workflow for Reproducible DANTE Gravitational-Wave Anomaly Analysis

Unsupervised analysis of gravitational-wave detector data combines expensive scientific stages with data acquisition, provenance checks, statistical calibration, follow-up products, and reporting. Reproducing such an analysis requires more than preserving model weights or a final candidate table: the execution graph, scientific contracts, retry semantics, and verification boundary must also be explicit. We present a local workflow layer for the Domain-Adaptive Network for Transient Evaluation (DANTE) that represents a corrected O4a analysis as a content-addressed, 15-stage directed acyclic graph. The architecture separates immutable stage evidence from mutable operational progress, binds command-line and browser interfaces to the same run identity, and resumes only work compatible with the frozen contract. The tagged v1 release records all 15 stages as verified in a machine-readable receipt and passed a documented human usability gate. The release verifies adopted scientific artifacts rather than recomputing them, and therefore makes no new claim of global significance, astrophysical discovery, or public real-time operation. This paper describes the architecture, its fail-closed controls, and the bounded evidence supporting reproducible local use.
Read the full abstract

High Energy Astrophysical Phenomena

Constraints on the Low-frequency Radio Emission of the Galactic FRB Source SGR 1935+2154

We present a search for radio pulses from the Galactic magnetar SGR 1935+2154, a well-known source of fast radio bursts (FRBs), at $\sim$110 MHz using the Large Phased Array (LPA) of the Pushchino Radio Astronomy Observatory. Data from two active periods in 2020 (March – May and September – November, with $\sim 3.5$ minutes of daily coverage) were analyzed with new methods tailored to both FRB-like single pulses and pulsar-like periodic signals. No significant FRB-like pulses were found. Using Monte Carlo simulations, $3\sigma$ upper limits were derived for the burst rate: for a log-normal energy distribution the limit is $\sim$${10}^{1.5}~{\rm{d}}^{-1}$ for a mean of average monochromatic isotropic luminosity $L_{\nu{\rm ,mean}}\sim1.3\times{10}^{29}~{\rm{erg~s^{-1}~ {Hz}^{-1}}}$ and a natural log-space scatter of $\sigma\sim0.85$; while for a power-law distribution it is $\sim$${10}^{1.8}~{\rm{d}}^{-1}$ for an index $\beta\lesssim3.0$ and a minimum average monochromatic isotropic luminosity $L_{\nu{\rm{,min}}}\lesssim0.7\times{10}^{25}~{\rm{erg~s^{-1}~{Hz}^{-1}}}$. When folded at the known 3.24781628 s period of SGR 1935+2154, a weak pulse was noted (S/N $<$ 3.16), but the significance is insufficient for a secure detection of the pulsar-like emission signal. A conservative upper limit on the average monochromatic isotropic luminosity of any possible periodic emission is $2.08\times{10}^{19}~{\rm{erg~s^{-1}~{Hz}^{-1}}}$. Our results offer meaningful low-frequency upper limits on the burst rate of SGR 1935+2154, and hint for very faint pulsar-like radiation at meter wavelengths.
Read the full abstract

All authors Chen-Ran Hu, Jinhuang Cao, S.A. Tyul'bashev, Pei Wang, Yong-Feng Huang, E.A. Brylyakova, G.E. Tyul'basheva, Jin-Jun Geng, Orkash Amat, Ze-Cheng Zou, Chen Du, Nurimangul Nurmamat, Lang Cui, Fan Xu, Xiao-Fei Dong, Chen Deng