Rovibrational energy levels of H$_2$O by quantum computing

In dit artikel wordt beschreven hoe met een quantumcomputer op basis van gevangen ionen de ro-vibratie-energieniveaus van watermoleculen worden berekend, waarbij een hybride quantum-klassieke aanpak wordt gebruikt om een nauwkeurigheid van enkele cm⁻¹ te bereiken.

De kern

🌊 De Dans van het Watermolecuul: Hoe een Quantumcomputer de Beweging van Water Begrijpt

Stel je een watermolecuul ($H_2O$) voor als een klein, levendig balletje met twee armen (de waterstofatomen) en een hoofd (het zuurstofatoom). Dit balletje zit nooit stil. Het trilt (de armen bewegen in en uit), het buigt (de hoek verandert) en het draait als een tol.

In de natuurkunde noemen we deze combinatie van trillen en draaien ro-vibratie. Om te begrijpen hoe water zich gedraagt (bijvoorbeeld in wolken of in je lichaam), moeten we precies weten hoeveel energie deze dans kost.

Vroeger was dit een onmogelijke taak voor onze gewone computers. Het is alsof je probeert de exacte beweging van duizenden dansers tegelijk te voorspellen terwijl ze allemaal een beetje uit hun ritme raken. De berekeningen worden zo complex dat zelfs de krachtigste supercomputers er van duizelig worden.

Maar in dit nieuwe onderzoek hebben wetenschappers van RIKEN (Japan) en de ELTE (Hongarije) een nieuwe aanpak geprobeerd: ze gebruikten een Quantumcomputer.

1. Het Grote Muziekstuk (De Hamiltoniaan)

Om de energie van het watermolecuul te berekenen, schrijven wetenschappers een enorme vergelijking op, een soort muziekpartituur die alle mogelijke bewegingen beschrijft. Deze heet de Watson-Hamiltoniaan.

• Het probleem: Deze partituur is gigantisch. Het bevat duizenden noten (termen) die allemaal met elkaar verweven zijn.
• De oplossing: De onderzoekers hebben deze partituur vertaald naar de taal van de quantumcomputer: qubits. In plaats van gewone bits (0 of 1), gebruiken qubits de vreemde eigenschappen van de quantumwereld om deze complexe partituur te "bevatten".

2. De Muzikale Zolder (De QSCI-methode)

Normaal gesproken zou een quantumcomputer proberen om alle noten in de partituur tegelijk te spelen om het antwoord te vinden. Maar omdat de huidige quantumcomputers nog een beetje "ruis" hebben (zoals een slecht gestemde piano), is dat te lastig.

De onderzoekers gebruikten een slimme truc, genaamd QSCI (Quantum-Selected Configuration Interaction).

• De analogie: Stel je voor dat je een enorme zolder vol met oude dozen hebt (alle mogelijke bewegingen). Je wilt weten welke dozen de belangrijkste zijn voor het liedje.
• De methode: In plaats van elke doos te openen, laat je de quantumcomputer even snel door de zolder "glijden" (tijdsevolutie). De computer zegt dan: "Hey, deze dozen hier hebben de meeste geluid!"
• De quantumcomputer selecteert dus alleen de belangrijkste dozen (de meest waarschijnlijke bewegingen).
• Vervolgens nemen de onderzoekers die selectie mee naar een gewone computer. Die gewone computer is heel goed in het ordenen van een kleine, geselecteerde lijst. Die berekent dan de exacte energie van die specifieke bewegingen.

Het is alsof de quantumcomputer de kandidaten selecteert voor een auditie, en de gewone computer de jury is die de winnaars kiest.

3. Het Experiment: Water op de Quantumcomputer

De onderzoekers hebben dit getest op de Quantinuum Reimei, een quantumcomputer die gebruikmaakt van gevangen ionen (atomen die zweven in een magnetisch veld).

• Ze lieten de computer de energie berekenen voor de laagste danspassen van het watermolecuul.
• Het resultaat? Ze haalden een nauwkeurigheid van ongeveer 1 cm⁻¹ (een heel klein stukje energie).
• Voor de wetenschap is dit "spectroscopische nauwkeurigheid". Het betekent dat hun berekening bijna perfect overeenkomt met wat we in het echt met lasers meten in een laboratorium.

4. Waarom is dit belangrijk?

Vroeger konden we alleen de elektronen (de kleine deeltjes rondom de atomen) goed berekenen. Maar de beweging van de atomen zelf (trillen en draaien) was te lastig.

