I Eksperimento tikslai
  • Suprasti Rabi osciliacijos fizikini reikinį — periodiskus perėjimus tarp dviejų kvantinių būsenų
  • Ismatuoti Rabi daznį ir suprasti jo priklausomybę nuo RD impulso amplitudės
  • Issiaiskinti, kaip Rabi osciliacijos naudojamos kvantiniu vartu kalibravimui
  • Nustatyti tikslią impulso trukmę, reikalingą π/2 (90 laipsnių) ir π (180 laipsnių) pasūkimams
  • Stebėti Rabi osciliacijas SpinQ Gemini platformoje abiem kubitais (¹H ir ¹⁹F)
  • Suprasti dekoherencijos itaką — kodėl Rabi osciliacijos gęsta
II Teorinis pagrindas

Dviejų lygmenų kvantinė sistema

Branduolinis sukinys I = 1/2 magnetiniame lauke sudaro dviejų lygmenų kvantinę sistemą. Zemesnė energijos būsena |0⟩ (spin-up, lygiagretis laukui) ir aukstesnė būsena |1⟩ (spin-down, priesingai laukui) yra atskirtos energijos tarpu ΔE = ℏω₀, kur ω₀ yra Larmor daznis.

Bendroji kubito būsena bet kuriuo momentu yra superpozicija:

|ψ(t)⟩ = c₀(t)|0⟩ + c₁(t)|1⟩ kur |c₀|² + |c₁|² = 1 — tikimybiu normavimo sąlyga

Rabi osciliacijos

Kai dviejų lygmenų sistemą veikiame rezonansiniu RD impulsu (daznis ω = ω₀), populiacijos tikimybės periodiskai osciliuoja tarp |0⟩ ir |1⟩ būsenų. Sis reiskinys vadinamas Rabi osciliacijomis.

Populiacijos tikimybė rasti sistemą būsenoje |1⟩ po laiko t:

P|1⟩(t) = sin²(ΩR · t / 2) kur ΩR = γ · B₁ — Rabi daznis, B₁ — RD impulso magnetinio lauko amplitudė

Rabi daznis ΩR yra tiesiogiai proporcingas impulso stiprumui B₁. Didesnis impulso stiprumas — greitesnės osciliacijos — trumpesnis laikas, reikalingas kubito būsenos apvertimui.

|0⟩ |1⟩ x
t = 0
Pradinė būsena |0⟩
|0⟩ |1⟩ x
π/2 impulsas
Superpozicija (|0⟩+|1⟩)/√2
|0⟩ |1⟩ x
π impulsas
Apverstoji būsena |1⟩
|0⟩ |1⟩ x
3π/2 impulsas
Superpozicija (|0⟩-|1⟩)/√2

Rysys su kvantiniais vartais

Rabi osciliacijos yra pagrindinis mechanizmas kvantiniams vartams realizuoti BMR kvantiniame kompiuteryje. Tiksliai kontroliuodami RD impulso trukmę, galime isgyvendinti bet kuri vieno kubito vartų operaciją:

Impulso kampas Impulso trukmė Kvantiniai vartai Poveikis
π/2 (90 laipsnių) t = π / (2ΩR) Hadamardo tipo |0⟩ → (|0⟩+|1⟩)/√2 superpozicija
π (180 laipsnių) t = π / ΩR Pauli-X (NOT) |0⟩ → |1⟩ bitų apvertimas
2π (360 laipsnių) t = 2π / ΩR Identitetas |0⟩ → |0⟩ pilnas ciklas

Gesimas ir dekoherencija

Realybėje Rabi osciliacijos nėra ideolios — jos pamazu gęsta dėl dekoherencijos. Pagrindinės gesimo priezastys:

T2 relaksacija

Skersinė relaksacija sukelia fazės nuostolius. Kuo trumpesnis T2, tuo greičiau gęsta osciliacijos. Rabi osciliacijos su gesimu: P(t) = sin²(ΩRt/2) · e-t/T2

B₁ lauko nehomogeniskumas

Jei RD impulso laukas nėra vienodas visame mėginyje, skirtingi branduoliai patiria skirtingą Rabi dazni, dėl ko bendra osciliacija issiplauna.

P|1⟩(t) = sin²(ΩR · t / 2) · e-t/T2 Rabi osciliacijos su eksponentiniu gesimu dėl dekoherencijos
III Interaktyvi Rabi osciliacijos simuliacija

Keiskite RD impulso parametrus ir stebėkite, kaip kubito būsena evoliucionuoja laike. Eksperimentuokite su Rabi dazniu ir dekoherencija.

Rabi osciliacijos simuliatorius

