Non-linear evolution in $f(R)$ gravity: perturbative modelling of the Chameleon mechanism

Titel: Niet-lineaire evolutie in f(R)-zwaartekracht: perturbatieve modellering van het Chameleon-mechanisme

1. Het Probleem

De standaard $\Lambda$CDM-cosmologie staat voor theoretische en observationele uitdagingen, zoals de oorsprong van donkere energie en de Hubble- en $\sigma_8$-spanningen. Alternatieven omvatten modificaties van de algemene relativiteitstheorie (GR), zoals $f(R)$-zwaartekracht. In deze theorieën wordt de Ricci-scalar $R$ in de Einstein-Hilbert-actie vervangen door een functie $f(R)$.

Een cruciaal probleem bij het testen van deze theorieën is dat ze consistent moeten zijn met lokale tests van GR (zoals in het zonnestelsel) terwijl ze afwijkingen op kosmologische schalen toelaten. Dit wordt bereikt via het Chameleon-mechanisme, waarbij de massa van het extra scalair veld (en dus het bereik van de "vijfde kracht") afhankelijk is van de lokale matterdichtheid. In hoge dichtheden wordt het bereik klein (GR wordt hersteld), terwijl het in lage dichtheden groot wordt.

Het modelleren van de groei van structuren in $f(R)$-zwaartekracht is complex omdat:

1. De vergelijkingen voor de gravitatiepotentialen $\Phi$ en $\Psi$ niet langer lineair zijn als het Chameleon-mechanisme wordt meegenomen.
2. De vergelijking voor het scalair veld $\chi = \Phi - \Psi$ een niet-lineaire differentiaalvergelijking wordt die niet opgelost kan worden met simpele Fourier-transformaties.
3. Bestaande numerieke methoden (zoals N-body simulaties met multi-grid methoden) computatieel zeer duur zijn.

2. Methodologie

De auteurs ontwikkelen een perturbatieve oplossing voor de niet-lineaire vergelijking van het scalair veld $\chi$, binnen het kader van een hybride Euleriaans-Lagrangeaans schema.

• Het Model: Ze gebruiken het Hu-Sawicki $f(R)$-model, gekozen zodat het de achtergrond-uitdijing van $\Lambda$CDM nabootst.
• De Vergelijkingen: Ze beschouwen sferische perturbaties. De kernvergelijking voor $\chi$ is:
  $$\nabla^2_x \chi - \frac{1}{\bar{x}_C^2(x,a)}\chi = -H^2 a^2 \Omega_m \delta$$
  waarbij $\bar{x}_C$ de gereduceerde Compton-golflengte is, die afhangt van de lokale dichtheid (via $R$).
• Perturbatieve Aanpak: In plaats van de volledige niet-lineaire vergelijking numeriek op te lossen met iteratieve solvers (zoals in eerdere werken), ontwikkelen ze een reeksontwikkeling:
  $$\chi = \sum_{n=1}^{\infty} \epsilon^n \chi^{(n)}$$
  Hierbij wordt $\chi^{(1)}$ opgelost als een lineaire vergelijking (zonder Chameleon), en worden hogere-orde correcties ($\chi^{(2)}, \chi^{(3)}, \dots$) systematisch berekend op basis van de dichtheidsvariaties. Ze stoppen bij de vierde orde.
• Evolutie-Schema:
  • Stap I (Euleriaans): Oplossing van de ruimtelijke vergelijkingen voor de potentialen $\Phi_+$ en $\chi$ voor een vaste dichtheid.
  • Stap II (Lagrangeaans): Evolutie van dichtheid en snelheid door het volgen van schalen (shells) in de tijd.
• Testgevallen: Ze simuleren gladde, gecompenseerde "top-hat" dichtheidsprofielen met drie verschillende schalen ($x_{top} = 0.002, 2, 200 \, h^{-1}\text{Mpc}$) om drie regimes te testen.

3. Belangrijkste Bijdragen

