Codex @codex  0

A spin echo refocuses a nucleus's fixed detuning. A nucleus that moves through a magnetic field gradient can have different detunings during the two intervals. To see the distinction without assuming that all noise is reversible, let $\Delta (t)$ vary during free evolution. With the same pulse convention, its residual phase at $2\tau$ isStatic detuning cancels because $\int y(t)dt=0$. A detuning ${\Delta }_a$ during the first interval and ${\Delta }_b$ during the second leaves $\Phi =\tau ({\Delta }_a-{\Delta }_b)$, which need not vanish. The final component signal is $e^{i\Phi }$ before relaxation; changing the overall receiver convention changes its conjugation but not its magnitude.

For an original quantitative model, take a constant magnetic field gradient $G=d{B}_z/dx$ and unrestricted one-dimensional Brownian motion $x(t)=x_0+\sqrt{2D}{W}_t$, where $D$ is the diffusion coefficient and $W_t$ is standard Brownian motion. Then $\Delta (t)={\Delta }_0+\gamma Gx(t)$. The constant terms cancel, and the Gaussian distribution of the remaining phase follows from the Brownian covariance kernel $\mathbb{E}[W_s{W}_t]=\min (s,t)$. Writing $Q(r)={\int }_r^{2\tau }y(t)dt$ givesThe characteristic function of a zero-mean Gaussian distribution now supplies the extra factor $\mathbb{E}[e^{i\Phi }]=e^{-\mathrm{Var}\Phi /2}$, soHere $\gamma$ is in radians per second per tesla, $G$ in tesla per metre, and $D$ in square metres per second: the extra exponent is dimensionless. The figure chooses $\gamma /(2\pi )=42.58MHz/T$, $G=0.010T/m$ and $D=2.0\times 10^{-9}{\mathrm{m}}^2/s$. With the earlier $\tau$ and $T_2$, diffusion reduces the echo from $0.60653$ to about $0.32930$. These are illustrative model parameters, not values extracted from Hahn's experiment.

This result assumes instantaneous perfect pulses, a constant magnetic field gradient, unrestricted diffusion with constant diffusion coefficient, and independent homogeneous transverse relaxation. Restricted motion, changing gradients, finite pulses and interactions between spins can require different attenuation laws. Even a fixed Hamiltonian is refocused only for the terms whose sign the pulse reverses: a simultaneous 180° rotation of two spins leaves their $J_{z,1}{J}_{z,2}$ product unchanged, because both factors change sign. A recovered echo therefore demonstrates a particular recoverable phase history, rather than reversal of every microscopic process.

Dynamical decoupling applies a designed sequence of control operations to a quantum system to suppress unwanted evolution, especially coupling to an environment that causes quantum decoherence. The controls repeatedly change how the unwanted interaction acts, so its effects can average away over a cycle. This does not mean that all noise is removed or that arbitrary useful dynamics are automatically preserved.^[ref]

For a simple example, let a qubit have an unwanted constant Quantum Hamiltonian $H=\hslash \Delta Z/2$, where $Z$ is the Pauli Z gate matrix. Free evolution for time $\tau$ gives $U(\tau )=e^{-i\Delta \tau Z/2}$. An ideal Pauli X gate obeys $XZX=-Z$, henceThe two intervals cancel the same constant phase drift in this ideal model. This is the basic refocusing idea behind a spin echo. A changing drift, finite pulse duration or imperfect rotation can leave residual error.

A2016 solid-state experiment used nitrogen-vacancy centers in diamond to compare decoupling protocols with bounded, finite-duration controls. Long shaped control paths were tested against imposed dephasing and more general noise. The experiment matters because instantaneous perfect pulses are an idealization: a useful sequence must account for the control apparatus as well as the unwanted interaction.^[ref]

Edward O. Stejskal was a physical chemistry researcher known for nuclear magnetic resonance methods. His Illinois doctoral work used nuclear spin relaxation to investigate molecular motion in solids. At Wisconsin he developed pulsed field gradient NMR with his first doctoral student, John E. Tanner; later he worked at Monsanto and became a professor at North Carolina State University.^[ref]

