An efficient one-loop EFTofLSS framework for Vainshtein-screened Horndeski gravity

Technische Samenvatting: Een efficiënt one-loop EFTofLSS-framework voor Vainshtein-afgeschermde Horndeski-zwaartekracht

Probleemstelling
Huidige en toekomstige Stage IV-galactische surveys (bijv. Euclid, DESI) breiden hun observationele bereik uit naar de quasi-lineaire en niet-lineaire regimes van de grootschalige structuur van het Universum (LSS). Om het beperkende vermogen van deze surveys volledig te benutten, zijn theoretische voorspellingen nodig die verder gaan dan de lineaire perturbatietheorie, specifiek op schalen $k \lesssim 0.3 \, h \, \text{Mpc}^{-1}$. Hoewel de Effective Field Theory of Large-Scale Structure (EFTofLSS) een robuust framework biedt voor het modelleren van deze schalen in $\Lambda$CDM, vormt het uitbreiden hiervan naar gemodificeerde zwaartekrachttheorieën (MG) een uitdaging. Specifiek vereist Horndeski-zwaartekracht — een algemene klasse van scalaar-tensor theorieën die de meeste dynamische donkere energie-modellen omvat — een nauwkeurige behandeling van niet-lineaire dynamica, Vainshtein-afscherming en tijdafhankelijke achtergrondevolutie. Bestaande implementaties vertrouwen vaak op de Einstein–de Sitter (EdS) benadering voor perturbatiekernen of computationeel dure Green's functie-methoden die moeite hebben met de schaal-onafhankelijke groei en de complexe tijdsafhankelijke evolutie die kenmerkend zijn voor MG-modellen.

Methodologie
De auteurs presenteren een uitbreiding van de publieke code PyBird, ontworpen voor het berekenen van one-loop EFTofLSS-voorspellingen voor gebiaste tracers in roodverschuivingsruimte. De methodologie omvat drie primaire technische vooruitgangen:

1. Framework Uitbreiding: De pipeline wordt uitgebreid om quasi-statische, Vainshtein-afgeschermde, luminale Horndeski-modellen te ondersteunen (waarbij de snelheid van zwaartekrachtgolven $c_T^2 = 1$) in zowel het covariante formalisme als de taal van de Effective Field Theory of Dark Energy (EFTofDE). Het ondersteunt ook het nDGP-model. De implementatie maakt gebruik van de $\alpha$-basis ($\alpha_B, \alpha_M, \alpha_K$) en koppelt met H–EFTCAMB om achtergrondevolutie en lineaire perturbatiehoeveelheden aan te leveren.
2. Directe ODE-benadering: De auteurs vervangen de standaard Green's functie-aanpak voor het berekenen van tijdsafhankelijke functies in perturbatiekernen door een directe Gewone Differentiaalvergelijking (ODE) methode. In plaats van Green's functies te convolueren met brontermen, lost de code scalaire ODE's iteratief op om de tijdsafhankelijke coëfficiënten ($G, U, V$) direct te verkrijgen. Deze aanpak vermijdt de noodzaak van een expliciete decompositie in groeiende en afnemende modi, wat moeilijk is in scenario's met gemodificeerde zwaartekracht waar de definitie van de groeiende modus niet triviaal is.
3. Exacte Tijdsevolutie: In tegen tegenstelling tot standaard EFTofLSS-analyses die vaak de EdS-benadering gebruiken voor kernen, berekent dit framework kernen met exacte tijdsafhankelijkheid. De auteurs valideren dat dit met een bescheiden computationele kost kan worden bereikt, waardoor de noodzaak verdwijnt om nauwkeurigheid in te ruilen voor efficiëntie.

De pipeline berekent de one-loop galaxy power spectrum multipolen (monopool en kwadrupool) inclusief counterterms en stochastische bijdragen. De auteurs valideren de numerieke consistentie van hun ODE-methode tegenover de standaard Green's functie-aanpak voor $\Lambda$CDM, $w$CDM en Quintessence, en tegenover onafhankelijke implementaties voor nDGP en Cubic Galileon modellen. Verder wordt de one-loop materie power spectrum gevalideerd tegen N-body simulaties (PySCo-EFT en ECOSMOG-EFT) voor de $\alpha_i \propto \Omega_{DE}$ parametrisatie.