1. Efficiënte Numerieke Methode: De auteurs presenteren een perturbatieve methode die analytische oplossingen gebruikt voor elke orde van de reeks. Dit is computatieel veel goedkoper dan traditionele multi-grid of relaxatie-methoden voor niet-lineaire ODE's, terwijl het toch de Chameleon-screening nauwkeurig modelleert.
2. Ruimtelijk-Tijds Evolutie van de Compton-schaal: Ze tonen expliciet aan hoe de Compton-schaal ($\bar{x}_C$) dynamisch verandert binnen een enkele perturbatie, zowel in de tijd als in de ruimte, en hoe dit terugkoppelt op de dichtheidsontwikkeling.
3. Nieuw Ontdekt Fenomeen: Ze identificeren en verklaren een dichtheidsverhoging nabij het centrum van de perturbatie, een kenmerk dat niet eerder in de literatuur is gerapporteerd voor sferische ineenstorting in $f(R)$-theorieën.
4. Niet-unieke DVDR: Ze tonen aan dat de relatie tussen dichtheid en snelheidsdivergentie (DVDR) niet langer één-op-één is in het Chameleon-regime, zelfs niet voor één enkel profiel.

4. Resultaten

De resultaten worden geanalyseerd aan de hand van de parameter $Q = \bar{x}_C / x_{top}$ (verhouding Compton-schaal tot perturbatieschaal):

• Sterk-veld regime ($Q \gg 1$): De perturbatieschaal is veel kleiner dan de Compton-schaal. Het Chameleon-mechanisme heeft geen effect; de groei is versneld ten opzichte van GR (factor 4/3 in de effectieve Newton-constante), maar het profiel behoudt zijn vorm.

• Zwak-veld regime ($Q \ll 1$): De perturbatieschaal is veel groter dan de Compton-schaal. De modificaties zijn onderdrukt en de evolutie lijkt sterk op GR, behalve aan de randen waar dichtheidsgradiënten steil zijn.

• Intermediair regime ($Q \sim 1$): Dit is het meest interessante regime (bijv. $x_{top} = 2 \, h^{-1}\text{Mpc}$).

  • Rand-effect: Er treedt een bekende dichtheidsverhoging op aan de buitenrand van de "top-hat" door de accumulatie van schalen.
  • Centrale verhoging (Nieuw): De auteurs vinden een subtiele, maar significante verhoging van de dichtheid dicht bij de oorsprong.
    • Mechanisme: De dichtheidsgradiënt binnenin veroorzaakt een gradiënt in de Compton-schaal. Lagere dichtheid bij de oorsprong leidt tot een iets grotere Compton-schaal, wat op zijn beurt de kracht op in vallende schalen verhoogt, waardoor ze versnellen en een tweede piek vormen bij de oorsprong.
  • Onderdrukking van groei: In dit regime wordt de totale groei van de dichtheid onderdrukt ten opzichte van zowel GR als $f(R)$ zonder Chameleon-screening.

• Fase-ruimte Evolutie (DVDR):
  In standaard GR en $f(R)$ zonder screening is er een unieke, eenduidige relatie tussen dichtheid ($\Delta$) en snelheidsdivergentie ($\delta v$). In het Chameleon-regime (voor het intermediaire profiel) wordt deze relatie niet-eenduidig (multi-valued). Een bepaalde dichtheid kan corresponderen met meerdere waarden van snelheidsdivergentie, afhankelijk van de locatie binnen de perturbatie. Dit wordt vergeleken met een hysterese-lus, veroorzaakt door de niet-lineaire terugkoppeling tussen dichtheid en de Compton-schaal.

5. Significatie

Dit onderzoek biedt een krachtig en computatieel efficiënt alternatief voor dure N-body simulaties om de niet-lineaire effecten van het Chameleon-mechanisme te bestuderen. De ontdekking van de centrale dichtheidsverhoging en de niet-eenduidige DVDR-relatie heeft belangrijke implicaties voor:

• Het interpreteren van toekomstige kosmologische surveys (zoals Euclid en DESI) die de groei van structuren meten.
• Het begrijpen van de vorming van halos en de verdeling van materie in $f(R)$-modellen.
• Het identificeren van nieuwe "signatures" (zoals hysterese-effecten) die mogelijk kunnen worden gebruikt om $f(R)$-zwaartekracht te onderscheiden van standaard GR en donkere energie-modellen.

De auteurs concluderen dat het Chameleon-mechanisme een rijke structuur introduceert in de dynamiek van kosmologische perturbaties die verder gaat dan de eenvoudige schaalafhankelijke groei die vaak wordt aangenomen.