Difference between revisions of "Overtone cp.m"

From Spinach Documentation Wiki
Jump to: navigation, search
(Created page with "Cross-polarization overtone experiment. Syntax: spectrum=overtone_cp(spin_system,parameters,H,R,K) where H is the hamiltonian commutation superoperator, R is the ...")
 
(Update function See also links and function index membership)
 
(13 intermediate revisions by 2 users not shown)
Line 1: Line 1:
−
Cross-polarization overtone experiment. Syntax:
+
{{DISPLAYTITLE:overtone_cp.m}} __NOTOC__
 +
Overtone cross-polarisation experiment that applies a cross-polarisation pulse pair and then performs frequency-domain acquisition at the overtone frequency; because time-domain overtone spectroscopy is difficult (see http://dx.doi.org/10.1039/C4CP03994G), this acquisition mode is preferable in practice, and simulations [[assume.m|assumptions]] should be set to 'qnmr'.
  
−
          spectrum=overtone_cp(spin_system,parameters,H,R,K)
+
==Syntax==
  
−
where H is the hamiltonian commutation superoperator, R is the relaxation
+
    spectrum=overtone_cp(spin_system,parameters,H,R,K)
−
superoperator and K is the chemical kinetics superoperator. The following
 
−
parameters are required:
 
  
−
    parameters.sweep          vector with two elements giving
+
==Parameters==
−
                              the spectrum frequency extents
 
−
                              in Hz around the overtone frequency
 
  
−
    parameters.npoints         number of points in the spectrum
+
    parameters.sweep        -  vector with two elements giving the spectrum frequency extents
 +
                                in Hz around the overtone frequency
 +
 +
    parameters.npoints     -  number of points in the spectrum
 +
 +
    parameters.spins        -  overtone-active nucleus, specified as a
 +
                                single-element cell array
 +
 +
    parameters.spc_dim      -  Fokker-Planck spatial dimension
 +
 +
    parameters.method      -  pulse simulation method, either
 +
                                'average' or 'fplanck'
 +
 +
    parameters.rho0        -  initial state
 +
 +
    parameters.coil        -  detection state
 +
 +
    parameters.Nx          -  X Zeeman operator on the quadrupolar nucleus
 +
 +
    parameters.Hx          -  X Zeeman operator on the spin-1/2 nucleus
 +
 +
    parameters.rf_frq      -  spin-lock frequency offset from the overtone frequency on the quadrupolar nucleus, Hz
 +
 +
    parameters.rf_pwr      -  a vector of spin-lock powers on the quadrupolar nucleus (first element) and the
 +
                                spin-1/2 nucleus (second element), rad/s
 +
 +
    parameters.rf_dur      -  spin-lock pulse duration, seconds
 +
 +
    parameters.method      -  'average' uses the average Hamiltonian theory, 'fplanck' uses Fokker-Planck formalism
 +
                                for the calculation of the pulse evolution.
 +
 +
    H                      -  Hamiltonian commutation superoperator
 +
 +
    R                      -  unthermalised relaxation superoperator
 +
 +
    K                      -  chemical kinetics superoperator
  
−
    parameters.rho0            initial state
+
==Returns==
 +
The function returns the resulting spectrum.
  
−
    parameters.coil            detection state
+
==Examples==
 +
The following 15N overtone spectrum is produced by examples/nmr_overtone/cpmas_glycine_plain.m example file:
  
−
    parameters.Nx
+
[[File:ot_example_3.png]]
−
    parameters.Ny              X and Y Zeeman operators on the
 
−
                              quadrupolar nucleus
 
  
−
    parameters.Hx
+
==Notes==
−
    parameters.Hy              X and Y Zeeman operators on the
+
# Relaxation must be present in the system dynamics, or the matrix inverse-times-vector operation performed by the frequency domain detection module would fail to converge. The relaxation superoperator should ''not'' be thermalised.
−
                              spin-1/2 nucleus
+
# Relaxation theory is not applied during the pulses.
 +
# Average Hamiltonian and Fokker-Planck pulses produce signals in different phases. We are trying to figure out why, but it looks non-trivial. The average Hamiltonian theory option is faster.
 +
# Irrespectively of the pulse algorithm option selection, the magic angle spinning is always handled with the Fokker-Planck formalism.
 +
# There's some really weird difference between the experimental data obtained on Bruker and Agilent spectrometers that goes beyond the spinning direction difference. While this is of no direct relevance to the simulations, at the moment Spinach seems to agree with the Agilent output.
  
−
    parameters.rf_frq          spin-lock frequency offset from
+
==See also==
−
                              the overtone frequency on the
+
[[overtone_pa.m]], [[overtone_dante.m]], [[overtone_a.m]], [[slowpass.m]], [[assume.m]], [[Built-in_experiments]]
−
                              quadrupolar nucleus, Hz
 
  
−
    parameters.rf_pwr          a vector of spin-lock powers on
+
''Version 2.9, authors: [[Ilya Kuprov]]''
−
                              the quadrupolar nucleus (first
 
−
                              element) and the spin-1/2 nucleus
 
−
                              (second element), Hz
 
−
 
 
−
    parameters.rf_dur          spin-lock pulse duration, seconds
 
−
 
 
−
Relaxation must be present in the system dynamics, or the matrix
 
−
inversion operation below would fail to converge. The relaxation
 
−
matrix should *not* be thermalized.
 

Latest revision as of 19:39, 6 June 2026

Overtone cross-polarisation experiment that applies a cross-polarisation pulse pair and then performs frequency-domain acquisition at the overtone frequency; because time-domain overtone spectroscopy is difficult (see http://dx.doi.org/10.1039/C4CP03994G), this acquisition mode is preferable in practice, and simulations assumptions should be set to 'qnmr'.

Syntax

    spectrum=overtone_cp(spin_system,parameters,H,R,K)

Parameters

    parameters.sweep        -  vector with two elements giving the spectrum frequency extents
                               in Hz around the overtone frequency

    parameters.npoints      -  number of points in the spectrum

    parameters.spins        -  overtone-active nucleus, specified as a
                               single-element cell array

    parameters.spc_dim      -  Fokker-Planck spatial dimension

    parameters.method       -  pulse simulation method, either
                               'average' or 'fplanck'

    parameters.rho0         -  initial state

    parameters.coil         -  detection state

    parameters.Nx           -  X Zeeman operator on the quadrupolar nucleus

    parameters.Hx           -  X Zeeman operator on the spin-1/2 nucleus

    parameters.rf_frq       -  spin-lock frequency offset from the overtone frequency on the quadrupolar nucleus, Hz

    parameters.rf_pwr       -  a vector of spin-lock powers on the quadrupolar nucleus (first element) and the 
                               spin-1/2 nucleus (second element), rad/s

    parameters.rf_dur       -  spin-lock pulse duration, seconds

    parameters.method       -  'average' uses the average Hamiltonian theory, 'fplanck' uses Fokker-Planck formalism
                               for the calculation of the pulse evolution.

    H                       -  Hamiltonian commutation superoperator

    R                       -  unthermalised relaxation superoperator 

    K                       -  chemical kinetics superoperator

Returns

The function returns the resulting spectrum.

Examples

The following 15N overtone spectrum is produced by examples/nmr_overtone/cpmas_glycine_plain.m example file:

Ot example 3.png

Notes

  1. Relaxation must be present in the system dynamics, or the matrix inverse-times-vector operation performed by the frequency domain detection module would fail to converge. The relaxation superoperator should not be thermalised.
  2. Relaxation theory is not applied during the pulses.
  3. Average Hamiltonian and Fokker-Planck pulses produce signals in different phases. We are trying to figure out why, but it looks non-trivial. The average Hamiltonian theory option is faster.
  4. Irrespectively of the pulse algorithm option selection, the magic angle spinning is always handled with the Fokker-Planck formalism.
  5. There's some really weird difference between the experimental data obtained on Bruker and Agilent spectrometers that goes beyond the spinning direction difference. While this is of no direct relevance to the simulations, at the moment Spinach seems to agree with the Agilent output.

See also

overtone_pa.m, overtone_dante.m, overtone_a.m, slowpass.m, assume.m, Built-in_experiments

Version 2.9, authors: Ilya Kuprov