QwaveMPS: An efficient open-source Python package for simulating non-Markovian waveguide-QED using matrix product states

Probleemstelling

Waveguide-kwantumelektrodynamica (waveguide-QED) bestudeert de interactie tussen kwantumemitters (zoals atomen of supergeleidende transmons) en een continuüm van gekwantiseerde veldmodi in quasi-één-dimensionale systemen. Deze systemen vertonen rijke fenomenen zoals tijdsvertraagde feedback (niet-Markoviaanse effecten), chirale koppeling en sterke niet-lineariteiten.

Het modelleren van deze systemen vormt echter een aanzienlijke uitdaging:

1. Explosie van de Hilbertruimte: De grootte van de Hilbertruimte groeit exponentieel met het aantal fotonen en emitters, wat volledige simulaties onmogelijk maakt voor grotere systemen.
2. Beperkingen van bestaande methoden: Traditionele benaderingen, zoals de master-vergelijking en input-output theorieën, maken vaak gebruik van de Markov-benadering. Hierbij worden tijdsvertraagde feedback-effecten genegeerd en worden koppelingsraten als instantane beschouwd. Dit faalt in regimes met sterke niet-lineariteiten of lange vertragingstijden.
3. Beperkingen van geavanceerdere methoden: Methoden die verder gaan dan de Markov-benadering (zoals verstoringstheorie of stochastische methoden) zijn vaak beperkt tot zwakke excitaties of lage fotonnummers, waardoor ze ongeschikt zijn voor sterk niet-lineaire dynamica.
4. Gebrek aan toegankelijke tools: Hoewel Matrix Product State (MPS) methoden theoretisch geschikt zijn, ontbreekt er vaak specifieke documentatie en gebruiksvriendelijke software binnen de kwantumoptica-gemeenschap voor deze specifieke toepassingen.

Methodologie

De auteurs introduceren QwaveMPS, een open-source Python-pakket dat gebruikmaakt van Matrix Product States (MPS) gebaseerd op tensornetwerktheorieën om waveguide-QED-systemen exact en schaalbaar te simuleren.

Kernprincipes van de methode:

• Discretisatie in tijd: Het elektromagnetische veld wordt gediskretiseerd in "tijdbins" (time bins). Dit beperkt de groei van de Hilbertruimte door het veld te behandelen als een keten van discrete tijdstappen in plaats van een continuüm.
• Tensornetwerk representatie:
  • De kwantumtoestand wordt weergegeven als een MPS, waarbij de emitters en de gediskretiseerde veldtijdstappen worden gemodelleerd als tensoren.
  • Operatoren (zoals de Hamiltoniaan) worden omgezet in Matrix Product Operators (MPO).
• Tijdevolutie: De systeemevolutie wordt berekend door de tijdevolutie-operator (U) toe te passen op relevante bins.
  • In het Markoviaanse regime werkt U op het systeem en de huidige tijdbin (2-site MPO).
  • In het niet-Markoviaanse regime (bijv. met spiegels of feedback) wordt een extra "feedback-bin" (tijdsvertraging) meegenomen. De operator wordt dan een 3-site MPO, waarbij een swap-operator wordt gebruikt om de feedback-bin naast het systeem te plaatsen voor de evolutie.
• SVD en Orthogonaliteit: Na elke tijdstap wordt de toestand opnieuw gedecomposeerd via Singular Value Decomposition (SVD) om de MPS-vorm te behouden en de bond-dimensie (complexiteit) te beheersen.

Software-architectuur (QwaveMPS):
Het pakket bestaat uit zes scripts die een modulaire workflow bieden:

• states.py en hamiltonians.py: Definieren initiële toestanden (bijv. TLS, vacuüm, Fock-pulsen) en Hamiltonianen.
• parameters.py: Bevat input-parameters (tijdstap, maximale bond-dimensie, koppelingssterktes).
• simulation.py: Voert de tijdevolutie uit (zowel voor Markoviaanse als niet-Markoviaanse gevallen).
• operators.py en correlation.py: Berekenen observabelen zoals populaties, entanglement-entropie en twee-tijds correlatiefuncties.

Belangrijkste Bijdragen

