Annexe B. Calcul de la diffusion spatiale à basse température par le théorème de régression

L'équation (86) donne un équivalent à basse température du coefficient de diffusion spatiale de la quasi-particle $\gamma$ dans un gaz de phonons. Dans le cas où la relation de dispersion $\epsilon_k$ de la quasi-particle atteint son minimum en un nombre d'onde $k_0 > 0$, la justification qu'en fournit le texte principal, à l'aide d'un modèle stochastique de Langevin minimal, n'est pas pleinement convaincante. Nous présentons ici une démonstration basée sur le théorème de régression et sur un développement systématique à basse température de la solution de l'équation de Fokker-Planck correspondante. Nous justifions aussi l'expression (86) dans la zone de transition entre le cas $k_0 \equiv$ et le cas $k_0 > 0$, qui doit faire l'objet d'un traitement à part.

**Théorème de régression.** La connaissance du coefficient de diffusion spatiale passe par celle de la fonction de corrélation $C_{ij}(\tau)$ des coordonnées du vecteur vitesse dans l'état stationnaire, voir les équations (80,81,82,83). Par isotropie, nous nous limitons au calcul de $C_{zz}(\tau)$, $\tau \ge 0$. Dans le cadre de Fokker-Planck, nous disposons alors du théorème de régression suivant :

$$C_{zz}(\tau) = \int \frac{d^3k}{(2\pi)^3} \frac{1}{\hbar} \frac{d\epsilon_k}{dk} \hat{k}_z \Pi(\mathbf{k}, \tau) \quad (B.1)$$

où la distribution $\Pi(\mathbf{k}, \tau)$ a évolué pendant $\tau$ selon l'équation de Fokker-Planck (50) à partir de la condition initiale

$$\Pi(\mathbf{k}, 0) = \frac{1}{\hbar} \frac{d\epsilon_k}{dk} \hat{k}_z \Pi_0(k) \quad (B.2)$$

avec $\Pi_0(k)$ la solution stationnaire normalisée, invariante par rotation, de cette même équation (50), et $\hat{k}_z = k_z/k$. Comme $\Pi(\mathbf{k}, \tau)$ est de moment cinétique un, porté par l'axe $Oz$, nous pouvons poser $\Pi(\mathbf{k}, \tau) = \Pi_1(k, \tau)\hat{k}_z$ et obtenir l'équation du mouvement

$$\frac{\partial}{\partial \tau} \Pi_1(k, \tau) = -\frac{2F(k)\Pi_1(k, \tau)}{\hbar k} - \frac{1}{\hbar} \frac{\partial}{\partial k} [F(k)\Pi_1(k, \tau)] + \frac{1}{\hbar^2} \frac{\partial^2}{\partial k^2} [D_{\parallel}(k)\Pi_1(k, \tau)] \\ + \frac{1}{\hbar^2 k} \frac{\partial}{\partial k} \{[4D_{\parallel}(k) - 2D_{\perp}(k)]\Pi_1(k, \tau)\} + \frac{1}{\hbar^2 k^2} [2D_{\parallel}(k) - 4D_{\perp}(k)]\Pi_1(k, \tau) \quad (B.3)$$

L'équation (72) permettrait d'exprimer $\Pi_0(k)$ en termes de la force moyenne $F(k)$ et des coefficients de diffusion en impulsion $D_{\parallel,\perp}(k)$. Nous trouvons plus éclairant de l'écrire sous la forme $\Pi_0(k) = \exp(-\beta\tilde{\epsilon}_k)/Z$, où $Z$ est un facteur de normalisation, et d'éliminer au contraire la force $F(k)$ au profit de la relation de dispersion effective $\tilde{\epsilon}_k$. Comme nous l'apprend d'ailleurs la comparaison des équations (72,73), celle-ci est soumise à la contrainte $\tilde{\epsilon}_k - \epsilon_k = O(T^2)$ lorsque $T \to 0$ avec $k - k_0 = O(T^{1/2})$, donc on peut, comme la vraie relation de dispersion $\epsilon_k$, la développer autour de la position $\tilde{k}_0$ de son minimum (par convention égal à $\Delta_*$), nulle si $k_0 \equiv 0$, dépendant de la température mais très proche de $k_0$ si $k_0 > 0$, avec des coefficients dépendant eux aussi faiblement de la température :$^{62}$

$$\tilde{\epsilon}_k = \Delta_* + \frac{\hbar^2(k - \tilde{k}_0)^2}{2\tilde{m}_*} + \frac{\hbar^2(k - \tilde{k}_0)^3\tilde{b}}{3\tilde{m}_*} + \frac{\hbar^2(k - \tilde{k}_0)^4\tilde{\ell}^2}{4\tilde{m}_*} + O(k - \tilde{k}_0)^5 \quad (B.4)$$

