Phonon number relaxation in a 3D superfluid with a concave acoustic branch

Dit artikel beschrijft hoe een superfluïde gas van fononen in drie dimensies via vijf-fononbotsingen, die veel trager zijn dan vier-fononbotsingen, langzaam evolueert naar volledige thermodynamische en chemische evenwichtstoestand, waarbij de fugaciteit en entropie specifieke tijdsafhankelijke wetten volgen.

De kern

De Dans van de Geluidsgolven: Waarom het evenwicht zo lang duurt

Stel je een supergeleidende vloeistof voor (een "superfluide") als een perfect georganiseerd dansfeest. In dit feestje zijn de gasten geen mensen, maar fononen: de kleinste deeltjes van geluid die door het materiaal reizen. Normaal gesproken willen deze gasten zich gedragen volgens de regels van de thermodynamica: ze willen een evenwichtige temperatuur bereiken en zich willekeurig verspreiden, net als mensen op een drukke dansvloer.

Maar in dit specifieke artikel kijken we naar een heel speciale situatie: een superfluide waar de "geluidssnelheid" op een rare manier afhangt van de energie. De wetenschappers Yvan Castin en Mariia Tsimokha hebben ontdekt dat in zo'n systeem het bereiken van rust (equilibrium) een enorm langzaam en lastig proces is.

Hier is hoe het werkt, stap voor stap:

1. Het Probleem: De "Verboden" Danspas

In de meeste superfluiden kunnen fononen makkelijk met elkaar botsen. Stel je voor dat twee gasten (fononen) samenkomen en direct in drie nieuwe gasten kunnen veranderen, of andersom. Dit is een drie-phonon botsing ($1 \leftrightarrow 2$). Dit gaat heel snel en zorgt ervoor dat het systeem snel tot rust komt.

Maar in dit specifieke geval (waar de "acoustische tak" hol is, ofwel concave) is deze danspas verboden. De wetten van behoud van energie en impuls staan het niet toe dat twee fononen simpelweg in drie veranderen. Het is alsof je op een dansfeest probeert twee mensen te laten samensmelten tot drie, maar de dansvloer (de natuurwetten) zegt: "Nee, dat kan niet, de geometrie klopt niet."

2. De Tussenstap: De Vier-Persoons Dans

Omdat de snelle drie-persoons botsing verboden is, moeten de fononen een langzamere route nemen. Ze proberen eerst een vier-phonon botsing ($2 \leftrightarrow 2$). Twee gasten botsen en gaan weg als twee andere gasten.

• Wat dit doet: Dit zorgt ervoor dat de gasten zich wel verdelen over de dansvloer (thermisch evenwicht), maar het verandert niet het totale aantal gasten.
• Het probleem: Het systeem heeft nu een "potentiaal" (een soort druk om meer gasten te hebben of minder), maar omdat het aantal gasten vastzit, kan het niet volledig tot rust komen. Het is alsof je een kamer hebt waar de temperatuur gelijk is, maar de luchtvochtigheid niet. Het is een half-weg evenwicht.

3. De Echte Oplossing: De Vijf-Persoons Dans

Om het systeem echt in evenwicht te krijgen (zodat het totale aantal fononen kan veranderen en het chemisch evenwicht bereikt), moeten ze wachten op een vijf-phonon botsing ($2 \leftrightarrow 3$).

• Hoe het werkt: Twee fononen botsen en worden drie, of drie worden twee. Dit is de enige manier om het totale aantal fononen te veranderen en het systeem volledig te laten "ademen".
• De snelheid: Dit is extreem zeldzaam. Het is alsof je wacht tot vijf mensen op een dansfeest tegelijkertijd in een perfecte cirkel staan en dan plotseling van vorm veranderen.
• De tijd: De wetenschappers hebben berekend dat dit proces zo langzaam is dat de tijd die het kost evenredig is met $T^{-9}$ (waarbij $T$ de temperatuur is). Bij heel lage temperaturen duurt dit eeuwen in menselijke termen, of in elk geval veel langer dan de vier-phonon botsingen.

4. De Reis van de "Fugaciteit" (De Druk)

De auteurs volgen de "fugaciteit" ($z_\phi$). Je kunt dit zien als een drukmeter voor het aantal fononen.

• Start: Het systeem begint met een druk van 0 (zeer weinig fononen, een lege dansvloer).
• Einde: Het systeem wil naar een druk van 1 (volledig evenwicht).
• De reis:
  • Korte tijd: Aan het begin groeit de druk heel snel, maar niet lineair. Het volgt een vreemde wet: $t^{4/5}$. Alsof de dansvloer eerst snel volloopt met een paar groepjes.
  • Lange tijd: Als het systeem bijna vol is, vertraagt het proces enorm en volgt het een exponentiële afname naar het eindpunt. Het is alsof de laatste paar gasten heel langzaam de juiste plek zoeken.

