⚛️ phenomenology

Transient PULSE: Code for the Transient Propagation of UltraLight (pseudo)Scalars to Earth

The paper introduces Transient PULSE, a software package that numerically models the propagation of transient ultralight bosonic signals from astrophysical bursts to Earth by discretizing source spectra, simulating their evolution through matter density profiles, and calculating time-frequency energy distributions while accounting for dispersion and screening effects.

Original authors: Jason Arakawa, Charles Cheung, Muhammad Hani Zaheer, Volodymyr Takhistov, Marianna Safronova

Published 2026-10-15

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Jason Arakawa, Charles Cheung, Muhammad Hani Zaheer, Volodymyr Takhistov, Marianna Safronova

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The universe is filled with invisible substances that scientists call dark matter, a mysterious substance that outweighs all the stars and planets combined. While we know it exists because of how it pulls on galaxies, we have never directly seen a single particle of it. For decades, the leading theory has been that this dark matter is cold and slow-moving, drifting quietly through our solar system like a silent fog. However, there is another possibility that has gained attention: that dark matter might be made of ultra-light particles that behave more like waves than solid objects. These particles are so light that a single one would be billions of times lighter than an electron. If they exist, they could be created in violent cosmic events, such as the collapse of a star or the explosion of a dense cluster of particles, and then race across the galaxy toward Earth. The challenge for scientists is that if these particles do arrive, they would not hit our detectors all at once. Instead, because they travel at slightly different speeds depending on their energy, a single burst of particles would stretch out over time, arriving as a long, shifting signal rather than a sharp flash.

A team of researchers has now released a new computer tool designed to predict exactly how these stretched-out signals would look when they reach Earth. This software, named Transient PULSE, acts as a bridge between the violent explosions happening in deep space and the sensitive instruments waiting for them here. The program takes a description of how a cosmic event emits these ultra-light particles and simulates their journey across the vast distances of space. It accounts for the fact that the particles must travel through clouds of gas and dust, which can slow them down or even block some of them entirely. By breaking the signal down into its individual components and tracking how each one moves, the software constructs a detailed map of what the signal would look like to a detector on our planet. This map shows how the signal's frequency changes over time, revealing a specific pattern where the highest energy particles arrive first, followed by a steady drift toward lower energies, creating a distinct signature that astronomers can look for.

The researchers tested their tool using several different scenarios to ensure it works correctly. They simulated signals from idealized bursts and compared them against a more complex, realistic model of an explosion from a dense cluster of particles known as a boson star. In one example, they modeled a burst occurring about eight thousand light-years away, a distance similar to the size of our own galaxy. The simulation showed that while the original burst might last for a few days, the signal stretching out by the time it reaches Earth could last for nearly two weeks. In another test involving a closer event just a few light-years away, the signal stretched out over a period of years. These results confirm that the software can accurately handle the complex physics of how these particles spread out as they travel. The tool also checks for "screening," a phenomenon where the density of matter in space acts like a filter, preventing certain types of particles from passing through at all. If the particles are too heavy or the matter they travel through is too dense, the signal simply stops, and the software correctly identifies which parts of the signal would be lost.

Beyond just showing what the signal looks like, the software helps scientists understand what they need to find it. It calculates how strong the connection between these particles and normal matter must be for a detector to notice them. This allows researchers to look at their current instruments and determine if they are sensitive enough to catch such a signal, or if they need to build new, more powerful detectors. The program can also predict how much time would pass between the arrival of a gravitational wave, which travels at the speed of light, and the arrival of these ultra-light particles. This time difference is a crucial clue; if two signals from the same cosmic event arrive at different times, it would provide strong evidence for the existence of these particles. The researchers have made their code available to the public, allowing other scientists to plug in their own data and explore different types of cosmic events. By turning complex equations into a visual map of time and energy, this work gives the scientific community a practical way to hunt for a new kind of dark matter that has so far remained hidden in the noise of the universe.

Technical Summary: Transient PULSE

Problem Statement

The detection of ultralight bosonic (ULB) fields, which are candidates for dark matter and potential new physics beyond the Standard Model, has traditionally focused on non-relativistic, persistent local backgrounds. However, astrophysical events (e.g., boson star collapses, black hole superradiance, supernovae) can emit energetic, transient bursts of relativistic ULBs. Detecting these bursts requires understanding how the emitted spectrum is reshaped during propagation to Earth.

Key challenges in modeling these signals include:

