Suppressed correlation-spreading in a one-dimensional Bose-Hubbard model with strong interactions

Dit artikel onderzoekt onderdrukte correlatieverspreiding in een één-dimensionaal Bose-Hubbard-model met sterke interacties, waarbij wordt aangetoond dat de dynamiek wordt gedomineerd door langzame domeinwand-excitaties en dat een parabolische val de correlatievoortplanting zelfs bij zwakke valsterktes verder onderdrukt.

De kern

Titel: De "Bevroren" Dans van Atomen: Waarom sommige quantum-systemen nooit tot rust komen

Stel je voor dat je een grote zaal hebt vol met dansers (atomen). Normaal gesproken, als je de muziek aanzet, beginnen ze te dansen, te wisselen van partner en uiteindelijk een perfecte, chaotische dans te vormen waarbij iedereen evenveel beweegt. Dit noemen wetenschappers "thermisch evenwicht" of "ergodisch gedrag". Alles mengt zich en de oorspronkelijke volgorde is vergeten.

Maar wat als je een speciale dansstijl kiest waarbij de dansers in een streng patroon staan: twee dansers op de ene plek, niemand op de volgende, weer twee, weer niemand? En wat als de dansers zo bang zijn om aan te raken (ze stoten elkaar af) dat ze bijna niet kunnen bewegen?

Dat is precies wat deze wetenschappers hebben onderzocht. Ze keken naar een heel specifiek quantum-systeem (het Bose-Hubbard-model) en ontdekten iets verrassends: in sommige gevallen stoppen de atomen met het vergeten van hun oorspronkelijke dans. Ze blijven "gevangen" in hun oude patroon.

Hier is hoe het werkt, vertaald in alledaagse taal:

1. De Dans van de "Dubbele" en de "Lege"

In hun experiment beginnen ze met een rij atomen die eruitziet als: 2 - 0 - 2 - 0 - 2 - 0.

• 2 betekent: twee atomen op dezelfde plek.
• 0 betekent: een lege plek.

Normaal gesproken zou een atoom van plek 2 naar plek 0 springen. Maar omdat de atomen elkaar zo sterk afstoten (sterke interactie), is het bijna onmogelijk voor één atoom om alleen te springen. Het is alsof ze in een drukke menigte zitten waar je niet doorheen kunt lopen zonder iemand te raken.

In plaats van één atoom te laten springen, gebeurt er iets slim: twee atomen springen samen als een koppel naar de lege plek ernaast.

• Van 2-0 wordt 0-2.
• Dit koppel (een "doublon") en de lege plek (een "holon") wisselen van plaats.

2. De Muur die Loopt (De "Domain Wall")

Stel je voor dat je een lange rij hebt met deze paren. Als je aan het ene einde begint met het wisselen van een koppel, ontstaat er een muur tussen het oude patroon en het nieuwe patroon.

• Links van de muur: 2-0-2-0 (het oude patroon).
• Rechts van de muur: 0-2-0-2 (het nieuwe patroon).

Deze muur beweegt door de rij als een golf. In een perfecte, lege ruimte (zonder obstakels) zou deze muur snel door de hele rij lopen, totdat alles is omgedraaid. Dit is de "relaxatie": het systeem verandert van vorm.

3. De Valstrik: De "Parabolische Loopt"

Nu komt het spannende deel. De onderzoekers voegden een "valstrik" toe: een zachte, ronde helling in het midden van de zaal. De atomen aan de randen zitten in een soort "diepe kuil" en voelen zich daar veilig en vastgeklemd.

Het verrassende resultaat:
Zelfs als deze valstrik heel zacht is (zwakker dan de kracht waarmee de atomen normaal kunnen springen), bevriest het systeem aan de randen.

• De "muur" die door de rij loopt, botst tegen de randen en kan niet verder.
• De atomen aan de randen kunnen niet bewegen.
• Hierdoor kan de informatie (de correlatie) zich niet verspreiden door het hele systeem. Het systeem "vergeet" niet wat het was. Het blijft in een staat van niet-ergodisch gedrag: het is niet volledig gemengd.

4. De Analogie van de Spiegels en de IJsklomp

Om dit beter te begrijpen, gebruiken de auteurs een slimme truc. Ze zeggen: "Laten we dit systeem niet zien als atomen, maar als een rij magneetjes (spins)."

• 2-0 wordt een pijl die naar boven wijst (↑).
• 0-2 wordt een pijl die naar beneden wijst (↓).

In dit nieuwe beeld gedraagt het systeem zich precies als een IJsklomp met een magneetveld. De "muur" die we eerder zagen, is nu gewoon een omgekeerde magneet die door de rij loopt.

