Impact of molecular orbital localization on quantum computational resources for Hamiltonian simulation: A benchmark study of hydrogen chain systems

Technische Samenvatting: Impact van de lokalisatie van moleculaire orbitalen op kwantumcomputationele middelen voor Hamiltoniaan-simulatie

Probleemstelling
Terwijl er aanzienlijk onderzoek is gedaan naar het optimaliseren van kwantumsimulatie-algoritmen (bijv. Trotterisatie, qubitisatie) en compressiestrategieën (bijv. tensorhypercontractie, QCELS), is er relatief weinig aandacht besteed aan de vraag hoe de keuze van een moleculaire orbitaalbasis de kwantumcomputationele middelen beïnvloedt. In de conventionele kwantumchemie zijn canonieke moleculaire orbitalen (CMO's), afgeleid van Hartree–Fock (HF) theorie, gedelokaliseerd, wat leidt tot "all-to-all" koppelingen in de gequantiseerde Hamiltoniaan. Deze delokalisatie vereist een groot aantal kwantumpoorten voor simulatie. Daartegenover staan gelokaliseerde moleculaire orbitalen (LMO's), die breed worden gebruikt in klassieke elektronische structuurmethoden om ruimtelijke lokaliteit uit te buiten, waarbij interacties tussen verre orbitalen verwaarloosbaar zijn. De centrale vraag die in deze studie wordt behandeld, is of orbitaal-lokalisatie op vergelijkbare wijze de complexiteit en de benodigde middelen voor kwantum Hamiltoniaan-simulatie kan verminderen, specifiek voor op Trotter-decompositie gebaseerde benaderingen.

Methodologie
De auteurs voerden een benchmarkstudie uit met ééndimensionale lineaire waterstofketensystemen ($H_n$) met $n$ variërend van 8 tot 100 atomen, bij een interatomaire afstand van 1,0 Å. De studie vergeleek twee orbitaalrepresentaties:

1. CMO's: Standaard Hartree–Fock canonieke orbitalen.
2. LMO's: Pipek–Mezey gelokaliseerde orbitalen, geherordend volgens hun relatieve ruimtelijke posities.

De Hamiltoniaan werd geconstrueerd met behulp van de Jordan–Wigner (JW) transformatie. De auteurs onderzochten twee verschillende strategieën voor Hamiltoniaan-truncatie om het aantal termen te verminderen:

• Op Coëfficiënt gebaseerde Truncatie ($c_{thre}$): Het weglaten van termen waarvan de absolute waarde onder een drempelwaarde ligt.
• Op Operator-lokaliteit gebaseerde Truncatie ($k_{max}$): Het weglaten van termen waarbij de operator-lokaliteit $k$ (het aantal niet-identiteit Pauli-operatoren in een Pauli-string) een drempelwaarde overschrijdt. Voor termen met twee elektronen werd de lokaliteit gedefinieerd op basis van de excitatie van één $\alpha$- en één $\beta$-elektron om spin-bias te voorkomen.

Nauwkeurigheid werd geëvalueerd door de grondtoestandsenergie-fout ($\Delta E$) en de fideliteit ($F$) te berekenen ten opzichte van niet-getrunceerde referentieberekeningen (Full-CI voor kleine ketens, DMRG voor grotere ketens). Empirische formules werden gefit aan de relatie tussen de truncatie-drempelwaarden en de fideliteit om de drempelwaarden te schatten die nodig zijn om een hoge fideliteit te bereiken ($F \geq 0,99$). Deze geschatte drempelwaarden werden vervolgens gebruikt om de kwantumpoort-aantallen (Clifford, CNOT, en niet-Clifford $R_z$ poorten) te projecteren voor ketens tot $H_{100}$.

Belangrijkste Bijdragen en Resultaten

1. Optimale Truncatiestrategie per Basis:

   • Voor CMO-gebaseerde expansies is coëfficiënt-gebaseerde truncatie ($c_{thre}$) effectief. Zelfs met ruime drempelwaarden blijven de grondtoestandsenergie-fout en de fideliteit dicht bij de ideale waarden, omdat de gedelokaliseerde aard van CMO's betekent dat termen met een hoge lokaliteit vaak noodzakelijk zijn om de golffunctie te beschrijven.
   • Voor LMO-gebaseerde expansies is lokaliteit-gebaseerde truncatie ($k_{max}$) superieur. Coëfficiënt-gebaseerde truncatie faalt onmiddellijk voor LMO's en veroorzaakt een snelle degradatie in nauwkeurigheid. Echter, lokaliteit-gebaseerde truncatie behoudt een hoge fideliteit zelfs bij lage $k_{max}$-waarden, omdat ruimtelijk gescheiden orbitalen in de LMO-basis zwakke correlaties vertonen.

