Simulation of SQ and MQ coherence line intensities of a powder.
Contributor: R. Hajjar

Home and Applets > Quadrupole Interaction > One-Pulse MAS > New Version for Rotating Powder

One pulse applied to MAS powder

AIM: We provide a new Mathematica-5 notebook to simulate single-quantum and multiple-quantum coherence line intensities for MQMAS NMR applied to half-integer quadrupole spin in a powder.

Method: We simulate single-quantum and multiple-quantum coherence line intensities of a spin I = 3/2 with increasing pulse duration in a powder rotating at the magic angle.

The parameters for these simulations are:

(A) Mathematica-5 notebook: onePowderMAS.nb

Get["QUADRUPOLE"];
(*------------- Nucleus ----------*)
quadrupoleSpin = 1.5;
larmorFrequencyMhz = 105.8731007;
(*----- Quadrupole interaction ----*)
quadrupoleOrder = 2; QCCMHz = 8; η = -1;
(*--- Rotor Euler angles in PAS ---*)
αPR = 30; βPR = 80; γPR = 120;
(*----------- Parameters ----------*)
startOperator = 0.4*Iz;
ωRFkHz = 100;
spinRatekHz = 15;
powderFile = "rep100_simp";
numberOfGammaAngles = 3;
t1 = 20;
Δt = 1;
np = t1/Δt;
(*--------- Pulse sequence ---------*)
detectelt = {{3, 2}};
fsimulation := (
  acq0;
  For [p = 1, p <= np, p++, {
      pulse[Δt, ωRFkHz];
      ;acq[p];
  }];
);
(*---Execute, plot, and save simulation 
  in "onePowderMAS" file--------------*)
run;
tabgraph["onePowderMAS"];
	

(1) Preliminary

  1. Download Mathematica-5 notebook onePowderMAS.nb (the corresponding PDF file), that for MAS NMR utilities QUADRUPOLE_1_0.nb (the corresponding PDF file), and the crystal file rep100_simp.
  2. Save these three files into Mathematica-5 folder. Forbidden the Operating System of your computer to include extra file extension to rep100_simp by providing the file name with double quotes such as "rep100_simp".
  3. Open QUADRUPOLE_1_0.nb file with Mathematica-5.
  4. Press "Ctrl-A" to select the notebook, then press "Shift-enter" to start the notebook. (Some warning messages appear but they have no consequences on the results.) A new file called QUADRUPOLE is created in Mathematica-5 folder.

(2) Simulation

  1. Open onePowderMAS.nb file with Mathematica-5.
  2. Press "Ctrl-A" to select the notebook, then press "Shift-enter" to start simulation. (Some warning messages precede the simulation.) At the end a data file called onePowderMAS is created in Mathematica-5 folder. MS Excel can open this data file for graphic representation.

(B) Result

The simulated line intensities for two opposite rotor spinning speeds are gathered in the following table.

t
(μs)
Rotor spinning speed:
15 kHz
Rotor spinning speed:
-15 kHz
0
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
0              
0.1966516113
0.1407515005
-0.07013333075
-0.1418308959
0.004308044472
0.1484301615
0.09519135239
-0.06704005251
-0.1129497492
0.01198907045
0.124670265
0.07777468433
-0.04395928328
-0.07373481123
0.01932615899
0.08592484193
0.03816589776
-0.04505499021
-0.04325721476
0.03423029552
0              
0.1966516113
0.1405010943
-0.07193143925
-0.1464688148
0.008288630193
0.1550069159
0.08993224876
-0.08208032528
-0.121105739
0.02564040147
0.1439525973
0.07957656183
-0.06646607056
-0.08892222642
0.02320356206
0.10851410743
0.05282811978
-0.05239192631
-0.05929601348
0.02762217095

(C) Conclusion

Powder averaging using crystal file (or powderFile) for line intensity simulation provides different results for two opposite rotor spinning speeds. However this discrepancy disappears when the asymmetry parameter is zero, or when one of the two angles α and β is zero in the crystal file (or powderFile).

Solid-state NMR bibliography for:

Aluminum-27
Antimony-121/123
Arsenic-75
Barium-135/137
Beryllium-9
Bismuth-209
Boron-11
Bromine-79/81
Calcium-43
Cesium-133
Chlorine-35/37
Chromium-53
Cobalt-59
Copper-63/65
Deuterium-2
Gallium-69/71
Germanium-73
Gold-197
Hafnium-177/179
Indium-113/115
Iodine-127
Iridium-191/193
Krypton-83
Lanthanum-139
Lithium-7
Magnesium-25
Manganese-55
Mercury-201
Molybdenum-95/97
Neon-21
Nickel-61
Niobium-93
Nitrogen-14
Osmium-189
Oxygen-17
Palladium-105
Potassium-39/41
Rhenium-185/187
Rubidium-85/87
Ruthenium-99/101
Scandium-45
Sodium-23
Strontium-87
Sulfur-33
Tantalum-181
Titanium-47/49
Vanadium-51
Xenon-131
Zinc-67
Zirconium-91
[Contact me] - Last updated August 30, 2020
Copyright © 2002-2022 pascal-man.com. All rights reserved.