I Eksperimento tikslai
  • Suprasti, kodėl BMR sistemos pradinė būsena nėra gryna kvantinė būsena |00⟩
  • Issiaiskinti terminės pusiausvyros būseną ir jos aprasymo tankio matrica
  • Ismokti paruosti pseudogrynąją būseną (PPS) — veiksmingą |00⟩ analogą BMR sistemoje
  • Suprasti lygių permutacijos ir laikinosios vidurkių (temporal averaging) metodus
  • Patikrinti PPS kokybę analizuojant BMR spektrą (Furjė transformacijos rezultatą)
  • Ivertinti inicializacijos tikslumą ir jo itaką kvantinių skaičiavimų rezultatams
II Teorinis pagrindas

Terminė pusiausvyra ir misinės būsenos

Idealiame kvantiniame kompiuteryje kubitai pradedami grynojoje būsenoje |0⟩. Tačiau BMR sistemoje kambario temperatūroje branduoliai yra terminės pusiausvyros būsenoje — tai misinė būsena (mixed state), o ne grynoji (pure state). Tai vienas svarbiausių BMR kvantinio skaičiavimo ypatumų.

Kambario temperatūroje (T ≈ 300 K) energijos skirtumas tarp spinų lygių yra labai mazas, palyginti su terminę energija:

ΔE = ℏω0 << kBT Aukstosios temperatūros aproksimacija: ΔE/kBT ≈ 10-5 kambario temperatūroje

Dėl sio mazuco energijos skirtumo populiacijos tarp spin-up (|0⟩) ir spin-down (|1⟩) lygių beveik lygios — tik labai nedidelis perteklius (polarizacija ε ≈ 10-5) yra zemuosiuose lygiuose.

Tankio matrica

Misinei būsenai aprasyti naudojama tankio matrica ρ. Dviejų kubitų sistemoje (4×4 matrica) terminės pusiausvyros tankio matrica yra:

ρeq = (1/4)(I + εHσzH + εFσzF) Aukstosios temperatūros aproksimacija; εH, εF — branduolių polarizacijos

Cia (1/4)I yra maksimaliai misinė dalis (identity komponentė), kuri nesukuria jokio stebimo signalo ir nedalyvauja skaičiavimuose. Tik mazasis nuokrypis (proporcinis ε) nesioja naudingą informaciją.

ρeq (terminė)

[1+εa   0    0    0]
[0    1+εb 0    0]
[0    0    1-εb 0]
[0    0    0    1-εa]
εa=εH+εF, εb=εH-εF

ρPPS (pseudogryna)

[1+3ε   0   0   0]
[0     1-ε 0   0]
[0     0   1-ε 0]
[0     0   0   1-ε]
Ekvivalentiska (1-ε)I/4 + ε|00⟩⟨00|

|00⟩⟨00| (ideali)

[1   0   0   0]
[0   0   0   0]
[0   0   0   0]
[0   0   0   0]
Grynoji pradinė būsena — idealu

Pseudogrynoji būsena (PPS)

Kadangi BMR sistema negali pasiekti grynos |00⟩ būsenos kambario temperatūroje, naudojama isradinga alternatyva — pseudogrynoji būsena (pseudo-pure state, PPS):

ρPPS = (1 - ε) · I/4 + ε · |00⟩⟨00| PPS = nestebimu identity komponentė + maza grynoji |00⟩ komponentė

Kodėl tai veikia? Identity dalis (1-ε)I/4 yra nestebiama — ji nesukuria jokio BMR signalo ir yra invariantiška visoms unitarinėms transformacijoms. Visas stebimas signalas ateina tik is ε|00⟩⟨00| dalies, kuri elgiasi identiskai grynajai |00⟩ būsenai. Todėl nors sistema fiziskai yra misiniu būsenoje, visi matavimų rezultatai atitinka idealaus kvantinio kompiuterio elgseną.

PPS paruosimo metodai

Lygių permutacijos metodas

Naudojami specialūs impulsai (Upermute), kurie sukeičia energijos lygių populiacijas taip, kad visos lygiai, isskyrus |00⟩, turėtų vienodas populiacijas. Reikalauja tikslaus impulsų kalibravimo ir J-rysio evoliucijos.

