DSWIM:Efficient and Stable Deterministic Computation of Warm Inflation Perturbations

Résumé technique : DSWIM – Calcul déterministe efficace et stable des perturbations de l'inflation chaude

Énoncé du problème
Les perturbations de l'inflation chaude (WI - Warm Inflation) sont pilotées par des fluctuations à la fois quantiques et thermiques, nécessitant la résolution d'équations différentielles stochastiques couplées. Bien que le code SWIM, disponible publiquement, implémente une approche stochastique robuste, il nécessite une moyenne d'ensemble sur de nombreuses réalisations pour atteindre une convergence statistique, ce qui le rend extrêmement coûteux en termes de calcul pour l'inférence de paramètres et les évaluations de vraisemblance.

Une alternative déterministe basée sur l'approche de Fokker-Planck, qui fait évoluer la matrice de corrélation des variables de perturbation via des équations différentielles ordinaires (ODE), offre une réduction significative du coût de calcul. Cependant, les auteurs identifient que l'évolution déterministe directe souffre d'un mauvais conditionnement numérique sévère dans de nombreux modèles de WI. Cette instabilité provient du fait que les variables de perturbation évoluent sur des échelles très différentes, provoquant des éléments de la matrice de corrélation couvrant de nombreux ordres de grandeur (par exemple, de 13 à 90 ordres dans les modèles testés). Cette hiérarchie entraîne des erreurs d'arrondi de virgule flottante, des distorsions artificielles dans le facteur de correction $G(Q)$, et une rupture finale de l'évolution numérique, particulièrement dans le régime de dissipation fort. De plus, les implémentations déterministes précédentes (par exemple, WI2easy) présentaient des écarts systématiques par rapport aux résultats stochastiques dans le régime de dissipation modéré en raison d'un traitement incorrect du bruit thermique corrélé.

Méthodologie
Pour remédier à ces problèmes, les auteurs introduisent DSWIM (Deterministic SWIM), un nouveau module au sein du cadre SWIM qui implémente un cadre de matrice d'échelle (scaling matrix framework) pour les équations d'évolution déterministe.

1. Transformation d'échelle : Les auteurs définissent un vecteur d'état mis à l'échelle $\tilde{\Phi} = S\Phi$, où $S$ est une matrice de mise à l'échelle diagonale composée de puissances du paramètre de Hubble $H$. La mise à l'échelle est physiquement motivée par le comportement de mise à l'échelle effectif des variables de perturbation pendant l'inflation :

   • La perturbation métrique $\psi$ (adimensionnelle) reste non mise à l'échelle.
   • Les fluctuations de l'inflaton $\delta\phi$ et leurs dérivées sont mises à l'échelle par $H^{-1}$.
   • La perturbation de la densité d'énergie du rayonnement $\delta\rho_r$ est mise à l'échelle par $H^{-2}$.
   • La perturbation de la quantité de mouvement du rayonnement $\delta q_r$ est mise à l'échelle par $H^{-1}$ en fonction de son apparition dans les équations de perturbation.
     Cette transformation normalise les variables, garantissant que les éléments de la matrice de corrélation évoluée restent numériquement comparables tout au long de l'intégration, préservant ainsi la plage dynamique et évitant la perte de précision.

2. Préservation des observables physiques : Les auteurs démontrent que cette transformation d'échelle laisse la perturbation de courbure comobile $R$ et, par conséquent, le spectre de puissance scalaire primordial $P_R$, exactement invariantes. La matrice de dérive $\tilde{A}$ et la matrice de diffusion $\tilde{D}$ transformées sont dérivées pour maintenir la cohérence physique du système.

3. Correction des corrélations de bruit : Le formalisme déterministe est mis à jour pour incorporer explicitement les corrélations hors-diagonales dans la matrice de diffusion $D$. Les auteurs montrent que les termes de bruit thermique dans les équations de l'inflaton et du rayonnement ne sont pas indépendants mais sont corrélés en raison de la conservation de l'énergie et de la quantité de mouvement. L'incorporation de ces termes hors-diagonaux non nuls résout les écarts systématiques précédemment observés entre les codes stochastiques et déterministes dans le régime de dissipation modérée.

Résultats clés

• Stabilité numérique : Le cadre déterministe mis à l'échelle supprime avec succès les artefacts numériques observés dans les implémentations non mises à l'échelle. Pour un modèle de potentiel exponentiel de type "runaway", la transformation d'échelle réduit la hiérarchie de la matrice de covariance d'environ 40 à 50 ordres de grandeur, étendant la plage d'évolution stable du solveur déterministe de près d'un ordre de grandeur en $Q$ par rapport aux méthodes non mises à l'échelle.
• Accord avec les méthodes stochastiques : Les résultats de DSWIM montrent un excellent accord avec l'implémentation stochastique de SWIM à travers les régimes de dissipation faible, intermédiaire et fort, validant la précision physique de l'approche déterministe mise à l'échelle.
• Efficacité de calcul : L'approche déterministe produit des accélérations substantielles. Pour des modèles représentatifs, DSWIM réduit le temps d'exécution pour calculer le facteur de correction $G(Q)$ d'environ deux ordres de grandeur par rapport au mode stochastique de SWIM et surpasse nettement le code déterministe existant WI2easy. Par exemple, le temps d'exécution pour un modèle de potentiel quartique est passé d'environ 116 secondes (stochastique) à environ 1,56 seconde (déterministe).
• Résolution des écarts : En implémentant correctement la structure du bruit thermique corrélé dans la matrice de diffusion, DSWIM élimine les déviations systématiques entre les résultats stochastiques et déterministes dans le régime de dissipation modérée ($0,01 \lesssim Q \lesssim 10$).

Signification et revendications
L'article établit DSWIM comme un cadre numériquement robuste et efficace pour calculer les perturbations scalaires de l'inflation chaude. Les auteurs affirment que DSWIM est actuellement le seul solveur de perturbations de WI disponible publiquement qui fournit un formalisme déterministe numériquement stabilisé directement interfacé avec le cadre d'analyse de vraisemblance Cobaya. Cette capacité permet des études d'inférence de paramètres efficaces pour les modèles d'inflation chaude qui étaient auparavant prohibitives en raison du coût de la moyenne stochastique ou de l'instabilité numérique des méthodes déterministes non mises à l'échelle. Ce travail ne propose pas de nouveaux modèles physiques, mais fournit un outil de calcul nécessaire pour explorer avec précision et efficacité l'espace des paramètres des scénarios d'inflation chaude existants.