En d'autres termes, il ne faut pas croire aux prédictions tirées de l'équation de Fokker-Planck (une approximation de l'équation maîtresse (47) dans la limite (3)) si elles dépendent du faible écart entre les masses effectives $m_*$ et $\tilde{m}_*$, entre les longueurs $b$ et $\tilde{b}$, entre les nombres d'onde minimiseurs $k_0$ et $\tilde{k}_0$, ou de la longueur $\tilde{\ell}$ ou des coefficients d'ordre supérieur dans le développement (B.4). Effectuons une dernière transformation : l'opérateur sur les fonctions de $k$ défini par le second membre de l'équation (B.3) n'est pas autoadjoint, mais peut l'être rendu par le changement de fonction $\Pi_1(k) = k^{-1} \exp(-\beta\tilde{\epsilon}_k/2)\psi(k)$. Après intégration temporelle formelle, nous aboutissons pour $\tau \ge 0$ à la belle écriture (le facteur 3 au dénominateur vient de l'intégration angulaire) :

$$C_{zz}(\tau) = \frac{\langle \psi_S | \exp(-\mathcal{H}_1 \tau) | \psi_S \rangle}{3\hbar^2 \int_0^{+\infty} dk \, k^2 e^{-\beta\tilde{\epsilon}_k}} = \sum_{n \in \mathbb{N}} \frac{|\langle \psi_S | \phi_n \rangle|^2 e^{-\omega_n t}}{3\hbar^2 \int_0^{+\infty} dk \, k^2 e^{-\beta\tilde{\epsilon}_k}} \quad \text{avec} \quad \left\{ \begin{array}{l} \mathcal{H}_1 = -\frac{d}{dk} \frac{D_{\parallel}(k)}{\hbar^2} \frac{d}{dk} + U_1(k) \\ \langle k | \psi_S \rangle = k \frac{d\epsilon_k}{dk} e^{-\beta\tilde{\epsilon}_k/2} \end{array} \right. \quad (B.5)$$

<footnote n="62">62. Si l'on admet que l'amplitude de diffusion $\phi - \gamma$ (31) possède un développement en puissances entières de $q$ à $\mathbf{k}$ fixé, alors $F(k)$, $D_{\parallel,\perp}(k)$ et $\tilde{\epsilon}_k$ possèdent un développement en puissances entières de $T$ à $k$ fixé, et l'on a $\tilde{k}_0/k_0 = 1 + O(T^2)$, $\tilde{m}_*/m_* = 1 + O(T)$, $\tilde{b}/b = 1 + O(T)$, $\tilde{\ell}/\ell = O(1)$.</footnote>
Nous avons introduit la décomposition spectrale du hamiltonien fictif $\mathcal{H}_1$ autoadjoint pour le produit scalaire unidimensionnel $\langle\phi_a|\phi_b\rangle = \int_0^{+\infty} dk \phi_a^*(k)\phi_b(k)$, l'état propre $|\phi_n\rangle$ de valeur propre $\omega_n$ étant normalisé. Le potentiel vaut

$$U_1(k) = \frac{1}{\hbar^2} \left[ \frac{1}{4} D_{\parallel}(k) \left( \beta \frac{d\tilde{\epsilon}_k}{dk} \right)^2 - \frac{1}{2} D_{\parallel}(k) \beta \frac{d^2\tilde{\epsilon}_k}{dk^2} - \frac{1}{2} \frac{dD_{\parallel}(k)}{dk} \beta \frac{d\tilde{\epsilon}_k}{dk} - \frac{D_{\parallel}(k)}{k} \beta \frac{d\tilde{\epsilon}_k}{dk} + \frac{1}{k} \frac{dD_{\parallel}(k)}{dk} + \frac{2D_{\perp}(k)}{k^2} \right] \quad (B.6)$$

Dans le cas $k_0 > 0$. Pour prendre la limite de basse température dans le cas $k_0 > 0$, nous effectuons un changement d'échelle adapté sur le vecteur d'onde et le temps en posant $k - \tilde{k}_0 = xk_*$, $\tau = \bar{\tau}/\bar{\Gamma}$, et donc $\mathcal{H}_1 = \tilde{\mathcal{H}}_1\bar{\Gamma}$, $\omega_n = \bar{\omega}_n\bar{\Gamma}$ et $\phi_n(k) = k_*^{-1/2}\tilde{\phi}_n(x)$, avec le nombre d'onde thermique $k_* = (\tilde{m}_*k_B T)^{1/2}/\hbar$ et le taux $\bar{\Gamma} = D_{\parallel}(\tilde{k}_0)/(\tilde{m}_*k_B T)$ construit sur le modèle du taux $\Gamma$ de l'équation (88). Il reste à développer $\tilde{\mathcal{H}}_1$ en puissances de $k_*$ à x fixé (x décrit $\mathbb{R}$ dans cette limite), puis à traiter les ordres successifs par la théorie des perturbations :

$$\tilde{\mathcal{H}}_1 = - \frac{d}{dx} \frac{D_{\parallel}(\tilde{k}_0 + xk_*)}{D_{\parallel}(\tilde{k}_0)} \frac{d}{dx} + \frac{\hbar^2 k_*^2}{D_{\parallel}(\tilde{k}_0)} U(\tilde{k}_0 + xk_*) \stackrel{x \text{ fixé}}{=} \tilde{\mathcal{H}}_1^{(0)} + \tilde{\mathcal{H}}_1^{(1)} + \tilde{\mathcal{H}}_1^{(2)} + \dots \quad (B.7)$$

On aura besoin aussi de développer la fonction d'onde source à l'ordre sous-dominant et le dénominateur de $C_{zz}(\tau)$ dans l'équation (B.5) à l'ordre dominant :

$$\langle k | \psi_S \rangle \stackrel{x \text{ fixé}}{=} \frac{\hbar^2 \tilde{k}_0 k_*}{\tilde{m}_*} e^{-\beta \Delta_*/2} x e^{-x^2/4} \left[ 1 + x \left( \frac{k_*}{\tilde{k}_0} + k_* \tilde{b} \right) - \frac{1}{6} k_* \tilde{b} x^3 + O(k_*^2) \right], \quad \int_0^{+\infty} dk \, k^2 e^{-\beta \tilde{\epsilon}_k} \stackrel{k_* \to 0}{\sim} \tilde{k}_0^2 k_* e^{-\beta \Delta_*} \quad \int_{-\infty}^{+\infty} dx \, x e^{-x^2/2} \quad (B.8)$$

À l'ordre zéro en $k_*$, on tombe sur un hamiltonien d'oscillateur harmonique unidimensionnel réduit (masse égale à 1/2, constante de Planck réduite égale à 1, pulsation d'oscillation égale à 1) et de spectre $\tilde{\omega}_n^{(0)} = n$ :

$$\tilde{\mathcal{H}}_1^{(0)} = - \frac{d^2}{dx^2} + \frac{1}{4} x^2 - \frac{1}{2} = \hat{a}^\dagger \hat{a} \quad \text{avec} \quad \hat{a} = \frac{1}{2} x + \frac{d}{dx}, \quad \hat{a}^\dagger = \frac{1}{2} x - \frac{d}{dx}, \quad [\hat{a}, \hat{a}^\dagger] = 1 \quad (B.9)$$

À l'ordre dominant, la fonction d'onde source est donc proportionnelle à la fonction d'onde $\tilde{\phi}_1^{(0)}(x)$ du premier état excité de l'oscillateur harmonique. La contribution des états excités $n > 0$ dans l'équation (B.5) est ainsi dominée par $n = 1$ et l'on garde

$$C_{zz}^{n \neq 0}(\tau) \stackrel{\bar{\tau} \text{ fixé}}{\sim} \frac{k_B T}{3 \tilde{m}_*} e^{-\bar{\tau}} \quad (B.10)$$

On justifie ainsi le premier terme entre crochets dans l'équation (85), compte tenu de la note 62. La contribution de $n = 0$ est plus délicate à obtenir. En effet, le recouvrement de $\psi_S$ avec l'état fondamental $\tilde{\phi}_0$ est nul à l'ordre de $\tilde{\mathcal{H}}_1^{(0)}$, comme la pulsation propre $\omega_0$ de l'état fondamental, ce qui conduit à une forme indéterminée 0/0 après intégration de $C_{zz}(\tau)$ sur le temps comme dans l'équation (83). Il faut aller à l'ordre un en $k_*$ pour obtenir un recouvrement non nul,

$$\begin{align*}
\tilde{\mathcal{H}}_1^{(1)} &= k_* \frac{\frac{dD_{\parallel}}{dk}(\tilde{k}_0)}{D_{\parallel}(\tilde{k}_0)} \left( -\frac{d}{dx}x \frac{d}{dx} + \frac{1}{4}x^3 - x \right) + k_* \tilde{b} \left( \frac{1}{2}x^3 - x \right) - \frac{k_*}{\tilde{k}_0} x = \left( \frac{k_* \frac{dD_{\parallel}}{dk}(\tilde{k}_0)}{D_{\parallel}(\tilde{k}_0)} + \frac{3}{2} k_* \tilde{b} \right) (\hat{a}^{\dagger 2} \hat{a} + \hat{a}^{\dagger} \hat{a}^2) \\
&\quad + \frac{1}{2} k_* \tilde{b} (\hat{a}^{\dagger 3} + \hat{a}^3) + \left( \frac{1}{2} k_* \tilde{b} - \frac{k_*}{\tilde{k}_0} \right) (\hat{a} + \hat{a}^{\dagger}) \implies |\tilde{\phi}_0\rangle_{k_* \to 0} = |\tilde{\phi}_0^{(0)}\rangle - \left( \frac{1}{2} k_* \tilde{b} - \frac{k_*}{\tilde{k}_0} \right) |\tilde{\phi}_1^{(0)}\rangle - \frac{k_* \tilde{b}}{\sqrt{6}} |\tilde{\phi}_3^{(0)}\rangle + O(k_*^2) \\
&\implies \frac{|\langle \psi_S | \phi_0 \rangle|^2}{3 \hbar^2 \int_0^{+\infty} dk \, k^2 e^{-\beta \tilde{\epsilon}_k}} \stackrel{k_* \to 0}{\sim} \frac{k_B T}{3 \tilde{m}_*} \left( \frac{2k_*}{\tilde{k}_0} \right)^2 \quad (B.11)
\end{align*}$$

et à l'ordre deux en $k_*$ pour obtenir une pulsation propre non nulle dans l'état fondamental :

$$\begin{align}
\tilde{\mathcal{H}}_1^{(2)} = \frac{k_*^2 \frac{d^2 D_{\parallel}}{dk^2}(\tilde{k}_0)}{2 D_{\parallel}(\tilde{k}_0)} \left[ -\frac{d}{dx} x^2 \frac{d}{dx} + \frac{1}{4} x^4 - \frac{3}{2} x^2 \right] &+ \frac{k_*^2 \frac{dD_{\parallel}}{dk}(\tilde{k}_0)}{D_{\parallel}(\tilde{k}_0)} \left[ \frac{\tilde{b}}{2}(x^4 - 3x^2) + \frac{1-x^2}{\tilde{k}_0} \right] \nonumber \\
&+ \frac{1}{4} k_*^2 \tilde{b}^2 x^4 + \frac{1}{2} k_*^2 \tilde{\ell}^2 (x^4 - 3x^2) \nonumber \\
&+ \frac{k_*^2}{\tilde{k}_0^2} (1 - \tilde{k}_0 \tilde{b}) x^2 + \frac{2 k_*^2 D_{\perp}(\tilde{k}_0)}{\tilde{k}_0^2 D_{\parallel}(\tilde{k}_0)} \implies \tilde{\omega}_0 \stackrel{k_* \to 0}{\sim} \langle \tilde{\phi}_0^{(0)} | \tilde{\mathcal{H}}_1^{(2)} | \tilde{\phi}_0^{(0)} \rangle + \sum_{n>0} \frac{|\langle \tilde{\phi}_n^{(0)} | \tilde{\mathcal{H}}_1^{(1)} | \tilde{\phi}_0^{(0)} \rangle|^2}{-n} = \frac{2 k_*^2 D_{\perp}(\tilde{k}_0)}{\tilde{k}_0^2 D_{\parallel}(\tilde{k}_0)} \quad (B.12)
\end{align}$$
sachant que $\langle x^2 \rangle = 1$, $\langle x^4 \rangle = 3$ et $\langle \frac{d}{dx} x^2 \frac{d}{dx} \rangle = -3/4$ où la moyenne est prise dans $|\tilde{\phi}_0^{(0)}\rangle$. On justifie ainsi le second terme entre crochets dans l'équation (85) et donc, au total, l'équivalent (86) du coefficient de diffusion spatiale 63 64.

Entre les cas $k_0 \equiv 0$ et $k_0 > 0$. Au point de transition entre ces deux cas, la masse effective $m_*$ est infinie, $1/m_* = 0$, la relation de dispersion $\epsilon_k$ varie quartiquement en nombre d'onde au voisinage de son minimum, et aucun des calculs du coefficient de diffusion spatiale présentés jusqu'ici dans ce travail ne s'applique. Pour raccorder les deux cas limites (par exemple lorsque la densité $\rho$ varie), il faut utiliser une approximation de degré 4 pour la relation de dispersion,

$$ \epsilon_k \approx \Delta_0 + \frac{\hbar^2 k^2}{2m_{\text{eff}}} + \frac{\hbar^2 k^2}{4m} (kl)^2 \quad (B.13) $$

où la longueur $l$ est essentiellement constante mais où la variation de la masse effective $m_{\text{eff}}$ en $k=0$ (à ne pas confondre avec celle $m_*$ à l'endroit du minimum) décrit le passage d'une relation de dispersion atteignant son minimum en $k=0$ (cas $1/m_{\text{eff}} > 0$) à une relation de dispersion atteignant son minimum en un nombre d'onde non nul (cas $1/m_{\text{eff}} < 0$). On fait ensuite un grossissement autour du point de transition, en regardant dans l'espace des $1/m_{\text{eff}}$ et des $k$ à des échelles $T^\gamma$ telles que les deux derniers termes de l'équation (B.13) soient du même ordre de grandeur $\approx k_B T$. Nous introduisons donc le nombre d'onde adimensionné $\kappa$ et le paramètre sans dimension $\delta$ qui pilote la transition :

$$ \kappa = k \left( \frac{\hbar^2 l^2}{mk_B T} \right)^{1/4} \quad \text{et} \quad \delta = \frac{m}{m_{\text{eff}}} \left( \frac{\hbar^2}{mk_B T l^2} \right)^{1/2} \quad (B.14) $$

Il reste à passer à la limite $T \to 0$ à $\kappa$ et $\delta$ fixés. On approxime la force moyenne et les coefficients de diffusion en impulsion par leur ordre dominant au point de transition, c'est-à-dire sous la forme réduite (53) par $\mathcal{F}(k) \approx -\alpha_0 e_k$ et $\mathcal{D}_{\perp}(k) \approx \mathcal{D}_{\parallel}(k) \approx \alpha_0$, où $\alpha_0$ est le coefficient de frottement réduit à la transition 65. Il reste à introduire le temps réduit

$$ \theta = t/t_0 \quad \text{avec} \quad t_0 \equiv \frac{\hbar}{mc^2} \frac{15}{\pi^5} \left( \frac{mc^2}{k_B T} \right)^{17/2} \frac{(\hbar\rho^{1/3}/mc)^6}{\alpha_0 mcl/\hbar} \quad (B.15) $$

pour que la marche au hasard de la quasi-particle $\gamma$ dans l'espace de Fourier soit décrite par un processus stochastique adimensionné universel à un paramètre $\delta$ et de même bruit gaussien que l'équation (59),

$$ d\kappa = -d\theta \mathbf{grad}_{\kappa} \varepsilon(\kappa) + \sqrt{2d\theta} \eta \quad \text{avec} \quad \varepsilon(\kappa) = \frac{1}{2}\kappa^2\delta + \frac{1}{4}\kappa^4 \quad (B.16) $$

<footnote n="63">63. Il existe une façon rapide et élégante de retrouver ces résultats. En restreignant l'équation de Fokker-Planck dépendant du temps aux distributions invariants par rotation, puis en appliquant le même changement de fonction qu'à $\Pi_1(k)$, on fait apparaître le hamiltonien fictif autoadjoint $\mathcal{H}_0$. Dans ce secteur isotrope, l'équation de Fokker-Planck admet une solution stationnaire non nulle $\Pi_0$. L'état fondamental $|\psi_0\rangle$ de $\mathcal{H}_0$ est donc connu, $\psi_0(k) = N_0 k \exp(-\beta \tilde{\epsilon}_k / 2)$ avec $N_0 = \left[ \int_0^{+\infty} dk \, k^2 \exp(-\beta \tilde{\epsilon}_k) \right]^{1/2}$, et sa valeur propre est exactement nulle. Or, le calcul le montre, on a simplement $\mathcal{H}_1 = \mathcal{H}_0 + 2D_{\perp}(k)/(\hbar k)^2$, si bien que $\tilde{\mathcal{H}}$ diffère de $\tilde{\mathcal{H}}_0$ par un terme du second ordre $2k_*^2 D_{\perp}(\tilde{k}_0 + xk_*)/[D_{\parallel}(\tilde{k}_0)(\tilde{k}_0 + xk_*)^2]$. L'état fondamental de $\tilde{\mathcal{H}}_1$ coïncide donc avec $|\psi_0\rangle$ à l'ordre sous-dominant, et l'on a $\langle \psi_S | \phi_0 \rangle \sim \langle \psi_S | \psi_0 \rangle$; au même ordre, on peut remplacer $d\tilde{\epsilon}_k/dk$ par $d\tilde{\epsilon}_k/dk$ dans $\langle k | \psi_S \rangle$. Une simple intégration par parties (il faut intégrer $d\tilde{\epsilon}_k/dk \exp(-\beta \tilde{\epsilon}_k)$ et dériver $k^2$), puis le remplacement des facteurs lentement variables $k$ au numérateur et $k^2$ au dénominateur par $\tilde{k}_0$ et $\tilde{k}_0^2$, redonnent alors l'équivalent du recouvrement de l'équation (B.11). De même, à l'ordre dominant, la pulsation propre $\omega_0$ de l'état fondamental de $\tilde{\mathcal{H}}_1$ s'obtient en traitant l'écart de $\tilde{\mathcal{H}}_1$ à $\mathcal{H}_0$ au premier ordre de la théorie des perturbations, $\omega_0 \sim \langle \psi_0 | 2D_{\perp}(k)/(\hbar k)^2 | \psi_0 \rangle = \int_0^{+\infty} dk \, 2D_{\perp}(k) \exp(-\beta \tilde{\epsilon}_k) / \int_0^{+\infty} dk \, (\hbar k)^2 \exp(-\beta \tilde{\epsilon}_k)$. En approximant les facteurs lentement variables $2D_{\perp}(k)$ et $(\hbar k)^2$ par leur valeur en $\tilde{k}_0$, on retrouve l'équivalent de $\omega_0$ donné dans l'équation (B.12).</footnote>

<footnote n="64">64. Nous avons aussi vérifié avec succès la validité à basse température du coefficient de diffusion spatiale prédit par l'approximation de Fokker-Planck en résolvant numériquement l'équation maîtresse (47) avec l'amplitude de diffusion $\phi - \gamma$ microréversible (74) pour le modèle simple $\epsilon_k = \Delta_k + (\hbar^2/2m_*)[(k^2 - k_0^2)/2k_0]^2$ ($m_* > 0, k_0 > 0$) et $R_k(u, u', w) = (1 - \eta u^2)(1 - \eta u'^2)$, dans lequel $\alpha = \frac{32}{1575}$ ($3\eta^2 - 10\eta + 15$)($15\eta^2 - 42\eta + 35$), $\mathcal{D}_{\perp}(k_0) = \frac{32}{1575}$ ($3\eta^2 - 10\eta + 15$)($3\eta^2 - 14\eta + 35$) $\neq \alpha$ et $\frac{d}{dk} \mathcal{D}_{\parallel}(k_0) = 0$ donc la force moyenne $F(k_0)$ est non nulle à l'ordre $T^9$ (sauf si $\eta = 0$). L'erreur commise par Fokker-Planck sur $\mathcal{D}^{\text{spa}}$ est < 10% pour tout $\eta \in [0, 3]$ dès que $k_*/k_0 < 1/20$. Le calcul utilise le théorème de régression quantique pour accéder à la fonction de corrélation de la vitesse, que l'on reporte ensuite dans l'expression (83) de $\mathcal{D}^{\text{spa}}$. Dans le secteur de moment cinétique nul, l'équation maîtresse s'écrit $\partial_t \Pi_0(k, t) = -\gamma(k)\Pi_0(k, t) + \int_0^{+\infty} dk' \gamma_0(k' \to k)\Pi_0(k', t)$; en pratique, on assure la conservation de la probabilité totale en calculant numériquement tous les taux d'alimentation $\gamma_0(k' \to k)$ puis en déduisant les taux de départ de la relation $\gamma(k) = \int_0^{+\infty} dk' (k'/k)^2 \gamma_0(k \to k')$. Dans le secteur de moment cinétique un, l'équation maîtresse s'écrit $\partial_t \Pi_1(k, t) = -\gamma(k)\Pi_1(k, t) + \int_0^{+\infty} dk' \gamma_1(k' \to k)\Pi_1(k', t)$ (les taux de départ sont inchangés). La microréversibilité impose $k^2 \exp(\beta \epsilon_k) \gamma_n(k' \to k) = k'^2 \exp(\beta \epsilon_k') \gamma_n(k \to k')$, $n = 0$ ou 1, ce qui permet de se ramener à des opérateurs intégraux à noyaux symétriques (numériquement avantageux) par le changement de fonction $\Pi_n(k) = k^{-1} \exp(-\beta \epsilon_k / 2) \psi_n(k)$. L'expression (B.5) s'applique alors, pourvu que l'on y remplace $\tilde{\epsilon}_k$ par $\epsilon_k$ et $-\mathcal{H}_1$ par l'opérateur intégral reliant $\partial_t \psi_1$ à $\psi_1$.</footnote>

<footnote n="65">65. On vérifie en effet dans l'équation (40) que, pour la relation de dispersion (B.13) et pour une limite finie de $d(m/m_{\text{eff}})/d\rho$ à la transition $\delta = 0$, $R_k(u, u', w) \to (1/2)(\tilde{\epsilon}_{x,0} + w\tilde{\epsilon}_{\rho,0})$ où $\tilde{\epsilon}_{x,0}$ et $\tilde{\epsilon}_{\rho,0}$ sont les valeurs en $\delta = 0$ de $\tilde{\epsilon}_x$ et $\tilde{\epsilon}_\rho$ définis par l'équation (45). C'est bien ce qu'on obtient en prenant $m/m_* = 0$ dans l'équation (46).</footnote>
Pour en extraire le coefficient de diffusion spatiale au sens de l'équation (83), nous passons par le théorème de régres-
sion et l'équation de Fokker-Planck comme au début de cette Annexe B et obtenons (en posant $\epsilon' = d\epsilon(\kappa)/d\kappa$)

$$
\mathcal{D}^{\text{spa}}_{\substack{\text{\delta fixé} \\ T \to 0}} \frac{\hbar}{m} \frac{\alpha_0^{-1} \Delta^{\text{spa}}(\delta)}{\frac{\pi^5}{15} \left(\frac{mc}{\hbar\phi^{1/3}}\right)^6 \left(\frac{\hbar c T}{mc^2}\right)^7} \quad \text{avec} \quad \Delta^{\text{spa}}(\delta) = \frac{\langle s | \hat{h}^{-1} | s \rangle}{3 \int_0^{+\infty} d\kappa \, \kappa^2 e^{-\epsilon(\kappa)}} \quad \text{et} \quad
\left\{
\begin{array}{@{}l@{}}
\hat{h} = -\frac{d^2}{d\kappa^2} + \frac{2}{\kappa^2} + \frac{1}{4}\epsilon'^2 - \frac{1}{\kappa}\epsilon' - \frac{1}{2}\epsilon'' \\
\langle \kappa | s \rangle = \kappa\epsilon' e^{-\epsilon(\kappa)/2}
\end{array}
\right.
\tag{B.17}
$$

Nous avons ici intégré formellement sur le temps la fonction de corrélation de la vitesse écrite comme dans l'équation (B.5), ce qui a fait apparaître l'inverse de l'opérateur hamiltonien fictif $\hat{h}$. Il se trouve que l'on sait déterminer analytiquement l'action de $\hat{h}^{-1}$ sur le vecteur source $|s\rangle$ (c'est facile à vérifier en faisant agir $\hat{h}$ sur le résultat) :

$$
\langle \kappa | \hat{h}^{-1} | s \rangle = \kappa^2 e^{-\epsilon(\kappa)/2} \implies \Delta^{\text{spa}}(\delta) = \frac{\int_0^{+\infty} d\kappa \, \kappa^3 \epsilon' e^{-\epsilon(\kappa)}}{3 \int_0^{+\infty} d\kappa \, \kappa^2 e^{-\epsilon(\kappa)}} \quad (\text{B.18})
$$

Une simple intégration par parties au numérateur de (B.18) conduit alors au résultat remarquable $\Delta^{\text{spa}}(\delta) \equiv 1,$$^{66}$ en accord parfait avec la prédiction (86), qui vaut donc partout dans la zone de transition entre les cas $k_0 \equiv 0$ et $k_0 > 0$.

Références

[1] T.J. Greytak, R.L. Woerner, « The two roton bound state », Journal de Physique Colloques, **33** (C1), C1-269 (1972).

[2] L.P. Pitaevskii, I.A. Fomin, « Structure of the bound-state spectrum of two rotons in superfluid helium », Zh. Eksp. Teor. Fiz. **65**, 2516 (1973).

[3] A. Nicolis, R. Penco, « Mutual interactions of phonons, rotons and gravity », Phys. Rev. B **97**, 134516 (2018).

[4] Y. Castin, A. Sinatra, H. Kurkjian, « Landau Phonon-Roton Theory Revisited for Superfluid He-4 and Fermi Gases », Phys. Rev. Lett. **119**, 260402 (2017).

[5] A.L. Gaunt, T.F. Schmidutz, I. Gotlibovych, R.P. Smith, Z. Hadzibabic, « Bose-Einstein condensation of atoms in a uniform potential », Phys. Rev. Lett. **110**, 200406 (2013).

[6] B. Mukherjee, Zhenjie Yan, P.B. Patel, Z. Hadzibabic, T. Yefsah, J. Struck, M.W. Zwierlein, « Homogeneous Atomic Fermi Gases », Phys. Rev. Lett. **118**, 123401 (2017).

[7] K. Hueck, N. Luick, L. Sobirey, J. Siegl, T. Lompe, H. Moritz, « Two-Dimensional Homogeneous Fermi Gases », Phys. Rev. Lett. **120**, 060402 (2018).

[8] P.B. Patel, Zhenjie Yan, B. Mukherjee, R.J. Fletcher, J. Struck, M.W. Zwierlein, « Universal Sound Diffusion in a Strongly Interacting Fermi Gas », arXiv:1909.02555 (2019).

[9] K.M. O'Hara, S.L. Hemmer, M.E. Gehm, S.R. Granade, J.E. Thomas, « Observation of a strongly interacting degenerate Fermi gas of atoms », Science **298**, 2179 (2002).

[10] T. Bourdel, J. Cubizolles, L. Khaykovich, K.M. Magalhães, S.J.J.M.F. Kokkelmans, G.V. Shlyapnikov, C. Salomon, « Measurement of the interaction energy near a Feshbach resonance in a $^6$Li Fermi gas », Phys. Rev. Lett. **91**, 020402 (2003).

[11] M. Bartenstein, A. Altmeyer, S. Riedl, S. Jochim, C. Chin, J.H. Denschlag, R. Grimm, « Collective excitations of a degenerate gas at the BEC-BCS crossover », Phys. Rev. Lett. **92**, 203201 (2004).

[12] M.W. Zwierlein, C.A. Stan, C.H. Schunck, S.M.F. Raupach, A.J. Kerman, W. Ketterle, « Condensation of pairs of fermionic atoms near a Feshbach resonance », Phys. Rev. Lett. **92**, 120403 (2004).

[13] S. Nascimbène, N. Navon, K.J. Jiang, F. Chevy, C. Salomon, « Exploring the thermodynamics of a universal Fermi gas », Nature **463**, 1057 (2010).

[14] M.J.H. Ku, A.T. Sommer, L.W. Cheuk, M.W. Zwierlein, « Revealing the superfluid lambda transition in the universal thermodynamics of a unitary Fermi gas », Science **335**, 563 (2012).

[15] W. Ketterle, « Making, probing and understanding ultracold Fermi gases », section 2, in *Ultra-cold Fermi gases*, Cours de l'école de physique Enrico Fermi 2006 de Varenne, édité par M. Inguscio, W. Ketterle, C. Salomon (SIF, Bologne, 2007).

[16] A. Schirotzek, Yong-il Shin, Ch.H. Schunck, W. Ketterle, « Determination of the Superfluid Gap in Atomic Fermi Gases by Quasiparticle Spectroscopy », Phys. Rev. Lett. **101**, 140403 (2008).

[17] J. Cubizolles, T. Bourdel, S.J.J.M.F. Kokkelmans, G.V. Shlyapnikov, C. Salomon, « Production of Long-Lived Ultracold Li₂ Molecules from a Fermi Gas », Phys. Rev. Lett. **91**, 240401 (2003).

[18] Y. Castin, A. Sinatra, H. Kurkjian, « Erratum: Landau Phonon-Roton Theory Revisited for Superfluid He-4 and Fermi Gases [Phys. Rev. Lett. **119**, 260402 (2017)] », Phys. Rev. Lett. **123**, 239904(E) (2019).

[19] B. Fåk, T. Keller, M. E. Zhitomirsky, A. L. Chernyshev, « Roton-phonon interaction in superfluid $^4$He », Phys. Rev. Lett. **109**, 155305 (2012).

[20] R. Combescot, M. Yu. Kagan, S. Stringari, « Collective mode of homogeneous superfluid Fermi gases in the BEC-BCS crossover », Phys. Rev. A **74**, 042717 (2006).

<footnote n="66">66. On l'aura compris, cette conclusion est indépendante de la forme particulière de la relation de dispersion réduite $\epsilon(\kappa)$.</footnote>
[21] H. Kurkjian, Y. Castin, A. Sinatra, « Concavity of the collective excitation branch of a Fermi gas in the BEC-BCS crossover », Phys. Rev. A **93**, 013623 (2016).

[22] Y. Castin, I. Ferrier-Barbut, C. Salomon, « La vitesse critique de Landau d'une particule dans un superfluide de fermions », Comptes Rendus Physique **16**, 241 (2015).

[23] L.P. Pitaevskii, « Properties of the Spectrum of Elementary Excitations near the Disintegration Threshold of the Excitations », Zh. Eksp. Teor. Fiz. **36**, 1168 (1959).

[24] L. Landau, I. Khalatnikov, « Teoriya vyazkosti Geliya-II », Zh. Eksp. Teor. Fiz. **19**, 637 (1949) [en anglais dans *Collected papers of L.D. Landau*, chapitre 69, p. 494-510, édité par D. ter Haar (Pergamon, New York, 1965)].

[25] N. Lerch, L. Bartosch, P. Kopietz, « Absence of Fermionic Quasiparticles in the Superfluid State of the Attractive Fermi Gas », Phys. Rev. Lett. **100**, 050403 (2008).

[26] H. Kurkjian, Y. Castin, A. Sinatra, « Three-phonon and four-phonon interaction processes in a pair-condensed Fermi gas », Annalen der Physik (Berlin) **529**, 1600352 (2017).

[27] R. Haussmann, M. Punk, W. Zwerger, « Spectral functions and rf response of ultracold fermionic atoms », Phys. Rev. A **80**, 063612 (2009).

[28] I.M. Khalatnikov, « Pogloshcheniye zvuka v gelii II », Zh. Eksp. Teor. Fiz. **20**, 243 (1950).

[29] I.M. Khalatnikov, D.M. Chernikova, « Relaxation phenomena in superfluid helium », Zh. Eksp. Teor. Fiz. **49**, 1957 (1965).

[30] D.T. Son, M. Wingate, « General coordinate invariance and conformal invariance in nonrelativistic physics: Unitary Fermi gas », Ann. Physics **321**, 197 (2006).

[31] G. Bighin, L. Salasnich, P.A. Marchetti, F. Toigo, « Beliaev damping of the Goldstone mode in atomic Fermi superfluids », Phys. Rev. A **92**, 023638 (2015).

[32] H. Kurkjian, Y. Castin, A. Sinatra, « Brouillage thermique d'un gaz cohérent de fermions », Comptes Rendus Physique **17**, 789 (2016).

[33] C. Cohen-Tannoudji, J. Dupont-Roc, G. Grynberg, *Processus d'interaction entre photons et atomes*, collection Savoirs Actuels (InterEditions/Editions du CNRS, Paris, 1988).

[34] Y. Castin, A. Sinatra, H. Kurkjian, « Le couplage de Landau phonon-roton revisité pour l'hélium 4 liquide et étendu aux gaz de fermions superfluides », hal-01570314v3 (2019).

[35] Y. Castin, « Étalement de la phase et cohérence temporelle d'un gaz de fermions condensé par paires à basse température », Comptes Rendus Physique **20**, 540 (2019).

[36] M. Marini, F. Pistolesi, G.C. Strinati, « Evolution from BCS superconductivity to Bose condensation: Analytic results for the crossover in three dimensions », Eur. Phys. J. B **1**, 151 (1998).

[37] P. Langevin, « Sur la théorie du mouvement brownien », Comptes Rendus de l'Académie des Sciences **146**, 530 (1908).

[38] Y. Nishida, D.T. Son, « *ε* Expansion for a Fermi gas at infinite scattering length », Phys. Rev. Lett. **97**, 050403 (2006).

[39] D.S. Petrov, C. Salomon, G.V. Shlyapnikov, « Scattering properties of weakly bound dimers of fermionic atoms », Phys. Rev. A **71**, 012708 (2005).

[40] I.V. Brodsky, A.V. Klaptsov, M. Yu. Kagan, R. Combescot, X. Leyronas, « Bound states of three and four resonantly interacting particles », Pis'ma v ZhETF **82**, 306 (2005) [JETP Letters **82**, 273 (2005)].

[41] G.V. Skorniakov, K.A. Ter-Martirosian, « Three-body problem for short range forces. I. Scattering of low energy neutrons by deuterons », Zh. Eksp. Teor. Phys. **31**, 775 (1957) [Sov. Phys. JETP **4**, 648 (1957)].

[42] J. Levinsen, D.S. Petrov, « Atom-dimer and dimer-dimer scattering in fermionic mixtures near a narrow Feshbach resonance », Eur. Phys. J. D **65**, 67 (2011).

[43] S. Giorgini, L.P. Pitaevskii, S. Stringari, « Theory of ultracold atomic Fermi gases », Rev. Mod. Phys. **80**, 1215 (2008).

[44] F. Alzetto, X. Leyronas, « Equation of state of a polarized Fermi gas in the Bose-Einstein-condensate limit », Phys. Rev. A **81**, 043604 (2010).