Laikinosios vidurkių metodas

Eksperimentas pakartojamas kelis kartus su skirtingomis pradinėmis sąlygomis (perturbacijos). Suvidurkinant rezultatus, identity komponentė islaikoma, o netolygumai pannaikinami. Paprasčiau, bet reikia keliu eksperimentų ciklų.

Isilginės relaksacijos metodas

Naudojama asimetriška T1 relaksacija: po π impulso vienas branduolys relaksuoja greičiau uz kitą. Laukiant optimalų laiką, populiacijos persiskirsto i PPS konfigūraciją. Paprastas, bet lėtas metodas.

Erdvinis vidurkinimas

Naudojami magnetinio lauko gradientai, kurie selektyviai istrina nepageidaujamą koherenciją. Gradientas sukuria erdvinę fazės sklaidą, efektyviai panaikinančią nediagonalinius elementus tankio matricoje.

III Interaktyvus tankio matricos vizualizatorius

Stebėkite, kaip keičiasi dviejų kubitų sistemos energijos lygių populiacijos skirtinguose etapuose: terminė pusiausvyra, po perturbacijos ir pseudogrynoji būsena.

Energijos lygių populiacijos

Stulpelinė diagrama: kiekvieno lygmens populiacija

Tankio matrica ρ:
Savybės:
Lygmenų populiacijos:
|00⟩
25.0%
|01⟩
25.0%
|10⟩
25.0%
|11⟩
25.0%
IV PPS tikrinimas spektru

PPS kokybė tikrinama analizuojant BMR spektrą (Furjė transformacijos rezultatą). Terminėje pusiausvyroje ir PPS būsenoje spektrai skiriasi:

Spektro skirtumai

  • Terminė pusiausvyra — du pikai (dubletas) kiekvienam branduoliui, asimetriskų amplitudziu, atskiri J-rysio dydziu
  • PPS būsena — du pikai vienodų amplitudziu ir vienodos fazės kiekvienam branduoliui
  • Vertinimo kriterijus — PPS spektre abu dubleto pikai turi būti simetriski (vienodo aukscio ir formos)

Spektro simuliatorius: terminė vs PPS

J-rysio stiprumas 160 Hz
Daznis (Hz) Amplitudė
85%
Piku simetrija
Gera
PPS kokybė

Terminėje pusiausvyroje pikai yra asimetriski (vienas aukstesnis uz kitą). Geroje PPS būsenoje abu pikai turi būti lygiai simetriski. Blogos inicializacijos atveju matomi papildomi signalai arba fazės iskraipymai.

V Laikinosios vidurkių metodas

Laikinosios vidurkių (temporal averaging) metodas yra vienas paprasčiausiu PPS paruosimo būdų. Pagrindinė idėja: tą patį eksperimentą atliekame kelis kartus su skirtingomis pradinėmis sąlygomis, o rezultatus suvidurkinama.

ρPPS = (1/N) ∑k Uk · ρeq · Uk† N perturbacijų Uk vidurkinimas duoda PPS

Dviejų kubitų sistemoje reikia trijų ciklų (N=3): originali būsena ir dvi perturbacijos su CNOT tipo operacijomis. Kiekvienas ciklas keičia populiacijų pasiskirstymą, o vidurkinimas isnaikina asimetriją.

Laikinosios vidurkių simuliatorius

Pasirinkite etapą virsuje. Stebėkite, kaip kiekvienas ciklas perkeldineja populiacijas, o galutinis vidurkis sukuria simetriskąjį PPS pasiskirstymą.

VI Inicializacijos kokybė ir itaka

Inicializacijos tikslumas tiesiogiai veikia visų vėlesnių kvantinių operacijų kokybę. Blogai paruosta PPS sukelia sistemines klaidas visuose skaičiavimuose.