• De snelheid waarmee deze magneet omklapt, hangt af van hoe sterk de atomen elkaar afstoten. Hoe sterker de afstoting, hoe trager de magneet beweegt.
• De onderzoekers ontdekten dat de snelheid waarmee de "informatie" door het systeem reist, precies overeenkomt met de snelheid van een golf in deze magneetrij.

Waarom is dit belangrijk?

In de wereld van quantumcomputers en nieuwe materialen willen we vaak weten: Hoe snel verspreidt informatie zich?

• Normaal: Informatie verspreidt zich snel als een lichtstraal (zoals een rimpeling in een meer).
• In dit experiment: Informatie blijft lokaal. Het verspreidt zich nauwelijks, vooral als er een zachte valstrik is.

Dit betekent dat we systemen kunnen bouwen die niet tot rust komen, maar hun quantum-informatie behouden. Het is alsof je een danspartij hebt waar de dansers aan de randen vastzitten in hun stoelen, waardoor de rest van de zaal ook niet meer kan dansen.

Kort samengevat:
De onderzoekers laten zien dat als atomen elkaar heel hard afstoten en er een zachte helling is, ze in een "sluimerende" toestand kunnen blijven. Ze wisselen wel af en toe van plaats (als paren), maar de grote verandering (thermisch evenwicht) wordt geblokkeerd. Het systeem "vergeet" zijn verleden niet, en dat is een heel zeldzaam en interessant fenomeen in de quantumwereld.

Technische samenvatting

Titel: Onderdrukte correlatieverspreiding in een eendimensionaal Bose-Hubbard-model met sterke interacties

Auteurs: Jose Carlos Pelayo en Ippei Danshita (Kindai University, Japan)

---

1. Probleemstelling en Context

Het artikel onderzoekt de dynamica van niet-ergodisch gedrag in geïsoleerde kwantumveeldeeltjessystemen, specifiek in het kader van het eendimensionale (1D) Bose-Hubbard-model.

• Ergoditeit en ETH: In de meeste systemen zorgt de Eigenstate Thermalization Hypothesis (ETH) ervoor dat systemen thermisch evenwicht bereiken. Echter, bij sterke interacties of specifieke beperkingen kan deze hypothesis worden geschonden, wat leidt tot vertraagde relaxatie of localisatie.
• Specifiek scenario: De auteurs focussen op de tijdevolutie van een systeem dat start in een dubbel bezette dichtheidsgolf-toestand (doubly occupied density-wave state, $|2, 0, 2, 0, \dots\rangle$).
• Vraagstelling: Waar eerdere studies aantoonden dat een enkel bezette dichtheidsgolf snel relaxeert, en dat een dubbel bezette golf in een val kan leiden tot Many-Body Localization (MBL), is het onduidelijk hoe de verspreiding van kwantumcorrelaties zich gedraagt in een systeem met sterke interacties ($U/J \gg 1$) en een zwakke valpotentiaal. De centrale vraag is of de correlatieverspreiding wordt onderdrukt in dit specifieke regime.

2. Methodologie

De auteurs gebruiken een combinatie van numerieke simulaties en analytische effectieve modellen:

• Model: Een 1D Bose-Hubbard keten met open randvoorwaarden en een parabolische valpotentiaal. De Hamiltoniaan bevat tunneling ($J$), on-site interactie ($U$) en een valterm ($\Omega$).
• Numerieke Simulatie:
  • Gebruik van Matrix Product States (MPS) met het Time-Evolving Block Decimation (TEBD) algoritme.
  • Parameters: Maximale bandbreedte $\chi_{max} = 1600$, maximale bezetting per site $n_{max} = 4$.
  • Grootte: $M=16$ sites (voor dynamiek) en $M=12$ (voor langetijdsanalyse).
  • Interactiestrategie: Sterke interacties ($U/J = 40$) en een zwakke val ($\Omega/J = 0.02$).
• Analytische Benadering: Afleiding van een effectieve Hamiltoniaan via de Schrieffer-Wolff-transformatie om het systeem te mappen naar een effectief spin-model in de limiet van sterke interacties.

3. Belangrijkste Bijdragen en Resultaten

A. Dynamica van de Site-bezetting (Site Occupation Dynamics)

