Relativistic Quantum Simulation under Periodic and Dirichlet Boundary Conditions: A First-Quantised Framework for Near-Term Devices

Technische Samenvatting: Relativistische Kwantumsimulatie onder Periodieke en Dirichlet-randvoorwaarden

Probleemstelling
Het simuleren van kwantumsystemen is computationeel veeleisend vanwege de exponentiële groei van complexiteit met de systeemgrootte, een uitdaging die wordt verergerd wanneer relativistische effecten moeten worden geïntegreerd om fysische verschijnselen nauwkeurig te beschrijven. Hoewel de perturbatietheorie een standaardmethode is voor het bestuderen van relativistische effecten buiten het niet-relativistische kader, blijft het implementeren van deze simulaties op nabijgelegen (near-term) kwantumtoestellen moeilijk. Bestaande methoden vertrouwen vaak op tweede kwantisatie, wat standaard is voor kwantumchemie, maar mogelijk tekortschiet in de intuïtieve afhandeling van systemen met niet-lineariteiten. Er is behoefte aan een praktisch, eerste-kwantisatie kader dat in staat is om relativistische grondtoestandsenergieën van systemen onder verschillende randvoorwaarden te schatten op nabijgelegen hardware.

Methodologie
De auteurs stellen een eerste-kwantisatie kader voor voor relativistische kwantumsimulatie dat de golffunctie discretiseert over een eindig rooster weergegeven door systeem-qubits. De kernmethodologie bestaat uit de volgende componenten:

1. Discretisatie en Translatie-operatoren:
   De continue 1D ruimtelijke golffunctie wordt gemapt naar een discreet rooster van $2^L$ punten met behulp van $L$ qubits, waarbij de positie $x$ wordt geschaald naar het domein $0 \le x < 1$. De gekwadrateerde impulsoperator $\hat{p}^2$ wordt benaderd met de eindige-verschil methode gebaseerd op kwantum translatie-operatoren. De translatie-operator $\hat{A}_x$ (kwantum-adder) verschuift de positietoestand $|j\rangle$ naar $|j+\delta x\rangle$, terwijl zijn conjunct $\hat{A}_x^\dagger$ (kwantum-subtractor) in de tegenovergestelde richting verschuift. De gekwadrateerde impuls wordt uitgedrukt als $\hat{p}^2 = -(mc L_m)^2 (\hat{A}^{(1)} - 2\hat{1})$, waarbij $L_m = \lambda_m / \delta x$ een dimensieloze parameter is die de gereduceerde Compton-golflengte omvat.

2. Perturbatieve Expansie van Kinetische Energie:
   De relativistische kinetische energie Hamiltoniaan $\hat{H}_{re}^K$ wordt afgeleid van de wortelvorm van de relativistische energie-impuls relatie. Deze wordt geëxpandeerd in een perturbatieve reeks van hogere-orde impulstermen ($\hat{p}^2, \hat{p}^4, \hat{p}^6, \dots$) die geldig is in het regime waar $\langle \hat{p}^2 \rangle \ll m^2 c^2$. De totale relativistische kinetische energie wordt benaderd door deze reeks te af te snijden, waarbij coëfficiënten worden opgenomen die zijn afgeleid van de Klein-Gordon vergelijking.

3. Afhandeling van Randvoorwaarden:
   Het kader adresseert twee onderscheidende randvoorwaarden:

• Periodieke Randvoorwaarden (PBC): De translatie-operatoren behouden unitariteit door het rooster te laten 'wrappen' ($\hat{A}_x |2^L-1\rangle = |0\rangle$). De verwachtingswaarden van even macht van de impuls worden uitgedrukt als combinaties van $\langle \hat{A}^{(l)} \rangle$.
• Dirichlet Randvoorwaarden (DBC): Hier 'wrapt' de golffunctie niet ($\psi(0) \neq \psi(1)$). De auteurs introduceren een correctie-operator $\hat{E}_0 = |2^L-1\rangle\langle 0| + |0\rangle\langle 2^L-1|$ om de rand te compenseren. De verwachtingswaarden voor DBC worden geformuleerd als de PBC-termen plus aanvullende termen die betrokken zijn bij $\hat{E}_0$ en gerelateerde operatoren ($\hat{E}_1, \hat{E}_2$), welke afhankelijk zijn van de waarschijnlijkheidsamplituden op de roostergrenzen.

