Processing math: 100%

Mathematics

The hamiltonian ˆH we use for NMR systems is of the form

ˆH=γBz(1+δ)ˆIz+2πk<lJklˆIkˆIl

with gyromagnetic factors γ, chemical shifts δ of nuclear spin , coupling between spins k and l denoted as Jkl, and ˆIα=ˆSα/, with ˆSα being the usual su(2) spin operators.

Within NMR, we have a strong magnetic field Bz in the z-direction, and electromagnetic pulses / oscillating fields are applied to flip the spins into the x/y plane. Since typically Bz is of the order of 500Mhz, the pulses of 10kHz bandwidth, and the required resolution is sub 1Hz, we refrain from modeling the explicit time dependence of the pulses. Instead, we model the pulses directly by calculating the spectral function, i.e., time-dependent correlations between the corresponding ˆM± operators.

The spectrum measured in an NMR experiment corresponds to the spectral function, which is the Fourier transform of the correlation function of the operators ˆM±, calculated by the quantum program created by the HQS Qorrelator App.

Calculation of the spectral function

The spectral function of the NMR problem is given by the Fourier transform

G(ω)=FTG(t)G(t)=iˆM(t)ˆM+(0)ˆI±=ˆIx±iˆIyˆMα=γˆIα

where the ˆMα operators, αx,y,z,+,, contain the gyromagnetic factors for convenience, and t is the real time dependence.

The contribution of an individual nuclear spin to the full NMR spectrum is obtained via

G(t)=iγˆI(t)ˆM+(0),

while the full NMR signal is the sum of individual contributions.