• Mechanisme: Bij sterke interacties ($U \gg J$) is het direct hoppen van een enkel deeltje energetisch ongunstig. In plaats daarvan domineert paar-tunneling (doublon-holon uitwisseling). Een doublon ($|2\rangle$) en een holon ($|0\rangle$) wisselen van positie: $|2, 0\rangle \to |0, 2\rangle$.
• Domeinwandexcitatie: Deze uitwisseling creëert een "domeinwand" die zich voortplant door het systeem. In afwezigheid van een val ($\Omega=0$) beweegt deze wand van de randen naar het midden en keert het systeem uiteindelijk om.
• Invloed van de Val: Zelfs bij een zeer zwakke val ($\Omega/J = 0.02$) wordt de dynamica aan de randen "bevroren". De energie-kost voor paar-tunneling aan de randen is te hoog, waardoor de bezetting daar vastzit. Dit onderdrukt de voortplanting van de domeinwand aanzienlijk.

B. Correlatieverspreiding

De auteurs analyseren drie soorten correlatiefuncties:

1. Een-deeltjescorrelatie ($C_{sp}$): Toont een sterke localisatie; correlaties vervallen snel buiten de naaste buren.
2. Paar-correlatie ($C_{pair}$): Toont iets meer verspreiding dan $C_{sp}$, maar blijft gedomineerd door naaste-buren-correlaties, wat wijst op localisatie.
3. Gekoppelde dichtheids-dichtheids correlatie ($C_{dd}$):
   • Toont een negatief gestreept patroon (negated staggered pattern) op grotere afstanden. Dit is een direct signaal van de voortplantende domeinwandexcitatie.
   • Onderdrukking: In aanwezigheid van een val wordt de verspreiding van deze correlaties sterk onderdrukt. De correlaties blijven beperkt tot kortere afstanden, zelfs op lange tijdschalen.

C. Tijdschalen en Relaxatie

• De tijdschaal voor de excitatie en voortplanting wordt bepaald door de inverse van de interactiestrakteerheid ($\sim U/J$).
• De effectieve tunnelingssterkte is $\tilde{J} = 2J^2/U$. Omdat $\tilde{J}$ omgekeerd evenredig is met $U$, wordt de relaxatie bij sterke interacties extreem traag.
• In de val ($\Omega > 0$) blijft het systeem niet-ergodisch op lange tijden (de correlatiewijdte $\sigma(t)$ neemt af of blijft klein), terwijl het zonder val ($\Omega=0$) uiteindelijk toch naar een relaxatie-regime neigt, zij het traag.

D. Effectief Transvers-veld Ising (TFI) Model

• De auteurs mappen het Bose-Hubbard-model in de limiet $U \gg J$ af op een antiferromagnetisch transvers-veld Ising-model.
• Mapping: De configuraties $|2, 0\rangle$ en $|0, 2\rangle$ worden gemapt op spin-up ($|\uparrow\rangle$) en spin-down ($|\downarrow\rangle$).
• Validatie: Het TFI-model vat de kwantitatieve dynamica van de site-bezetting en de dichtheids-correlaties uitstekend samen.
• Voorspellende kracht: De voortplantingssnelheid van de correlaties in het oorspronkelijke model komt overeen met de groepssnelheid ($v_g$) van de spin-golfexcitatie in het effectieve spin-model. Dit biedt een intuïtief inzicht in de onderliggende fysica.

4. Significatie en Conclusie

• Nieuw inzicht in niet-ergodisch gedrag: Het artikel toont aan dat zelfs bij interacties die niet extreem sterk zijn (maar wel $U \gg J$), de combinatie van een dubbel bezette initiële toestand en een zwakke valpotentiaal leidt tot een sterke onderdrukking van correlatieverspreiding.
• Rol van de Val: Een verrassend resultaat is dat een valpotentiaal die zwakker is dan de tunnelingssterkte ($\Omega < J$) toch voldoende is om de dynamica aan de randen te bevriezen en de thermalisatie te onderdrukken.
• Methodologische bijdrage: De succesvolle mapping naar het TFI-model biedt een krachtige analytische tool om complexe bosonische dynamica te begrijpen zonder zware numerieke simulaties, en verklaart de waargenomen tijdschalen en snelheden.
• Implicaties: Deze bevindingen zijn relevant voor het begrijpen van kwantuminformatie-voortplanting en thermalisatie in ultrakoude atoomsystemen, en suggereren dat Hilbert-ruimte-fragmentatie (Hilbert Space Fragmentation) een cruciale rol speelt in het voorkomen van thermalisatie in deze systemen.

Samenvattend demonstreert dit werk dat de dynamiek van een 1D Bose-Hubbard-systeem met sterke interacties en een dubbel bezette starttoestand wordt gedomineerd door langzame, domeinwand-gedreven processen die gevoelig zijn voor randvoorwaarden en valpotentiaties, wat resulteert in een onderdrukte verspreiding van kwantumcorrelaties.