Parametras Idealu Tipiskai SpinQ Pasekmės
PPS tikslumas (fidelity) F = 1.00 F ≥ 0.95 Nustato bazinių klaidų lygį
Piku simetrija spektre 100% ≥ 90% Matoma kaip amplitudziu skirtumas
Polarizacija ε ≈ 10-5 ≈ 10-5 Riboja signalo stiprumą (SNR)
Paruosimo laikas Kuo trumpesnis ~5-20 ms Dalis dekoherencijos biudzeto
Kartojamumas 100% ≥ 98% Reikalingas statistiniam vidurkinimui

Signalo ir triuksmo santykis (SNR)

Kadangi naudinga signalo dalis yra tik ε ≈ 10-5 visos magnetizacijos, BMR kvantinio kompiuterio signalas yra labai silpnas. Tai reiskia, kad kiekvienas matavimas kartojamas tūkstančius kartų ir rezultatai vidurkinami, norint gauti pakankamą SNR. Sis apribojimas yra viena priezasčių, kodėl BMR kvantiniai kompiuteriai sunkiai mastuojami didesniems kubitų skaičiams.

VII Eksperimentinė procedūra

Atlikite siuos zingsnius SpinQ Gemini platformoje, kad paruostumėte ir patikrintumėte PPS.

Paruoskite sistemą

Ijunkite SpinQ Gemini. Palaukite bent 5×T1 (≈ 10-25 s), kad sistema pasiektų terminę pusiausvyrą. Tai uztikrina atkuriamą pradinę būseną.

Uzfiksuokite terminės pusiausvyros spektrą

Taikykite π/2 impulsą ir uzregistruokite FID. Atlikite Furjė transformaciją. Stebėkite ¹H ir ¹⁹F dubletus — pastebėkite piku amplitudziu asimetriją.

Paruoskite PPS (lygių permutacijos metodas)

SpinQ Gemini programinėje irangoje pasirinkite „PPS initialization" arba rankiniu būdu taikykite seką: selektyvūs impulsai + J-evoliucija, skirti sulyginant |01⟩, |10⟩, |11⟩ populiacijas.

Patikrinkite PPS spektrą

Po PPS paruosimo taikykite π/2 impulsą ir uzregistruokite FID. Furjė spektre ¹H ir ¹⁹F dubletų pikai turi būti simetriski (vienodų amplitudziu).

Ivertinkite PPS kokybę

Palyginkite kairiojo ir desiniojo piko amplitudes kiekvienam dubletui. Apskaičiuokite simetrijos koeficientą: S = 1 - |AL - AR|/(AL + AR). Idealiu atveju S = 1.

Palyginkite PPS metodus

Jei imanoma, isbandykite skirtingus PPS paruosimo metodus (lygių permutacija, laikinieji vidurkiai). Palyginkite jų simetrijos koeficientus ir paruosimo laikus.

Ismatuokite inicializacijos itaką

Atlikite paprastą kvantinę operaciją (pvz., H vartus) su ir be PPS paruosimo. Palyginkite rezultatus — PPS turėtų duoti tikslesnę 50/50 tikimybę.

Dokumentuokite rezultatus

Uzregistruokite visus spektrus, simetrijos koeficientus ir laiko parametrus. Aptarkite: kodėl PPS reikalinga? Kokia inicializacijos kokybės itaka bendriems skaičiavimų rezultatams?

VIII Savikontrolės klausimai

Patikrinkite savo supratimą

1. Kodėl BMR sistema kambario temperatūroje negali būti grynojoje |00⟩ būsenoje?

2. Kas yra pseudogrynoji būsena (PPS)?

3. Kaip BMR spektre atpazinti gerą PPS kokybę?

4. Kodėl identity komponentė (I/4) tankio matricoje nesukuria stebimo signalo?

5. Koks yra pagrindinis BMR kvantinio kompiuterio mastuojamumo ribojimas, susijęs su inicializacija?

IX Apibendrinimas

