ERMES 20.0: Dual Applications in Seismology & Plasma

• ERMES 20.0 is a dual-use concept, representing a muon-beam seismic monitoring system in geophysics and a finite element solver in plasma physics.
• In seismology, it employs high-energy muon beams to detect stress-induced piezoelectric fields in quartz-rich rocks near faults for early earthquake detection.
• In plasma physics, it validates cold-plasma Maxwell equations and O–X mode conversion using both standard and regularized finite element formulations.

Searching arXiv for the cited papers to ground the article in current records. arXiv_search({"query":"id:(Serafini et al., 15 Jan 2026) OR id:(Otin et al., 30 Sep 2025)","max_results":5,"sort_by":"submittedDate","sort_order":"descending"}) arXiv_search returned 2 papers:

• (Serafini et al., 15 Jan 2026): "Active interrogation of underground piezoelectric fabrics using high energy muon beams propagating across seismogenic faults"
• (Otin et al., 30 Sep 2025): "Validation of ERMES 20.0 finite element code for MAST Upgrade O-X mode conversion" ERMES 20.0 is a term with two distinct technical uses in current arXiv literature. In geophysics, it denotes a conceptual mature version of ERMES, “Earthquake Reconnaissance using Muon beam Evolution in Silicon dioxide,” an active interrogation system in which a high-energy, collimated artificial muon beam traverses quartz-rich rocks across a seismogenic fault to monitor stress-dependent internal piezoelectric fields as a possible earthquake precursor (Serafini et al., 15 Jan 2026). In plasma physics and computational electromagnetics, ERMES 20.0 denotes an open-source, frequency-domain finite element code validated for cold-plasma Ordinary–Extraordinary mode conversion in an Electron Bernstein Wave benchmark representative of the MAST Upgrade EBW system (Otin et al., 30 Sep 2025).

1. Dual usage and nomenclature

The name “ERMES 20.0” does not identify a single research object. It is used for two technically unrelated systems: a proposed accelerator-based geophysical monitoring architecture, and a finite element plasma-wave solver.

Usage	Description	Domain
ERMES	“Earthquake Reconnaissance using Muon beam Evolution in Silicon dioxide”; “ERMES 20.0” is the fully developed concept	Seismology, muon beam interrogation
ERMES 20.0	Open-source, frequency-domain finite element code for computational electromagnetics in plasmas	Fusion plasma modeling

In the geophysical usage, “ERMES 20.0” is explicitly described as the fully developed version of a previously introduced concept: multi-TeV muon beamlines aligned with major faults, long rock traversals up to about $3\ \mathrm{km}$, muonic lenses, precision tracking detectors, and continuous integration into seismic hazard systems (Serafini et al., 15 Jan 2026). In the plasma-physics usage, ERMES 20.0 is a validated numerical platform for solving the frequency-domain cold-plasma Maxwell problem in the O–X stage of O–X–B conversion, with benchmark conditions representative of MAST Upgrade (Otin et al., 30 Sep 2025).

A common misconception would be to treat the shared name as evidence of a shared methodology. The published descriptions instead place the two usages in different physical regimes, with different governing equations, infrastructures, and validation criteria.

2. ERMES as active interrogation of stressed quartz near seismogenic faults

The geophysical ERMES concept is an active geophysical interrogation system rather than a passive precursor detector. Its central aim is to monitor and detect stable and reliable precursor signals on an adequate time scale, well before an earthquake event, by continuously measuring the time evolution of tectonic stress inside quartz-rich granite-like rocks surrounding a known seismogenic fault (Serafini et al., 15 Jan 2026).

The underlying mechanism is piezoelectric. Quartz crystals under tectonic stress of order $1$–$5\ \mathrm{kbar}$ develop electric fields of order a few $\mathrm{MV/m}$, parameterized in the analysis as

$E_p \simeq 5~\mathrm{MV/m\ per\ kbar\ of\ stress},$

so that for stresses of order $1$–$5\ \mathrm{kbar}$,

$E_p \sim 5\text{–}25~\mathrm{MV/m}.$

At the constitutive level, the crystal response is written as

$P_i=d_{ijk}\sigma_{jk},$

and at coarse-grained level the local electric field scales linearly with applied stress, $E_p\propto \sigma$ (Serafini et al., 15 Jan 2026).