• De toekomst: Als we dit kunnen doen voor water, kunnen we het ook doen voor grotere, complexere moleculen. Denk aan medicijnen die in je lichaam bewegen, of materialen die we nog niet hebben ontdekt.
• De "Killer-app": Dit onderzoek toont aan dat zelfs de huidige, nog wat onvolmaakte quantumcomputers (die ruis hebben) nuttig kunnen zijn voor chemie. Ze kunnen namelijk taken doen die voor gewone computers te zwaar zijn, maar die we wel kunnen controleren met echte metingen.

Samenvatting in één zin

De onderzoekers hebben een slimme samenwerking tussen een quantumcomputer (die de beste bewegingen selecteert) en een gewone computer (die de exacte energie berekent) gebruikt om de dans van een watermolecuul met extreme precisie te voorspellen, een stap dichter bij het oplossen van de grootste mysteries in de chemie.

Technische samenvatting

Hier is een gedetailleerde technische samenvatting van het artikel "Rovibrational energy levels of H2O by quantum computing" van Lötstedt en Szidarovszky, geschreven in het Nederlands.

Titel: Rovibratie-energieniveaus van H2O berekend met kwantumcomputing

Auteurs: Erik Lötstedt en Tamás Szidarovszky
Instituut: RIKEN (Japan) en ELTE (Hongarije)
Datum: 9 maart 2026

---

1. Het Probleem

Traditioneel ligt de focus van kwantumcomputing in de chemie op het oplossen van het elektronische structuurprobleem (de golffuncties en energieën van elektronen). Echter, om volledige moleculaire energieniveaus te bepalen, moet ook de Schrödinger-vergelijking voor de gekoppelde rotatie- en vibratiebeweging (rovibratie) van het molecuul worden opgelost.

• Complexiteit: Het aantal variabelen in de rovibratie-golffunctie groeit lineair met het aantal atomen ($N_{var} = 3N_a - 3$), wat leidt tot een exponentiële groei van de Hilbertruimte bij klassieke methoden.
• Beperkingen van bestaande methoden: Hoewel er veel literatuur is over elektronische structuur, zijn er weinig rapporten over rovibratie. Eerdere studies met kwantumcomputers hebben vaak alleen vibratie (zonder rotatie) behandeld of vereisten fouttolerante hardware.
• Noodzaak van koppeling: Voor spectroscopische nauwkeurigheid (binnen 1 cm⁻¹) is het essentieel om de koppeling tussen rotatie en vibratie (rovibratie-koppeling) en Coriolis-koppeling correct te modelleren. Bestaande benaderingen op huidige "noisy" (ruisgevoelige) kwantumcomputers zijn vaak te beperkt door de enorme grootte van de Hamilton-operator.

2. Methodologie

De auteurs presenteren een hybride klassiek-kwantum aanpak om de rovibratie-energieniveaus van het watermolecuul (H₂O) te berekenen.

A. Theoretische Basis: De Watson-Hamiltoniaan

• Ze leiden de qubit-vorm af van de Watson-Hamiltoniaan, die de volledige rovibratie-dynamiek beschrijft.
• De Hamiltoniaan ($H$) wordt opgesplitst in termen:
  • $H_{RR}$: Rigid-rotor (stijf rotor) term.
  • $H_{HO}$: Harmonische oscillator term.
  • $V_{anharm}$: Anharmonische potentiaal (Taylor-reeks tot 4e orde).
  • $H_{vibCor}$: Vibratoire Coriolis-koppeling.
  • $H_{rovib}$: Koppeling tussen rotatie en vibratie.
• Qubit-encoding: Ze gebruiken een binair mapping-schema om de bosonische vrijheidsgraden (vibratie- en rotatiekwantumgetallen) te coderen op qubits. Dit is compact en vereist minder qubits dan andere methoden.
• Complexiteit: De resulterende qubit-Hamiltoniaan bestaat uit duizenden Pauli-strings (bijv. >200.000 termen voor hogere rotatietoestanden), wat standaard methoden zoals VQE (Variational Quantum Eigensolver) onpraktisch maakt vanwege de enorme meetkosten.

B. Algoritme: Quantum-Selected Configuration Interaction (QSCI)

Om de beperkingen van ruis en het groot aantal termen te omzeilen, gebruiken ze een variant van de QSCI-methode (ook wel Sample-based Quantum Diagonalization genoemd):

1. Kwantumstap: Een proef-golffunctie $|\psi(\tau, N_{ST})\rangle$ wordt voorbereid op de kwantumcomputer via tijdevolutie (Suzuki-Trotter-decompositie) van een initiële toestand.
2. Sampling: De golffunctie wordt gemeten in de computationele basis om een kansverdeling $p(x)$ te verkrijgen.
3. Selectie: Een subset van basisstaten ($\Omega$) wordt geselecteerd op basis van deze verdeling (staten met hoge waarschijnlijkheid).
4. Klassieke stap: De Hamiltoniaan wordt beperkt tot deze geselecteerde basis $\Omega$. De matrix wordt op een klassieke computer gegenereerd en diagonaal gemaakt om de energieniveaus te vinden.

