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

Dit artikel presenteert een nauwkeurige eindige-elementen-simulatie van de Actieve Magnetische Schild van het n2EDM-experiment, waarbij de bruikbaarheid ervan wordt aangetoond voor het optimaliseren van de plaatsing en het aantal feedbacksensoren via genetische algoritmen om magnetische stabiliteit binnen de magnetisch afgeschermde ruimte te waarborgen.

De kern

Het Grote Geheel: Een Neutron Rustig Houden

Stel je voor dat je probeert een zeer delicaat, draaiend tolletje (de neutron) op een tafel in evenwicht te houden. Als de kamer schudt, of als er vlakbij een enorme ventilator aanzet, gaat de tol wiebelen en valt hij om. Wetenschappers willen dit tolletje bestuderen om te zien of het een klein, verborgen "schuine stand" heeft (een elektrisch dipoolmoment genoemd) die geheimen over het universum zou kunnen verklaren.

Om dit te doen, moeten ze de kamer perfect stil houden en moet de magnetische "wind" perfect rustig zijn. Het n2EDM-experiment bij het Paul Scherrer Instituut is deze hoog-risico kamer.

Het Probleem: Een Luidruchtige Buurt

Het experiment bevindt zich in een drukke wetenschappelijke buurt. In de buurt staan enorme supergeleidende magneten (zoals de SULTAN- en COMET-machines) die fungeren als gigantische elektromagneten. Wanneer deze machines op- of afregelen, creëren ze enorme magnetische "stormen" die de delicate meting van de neutron volledig zouden verpesten.

De Oplossing: Een Dubbel-laags Verdedigingssysteem

Om de kamer rustig te houden, bouwden de wetenschappers een tweeledig verdedigingssysteem:

1. Het Passieve Schild (Het Fort): Ze bouwden een speciale kamer, een Magnetically Shielded Room (MSR). Denk hierbij aan een fort gemaakt van zeven lagen van een supermagnetisch metaal genaamd mu-metall. Het werkt als een dik, zwaar deken dat het grootste deel van de magnetische ruis van buitenaf absorbeert.
2. Het Actieve Schild (De Ruisonderdrukkende Koptelefoon): Zelfs het beste deken heeft kleine lekken. Om dit op te lossen, voegden ze een Actief Magnetisch Schild (AMS) toe.
   • Hoe het werkt: Stel je voor dat de MSR wordt omringd door acht gigantische, onzichtbare "magnetische handen" (spoelen).
   • De Sensoren: Kleine apparaten genaamd fluxgate-sensoren (zoals tiny magnetische oren) zijn rond de kamer geplaatst. Ze luisteren naar de magnetische ruis.
   • De Feedbacklus: Wanneer de "oren" een verstoring horen (zoals een nabijgelegen magneet die opregelt), vertelt een computer de "handen" om direct terug te duwen met een gelijke en tegengestelde magnetische kracht. Het is precies zoals ruisonderdrukkende koptelefoons: ze horen de buitenste ruis en genereren een "anti-ruis" om deze perfect op te heffen.

De Uitdaging: Het Schild Verandert het Geluid

De wetenschappers realiseerden zich dat het "fort" (de mu-metall kamer) niet alleen ruis blokkeert; het vervormt het ook.

• De Analogie: Stel je voor dat je schreeuwt in een grot. De muren van de grot kaatsen het geluid rond, waardoor het in de hoeken harder echoot en in het midden stiller is.
• De Realiteit: De mu-metall wanden van de MSR buigen de magnetische velden. Dit betekent dat de magnetische "ruis" niet uniform is; het wordt versterkt in de hoeken van de kamer. Als de wetenschappers gewoon zouden gokken waar ze hun "oren" (sensoren) moeten plaatsen, zouden ze misschien de luidste plekken missen of proberen ruis op te heffen die er eigenlijk niet is.

De Simulatie: Een Virtuele Tweeling

