Publications

2026

• Formation of gradients of atomic oxygen in nanosecond plasma for plasma-assisted detonation: experimental and numerical study
  
  • Lafaurie Victor
  • Shu Zhan
  • Zhang B
  • Terentjeviene M
  • Billeau Jean-Baptiste
  • Orel Inna
  • Hoyos Aristizabal Samuel
  • Vidal Pierre
  • Starikovskaia Svetlana
  Plasma Sources Science and Technology, IOP Publishing, 2026, 35 (2), pp.025022. This work aims at producing a gradient of atomic oxygen on a scale of 10 cm in a plane-to-plane nanosecond discharge in 150 mbar of air with a varying gap size for applications in combustion and ignition of detonation waves. Local measurements of atomic oxygen density along the discharge span, at varying heights between high-voltage and grounded electrode, are performed with Xe calibrated O-TALIF and validated by 2D numerical modelling. They both show existence of a gradient of atomic density of oxygen along the span. Reduced electric field is measured with two experimental techniques: optical emission spectroscopy by a spectral band intensity ratio of the first negative system and the second positive system of nitrogen, and E-FISH. It is also compared with numerical modelling. All techniques show existence of a gradient of reduced electric field along the span. This distribution of reduced electric field, in combination with the non-uniform energy deposition in the plasma, is shown to explain the measured gradient of density of atomic oxygen. (10.1088/1361-6595/ae3f53)
  DOI : 10.1088/1361-6595/ae3f53
• Benchmark for two-dimensional large scale coherent structures in partially magnetized E × B plasmas—community collaboration & lessons learned
  
  • Powis Andrew T
  • Ahedo Eduardo
  • Álvarez Laguna Alejandro
  • Barléon Nicolas
  • Bello-Benítez Enrique
  • Beving Lucas
  • Boeuf Jean-Pierre
  • Bogopolsky Guillaume
  • Bourdon Anne
  • Cichocki Filippo
  • Cuenot Bénédicte
  • Denig Andrew
  • Donkó Zoltán
  • Elias Paul-Quentin
  • Encinar Miguel
  • Eremin Denis
  • Fajardo Pablo
  • Faraji Farbod
  • Fubiani Gwenael
  • Garrigues Laurent
  • Hara Kentaro
  • Hartmann Peter
  • Hopkins Matthew
  • Kaganovich Igor D
  • Knoll Aaron
  • Lapenta Giovanni
  • Magin Thierry
  • Marín-Cebrián Alberto
  • Merino Mario
  • Minelli Pierpaolo
  • Papahn Zadeh Mina
  • Parodi Pietro
  • Petronio Federico
  • Reza Maryam
  • Smolyakov Andrei I
  • Sydorenko Dmytro
  • Taccogna Francesco
  • Turner Miles M
  • Vermorel Olivier
  • Villafana Willca
  • Xu Liang
  Plasma Sources Science and Technology, IOP Publishing, 2026, 35 (2), pp.025002. Abstract Low-temperature plasmas (LTPs) are essential to both fundamental scientific research and critical industrial applications. As in many areas of science, numerical simulations have become a vital tool for uncovering new physical phenomena and guiding technological development. Code benchmarking remains crucial for verifying implementations and evaluating performance. This work continues the Landmark benchmark initiative, a series specifically designed to support the verification of LTP codes. In this study, seventeen simulation codes from a collaborative community of nineteen international institutions modeled a partially magnetized E × B Penning discharge. The emergence of large scale coherent structures, or rotating plasma spokes, endows this configuration with an enormous range of time scales, making it particularly challenging to simulate. The codes showed excellent agreement on the rotation frequency of the spoke as well as key plasma properties, including time-averaged ion density, plasma potential, and electron temperature profiles. Achieving this level of agreement came with challenges, and we share lessons learned on how to conduct future benchmarking campaigns. Comparing code implementations, computational hardware, and simulation runtimes also revealed interesting trends, which are summarized with the aim of guiding future plasma simulation software development. (10.1088/1361-6595/ae3985)
  DOI : 10.1088/1361-6595/ae3985