Dispersion: Different momentum modes travel at different velocities, causing the signal to spread in time. High-momentum modes arrive first, creating an "anti-chirp" signal where frequency decreases over time.
In-Medium Effects: Propagation through matter (e.g., interstellar gas, Earth's atmosphere) introduces effective mass terms that can slow waves or completely screen (block) low-momentum modes if the effective mass exceeds the wave energy.
Non-Gaussian Spectra: Realistic astrophysical sources often produce complex, non-Gaussian emission spectra that are difficult to model analytically.
Signal Reconstruction: Connecting the source emission spectrum to the time-frequency energy density observable at a terrestrial detector requires numerically tracking the arrival time and energy contribution of each momentum mode through a varying density profile.

Methodology

The paper introduces Transient PULSE, a Python-based software package that numerically implements a source-to-detector framework for transient ULB signals. The methodology is divided into two complementary workflows:

1. Numerical Propagation Workflow

This workflow simulates the full propagation of a signal from source to detector:

Input: Users define a source emission spectrum ($dE/dp$) via SpectrumSource objects (supporting CSV data, analytical functions, or time-varying signals).
Discretization: The continuous spectrum is discretized into momentum modes ( $p_i$ ).
Propagation: For each mode, the code numerically integrates the inverse group velocity along the line-of-sight through a specified matter density profile ( $\rho_{matter}$ $ρ_{ma tt er}$ ). The group velocity $v_\phi$ $v_{ϕ}$ is determined by the dispersion relation, accounting for the particle mass $m$ $m$ and in-medium contributions $\beta(\ell) \propto \rho_{matter}$ $β (ℓ) \propto ρ_{ma tt er}$ .
- If $m^2 + \beta(\ell) \geq \omega_p^2$ at any point, the mode is screened and excluded.
Signal Construction: The arrival time $t_\phi(p_i)$ for each mode is calculated. The code constructs a time-frequency energy density distribution (spectrogram) at Earth by summing the contributions of all modes, accounting for the "anti-chirp" effect where frequency drifts downward as lower-momentum modes arrive later.
Output: The primary output is a spectrogram (energy density vs. time and frequency) and integrated projections (energy density vs. time, and time-integrated frequency distribution).

2. Projected Sensitivity Workflow

This workflow provides analytical estimates for detection capabilities without running full numerical simulations:

Inputs: Users specify source parameters (total energy, mass, burst duration, distance) and detector capabilities (integration time, frequency sensitivity, time-delay thresholds).
Calculations: The code computes:
- Detection threshold couplings (rescaling existing Dark Matter sensitivities).
- Time-delay couplings (required to produce measurable delays relative to gravitational waves).
- Critical screening couplings (where environmental screening suppresses the signal).
- Astrophysical constraints and uncertainty principle limits.
Output: Visual projections of viable parameter space in coupling strength vs. energy/mass plots, including exclusion regions and sensitivity curves.

Key Contributions

Software Implementation: The release of Transient PULSE, a modular Python package that bridges the gap between theoretical source models and terrestrial detector observables.
Numerical Treatment of Dispersion and Screening: The code explicitly handles non-Gaussian source spectra and spatially varying propagation media, including the complex interplay of dispersion and in-medium screening which can alter signal duration and spectral content.
Time-Dependent Source Modeling: The framework supports time-varying emission amplitudes, allowing for the simulation of complex burst structures beyond simple static bursts.
Bridge Functionality: The package includes "bridge functions" that allow users to export parameters from a detailed numerical propagation run and feed them directly into the analytical sensitivity workflow, enabling a seamless transition from signal simulation to experimental feasibility studies.
Validation: The code is validated against energy conservation principles (showing <1% deviation in energy integrals for tested resolutions) and analytical limits for propagation time in constant density media.

Results and Examples

The paper demonstrates the package's utility through several examples:

Gaussian Signals: Simulations of time-independent and time-dependent Gaussian bursts illustrate the formation of the "anti-chirp" spectrogram, where high frequencies arrive first and the signal duration at Earth ( $\tilde{t}^*$ ) significantly exceeds the intrinsic source duration ( $t^*$ ).
Screening Effects: A scenario with increased coupling strength demonstrates how specific momentum modes are screened out, altering the arrival time distribution and reducing the total signal energy reaching the detector.
Bosenova Spectrum: The code is applied to a realistic bosenova emission spectrum (derived from numerical simulations of boson star collapse). The results show a multi-peak structure and a very long signal duration (order of years) due to wave spreading, consistent with analytical approximations.
Convergence: Tests confirm that the numerical results converge as the number of time and momentum grid points increases, with energy conservation maintained to high precision.

Significance and Claims

The paper positions Transient PULSE as a necessary tool for the emerging field of transient ULB detection. Its significance lies in:

Enabling Direct Numerical Treatment: It allows researchers to move beyond idealized analytical approximations to handle realistic, non-Gaussian source spectra and complex propagation environments.
Facilitating Multi-Messenger Astronomy: By calculating precise arrival time delays relative to gravitational waves, the tool aids in identifying coincident signals from transient astrophysical events.
Guiding Experimental Searches: The sensitivity workflow provides a rapid method to project detection prospects for various coupling types (photon, electron, gluon, etc.) and experimental configurations, helping to identify viable parameter spaces for future searches.
Complementing Local DM Searches: It offers a complementary approach to persistent local dark matter searches by probing fields that may be suppressed locally but produced in energetic astrophysical events.

The authors emphasize that the detection treatment is a "simplified rescaling of DM sensitivities" rather than a full experimental data-analysis pipeline, and that the code serves as a generator of inputs for detector-specific analyses.