Publications

2026

• 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>
• 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