ERMES assumes the least favorable case of totally incoherent quartz orientations. Under that assumption, each traversed quartz crystal gives a very small transverse kick to a relativistic muon, and the aggregate effect over many crystals is a piezoelectric random walk superimposed on standard multiple Coulomb scattering. The measurable signature is not a mean deflection but a variance increase in transverse phase space, expressed through the relative increase of the RMS beam spot size,

$1$0

For incoherent fabrics, the model gives

$1$1

with $1$2 in $1$3 and $1$4 in meters. Because $1$5, the signal scales as $1$6 (Serafini et al., 15 Jan 2026).

For the representative parameters stated in the paper—tectonic stress $1$7–$1$8, $1$9–$5\ \mathrm{kbar}$0, and centimeter-scale crystals $5\ \mathrm{kbar}$1—the predicted relative spot-size increase is

$5\ \mathrm{kbar}$2

The system therefore hinges on resolving beam-width changes from tens of ppm to below $5\ \mathrm{kbar}$3 against a much larger background broadening from rock traversal (Serafini et al., 15 Jan 2026).

The conceptual significance of this architecture lies in where it probes the piezoelectric process. Conventional electromagnetic precursor searches attempt to detect far-field electromagnetic emissions outside the crust. ERMES instead interrogates the near field inside the source region, namely quartz grains embedded in the stressed rock volume.

3. Beam propagation, muonic optics, and feasibility constraints

ERMES imposes unusually demanding beam requirements because the relevant pre-rupture “plastic loading” phase is taken to last tens of minutes. The proposed operational target is one measurement of $5\ \mathrm{kbar}$4 every $5\ \mathrm{kbar}$5, with 1–2 minute sampling, using $5\ \mathrm{kbar}$6 detected muons per measurement. After accounting for survival fraction and acceptance, this translates into an average beam intensity at rock entry of $5\ \mathrm{kbar}$7 muons/s (Serafini et al., 15 Jan 2026).

Muon range places a hard limit on accessible depth. The paper summarizes that below a few TeV the energy loss is dominated by ionization, roughly $5\ \mathrm{kbar}$8–$5\ \mathrm{kbar}$9 in rock, whereas above $\mathrm{MV/m}$0 radiative losses dominate. The Continuous Slowing Down Approximation range is stated as $\mathrm{MV/m}$1 for $\mathrm{MV/m}$2 muons in granite-like rock and $\mathrm{MV/m}$3 for $\mathrm{MV/m}$4 muons in $\mathrm{MV/m}$5, with limited benefit in exceeding $\mathrm{MV/m}$6 because radiative losses saturate the range (Serafini et al., 15 Jan 2026).

The FLUKA study for $\mathrm{MV/m}$7 muons crossing $\mathrm{MV/m}$8 of $\mathrm{MV/m}$9 starts from a monochromatic beam with transverse RMS $E_p \simeq 5~\mathrm{MV/m\ per\ kbar\ of\ stress},$0. The reported results are a survival fraction of $E_p \simeq 5~\mathrm{MV/m\ per\ kbar\ of\ stress},$1, an exit energy spectrum with average $E_p \simeq 5~\mathrm{MV/m\ per\ kbar\ of\ stress},$2, and a final transverse RMS of order $E_p \simeq 5~\mathrm{MV/m\ per\ kbar\ of\ stress},$3 in $E_p \simeq 5~\mathrm{MV/m\ per\ kbar\ of\ stress},$4, with $E_p \simeq 5~\mathrm{MV/m\ per\ kbar\ of\ stress},$5. A separate Geant4 study for $E_p \simeq 5~\mathrm{MV/m\ per\ kbar\ of\ stress},$6 muons through $E_p \simeq 5~\mathrm{MV/m\ per\ kbar\ of\ stress},$7 of $E_p \simeq 5~\mathrm{MV/m\ per\ kbar\ of\ stress},$8 finds mean energy loss $E_p \simeq 5~\mathrm{MV/m\ per\ kbar\ of\ stress},$9, primary survival $1$0, strong spectrum broadening, and negligible secondary muons at the measurement depth: only $1$1 of primary count survive beyond $1$2, and they are negligible at $1$3 (Serafini et al., 15 Jan 2026).

Because the post-rock beam has meter-scale spot size and significant divergence, the paper introduces a muonic lens: a large-radius, current-carrying cylindrical conductor through which the muons propagate in solid aluminum rather than vacuum. For a cylinder of radius $1$4 and current density $1$5, the magnetostatic field is

$1$6

Inside the lens, the field is linear in radius and provides simultaneous cylindrical focusing in both transverse directions. The focusing gradient is

