Simulation and optimization of the Active Magnetic Shield of the n2EDM experiment
Pith reviewed 2026-05-16 09:31 UTC · model grok-4.3
The pith
Finite element simulation of the active magnetic shield matches measurements and optimizes feedback sensor placement via genetic algorithms.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We present a full finite element simulation of magnetic fields generated by the AMS in the presence of the MSR. This simulation is of sufficient accuracy to approach our measurements. We demonstrate how the simulation can be used with an example, obtaining an optimal number and placement of feedback sensors using genetic algorithms.
What carries the argument
Finite element simulation of the AMS coils that includes the MSR shielding, paired with genetic-algorithm search for optimal feedback sensor count and locations.
Load-bearing premise
The finite element simulation accurately captures the magnetic field behavior in the presence of the MSR without needing additional corrections for real-world imperfections.
What would settle it
If measured residual fields after applying the genetically optimized sensor placement differ by more than a few microtesla from the simulation predictions inside the n2EDM volume, the claimed accuracy would be falsified.
read the original abstract
The n2EDM experiment at the Paul Scherrer Institute aims to conduct a high-sensitivity search for the electric dipole moment of the neutron. Magnetic stability and control are achieved through a combination of passive shielding, provided by a magnetically shielded room (MSR), and a surrounding active field compensation system by an Active Magnetic Shield (AMS). The AMS is a feedback-controlled system of eight coils spanned on an irregular grid, designed to provide magnetic stability to the enclosed volume by actively suppressing external magnetic disturbances. It can compensate static and variable magnetic fields up to $\pm 50$ $\mu$T (homogeneous components) and $\pm 5$ $\mu$T/m (first-order gradients), suppressing them to a few $\mu$T in the sub-Hertz frequency range. We present a full finite element simulation of magnetic fields generated by the AMS in the presence of the MSR. This simulation is of sufficient accuracy to approach our measurements. We demonstrate how the simulation can be used with an example, obtaining an optimal number and placement of feedback sensors using genetic algorithms.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents a full finite element simulation of the magnetic fields generated by the Active Magnetic Shield (AMS) in the presence of the magnetically shielded room (MSR) for the n2EDM experiment. It claims this simulation is accurate enough to approach experimental measurements and demonstrates its application by using genetic algorithms to optimize the number and placement of feedback sensors for compensating static and variable fields up to ±50 μT (homogeneous) and ±5 μT/m (gradients) in the sub-Hz range.
Significance. If the simulation is shown to reproduce measurements with quantified accuracy, the work would supply a practical design tool for active compensation systems in precision neutron EDM searches, allowing simulation-driven optimization of sensor arrays to maintain field stability against external disturbances.
major comments (2)
- [Abstract] Abstract: the claim that the simulation 'is of sufficient accuracy to approach our measurements' is unsupported by any quantitative metrics (RMS error, peak deviation, or residual spectra for homogeneous and gradient components). This is load-bearing for the GA optimization, as unquantified mismatch near coil-MSR interfaces could shift the reported optimal sensor count and placement.
- [Validation section] Validation section (inferred from abstract description of FEM model): no error analysis, cross-validation on gradient components, or details on model adjustments for real-world imperfections (e.g., permeability variations or eddy-current effects) are provided, leaving the central fidelity claim unverified.
minor comments (2)
- [Abstract] Abstract: the sub-Hertz frequency range should be quantified (e.g., 0.01–1 Hz) to allow direct comparison with the compensation performance figures.
- [Optimization section] Optimization example: the genetic algorithm fitness function, population size, and convergence criteria are not specified, reducing reproducibility of the sensor-placement result.
Simulated Author's Rebuttal
We thank the referee for their thorough review and valuable comments on our manuscript. We address each of the major comments below and have updated the manuscript to include the requested quantitative validations and details.
read point-by-point responses
-
Referee: [Abstract] Abstract: the claim that the simulation 'is of sufficient accuracy to approach our measurements' is unsupported by any quantitative metrics (RMS error, peak deviation, or residual spectra for homogeneous and gradient components). This is load-bearing for the GA optimization, as unquantified mismatch near coil-MSR interfaces could shift the reported optimal sensor count and placement.
Authors: We agree with this observation. The current manuscript relies on visual agreement in figures for the accuracy claim, without explicit quantitative metrics. In the revised manuscript, we will add RMS errors, peak deviations, and residual spectra for both homogeneous and gradient components. This will provide the necessary support for the simulation's fidelity and its use in the genetic algorithm optimization. We do not anticipate changes to the optimal sensor number and placement based on preliminary checks, but the added metrics will allow readers to assess this independently. revision: yes
-
Referee: [Validation section] Validation section (inferred from abstract description of FEM model): no error analysis, cross-validation on gradient components, or details on model adjustments for real-world imperfections (e.g., permeability variations or eddy-current effects) are provided, leaving the central fidelity claim unverified.
Authors: We acknowledge the need for more rigorous validation details. The revised version will include a dedicated error analysis section with cross-validation on gradient components. Additionally, we will provide details on how the FEM model was adjusted for real-world imperfections, including any considerations for permeability variations and eddy-current effects in the sub-Hz range. If certain effects were neglected, we will justify the approximation and its expected impact on the results. revision: yes
Circularity Check
No significant circularity: FEM model validated externally and GA optimization is downstream
full rationale
The paper's central steps are (1) construction of a finite-element model of the AMS coils inside the MSR and (2) use of that model to run a genetic algorithm that selects sensor number and placement. The abstract states the simulation 'is of sufficient accuracy to approach our measurements,' but this is presented as an empirical validation claim rather than a self-referential definition. No equation or procedure is shown to reduce to its own fitted parameters by construction, no uniqueness theorem is imported from prior self-citations, and the GA output is a downstream application of the model rather than a re-derivation of the model's inputs. The derivation chain therefore remains self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Finite element method accurately models the magnetostatic fields generated by the AMS coils in the presence of the MSR
Reference graph
Works this paper leans on
-
[1]
Pignol,A global perspective on searches for Electric Dipole Moments,arXiv(2019)
G. Pignol,A global perspective on searches for Electric Dipole Moments,arXiv(2019)
work page 2019
-
[2]
R. Alarcon, J. Alexander, V. Anastassopoulos, T. Aoki, R. Baartman, S. Baeßler et al.,Electric dipole moments and the search for new physics,arXiv(2022)
work page 2022
-
[3]
M.S. Safronova, D. Budker, D. DeMille, D.F.J. Kimball, A. Derevianko and C.W. Clark,Search for new physics with atoms and molecules,Rev. Mod. Phys.90(2018) 025008
work page 2018
-
[4]
Jungmann,Searching for electric dipole moments,Ann
K. Jungmann,Searching for electric dipole moments,Ann. Phys.525(2013) 550
work page 2013
-
[5]
C. Abel, N. Ayres, G. Ban, G. Bison, K. Bodek, V. Bondar et al.,A search for neutron to mirror-neutron oscillations using the nEDM apparatus at PSI,Physics Letters B812(2021) 135993
work page 2021
-
[6]
V. Santoro, D. Milstead, P. Fierlinger, W.M. Snow, J. Barrow, M. Bartis et al.,The HIBEAM program: search for neutron oscillations at the European Spallation Source,J. Phys. G: Nucl. Part. Phys. (2025) 040501. – 16 –
work page 2025
- [7]
-
[8]
K. Kobayashi, A. Kon, M. Yoshizawa and Y. Uchikawa,Active magnetic shielding using symmetric magnetic field sensor method,IEEE Transactions on Magnetics48(2012) 4554
work page 2012
-
[9]
T. Bryś, S. Czekaj, M. Daum, P. Fierlinger, D. George, R. Henneck et al.,Magnetic field stabilization formagneticallyshieldedvolumesbyexternalfieldcoils,NuclearInstrumentsandMethodsinPhysics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment554(2005) 527
work page 2005
-
[10]
H.J.M. ter Brake, H.J. Wieringa and H. Rogalla,Improvement of the performance of a mu -metal magnetically shielded room by means of active compensation (biomagnetic applications),Meas. Sci. Technol.2(1991) 596
work page 1991
-
[11]
T. Spemann, D.and Reinert, J. Vogt, J. Wassermann and B. T.,Active compensation of stray magnetic fields at LIPSION,Nucl. Instrum. Methods. Phys. Res. B210(2003) 79
work page 2003
- [12]
-
[13]
Y. Liu, J. Yang, F. Cao, X. Zhang and S. Zheng,Enhancement of magnetic shielding based on low-noise materials, magnetization control, and active compensation: A review,Materials17(2024) 5469
work page 2024
- [14]
-
[15]
B. Lauss and on behalf of the PSI UCN Project Team,Startup of the high-intensity ultracold neutron source at the Paul Scherrer Institute,Hyperfine Interact.211(2012) 21
work page 2012
- [16]
-
[17]
B. Lauss and B. Blau,UCN, the ultracold neutron source – neutrons for particle physics, p. 004, SciPost, 2021, DOI
work page 2021
- [18]
- [19]
-
[20]
J. Elen, W. Franken, I. Horvath, G. Pasotti, M. Ricci, J. Roeterdink et al.,The superconductor test facility sultan,IEEE Transactions on Magnetics17(1981) 490
work page 1981
-
[21]
J. Schippers, R. Dölling, J. Duppich, G. Goitein, M. Jermann, A. Mezger et al.,The SC cyclotron and beam lines of PSI’s new protontherapy facility PROSCAN,Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms261(2007) 773
work page 2007
-
[22]
S. Emmenegger,Next Generation Active Magnetic Shielding for n2EDM and Axion-Like Particle Search, Doctoral Thesis, ETH Zurich, 2021. https://doi.org/10.3929/ethz-b-000515206
-
[23]
Rawlik, A.Eggenberger, J.Krempel, C.Crawford, K
M. Rawlik, A.Eggenberger, J.Krempel, C.Crawford, K. Kirch, F.M.Piegsaet al.,A simplemethodof coil design,Am. J. Phys.86(2018) 602. – 17 –
work page 2018
-
[24]
Rawlik,Active Magnetic Shielding and Axion-Dark-Matter Search, Doctoral Thesis, ETH Zurich,
M. Rawlik,Active Magnetic Shielding and Axion-Dark-Matter Search, Doctoral Thesis, ETH Zurich,
-
[25]
https://doi.org/10.3929/ethz-b-000273039
-
[26]
B. Franke,Investigations of the internal and external magnetic fields of the neutron electric dipole moment experiment at the Paul Scherrer Institute, Doctoral Thesis, ETH Zurich, 2013. https://doi.org/10.3929/ethz-a-010144564
- [27]
- [28]
-
[29]
Stenderm,Python-MPH: A python interface to COMSOL’s Java API, 2023
N. Stenderm,Python-MPH: A python interface to COMSOL’s Java API, 2023
work page 2023
-
[30]
E. Cartan and R. Hermann,Geometry of Riemannian Spaces, Math Sci Press, Massachusetts, translation or edition ed. (1983)
work page 1983
-
[31]
D.E. Goldberg and J.H. Holland,Genetic algorithms and machine learning.,Machine Learning 3 (1988)
work page 1988
-
[32]
J.H. Holland,Adaptation in natural and artificial systems: An introductory analysis with applications to biology, control, and artificial intelligence,The MIT Press(1992)
work page 1992
-
[33]
R.L. Haupt and S.E. Haupt,Practical genetic algorithms,John Wiley & Sons, Inc.(2003)
work page 2003
-
[34]
G. Heeralal, S.K. Chaturvedi and A.K. Thakur,Design and optimization of multilayered electromagnetic shield using a real-coded genetic algorithm,Progress In Electromagnetics Research B 39(39)(2012)
work page 2012
-
[35]
D. Pais,Development of the caesium magnetometer array for the n2EDM experiment, Doctoral Thesis, ETH Zurich, 2018. https://doi.org/10.3929/ethz-b-000511496
-
[36]
N.R. Shaw and R.E. Ansorge,Genetic algorithms for MRI magnet design,IEEE Transactions on Applied Superconductivity(2002)
work page 2002
- [37]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.