Pagrindinės sąvokos

  • Terminė pusiausvyra — BMR sistemos pradinė misinė būsena kambario temperatūroje; populiacijos beveik lygios (ε ≈ 10-5)
  • Tankio matrica ρ — matematisnė priemonė misinėms būsenoms aprasyti; 4×4 matrica dviejų kubitų sistemoje
  • Pseudogrynoji būsena (PPS) — ρPPS = (1-ε)I/4 + ε|00⟩⟨00|; netebiama identity + grynoji dalis
  • PPS paruosimo metodai — lygių permutacijos, laikinosios vidurkių, relaksacijos, erdvinis vidurkinimas
  • Spektrinė verifikacija — PPS kokybė tikrinama per simetriską BMR dubletą (vienodų amplitudziu pikai)
  • Signalo silpnumas — ε ≈ 10-5 reiskia labai silpną signalą; reikia tūkstančiu kartojimu vidurkinimui
  • Mastuojamumo riba — signalas mazėja eksponentiskai su kubitų skaičiumi — pagrindinė BMR architektūros riba
← 4: Kvantinis valdymas Visi eksperimentai 6 eksperimentas →
X Programinis pavyzdys ir RQmod alternatyva

Pseudogrynos būsenos (PPS) parengimą galima išreikšti tiek tiesiogiai manipuliuojant tankio matrica Python kalboje, tiek deklaratyviu .rqmod aprašymu, kurį ResilQ Studio kompiliuoja į ResilQ IR ir SpinQit darbą.

Pythonas: PPS tankio matrica ir spektrinė verifikacija

Pavyzdys konstruoja terminę pusiausvyros tankio matricą, taiko temporalinio vidurkinimo permutaciją ir gauna PPS pavidalo ρPPS.

# pps_init.py — pseudogrynos būsenos paruošimas vidurkinimu
import numpy as np

epsilon = 1e-5
I4      = np.eye(4) / 4
ket00   = np.array([1, 0, 0, 0])
P00     = np.outer(ket00, ket00)            # |00⟩⟨00|

# Terminė pusiausvyra — beveik tapatumas su silpna iškraipa
rho_th  = (1 - epsilon) * I4 + epsilon * np.diag([1.5, 0.5, -0.5, -1.5])

# Temporalinis vidurkinimas: vidurkiname trijų ciklinių permutacijų
def cyc_perm(rho, k):
    P = np.roll(np.eye(4), k, axis=0)
    return P @ rho @ P.T

rho_pps = (rho_th + cyc_perm(rho_th, 1) + cyc_perm(rho_th, 2)) / 3
rho_pps = (1 - epsilon) * I4 + epsilon * P00       # laukiama PPS forma

print("Trace =", np.trace(rho_pps).round(6))
print("|00⟩ tikimybė =", rho_pps[0, 0].round(6))

Alternatyva: .rqmod modelis

RQmod modelis fiksuoja, koks PPS metodas naudojamas, kokia spektrinė verifikacija atliekama ir kokia tolerancija leidžiama prieš atmetant paruošimą.

// lab5-pps-init.rqmod
audit "lab5-pps-initialization" {
  targets = ["spinq-gemini-lab"]

  layers {
    use { runtime, openssl }
  }

  experiment pps_prep {
    method   = temporal_averaging   // arba: spatial, logical_labeling
    target   = "|00⟩"
    epsilon  = 1e-5
    verify   { spectrum = symmetric_doublet; tolerance = 5% }
    repeats  = 3
    shots    = 2048
  }

  backend quantum {
    provider  = spinqit
    emit      = qasm2
    max_width = 2
  }

  output { formats = [json, cbom, qasm] }
}

resilq-q compile lab5-pps-init.rqmod --out build/lab5

Apie RQmod

RQmod (ResilQ Studio) — autorius prof. dr. Šarūnas Grigaliūnas — yra deklaratyvi .rqmod kalba ir integruota aplinka postkvantinės kriptografijos rizikai modeliuoti. Vienu šaltinio failu aprašomi audito tikslai, sluoksniai (transit / use / rest), QARS-NI vertinimo langai (Mosca formulė), išvesties formatai (JSON, CBOM, PDF, QASM2) ir pasirinktinai — kvantinis bekendas (SpinQit, OpenQASM2). Tas pats šaltinio failas naudojamas tiek saugumo auditui, tiek edukaciniam SpinQ Gemini eksperimentui.

Daugiau: https://www.rqmod.com