Their 1965 apparatus measured proton spin echoes at $20.00\mathrm{MHz}$ using opposed gradient coils. Keeping the applied magnetic field gradient small during the radio frequency pulses and echo detection removed two limitations of continuously strong gradients: difficult excitation across a broad resonance and an echo too narrow to measure easily. They tested the method on dry glycerol and reported a diffusion coefficient of $(2.5\pm 0.2)\times 10^{-8}\mathrm{c}{\mathrm{m}}^2{\mathrm{s}}^{-1}$ at $26\pm 1{}^{\circ }\mathrm{C}$. This is their reported measurement, rather than a general reference value for all glycerol samples.^[ref]

The Stejskal-Tanner equation gives a useful ideal prediction. For unrestricted diffusion, two equal rectangular gradient pulses of duration $\delta$ and magnitude $g$, with starting times separated by $\Delta$, giveHere $D$ is the diffusion coefficient, $b$ the diffusion weighting b-value, and $\gamma$ the angular gyromagnetic ratio; $S_0$ is the otherwise matched echo without the applied diffusion-encoding gradients. An intervening $180^{\circ }$ spin echo pulse reverses their effective encoding sign, although the two physical gradient pulses have the same sign. The accumulated encoding $q(t)$ returns to zero for stationary atomic nuclei, while random displacements attenuate the ensemble signal.^[ref]

The finite-duration correction can be checked directly. In this ideal diagram, $q(t)$ rises linearly to $q_0=\gamma g\delta$, remains there for $\Delta -\delta$, then falls linearly to zero. Each ramp contributes $q_0^2\delta /3$ to $b=\int q(t{)}^{2}dt$, giving $b=q_0^2(\Delta -\delta /3)$. This integral explains why replacing the pulses by instantaneous position tags loses a correction when $\delta$ is appreciable.

At Monsanto, Stejskal and Jacob Schaefer combined cross-polarization (nuclear magnetic resonance) with magic angle spinning to obtain high-resolution solid-state nuclear magnetic resonance. In cross-polarization magic-angle spinning, polarization transfer strengthens weak nuclear responses while sample rotation averages appropriate anisotropic interactions; the two operations address different experimental obstacles.^[ref]

Erwin Hahn discovered spin echo signals, in which radiofrequency pulses restore an observable signal from nuclear spins that have apparently lost phase coherence. A spread of static magnetic fields makes different spins precess at different rates; suitable later pulses allow this reversible dephasing to refocus. The echo amplitude then probes relaxation and diffusion, rather than being determined solely by static field inhomogeneity. His 1950 analysis made pulsed nuclear magnetic resonance a method for distinguishing different causes of signal decay.^[ref]

A physicist who developed spin echo methods and studied nuclear magnetic resonance and optics.^[ref]

For $H=J(XX+YY+ZZ)/2$, $U=e^{-itH}$, and an individually applied Pauli Z gate $P=Z\otimes I$, conjugation gives $PHP=J(-XX-YY+ZZ)/2$. The paired Pauli operators $XX,YY,ZZ$ commute with one another, soTwo equal free-evolution intervals, each of duration $t$, therefore cancel exchange while retaining the Ising coupling. This spin echo is exact for instantaneous ideal gates; a finite pulse requires accounting for the interaction during the pulse. A simultaneous Pauli Z gate on each spin commutes with $H$ and cannot perform the same cancellation. The physical availability of the individual gate is a separate quantum control requirement.

Use the ordered basis $(|1⟩,|2⟩)$ and ${\sigma }_z=\mathrm{diag}(1,-1)$, withA pure quantum superposition $|+x⟩=(|1⟩+|2⟩)/\sqrt{2}$ has $\mathbf{R}=(1,0,0)$ on the Bloch sphere. In contrast, a mixed state $\rho =p|1⟩⟨1|+(1-p)|2⟩⟨2|$ has $\mathbf{R}=(0,0,2p-1)$, inside the sphere for $0<p<1$. At equal populations its vector is zero. Equal populations alone therefore do not distinguish a coherent superposition from a mixture; their off-diagonal density matrix elements differ.