1. Open-source implementatie: QwaveMPS is de eerste gebruiksvriendelijke, open-source Python-library die MPS-methoden toegankelijk maakt voor waveguide-QED, inclusief complexe niet-Markoviaanse scenario's.
2. Behandeling van niet-lineariteit en feedback: Het pakket kan systemen simuleren met sterke niet-lineariteiten, tijdsvertraagde feedback en meerdere propagende quanta, zonder de beperkingen van de Markov-benadering.
3. Efficiëntie: Door zich te richten op de relevante delen van de Hilbertruimte via MPS, worden berekeningen mogelijk gemaakt op standaard werkstations (zelfs laptops), terwijl volledige Hilbertruimte-methoden onuitvoerbaar zouden zijn.
4. Veelzijdigheid: Het ondersteunt diverse scenario's:
   • Lineaire en niet-lineaire regimes.
   • Klassieke aandrijvingen (CW en gepulseerd) en kwantumpulsen (Fock-toestanden).
   • Enkele en meerdere emitters (TLS's).
   • Chirale en symmetrische koppelingen.

Resultaten

De auteurs presenteren een reeks simulaties die de kracht van het pakket demonstreren:

• Lineair regime (Vacuüm dynamica):

  • Simulatie van het verval van een enkele TLS in een oneindige en semi-oneindige waveguide. De resultaten komen perfect overeen met analytische oplossingen en ODE-oplossers.
  • Demonstratie van destructieve en constructieve feedback (via een spiegel), waarbij de populatiedynamica en fotonflux in de feedbacklus nauwkeurig worden berekend.
  • Simulatie van twee TLS's, waarbij de overgang van Markoviaanse naar niet-Markoviaanse dynamica (met vertraging) wordt getoond.

• Niet-lineair regime (Meerdere quanta):

  • Simulatie van twee tegelijkertijd geëxciteerde TLS's. Hier wordt entanglement-entropie berekend tussen de emitters en het veld, wat in het niet-Markoviaanse regime afwijkt van het Markoviaanse geval.
  • Berekening van de dynamica duurt minder dan een seconde, zelfs met hoge bond-dimensies.

• Klassieke aandrijvingen:

  • CW-pomp: Berekening van de Mollow-triplet in het spectrum en steady-state correlatiefuncties ($g^{(1)}$ en $g^{(2)}$). Het pakket kan sterke aandrijvingen hanteren waar methoden met lage fotonnummers falen.
  • Gepulseerde aandrijving: Simulatie van tijdsafhankelijke klassieke pulsen (bijv. Gaussische $\pi$-pulsen) en de daaruit voortvloeiende twee-tijds correlaties.

• Kwantum-pulsen (Fock-toestanden):

  • Directe injectie van 1- en 2-foton Fock-pulsen als initiële toestand.
  • Berekening van $G^{(1)}$ en $G^{(2)}$ correlaties voor deze kwantumtoestanden, waarbij kenmerkende patronen (zoals de "bird-shape" in $G^{(2)}$ voor 2-fotonen) worden waargenomen.

Berekeningstijden:
Alle voorbeelden worden uitgevoerd op een standaard workstation (18 cores, 527 GB RAM, maar slechts ~180 MB gebruikt). De meeste simulaties duren minder dan een seconde; zelfs complexe niet-Markoviaanse berekeningen met hoge bond-dimensies blijven binnen enkele seconden tot minuten.

Betekenis en Conclusie

QwaveMPS vormt een brug tussen geavanceerde theoretische technieken (tensornetwerken) en praktische, reproduceerbare studies in kwantumoptica. Het pakket lost het probleem van de exponentiële groei van de Hilbertruimte op voor waveguide-QED-systemen en stelt onderzoekers in staat om:

• Niet-Markoviaanse effecten exact te modelleren zonder benaderingen.
• Sterk niet-lineaire dynamica te bestuderen met meerdere fotonen.
• Kwantumcorrelaties en entanglement in complexe circuits te analyseren.

De auteurs benadrukken dat het pakket nog uitbreidbaar is (bijv. naar meer dan twee emitters, hogere energieniveaus, en het toevoegen van verliesprocessen zoals pure dephasing). QwaveMPS positioneert zich hiermee als een essentieel open-source hulpmiddel voor de kwantumoptica-gemeenschap, vooral in het licht van de opkomst van supergeleidende kwantumcircuits en nanofotonica.