$1$7

and the thin-lens focal length is

$1$8

For $1$9, $5\ \mathrm{kbar}$0, $5\ \mathrm{kbar}$1, and $5\ \mathrm{kbar}$2, the focal length is $5\ \mathrm{kbar}$3 (Serafini et al., 15 Jan 2026).

The lens is presented as a practical necessity for any “ERMES 20.0” infrastructure. Simulated performance includes a focal point $5\ \mathrm{kbar}$4 downstream of the $5\ \mathrm{kbar}$5 exit and minimum RMS radius $5\ \mathrm{kbar}$6, as well as a muonic FODO-like lattice with 8 lenses transporting a $5\ \mathrm{kbar}$7 beam over $5\ \mathrm{kbar}$8 (Serafini et al., 15 Jan 2026). At the same time, the feasibility limits are explicit: collider-class infrastructure, underground caverns, geological characterization, high stability in beam and detector systems, and long-term calibration against systematic drift.

4. ERMES 20.0 as a finite element code for cold-plasma full-wave modeling

In plasma modeling, ERMES 20.0 is an open-source, frequency-domain finite element code for computational electromagnetics in plasmas. The validation study uses it for Ordinary–Extraordinary mode conversion in the Electron Bernstein Wave regime, specifically the first O$5\ \mathrm{kbar}$9X stage of the O–X–B scheme, under a cold-plasma slab benchmark representative of the MAST Upgrade EBW system (Otin et al., 30 Sep 2025).

The solved continuous problem is the time-harmonic Maxwell system with fields $E_p \sim 5\text{–}25~\mathrm{MV/m}.$0, $E_p \sim 5\text{–}25~\mathrm{MV/m}.$1:

$E_p \sim 5\text{–}25~\mathrm{MV/m}.$2

which yields the curl–curl equation

$E_p \sim 5\text{–}25~\mathrm{MV/m}.$3

Here $E_p \sim 5\text{–}25~\mathrm{MV/m}.$4, and $E_p \sim 5\text{–}25~\mathrm{MV/m}.$5 is the cold electron plasma tensor in a uniform magnetic field, with collision frequency incorporated via the complex-frequency substitution $E_p \sim 5\text{–}25~\mathrm{MV/m}.$6 in the electron contribution (Otin et al., 30 Sep 2025).

The benchmark geometry is quasi-2D in $E_p \sim 5\text{–}25~\mathrm{MV/m}.$7–$E_p \sim 5\text{–}25~\mathrm{MV/m}.$8, implemented as a thin $E_p \sim 5\text{–}25~\mathrm{MV/m}.$9-slab of thickness $P_i=d_{ijk}\sigma_{jk},$0 in a 3D solver. The stated dimensions are $P_i=d_{ijk}\sigma_{jk},$1, $P_i=d_{ijk}\sigma_{jk},$2, with $P_i=d_{ijk}\sigma_{jk},$3 and $P_i=d_{ijk}\sigma_{jk},$4. The electron density profile is vacuum for $P_i=d_{ijk}\sigma_{jk},$5, and for $P_i=d_{ijk}\sigma_{jk},$6,

$P_i=d_{ijk}\sigma_{jk},$7

with the dimensionless gradient parameter $P_i=d_{ijk}\sigma_{jk},$8 scanned from $P_i=d_{ijk}\sigma_{jk},$9 to $E_p\propto \sigma$0 (Otin et al., 30 Sep 2025).

Two formulations are used. The first is a standard double-curl electric-field formulation with Nédélec edge elements, denoted EDG. The second is a regularized Maxwell formulation with nodal elements, denoted RME, which adds a divergence-penalty term and is described as unique to ERMES. The study uses second-order elements throughout. The unstructured 3D mesh has characteristic element size $E_p\propto \sigma$1 and about $E_p\propto \sigma$2 million elements; each case is solved on an HPC cluster with $E_p\propto \sigma$3–$E_p\propto \sigma$4 memory, $E_p\propto \sigma$5–$E_p\propto \sigma$6 hours runtime, and a direct MUMPS solver via PETSc (Otin et al., 30 Sep 2025).

This formulation choice is central to the validation. EDG is standard but becomes poorly conditioned near resonances in anisotropic media, while RME is introduced to improve conditioning and produce stable solutions even near resonances, though the paper explicitly notes an open question as to whether regularization may numerically smooth the solution in some regimes.