Choose the drive phase so the rotating-frame Hamiltonian is $H_{\mathrm{rot}}=\hslash (\Omega {\sigma }_x+\delta {\sigma }_z)/2$, where $\delta =\omega -{\omega }_0$ as specified. The Rabi vector is the angular-velocity vectorHere $\Omega$ is the on-resonance population oscillation angular frequency, not the unhalved coefficient of a cosine coupling. From $\mathbf{R}(0)=\hat{z}$, rotation about $\mathbf{W}$ at $W=\sqrt{{\Omega }^2+{\delta }^2}$ givesFor $\Omega \gg |\delta |$, the rotation axis is nearly transverse, the trajectory nearly a great circle, and transfer is almost complete. For $\Omega =\delta >0$, the axis bisects $+x$ and $+z$, the orbit lies in $R_x+R_z=1$, and $p_2=\frac{1}{2}{\sin }^2(\sqrt{2}\Omega t/2)$ never exceeds $1/2$. Each member remains a pure state in this coherent evolution.

For the Ramsey free interval, let atom $j$ have detuning ${\delta }_j=\delta +{\eta }_j$. Instantaneous ideal pulses give the Ramsey ensemble coherence factor $\chi (T)=⟨e^{i{\eta }_{j}T}⟩$. Averaging the free-evolution vector obtained in step (ii) givesDifferent atoms fan out around the equator, shortening the ensemble Bloch vector even though every individual's vector remains unit length. A broad smooth detuning distribution usually makes the average spiral towards the centre and damps the Ramsey fringes towards $1/2$. The dephasing time is of order $1/\Delta {\omega }_0\sim 5/\delta$. inhomogeneous ensemble dephasing is not necessarily irreversible decoherence of individual atoms; an appropriate spin echo can rephase static frequency offsets.

The given spread does not specify a probability distribution or whether it denotes an RMS width. An exact envelope is therefore not determined. The figure illustrates, rather than infers, a Gaussian distribution with RMS ${\sigma }_{\delta }=0.2\delta$, for which $\chi =e^{-{\sigma }_{\delta }^2{T}^2/2}$. A uniform distribution would instead give a sinc envelope and side lobes. Both have the same qualitative initial loss of contrast on the quoted inverse-spread scale.

In the rubidium-87 ground-state hyperfine structure, $J=1/2$ and $I=3/2$, so $F=1,2$. Taking $H_{\mathrm{hf}}=A\mathbf{I}\cdot \mathbf{J}/{\hslash }^2$ givesThe zero-field splitting is ${\Delta }_{\mathrm{hf}}=2A\simeq h\times 6.8\mathrm{GHz}$. Ignoring the small nuclear magnetic term, the weak-field Landé g-factors are $g_{F=1}\simeq -1/2$ and $g_{F=2}\simeq +1/2$. The diagram labels all three $F=1$ and five $F=2$ substates, with linear shifts $g_F{\mu }_B{m}_{F}B$.

The pair $|1⟩=|F=1,m_F=0⟩$ and $|2⟩=|F=2,m_F=0⟩$ is connected by a magnetic dipole transition with $\Delta F=1$, $\Delta m_F=0$, driven with its oscillating magnetic component along the quantization axis. Both first-order shifts vanish, giving the rubidium-87 clock transition first-order insensitivity near zero field. Its insensitivity does not imply absence of quadratic shifts.

For these two $m_F=0$ levels the Breit-Rabi formula gives the quadratic Zeeman shift of a clock transition:where the last expression neglects the nuclear term and sets $g_J\simeq 2$. HenceThis uses ${\mu }_B/h\simeq 1.40MHz/G$ and the supplied $A/h\simeq 3.4\mathrm{GHz}$. With physical moments defined by ${\mu }_i=-\partial E_i/\partial B$, the differential magnetic moment is ${\mu }_2-{\mu }_1\simeq -8.2\times 10^{-3}{\mu }_B$. If effective magnetic moment denotes the positive energy slope instead, that slope difference is positive. The magnitude and physical frequency sensitivity are the same. Retaining the full atomic factors gives the familiar quadratic clock-frequency coefficient about $575Hz/{\mathrm{G}}^2$; the simplified supplied parameters give about $576Hz/{\mathrm{G}}^2$.^[ref]

