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.