Rabi daznis ΩR 5.0 Hz
T2 relaksacijos laikas 2.00 s
Impulso trukmė t 0.00 s
t (s) P(|1⟩)
100%
P(|0⟩)
0%
P(|1⟩)
|0⟩
100%
|1⟩
0%
|0⟩ |1⟩ x
|0⟩

Blocho sfera: magnetizacijos vektoriaus kryptis

Kairė: P(|1⟩) tikimybės osciliacijos laike su gesimu. Desinė: kubito būsena Blocho sferoje. Pastebėkite, kaip vektorius sukasi is |0⟩ polaus per ekvatorių i |1⟩ ir atgal.

IV Impulso kalibravimas

Rabi osciliacijos yra esminis irankis kvantinių vartų kalibravimui. Norint tiksliai isgyvendinti, pavyzdziui, π/2 vartus (Hadamardo tipo), reikia zinoti tikslią RD impulso trukmę, kuri priklauso nuo Rabi dazno:

tπ/2 = π / (2ΩR)     |     tπ = π / ΩR Impulso trukmės π/2 ir π pasūkimams

Impulso trukmės skaičiuoklė

Iveskite ismatuotą Rabi dazni ir gaukite reikalingas impulso trukmes.

Ismatuotas Rabi daznis ΩR 5.0 Hz
Pasukimo kampas Vartai Impulso trukmė
π/4 (45 laipsnių) T vartai 25.00 ms
π/2 (90 laipsnių) Hadamardo tipo / √X 50.00 ms
π (180 laipsnių) Pauli-X (NOT) 100.00 ms
2π (360 laipsnių) Identitetas (pilnas ciklas) 200.00 ms
V Nerezonansinės Rabi osciliacijos

Kai RD impulso daznis ω neatitinka tikslaus Larmor dazno ω₀, osciliacijos vis dar vyksta, bet su pakeistomis savybėmis. Efektyvusis Rabi daznis tampa:

Ωeff = √(ΩR² + Δω²) kur Δω = ω - ω₀ — daznis dekonfigūracija (detuning)

Nerezonansiniu atveju osciliacijos yra greitesnės (didesnis efektyvusis daznis), bet mazesnės amplitudės — sistema niekada nepasiekia visiskai |1⟩ būsenos. Maksimali pereigos tikimybė:

Pmax = ΩR² / Ωeff² = ΩR² / (ΩR² + Δω²) Pmax < 1 kai Δω ≠ 0

Rezonansinio ir nerezonansinio atvejo palyginimas

Daznis dekonfigūracija Δω 0.0 Hz
t (s) P(|1⟩)
5.0 Hz
Ωeff
100%
Pmax

Melyna linija — rezonansinis atvejis (Δω = 0). Geltona linija — nerezonansinis atvejis su pasirinkta Δω. Didėjant dekonfigūracijai, osciliacijos greitėja, bet amplitudė mazėja.

VI Eksperimentinė procedūra

Atlikite siuos zingsnius SpinQ Gemini platformoje, kad ismatuotumėte Rabi osciliacijas.

Paruoskite sistemą

Ijunkite SpinQ Gemini ir paleiskite programinę irangą. Isitikinkite, kad sistema susieta ir veikia stabilia temperatūrą.

Pasirinkite branduolį

Pasirinkite ¹H (protono) kanalą. Nustatykite RD dazni ties Larmor rezonansiniu dazniu.

Nustatykite impulso parametrus

Pasirinkite fiksuotą impulso amplitudę (B₁ stiprumą). Nustatykite kintamos trukmės impulso eksperimentą su zingsniu, pvz., nuo 0 iki 500 μs kas 10 μs.

Paleiskite Rabi eksperimentą

Paleiskite matavimų seriją. Kiekvienam impulso trukmės taskui sistema atliks: inicializaciją → RD impulsą → signalo matavimą. Stebėkite, kaip signalas osciliuoja.

Analizuokite Rabi kreive

Nubręzkite signalo amplitudę priesais impulso trukmę. Pritaikykite sinusinę funkciją su eksponentiniu gesimu. Is pritaikytos kreivės nustatykite Rabi daznį ΩR.

Nustatykite π/2 ir π impulsų trukmes

Is Rabi kreivės nustatykite impulsų trukmes: pirmasis kreivės maksimumas atitinka π/2, o pirmasis nulis po maksimumo — π impulsą. Uzsirasykite sias vertes.

Pakartokite su ¹⁹F branduoliu

Perjunkite i ¹⁹F kanalą ir pakartokite Rabi eksperimentą. Atkreipkite dėmesį, kad Rabi daznis gali skirtis dėl kito giromagnetinio santykio.

