Classical simulation and quantum resource theory of non-Gaussian optics

Titel: Klassieke simulatie en kwantumbrontheorie van niet-Gaussische optica

Auteurs: Oliver Hahn, Ryuji Takagi, Giulia Ferrini, en Hayata Yamasaki.

---

1. Het Probleem

Kwantumcomputing met continue variabelen (CV-systemen), zoals optische systemen, biedt veelbelovende voordelen zoals inherentere weerstand tegen ruis. Echter, om universele kwantumcomputing te bereiken, zijn niet-Gaussische bronnen (zoals Fock-toestanden of GKP-toestanden) noodzakelijk, aangezien puur Gaussische operaties en toestanden efficiënt klassiek te simuleren zijn en geen exponentiële kwantumvoordeel bieden.

De uitdaging ligt in het evalueren en ontwerpen van protocollen voor deze systemen:

• Klassieke simulatie: Het is algemeen moeilijk om kwantumsystemen klassiek te simuleren, vooral CV-systemen. Bestaande methoden hebben beperkingen:
  • Wigner-functie: Simulatiekosten schalen met de hoeveelheid negativiteit in de Wigner-functie. Dit werkt goed voor lage negativiteit, maar faalt bij complexe niet-Gaussische toestanden.
  • Fock-basis decompositie: Methodes die toestanden ontleden in een Fock-basis (zoals de "stellar rank") lopen vast bij toestanden die oneindig veel Fock-toestanden vereisen (zoals GKP- en kat-toestanden), wat leidt tot schaalingsproblemen.
• Brontheorie: Er ontbreekt een operationele maatstaf die direct de kosten van klassieke simulatie koppelt aan de hoeveelheid "niet-Gaussischheid" in een toestand.

2. Methodologie

De auteurs introduceren twee nieuwe klassieke simulatiealgoritmen die gebaseerd zijn op het ontleden van niet-Gaussische toestanden in een lineaire superpositie van Gaussische toestanden.

A. Fase-gevoelige Gaussische simulator
Om de relatieve fasen in superposities van Gaussische toestanden correct te volgen (wat de standaard covariantiematrix-formalisme niet doet), ontwikkelen de auteurs een uitgebreide methode:

• Ze gebruiken een referentiestaat $|G_0\rangle$ en berekenen de overlap tussen twee Gaussische toestanden $|G_1\rangle$ en $|G_2\rangle$ via de trace $\text{Tr}(G_0 G_1 G_2)$.
• Hierdoor kunnen ze de globale fase $\theta$ bepalen in $\langle G_1 | G_2 \rangle = e^{i\theta} |\langle G_1 | G_2 \rangle|$, essentieel voor het correct simuleren van interferentie in niet-Gaussische superposities.

B. Exacte Simulatie Algoritme

• Principe: Een niet-Gaussische toestand $|\psi\rangle$ wordt geschreven als $|\psi\rangle = \sum_{i=1}^\chi c_i |G_i\rangle$.
• Kosten: De simulatie vereist het berekenen van alle kruis-overlappen in de Born-waarschijnlijkheid. De complexiteit schaalt kwadratisch met het aantal termen: $O(\chi^2)$, waarbij $\chi$ het aantal Gaussische toestanden in de decompositie is.

C. Benaderende Simulatie Algoritme (Sparsificatie)

• Principe: Om de kosten te verlagen, wordt de toestand gesparseerd (verdund) door termen in de decompositie te stochastisch te selecteren op basis van hun gewicht $|c_i|$.
• Snelheidsnorm-schatting: Omdat de gesparseerde toestand niet genormaliseerd is, ontwikkelen de auteurs een "fast norm estimation" algoritme. Dit gebruikt steekproeven uit een Gaussische ensemble van verplaatsingsoperatoren om de norm van de toestand lineair te schatten in plaats van kwadratisch.
• Kosten: De complexiteit schaalt lineair met de $l_1$-norm van de coëfficiënten ($\|c\|_1$), wat aanzienlijk efficiënter is voor toestanden met een groot aantal termen.

3. Belangrijkste Bijdragen

A. Definities van Nieuwe Maatstaven
De auteurs definiëren twee maatstaven binnen de kwantumbrontheorie van niet-Gaussischheid die direct corresponderen met de simulatiekosten:

