Dynamical structure factor with a pumping approach on a trapped-ion quantum computer

Technische Samenvatting: Dynamische Structuurfactor met een Pumping-benadering op een Trapped-Ion Kwantumcomputer

Probleemstelling
De dynamische structuurfactor (DSF), aangeduid als $S(q, \omega)$, is een cruciale grootheid in de gecondenseerde materie natuurkunde die de respons van een materiaal op neutronenverstrooiing beschrijft. Het wordt gedefinieerd als de ruimte-tijd Fourier-transformatie van dynamische correlaties $\langle O(x, t)O(0, 0) \rangle$. Het berekenen van de DSF is uitdagend voor klassieke computers omdat het zowel het vinden van een representatie van de grondtoestand of thermische toestand als het simuleren van Hamiltoniaanse dynamica over de tijd vereist, een taak die over het algemeen als BQP-hard wordt beschouwd. Hoewel kwantumcomputers een pad bieden om eindtemperatuur-toestanden voor te bereiden en efficiënte Hamiltoniaanse dynamica uit te voeren, vertrouwen bestaande kwantumalgoritmen voor het berekenen van $S(q, \omega)$ doorgaans op het berekenen van dynamische correlaties op veel discrete tijdstippen en het vervolgens uitvoeren van een Fourier-transformatie. Deze aanpak brengt aanzienlijke shot-overhead met zich mee, vooral wanneer slechts enkele specifieke frequentiewaarden van belang zijn.

Methodologie: De Pumping-benadering
De auteurs introduceren een "pumping-benadering" die direct specifieke frequentiewaarden $\omega$ target zonder de noodzaak om tijd-domein correlaties te berekenen en een Fourier-transformatie uit te voeren. De kern van de methode omvat het uit evenwicht drijven van het systeem met behulp van een tijdsafhankelijke Hamiltoniaan die wordt verstoord door een bronterm die oscilleert op de doel frequentie.

1. Hamiltoniaan Opzet: Gegeven een systeem-Hamiltoniaan $H$ en een observeerbare $O$, definiëren de auteurs een tijdsafhankelijke Hamiltoniaan:
   $$H(t) = H - \epsilon O(0) \sin(\omega(T - t))$$
   waarbij $\epsilon$ een koppelingsconstante is en $T$ de uiteindelijke evolutietijd.
2. Meetprotocol: Het systeem wordt geëvolueerd vanuit een initiële evenwichtstoestand (bijv. grondtoestand) met behulp van deze tijdsafhankelijke Hamiltoniaan. De verwachtingswaarde van een momentum-ruimte operator $\hat{O}_q$ wordt gemeten op tijdstip $T$.
3. Extractie van DSF: Door het systeem te meten met koppeling $+\epsilon$ en $-\epsilon$, wordt de DSF geëxtraheerd via de relatie:
   $$S(q, \omega) \approx \frac{\langle \hat{O}_q(T) \rangle_\epsilon - \langle \hat{O}_q(T) \rangle_{-\epsilon}}{\epsilon(1 - e^{-\beta\omega})}$$
   Deze formule is afgeleid van een eerste-orde expansie in $\epsilon, waarbij gebruik wordt gemaakt van de gedetailleerde balans-relatie voor thermische toestanden.
4. Implementatie details:
   • Toestandsvoorbereiding: Voor het 1D Heisenberg-model wordt de grondtoestand voorbereid met een adiabatisch-geïnspireerde ansatz. Het systeem wordt geïnitialiseerd in een Néel-toestand en geëvolueerd via een lineaire interpolatie tussen een gestaffeld-veld Hamiltoniaan en de Heisenberg Hamiltoniaan. De parameters worden klassiek geoptimaliseerd met statevector-simulaties op een keten van 20 sites en vervolgens toegepast op de hardware.
   • Tijds-evolutie: De tijds-evolutie operator wordt benaderd met behulp van een Trotter-decompositie. Het circuit bestaat uit opeenvolgende single-qubit en two-qubit gates (specifiek TK2 gates op de Quantinuum hardware).
   • Regularisatie: Omdat $T \to \infty$ niet haalbaar is, wordt een eindige $T$ gebruikt, wat effectief de ware DSF convolueert met een sinc-functie, wat het resultaat gladstrijkt maar potentieel spurieuze kenmerken introduceert.

