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.
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.
Aido amplitude gesta pagal eksponentini desni:
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.
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.
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 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:
UDD seka
Uhrig dinaminis atsiejimas (UDD) naudoja netolygiai isdestymus impulsus, optimizuotus konkretiems triuksmo spektrams:
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.
VII DD impulsu isdestymo vizualizacija
Impulsu laiko pozicijos
Palyginkite kaip CPMG (vienodai isdestyti) ir UDD (netolygiai isdestyti) impulsai pasiskirsto laike.
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.
Tikslo funkcija (fidelity) apibreziama kaip:
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:
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:
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.
X Robustumo analize
Impulso tikslumas esant paklaidoms
Palyginkite standartinio ir robusto impulso tiksluma esant daznio ir amplitudes paklaidoms.
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).
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?
2. Kuo CPMG seka pranasesne uz vienkartini Hahn aida?
3. Kuo UDD impulsu isdestymas skiriasi nuo CPMG?
4. Koks GRAPE algoritmo privalumas pries SMP?
5. Ka reiskia robusto impulsu optimizavimas?
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