1. Gaussische Rang (Gaussian Rank, $\chi$): Het minimale aantal Gaussische toestanden nodig om een toestand exact te decomponeren. Dit correleert met de kosten van het exacte algoritme.
2. Gaussische Extent (Gaussian Extent, $\xi$): Het kwadraat van de minimale $l_1$-norm van de coëfficiënten in een decompositie ($\xi = \inf (\sum |c_i|)^2$). Dit correleert met de kosten van het benaderende algoritme.

B. Theoretische Eigenschappen

• Monotonie: Beide maatstaven zijn monotoon onder Gaussische operaties (unitaire transformaties, verplaatsingen, homodyne/heterodyne detectie). Ze nemen niet toe door gratis operaties.
• Verbinding met Robuustheid: De auteurs tonen aan dat de Gaussische extent gerelateerd is aan de lower semicontinuous robustness van niet-Gaussischheid. Dit stelt hen in staat om optimale decomposities te vinden voor specifieke toestanden door gebruik te maken van getuige-operatoren (witness operators).
• Niet-multiplicativiteit: In tegenstelling tot wat vaak wordt verondersteld, is de Gaussische extent niet multiplicatief voor tensorproducten. Dit betekent dat het simuleren van multi-modale systemen efficiënter kan zijn dan het simuleren van individuele modi apart (sub-multiplicatieve schaling).

4. Resultaten en Toepassingen

A. Optimale Decomposities
De auteurs berekenen optimale decomposities voor toestanden die cruciaal zijn voor CV-kwantumcomputing:

• Fock-toestanden ($|1\rangle$): De optimale decompositie is een superpositie van gecomprimeerde coherent-toestanden met specifieke parameters ($\alpha = 2/3, \xi = \ln\sqrt{3}$). De auteurs tonen aan dat de Gaussische extent voor $|1\rangle$ gelijk is aan $4e / (3\sqrt{3})$.
• GKP-toestanden: Deze kunnen worden geschreven als een som van gecomprimeerde toestanden. De auteurs tonen aan dat de Gaussische extent hier veel lager is dan de stellar rank, wat de simulatie haalbaarder maakt.
• Kat-toestanden (Cat states): Voor grote amplitudes convergeren de decomposities naar een optimale vorm met $\|c\|_1 \approx \sqrt{2}$.

B. Toepassingen in Protocollen

• Gaussische Boson Sampling (GBS): De auteurs geven een bovengrens voor de simulatiekosten van GBS. Hoewel de dynamica volledig binnen de "gratis" operaties van de theorie van niet-klassicaliteit valt, levert het gebruik van de bredere Gaussische decompositie (in plaats van alleen coherent-toestanden) een betere simulatiecapaciteit op.
• Kat-toestanden Kweken (Cat State Breeding): Ze gebruiken de Gaussische extent om een fundamentele ondergrens af te leiden voor het aantal kat-toestanden dat nodig is om een GKP-toestand (sensor toestand) te genereren. Deze grenzen zijn onafhankelijk van de grootte van de input-kat-toestanden en zijn strikter dan die gebaseerd op Wigner-negativiteit of stellar rank.

5. Betekenis en Conclusie

Dit werk vormt een brug tussen praktische klassieke simulatie en fundamentele kwantumbrontheorie:

1. Efficiëntie: De voorgestelde algoritmes bieden een nieuwe weg om complexe niet-Gaussische CV-systemen te simuleren, waarbij de kosten direct gekoppeld zijn aan de hoeveelheid niet-Gaussische hulpbronnen.
2. Operationele Interpretatie: De "Gaussische extent" krijgt een directe operationele betekenis als maatstaf voor de simulatiekosten, wat een nieuw perspectief biedt op het kwantificeren van kwantumvoordeel.
3. Fundamentele Inzicht: De ontdekking dat de Gaussische extent niet multiplicatief is, suggereert dat multi-modale correlaties essentieel zijn voor het begrijpen van de schaalbaarheid van kwantumvoordeel in CV-systemen.
4. Toekomstperspectief: De methoden zijn niet beperkt tot input-toestanden; ze kunnen ook worden toegepast op circuits met niet-Gaussische poorten via "gadgetisatie", vergelijkbaar met methoden voor Clifford+T circuits in qubit-systemen.

Samenvattend bieden de auteurs een krachtig raamwerk om de grens tussen klassiek en kwantumrekenvermogen in optische systemen te verkennen, met directe toepassingen voor het ontwerpen van foutcorrectiecodes en het evalueren van de haalbaarheid van toekomstige kwantumexperimenten.