Here we present results from an experiment performed at the GSI Helmholtz Centre for Heavy Ion Research. A mono-energetic beam of chromium ions with initial energies of ∼ 450 MeV was fired through a magnetized interaction region formed by the collision of two counter-propagating laser-ablated plasma jets. While laser interferometry revealed the absence of strong fluid-scale turbulence, acceleration and diffusion of the beam ions was driven by wave-particle interactions. A possible mechanism is particle acceleration by electrostatic, short scale length kinetic turbulence, such as the lower-hybrid drift instability.

We present an extension of the OSIRIS particle-in-cell (PIC) code that introduces an axion macroparticle species and three axion-production channels commonly used in thermal-plasma axion phenomenology: screened Primakoff conversion (γ+ Z ↔ a + Z), Compton-like photoproduction on electrons in a blackbody photon bath (γ+ e → a + e), and thermal axion bremsstrahlung from electron-ion and electron-electron scattering (e + Z → e + Z + a and e + e → e + e + a). The package is integrated into the existing OSIRIS quantum-electrodynamics (QED) Monte Carlo infrastructure and provides Poisson macro-event sampling with unbiased weight rescaling for variance control. Optional modules implement conservative cell-local energy and momentum feedback and temperature-field evolution, and each channel includes an inverse absorption operator constructed to satisfy detailed balance with a thermal bath. We benchmark forward spectral emissivities for uniform plasmas at T_e = 1.3 keV, 3 keV, and 5 keV against screened analytic results based on Raffelt-style calculations, finding percent-level agreement in integrated power for all channels and good reproduction of spectral peak positions. In addition, homogeneous relaxation tests with forward and inverse operators enabled show that, for all three implemented channels, the axion population and total axion energy evolve toward stable steady-state values, providing an initial validation of detailed-balance recovery in the inverse-process implementation. These results establish a foundation for kinetic simulations of axion production, absorption, and transport in high-energy-density plasmas, while more extensive validation of feedback physics and fully dynamic multidimensional coupled scenarios remains future work.

A robust and flexible architecture capable of providing real-time analysis on diagnostic data is of crucial importance to physics experiments. In this paper, we present such an online framework, used in June 2025 as part of the HRMT-68 experiment, performed at the HiRadMat facility at CERN, using the Super Proton Synchrotron (SPS) beam line. HRMT-68 was a fixed-target laboratory astrophysics experiment aiming to identify plasma instabilities generated by a relativistic electron-positron beam during traversal of an argon plasma. This framework was essential for experimental data acquisition and analysis, and can be adapted for a broad range of similar-scale experiments with a variety of experimental diagnostics, even those without a standard direct network communication interface. The developed framework's customizable design enabled us to rapidly observe and extract emergent features from a diverse range of diagnostic data. Simultaneously, its modularity allowed for a quick introduction of new diagnostic devices and the modification of our analysis as features of interest were identified. As a result, we were able to effectively diagnose equipment malfunction, and infer the beam's response to varying bunch duration, beam intensity, and the plasma state without resorting to offline analysis, at which time adjustment or improvement would have been impossible. We present the features of this agile framework, whose codebase we have made publicly available so that it may be adapted for future experiments with minimal modification.

When an intense laser pulse strikes a flat solid target at normal incidence, only odd, equally spaced harmonics are typically generated. Using particle-in-cell simulations, we show that introducing transverse spatial modulation—either via two beams incident at small opposing angles or via a single beam reflecting from a corrugated surface—produces a two-dimensional harmonic spectrum. Selecting the emission through a narrow slit isolates an angle-dependent subset of harmonics, forming a tunable frequency comb controlled by both the viewing angle and the drive-laser frequency. Distinct combs arise for s- and p-polarized emission, while increasing the focal spot size reduces divergence and improves efficiency. Additional control is achieved by varying beam frequencies, target tilt, or beam number. Replacing the dual-beam configuration with a single beam and structured surface provides equivalent functionality with significantly reduced experimental complexity.