Om dit op te lossen, bouwde het team een digitale tweeling van hun hele experiment met behulp van computersoftware (COMSOL).

• Ze creëerden een virtuele versie van het fort en de acht magnetische handen.
• Ze testten hoe de "handen" zouden terugduwen tegen de "ruis" terwijl het "fort" de golven vervormde.
• Het Resultaat: De computersimulatie kwam bijna perfect overeen met hun experimenten in de echte wereld. Dit bewees dat hun wiskunde klopte en dat het systeem zich op een voorspelbare, lineaire manier gedraagt (zoals een simpele volumeknop: draai het harder, het geluid wordt luider; draai het zachter, het wordt stiller).

De Optimalisatie: De Perfecte Plek Vinden

Zodra ze een werkende digitale tweeling hadden, vroegen ze zich af: "Waar is de absolute beste plek om onze magnetische oren te plaatsen?"

• De Oude Manier: Ze gebruikten een standaard algoritme om de posities te raden.
• De Nieuwe Manier: Ze gebruikten Genetische Algoritmen. Denk hierbij aan "digitale evolutie".
  • De computer creëerde duizenden willekeurige opstellingen van sensoren.
  • Het testte welke opstellingen het beste werkten bij het opheffen van ruis.
  • Het hield de "fitste" opstellingen (diegenen die ruis het beste opheften) en mengde ze samen om nog betere generaties te creëren.
  • Het Doel: Ze wilden de "conditienummer" minimaliseren. In gewone taal is dit een score die aangeeft hoe stabiel en makkelijk te controleren het systeem is. Een lagere score betekent dat het systeem minder snel in de war raakt of instabiel wordt.

Het Resultaat:
Het genetische algoritme vond een nieuwe opstelling van sensoren die wiskundig superieur was. Echter, de perfecte plek was fysiek onmogelijk te bouwen (er was niet genoeg ruimte). Dus kozen de wetenschappers de beste mogelijke plek die in de echte kamer paste.

• Ze verplaatsten de sensoren naar deze nieuwe plekken.
• Het systeem werkte precies zoals de computer voorspelde. Het "conditienummer" verbeterde, wat betekent dat het systeem nu stabieler is en beter in staat is om de magnetische stormen van nabijgelegen machines op te heffen.

Samenvatting

Het artikel beschrijft hoe wetenschappers een high-tech "ruisonderdrukkend" systeem bouwden voor een neutron-experiment. Ze realiseerden zich dat de kamer zelf de magnetische velden vervormde, dus bouwden ze een super-accurate computersimulatie om de vervorming te begrijpen. Met die simulatie en een "digitale evolutie"-algoritme vonden ze de perfecte plekken om hun sensoren te plaatsen om ervoor te zorgen dat het systeem stabiel blijft en succesvol enorme magnetische verstoringen van nabijgelegen machines kan opheffen.

Technische samenvatting

1. Probleemstelling

Het n2EDM-experiment aan het Paul Scherrer Instituut (PSI) heeft als doel het elektrische dipoolmoment van de neutron (nEDM) met extreme precisie te meten. Deze meting vereist een uiterst stabiel en uniform magnetisch veld binnen de opslagkamers voor Ultra-Koude Neutronen (UCN).

