CyberSec Research

The May 2024 geomagnetic superstorm provided the opportunity to explore how strong wave-particle interactions affect energetic electron precipitation under intense driving. Using coordinated measurements from a balloon-borne Timepix-based X-ray detector, ground-based riometers and magnetometers, and Arase satellite observations, we identified quasi-periodic bursts of energetic electron precipitation coincident with Pc5 ultra low frequency (ULF) wave oscillations. Arase satellite data revealed energy-dispersed trapped energetic electron flux modulations in the 'seed' energy range, indicating that trapped electron flux was likely modulated by ULF waves. This letter reveals that these flux enhancements surpassed the Kennel-Petschek (K-P) limit, creating intense chorus waves and driving periodic electron precipitation. Drift-dispersion analysis traced these modulations back to a source in the post-noon magnetospheric sector, matching balloon and ground-based measurements. Here, we propose a novel indirect ULF wave-driven mechanism for modulated energetic electron precipitation, whereby periodic modulations of `seed' electron fluxes enhance electron losses.

High-energy heavy-ion particle accelerators have long served as a proxy for the harsh space radiation environment, enabling both fundamental life-science research and applied testing of flight components. Typically, monoenergetic high-energy heavy-ion beams are used to mimic the complex mixed radiation field encountered in low Earth orbit and beyond. However, synergistic effects arising from the spatial or temporal proximity of interactions of different radiation qualities in a mixed field cannot be fully assessed with such beams. Therefore, spearheaded by developments at the NASA Space Radiation Laboratory, the GSI Helmholtzzentrum fuer Schwerionenforschung, supported by ESA, has developed advanced space radiation simulation capabilities to support space radiation studies in Europe. Here, we report the design, optimization, and in-silico benchmarking of GSI's hybrid active-passive GCR simulator. Additionally, a computationally optimized phase-space particle source for Geant4 is presented, which will be made available to external users to support their own in-silico studies and experimental planning.

Space radiation is one of the major obstacles to space exploration. If not mitigated, radiation can interact both with biological and electronic systems, inducing damage and posing significant risk to space missions. Countermeasures can only be studied effectively with ground-based accelerators that act as a proxy for space radiation, typically with a harsher radiation field that worsen the effects of space radiation. Following an in-silico design and optimization process we have developed a galactic cosmic ray (GCR) simulator using a hybrid active-passive methodology. In this approach, the primary beam energy is actively switched and the beam interacts with specifically designed passive modulators. In this paper, we present the implementation of such a GCR simulator and its experimental microdosimetric characterization. Measuring the GCR field is of paramount importance, both before providing it to the user as a validated radiation field and for achieving the best possible radiation description. The issue is addressed in this paper by using a tissue-equivalent proportional counter to measure radiation quality and by comparing experimental measurements with Monte Carlo simulations. In conclusion, we will demonstrate the GCR simulator's capability to reproduce a GCR field.

The CUbesat Solar Polarimeter (CUSP) project is a CubeSat mission planned for a launch in low-Earth orbit and aimed to measure the linear polarization of solar flares in the hard X-ray band by means of a Compton scattering polarimeter. CUSP will allow us to study the magnetic reconnection and particle acceleration in the flaring magnetic structures of our star. CUSP is a project in the framework of the Alcor Program of the Italian Space Agency aimed at developing new CubeSat missions. It is undergoing a 12-month Phase B that started in December 2024. The Compton polarimeter on board CUSP is composed of two acquisition chains based on plastic scintillators read out by Multi-Anode PhotoMultiplier Tubes for the scatterer part and GAGG crystals coupled to Avalanche PhotoDiodes for the absorbers. An event coincident between the two readout schemes will lead to a measurement of the incoming X-ray's azimuthal scattering angle, linked to the polarization of the solar flare in a statistical manner. The current status of the CUSP mission design, mission analysis, and payload scientific performance will be reported. The latter will be discussed based on preliminary laboratory results obtained in parallel with Geant4 simulations.

Magnetic reconnection is a fundamental and omnipresent energy conversion process in plasma physics. Novel observations of fields and particles from Parker Solar Probe (PSP) have shown the absence of reconnection in a large number of current sheets in the near-Sun solar wind. Using near-Sun observations from PSP Encounters 4 to 11 (Jan 2020 to March 2022), we investigate whether reconnection onset might be suppressed by velocity shear. We compare estimates of the tearing mode growth rate in the presence of shear flow for time periods identified as containing reconnecting current sheets versus non-reconnecting times, finding systematically larger growth rates for reconnection periods. Upon examination of the parameters associated with reconnection onset, we find that 85% of the reconnection events are embedded in slow, non-Alfvenic wind streams. We compare with fast, slow non-Alfvenic, and slow Alfvenic streams, finding that the growth rate is suppressed in highly Alfvenic fast and slow wind and reconnection is not seen in these wind types, as would be expected from our theoretical expressions. These wind streams have strong Alfvenic} flow shear, consistent with the idea of reconnection suppression by such flows. This could help explain the frequent absence of reconnection events in the highly Alfvenic, near-Sun solar wind observed by PSP. Finally, we find a steepening of both the trace and magnitude magnetic field spectra within reconnection periods in comparison to ambient wind. We tie this to the dynamics of relatively balanced turbulence within these reconnection periods and the potential generation of compressible fluctuations.

Near-Earth asteroid 2024 YR4 was discovered on 2024-12-27 and its probability of Earth impact in December 2032 peaked at about 3% on 2025-02-18. Additional observations ruled out Earth impact by 2025-02-23. However, the probability of lunar impact in December 2032 then rose, reaching about 4% by the end of the apparition in May 2025. James Webb Space Telescope (JWST) observations on 2025-03-26 estimated the asteroid's diameter at 60 +/- 7 m. Studies of 2024 YR4's potential lunar impact effects suggest lunar ejecta could increase micrometeoroid debris flux in low Earth orbit up to 1000 times above background levels over just a few days, possibly threatening astronauts and spacecraft. In this work, we present options for space missions to 2024 YR4 that could be utilized if lunar impact is confirmed. We cover flyby & rendezvous reconnaissance, deflection, and robust disruption of the asteroid. We examine both rapid-response and delayed launch options through 2032. We evaluate chemical and solar electric propulsion, various launch vehicles, optimized deep space maneuvers, and gravity assists. Re-tasking extant spacecraft and using built spacecraft not yet launched are also considered. The best reconnaissance mission options launch in late 2028, leaving only approximately three years for development at the time of this writing in August 2025. Deflection missions were assessed and appear impractical. However, kinetic robust disruption missions are available with launches between April 2030 and April 2032. Nuclear robust disruption missions are also available with launches between late 2029 and late 2031. Finally, even if lunar impact is ruled out there is significant potential utility in deploying a reconnaissance mission to characterize the asteroid.

Large-amplitude electrostatic fluctuations are routinely observed by spacecraft upon traversal of collisionless shocks in the heliosphere. Kinetic simulations of shocks have struggled to reproduce the amplitude of such fluctuations, complicating efforts to understand their influence on energy dissipation and shock structure. In this paper, 1D particle-in-cell simulations with realistic proton-to-electron mass ratio are used to show that in cases with upstream electron temperature $T_e$ exceeding the ion temperature $T_i$, the magnitude of the fluctuations increases with the electron plasma-to-cyclotron frequency ratio $\omega_{pe}/\Omega_{ce}$, reaching realistic values at $\omega_{pe}/\Omega_{ce} \gtrsim 30$. The large-amplitude fluctuations in the simulations are shown to be associated with electrostatic solitary structures, such as ion phase-space holes. In the cases where upstream temperature ratio is reversed, the magnitude of the fluctuations remains small.

We investigate energetic particle diffusion in the inner heliosphere (approximately 0.06-0.3 AU) explored by Parker Solar Probe (PSP). Parallel (kappa_parallel) and perpendicular (kappa_perp) diffusion coefficients are calculated using second-order quasi-linear theory (SOQLT) and unified nonlinear transport (UNLT) theory, respectively. PSP's in situ measurements of magnetic turbulence spectra, including sub-Alfvenic solar wind, are decomposed into parallel and perpendicular wavenumber spectra via a composite two-component turbulence model. These spectra are then used to compute kappa_parallel and kappa_perp across energies ranging from sub-GeV to GeV. Our results reveal a strong energy and radial distance dependence in kappa_parallel. While kappa_perp remains much smaller, it can increase in regions with relatively high turbulence levels delta B / B0. To validate our results, we estimate kappa_parallel using the upstream time-intensity profile of a solar energetic particle event observed by PSP and compare it with theoretical values from different diffusion models. Our results suggest that the SOQLT-calculated parallel diffusion generally shows better agreement with SEP intensity-derived estimates than the classic QLT model. This indicates that the SOQLT framework, which incorporates resonance broadening and nonlinear corrections and does not require an ad hoc pitch-angle cutoff, may provide a more physically motivated description of energetic particle diffusion near the Sun.

Optically trapped Silica nanoparticles are a promising tool for precise sensing of gravitational or inertial forces and fundamental physics, including tests of quantum mechanics at 'large' mass scales. This field, called levitated optomechanics can greatly benefit from an application in weightlessness. In this paper we demonstrate the feasibility of such setups in a microgravity environment for the first time. Our experiment is operated in the GraviTower Bremen that provides up to 2.5 s of free fall. System performance and first release-recapture experiments, where the particle is no longer trapped are conducted in microgravity. This demonstration should also be seen in the wider context of preparing space missions on the topic of levitated optomechanics.

Two fluid simulations using local Landau-fluid closures derived from linear theory provide an efficient computational framework for plasma modelling, since they bridge the gap between computationally intensive kinetic simulations and fluid descriptions. Their accuracy in representing kinetic effects depends critically on the validity of the linear approximation used in the derivation: the plasma should not be too far from local thermodynamic equilibrium, LTE. However, many of the problems where these models are of particular interest (such as plasma turbulence and instabilities) are in fact quite far from LTE. The question then arises, if kinetic scale processes are still sufficiently well captured outside of the theoretical regime of applicability of the closure. In this paper, we show that two fluid simulations with Landau fluid closures can effectively reproduce the energy spectra obtained with fully kinetic Vlasov simulations, used as references, as long as the local closure parameter is appropriately chosen. Our findings validate the usage of two fluid simulations with Landau-fluid closure as a possible alternative to fully kinetic simulations of turbulence, in cases where being able to simulate extremely large domains is of particular interest.

Magnetic reconnection is an explosive process that accelerates particles to high energies in Earth's magnetosphere, offering a unique natural laboratory to study this phenomenon. We performed fully kinetic 2D simulations of a reconnection event observed by the Magnetospheric Multiscale mission and compared the resulting ion and electron energy distributions with observations. The simulations capture the overall shape and evolution of non-thermal energy distributions for both species, but generally underestimate the very-high-energy tail of the electron spectrum. Variations in numerical parameters have negligible effects on the resulting spectra, while the initial upstream temperatures instead play a critical role in reproducing the observed distributions. This work presents a novel analysis of particle acceleration in fully kinetic modeling of reconnection directly informed by observed, realistic parameters; highlights the limitations of 2D simulations and underlines the need for more realistic simulations (e.g. employing 3D setups) to capture the observed particle energization more accurately.

This paper introduces a versatile approach for computing the risk of collision specifically tailored for scenarios featuring low relative encounter velocities, but with potential applicability across a wide range of situations. The technique employs Differential Algebra (DA) to express the non-linear dynamical flow of the initial distribution in the primary-secondary objects relative motion through high-order Taylor polynomials. The entire initial uncertainty set is subdivided into subsets through Automatic Domain Splitting (ADS) techniques to control the accuracy of the Taylor expansions. The methodology samples the initial conditions of the relative state and evaluates the polynomial expansions for each sample while retaining their temporal dependency. The classical numerical integration of the initial statistics over the set of conditions for which a collision occurs is thus reduced to an evaluation of mono-dimensional time polynomials. Specifically, samples reaching a relative distance below a critical value are identified along with the time at which this occurs. The approach is tested against a Monte Carlo (MC) simulation for various literature test cases, yielding accurate results and a consistent gain in computational time.

The Integrating Miniature Piggyback for Impulsive Solar Hard X-rays (IMPISH) is a piggyback mission originally designed for the second flight of the Gamma-Ray Imager/Polarimeter for Solar flares (GRIPS-2) Antarctic balloon. IMPISH will take measurements of collimated, full-Sun X-ray spectra with the goal of detecting sub-second variations (order of tens of milliseconds) of nonthermal X-ray emission during the impulsive phase of large solar flares to probe particle acceleration mechanisms driven by magnetic reconnection. The IMPISH detector system, made up of four identical detectors totaling 64 cm$^2$ effective area, is capable of measuring from ~10 keV to over 200 keV through the use of silicon photomultipliers (SiPMs) and LYSO scintillators. At the stratospheric altitude of GRIPS-2, the effective lower energy limit is ~20-30 keV. The geometry of the LYSO crystal has been optimized to balance the light collection efficiency with the effective area required for stratospheric X-ray measurements. Development of the IMPISH detectors has introduced a path for a low-cost solution to fast solar X-ray measurements across a large energy range utilizing commercially available components. The payload has a 3U form factor and has been designed so that both the electronics and detectors may be easily adaptable for space-based missions.

The Martian magnetosheath acts as a conduit for mass and energy transfer between the upstream solar wind and its induced magnetosphere. However, our understanding of its global properties remains limited. Using nine years of data from NASA's Mars Atmosphere and Volatile EvolutioN (MAVEN) mission, we performed a quantitative statistical analysis to explore the spatial distribution of the magnetic fields, solar wind and planetary ions in the magnetosheath. We discovered significant asymmetries in the magnetic field, solar wind protons, and planetary ions between the quasi-perpendicular and quasi-parallel magnetosheaths. The asymmetries in the Martian magnetosheath exhibit both similarities and differences compared to those in the Earth's and Venus' magnetosheaths. These results indicate that the Martian magnetosheath is distinctly shaped by both shock geometry and planetary ions.

Understanding how shocks interact with coronal structures is crucial for understanding the mechanisms of particle acceleration in the solar corona and inner heliosphere. Using simultaneous radio and white-light observations, we investigate the interaction between a CME-driven shock and a plasmoid. LASCO and STEREO-A COR-2 white-light images are analyzed to track the evolution of the plasmoid, CME and its associated shock, while the Wind/WAVES and STEREO/WAVES dynamic spectra provide complementary radio signatures of the shock-plasmoid interaction at $\approx$7 R$_\odot$. An interplanetary Type II radio burst was detected as the shock propagated through the plasmoid. The merging of the plasmoid into the CME was accompanied by interplanetary Type III radio bursts, suggesting escaping electron beams during the reconnection process. These observations clearly demonstrate that shock-plasmoid interactions can enhance the efficiency of particle acceleration associated with CMEs, with implications for electron acceleration in flare and heliospheric current sheets as well.

Recent years have seen many arguments for cosmic rays (CRs) as an important influence on galactic and circumgalactic (CGM) physics, star and galaxy formation. We present a pedagogical overview of state-of-the-art modeling of CR-magnetohydrodynamics (CR-MHD) on macro scales (~kpc), highlighting their fundamental dependence on the micro (< au) scales of CR gyro orbits and meso (~pc) scales of CR mean-free-paths, intended to connect the extragalactic, Galactic, and plasma CR transport modeling communities. We note the pitfalls and systematic errors that arise from older assumptions in CR modeling, including: use of a simple Fokker-Planck equation or ad-hoc two-moment formalisms for transport; assumption of leaky boxes or plane-parallel or shear-periodic boundaries for comparison to local interstellar medium (LISM) observations; ignoring detailed LISM constraints on CR spectra (e.g. focusing only on extragalactic observables or spectrally integrated models); assuming CR transport is mediated by classical models of advection, streaming from self-confinement (super-Alfvenic or Alfvenic), or extrinsic turbulence. We emphasize recent progress addressing these: development of rigorously-derived CR-MHD equations; use of global, 3D galaxy+halo models for LISM comparisons; new methods for full-spectrum dynamics; novel models for intermittent scattering and/or new drivers. We compile extragalactic+LISM observations to show how ~GeV CR transport is being rapidly constrained in the CGM, and present phenomenological models which can be used in future simulations. We conclude by highlighting critical open questions for micro, meso, and macro-scale CR-MHD simulations.

The Alfv\'en surface -- where the solar wind exceeds the local Alfv\'en speed as it expands into interplanetary space -- is now routinely probed by NASA's Parker Solar Probe (PSP) in the near-Sun environment. The size of the Alfv\'en surface governs how efficiently the solar wind braking torque causes the Sun to spin-down. We aimed to characterise the size and evolution of the Alfv\'en surface as magnetic activity increased during solar cycle 25. The Alfv\'en surface was extrapolated from the solar wind mass and magnetic flux measured by the SWEAP and FIELDS instrument suites onboard PSP. We accounted for the acceleration of the solar wind along Parker spiral magnetic field lines and used potential field source surface modelling to determine the sources of the solar wind. The longitudinally averaged Alfv\'en radius measured by PSP grew from 11 to 16 solar radii as solar activity increased. Accordingly, the solar wind angular momentum-loss rate grew from $\sim$1.4$\times 10^{30}$ erg to 3$\times 10^{30}$ erg. Both the radial and longitudinal scans of the solar wind contained fluctuations of 10-40\% from the average Alfv\'en radius in each encounter. Structure in the solar corona influenced the morphology of the Alfv\'en surface, which was smallest around the heliospheric current sheet and pseudostreamers. The Alfv\'en surface was highly structured and time-varying however, at large-scales, organised by the coronal magnetic field. The evolution of the solar corona over the solar cycle systematically shifted the magnetic connectivity of PSP and influenced our perception of the Alfv\'en surface. The Alfv\'en surface was 30\% larger than both thermally-driven and Alfv\'en wave-driven wind simulations with the same mass-loss rate and open magnetic flux, but had a similar dependence on the wind magnetisation parameter.

The Wang-Sheeley-Arge (WSA) model has been the cornerstone of operational solar wind forecasting for nearly two decades, owing to its simplicity and physics-based formalism. However, its performance is strongly dependent on several empirical parameters that are typically fixed or tuned manually, limiting its adaptability across varying solar conditions. In this study, we present a neural enhancement to the WSA framework (referred to as WSA+) that systematically optimizes the empirical parameters of the WSA solar wind speed relation using in-situ observations within a differentiable physics-constrained pipeline. The approach operates in two stages: first, a neural optimizer adjusts WSA parameters independently for each Carrington Rotation to better match the observed solar wind data. Then, a neural network learns to predict these optimized speed maps directly from magnetogram-derived features. This enables generalization of the optimization process and allows inference for new solar conditions without manual tuning. The developed WSA+ preserves the interpretability of the original relation while significantly improving the match with OMNI in-situ data across multiple performance metrics, including correlation and error statistics. It consistently outperforms the traditional WSA relation across both low and high solar activity periods, with average improvements of approximately 40 percent. By integrating data-driven learning with physical constraints, WSA+ offers a robust and adaptable enhancement, with immediate utility as a drop-in replacement in global heliospheric modeling pipelines.