We report on a measurement of collective behavior in a relativistic electron-positron pair plasma produced in the laboratory. Using the Fireball platform at CERN's HiRadMat facility, 440 GeV protons were used to generate an ultra-relativistic, charge-neutral electron-positron pair beam that propagated through an ambient RF discharge plasma. Magnetic-field amplification due to a beam-plasma instability was diagnosed using a high-sensitivity Faraday-rotation probe, supported by detailed characterization of the diagnostic impulse response. The measured path-integrated magnetic field agrees quantitatively with predictions from particle-in-cell simulations. The results provide a critical benchmark for models of relativistic beam-plasma interactions in astrophysical contexts such as blazar jets and pulsar-wind nebulae.

Axions are hypothetical particles, proposed to account for the invariance of CP symmetry in quantum chromodynamics. While axions and axion-like-particles are well-motivated by string theory and beyond-Standard-Model extensions, they have remained elusive to experimental searches even after significant effort over many decades. Building on a recent development using an X-ray free electron laser to search for cosmologically favoured axions of mass m_a≲ 0.01 eV, we extend previous bounds on the ALP-photon coupling, g_aγγ, by over an order of magnitude. We exploit the Bormann effect of Laue crystals in a light-shining-through-wall experiment, with broad sensitivity to m_a ≲ 22 eV. Moreover for m_a∈ (3460, 3480) eV our sensitivity reaches down to the QCD axion coupling prediction, providing the most stringent laboratory constraints in this mass range.

High-resolution fluid simulations for plasma physics and astrophysics rely on Particle in cell (PIC) and hydrodynamic solvers (e.g., FLASH) to resolve shock dominated, multiscale phenomena, but their high computational cost severely limits scalability. This motivates the development of learning based surrogate models, which offer a promising route to accelerate these simulations while preserving physical fidelity. In this work, we study the Sedov Taylor shock propagation problem using a physics informed graph based surrogate model, Physics Informed MeshGraphNet (PhyMGN), designed for grid-based hydrodynamics. By incorporating weak physics constraints derived from the Euler equations using finite difference method, the model captures the self similar shock evolution and associated flow structures without explicitly solving the full hydrodynamic equations at each timestep. Comparing to the baseline MeshGraphNet model, PhyMGN is able to generalize beyond the training regime with a higher accuracy and preserves differentiability in parameter space while achieving a substantial reduction in computational cost relative to conventional numerical solvers.

We extend the beat-wave framework for laser harmonic generation - where spectra form regular lattices in Fourier space - to the nonlinear response of isotropic chiral media driven by locally chiral light. We represent the enantio-sensitive response of the medium by a chiral zero-frequency (DC) mode derived from the transverse spin density induced by structured or focused fields. Beating between this DC mode and the driving electromagnetic modes yields alternating chiral and achiral contributions on a regular harmonic lattice. We derive a general criterion for when chiral and achiral pathways overlap at the same harmonic and generate enantio-sensitive interference that survives spatial or angular integration (global chirality), versus when enantio-sensitivity remains confined to spatially varying patterns (local chirality). We apply the criterion to published configurations of synthetic chiral light, including OAM-carrying bicircular fields and crossed multicolour beams, and show that it reproduces and clarifies their reported global-chirality and beam-bending regimes.

We discuss the emission of pairs of photons by charges with generic worldlines in the Minkowski vacuum from the viewpoint of inertial observers and interpret them from the perspective of Rindler observers. We show that the emission of pairs of Minkowski photons-commonly referred to as Unruh radiation-corresponds, in general, to three distinct processes according to Rindler observers: scattering, and emission and absorption of pairs of Rindler photons. In the special case of uniformly accelerated charges, the radiation observed in the inertial frame can be fully described by the scattering channel in the Rindler frame. Therefore, the emission of pairs of Minkowski photons can be seen as further evidence supporting the Unruh effect.

We discuss the emission of radiation from general sources in quantum scalar, electromagnetic, and gravitational fields using the Rindler coordinate frame, which is suitable for uniformly accelerated observers, in the Minkowski vacuum. In particular, we point out that to recover, from the point of view of uniformly accelerated observers in the interaction picture, the usual Larmor radiation, which is independent of the choice of the vacuum state, it is necessary to incorporate the Unruh effect assuming the Minkowski vacuum. Thus, the observation of classical Larmor radiation in the Minkowski vacuum could be seen as vindicating the Unruh effect in the sense that it is not correctly recovered in the uniformly accelerated frame unless the Unruh effect is taken into account.