4. Kwantumsimulatie Algoritme:
   De auteurs ontwerpen kwantumcircuits om de vereiste verwachtingswaarden te schatten:

• Kinetische Energie: Gecontroleerde-translatie poorten worden iteratief toegepast op een controle-qubit. Het meten van de controle-qubit in de Z-basis levert het reële deel van $\langle (\hat{A}_x^\dagger)^l \rangle$ op, wat overeenkomt met de noodzakelijke impuls-termen.
• Potentiële Energie: Een dichtheidsmatrix $\hat{\rho}_V$ die de potentiaal representeert wordt gebruikt in een circuit om $\langle \hat{V} \rangle$ te schatten via de waarschijnlijkheden van de systeemtoestand.
• Variational Optimalisatie: De totale energie wordt geschat door de som van de kinetische en potentiële energie verwachtingswaarden te minimaliseren over een verzameling variatiemparameters $\{\bar{\kappa}\}$ in een ansatz-toestand $|\psi(\{\bar{\kappa}\}, x)\rangle$.

Belangrijkste Bijdragen

• Eerste-kwantisatie Recept: Het artikel biedt een specifiek recept voor relativistische simulatie met behulp van eerste kwantisatie, wat een alternatief biedt voor de veelgebruikte tweede-kwantisatie methoden, met name gunstig voor systemen met niet-lineariteiten.
• Verenigd Randvoorwaarden Kader: Het presenteert een verenigde aanpak voor het afhandelen van zowel Periodieke als Dirichlet randvoorwaarden. Cruciaal is dat het aantoont dat de verkregen data voor PBC kan worden hergebruikt om DBC-energieën te schatten door specifieke randcorrectie-termen toe te voegen, wat de computationele overhead vermindert.
• Compatibiliteit met Nabijgelegen Toestellen: De methode maakt gebruik van standaard kwantumcircuit-primitieven (gecontroleerde-adders, Hadamard-poorten en metingen) en variationele optimalisatie, waardoor het geschikt is voor implementatie op nabijgelegen kwantumtoestellen.
• Expliciete Circuitontwerpen: De auteurs leveren gedetailleerde circuitontwerpen voor het berekenen van de verwachtingswaarden van de translatie-operatoren en de randcorrectie-termen die nodig zijn voor de perturbatieve expansie.

Resultaten en Claims
Het artikel beweert een "praktische route" te bieden voor het simuleren van relativistische effecten op nabijgelegen kwantumtoestellen. De theoretische resultaten tonen aan dat:

• Relativistische kinetische energietermen tot een specifieke orde van perturbatie kunnen worden gedecomposeerd in lineaire combinaties van verwachtingswaarden van translatie-operatoren ($\langle \hat{A}^{(l)} \rangle$) en rand-operatoren ($\langle \hat{E}_d \rangle$).
• De coëfficiënten voor deze expansies vooraf kunnen worden berekend op basis van fysische parameters ($m, c, \hbar, \delta x$).
• De methode de variationele optimalisatie van ansatz-toestanden mogelijk maakt om zowel niet-relativistische als relativistische grondtoestandsenergieën te schatten.

Betekenis
De auteurs positioneren dit werk als een fundamentele stap naar het ondersteunen van toekomstige ontwikkelingen in de kwantumfysica en kwantumchemie. Door de simulatie van relativistische effecten op nabijgelegen hardware mogelijk te maken, adresseert het kader een gat in de huidige capaciteiten van kwantumsimulatie. Het artikel benadrukt dat hoewel het verhogen van het aantal qubits de roosterresolutie verbetert, men voorzichtig moet zijn om ervoor te zorgen dat de dimensieloze parameter $L_m$ klein genoeg blijft zodat de perturbatieve expansie geldig blijft. Het werk claimt niet de volledige Dirac-vergelijking op te lossen of ultra-relativistische limieten te behandelen, maar richt zich specifiek op het regime waar relativistische correcties aan het niet-relativistische kader relevant zijn. Toekomstig werk wordt gesuggereerd om hogere-orde Laplacian-operatoren en nauwkeurigere ansatz-constructiemethoden te verkennen.