I Tikslai

Eksperimento tikslai

  • Suprasti T2* ir T2 dekoherencijos laiku skirtuma ir sukinio aido mechanizma
  • Istirti Hahn aido seka ir jos taikyma T2 matavimui
  • Isnagrineti dinaminio atsiejimo sekas (CPMG, XY-4, UDD) ir ju efektyvuma
  • Suprasti skaitmeninio impulsu optimizavimo metodus (SMP ir GRAPE)
  • Isnagrineti robusto impulsu optimizavima ir eksperimentiniu paklaidu kompensavima
  • Atlikti integruotos kvantines technologijos eksperimenta SpinQ Gemini platformoje

II Teorija: Sukinio aidas

T2* ir T2 dekoherencijos laikai

BMR sistemoje, po 90 laipsniu impulso, magnetizacijos vektoriaus precesija sukurianciame magnetiniame lauke nera vienoda visiems sukiniams. Del magnetinio lauko nehomogeniskumo skirtingos sukiniu grupes patiria skirtingus Larmor daznius, todel precesijos fazeje susidaro issklaidymas. Sis procesas vadinamas laisvu indukciniu irimo signalu (FID — Free Induction Decay).

T2* — Efektyvusis dekoherencijos laikas

Apima visus defazavimo saltinius: tikraja T2 relaksacija ir nehomogenini issklaidyma. FID signalas gesta greiciau nei tikroji T2 relaksacija.

T2 — Tikrasis dekoherencijos laikas

Atspindi tik negrisztamaja (stochastine) relaksacija, nesusijusia su statiniu lauku nehomogeniskumu. Imatuojamas naudojant sukinio aido technika.

1/T2* = 1/T2 + 1/T2' T2' — nehomogeninis defazavimo indelas, T2 — tikroji spin-spin relaksacija

Hahn aido seka

Hahn aido seka leidzia atskirti T2 nuo T2* poveikio. Po pradinio 90 laipsniu impulso sukiniai ima precessuoti ir issibarstyt. Po laiko tau pritaikomas 180 laipsniu impulso, kuris atvercia sukiniu faziu tvarka. Po dar vieno tau periodo sukiniai vel susifokusuoja ir susidaro aido signalas.

Hahn aido seka: |0> ——[90x]—— tau ——[180x]—— tau ——{Aidas}—— Laiko asis: t=0 t=tau t=2*tau 90x 180x Aidas

Aido amplitude gesta pagal eksponentini desni:

M(2*tau) = M(0) * exp(-2*tau / T2) Matuojant aido amplitude prie skirtingu tau verciu, galima nustatyti T2

Refokusavimo mechanizmas

180 laipsniu impulsas veikia kaip "laiko apvertimas" statiniu lauku atvilgiu. Sukiniai, kurie precessavo greiciau del didesnes Larmor daznio, po 180 laipsniu impulso atsilieka, ir atvirksciai. Todel po laiko tau visi sukiniai vel susitinka fazeje — susidaro aidas. Taciau stochastine relaksacija (T2) yra negrisztama ir negali buti atversta.

III Sukinio aido simuliatorius

Hahn aido vizualizacija

Stebekite sukiniu defazavima ir refokusavima. Keiskite tau ir T2 vertes, kad pamatytumete aido formavimasi.

Delsos laikas tau 50 us
T2 relaksacijos laikas 300 us
T2* efektyvusis laikas 80 us
FID amplitude t=tau
—
Aido amplitude t=2*tau
—
T2/T2* santykis
—

IV Sukiniu faziu refokusavimo vizualizacija

Sukiniu fazinis issibarstimas ir aidas

Stebekite 12 sukiniu precesija XY plokstumoje. Pries 180 laipsniu impulsa sukiniai issisklaido, po jo — vel susifokusuoja.

Etapas
Pradzia
Faziu sklaida
0.00
Signalo stiprumas
1.000

V Teorija: Dinaminis atsiejimas

Vienas 180 laipsniu impulsas Hahn aido sekoje kompensuoja tik statinius lauku nehomogeniskumus. Taciau realioswe sistemose egzistuoja ir laikui besikeiciantis triuksmas. Dinaminis atsiejimas (DD) naudoja daugeli 180 laipsniu impulsu, kad efektyviau slopintu triuksma ir pratestum koherencijos laika.

CPMG seka

Carr-Purcell-Meiboom-Gill (CPMG) seka kartoja 180 laipsniu Y impulsus kas 2*tau intervalais:

CPMG seka (N kartojimas): |0> —[90x]— tau —[180y]— tau — tau —[180y]— tau — ... —{Aidas}— Bendra trukme: T = 2*N*tau

CPMG seka veikia kaip aukstu dazniu filtras triuksmui — slopina zemo daznio triuksmo komponentes ir leidzia auksto daznio komponentems praeiti.

