pith. sign in

arxiv: 2601.22960 · v1 · submitted 2026-01-30 · ⚛️ physics.ins-det · nucl-ex

Simulation and optimization of the Active Magnetic Shield of the n2EDM experiment

Pith reviewed 2026-05-16 09:31 UTC · model grok-4.3

classification ⚛️ physics.ins-det nucl-ex
keywords active magnetic shieldfinite element simulationmagnetic shielding roomgenetic algorithmfeedback sensorsneutron EDMn2EDM experimentfield compensation
0
0 comments X

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.

The n2EDM experiment needs tight magnetic field control inside its shielded volume to reach high sensitivity in the search for the neutron electric dipole moment. An active magnetic shield of eight coils provides feedback compensation for external disturbances, working together with a passive magnetically shielded room. This paper builds a full finite element model that calculates the fields from the coils while including the effect of the room. The model reproduces measured fields closely enough that it can replace physical trial-and-error. The authors then use the model inside a genetic algorithm to find the best number and positions for the feedback sensors.

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.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

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)
  1. [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.
  2. [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)
  1. [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.
  2. [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

2 responses · 0 unresolved

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
  1. 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

  2. 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

0 steps flagged

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

0 free parameters · 1 axioms · 0 invented entities

Central claim rests on standard assumptions of finite element magnetostatics and material properties for the MSR; no explicit free parameters or invented entities are described in the abstract.

axioms (1)
  • domain assumption Finite element method accurately models the magnetostatic fields generated by the AMS coils in the presence of the MSR
    Invoked implicitly when claiming the simulation approaches measurements.

pith-pipeline@v0.9.0 · 5771 in / 1184 out tokens · 36346 ms · 2026-05-16T09:31:27.906234+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

37 extracted references · 37 canonical work pages

  1. [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)

  2. [2]

    Alarcon, J

    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)

  3. [3]

    Safronova, D

    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

  4. [4]

    Jungmann,Searching for electric dipole moments,Ann

    K. Jungmann,Searching for electric dipole moments,Ann. Phys.525(2013) 550

  5. [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

  6. [6]

    Santoro, D

    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 –

  7. [7]

    Lin, Y.-L

    H.-Y. Lin, Y.-L. Song and L.-M. Chang,A magnetic field cancelling system design for mitigating extremely low frequency magnetic field in a high tech fab, in2021 Asia-Pacific International Symposium on Electromagnetic Compatibility (APEMC), pp. 1–4, 2021, DOI

  8. [8]

    Kobayashi, A

    K. Kobayashi, A. Kon, M. Yoshizawa and Y. Uchikawa,Active magnetic shielding using symmetric magnetic field sensor method,IEEE Transactions on Magnetics48(2012) 4554

  9. [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

  10. [10]

    ter Brake, H.J

    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

  11. [11]

    Spemann, D.and Reinert, J

    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

  12. [12]

    Holmes, M

    N. Holmes, M. Rea, J. Chalmers, J. Leggett, L.J. Edwards, P. Nell et al.,A lightweight magnetically shielded room with active shielding,Sci. Rep.12(2022) 13561

  13. [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

  14. [14]

    Ayres, G

    N.J. Ayres, G. Ban, L. Bienstman, G. Bison, K. Bodek, V. Bondar et al.,The design of the n2EDM experiment,Eur. Phys. J. C81(2021) 512

  15. [15]

    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

    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

  16. [16]

    Bison, M

    G. Bison, M. Daum, K. Kirch, B. Lauss, D. Ries, P. Schmidt-Wellenburg et al.,Ultracold neutron storage and transport at the PSI UCN source,Eur. Phys. J. A58(2022) 103

  17. [17]

    Lauss and B

    B. Lauss and B. Blau,UCN, the ultracold neutron source – neutrons for particle physics, p. 004, SciPost, 2021, DOI

  18. [18]

    Ayres, G

    N.J. Ayres, G. Ban, G. Bison, K. Bodek, V. Bondar, T. Bouillaud et al.,The very large n2EDM magnetically shielded room with an exceptional performance for fundamental physics measurements, Rev. Sci. Instrum.93(2022) 095105

  19. [19]

    Abel, N.J

    C. Abel, N.J. Ayres, G. Ban, G. Bison, K. Bodek, V. Bondar et al.,A large ‘Active Magnetic Shield’ for a high-precision experiment,Eur. Phys. J. C83(2023) 1061

  20. [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

  21. [21]

    Schippers, R

    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

  22. [22]

    Emmenegger,Next Generation Active Magnetic Shielding for n2EDM and Axion-Like Particle Search, Doctoral Thesis, ETH Zurich, 2021

    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. [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 –

  24. [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. [25]

    https://doi.org/10.3929/ethz-b-000273039

  26. [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. [27]

    Afach, G

    S. Afach, G. Bison, K. Bodek, F. Burri, Z. Chowdhuri, M. Daum et al.,Dynamic stabilization of the magnetic field surrounding the neutron electric dipole moment spectrometer at the Paul Scherrer Institute,J. Appl. Phys.116(2014) 084510

  28. [28]

    COMSOL Multiphysics©

    COMSOL AB, “COMSOL Multiphysics©.”

  29. [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

  30. [30]

    Cartan and R

    E. Cartan and R. Hermann,Geometry of Riemannian Spaces, Math Sci Press, Massachusetts, translation or edition ed. (1983)

  31. [31]

    Goldberg and J.H

    D.E. Goldberg and J.H. Holland,Genetic algorithms and machine learning.,Machine Learning 3 (1988)

  32. [32]

    Holland,Adaptation in natural and artificial systems: An introductory analysis with applications to biology, control, and artificial intelligence,The MIT Press(1992)

    J.H. Holland,Adaptation in natural and artificial systems: An introductory analysis with applications to biology, control, and artificial intelligence,The MIT Press(1992)

  33. [33]

    Haupt and S.E

    R.L. Haupt and S.E. Haupt,Practical genetic algorithms,John Wiley & Sons, Inc.(2003)

  34. [34]

    Heeralal, S.K

    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)

  35. [35]

    Pais,Development of the caesium magnetometer array for the n2EDM experiment, Doctoral Thesis, ETH Zurich, 2018

    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. [36]

    Shaw and R.E

    N.R. Shaw and R.E. Ansorge,Genetic algorithms for MRI magnet design,IEEE Transactions on Applied Superconductivity(2002)

  37. [37]

    Akiba, S

    T. Akiba, S. Sano, T. Yanase, T. Ohta and M. Koyama,Optuna: A next-generation hyperparameter optimization framework, inProceedings of the 25th ACM SIGKDD International Conference on Knowledge Discovery & Data Mining, pp. 2623–2631, ACM, 2019, DOI. – 18 –