J. 5117
SESSION 2001
Filière Physique‐Chimie
PHYSIQUE
(ENS : Ulm)
Durée: 6 heures
Il est conseillé d’border les différentes parties dans leur ordre d’apparition, certaines questions ayant une résolution utile pour la suite \(du\) problème.
L’usage d’une calculatrice électronique de \(\rho oche\) à alimentation autonome, sans imprimante, est autorisé. Une seule calculatrice à la fois est admise sur la table et aucun échange n’est autorisé entre les candidats. Le candidat est prié d’accorder une importance particulière aux applications numériques.
Tournez la page S.V.P.
Des progrès importants effectués dans les techniques de refroidissement des gaz d’atomes neutres permettent depuis 1995 d’atteindre des températures tellement basses que ces gaz perdent toute viscosité. Par analogie avec l’hélium liquide, qui perd lui aussi toute viscosité à très basse température, on parle de gaz superfluides.
On se propose d’établir ici quelques propriétés des gaz superfluides et de comparer les prédictions obtenues aux résultats expérimentaux. La partie I est consacrée à l’étude des propriétés à l’équilibre thermodynamique d’un gaz superfluide en l’absence de potentiel extérieur. La partie II considère la situation réalisée en pratique d’un gaz superfluide à l’équilibre dans un potentiel harmonique créé par un champ magnétique. La partie III détermine la réponse du gaz superfluide à une faible perturbation du potentiel harmonique. On rappelle la valeur du nombre d’Avogadro, ${{\mathcal{N}}_{a}}=6,02\times {{10}^{23}}.$
1 Gaz superfluide en l’absence de potentiel extérieur
1.1 Équation d’état d’un gaz superfluide
On considère un gaz superfluide de \(N\) particules dans une enceinte de volume\(V\). Les parois de l’enceinte constituent un thermostat de température\(T\), avec lequel le gaz est à l’équilibre thermodynamique. Les particules interagissent deux à deux avec un potentiel d’interaction \(\mathcal{V}\left( {\overrightarrow {{r_i}} - \vec r} \right)\) , où \(\overrightarrow {{r_i}} ,\) \(\overrightarrow {{r_j}} \) sont les vecteurs positions de deux particules quelconques \(i\) et \(j\) du gaz. On suppose que le potentiel d’interaction est négligeable pour deux particules séparées par une distance macroscopique de l’ordre de la taille de l’enceinte. En d’autres termes, la portée du potentiel d’interaction est négligeable devant la taille de l’enceinte. De plus, la température \(T\) du gaz est si basse que l’énergie cinétique des particules est négligeable devant l’énergie d’interaction. On suppose donc que les particules sont au repos, avec des positions réparties aléatoirement avec une densité moyenne uniforme \(\rho = \frac{N}{V}.\)
(a) Calculer l’énergie moyenne \(Udu\) gaz en fonction du nombre de particules \(N,\) \(du\) volume \(V\) de l’enceinte et de la constante de couplage
\(\gamma = \mathop \smallint \limits_{ - \infty }^{ + \infty } d{r_x}\mathop \smallint \limits_{ - \infty }^{ + \infty } d{r_y}\mathop \smallint \limits_{ - \infty }^{ + \infty } d{r_z}\mathcal{V}\left( {\vec r} \right)\) (1)
Où \({r_x},\) \({r_y},\) \({r_z}\) sont les composantes du vecteur \({r^ \to }\) dans une base orthonormée \(\left( {{{\vec e}_x},{{\vec e}_y},{{\vec e}_z}} \right)\) .
(b) On admet que l’entropie du gaz tend vers zéro lorsque la température tend vers zéro. En déduire l’énergie libre\(Fdugaz\), dans la limite\(T \simeq O\), en fonction des quantités \(N,\) \(V,\) \(\gamma \). Montrer que la pression \(P\) du gaz est donnée par
\(P = \frac{1}{2}\gamma {\rho ^2}\) (2)
Exprimer le coefficient de compressibilité isotherme \(\chi \) du gaz et son potentiel chimique \(\mu \) en fonction de \(\gamma \) et \(\rho .\)
1.2 Quelques contraintes sur l’équilibre thermodynamique
Il s’agit de vérifier si l’état d’équilibre considéré à la partie 1.1 précédente est bien admissible d’un point de vue thermodynamique. Pour cela, on sépare par la pensée l’enceinte de volume \(V\) en deux parties 1 et 2 de volumes respectifs \({V_1}\) et\({V_2} = V - {V_1}\), et l’on partage le gaz en \({N_1}\) particules dans la zone 1 et \({N_2} = N - {N_1}\) particules dans la zone 2. Chaque zone est supposée être à l’équilibre avec le thermostat à la température\(T\), avec une densité uniforme de particules \({\rho _i} = \frac{{{N_i}}}{{{V_i}}},\) \(i = 1,2\). Les densités \({\rho _1}\) et \({\rho _2}\) peuvent a priori être différentes. Le raisonnement va s’appuyer sur la fonction thermodynamique ${{F}_{e}}\left( N,~T,~V \right)$ donnant l’énergie libre du gaz à l’équilibre thermodynamique pour \(N\) particules dans un volume \(V\) à la température \(T.\)
(a) Rappeler pourquoi l’énergie libre ${{F}_{e}}\left( N,~T,~V \right)$ est inférieure à l’énergie libre de tout état du système hors d’équilibre à \(N,\) \(T,\) \(V\) fixés.
(b) Exprimer l’énergie libre totale du gaz en fonction des énergies libres des parties 1 et 2 du gaz. En déduire une équation fonctionnelle vérifiée par \({F_e}.\)
(c) Traduire le fait que l’énergie libre du gaz est un minimum vis‐à‐vis d’un changement de \({V_1}\) et \({V_2}\) à \({N_1},\) \({N_2},\) \(V\) constants. Interpréter physiquement les résultats.
(d) Mêmes questions pour un changement de \({N_1}\) et \({N_2}\) à \({V_1},\) \({V_2},\) \(N\) constants.
(e) Appliquer les conditions obtenues en (c) et (d) à l’expression de l’énergie libre obtenue dans la partie 1.1. En déduire que le gaz doit avoir une densité uniforme à l’équilibre thermodynamique, et que seul un signe de la constante de couplage \(\gamma ,\) que l’on précisera, est acceptable.
1.3 Mélange de deux gaz superfluides
On mélange deux gaz superfluides d’espèces atomiques différentes, \({N_a}\) atomes de l’espèce \(a\) et \({N_b}\) atomes de l’espèce\(b\), dans la même enceinte de volume\(V\). Les deux gaz sont à la même température\(T\), que l’on supposera très basse au sens de la partie 1.1. Les atomes de l’espèce \(a\) (respectivement b) interagissent par un potentiel d’interaction \({\mathcal{V}_{aa}}\left( {\vec r} \right)\) (respectivement \({\mathcal{V}_{bb}}\left( {\vec r} \right)\) ), où \({r^ \to }\) est le rayon vecteur entre deux atomes. De plus, chaque atome de l’espèce \(a\) interagit avec chaque atome de l’espèce \(b\) par le potentiel \({\mathcal{V}_{ab}}\left( {\vec r} \right)\) . Comme dans la partie 1.1 on néglige l’énergie d’interaction de deux atomes quelconques séparés par une distance macroscopique de l’ordre de la taille de l’enceinte. Il s’agit de déterminer l’état d’équilibre du mélange, en choisissant parmi deux configurations possibles.
(a) Dans la configuration I, chaque gaz superfluide remplit toute l’enceinte avec une densité moyenne uniforme, \({\rho _a} = {N_a}/V\) pour l’espèce \(a\) et \({\rho _b} = {N_b}/V\) pour l’espèce\(b\). Calculer l’énergie libre \({F_{I,e}}\) de cette configuration en fonction des nombres de particules \({N_a},\) \({N_b},\) \(du\) volume \(V\) et des constantes de couplage
${{\gamma }_{aa}}~=~\underset{-\infty }{\overset{+\infty }{\mathop \int }}\,d{{r}_{x}}\underset{-\infty }{\overset{+\infty }{\mathop \int }}\,d{{r}_{y}}\underset{-\infty }{\overset{+\infty }{\mathop \int }}\,d{{r}_{z}}{{\mathcal{V}}_{aa}}\left( {\vec{r}} \right)$ , (3)
${{\gamma }_{bb}}~=~\underset{-\infty }{\overset{+\infty }{\mathop \int }}\,d{{r}_{x}}\underset{-\infty }{\overset{+\infty }{\mathop \int }}\,d{{r}_{y}}\underset{-\infty }{\overset{+\infty }{\mathop \int }}\,d{{r}_{z}}{{\mathcal{V}}_{bb}}\left( {\vec{r}} \right)$ , (4)
${{\gamma }_{ab}}~=~\underset{-\infty }{\overset{+\infty }{\mathop \int }}\,d{{r}_{x}}\underset{-\infty }{\overset{+\infty }{\mathop \int }}\,d{{r}_{y}}\underset{-\infty }{\overset{+\infty }{\mathop \int }}\,d{{r}_{z}}{{\mathcal{V}}_{ab}}\left( {\vec{r}} \right)$ . (5)
(b) Dans la configuration II, les deux gaz superfluides occupent des régions de l’enceinte disjointes et complémentaires, chaque gaz ayant une densité moyenne uniforme dans son domaine respectif. L’espèce \(a\) occupe une fraction \(\alpha du\) volume\(V\), l’espèce \(b\) une fraction \(1 - \alpha \). Calculer l’énergie libre \({F_{II}}\) de cette configuration, en fonction de \({N_a},\) \({N_b},\) \(\alpha ,\) \({\gamma _{aa}},\) \({\gamma _{bb}}\) et du volume \(V.\)
(c) À l’équilibre, le paramètre \(\alpha \) de la configuration II prend une valeur bien précise. À l’aide d’un des principes de la thermodynamique, déterminer cette valeur \({\alpha _e}\) en fonction de \({N_{a)}}{N_b},\) \({\gamma _{aa}}\) et\({\gamma _{bb}}\). En déduire\({F_{II,e}}\), énergie libre de la configuration II à l’équilibre, en fonction de \({N_a},\) \({N_b},\) \({\gamma _{aa}},\) \({\gamma _{bb}}\) et du volume \(V.\)
(d) Parmi les configurations I et II, laquelle est la plus favorable d’un point de vue thermodynamique? On précisera les différents cas possibles suivant les valeurs des constantes de couplage. Interpréter physiquement le résultat.
(e) Le groupe de Wolfgang Ketterle, au Massachusetts Institute of Technology de Boston (MIT), a étudié des mélanges de trois espèces différentes, appelées \(a,\) \(b,\) \(c\). Les constantes de couplage de chaque paire d’espèces sont données dans la table 1. À l’aide de la théorie développée ici, prédire si les paires d’espèces \(a - b,\) \(b - c\) et \(a - c\) sont miscibles ou pas.
TAB. \(1 - \) Dans l’expérience \(du\) MIT sur les mélanges de trois espèces \(a,\) \(b,\) \(c\) de \(gaz\) superfluides, constantes de couplage \({\gamma _{ij}}\) en \(J \cdot {m^3}\) entre les espèces \(i\) et \(j\), avec \(i,\) \(j = a,\) \(b\) ou \(c.\)
2 Gaz superfluide au repos dans un piège
2.1 Équilibre dans un piège de forme arbitraire
Dans les expériences, le gaz superfluide est piégé, c’est‐à‐dire qu’il est confiné par un champ de forces dérivant du potentiel \({\Phi _e}\left( {\vec r} \right)\) dépendant de la position \(\vec r\) des particules dans le gaz, mais indépendant du temps. La force de pesanteur est incluse dans ce champ de forces. Contrairement au modèle spatialement homogène de la partie 1, la densité \({\rho _e}\left( {\vec r} \right)\) et la pression \({P_e}\left( {\vec r} \right)du\) gaz à l’équilibre thermodynamique dépendent de la position \(\vec r.\)
(a) Rappeler la relation fondamentale de l’hydrostatique portant sur la pression à l’équilibre d’un fluide soumis à un champ de forces extérieures.
(b) A l’aide de l’équation d’état obtenue dans la partie 1, exprimer la densité de particules \({\rho _e}\left( {{r^ \to }} \right)\) en fonction de la constante de couplage \(\gamma \) définie par l’équation (1) et du potentiel de piégeage \({\Phi _e}\left( {\vec r} \right)\) , à une constante additive près.
(c) Indiquer comment déterminer la valeur de cette constante additive, sans chercher ici à la calculer explicitement.
2.2 Équilibre dans un piège harmonique
Le plus souvent, le gaz superfluide reste confiné suffisamment près du point où le potentiel de piégeage est minimum pour que le potentiel de piégeage soit très bien représenté par son approximation harmonique. Pour simplifier, on suppose de plus dans cette partie que le potentiel \({\Phi _e}\left( {\vec r} \right)\) est isotrope:
\({\Phi _e}\left( {\vec r} \right) = \frac{1}{2}m{\omega ^2}{r^2}\), (6)
où \(m\) est la masse d’une particule, \(r\) est le module du vecteur position $\vec{r}$ et \(\omega \) une constante.
(a) Rappeler la signification physique de la quantité \(\omega .\)
(b) Montrer que le gaz superfluide reste effectivement confiné au voisinage de \(\vec r = \vec 0.\) Quelle est la forme d’équilibre de la surface du gaz ? On appelle \(R\) la distance maximale au centre du piège accessible au gaz.
(c) Exprimer la densité du gaz \({\rho _e}\left( {\vec r} \right)\) en fonction de\(\vec r\), de \(R\) et de la densité \(\rho _e^0\) au centre du piège.
(d) On connait le nombre total \(N\) de particules dans le superfluide. En déduire \(\rho _e^0\) en fonction de \(N\) et de \(R.\)
(e) Finalement, exprimer \(\rho _e^0\) en fonction de \(m,\) \(\omega ,\) \(N\) et de la constante de couplage \(\gamma \) définie par l’équation (1).
2.3 Considérations thermodynamiques
On souhaite obtenir une interprétation thermodynamique de l’état d’équilibre du gaz obtenu à la partie précédente 2.2.
(a) Calculer l’énergie potentielle harmonique \({E_{harm}}\) du gaz en fonction de \(N\) et \(\gamma \rho _e^0.\)
(b) Faire de même pour l’énergie d’interaction \({E_{int}}\) du gaz due au potentiel d’interaction \(\mathcal{V}\left( {\vec r} \right)\) de la partie 1.1, en supposant que \(\mathcal{V}\left( {\vec r} \right)\) est négligeable pour deux particules séparées d’une distance \(r\) macroscopique de l’ordre de \(R.\)
(c) Exprimer l’énergie libre du gaz à \(T \simeq 0\) en fonction de \(\gamma \rho _e^0\) et\(N\). Montrer que la quantité \(\gamma \rho _e^0\) n’est autre que le potentiel chimique du gaz.
(d) Retrouver la condition d’équilibre obtenue par l’hydrostatique à l’aide d’un argument thermodynamique.
2.4 Comparaison aux expériences
On compare les prédictions théoriques précédentes à quelques résultats expérimentaux obtenus avec un gaz superfluide d’atomes de sodium \({}^{23}Na\). Ces atomes interagissent avec une constante de couplage \(\gamma = 1,0 \times {10^{ - 50}}J \cdot {m^3}.\)
(a) Le groupe de Lene Hau, au Rowland Institute (Boston), a mesuré un diamètre maximal de 73 \(\mu m\) pour un nuage superfluide de \(N = 1,6 \times {10^6}\) atomes de sodium dans un piège harmonique de paramètre \(\omega = 87\) rad\( \cdot {s^{ - 1}}\). Comparer à la valeur attendue théoriquement.
(b) On coupe brutalement le potentiel de piégeage. On constate alors que le gaz entre en expansion, bien que les particules soient initialement au repos, à\(T \simeq 0\). Comment expliquer ce phénomène?
(c) Après un temps assez long pour que la vitesse d’expansion du gaz ait atteint un régime stationnaire, on mesure l’énergie cinétique d’expansion\({E_{cin}}\)du gaz. Les valeurs obtenues par le groupe de Wolfgang Ketterle, au MIT, de l’énergie d’expansion par particule pour différentes valeurs du nombre de particules dans le gaz, sont données dans la table 2. À l’aide d’une régression linéaire des résultats expérimentaux en échelle log‐log, montrer que \({E_{cin}}\) varie approximativement comme une puissance de \(N\) avec un exposant \(\delta \) que l’on précisera.
(d) Donner l’expression de \({E_{cin}}\) prédite par la théorie. Comparer la valeur de \(\delta \) obtenue expérimentalement à la prédiction théorique.
TAB. \(2\) —Énergie d’expansion par particule en Joule, mesurée au MIT en fonction \(du\) nombre de particules \(N\) dans le \(gaz\) superfluide.
3 Gaz superfluide en mouvement
3.1 Équations du mouvement
Le gaz superfluide, initialement à l’équilibre, est maintenant soumis à un potentiel extérieur $\Phi \left( \vec{r},~t \right)$ pouvant dépendre du temps\(t\). Pour décrire le gaz superfluide en mouvement, on l’assimile à un fluide sans viscosité. On note $\rho \left( \vec{r},~t \right),$ $P\left( \vec{r},~t \right),\vec{v}\left( \vec{r},~t \right)$ la densité volumique, la pression et le champ de vitesse du fluide.
(a) Rappeler l’équation de continuité portant sur la densité de particules $\rho \left( \vec{r},~t \right)$ et le champ de vitesse $\vec{v}\left( \vec{r},~t \right)$ , traduisant la conservation de la matière.
(b) Rappeler l’équation d’Euler donnant l’évolution du champ de vitesse $\vec{v}\left( \vec{r},~t \right)$ sous l’effet du potentiel extérieur $\Phi \left( \vec{r},~t \right)$ pour un fluide sans viscosité de pression $P\left( \vec{r},~t \right)$ .
(c) À l’aide de l’équation d’état obtenue dans la partie 1 et sous l’hypothèse d’un équilibre local du gaz, transformer l’équation d’Euler en une équation sur $\vec{v}\left( \vec{r},~t \right)$ et $p\left( \vec{r},~t \right)$ .
3.2 Régime de réponse linéaire
On applique au gaz superfluide une perturbation de potentiel \(\delta \Phi \) sur l’intervalle de temps $\left[ 0,~\tau \right]$:
$\Phi \left( \vec{r},~t \right)={{\Phi }_{e}}\left( {\vec{r}} \right)+\delta \Phi \left( \vec{r},~t \right)$ (7)
où \({\Phi _e}\left( {\vec r} \right)\) décrit le piège statique de la partie 2 et $\left| \delta \Phi \left( \vec{r},~t \right) \right|\ll {{\Phi }_{e}}\left( {\vec{r}} \right)$ . Cette perturbation entraîne des déviations $\delta p\left( \vec{r},~t \right)$ et $\delta \vec{v}\left( \vec{r},~t \right)$ de la densité et du champ de vitesse de leurs valeurs à l’équilibre \({\rho _e}\left( {\vec r} \right)\) et \(\vec v\left( {\vec r} \right) = \vec 0.\)
(a) En négligeant les termes quadratiques en les écarts à l’équilibre, obtenir des équations d’évolution linéaires pour \(\delta p\) et\(\delta \vec v\), sans chercher à les résoudre. Pour simplifier, on se limitera aux instants ultérieurs à la perturbation, \(t > \tau .\)
(b) Montrer qu’après avoir dérivé une fois par rapport au temps l’équation linéaire sur\(\delta p\), il est possible d’éliminer\(\delta \vec v\). En déduire une équation d’évolution portant seulement sur \(\delta \rho .\)
3.3 Modes propres d’un gaz superfluide homogène
On suppose que le gaz superfluide à l’équilibre remplit tout l’espace avec une densité uniforme\({\rho _e}\). On souhaite déterminer les modes propres de l’équation d’évolution sur \(\delta \rho .\) Pour cela, on considère une solution de la forme:
$\delta \rho \left( \vec{r},~t \right)=A\text{ }\!\!~\!\!\text{ cos }\!\!~\!\!\text{ }\left( \vec{k}\cdot \vec{r}-\Omega t \right)$ (8)
où \(A,\) \(\Omega \) et le vecteur réel \(\vec k\) sont des constantes.
(a) Vérifier que la forme proposée convient à condition que \({\Omega ^2}\) soit une certaine fonction de \(\vec k\) que l’on précisera.
(b) Interpréter physiquement le résultat obtenu. On précisera à quelle vitesse une perturbation appliquée en un point du gaz superfluide se propage. Comment cette vitesse dépend‐elle de la forme de la perturbation?
(c) Après envoi d’une impulsion laser dans un gaz superfluide d’atomes de sodium \({}^{23}Na\)initialement au repos, l’équipe du MIT a constaté la propagation d’ondes de densité à une vitesse de\(1 \times {10^{ - 2}}m/s\). On suppose que le gaz était initialement quasi homogène avec une densité moyenne de \(3,8 \times {10^{20}}\) atomes par\({m^3}\). Comparer la vitesse mesurée à la prédiction théorique.
3.4 Modes propres dans un piège harmonique
Le gaz superfluide est initialement au repos dans un potentiel de piégeage \({\Phi _e}\left( {{r^ \to }} \right)\) harmonique anisotrope:
\({\Phi _e}\left( {\vec r} \right) = \frac{1}{2}m\mathop \sum \limits_{\alpha = x,y,z}^{} \omega _\alpha ^2r_\alpha ^2\) (9)
où \({r_\alpha }\) est la composante du vecteur position $\vec{r}$ sur le vecteur de base \(\vec e,\) \(\alpha = x,\) \(y,\) \(z\). On excite le gaz en modifiant faiblement l’un des paramètres \({\omega _a}\) pendant une durée\(\tau \). On admet qu’on excite ainsi des modes propres du gaz superfluide de la forme
$\delta \rho \left( {{r}^{\to }},~t \right)=\left[ B+\underset{\alpha =x,y,z}{\overset{{}}{\mathop \sum }}\,{{A}_{\alpha }}r_{\alpha }^{2} \right]\text{ }\!\!~\!\!\text{ cos }\!\!~\!\!\text{ }\Omega t$ (10)
où \(B\), les coefficients \({A_\alpha }\) et la pulsation \(\Omega \) sont des constantes.
(a) Exprimer la densité \({\rho _e}\left( {\vec r} \right)\) du gaz à l’équilibre à une constante additive près, en fonction des \({r_\alpha },\) \({\omega _\alpha }\), de \(m\) et de \(\gamma .\)
(b) À quelles conditions sur les coefficients \({A_\alpha }\) la forme proposée pour \(\delta \rho \) satisfait‐elle l’équation d’évolution obtenue à la partie 3.2? Ces conditions peuvent‐elles être satisfaites quelle que soit la pulsation \(\Omega \)?
(c) Au MIT, le gaz superfluide, initialement piégé dans un potentiel à symétrie de révolution d’axe\(z\), est excité pendant la durée \(\tau \) par une petite modification du paramètre \({\omega _z}\) seulement. Expliquer pourquoi on a alors\({A_x} = {A_y}\), et écrire les conditions que doivent satisfaire \({A_x}\) et \({A_z}.\)
(d) En déduire que \(X = \frac{{{\Omega ^2}}}{{\omega _z^2}}\) doit vérifier une équation du second degré que l’on précisera. Donner les solutions de cette équation en fonction du paramètre \(\eta = \frac{{{\omega _x}}}{{{\omega _z}}}.\)
(e) Dans l’expérience du MIT, le paramètre \(\eta \) vaut 13,60. Donner les pulsations propres des deux modes excités, ainsi que les valeurs correspondantes du rapport \(\frac{{{A_z}}}{{{A_x}}}.\) Expérimentalement, on constate que le gaz bat à la pulsation \(1,569 \pm 0,004{\omega _z}.\) Quel est l’écart relatif entre la prédiction théorique et le résultat expérimental?
(f) À quelle fréquence devrait battre le gaz en l’absence d’interaction entre les particules? Qu’en concluez‐vous?
SESSION 2001
Filière Physique‐Chimie
PHYSIQUE
(ENS : Ulm)
Durée: 6 heures
Il est conseillé d’border les différentes parties dans leur ordre d’apparition, certaines questions ayant une résolution utile pour la suite \(du\) problème.
L’usage d’une calculatrice électronique de \(\rho oche\) à alimentation autonome, sans imprimante, est autorisé. Une seule calculatrice à la fois est admise sur la table et aucun échange n’est autorisé entre les candidats. Le candidat est prié d’accorder une importance particulière aux applications numériques.
Tournez la page S.V.P.
Des progrès importants effectués dans les techniques de refroidissement des gaz d’atomes neutres permettent depuis 1995 d’atteindre des températures tellement basses que ces gaz perdent toute viscosité. Par analogie avec l’hélium liquide, qui perd lui aussi toute viscosité à très basse température, on parle de gaz superfluides.
On se propose d’établir ici quelques propriétés des gaz superfluides et de comparer les prédictions obtenues aux résultats expérimentaux. La partie I est consacrée à l’étude des propriétés à l’équilibre thermodynamique d’un gaz superfluide en l’absence de potentiel extérieur. La partie II considère la situation réalisée en pratique d’un gaz superfluide à l’équilibre dans un potentiel harmonique créé par un champ magnétique. La partie III détermine la réponse du gaz superfluide à une faible perturbation du potentiel harmonique. On rappelle la valeur du nombre d’Avogadro, ${{\mathcal{N}}_{a}}=6,02\times {{10}^{23}}.$
1.1 Équation d’état d’un gaz superfluide
On considère un gaz superfluide de \(N\) particules dans une enceinte de volume\(V\). Les parois de l’enceinte constituent un thermostat de température\(T\), avec lequel le gaz est à l’équilibre thermodynamique. Les particules interagissent deux à deux avec un potentiel d’interaction \(\mathcal{V}\left( {\overrightarrow {{r_i}} - \vec r} \right)\) , où \(\overrightarrow {{r_i}} ,\) \(\overrightarrow {{r_j}} \) sont les vecteurs positions de deux particules quelconques \(i\) et \(j\) du gaz. On suppose que le potentiel d’interaction est négligeable pour deux particules séparées par une distance macroscopique de l’ordre de la taille de l’enceinte. En d’autres termes, la portée du potentiel d’interaction est négligeable devant la taille de l’enceinte. De plus, la température \(T\) du gaz est si basse que l’énergie cinétique des particules est négligeable devant l’énergie d’interaction. On suppose donc que les particules sont au repos, avec des positions réparties aléatoirement avec une densité moyenne uniforme \(\rho = \frac{N}{V}.\)
(a) Calculer l’énergie moyenne \(Udu\) gaz en fonction du nombre de particules \(N,\) \(du\) volume \(V\) de l’enceinte et de la constante de couplage
\(\gamma = \mathop \smallint \limits_{ - \infty }^{ + \infty } d{r_x}\mathop \smallint \limits_{ - \infty }^{ + \infty } d{r_y}\mathop \smallint \limits_{ - \infty }^{ + \infty } d{r_z}\mathcal{V}\left( {\vec r} \right)\) (1)
Où \({r_x},\) \({r_y},\) \({r_z}\) sont les composantes du vecteur \({r^ \to }\) dans une base orthonormée \(\left( {{{\vec e}_x},{{\vec e}_y},{{\vec e}_z}} \right)\) .
(b) On admet que l’entropie du gaz tend vers zéro lorsque la température tend vers zéro. En déduire l’énergie libre\(Fdugaz\), dans la limite\(T \simeq O\), en fonction des quantités \(N,\) \(V,\) \(\gamma \). Montrer que la pression \(P\) du gaz est donnée par
\(P = \frac{1}{2}\gamma {\rho ^2}\) (2)
Exprimer le coefficient de compressibilité isotherme \(\chi \) du gaz et son potentiel chimique \(\mu \) en fonction de \(\gamma \) et \(\rho .\)
1.2 Quelques contraintes sur l’équilibre thermodynamique
Il s’agit de vérifier si l’état d’équilibre considéré à la partie 1.1 précédente est bien admissible d’un point de vue thermodynamique. Pour cela, on sépare par la pensée l’enceinte de volume \(V\) en deux parties 1 et 2 de volumes respectifs \({V_1}\) et\({V_2} = V - {V_1}\), et l’on partage le gaz en \({N_1}\) particules dans la zone 1 et \({N_2} = N - {N_1}\) particules dans la zone 2. Chaque zone est supposée être à l’équilibre avec le thermostat à la température\(T\), avec une densité uniforme de particules \({\rho _i} = \frac{{{N_i}}}{{{V_i}}},\) \(i = 1,2\). Les densités \({\rho _1}\) et \({\rho _2}\) peuvent a priori être différentes. Le raisonnement va s’appuyer sur la fonction thermodynamique ${{F}_{e}}\left( N,~T,~V \right)$ donnant l’énergie libre du gaz à l’équilibre thermodynamique pour \(N\) particules dans un volume \(V\) à la température \(T.\)
(a) Rappeler pourquoi l’énergie libre ${{F}_{e}}\left( N,~T,~V \right)$ est inférieure à l’énergie libre de tout état du système hors d’équilibre à \(N,\) \(T,\) \(V\) fixés.
(b) Exprimer l’énergie libre totale du gaz en fonction des énergies libres des parties 1 et 2 du gaz. En déduire une équation fonctionnelle vérifiée par \({F_e}.\)
(c) Traduire le fait que l’énergie libre du gaz est un minimum vis‐à‐vis d’un changement de \({V_1}\) et \({V_2}\) à \({N_1},\) \({N_2},\) \(V\) constants. Interpréter physiquement les résultats.
(d) Mêmes questions pour un changement de \({N_1}\) et \({N_2}\) à \({V_1},\) \({V_2},\) \(N\) constants.
(e) Appliquer les conditions obtenues en (c) et (d) à l’expression de l’énergie libre obtenue dans la partie 1.1. En déduire que le gaz doit avoir une densité uniforme à l’équilibre thermodynamique, et que seul un signe de la constante de couplage \(\gamma ,\) que l’on précisera, est acceptable.
On mélange deux gaz superfluides d’espèces atomiques différentes, \({N_a}\) atomes de l’espèce \(a\) et \({N_b}\) atomes de l’espèce\(b\), dans la même enceinte de volume\(V\). Les deux gaz sont à la même température\(T\), que l’on supposera très basse au sens de la partie 1.1. Les atomes de l’espèce \(a\) (respectivement b) interagissent par un potentiel d’interaction \({\mathcal{V}_{aa}}\left( {\vec r} \right)\) (respectivement \({\mathcal{V}_{bb}}\left( {\vec r} \right)\) ), où \({r^ \to }\) est le rayon vecteur entre deux atomes. De plus, chaque atome de l’espèce \(a\) interagit avec chaque atome de l’espèce \(b\) par le potentiel \({\mathcal{V}_{ab}}\left( {\vec r} \right)\) . Comme dans la partie 1.1 on néglige l’énergie d’interaction de deux atomes quelconques séparés par une distance macroscopique de l’ordre de la taille de l’enceinte. Il s’agit de déterminer l’état d’équilibre du mélange, en choisissant parmi deux configurations possibles.
(a) Dans la configuration I, chaque gaz superfluide remplit toute l’enceinte avec une densité moyenne uniforme, \({\rho _a} = {N_a}/V\) pour l’espèce \(a\) et \({\rho _b} = {N_b}/V\) pour l’espèce\(b\). Calculer l’énergie libre \({F_{I,e}}\) de cette configuration en fonction des nombres de particules \({N_a},\) \({N_b},\) \(du\) volume \(V\) et des constantes de couplage
${{\gamma }_{aa}}~=~\underset{-\infty }{\overset{+\infty }{\mathop \int }}\,d{{r}_{x}}\underset{-\infty }{\overset{+\infty }{\mathop \int }}\,d{{r}_{y}}\underset{-\infty }{\overset{+\infty }{\mathop \int }}\,d{{r}_{z}}{{\mathcal{V}}_{aa}}\left( {\vec{r}} \right)$ , (3)
${{\gamma }_{bb}}~=~\underset{-\infty }{\overset{+\infty }{\mathop \int }}\,d{{r}_{x}}\underset{-\infty }{\overset{+\infty }{\mathop \int }}\,d{{r}_{y}}\underset{-\infty }{\overset{+\infty }{\mathop \int }}\,d{{r}_{z}}{{\mathcal{V}}_{bb}}\left( {\vec{r}} \right)$ , (4)
${{\gamma }_{ab}}~=~\underset{-\infty }{\overset{+\infty }{\mathop \int }}\,d{{r}_{x}}\underset{-\infty }{\overset{+\infty }{\mathop \int }}\,d{{r}_{y}}\underset{-\infty }{\overset{+\infty }{\mathop \int }}\,d{{r}_{z}}{{\mathcal{V}}_{ab}}\left( {\vec{r}} \right)$ . (5)
(b) Dans la configuration II, les deux gaz superfluides occupent des régions de l’enceinte disjointes et complémentaires, chaque gaz ayant une densité moyenne uniforme dans son domaine respectif. L’espèce \(a\) occupe une fraction \(\alpha du\) volume\(V\), l’espèce \(b\) une fraction \(1 - \alpha \). Calculer l’énergie libre \({F_{II}}\) de cette configuration, en fonction de \({N_a},\) \({N_b},\) \(\alpha ,\) \({\gamma _{aa}},\) \({\gamma _{bb}}\) et du volume \(V.\)
(c) À l’équilibre, le paramètre \(\alpha \) de la configuration II prend une valeur bien précise. À l’aide d’un des principes de la thermodynamique, déterminer cette valeur \({\alpha _e}\) en fonction de \({N_{a)}}{N_b},\) \({\gamma _{aa}}\) et\({\gamma _{bb}}\). En déduire\({F_{II,e}}\), énergie libre de la configuration II à l’équilibre, en fonction de \({N_a},\) \({N_b},\) \({\gamma _{aa}},\) \({\gamma _{bb}}\) et du volume \(V.\)
(d) Parmi les configurations I et II, laquelle est la plus favorable d’un point de vue thermodynamique? On précisera les différents cas possibles suivant les valeurs des constantes de couplage. Interpréter physiquement le résultat.
(e) Le groupe de Wolfgang Ketterle, au Massachusetts Institute of Technology de Boston (MIT), a étudié des mélanges de trois espèces différentes, appelées \(a,\) \(b,\) \(c\). Les constantes de couplage de chaque paire d’espèces sont données dans la table 1. À l’aide de la théorie développée ici, prédire si les paires d’espèces \(a - b,\) \(b - c\) et \(a - c\) sont miscibles ou pas.
| a | b | c | |
|---|---|---|---|
| a b c |
\(1,0 \times {10^{ - 50}}\) \(1,0 \times {10^{ - 50}}\) \(9,4 \times {10^{ - 51}}\) |
\(1,0 \times {10^{ - 50}}\) \(9,7 \times {10^{ - 51}}\) \(1,0 \times {10^{ - 50}}\) |
\(9,4 \times {10^{ - 51}}\) \(1,0 \times {10^{ - 50}}\) \(1,0 \times {10^{ - 50}}\) |
2.1 Équilibre dans un piège de forme arbitraire
Dans les expériences, le gaz superfluide est piégé, c’est‐à‐dire qu’il est confiné par un champ de forces dérivant du potentiel \({\Phi _e}\left( {\vec r} \right)\) dépendant de la position \(\vec r\) des particules dans le gaz, mais indépendant du temps. La force de pesanteur est incluse dans ce champ de forces. Contrairement au modèle spatialement homogène de la partie 1, la densité \({\rho _e}\left( {\vec r} \right)\) et la pression \({P_e}\left( {\vec r} \right)du\) gaz à l’équilibre thermodynamique dépendent de la position \(\vec r.\)
(a) Rappeler la relation fondamentale de l’hydrostatique portant sur la pression à l’équilibre d’un fluide soumis à un champ de forces extérieures.
(b) A l’aide de l’équation d’état obtenue dans la partie 1, exprimer la densité de particules \({\rho _e}\left( {{r^ \to }} \right)\) en fonction de la constante de couplage \(\gamma \) définie par l’équation (1) et du potentiel de piégeage \({\Phi _e}\left( {\vec r} \right)\) , à une constante additive près.
(c) Indiquer comment déterminer la valeur de cette constante additive, sans chercher ici à la calculer explicitement.
2.2 Équilibre dans un piège harmonique
Le plus souvent, le gaz superfluide reste confiné suffisamment près du point où le potentiel de piégeage est minimum pour que le potentiel de piégeage soit très bien représenté par son approximation harmonique. Pour simplifier, on suppose de plus dans cette partie que le potentiel \({\Phi _e}\left( {\vec r} \right)\) est isotrope:
\({\Phi _e}\left( {\vec r} \right) = \frac{1}{2}m{\omega ^2}{r^2}\), (6)
où \(m\) est la masse d’une particule, \(r\) est le module du vecteur position $\vec{r}$ et \(\omega \) une constante.
(a) Rappeler la signification physique de la quantité \(\omega .\)
(b) Montrer que le gaz superfluide reste effectivement confiné au voisinage de \(\vec r = \vec 0.\) Quelle est la forme d’équilibre de la surface du gaz ? On appelle \(R\) la distance maximale au centre du piège accessible au gaz.
(c) Exprimer la densité du gaz \({\rho _e}\left( {\vec r} \right)\) en fonction de\(\vec r\), de \(R\) et de la densité \(\rho _e^0\) au centre du piège.
(d) On connait le nombre total \(N\) de particules dans le superfluide. En déduire \(\rho _e^0\) en fonction de \(N\) et de \(R.\)
(e) Finalement, exprimer \(\rho _e^0\) en fonction de \(m,\) \(\omega ,\) \(N\) et de la constante de couplage \(\gamma \) définie par l’équation (1).
2.3 Considérations thermodynamiques
On souhaite obtenir une interprétation thermodynamique de l’état d’équilibre du gaz obtenu à la partie précédente 2.2.
(a) Calculer l’énergie potentielle harmonique \({E_{harm}}\) du gaz en fonction de \(N\) et \(\gamma \rho _e^0.\)
(b) Faire de même pour l’énergie d’interaction \({E_{int}}\) du gaz due au potentiel d’interaction \(\mathcal{V}\left( {\vec r} \right)\) de la partie 1.1, en supposant que \(\mathcal{V}\left( {\vec r} \right)\) est négligeable pour deux particules séparées d’une distance \(r\) macroscopique de l’ordre de \(R.\)
(c) Exprimer l’énergie libre du gaz à \(T \simeq 0\) en fonction de \(\gamma \rho _e^0\) et\(N\). Montrer que la quantité \(\gamma \rho _e^0\) n’est autre que le potentiel chimique du gaz.
(d) Retrouver la condition d’équilibre obtenue par l’hydrostatique à l’aide d’un argument thermodynamique.
On compare les prédictions théoriques précédentes à quelques résultats expérimentaux obtenus avec un gaz superfluide d’atomes de sodium \({}^{23}Na\). Ces atomes interagissent avec une constante de couplage \(\gamma = 1,0 \times {10^{ - 50}}J \cdot {m^3}.\)
(a) Le groupe de Lene Hau, au Rowland Institute (Boston), a mesuré un diamètre maximal de 73 \(\mu m\) pour un nuage superfluide de \(N = 1,6 \times {10^6}\) atomes de sodium dans un piège harmonique de paramètre \(\omega = 87\) rad\( \cdot {s^{ - 1}}\). Comparer à la valeur attendue théoriquement.
(b) On coupe brutalement le potentiel de piégeage. On constate alors que le gaz entre en expansion, bien que les particules soient initialement au repos, à\(T \simeq 0\). Comment expliquer ce phénomène?
(c) Après un temps assez long pour que la vitesse d’expansion du gaz ait atteint un régime stationnaire, on mesure l’énergie cinétique d’expansion\({E_{cin}}\)du gaz. Les valeurs obtenues par le groupe de Wolfgang Ketterle, au MIT, de l’énergie d’expansion par particule pour différentes valeurs du nombre de particules dans le gaz, sont données dans la table 2. À l’aide d’une régression linéaire des résultats expérimentaux en échelle log‐log, montrer que \({E_{cin}}\) varie approximativement comme une puissance de \(N\) avec un exposant \(\delta \) que l’on précisera.
(d) Donner l’expression de \({E_{cin}}\) prédite par la théorie. Comparer la valeur de \(\delta \) obtenue expérimentalement à la prédiction théorique.
| \(N\) | \({E_{cin}}/N\) |
|---|---|
| $6,~7\times {{10}^{4}}$ $1,~3\times {{10}^{5}}$ $1,~1\times {{10}^{6}}$ $1,~3\times {{10}^{6}}$ $2,~8\times {{10}^{6}}$ $3,~3~\times {{10}^{6}}$ \(3,9 \times {10^6}\) $4,~3\times {{10}^{6}}$ |
$3,~3\times {{10}^{-31}}$ \(4,6 \times {10^{ - 31}}\) $1,~1\times {{10}^{-30}}$ $1,~1\times {{10}^{-30}}$ $1,~7\times {{10}^{-30}}$ $1,~8\times {{10}^{-30}}$ $1,~9\times {{10}^{-30}}$ $2,~0\times {{10}^{-30}}$ |
3 Gaz superfluide en mouvement
3.1 Équations du mouvement
Le gaz superfluide, initialement à l’équilibre, est maintenant soumis à un potentiel extérieur $\Phi \left( \vec{r},~t \right)$ pouvant dépendre du temps\(t\). Pour décrire le gaz superfluide en mouvement, on l’assimile à un fluide sans viscosité. On note $\rho \left( \vec{r},~t \right),$ $P\left( \vec{r},~t \right),\vec{v}\left( \vec{r},~t \right)$ la densité volumique, la pression et le champ de vitesse du fluide.
(a) Rappeler l’équation de continuité portant sur la densité de particules $\rho \left( \vec{r},~t \right)$ et le champ de vitesse $\vec{v}\left( \vec{r},~t \right)$ , traduisant la conservation de la matière.
(b) Rappeler l’équation d’Euler donnant l’évolution du champ de vitesse $\vec{v}\left( \vec{r},~t \right)$ sous l’effet du potentiel extérieur $\Phi \left( \vec{r},~t \right)$ pour un fluide sans viscosité de pression $P\left( \vec{r},~t \right)$ .
(c) À l’aide de l’équation d’état obtenue dans la partie 1 et sous l’hypothèse d’un équilibre local du gaz, transformer l’équation d’Euler en une équation sur $\vec{v}\left( \vec{r},~t \right)$ et $p\left( \vec{r},~t \right)$ .
On applique au gaz superfluide une perturbation de potentiel \(\delta \Phi \) sur l’intervalle de temps $\left[ 0,~\tau \right]$:
$\Phi \left( \vec{r},~t \right)={{\Phi }_{e}}\left( {\vec{r}} \right)+\delta \Phi \left( \vec{r},~t \right)$ (7)
où \({\Phi _e}\left( {\vec r} \right)\) décrit le piège statique de la partie 2 et $\left| \delta \Phi \left( \vec{r},~t \right) \right|\ll {{\Phi }_{e}}\left( {\vec{r}} \right)$ . Cette perturbation entraîne des déviations $\delta p\left( \vec{r},~t \right)$ et $\delta \vec{v}\left( \vec{r},~t \right)$ de la densité et du champ de vitesse de leurs valeurs à l’équilibre \({\rho _e}\left( {\vec r} \right)\) et \(\vec v\left( {\vec r} \right) = \vec 0.\)
(a) En négligeant les termes quadratiques en les écarts à l’équilibre, obtenir des équations d’évolution linéaires pour \(\delta p\) et\(\delta \vec v\), sans chercher à les résoudre. Pour simplifier, on se limitera aux instants ultérieurs à la perturbation, \(t > \tau .\)
(b) Montrer qu’après avoir dérivé une fois par rapport au temps l’équation linéaire sur\(\delta p\), il est possible d’éliminer\(\delta \vec v\). En déduire une équation d’évolution portant seulement sur \(\delta \rho .\)
3.3 Modes propres d’un gaz superfluide homogène
On suppose que le gaz superfluide à l’équilibre remplit tout l’espace avec une densité uniforme\({\rho _e}\). On souhaite déterminer les modes propres de l’équation d’évolution sur \(\delta \rho .\) Pour cela, on considère une solution de la forme:
$\delta \rho \left( \vec{r},~t \right)=A\text{ }\!\!~\!\!\text{ cos }\!\!~\!\!\text{ }\left( \vec{k}\cdot \vec{r}-\Omega t \right)$ (8)
où \(A,\) \(\Omega \) et le vecteur réel \(\vec k\) sont des constantes.
(a) Vérifier que la forme proposée convient à condition que \({\Omega ^2}\) soit une certaine fonction de \(\vec k\) que l’on précisera.
(b) Interpréter physiquement le résultat obtenu. On précisera à quelle vitesse une perturbation appliquée en un point du gaz superfluide se propage. Comment cette vitesse dépend‐elle de la forme de la perturbation?
(c) Après envoi d’une impulsion laser dans un gaz superfluide d’atomes de sodium \({}^{23}Na\)initialement au repos, l’équipe du MIT a constaté la propagation d’ondes de densité à une vitesse de\(1 \times {10^{ - 2}}m/s\). On suppose que le gaz était initialement quasi homogène avec une densité moyenne de \(3,8 \times {10^{20}}\) atomes par\({m^3}\). Comparer la vitesse mesurée à la prédiction théorique.
3.4 Modes propres dans un piège harmonique
Le gaz superfluide est initialement au repos dans un potentiel de piégeage \({\Phi _e}\left( {{r^ \to }} \right)\) harmonique anisotrope:
\({\Phi _e}\left( {\vec r} \right) = \frac{1}{2}m\mathop \sum \limits_{\alpha = x,y,z}^{} \omega _\alpha ^2r_\alpha ^2\) (9)
où \({r_\alpha }\) est la composante du vecteur position $\vec{r}$ sur le vecteur de base \(\vec e,\) \(\alpha = x,\) \(y,\) \(z\). On excite le gaz en modifiant faiblement l’un des paramètres \({\omega _a}\) pendant une durée\(\tau \). On admet qu’on excite ainsi des modes propres du gaz superfluide de la forme
$\delta \rho \left( {{r}^{\to }},~t \right)=\left[ B+\underset{\alpha =x,y,z}{\overset{{}}{\mathop \sum }}\,{{A}_{\alpha }}r_{\alpha }^{2} \right]\text{ }\!\!~\!\!\text{ cos }\!\!~\!\!\text{ }\Omega t$ (10)
où \(B\), les coefficients \({A_\alpha }\) et la pulsation \(\Omega \) sont des constantes.
(a) Exprimer la densité \({\rho _e}\left( {\vec r} \right)\) du gaz à l’équilibre à une constante additive près, en fonction des \({r_\alpha },\) \({\omega _\alpha }\), de \(m\) et de \(\gamma .\)
(b) À quelles conditions sur les coefficients \({A_\alpha }\) la forme proposée pour \(\delta \rho \) satisfait‐elle l’équation d’évolution obtenue à la partie 3.2? Ces conditions peuvent‐elles être satisfaites quelle que soit la pulsation \(\Omega \)?
(d) En déduire que \(X = \frac{{{\Omega ^2}}}{{\omega _z^2}}\) doit vérifier une équation du second degré que l’on précisera. Donner les solutions de cette équation en fonction du paramètre \(\eta = \frac{{{\omega _x}}}{{{\omega _z}}}.\)
(e) Dans l’expérience du MIT, le paramètre \(\eta \) vaut 13,60. Donner les pulsations propres des deux modes excités, ainsi que les valeurs correspondantes du rapport \(\frac{{{A_z}}}{{{A_x}}}.\) Expérimentalement, on constate que le gaz bat à la pulsation \(1,569 \pm 0,004{\omega _z}.\) Quel est l’écart relatif entre la prédiction théorique et le résultat expérimental?
(f) À quelle fréquence devrait battre le gaz en l’absence d’interaction entre les particules? Qu’en concluez‐vous?