C. Experimentele Setup

• Hardware: Berekeningen zijn uitgevoerd op de Quantinuum Reimei (een 20-qubit ionenval-kwantumcomputer) op het RIKEN-campus.
• Foutcorrectie: Er wordt gebruik gemaakt van post-selectie om de pariteit (symmetrie onder protonuitwisseling) te behouden, wat een eenvoudige vorm van foutmitigatie is.
• Doel: Berekening van de laagste rovibratie-niveaus van H₂O voor verschillende rotatiekwantumgetallen ($J$).

3. Belangrijkste Bijdragen

1. Volledige Hamiltoniaan: Voor het eerst wordt de volledige Watson-Hamiltoniaan, inclusief alle rovibratie-koppelings- en Coriolis-termen, succesvol omgezet naar qubit-vorm en toegepast op een echte kwantumcomputer.
2. Hybride QSCI-toepassing: Het aantonen dat QSCI zeer effectief is voor rovibratieproblemen, zelfs op ruisgevoelige hardware, door de zware diagonalisatie op een klassieke computer te verplaatsen.
3. Nauwkeurigheid: Het bereiken van spectroscopische nauwkeurigheid (enkele cm⁻¹) voor lage energieniveaus, wat een mijlpaal is voor de toepassing van NISQ (Noisy Intermediate-Scale Quantum) apparatuur op moleculaire dynamica.
4. Validatie: Het vergelijken van de resultaten met exacte klassieke diagonalisatie en perturbatietheorie, wat de betrouwbaarheid van de kwantumbenadering bevestigt.

4. Resultaten

• Nauwkeurigheid: De methode levert energieniveaus op met een foutmarge van enkele cm⁻¹ (vaak < 1 cm⁻¹) voor de laagste vibratie- en rovibratie-toestanden van H₂O.
• Invloed van koppelingen: De resultaten tonen aan dat het negeren van de rovibratie-koppeling ($H_{rovib}$) en Coriolis-koppeling leidt tot significante fouten (tot wel 30 cm⁻¹ verschuiving in energie). Zelfs voor de grondtoestand ($J=0$) is de Coriolis-koppeling cruciaal.
• Convergentie:
  • Voor $J=0$ en $v_{max}=3$ (6 qubits) convergeren de kwantumberekeningen snel naar de exacte waarden naarmate het aantal Suzuki-Trotter-stappen toeneemt.
  • Zelfs met ruis (Reimei-data) worden de resultaten binnen 1-3 cm⁻¹ van de exacte waarden, wat vaak beter is dan wat verwacht wordt van ruisgevoelige hardware.
  • Interessant genoeg bleek dat voor sommige toestanden de ruis op de kwantumcomputer soms hielp bij het vinden van geschikte basisfuncties (beter dan de ruisvrije simulatie in specifieke gevallen).
• Schalbaarheid: De studie toont aan dat de methode werkt voor $J=1$ en hogere vibratietoestanden, hoewel de grootte van de benodigde basisset toeneemt.

5. Betekenis en Toekomstperspectief

• Verifieerbaar Kwantumvoordeel: Rovibratie-berekeningen bieden een uniek platform om kwantumvoordeel te demonstreren. In tegenstelling tot elektronische structuurproblemen, waar experimentele data vaak minder precies is, kunnen rovibratie-overgangen met spectroscopie tot op $10^{-6}$ cm⁻¹ worden gemeten. Dit maakt het mogelijk om kwantumberekeningen direct en strikt te valideren.
• Toepassing op grotere systemen: Hoewel dit artikel H₂O behandelt, is de methode schaalbaar naar grotere moleculen (zoals benzeen) of moleculaire clusters. Voor deze systemen is de klassieke diagonalisatie onmogelijk, maar de kwantumcomputer kan de cruciale basisstaten selecteren.
• Hardware-onafhankelijkheid: De succesvolle toepassing op een ionenval-computer (Reimei) met all-to-all connectiviteit suggereert dat deze hybride aanpak veelbelovend is voor de huidige generatie kwantumhardware, zonder dat volledige foutcorrectie nodig is.

Conclusie: Dit werk markeert een belangrijke stap in de kwantumchemie door te bewijzen dat huidige ruisgevoelige kwantumcomputers, in combinatie met geavanceerde hybride algoritmen (QSCI), in staat zijn om complexe, gekoppelde moleculaire bewegingen met spectroscopische nauwkeurigheid te simuleren.