5. De Entropie (De Chaos)

Een mooi resultaat van hun berekening is dat de entropie (de chaos of wanorde in het systeem) altijd toeneemt, zoals de tweede wet van de thermodynamica voorspelt. Maar de snelheid waarmee de chaos toeneemt, hangt af van het kwadraat van de snelheid waarmee de druk verandert.

• Vergelijking: Stel je voor dat je een auto remt. Hoe harder je remt (snelheidswijziging), hoe meer hitte (entropie) er vrijkomt. In dit geval is de hitte evenredig met het kwadraat van de remkracht. Dit bevestigt een oude voorspelling van de legendarische natuurkundige Landau.

Waarom is dit belangrijk?

Dit is niet alleen maar droge theorie. De auteurs zeggen dat we dit in de praktijk kunnen testen:

1. Koude atoomgassen: In laboratoria met ultra-koude atomen (die gedragen zich als een superfluide) kunnen we de druk zo instellen dat we deze holle golfkromme krijgen.
2. Helium-4: Als we vloeibare helium onder zeer hoge druk zetten, zou dit ook kunnen gebeuren.

Samenvatting in één zin

In een speciaal soort superfluide kunnen geluidsgolven niet snel met elkaar botsen om in evenwicht te komen; ze moeten wachten op een extreem zeldzame gebeurtenis waarbij vijf golven tegelijk interageren, wat leidt tot een heel langzaam, maar wiskundig prachtig voorspelbaar proces van rust en orde.

Technische samenvatting

Probleemstelling

Het artikel onderzoekt de collisionele evolutie naar thermodynamisch evenwicht van een ruimtelijk homogeen en isotroop gas van fononen in een driedimensionaal superfluïde. Het specifieke scenario is een superfluïde met een concave akoestische excitatie-tak (gekenmerkt door een krommingsparameter $\gamma < 0$) bij een temperatuur $T$ die niet nul is, maar willekeurig laag.

De centrale uitdaging is het begrijpen van hoe dit systeem relaxeert naar een volledig thermisch en chemisch evenwicht (waarbij het chemische potentieel van de fononen $\mu_\phi = 0$ is).

• Drie-fonon processen ($1\phi \leftrightarrow 2\phi$): Deze zijn, vanwege de wetten van behoud van energie en impuls, verboden in een concave tak. Ze kunnen dus niet leiden tot evenwicht.
• Vier-fonon processen ($2\phi \leftrightarrow 2\phi$): Deze processen (Landau en Khalatnikov) zorgen wel voor thermische relaxatie van de golfvectorverdeling naar een Bose-verdeling, maar ze behouden het totale aantal fononen $N_\phi$. Hierdoor bereiken ze slechts een partieel evenwicht met een niet-nul chemisch potentieel ($\mu_\phi \neq 0$) op een tijdschaal van $\propto T^{-7}$.
• Het probleem: De relaxatie naar volledig thermochemisch evenwicht ($\mu_\phi \to 0$) vereist processen die het aantal fononen niet behouden. Dit wordt verzorgd door de veel langzamere vijf-fonon processen ($2\phi \leftrightarrow 3\phi$), die een karakteristieke tijdschaal van $\propto T^{-9}$ hebben.

Het doel van dit werk is om de relaxatiewet van het chemisch potentieel (of de fugaciteit $z_\phi$) kwantitatief te bepalen voor het geval van een concave tak, een studie die kwalitatief begon met Khalatnikov in 1950 maar nu volledig kwantitatief wordt gemaakt.

Methodologie

De auteurs gebruiken een combinatie van kinetische theorie en kwantumhydrodynamica:

1. Kinetische Vergelijkingen: Ze stellen de Boltzmann-vergelijkingen op voor de bezettingsgetallen $n_q$ van de fononmodi, specifiek voor de vijf-fonon botsingen ($2\phi \leftrightarrow 3\phi$). De verandering in het totale aantal fononen $(d/dt)N_\phi$ wordt uitgedrukt als een functioneel van de verdelingen.
2. Berekening van de Botsingsamplitude: Een cruciaal onderdeel is de expliciete berekening van de overgangsamplitude $A_{2\phi \to 3\phi}$ voor lage temperaturen.
   • Ze gebruiken de kwantumhydrodynamica van Landau en Khalatnikov (kubische interactie Hamiltoniaan $\hat{H}_3$).
   • Ze passen de veralgemeende regel van Fermi toe, waarbij ze zorgvuldig rekening houden met de subdominante orde en de aanwezigheid van een thermisch bad (in plaats van een vacuüm). Ze tonen aan dat de amplitude op de energieschil ($E_i = E_f$) niet beïnvloed wordt door de thermische bezetting van de interne modi.
   • De berekening omvat zes Feynman-diagrammen (derde orde in de verstoring) en maakt gebruik van schaaltransformaties en complexe variabelen om de hoekafhankelijkheid te vereenvoudigen.