2. Schaling van Truncatie-drempelwaarden:

   • Om $F \geq 0,99$ te handhaven, schaalt de coëfficiënt-drempelwaarde voor CMO's als $c_{thre} \propto n^{-2,07}$.
   • Voor LMO's schaalt de vereiste lokaliteit-drempelwaarde $k_{max}$ polylogaritmisch met de systeemgrootte, specifiek $k_{max} \propto (\log n)^{1,86}$.

3. Schatting van Kwantumpoort-middelen:

   • CMO Basis: Het aantal kwantumpoorten dat vereist is voor Hamiltoniaan-simulatie groeit polynomiaal met de systeemgrootte, specifiek $O(n^4)$ tot $O(n^5)$. Hoewel truncatie de prefactor vermindert, verandert het de polynomiale schalingorde niet.
   • LMO Basis: Wanneer gecombineerd met lokaliteit-gebaseerde truncatie, groeit het aantal poorten polylogaritmisch, specifiek $O((\log n)^{12})$ tot $O((\log n)^{14})$.
   • CNOT Poorten: In de CMO-basis blijven veel termen met een hoge lokaliteit bestaan, zelfs na truncatie, wat leidt tot een hoog CNOT-aantal. In de LMO-basis beperkt de truncatie de Hamiltoniaan tot termen met een lage lokaliteit, waardoor de groei van CNOT-poorten wordt onderdrukt. Voor $n \geq 50$ vereist de LMO-gebaseerde aanpak minder CNOT-poorten dan de CMO-gebaseerde aanpak, ondanks het feit dat de LMO-basis aanvankelijk meer Hamiltoniaan-termen heeft door de gebroken puntgroep-symmetrie.

4. Vergelijking van Encoding:

   • De studie vergeleek ook Jordan–Wigner (JW) en Bravyi–Kitaev (BK) encodings. Hoewel BK de CNOT-schaling voor CMO's vermindert van $O(n^5)$ naar $O(n^4 \log n)$, blijven de resulterende poort-aantallen vergelijkbaar met de LMO-gebaseerde aanpak. Bovendien verhoogt BK het aantal single-qubit Clifford-poorten vanwege het variabele aantal Pauli-X/Y-termen in de strings.

Significantie en Claims
Het artikel stelt dat de keuze van de orbitaalrepresentatie en de bijbehorende truncatiestrategie cruciale factoren zijn bij het bepalen van de haalbaarheid van grootschalige kwantumchemische simulaties. De auteurs concluderen dat het gebruik van LMO's in combinatie met operator-lokaliteit-gebaseerde truncatie een aanzienlijk exponentieel voordeel biedt ten opzien van traditionele CMO-gebaseerde benaderingen voor de Trotter-gebaseerde Hamiltoniaan-simulatie. Specifiek suggereert de overgang van polynomiale naar polylogaritmische schaling in de poort-vereisten dat LMO's de simulatie van veel grotere moleculaire systemen op toekomstige fouttolerante kwantumcomputers mogelijk kunnen maken.

De auteurs merken op dat hoewel hun studie zich richt op 1D waterstofketens waar de operator-lokaliteit direct correleert met de ruimtelijke afstand, de bevindingen nuttige richtlijnen bieden voor basis-ontwerptstrategieën. Ze suggereren dat voor complexe 3D-moleculen afstand-gebaseerde truncatie geschikter zou kunnen zijn dan strikte lokaliteit-gebaseerde truncatie, maar dat het algemene principe van het exploiteren van orbitaal-lokalisatie een veelbelovende weg blijft om computationele kosten te verlagen. De studie benadrukt ook dat staat-preparatie (state preparation) een potentiële bottleneck blijft en dat verdere middelen-inschatting voor staat-preparatie noodzakelijk is voor een volledige end-to-end kostenbeoordeling.