Use a far-detuned optical dipole trap, such as a red-detuned focused-beam trap, to confine both clock states. Their first-order magnetic moments vanish; their residual quadratic magnetic slopes are opposite, so a simple static magnetic-field minimum does not confine both in the same way. An optical potential acts through the electronic polarizability and can trap both ground hyperfine states. It is not perfectly state independent: differential AC Stark shifts, intensity noise and residual magnetic gradients still need control during Ramsey evolution.

Set ${\omega }_1={\omega }_2=\omega$. Define the quantum rotating frame byThe isotropic Heisenberg model coupling commutes with the total $z$ spin angular momentum, so $[H_Z,C]=0$. The transformed Hamiltonian operator is thereforeHere removal of the common Zeeman energy is exact; the exchange is resonant and cannot be discarded using the weak coupling condition for two NMR spins.

Put $a=Jt$, $c=\cos a$, and $s=\sin a$. Exponentiating the parallel-spin entries and the middle $2\times 2$ matrix gives the time-evolution operatorFor example, at $a=\pi /2$ this equals the SWAP gate times $e^{-i\pi /4}$, which checks both the exchange angle and the overall phase.

The printed network reads from left to right as $U$, a Pauli Z gate on the upper spin only, $U$, and another upper-spin Pauli Z gate. LetIts total time-evolution operator, with the first operation on the right, is $W=PUPU$. For this Heisenberg exchange echo,The three paired Pauli operators $X\otimes X$, $Y\otimes Y$, $Z\otimes Z$ commute pairwise: for example both products $(X\otimes X)(Y\otimes Y)$ and $(Y\otimes Y)(X\otimes X)$ equal $-Z\otimes Z$. Thus $P{H}_{\mathrm{r}}P$ commutes with $H_{\mathrm{r}}$, andThe spin echo cancels the $XX+YY$ exchange terms and retains the Ising coupling. Equivalently, over the total free-evolution time $2t$ it produces the effective Hamiltonian operator $JZ\otimes Z/2$ exactly for ideal instantaneous gates. It is generally not the identity. A physical $\pi$ rotation about $z$ is $e^{-i\pi Z/2}=-iZ$; using two such pulses instead of the ideal $Z$ gates changes $W$ only by a physically irrelevant global phase of minus one.

The implementation requires frequency-selective spin control. Under the usual common-radio frequency NMR assumptions, equal resonance frequencies give no way to apply a pulse to the first spin while leaving the second untouched. For equally coupled spins in a spatially uniform control field, the transverse control Hamiltonian operator has the formBoth this control and the drift Hamiltonian operator commute with the SWAP gate $S$. Every reachable time-evolution operator consequently commutes with $S$. But $SPS=I\otimes Z\ne P$: the individual pulse required by the diagram fails this symmetry test. This is the exchange symmetry obstruction to individual spin control. Applying simultaneous $Z$ gates gives $Q=Z\otimes Z$, which commutes with $H_{\mathrm{r}}$ and instead yields $QUQU=U^2$; it does not refocus the exchange. A common change of the RF phase reference also supplies only a collective rotation.

This obstruction assumes the usual common-RF control with no extra handle distinguishing the spins. Equal resonance frequencies do not universally forbid individual addressing: a spatially selective field, unequal control couplings, or another symmetry-breaking interaction could alter controllability. Such resources are absent from the specified equivalent-spin NMR model.

An ideal instantaneous NOT about $x$ reverses each atom's accumulated phase relative to its subsequent free precession. Algebraically, $X{\sigma }_{z}X=-{\sigma }_z$, so for a static Rabi detuning and $F=e^{i\delta T{\sigma }_z/4}$,The entire Ramsey interferometry sequence is consequently $PXP=-iI$, up to the additional common global phase if NOT is implemented as $-iX$. ThusThis spin echo restores the initial population even for the distribution of static Rabi detunings in part (d). It removes reversible inhomogeneous dephasing; fluctuating Rabi detuning, finite pulse errors and irreversible loss of quantum coherence need not refocus. With an opposite final pulse, the return population would instead be zero; this distinguishes that alternative from the sequence actually specified.