XY-4 seka

XY-4 seka keicia impulsu faze tarp X ir Y krypciu, taip kompensuojant impulsu klaidas:

XY-4 seka: |0> —[90x]— tau —[180x]— tau —[180y]— tau —[180x]— tau —[180y]— tau —{Aidas}— Impulsu tvarka: X, Y, X, Y (kartojama N kartu)

UDD seka

Uhrig dinaminis atsiejimas (UDD) naudoja netolygiai isdestymus impulsus, optimizuotus konkretiems triuksmo spektrams:

t_j = T * sin^2(pi * j / (2N + 2)), j = 1, 2, ..., N T — bendra sekos trukme, N — impulsu skaicius

UDD yra optimalus kai triuksmo spektras turi ryzku auksto daznio nukirpima, o CPMG geriau veikia su platesniu triuksmo spektru.

VI Dinaminio atsiejimo palyginimas

DD seku efektyvumo palyginimas

Palyginkite Hahn, CPMG, XY-4 ir UDD seku koherencijos islaikyma esant skirtingam triuksmo lygiui.

Impulsu skaicius N 4
Hahn koherencija
—
CPMG koherencija
—
XY-4 koherencija
—
UDD koherencija
—

VII DD impulsu isdestymo vizualizacija

Impulsu laiko pozicijos

Palyginkite kaip CPMG (vienodai isdestyti) ir UDD (netolygiai isdestyti) impulsai pasiskirsto laike.

Impulsu skaicius N 6

VIII Teorija: Skaitmeninis impulsu optimizavimas

Stipriai moduliuojami impulsai (SMP)

Stipriai moduliuojamas impulsas (SMP — Strongly Modulating Pulse) yra sudarytas is keliu nuosekliu standziu impulsu, kur kiekvienas impulsas apibudinamas keturiais parametrais: amplitude R, radijo dazniu omega_rf, faze phi_rf ir trukme t. Optimizuojant siuos parametrus galima pasiekti norimaji unitarini operatoriu su aukstu tikslumu.

U_SMP = U_N * U_(N-1) * ... * U_2 * U_1 kur U_k = exp(-i*(H_0 + H_1(R_k, omega_k, phi_k)) * t_k)

Tikslo funkcija (fidelity) apibreziama kaip:

F = |Tr(U_tar^dag * U_SMP)| / 2^n F = 1 reiskia tobula atitikima, 2^n — Hilberto erdves dimensija

GRAPE algoritmas

GRAPE (Gradient Ascent Pulse Engineering) dalyja bendra impulso trukme T i N lygiu zingsniu, kiekvienas su fiksuota trukme delta_t = T/N. Kiekviename zingsnyje optimizuojamos tik amplitudes R_k^x ir R_k^y. GRAPE yra efektyvesnis nei SMP, nes turi analizini gradiento sprendini:

dF/dR_k(j) = -2*Re[(P_j^dag * i*delta_t*H_k*X_j) * (X_j^dag * P_j)] kur X_j = U_j...U_1, P_j = U_(j+1)^dag...U_N^dag * U_tar

Robusto impulsu optimizavimas

Realiose sistemose eksperimentines paklaidos (rezonansinio daznio dreifas, RF lauko nehomogeniskumas, magnetinio lauko nehomogeniskumas) mazina impulsu tiksluma. Robusto optimizavimas itraukia sias paklaidas i tikslo funkcija:

F_robust = Sum_i p_i * |Tr(U_tar^dag * U_i)| / 2^n p_i — svoriniai koeficientai, U_i — unitarinis operatorius su i-taja paklaida

Daznio dreifas

Magnetinio lauko pokyciai kecia rezonansini dazni. Nuolatiniai magnetai jautrus temperaturos pokyciams, elektromagnetai — sroves nestabilumui.

RF lauko variacija

RF ritu galios netolygumas erdveje lemia, kad skirtingose vietose sukiniai patiria skirtinga Rabi dazni.

Magnetinio lauko nehomogeniskumas

Skirtingose erdves vietose sukiniai patiria skirtinga magnetini lauka, kas sukelia skirtingus rezonansinius daznius.

IX Impulsu optimizavimo simuliatorius

SMP vs GRAPE palyginimas

Stebakite kaip SMP ir GRAPE optimizavimas randa optimalius impulsu parametrus. Pasirinkite tikslini varta ir stebakite konvergencija.

Impulsu segmentu skaicius (GRAPE) 10
Tikslinis vartai
H
SMP tikslumas
—
GRAPE tikslumas
—
GRAPE segmentai
10

X Robustumo analize

Impulso tikslumas esant paklaidoms