The generation of dense electron-positron pair beams in the laboratory can enable direct tests of theoretical models of γ-ray bursts and active galactic nuclei. We have successfully achieved this using ultrarelativistic protons accelerated by the Super Proton Synchrotron at CERN. In the first application of this experimental platform, the stability of the pair beam is studied as it propagates through a meter-length plasma, analogous to TeV γ-ray-induced pair cascades in the intergalactic medium. It has been argued that pair beam instabilities disrupt the cascade, thus accounting for the observed lack of reprocessed GeV emission from TeV blazars. If true, this would remove the need for a moderate strength intergalactic magnetic field to explain the observations. We find that the pair beam instability is suppressed if the beam is not perfectly collimated or monochromatic, hence the lower limit to the intergalactic magnetic field inferred from γ-ray observations of blazars is robust.

A high-power laser pulse at normal incidence onto a plane solid target will generate odd harmonics of its frequency. However, the spacing of the harmonic lines in this configuration is fixed. Here, we study harmonic generation using two laser beams incident on a plane target at small, opposite angles to the target normal, via particle-in-cell simulations. When looking at the harmonic radiation in a specific direction via a narrow slit or pinhole, we select an angle-dependent subset of the harmonic spectrum. This way, we obtain a harmonic frequency comb that we control via the observation angle and the input laser frequency. The divergence of the harmonic radiation will be reduced by using wider laser spots, thus increasing the efficacy of the scheme. We will discuss extensions to this scheme, such as using beams with unequal frequencies, a slight tilt of the target, or employing more than two beams.

Resonant axion-plasmon conversion in the magnetospheres of magnetars may substantially impact the landscape of dark-matter axion detection. This Letter explores how resonant axion-plasmon conversion modify the expected radio signals from axion-photon conversions observed on Earth. Critically, the resonant conversion radius lies within the region expected for axion-photon conversion, introducing a nonradiative power loss that diminishes the anticipated photon flux. Our analysis demonstrates that this effect can reduce radio telescope sensitivities, shifting them into regions excluded by previous experiments. These findings compel a reassessment of experimental constraints derived from radio signatures of axion-photon conversions and highlight the necessity of accounting for plasmon effects in astrophysical axion searches. The presented corrections provide critical insights for refining the detection strategies of future telescope-based dark matter axion experiments.

We present new exclusion bounds obtained at the European X-ray Free Electron Laser facility (EuXFEL) on axion-like particles (ALPs) in the mass range 10^-3 eV < m_a < 10^4 eV. Our experiment exploits the Primakoff effect via which photons can, in the presence of a strong external electric field, decay into axions, which then convert back into photons after passing through an opaque wall. While similar searches have been performed previously at a 3^rd generation synchrotron, our work demonstrates improved sensitivity, exploiting the higher brightness of X-rays at EuXFEL.

In our earlier work on harmonic generation with complex light [Nature Communications 15, 6878 (2024)], we demonstrated how the harmonic spectrum of a complex laser beam in a nonlinear medium can be obtained through the judicious application of the ``beat wave'' concept. In this paper, we show how the same results can be obtained via the full set of symmetries of the initial laser-target configuration, and how this can be reconciled with the ``beat wave'' approach. We also highlight the connections between our work and existing theory for diffraction of EM waves from crystals: Laue equations, Mathieu equation, and theorems by Noether, Floquet and Bloch. The specific nature of our approach to harmonic spectra allows these connections to be revealed. We illustrate this with numerous examples taken from existing literature to show the wide applicability of our approach.

The axion is a hypothetical particle associated with a possible solution to the strong CP problem and is a leading candidate for dark matter. In this Letter we investigate the emission of axions by accelerated electrons. We find the emission probability and energy within the WKB approximation for an electron accelerated by an electromagnetic field. As an application, we estimate the number of axions produced by electrons accelerated using two counterpropagating high-intensity lasers and discuss how they would be converted to photons to be detected. We find that, under realistic experimental conditions, competitive model-independent bounds on the coupling between the axion and the electron could be achieved in such an experiment.