5. Boundary conditions, launch model, and quantitative validation

ERMES 20.0 uses generalized Robin boundary conditions for electric-field absorption and source injection,

$E_p\propto \sigma$7

with an additional scalar Robin condition for the RME formulation. On the top and side boundaries, $E_p\propto \sigma$8 and the admittance $E_p\propto \sigma$9 is chosen to absorb the appropriate cold-plasma mode. On the input and output planes, $1$00, while the Gaussian beam source is imposed at the input boundary through nonzero $1$01 and $1$02 built from the incident beam field $1$03 (Otin et al., 30 Sep 2025).

The injected beam is a pure O-mode Gaussian beam launched at the optimal O–X conversion angle

$1$04

which gives $1$05 for $1$06 and $1$07. The beam waist radii are $1$08, and the waist center is located at $1$09 (Otin et al., 30 Sep 2025).

The principal scalar validation metric is the reflection coefficient,

$1$10

with the paper stating that in this slab setup the reflection coefficient is equal to one minus the mode conversion efficiency (Otin et al., 30 Sep 2025).

A key numerical issue is damping. EDG is reported to be numerically unstable for $1$11, while for $1$12 it stabilizes and agrees with RME. RME yields smooth reflection curves for all tested $1$13. The authors adopt $1$14 as the minimal viable damping for meaningful cross-code comparison (Otin et al., 30 Sep 2025).

Validation is performed against four independent full-wave solvers: IPF-FDMC, EMIT-2D, CUWA-2D, and the Fourier-based FFW code. For $1$15, ERMES reproduces the beam structure, O-mode cutoff location, converted X-mode propagation, penetration depth, and interference pattern seen in IPF-FDMC and FFW. For $1$16, the magnitude of the Poynting vector agrees well with CUWA in the direction and focusing of the injected O-mode beam, the conversion region, and the reflected power pattern. Across the full scan $1$17, ERMES, IPF-FDMC, EMIT-2D, and FFW show excellent agreement in $1$18, while CUWA exhibits systematically higher reflection attributed to different damping models; the CUWA curve nonetheless falls within the ERMES band obtained by varying $1$19 (Otin et al., 30 Sep 2025).

6. Significance, limitations, and prospective development

The two meanings of ERMES 20.0 occupy different positions on the spectrum from validated tool to speculative infrastructure. The plasma-wave ERMES 20.0 is already presented as a validated numerical code with demonstrated agreement in electric-field structure, Poynting-vector distribution, and reflection coefficient against established FDTD and Fourier full-wave solvers (Otin et al., 30 Sep 2025). The geophysical ERMES 20.0, by contrast, is explicitly characterized as a long-term R&D direction rather than an immediately deployable technology, despite quantitatively developed beam, transport, and detection models (Serafini et al., 15 Jan 2026).

For the geophysical system, the principal limitations are infrastructural and metrological: the need for a high-energy muon accelerator up to $1$20, beamline infrastructure colocated with seismogenic faults, underground entry and exit caverns, stable detector systems able to resolve $1$21, and detailed characterization of quartz content, crystal size distribution, fabric orientation statistics, and stress fields along the beam path. Open questions identified by the authors include detector technology, full 3D modeling of muonic lenses with fringing fields, engineering feasibility of large-radius high-current aluminum conductors, geological modeling of non-incoherent piezoelectric scenarios, accelerator cost, and laboratory validation of the piezoelectric random-walk model (Serafini et al., 15 Jan 2026).

For the finite element code, the main technical limitations are numerical stability near cold-plasma resonances for the EDG formulation, the heavy computational cost of the quasi-2D 3D implementation, and the unresolved question of whether the RME divergence term may in some settings smooth the solution rather than only regularize it. At the same time, the code’s reported advantages are geometric flexibility, second-order unstructured 3D discretization, multiple formulations, and boundary conditions tailored to plasma modes and Gaussian-beam launch. The stated development direction is extension to warm and hot plasma effects, with improved predictive modeling of electromagnetic wave heating and current drive in next-generation fusion devices (Otin et al., 30 Sep 2025).

Taken together, the two usages of ERMES 20.0 illustrate a naming convergence rather than a scientific convergence. One denotes a proposed method for converting stress-dependent piezoelectric microphysics into a beam-optics observable in the Earth’s crust; the other denotes a validated finite element framework for anisotropic full-wave plasma simulation. Their commonality is terminological, whereas their technical content, evidence base, and maturity are distinct.