Palyginkite standartinio ir robusto impulso tiksluma esant daznio ir amplitudes paklaidoms.

Paklaidos tipas
Daznis
Standartinis F_avg
—
Robustas F_avg
—

XI SpinQ Gemini realizacija

SpinQ Gemini platformoje sis eksperimentas naudoja dimetilfosfato (CH3O2POH) megini su dviem branduoliais: 1H (vandenilyje) ir 31P (fosfore). Abieju sukiniu kvantinis skaicius yra 1/2.

Parametras Verte Aprasymas
Giromagnetiniai santykiai gamma_H/2pi = 42.58 MHz/T, gamma_P/2pi = 17.23 MHz/T Vandenilio ir fosforo branduoliu giromagnetiniai santykiai
Saistos konstanta J_12 = 696 Hz Spin-spin saisa tarp H-1 ir P-31
90 laipsniu impulsas 40 us (100% amplitude) Standartinio standaus impulso trukme
Rabi dazniai omega_1^H = omega_1^P = 6250 Hz Apskaiciuoti is 90 laipsniu impulso trukmes

Impulso failo formatas

Optimizuoti impulsai irasyti JSON formatu su trimis dalimis: aprasymas (pavadinimas, tikslumas, segmentu skaicius), parametrai (dazniu ir kadro postumiai dviem kanalams), ir impulso duomenys (derinimas, faze, amplitude, plotis kiekvienam segmentui).

{ "description": { "TITLE": "X_gate_robust", "FIDELITY": "0.9998", "TOTALPULSEWIDTH": "320", "SLICES": "8" }, "parameters": { "offset": [{"channel1_pulsefre_offset": 0, ...}] }, "pulse": { "channel1_pulse": [ {"detuning": 0, "phase": 90, "amplitude": 100, "width": 40}, ... ] } }

XII Eksperimento procedura

Sistemos paruosimas

Ijunkite SpinQ Gemini aparata ir palaukite kol sistema pasieks stabilia bukle. Atlikite magnetinio lauko kalibravima ir nustatykite rezonansini dazni abiems kanalams.

Sukinio aido eksperimentas

Ivykdykite Hahn aido seka su keliais skirtingais tau reiksmemis (pvz., 50, 100, 150, 200 us). Uzfiksuokite kiekvieno aido amplitude ir nubraizykite amplitude vs 2*tau grafika T2 nustatymui.

CPMG dinaminis atsiejimas

Ivykdykite CPMG seka su N=2, 4, 8 impulsais, isvengiant vienodos bendros sekos trukmes. Palyginkite koherencijos islaikyma su Hahn aidu.

XY-4 ir UDD sekos

Ivykdykite XY-4 ir UDD sekas su tuo paciu impulsu skaiciumi kaip CPMG. Palyginkite visas tris sekas to paties triuksmo salygomis.

Impulso optimizavimas

Naudojant isoriniu programavimo irankius (Python/MATLAB), atlikite SMP arba GRAPE optimizavima pasirinktam kvantiniam vartui. Eksportuokite optimizuota impulsa JSON formatu.

Optimizuoto impulso ikrovimas ir testavimas

Naudojant "Pulse Load" mygtuka SpinQ sasajoje, ikraukite optimizuota impulsa per USB. Palyginkite standartinio ir optimizuoto impulso fidelity rezultatus.

Robusto impulso tikrinimas

Atlikite robusta optimizavima su 3-5 diskreciais paklaidu taskais. Ikraukite robusta impulsa ir palyginkite jo stabiluma su standartiniu impulsu keiciant eksperimentines salygas.

XIII Savikontrol

1. Ka Hahn aido seka leidzia imatuoti?