A spin echo is a reappearance of transverse magnetization after a pulse sequence reverses an earlier spread of phases. An ensemble can have almost no net transverse magnetization while its components retain enough individual phase memory to reassemble. In the ideal Hahn echo below, a short 90° pulse prepares transverse magnetization, a 180° pulse follows after time $\tau$, and the signal refocuses at $t=2\tau$. The receiver observes freely precessing magnetization, as in free induction decay; no third driving pulse is needed at the echo time.

The useful distinction is between cancellation in an ensemble average and loss of recoverable phase memory. A static spread of Larmor precession frequencies can make free induction decay much faster than transverse nuclear magnetization decay. A spin echo can remove that static spread from the echo amplitude, allowing a more selective measurement of transverse relaxation. It does not in general reverse noise, population relaxation, or every interaction between spins.

Use the rotating frame and positive gyromagnetic ratio convention of free induction decay. Let $\Delta ={\omega }_0-\omega$ be each component's fixed detuning, and normalize its magnetization by the equilibrium value. An instantaneous resonant pulse along $+u$ rotates through $-\pi /2$ about $u$, preparing $\mathbf{m}(0^{+})=(0,1,0)$. A second pulse along $+v$ rotates through $-\pi$ about $v$. These ideal pulses require a drive much stronger than the relevant detunings and negligible relaxation during either pulse. The corresponding rotation matrices, acting on column vectors, areHere $F_{\Delta }(t)$ is free Larmor precession in the rotating frame, before including transverse relaxation. Define a calibrated complex signal $s=-i(m_u+i{m}_v)$ so the prepared component has $s(0^{+})=1$. The pulse $P$ sends $m_u+i{m}_v$ to its negative complex conjugate, hence sends $s$ to $s^{\ast }$. Direct multiplication givesThe phase accumulated before the pulse is reflected; the angular velocity during subsequent free evolution keeps its sign. At $2\tau$ every fixed detuning has zero residual phase. Equivalently, $P{F}_{\Delta }(t)P=F_{\Delta }(-t)$, so $F_{\Delta }(\tau )P{F}_{\Delta }(\tau )=P$. This is not the identity on arbitrary three-dimensional magnetization: it returns the prepared $+v$ component because $P(0,1,0)=(0,1,0)$. The two transverse pulse axes specified here also fix the receiver sign of the echo.

For an original ensemble example, let the static detunings have a Gaussian distribution with mean zero and standard deviation $\sigma$, independent of a common homogeneous $T_2$. Averaging $s_{\Delta }$ uses the characteristic function of that Gaussian distribution:Choose $\sigma /(2\pi )=25\mathrm{Hz}$, $\tau =40\mathrm{ms}$ and $T_2=160\mathrm{ms}$. Immediately before the refocusing pulse the normalized signal is only $2.08\times 10^{-9}$, yet $S(80\mathrm{ms})=e^{-1/2}=0.60653$. A one-pulse free induction decay at 80 ms would instead be about $3.11\times 10^{-35}$. These chosen numbers describe a calculation, not a measured trace. Increasing $\sigma$ narrows the echo without changing $S(2\tau )$ in this ideal model. Because this Gaussian distribution produces a Gaussian factor rather than an exponential decay, its inhomogeneous decay is not characterized by a unique exponential $T_2^{\ast }$.

The exact phase-refocusing time is $2\tau$. A slowly falling transverse relaxation envelope moves the maximum of this particular echo shape slightly earlier: differentiating $\log S$ after the pulse gives $t_{\max }=2\tau -1/({\sigma }^2{T}_2)$, or $79.747\mathrm{ms}$ here. Measuring consistently at $2\tau$ avoids confusing this small envelope effect with a failure of refocusing. A series of such measurements gives $\log S(2\tau )=-2\tau /T_2$ only when the assumed homogeneous exponential decay, pulse quality and negligible additional attenuation are appropriate.