Patikrinkite kalibravimą

Naudodami rastas π/2 ir π impulsų trukmes, paleiskite atitinkamus eksperimentus ir patikrinkite, ar gaunami laukiami rezultatai (superpozicija ir bito apvertimas).

VII Savikontrolės klausimai

Patikrinkite savo supratimą

1. Kas yra Rabi daznis?

2. Kokia impulso trukmė reikalinga kubito būsenai apversti is |0⟩ i |1⟩?

3. Kas atsitinka su Rabi osciliacijos amplitudę, kai RD impulso daznis nesutampa su Larmor dazniu?

4. Kodėl Rabi osciliacijos gęsta eksperimente?

5. Kaip padidinti Rabi daznį?

VIII Apibendrinimas

Pagrindinės sąvokos

  • Rabi osciliacijos — periodiskas kubito populiacijos perėjimas tarp |0⟩ ir |1⟩ būsenų veikiant rezonansiniam RD impulsui
  • Rabi daznis ΩR — osciliacijos daznis, proporcingas impulso stiprumui: ΩR = γ · B₁
  • π/2 impulsas — sukuria superpoziciją, atitinka Hadamardo vartus; trukmė t = π/(2ΩR)
  • π impulsas — apverčia kubito būseną, atitinka Pauli-X vartus; trukmė t = π/ΩR
  • Nerezonansinis atvejis — Ωeff = √(ΩR² + Δω²), mazesnė amplitudė
  • Dekoherencija — T2 relaksacija ir lauko nehomogeniskumas sukelia Rabi osciliacijos gesimą
  • Vartu kalibravimas — Rabi eksperimentas nustato tikslias impulsu trukmes kiekvieniems kvantiniams vartams
← 1: BMR reiskiniai Visi eksperimentai 3 eksperimentas →
IX Programinis pavyzdys ir RQmod alternatyva

Rabi osciliacijų eksperimentą galima atkartoti dviem būdais: tiesiogiai Python kodu (skaitmeninė imitacija arba SpinQit kreipinys į spektrometrą) arba aprašant .rqmod modelyje, kurį ResilQ Studio kompiliuoja į ResilQ IR ir pasirinktinai eksportuoja SpinQ/OpenQASM2 darbą.

Pythonas: Rabi osciliacijų skenavimas

Pavyzdys keičia π impulso trukmę nuo 0 iki 2π ir braižo |1⟩ tikimybę — būtent toks skenavimas kalibruoja X ir H vartus SpinQ Gemini.

# rabi_scan.py — Rabi osciliacijų ir π/2, π impulsų kalibravimas
import numpy as np
import matplotlib.pyplot as plt

Omega_R = 2 * np.pi * 1e3     # Rabi dažnis, rad/s (1 kHz)
T2      = 5e-3                # skersinė relaksacija, s
delta   = 0.0                 # nukrypimas nuo rezonanso, rad/s

t = np.linspace(0, 3e-3, 600)
Omega_eff = np.sqrt(Omega_R**2 + delta**2)

# P(|1⟩) — Rabi formulė su gesimu (T2)
P1 = (Omega_R**2 / Omega_eff**2) * np.sin(Omega_eff * t / 2)**2 \
     * np.exp(-t / T2)

# π/2 → Hadamardo analogas; π → Pauli-X
t_pi_2 = np.pi / (2 * Omega_R)
t_pi   = np.pi / Omega_R

print(f"π/2 impulsas: {t_pi_2*1e6:.1f} µs")
print(f"π impulsas:   {t_pi*1e6:.1f} µs")

plt.plot(t * 1e6, P1)
plt.xlabel("Impulso trukmė, µs"); plt.ylabel("P(|1⟩)")
plt.title("Rabi osciliacijos"); plt.grid(True); plt.show()

Alternatyva: .rqmod modelis

Tas pats kalibravimas išreikštas RQmod kalba. Modelyje fiksuojamos kanalų konfigūracijos, impulso trukmės skenavimas ir matavimo bazė; resilq-q emituoja SpinQit darbą bei CBOM įrodymus apie kalibravimo procedūrą.

// lab2-rabi.rqmod
audit "lab2-rabi-calibration" {
  targets = ["spinq-gemini-lab"]

  layers {
    use { runtime, openssl }
  }

  experiment rabi_scan {
    qubit       = q0
    drive       { freq = larmor; amp = 1.0 }
    sweep       { var = pulse_us; from = 0; to = 3000; step = 5 }
    measure     = z_basis
    shots       = 1024
  }

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

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

resilq-q compile lab2-rabi.rqmod --out build/lab2-rabi

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