• Eisen: Tijdelijke fluctuaties in het magnetische veld moeten onder de 10 pT worden gehouden, en de verticale veldgradiënt over de kamers moet binnen 0,6 pT cm⁻¹ blijven gedurende een meetcyclus van ongeveer 180 seconden.
• Uitdagingen: Het experiment is gelegen in de buurt van faciliteiten met sterke, fluctuerende magnetische velden, met name de SULTAN-supraconductieve magneet (frequente op- en afbouwen) en de COMET-cyclotron (sporadische op- en afbouwen).
• Bestaande oplossing: Een passieve Magnetically Shielded Room (MSR) bestaande uit zeven lagen (zes mu-metal, één aluminium) biedt hoge demping (schildfactor $10^5$ tot $10^8$). Echter, passieve afscherming alleen is onvoldoende voor stabiliteit bij lage frequenties en voorkomt niet dat de buitenste laag van de MSR van magnetisatie verandert, wat naar binnen propageert.
• Het gat: Er is een Actieve Magnetische Schild (AMS) geïnstalleerd om externe storingen actief te compenseren. De interactie tussen de AMS-spoelen en het hoog-permeabele mu-metal van de MSR vervormt echter de magnetische velden, waardoor het regelsysteem complex wordt. Een nauwkeurige simulatie was nodig om deze interacties te begrijpen en het feedbackregelsysteem te optimaliseren, met name de plaatsing van fluxgate-sensoren.

2. Methodologie

De auteurs ontwikkelden een uitgebreide workflow die eindige-elementensimulatie, experimentele verificatie en algoritmische optimalisatie combineert.

A. Systeemontwerp

• AMS-architectuur: Acht feedback-gereguleerde spoelen, gerangschikt op een onregelmatig rooster rondom de MSR.
  • Spoeltypes: 3 spoelen voor homogene velden (X, Y, Z) en 5 spoelen voor gradiënten van de eerste orde.
  • Capaciteiten: Ontworpen om statische/variabele velden tot ±50 µT (homogeen) en ±5 µT/m (gradiënten) te compenseren.
• Feedbackmechanisme: Gebruik van acht 3-assige fluxgasmagnetometers (Sensys FGM3D/125). Het systeem maakt gebruik van een integraal regelfunctie gebaseerd op een lineair model:
  $$\mathbf{B} = \mathbf{B}_0 + \mathbf{M}\mathbf{I}$$
  Waarbij $\mathbf{B}$ het gemeten veld is, $\mathbf{B}_0$ de achtergrond, $\mathbf{I}$ de spoelstromen en $\mathbf{M}$ de proportionaliteitsmatrix. De regelstroom wordt berekend met behulp van de Moore-Penrose-pseudoinverse ($\mathbf{M}^\dagger$).
• Optimalisatiemetric: De stabiliteit van de feedbacklus hangt af van het conditiegetal ($\sigma$) van de matrix $\mathbf{M}^\dagger$. Een lagere $\sigma$ impliceert een stabieler en beter gedefinieerd regelsysteem.

B. Simulatie-implementatie (COMSOL)

• Modellering: Een volledig 3D eindige-elementenmodel werd gemaakt in COMSOL Multiphysics 6.1 met behulp van het Python-MPh-pakket.
• Geometrie: Het AMS-rooster (ongeveer $10,3 \times 8,6 \times 9,8$ m) en de MSR werden gemodelleerd.
• Vereenvoudiging: De zeslaagse mu-metal MSR werd benaderd als een enkele 5 cm dikke wand met een effectieve relatieve permeabiliteit ($\mu_r$) van 2000 cm⁻¹ (waardoor $\mu_{eff} = 10.000$). Simulaties bevestigden dat waarden boven deze drempel het gedrag van het externe veld niet significant veranderen.
• Lineariteitsverificatie: De simulatie bevestigde dat het systeem lineair gedraagt, waardoor de magnetische veldrespons kan worden geschaald vanuit een referentiestroom in plaats van opnieuw te simuleren voor elke stroomwaarde.

C. Optimalisatiestrategie

• Algoritme: Genetische Algoritmen (GA) werden gebruikt om een hoog-dimensionaal meerdoelig optimalisatieprobleem op te lossen.
• Doelstellingen:
  1. Minimaliseren van het conditiegetal ($\sigma$) van de feedbackmatrix.
  2. Minimaliseren van de gemiddelde afstand van sensoren tot de hoeken van de MSR (waar veldvervorming het grootst is).