3. Thermodynamische Koppeling: Omdat het systeem geïsoleerd is, is de totale energie $E_\phi$ behouden. Dit biedt een tweede vergelijking die de temperatuur $T(t)$ koppelt aan de fugaciteit $z_\phi(t)$.
4. Numerieke Integratie: De auteurs berekenen een dimensieloze coëfficiënt $C_\phi(z_\phi)$ die voortkomt uit de integratie over alle golfvectoren en hoeken. Deze wordt numeriek geëvalueerd met geavanceerde integratietechnieken (inclusief analytische integratie over een hyperradius om singulariteiten te vermijden).

Belangrijkste Bijdragen

• Kwantitatieve Afleiding: Voor het eerst wordt de relaxatiesnelheid van het chemisch potentieel in een superfluïde met een concave tak volledig kwantitatief afgeleid, inclusief alle numerieke factoren.
• Expliciete Amplitude: De auteurs leveren een expliciete uitdrukking voor de amplitude van de vijf-fonon botsing ($2\phi \to 3\phi$) bij lage temperaturen, wat een ontbrekende schakel was in de theorie.
• Gedetailleerde Kinetiek: Ze beschrijven het volledige traject van de fugaciteit $z_\phi$ van het niet-gedegenereerde regime ($z_\phi \to 0$) tot het volledige evenwicht ($z_\phi \to 1$).
• Entropie-analyse: Ze tonen aan dat de entropietoename consistent is met de tweede wet van de thermodynamica en dat de snelheid van entropieverandering evenredig is met het kwadraat van de snelheid van verandering van de fugaciteit, in overeenstemming met Landau's voorspelling voor langzame adiabatische transformaties.

Resultaten

De analyse leidt tot de volgende specifieke resultaten voor de evolutie van de fugaciteit $z_\phi(t)$:

1. Kinetische Vergelijking: De tijdsafgeleide van de fugaciteit wordt beschreven door een differentiaalvergelijking van de vorm:
   $$\frac{1}{\Gamma^{eq}_\phi} \frac{d z_\phi}{dt} = F(z_\phi)$$
   waarbij $\Gamma^{eq}_\phi \propto T^9$ de relaxatiesnelheid nabij evenwicht is en $F(z_\phi)$ een functie is die afhangt van de fugaciteit en de numeriek berekende coëfficiënt $C_\phi(z_\phi)$.

2. Korte Tijden (Niet-gedegenereerd regime):
   Wanneer het systeem start vanuit een zeer lage dichtheid ($z_\phi \to 0$), volgt de fugaciteit een niet-geheel machtswet:
   $$z_\phi(t) \propto t^{4/5}$$
   Dit is een opvallend resultaat dat afwijkt van de gebruikelijke lineaire of exponentiële benaderingen.

3. Lange Tijden (Evenwichtsnabijheid):
   Naarmate het systeem het volledige evenwicht nadert ($z_\phi \to 1$), gedraagt de relaxatie zich exponentieel:
   $$z_\phi(t) \approx 1 - \exp(B - \Gamma^{eq}_\phi t)$$
   De relaxatiesnelheid is hier bepaald door $\Gamma^{eq}_\phi \propto T^9$.

4. Entropie: De snelheid van entropieverandering $(d/dt)S_\phi$ is altijd positief en asymptotisch evenredig met $[(d/dt)z_\phi]^2$.

5. Numerieke Coëfficiënten: De auteurs berekenen de waarden voor de coëfficiënt $C_\phi$ bij de grenzen:

   • $C_\phi(0) \approx 3.528 \times 10^4$
   • $C_\phi(1) \approx 5.422 \times 10^4$

Betekenis en Toepassingsmogelijkheden

Dit werk sluit een theoretisch onderzoek af dat in 1950 door Khalatnikov werd gestart. De resultaten zijn niet louter academisch; ze bieden voorspellingen die experimenteel getest kunnen worden in:

• Koude Fermi-gassen: Specifiek aan de BCS-kant van de BEC-BCS crossover, waar de collectieve excitatie-tak concave is ($\gamma < 0$).
• Superfluïde Helium-4: Bij voldoende hoge druk, waar de dispersierelatie van fononen ook concave wordt.

De studie bevestigt dat in superfluïde systemen met een concave dispersierelatie, de relaxatie naar chemisch evenwicht wordt gedomineerd door zeer zeldzame, maar cruciale, vijf-fonon botsingen, en levert een nauwkeurige wiskundige beschrijving van dit langzame proces.