The non-linear optical process of laser harmonic generation (HG) enables the creation of high quality pulses of UV or even X-ray radiation, which have many potential uses at the frontiers of experimental science, ranging from lensless microscopy to ultrafast metrology and chiral science. Although many of the promising applications are enabled by generating harmonic modes with orbital angular momentum (OAM), independent control of the harmonic frequency and OAM level remains elusive. Here we show, through a theoretical approach, validated with 3D simulations, how unique 2-D harmonic progressions can be obtained, with both frequency and OAM level tuned independently, from tailored structured targets in both reflective and transmissive configurations. Through preferential selection of a subset of harmonic modes with a specific OAM value, a controlled frequency comb of circularly polarised harmonics can be produced. Our approach to describe HG, which simplifies both the theoretical predictions and the analysis of the harmonic spectrum, is directly applicable across the full range of HG mechanisms and can be readily applied to investigations of OAM harmonics in other processes, such as OAM cascades in Raman amplification, or the analysis of harmonic progressions in nonlinear optics.

On the path towards high-gain inertial confinement fusion ignition, foams are being considered to tamp the hohlraum wall-motion, and mitigate laser backscattering from Stimulated Raman Scattering (SRS) and Stimulated Brillouin Scattering (SBS). Here we present the results from an experimental campaign on foam-filled hohlraums conducted at the OMEGA laser facility. SiO2 foam-fills, with densities as low as 1 mg/cm3, successfully reduce the gold wall expansion, with laser backscattering comparable to gas-fill.

We present numerical simulations used to interpret laser-driven plasma experiments at the GSI Helmholtz Centre for Heavy Ion Research. The mechanisms by which non-thermal particles are accelerated in astrophysical environments, e.g., the solar wind, supernova remnants, and gamma ray bursts, is a topic of intense study. When shocks are present, the primary acceleration mechanism is believed to be first-order Fermi, which accelerates particles as they cross a shock. Second-order Fermi acceleration can also contribute, utilizing magnetic mirrors for particle energization. Despite this mechanism being less efficient, the ubiquity of magnetized turbulence in the universe necessitates its consideration. Another acceleration mechanism is the lower-hybrid drift instability, arising from gradients of both density and magnetic field, which produce lower-hybrid waves with an electric field that energizes particles as they cross these waves. With the combination of high-powered laser systems and particle accelerators, it is possible to study the mechanisms behind cosmic-ray acceleration in the laboratory. In this work, we combine experimental results and high-fidelity three-dimensional simulations to estimate the efficiency of ion acceleration in a weakly magnetized interaction region. We validate the FLASH magneto-hydrodynamic code with experimental results and use OSIRIS particle-in-cell code to verify the initial formation of the interaction region, showing good agreement between codes and experimental results. We find that the plasma conditions in the experiment are conducive to the lower-hybrid drift instability, yielding an increase in energy Delta E of similar to 264 keV for 242 MeV calcium ions.

There is renewed interest in direct-drive inertial confinement fusion, following the milestone December 2022 3.15 MJ ignition result on the National Ignition Facility. A key obstacle is the control of the two-plasmon decay instability. Here, recent advances in inhomogeneous turbulence theory are applied to the broadband parametric instability problem for the first time. A novel dispersion relation is derived for the two-plasmon decay in a uniform plasma valid under broad-bandwidth laser fields with arbitrary power spectra. The effects of temporal incoherence on the instability are then studied. In the limit of large bandwidth, the well-known scaling relations for the growth rate are recovered, but it is shown that the result is more sensitive to the spectral shape of the laser pulse rather than to its coherence time. The range of wavenumbers of the excited plasma waves is shown to be substantially broadened, suggesting that the absolute instability is favoured in regions further away from the quarter critical density. The intermediate-bandwidth regime is explored numerically - the growth rate is reduced to half its monochromatic value for laser intensities of $10<^>{15} \ \text {W}\ \text {cm}<^>{-2}$ and relatively modest bandwidths of $5 \ \text {THz}$. The instability-quenching properties of a spectrum of discrete lines spread over some bandwidth have also been studied. The reduction in the growth rate is found to be somewhat lower compared with the continuous case but is still significant, despite the fact that, formally, the coherence time of such a laser pulse is infinite.