Belangrijkste Bijdragen

• Computationele Efficiëntie: De directe ODE-methode verbetert de computationele efficiëntie aanzienlijk, vooral bij het evalueren van meerdere roodverschuivingen. Terwijl de Green's functie-methode voor elke roodverschuiving de convolutie-integralen opnieuw moet berekenen, lost de ODE-methode de groeivergelijking één keer op en hergebruikt de oplossing. Dit resulteert in een versnelling van ongeveer 20–25% voor volledige one-loop spectrumberekeningen en schaalt lineair met het aantal doel-roodverschuivingen.
• Gegeneraliseerde MG Ondersteuning: De code biedt een modulaire interface voor elk schaal-onafhankelijk, Vainshtein-afgeschermd luminaal Horndeski-model, waardoor de kloof tussen fenomenologische parametrisaties en covariante theorieën wordt overbrugd.
• Validatie: Het framework is rigoureus getest tegen N-body simulaties en vertoont overeenstemming op sub-percent niveau voor de boost factor $R(k) = P_{\text{EFT}}/P_{\Lambda\text{CDM}}$ in het bereik $0.1 \lesssim k \lesssim 0.3 \, h \, \text{Mpc}^{-1}$.
• Exacte versus Benaderde Kernen: De auteurs tonen aan dat de EdS-benadering, hoewel geldig voor $\Lambda$CDM, in extreme gemodificeerde zwaartekrachtscenario's op kleine schalen kan afwijken van de exacte tijdsevolutie met meerdere standaarddeviaties ($\sim 4\sigma$). Ze laten zien dat het behouden van exacte tijdsafhankelijkheid computationeel haalbaar en noodzakelijk is voor precisie-kosmologie in MG-contexten.

Resultaten
De auteurs passen het framework toe om twee fenomenologische parametrisaties van de EFTofDE amplitudes te beperken: $\alpha_i \propto a^3$ en $\alpha_i \propto \Omega_{DE}$, met behulp van Planck CMB, BOSS full-shape (FS) en DESI DR2 BAO data.

• Beperkingen: De inclusie van BOSS full-shape data versterkt de beperkingen op de gemodificeerde zwaartekrachtparameters ($c_B, c_M$) aanzienlijk vergeleken met CMB-only analyses. Voor de $\alpha_i \propto \Omega_{DE}$ casus levert de gecombineerde data $c_B = -0.41^{+0.25}_{-0.19}$ en $c_M = 0.32^{+0.61}_{-0.70}$ (68% CL) op, wat consistent is met Algemene Relativiteit binnen het 95% geloofwaardigheidsinterval.
• Projectie-effecten: De analyse vindt dat projectie-effecten (biases in MAP-schattingen door nuisance parameter degeneraties) geen significante zorg vormen wanneer sterke CMB-priors op $A_s$ en $n_s$ worden opgenomen.
• Model Showcase: De pipeline genereert succesvol one-loop power spectrum voorspellingen voor covariante modellen zoals de Cubic Galileon en nDGP, wat de toepasbaarheid voorbij eenvoudige amplitude-parametrisaties aantoont.

Betekenis
Dit artikel vestigt een systematisch, flexibel en computationeel efficiënt framework voor het uitbreiden van EFTofLSS-analyses buiten $\Lambda$CDM. Door het combineren van exacte tijdsafhankende perturbatietheorie met een praktische implementatie in PyBird, maakt dit werk robuuste one-loop voorspellingen mogelijk voor een brede klasse van gemodificeerde zwaartekrachttheorieën. De auteurs stellen dat het vermogen om exacte tijdsafhankelijkheid te behouden tegen een bescheiden computationele kosten is cruciaal voor toekomstige precisie-surveys (zoals DESI en Euclid), waarbij de impact van benaderingen zoals EdS model-afhankelijk en significant kan worden. Dit framework plaveit de weg voor precisietesten van zwaartekracht op kosmologische schalen met huidige en toekomstige LSS-observaties.