Belangrijkste Bijdragen

• Directe Frequentie Targeting: De primaire bijdrage is een methode die kwantumcomputers in staat stelt om direct specifieke willekeurige frequenties $\omega$ te targeten, waardoor de computationele kosten van het samplen van het volledige tijd-domein worden vermeden.
• Shot Efficiëntie: De auteurs demonstreren dat voor scenario's waarin slechts enkele frequentiewaarden nodig zijn, deze aanpak een significante reductie in shot-overhead biedt vergeleken met methoden die tijd-domein correlatiemetingen vereisen.
• Hybride Implementatie: Het werk presenteert een volledige implementatie op de Quantinuum Reimei trapped-ion kwantumcomputer, geïntegreerd binnen een Fugaku–Reimei hybride kwantum–HPC omgeving. De klassieke optimalisatie van de toestand-voorbereidingsparameters wordt uitgevoerd op de Fugaku supercomputer, terwijl de kwantumcircuits worden uitgevoerd op Reimei.

Resultaten
De auteurs valideerden de methode met twee gevallen:

1. 1D Heisenberg Model (20 sites): De DSF werd berekend voor een keten van 20 sites en vergeleken met ruisloze simulaties. De hardware-resultaten vertoonden een uitstekende overeenkomst met de ruisloze simulaties, behaald zonder enige technieken voor ruis-mitigatie.
2. Koper Sulfaat Kristal ($CuSO_4$): Het 1D Heisenberg-model beschrijft bekendelijk de magnetische excitaties van koper sulfaat kristallen. De berekende DSF werd vergeleken met experimentele neutronenverstrooiingsdata. De resultaten van de kwantumhardware kwamen zeer goed overeen met de experimentele data.

Specifieke parameters die werden gebruikt omvatten een systeemgrootte van $L=20$, een Trotter stapgrootte van $dt=0.24$, een totale evolutietijd van $T=2.4$, en een koppeling van $\epsilon=2$. De grondtoestand-voorbereiding maakte gebruik van een enkele ronde ($R=1$) van de geoptimaliseerde ansatz. De studie merkte op dat het gebruik van een enkele laag van de ansatz feitelijk een kleinere afwijking opleverde van de exacte grondtoestand resultaten dan het gebruik van twee lagen, waarschijnlijk door de reductie in gate-aantal.

Betekenis en Claims
Het artikel claimt dat deze pumping-aanpak een praktische en efficiënte weg biedt voor het berekenen van dynamische structuurfactoren op nabije-term kwantumhardware. Door direct een "speelgoedmodel versie van een neutronenverstrooiing experiment" te simuleren, omzeilt de methode de noodzaak van complexe real-time correlatieberekeningen.

De auteurs benadrukken de potentiële relevantie van dit protocol buiten neutronenverstrooiing, specifiek voor Nuclear Magnetic Resonance (NMR) simulaties. In NMR zijn spectra vaak scherp gelokaliseerd rond specifieke frequenties, en deze methode zou aanzienlijke runtime-besparingen kunnen bieden door vaste frequenties te samplen in plaats van lange tijd-domein vervalprocessen.

Wat betreft schaalbaarheid zijn de auteurs bescheiden over de huidige kwantumvoorsprong. Ze erkennen dat de toestand-voorbereiding in dit werk steunde op klassiek geoptimaliseerde ansatzes afgeleid van exacte statevector-simulaties, een regime waar geen kwantumvoorsprong bestaat. Ze stellen echter voor dat het combineren van hun DSF-algoritme met schaalbare klassieke toestand-voorbereidingsmethoden (zoals Sparse Pauli Dynamics of Matrix Product States via DMRG) binnen een hybride framework uiteindelijk kan leiden tot een regime waar de gezamenlijke taak van toestand-voorbereiding en structuurfactor-berekening uitdagend wordt voor klassieke computers alleen.