• Parameter ruimte: De GA optimaliseerde het aantal sensoren ($n$), hun 3D-posities en hun oriëntatie (azimutale en polaire hoeken), waardoor een 6$n$-dimensionale zoekruimte ontstond.

3. Belangrijkste Resultaten

A. Simulatie versus Experiment

• Nauwkeurigheid: De COMSOL-simulatie toonde uitstekende overeenkomst met experimentele data.
  • Lineariteit: De relatie tussen spoelstroom en magnetisch veld werd bevestigd als lineair.
  • Conditiegetal: Het gesimuleerde conditiegetal ($\sigma_{sim} = 10,41$) kwam overeen met de experimentele waarde ($\sigma_{exp} = 10,18$) binnen 1,37 standaardafwijkingen, rekening houdend met onzekerheden in de sensorpositionering (±2,5 cm).
• Veldvervorming: De simulatie visualiseerde hoe de mu-metal MSR de homogene velden die door de AMS worden gegenereerd, vervormt, waardoor de veldsterkte bij de hoeken en randen van de kamer wordt versterkt (tot >10 µT residu-velden ondanks compensatie).

B. Optimalisatie-uitkomsten

• Pareto-front: De GA genereerde een Pareto-front die de afweging toont tussen het conditiegetal en de nabijheid van sensoren tot de hoeken van de MSR.
• Optimale configuratie:
  • De theoretische optimum suggereerde een conditiegetal van $\sigma \approx 6,3$.
  • Vanwege fysieke ruimtebeperkingen rondom het experiment werd een praktische configuratie geselecteerd met $\sigma = 7,85$.
  • Deze configuratie gebruikte 8 fluxgates die dichter bij de hoeken van de MSR waren geplaatst dan de vorige opstelling.
• Experimentele validatie: De nieuwe configuratie werd geïmplementeerd, wat resulteerde in een experimenteel conditiegetal van $\sigma_{exp} = 8,15$, waarmee de simulatie- en optimalisatieworkflow werd gevalideerd.

C. Prestaties

• De AMS onderdrukt externe magnetische storingen (bijv. SULTAN-rampen) succesvol tot het niveau van enkele µT in het sub-Hertz-frequentiebereik.
• Hoewel de nieuwe sensorconfiguratie resulteerde in iets hogere residu-velden op de sensorlocaties (door dichter bij de hoeken met hoge vervorming te zitten), verbeterde het aanzienlijk de regelbaarheid en stabiliteit van de feedbacklus.

4. Betekenis en Conclusie

• Systeembegrip: Het artikel biedt de eerste volledige eindige-elementensimulatie van de AMS-MSR-interactie, waarmee wordt bewezen dat het mu-metal velden vervormt maar de lineariteit behoudt die nodig is voor regelalgoritmen.
• Methodologische vooruitgang: Het demonstreert een succesvolle workflow voor het gebruik van Genetische Algoritmen om sensorplaatsing te optimaliseren in complexe magnetische omgevingen, en gaat voorbij aan trial-and-error inbedrijfstelling.
• Impact op n2EDM: De geoptimaliseerde AMS zorgt voor de magnetische stabiliteit die nodig is voor het n2EDM-experiment om zijn gevoeligheidsdoelen te bereiken, en minimaliseert systematische onzekerheden die gerelateerd zijn aan fluctuaties in het magnetische veld.
• Brede toepasbaarheid: Het ontwikkelde simulatie- en optimalisatiekader is overdraagbaar naar andere actieve magnetische afschermingssystemen in hoog-precisie natuurkunde-experimenten en medische instrumentatie.

Kortom, de auteurs hebben succesvol de kloof overbrugd tussen theoretisch ontwerp en experimentele realiteit voor de Actieve Magnetische Schild van n2EDM, waarbij geavanceerde simulatie werd gebruikt om het regelsysteem te verfijnen en de stabiliteit te garanderen die nodig is voor metingen van fundamentele natuurkunde van de volgende generatie.