• Evidence of an Extended Alfvén Wing System at Enceladus: Cassini's Multi‐Instrument Observations
  
  • Hadid Lina Z.
  • Chust Thomas
  • Wahlund Jan-Erik
  • Morooka Michiko W
  • Roussos Elias
  • Witasse Olivier
  • Rabia Jonas
  • Pisa David
  • Kim Konstantin
  • Edberg Niklas J T
  • Rymer Abigail M
  • Lamy Laurent
  • Kotsiaros Stavros
  • Aizawa Sae
  • Jeandet Alexis
  • Modolo Ronan
  • André Nicolas
  • Canu Patrick
  • Bowers Charles F
  • Jia Xianzhe
  • Coates Andrew J
  • Jones Geraint H
  • Parsec‐wallis Anna
  • Agiwal Omakshi
  • Holmberg Mika K G
  • Nénon Quentin
  • Cao Hao
  • Kurth William S
  • Dougherty Michele K
  Journal of Geophysical Research Space Physics, American Geophysical Union/Wiley, 2026, 131 (2). We report in situ evidence for Enceladus' Alfvén wing system and its coupling with Saturn's ionosphere, based on multi-instrument observations from the Cassini spacecraft. Analysis of 36 events, including 13 from non-flyby paths, confirms the existence of a Main Alfvén Wing (MAW) current system generated at Enceladus, and associated Reflected Alfvén Wings (RAWs) occurring both at Saturn's ionosphere and on the density gradient of Enceladus' plasma torus, extending longitudinally to at least ∼ 120°(∼2,000 moon radii) downstream of the moon. Additionally, the observations reveal the systematic existence of a filamentation process of these large-scale Alfvénic perturbations (MAW and RAWs) during their propagation at any distance from their source. These findings demonstrate a more extensive electrodynamic coupling than previously reported for Enceladus and more generally for any moon-magnetosphere interaction. Moreover, the observation of energetic electron depletions and water-group ion signatures at longitudes even further from the moon supports the interpretation of an extended and persistent interaction region. These results highlight Enceladus' role in shaping Saturn's magnetospheric environment and underscore the importance of future missions to exhaustively analyze this type of complex interaction between a moon and a planet. Plain Language Summary Saturn's small icy moon Enceladus interacts with the planet's magnetic field, generating intermittent aurora in Saturn's upper atmosphere and electromagnetic waves that travel along invisible magnetic connections between them. During its 13-year mission, the Cassini spacecraft repeatedly crossed these magnetic field lines linked to Enceladus. We used data from several Cassini instruments to study how energy and particles move between the moon and Saturn. We detected wave activity characteristic of Alfvén waves (similar to vibrations on a string), forming as Saturn's magnetic field flows past Enceladus. Due to a complex system of reflection at both Saturn's ionosphere and the boundary of Enceladus' torus, these waves were found not only near the moon but also trailing far behind it, extending more than 504,000 km (over 2,000 times the moon's radius) behind it. This is the first time that Alfvén waves have been observed to be directly linked to the charged particles associated with Enceladus. This shows that Enceladus plays a much bigger role in shaping Saturn's space environment than previously thought, and reveals how moons can influence their host planet across vast distances. (10.1029/2025ja034657)
  DOI : 10.1029/2025ja034657
• Magnetic reversals in a geodynamo model with a stably–stratified layer
  
  • Müller Nicolás P
  • Gissinger Christophe
  • Pétrélis François
  Physics of the Earth and Planetary Interiors, Elsevier, 2026, 371, pp.107502. We study the process of magnetic reversals in the presence of a stably-stratified layer below the core-mantle boundary using direct numerical simulations of the incompressible magnetohydrodynamics equations under the Boussinesq approximation in a spherical shell. We show that the dipolar-multipolar transition shifts to larger Rayleigh numbers in the presence of a stably-stratified layer, and that the dipolar strength of the magnetic field at the core-mantle boundary increases due to the skin effect. By imposing an heterogeneous heat flux at the outer boundary, we break the equatorial symmetry of the flow, and show that different heat flux patterns can trigger different dynamo solutions, such as hemispheric dynamos and polarity reversals. Using kinematic dynamo simulations, we show that the stably-stratified layer leads to similar growth rates of the dipole and quadrupole components of the magnetic field, playing the role of a conducting boundary layer, favouring magnetic reversals, and a dynamics predicted by low-dimensional models. (10.1016/j.pepi.2026.107502)
  DOI : 10.1016/j.pepi.2026.107502
• Benchmark for two-dimensional large scale coherent structures in partially magnetized E×B plasmas -Community collaboration & lessons learned
  
  • Powis Andrew T
  • Ahedo Eduardo
  • Laguna Alejandro Álvarez
  • Barléon Nicolas
  • Bello-Benítez Enrique
  • Beving Lucas
  • Boeuf Jean-Pierre
  • Bogopolsky Guillaume
  • Bourdon Anne
  • Cichocki Filippo
  • Cuenot Bénédicte
  • Denig Andrew
  • Donkó Zoltán
  • Elias Paul-Quentin
  • Encinar Miguel P
  • Eremin Denis
  • Fajardo Pablo
  • Faraji Farbod
  • Fubiani Gwenael
  • Garrigues Laurent
  • Hara Kentaro
  • Hartmann Peter
  • Hopkins Matthew
  • Kaganovich Igor D
  • Knoll Aaron
  • Lapenta Giovanni
  • Magin Thierry E
  • Marín-Cebrián Alberto
  • Merino Mario
  • Minelli Pierpaolo
  • Papahn Zadeh Mina
  • Parodi Pietro
  • Petronio Federico
  • Reza Maryam
  • Smolyakov Andrei I
  • Sydorenko Dmytro
  • Taccogna Francesco
  • Turner Miles M
  • Vermorel Olivier
  • Villafana Willca
  • Xu Liang
  , 2025. <div><p>Low-temperature plasmas are essential to both fundamental scientific research and critical industrial applications. As in many areas of science, numerical simulations have become a vital tool for uncovering new physical phenomena and guiding technological development. Code benchmarking remains crucial for verifying implementations and evaluating performance. This work continues the Landmark benchmark initiative, a series specifically designed to support the verification of low-temperature plasma codes. In this study, seventeen simulation codes from a collaborative community of nineteen international institutions modeled a partially magnetized E×B Penning discharge. The emergence of large scale coherent structures, or rotating plasma spokes, endows this configuration with an enormous range of time scales, making it particularly challenging to simulate. The codes showed excellent agreement on the rotation frequency of the spoke as well as key plasma properties, including time-averaged ion density, plasma potential, and electron temperature profiles. Achieving this level of agreement came with challenges, and we share lessons learned on how to conduct future benchmarking campaigns. Comparing code implementations, computational hardware, and simulation runtimes also revealed interesting trends, which are summarized with the aim of guiding future plasma simulation software development.</p></div>
• Hall thruster modeling with multiple simulation techniques: Model benchmarking, fluid–kinetic consistency, and experimental validation
  
  • Petronio Federico
  • Alvarez Laguna Alejandro
  • Bourdon Anne
  • Lafleur Trevor
  • Chabert Pascal
  Journal of Applied Physics, American Institute of Physics, 2026, 139 (9), pp.093303. Numerical plasma models are critical tools for aiding the design and understanding of electric propulsion systems, such as Hall thrusters, particularly when considering challenges associated with diagnostic access and reliable internal measurements. For complex plasma systems, such as Hall thrusters, theoretical verification solutions are often missing, and therefore, benchmarking represents an important element in assessing the correctness and consistency of the underlying mathematical model, and the computational performance of the numerical implementation. In this work, we benchmark three different numerical codes by simulating an SPT-100 Hall thruster under identical operating conditions. The codes include one-dimensional stationary and non-stationary fluid models describing the axial thruster direction, as well as a two-dimensional axial–azimuthal Particle-In-Cell/Monte Carlo Collision (PIC/MCC) model. A partial validation is performed with available experimental measurements of the discharge current, thrust, and anode specific impulse, showing good agreement. Overall, the fluid and PIC/MCC models compare favorably with each other, and several fluid approximations are found to be acceptable. For example, axial electron energy transport is relatively minor such that the electron temperature is reasonably determined by a local energy balance. Other approximations, however, require a more careful examination: particularly the assumption of Maxwellian electrons and the neglect of electron–wall collisions in the electron momentum balance equations. (10.1063/5.0305019)
  DOI : 10.1063/5.0305019
• Nonlinear phase synchronization and the role of spacing in shell models
  
  • Manfredini L.
  • Gürcan Ö D
  Physical Review E, American Physical Society (APS), 2026, 113 (1), pp.015101. A shell model can be considered as a self-similar chain of interacting triads, where each triad can be interpreted as a nonlinear oscillator that can be mapped to a spinning top. Investigating the relation between phase dynamics and intermittency in such a chain of nonlinear oscillators, it is found that synchronization is linked to increased energy transfer. In particular, our results indicate that the observed systematic increase of intermittency, as the shell spacing is decreased, is associated with strong phase alignment among consecutive triadic phases, facilitating the energy cascade. It is shown that while the overall level of synchronization can be quantied using a Kuramoto order parameter for the global phase coherence in the inertial range, a local, weighted Kuramoto parameter can be used for the detection of burst-like events propagating across shells in the inertial range. This novel analysis reveals how locally phase-locked states are associated with the passage of extreme events of energy ux. Applying this method to helical shell models ( i.e. for a class of helical interactions that couple the two helicities in a non separable topology) reveals that a reduction in phase coherence correlates with suppression of intermittency. When inverse cascade scenarios are considered using two dierent shell models including a non local helical shell model, and a local standard shell model with a modied conservation law, it was shown that a particular phase organization is needed in order to sustain the inverse energy cascade. It was also observed that the PDFs of the triadic phases were peaked in accordance with the basic considerations of the form of the ux, which suggests that a triadic phase of π/2 and -π/2 maximizes the forward and the inverse energy cascades respectively. (10.1103/2vxp-1k2t)
  DOI : 10.1103/2vxp-1k2t
• Properties of Magnetic Switchbacks in the Near-Sun Solar Wind
  
  • Badman Samuel T.
  • Fargette Naïs
  • Matteini Lorenzo
  • Agapitov Oleksiy V.
  • Akhavan-Tafti Mojtaba
  • Bale Stuart D.
  • Bharati Das Srijan
  • Bizien Nina
  • Bowen Trevor A.
  • Dudok de Wit Thierry
  • Froment Clara
  • Horbury Timothy
  • Huang Jia
  • Jagarlamudi Vamsee Krishna
  • Larosa Andrea
  • Madjarska Maria S.
  • Panasenco Olga
  • Pariat Etienne
  • Raouafi Nour E.
  • Rouillard Alexis P.
  • Ruffolo David
  • Sioulas Nikos
  • Soni Shirsh Lata
  • Sorriso-Valvo Luca
  • Suen Gabriel Ho Hin
  • Velli Marco
  • Verniero Jaye
  Space Science Reviews, Springer Verlag, 2026, 222. Magnetic switchbacks are fluctuations in the solar wind in which the interplanetary magnetic field sharply deflects away from its background direction so as to create folds in magnetic field lines while remaining of roughly constant magnitude. The magnetic field and velocity fluctuations are extremely well correlated in a way corresponding to Alfvénic fluctuations propagating away from the Sun. For a background field which is nearly radial this causes an outwardly propagating jet to form. Switchbacks and their characteristic velocity jets have recently been observed to be nearly ubiquitous by Parker Solar Probe with in situ measurements in the inner heliosphere within 0.3 AU. Their prevalence, substantial energy content, and potentially fundamental role in the dynamics of the outer corona and solar wind motivate the significant research efforts into their understanding. Here we review the in situ measurements of these structures (primarily by Parker Solar Probe). We discuss how they are identified and measured, and present an overview of the primary observational properties of these structures, both in terms of individual switchbacks and their collective arrangement into "patches". We identify both properties for which there is a strong consensus and those that have limited or qualified support and require further investigation. We identify and collate several open questions and recommendations for future studies. (10.1007/s11214-026-01267-w)
  DOI : 10.1007/s11214-026-01267-w