Teisingai! Hahn aido seka eliminuoja nehomogenini defazavima (T2' indeli) ir leidzia imatuoti tikraji T2.
Neteisingai. Hahn aido 180 laipsniu impulsas atvercia nehomogenini defazavima, todel aido amplitude priklauso tik nuo T2, ne T2*.

2. Kuo CPMG seka pranasesne uz vienkartini Hahn aida?

Teisingai! CPMG kartojami 180 laipsniu impulsai veikia kaip aukstu dazniu filtras, slopindami zemo daznio triuksmo komponentes ir pratesdami koherencijos laika.
Neteisingai. CPMG pranasumas — daugkartiniai refokusavimo impulsai efektyviau slopina laikui besikeicianti triuksma.

3. Kuo UDD impulsu isdestymas skiriasi nuo CPMG?

Teisingai! UDD pozicijos apskaiciuojamos pagal t_j = T*sin^2(pi*j/(2N+2)), kas optimizuoja apsauga nuo tam tikro triuksmo spektro.
Neteisingai. Pagrindinis UDD skirtumas — netolygus impulsu isdestymas, optimizuotas konkretiems triuksmo spektrams.

4. Koks GRAPE algoritmo privalumas pries SMP?

Teisingai! GRAPE turi analizini gradiento isveda, todel nereikia skaiciuoti baigtinio skirtumo aproksimacijos, kas yra daug efektyviau daugelio sukiniu sistemoms.
Neteisingai. GRAPE privalumas — analizinis gradiento sprendinys, kuris padaro optimizavima efektyvesni, ypac daugelio sukiniu sistemose.

5. Ka reiskia robusto impulsu optimizavimas?

Teisingai! Robusto optimizavimas itaukia i tikslo funkcija skirtingas paklaidu kombinacijas su atitinkamais svoriais: F = Sum p_i*|Tr(U_tar^dag*U_i)|/2^n.
Neteisingai. Robustas optimizavimas reiskia, kad tikslo funkcija apima ne tik idealias salygas, bet ir ivairias eksperimentines paklaidas.

XIV Santrauka

Pagrindiniai rezultatai

  • Sukinio aido technika leidzia atskirti T2 nuo T2*, naudojant 180 laipsniu refokusavimo impulsa statiniam defazavimui kompensuoti
  • Dinaminis atsiejimas (CPMG, XY-4, UDD) pratesia koherencijos laika naudodamas daugkartinius refokusavimo impulsus
  • CPMG naudoja vienodai isdestytus Y impulsus, XY-4 keicia X/Y fazes geresniam klaidu kompensavimui
  • UDD netolygiai isdestyti impulsai optimizuoti konkretiems triuksmo spektrams
  • SMP optimizavimas naudoja keliu standziu impulsu seka su 4 parametrais kiekvienam segmentui
  • GRAPE dalija impulsa i N lygiu segmentu ir naudoja analizini gradiento sprendini efektyviam optimizavimui
  • Robustas optimizavimas itraukia eksperimentines paklaidas (daznio dreifas, RF nehomogeniskumas) i tikslo funkcija

XV Programinis pavyzdys ir RQmod alternatyva

Dinaminio atsiejimo sekas (CPMG, XY-4, UDD) ir GRAPE optimizavimą galima modeliuoti tiek tiesiogiai Python kodu, tiek deklaratyviai .rqmod aprašyme. ResilQ Studio sugeneruoja impulsų tvarkaraščio metaduomenis ir SpinQit darbą.

Pythonas: CPMG koherencijos pailginimas

Pavyzdys lygina laisva FID (T2*) su CPMG seka N refokusavimo impulsų — sumodeliuojama koherencijos kreive priklausanti nuo impulsų skaičiaus.

# cpmg_decoupling.py — CPMG sekos efektyvumas
import numpy as np
import matplotlib.pyplot as plt

T2_star = 0.2     # s — laisvo FID gesimo laikas
T2      = 2.0     # s — vidinis T2 (limit)
t       = np.linspace(0, 3.0, 600)

def cpmg_envelope(t, N):
    """CPMG: efektyvus koherencijos laikas T2_eff ≈ T2 * (N/τ_corr)^β"""
    beta    = 2/3
    T2_eff  = T2_star + (T2 - T2_star) * (1 - np.exp(-N / 10))**beta
    return np.exp(-t / T2_eff)

fig, ax = plt.subplots(figsize=(8, 5))
ax.plot(t, np.exp(-t/T2_star), label="FID (T2*)")
for N in (1, 4, 16, 64):
    ax.plot(t, cpmg_envelope(t, N), label=f"CPMG-{N}")
ax.set_xlabel("t, s"); ax.set_ylabel("|M_xy|"); ax.legend(); plt.show()

Alternatyva: .rqmod modelis

RQmod modelyje sekos aprašomos auksto lygmens primityvais — fiksuojama atsiejimo schema, impulsų skaičius ir optimizavimo metodas. ResilQ Studio sugeneruoja tvarkaraščio metaduomenis ir SpinQit/QASM2 darbą.

// lab10-pulse-engineering.rqmod
audit "lab10-dynamical-decoupling" {
  targets = ["spinq-gemini-lab"]

  layers {
    use { runtime, openssl }
  }

  experiment dd_compare {
    qubit = q0

    sequences {
      hahn  { tau = 100ms }
      cpmg  { n_pulses = [1, 4, 16, 64]; tau = 25ms }
      xy4   { n_blocks = 8; tau = 50ms }
      udd   { order    = 5; total_time = 400ms }
    }
    measure = z_basis; shots = 1024
  }

  experiment grape_optim {
    target_gate = X
    segments    = 32
    optimizer   = grape
    fidelity    { min = 0.999; iters = 200 }
  }

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

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

resilq-q compile lab10-pulse-engineering.rqmod --out build/lab10

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