We present SimpleBounce, a C++ package for finding the bounce solution for the false vacuum decay. This package is based on a flow equation which is proposed by the author and solves Coleman-Glaser-Martin's reduced problem: the minimization problem of the kinetic energy while fixing the potential energy. The bounce configuration is obtained by a scale transformation of the solution of this problem. For models with 1--8 scalar field(s), the bounce action can be calculated with O(0.1) % accuracy in O(0.1) s. This package is available at https://github.com/rsato64/SimpleBounce.
- Papier-ID: 1908.10868
- Titel: SimpleBounce: ein einfaches Paket für den Zerfall des falschen Vakuums
- Autor: Ryosuke Sato (Deutsches Elektronen-Synchrotron (DESY))
- Klassifizierung: hep-ph astro-ph.CO hep-th
- Veröffentlichungszeit: DESY 19-148 (2019)
- Papierlink: https://arxiv.org/abs/1908.10868
Dieses Papier stellt SimpleBounce vor, ein C++-Softwarepaket zur Berechnung von Bounce-Lösungen beim Zerfall des falschen Vakuums. Das Paket basiert auf den vom Autor vorgeschlagenen Flussgleichungen und löst das Coleman-Glaser-Martin-Vereinfachungsproblem: Minimierung der kinetischen Energie bei festgelegtem Potenzial. Die Bounce-Konfiguration wird durch die Skalentransformation der Lösung dieses Problems erhalten. Für Modelle mit 1-8 Skalarfeldern kann das Paket die Bounce-Wirkung in etwa 0,1 Sekunden mit einer Genauigkeit von O(0,1)% berechnen.
Der Zerfall des falschen Vakuums ist ein wichtiges Phänomen in der Teilchenphysik und Kosmologie. Die Lebensdauer des falschen Vakuums kann durch das euklidische Pfadintegral berechnet werden, wobei die Bounce-Lösung als Sattelpunkt der Wirkung den dominanten Beitrag zur Zerfallsbreite liefert.
Es existieren bereits mehrere numerische Softwarepakete zur Berechnung von Bounce-Lösungen, darunter:
- CosmoTransitions
- AnyBubble
- BubbleProfiler
Allerdings gibt es bei diesen Methoden noch Verbesserungspotenzial in Bezug auf Recheneffizienz und numerische Stabilität.
Der Autor schlägt einen neuen Rahmen basierend auf Gradientenflussgleichungen vor, der kürzlich von Chigusa, Moroi und Shoji vorgestellt wurde. Dieses Papier implementiert die vom Autor in Referenz 1 vorgeschlagenen spezifischen Flussgleichungen mit dem Ziel, ein schnelles, stabiles und benutzerfreundliches Werkzeug zur Berechnung von Bounce-Lösungen bereitzustellen.
- Entwicklung des SimpleBounce-Pakets: C++-Implementierung basierend auf Gradientenflussgleichungen zur Berechnung von Bounce-Lösungen beim Zerfall des falschen Vakuums
- Implementierung der numerischen Lösung des Coleman-Glaser-Martin-Vereinfachungsproblems: Erhalt von Bounce-Lösungen durch Minimierung der kinetischen Energie bei festgelegtem Potenzial
- Bereitstellung hocheffizienter Rechenleistung: Für Modelle mit 1-8 Skalarfeldern kann eine Genauigkeit von O(0,1)% in O(0,1) Sekunden erreicht werden
- Validierung der Methodengenauigkeit: Überprüfung der Ergebnisse durch Vergleich mit bestehenden Werkzeugen wie CosmoTransitions
Gegeben sei eine euklidische Wirkung mit nφ Skalarfeldern:
SE[ϕ]=T[ϕ]+V[ϕ]
wobei:
T[ϕ]=∑i=1nϕ∫ddx21(∇ϕi)2V[ϕ]=∫ddx[V(ϕ)−V(ϕFV)]
Das Ziel ist es, die Bounce-Lösung zu finden, die die Bewegungsgleichungen und Randbedingungen erfüllt.
Anstatt die Bewegungsgleichungen direkt zu lösen:
−∇2ϕi+∂ϕi∂V=0,lim∣x∣→∞ϕi(x)=ϕFV,i
wird stattdessen das Minimierungsproblem der kinetischen Energie Tφ gelöst, während das negative Potenzial Vφ < 0 festgelegt wird.
Es wird eine Funktion φᵢ(r,τ) eingeführt, deren Fluss durch die folgende Gleichung beschrieben wird:
∂τ∂ϕi(r,τ)=∇2ϕi−λ[ϕ]∂ϕi∂V(ϕ)
wobei:
λ[ϕ]=∑i=1nϕ∫0∞drrd−1(∂ϕi∂V(ϕ))2∑i=1nϕ∫0∞drrd−1∂ϕi∂V(ϕ)∇2ϕi
Die Bounce-Lösung wird durch die Skalentransformation des Fixpunktes erhalten:
ϕB(r)=limτ→∞ϕ(λ−1/2r,τ)
- Numerische Stabilität: Da das Potenzial V festgelegt ist, ist der Fixpunkt kein Sattelpunkt, was die numerische Stabilität der Berechnung gewährleistet
- Adaptive Algorithmen: Durch Anpassung der Anfangsparameter r₀ und σ wird sichergestellt, dass die Randbedingungen erfüllt sind
- Effiziente Diskretisierung: Verwendung eines speziellen Laplace-Operator-Diskretisierungsschemas, das für kugelsymmetrische Probleme geeignet ist
Das Papier verwendet mehrere Benchmark-Modelle für Tests:
- Einfeldiges Skalarmodell: V(φ) = ½φ² - ⅓φ³
- Acht Benchmark-Modelle aus Referenz 10
- Vier Benchmark-Modelle aus Referenz 24
- Euklidische Wirkung SE: Euklidische Wirkung der Bounce-Lösung
- Rechenzeit: Laufzeit des Algorithmus
- Genauigkeit: Relativer Fehler zu Referenzergebnissen
Hauptsächlich Vergleich mit CosmoTransitions, dem weit verbreiteten Standardwerkzeug in diesem Bereich.
- Raumdimension: d = 3 oder 4
- Gitterpunkte: n = 100-800
- Zeitschrittweite: δτ ~ O(δr²) zur Gewährleistung numerischer Stabilität
- Randbedingungen: φᵢ,ₙ = φFV,i
Tabelle 1: Vergleich mit CosmoTransitions
| Modell | SE (SimpleBounce) | SE (CosmoTransitions) | Zeit(SB) | Zeit(CT) |
|---|
| #1 | 52,4 | 52,6 | 0,04s | 0,05s |
| #2 | 20,8 | 21,1 | 0,04s | 0,35s |
| #8 | 45,9 | 46,1 | 0,16s | 0,23s |
Die Ergebnisse zeigen, dass die Berechnungsergebnisse von SimpleBounce mit CosmoTransitions hochgradig übereinstimmen und in den meisten Fällen schneller sind.
Abbildung 2 zeigt den Einfluss verschiedener Gitterpunkte n und Flusszeiten τ₁ auf die Ergebnisse:
- SE-Werte konvergieren mit zunehmenden n und τ₁
- Mit n=100 kann eine Genauigkeit von O(0,1)% in O(0,1) Sekunden erreicht werden
- Die Laufzeit ist proportional zu n³
- Stabilität gegenüber der Wahl der Anfangskonfiguration
- Dünnwand-Bounce erfordert mehr Rechenzeit als Dickwand-Bounce
- Die Rechenzeit ist ungefähr proportional zur Anzahl der Skalarfelder
- CosmoTransitions: Basierend auf Finite-Differenzen-Methoden zur Lösung von Randwertproblemen
- AnyBubble: Hocheffizienter numerischer Löser für Mehrfeld-Fälle
- BubbleProfiler: Blasenwand-Berechnung mit halbanalytischen Techniken
- Coleman-Glaser-Martin-Theorie: Verbindung zwischen Bounce-Lösung und Minimierungsproblem
- Gradientenfluss-Methode: Neuer Rahmen von Chigusa et al.
- O(d)-Symmetrie: Kugelsymmetrische Eigenschaften der Bounce-Lösung
- SimpleBounce bietet eine schnelle und genaue Methode zur Berechnung von Bounce-Lösungen
- Die auf Gradientenfluss basierende Methode weist gute numerische Stabilität auf
- Das Softwarepaket ist recheneffizient gegenüber bestehenden Werkzeugen
- Dünnwand-Grenzfall: Für Dünnwand-Bounce nimmt die Rechenzeit erheblich zu
- Gitterabhängigkeit: Bei hohen Genauigkeitsanforderungen ist ein feineres Gitter erforderlich, was die Rechenkosten erhöht
- Beschränkung auf Skalarfelder: Bisherige Tests konzentrieren sich hauptsächlich auf 1-8 Skalarfelder
- Optimierung der Recheneffizienz für Dünnwand-Bounce
- Erweiterung auf mehr Skalarfelder
- Integration in größere physikalische Rechnerrahmen
- Methodische Innovation: Der auf Gradientenfluss basierende Ansatz bietet neue Perspektiven für numerische Lösungen
- Hoher praktischer Wert: Das Open-Source-Softwarepaket ermöglicht die Verwendung und Überprüfung durch die Gemeinschaft
- Ausgezeichnete Leistung: Signifikante Verbesserung der Recheneffizienz bei Gewährleistung der Genauigkeit
- Solide theoretische Grundlagen: Aufbau auf der Coleman-Glaser-Martin-Theorie
- Begrenzte Testabdeckung: Tests konzentrieren sich hauptsächlich auf relativ einfache Modelle
- Dünnwand-Leistung: Die Behandlung von Dünnwand-Bounce erfordert noch Verbesserungen in der Effizienz
- Dokumentationsvollständigkeit: Als Softwarepaket könnten Benutzerdokumente detaillierter sein
- Akademischer Beitrag: Bietet neue Werkzeuge für die Berechnung des Zerfalls des falschen Vakuums
- Praktischer Wert: Könnte zum Standardberechnungswerkzeug in diesem Bereich werden
- Reproduzierbarkeit: Open-Source-Code gewährleistet die Wiederholbarkeit der Ergebnisse
- Vakuumstabilitätsanalyse in der Teilchenphysik
- Forschung zu kosmologischen Phasenübergängen
- Numerische Berechnungen in der Skalarfeldtheorie
- Theoretische Forschung, die schnelle Bounce-Lösungsberechnungen erfordert
Die Hauptreferenzen umfassen klassische theoretische Arbeiten von Coleman et al., Papiere zu bestehenden numerischen Methoden sowie frühere theoretische Forschungen des Autors zu Gradientenflussgleichungen. Die vollständige Implementierung des Softwarepakets ist auf GitHub verfügbar.