Laboratorny_praktikum_po_YaMR_Chernov_Butakov_2024
.pdf
С установлением квазиравновесия dndt 0 , откуда можно найти
n |
n0 |
|
|
|
. |
(23) |
|
1 2wT |
|||
|
1 |
|
|
Если 2wT1 1, то n почти не отличается от равновесной величины n0. Поскольку w пропорционально квадрату амплитуды высокочастотного поля, это условие вы-
полняется для малых полей или быстрых релаксационных процессов (малое T1).
С ростом амплитуды переменного поля, когда 2wT1 становится порядка единицы, n начинает заметно падать. При дальнейшем увеличении поля скорость поглоще-
ния энергии стремится к постоянной величине
dE |
ωwn ω |
wn0 |
|
ωn0 |
. |
(24) |
|
1 2wT |
|
||||
dt |
|
|
2T |
|
||
|
|
1 |
|
1 |
|
|
Однако результирующий сигнал, который мы наблюдаем в эксперименте по ЯМР, пропорционален не поглощённой мощности, а избытку населенности на нижнем уровне — n. Поэтому, как это следует из (23), при неограниченном росте напряжённости радиочастотного поля значение n и, следовательно, полезный сиг-
нал, будут уменьшаться до нуля.
4.Поведение системы спинов в постоянном
ипеременном магнитном поле
Вполученном нами квантово-механическом условии резонанса (9) отсутству-
ет постоянная планка . Это указывает на возможность классической интерпрета-
ции явления, при которой ряд характерных особенностей магнитного резонанса удаётся изложить гораздо проще и нагляднее. Поэтому классическая теория ядер-
ного магнитного резонанса наряду с квантовой получила широкое распростране-
ние.
В классической механике доказывается, что изменение момента количества движения должно равняться моменту действующих сил. Применительно к спину
|
|
|
|
|
|
|
|
|
μ |
|
|
||||
с моментом J |
I |
|
|
, находящемуся в постоянном магнитном поле H |
0 |
, это даёт |
|
γ |
|||||||
|
|
|
|
|
|||
|
|
|
|
|
|
11
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dJ |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
μ H0 |
|
|
|
(25) |
||
|
|
|
|
|
dt |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dμ |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
γ μ |
H0 |
. |
|
(25а) |
|||
|
|
|
|
dt |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
есть вектор, направленный пер- |
|||
Поскольку векторное произведение μ H0 |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
пендикулярно плоскости μ |
и |
H0 |
, вектор μ будет описывать конус вокруг |
H0 |
с по- |
|||||||||
стоянным углом θ при вершине. Этот результат можно получить более строго, рас-
писав уравнение (25а) в проекциях по осям:
|
dμ x |
|
γ μ |
|
H |
|
μ |
H |
|
. |
(26) |
|
|
|
y |
z |
y |
||||||
|
dt |
|
|
|
z |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
(Выражения для dμy /dt и dμ z /dt |
аналогичны и отличаются только циклической |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
перестановкой координатных индексов). Если ось z выбрать параллельно H0 |
, то |
||||||||||
Hz = H0 и Hx = Hy = 0. Поэтому из (26) и аналогичных выражений для dμy /dt
и dμ z /dt получим
|
dμ |
x |
γH |
μ |
|
; |
dμ y |
γH |
μ |
|
; |
|
dμ |
z |
0 . |
(27) |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
y |
|
|
|
|
x |
|
|
|
|
|
|||||||||||||||
|
dt |
|
|
|
|
0 |
|
|
|
|
dt |
|
|
|
|
0 |
|
|
|
dt |
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Для μ x отсюда можно найти |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
d2μ |
x |
|
γH |
|
dμ y |
γ2 H2 |
μ |
|
, |
|
|
(28) |
||||||||||||
|
|
|
|
|
dt 2 |
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
dt |
|
|
|
|
|
0 |
|
|
x |
|
|
|
|
|
или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d2μ x |
γ2 H2μ |
|
0 . |
|
|
|
|
|
|
(28а) |
||||||||||||
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
dt 2 |
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Следовательно, μ x совершает гармонические колебания с частотой ω0 γH0 |
по за- |
||||||||||||||||||||||||||||
кону |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
μ x |
Asin ω0 t , |
|
|
|
|
|
|
|
||||||||||||||||
где А и φ — постоянные интегрирования. Аналогично можно получить для μ y |
|||||||||||||||||||||||||||||
|
|
|
|
|
μy |
Acos ω0 t . |
|
|
|
|
|
|
(29) |
||||||||||||||||
12
|
|
|
|
μ2x μ2y , |
|||
Отсюда следует, что проекция вектора μ на плоскость xy, т. е. μ |
|||
остаётся постоянной по величине и вращается с частотой ω0 против часовой стрел-
|
|
|
|
|
|
ки (если смотреть по направлению вектора |
H0 ) |
вокруг оси z . Таким образом, |
|||
|
dμ |
z |
|
|
|
с учётом условия |
|
0 , означающего, что |
z=const, мы видим, что вектор μ |
||
|
|
||||
|
dt |
|
|
||
|
|
|
|
|
|
вращается против часовой стрелки вокруг поля H0 |
с так называемой ларморовской |
||||
частотой 0 = H0, совпадающей с частотой, определяемой формулой (9), получен-
ной из квантово-механического условия резонанса.
Пусть теперь кроме постоянного поля H0 имеется ещё и переменное поле, дей-
ствующее в плоскости, перпендикулярной H0: Hx(t) = Hxmcos t. Это поле может быть представлено как поле, состоящее из двух компонент, вращающихся с частотой ω в
разные стороны. Вблизи резонанса (ω ≈ 0) с магнитным полем будет взаимодейство-
|
|
|
|
вать только компонента магнитного поля, вращающаяся в ту же сторону, что и μ : |
|
||
|
|
|
|
H1 |
i H1cos t jH1sin t , |
(30) |
|
в то время как действием компоненты, вращающейся в противоположную сторону,
можно пренебречь. При этом суммарное магнитное поле будет равно
|
|
|
|
|
|
|
|
|
|
H H1 |
i cos t jsin t kH0 |
, |
(31) |
||
где |
|
— орты координатных осей. Для выяснения действия поля H1(t) удобно |
|||||
i , j, k |
|||||||
ввести систему координат, вращающуюся с частотой в ту же сторону, что и H1(t)
вокруг оси z. В ней вектор H будет покоиться. Из классической механики извест-
но, что скорость изменения вектора во вращающейся системе координат
|
μ |
связана со скоростью изменения этого же вектора в лабораторной системе |
|
|
|
|
dt |
|
координат μ соотношением (вектор угловой скорости направлен в сторону отри- dt
цательного направления оси z)
|
μ |
|
μ |
μ ω . |
(32) |
|
|
|
|||
|
dt |
|
dt |
|
|
13
Если направить ось х вращающейся системы координат (ВСК) вдоль H1 , то вме-
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
сто (30) мы будем иметь H1 |
i H1, и, имея в виду (25а), можно записать: |
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ω |
|
|
|
|
|
|
||||
μ |
ω |
|
|
|
|
|
|
|||||||||||||||||||||
|
|
γ μ H μ ω |
μ |
γ H |
|
|
|
|
μ γ |
H |
|
|
|
|
k H |
i |
|
γ μ H |
|
(33) |
||||||||
|
|
|
|
|
|
|||||||||||||||||||||||
dt |
|
|
|
|
|
|
|
|
γ |
|
|
|
|
|
|
0 |
|
|
γ |
|
1 |
|
|
|
эфф |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ω |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Hэфф |
|
H0 |
|
|
k H1i |
. |
|
|
|
|
|
|
|
(34) |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
γ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Как видно из сравнения (33) с выражением (25а), во вращающейся системе коор-
динат магнитный момент движется так, как если бы на него действовало эффек-
|
|
тивное магнитное поле Hэфф , т. е. он прецессирует вокруг |
Hэфф с угловой частотой |
ωэфф γHэфф (рис. 2). |
|
Рис. 2. Движение магнитного момента ядра μ в постоянном и переменном магнитном полях во вращающейся с частотой ω системе координат: а — когда H0, б — когда = H0
Если частота переменного поля равна ларморовской частоте 0, то, поскольку
|
|
|
|
|
|
|
|
вектор |
ω |
антипараллелен полю H0 |
(см. рис. 2), ω γH0 |
и |
Hэфф |
iH1 . Поэтому |
|
при условии точного резонанса вектор магнитного момента прецессирует вокруг оси х' вращающейся системы координат с частотой
ω1 γH1 .
14
Заметим, что эта частота обычно много меньше 0 = H0, так как H1 имеет по-
рядок единиц эрстед, тогда как H0 104 Э.
Поведение вектора суммарного магнитного момента образца, содержащего
большое число спинов M μi несколько отличается от поведения индивиду-
i
ального спина μ i . Если действует только постоянное поле H0, то нетрудно понять,
суммируя проекции спинов на ось z и на плоскость xy, что величина Mz, пропор-
циональная разнице числа спинов, ориентированных «по» и «против» поля
Mz = kn, остаётся постоянной, в то время как Mx = My = M = 0 в отличие от соот-
ветствующих величин для отдельного спина. Это видно из того, что фазы прецес-
сии отдельных спинов произвольны, следовательно, при большом числе спинов в любой момент времени для любого спина, имеющего определённое направление проекции в плоскости xy, найдётся другой спин, имеющий прямо противополож-
ное направление проекции, лежащей в той же плоскости.
5.Уравнение Блоха
Всвоей первоначальной теории магнитного резонанса Блох исходил из систе-
мы феноменологических уравнений, описывающих поведение компонент суммар-
|
|
|
|
|
|
|
|
|
|
|
|
ного (макроскопического) вектора намагниченности образца |
M . Если каким-либо |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
образом, например путём подачи поля Н1 , вывести |
M из равновесного состояния, |
||||||||||
|
|
|
|
|
|
|
|
|
|||
то после выключения Н1 , как предположил Блох, по аналогии с уравнением (20) |
|||||||||||
|
|
dn |
|
n |
0 n |
, |
|
|
|
||
|
|
dt |
|
T1 |
|
|
|
||||
равновесное значение Mz будет устанавливаться по закону |
|
||||||||||
|
dMz |
|
M0 Mz |
, |
|
(35) |
|||||
|
|
|
|
|
|
||||||
|
dt |
|
|
|
T1 |
|
|
|
|
||
где M0 = 0H0 — равновесное значение намагниченности ( 0 — статическая ядер-
ная магнитная восприимчивость). С другой стороны, равновесное значение компо-
нент Mx и My, как выяснилось выше, равно нулю, поэтому Блох предположил, что эти компоненты стремятся к равновесию с характеристическим временем Т2, 15
которое он назвал временем поперечной релаксации. Соответствующие дифферен-
циальные уравнения запишутся в виде
dMx |
|
|
|
Mx |
, |
|
|
|
|
|
|
||||
dt |
|
|
|
T2 |
(36) |
||
dM y |
|
|
M y |
|
|||
|
. |
|
|||||
|
|
|
|||||
dt |
|
|
|
T2 |
|
||
Учитывая, кроме релаксационных членов (35), (36), движение магнитного момента под действием поля H, можно записать уравнение Блоха во вращающейся системе координат как
|
|
|
|
|
|
|
|
|
|
|
dM |
|
|
i Mx jMy |
Mz |
|
|
||||
γ M H |
|
|
k |
M0 |
. |
(37) |
||||
|
эфф |
|
|
|
||||||
dt |
|
|
|
T2 |
|
|
T1 |
|
||
|
|
|
|
|
|
|
||||
Решения этого уравнения для отдельных компонент момента в ВСК при опреде-
лённых условиях могут быть записаны в виде
Mx |
|
ΔωγH T2 |
|
|
|
(38а) |
||
1 T Δω 2 γ2 H2 T T M0 ; |
||||||||
|
|
|
1 |
2 |
|
|
|
|
|
2 |
|
1 |
1 |
2 |
|
|
|
My |
|
γH1T2 |
|
|
|
M0 ; |
(38б) |
|
|
1 T Δω 2 |
γ2 H2T T |
||||||
|
2 |
|
1 |
1 |
2 |
|
|
|
Mz |
|
1 T2 2 |
|
|
M0 ; |
(38в) |
||
1 T Δω 2 |
γ2 H2 T T |
|||||||
|
2 |
|
1 |
1 |
2 |
|
|
|
где = — 0.
На рис. 3 представлены графики величин My и Mx в зависимости от .
Рис. 3. Кривые поглощения (My ) и дисперсии (Mx ) в соответствии с уравнениями Блоха
16
В лабораторной системе координат компоненты My и Mx вращаются вокруг оси z с угловой частотой . Следовательно, если, как это обычно делается в экспе-
риментах по ЯМР, установить приёмную катушку в плоскости xy, то в ней будет наводиться ЭДС. При этом в зависимости от сдвига фаз между переменным полем
H1 и ЭДС, наводимой в катушке, можно наблюдать сигнал поглощения, пропорцио-
нальный My , или дисперсии, пропорциональный Mx . Вместо My для характе-
ристики поглощения часто пользуются так называемой нормализованной функцией формы линии f( ) kMy . По определению
6.Спектры ЯМР и их наблюдение в жидкостях и твёрдых телах
Втвёрдом теле на ядра помимо приложенного поля H0 действуют локальные магнитные поля соседних ядер и электронов Hлок. В диамагнитных веществах зна-
чительными оказываются вклады соседних магнитных ядер, а в парамагнитных и ферромагнитных — более значительным оказывается вклад электронов. Если для оценки взять расстояние между ядрами 1A0 и считать магнитный момент ядра равным одному ядерному магнетону, то
H |
лок |
μr 3 |
5 Э. |
(39) |
|
|
|
|
Это средняя величина, которая определяет порядок ширины линии в твёрдых
диамагнитных телах. Поскольку в жидкостях в результате быстрого движения мо-
лекул локальные поля усредняются почти до нуля, линии ЯМР в них очень узкие.
Измеренная во многих жидкостях ширина линии оказываются меньше 0,05 Гц или
10-5 Э.
При стандартной нормировке функции формы спектральной линии |
f( ), |
а именно: |
|
|
|
f (ω)dω 1, |
(40) |
|
|
17
получаем
f(ω) |
2Т2 |
|
|
. |
(41) |
|
1 4π2 Т |
(ω |
0 |
ω)2 |
|||
|
2 |
|
|
|
|
|
Функция (41) во многих случаях правильно отражает форму одиночной спектральной линии магнитного резонанса. Из формулы (41) можно вычислить ширину спектральной линии магнитного резонанса (исходя из стандартного опре-
деления Δω1 как интервала частот между точками, в которых функция f() падает до половины от своего максимального значения):
Δω |
1 |
. |
(42) |
|
|||
1 |
πT2 |
|
|
|
|
||
Врезультате различных взаимодействий ядер или электронов с окружением
имежду собой спектры ядерного магнитного и электронного парамагнитного ре-
зонансов расщепляются на некоторое (иногда очень большое) число линий. Оче-
видно, что структуру спектров магнитного резонанса можно наблюдать, если ши-
рина отдельных линий не превышает величины расщеплений в спектре.
Как видно из формулы (42), ширина спектральной линии магнитного резонан-
са определяется временем поперечной релаксации T2 , которое учитывает все фак-
торы, приводящие к затуханию поперечной относительно вектора H0 компоненты спиновой намагниченности (быстрое затухание означает потерю когерентности излучения и, как следствие, более широкое спектральное распределение излуче-
ния). Если отсутствуют аппаратурные причины, то ширина линии магнитного ре-
зонанса определяется временем спин-спиновой релаксации T2. Основной аппара-
турной причиной потери разрешения спектров является существование неодно-
родности статического магнитного поля H0. Действительно, при наличии разброса значений индукции магнитного поля в пределах образца на величину H0 получа-
ем «размывание» частот:
Δω2 γΔH0 . |
(43) |
Учитывая, что суммарная ширина линии даётся выражением Δω
18
Δω Δωi ,
i
где i нумерует каждый отдельный механизм уширения линии, и, принимая во вни-
мание формулы (42) и (43), имеем:
1 |
|
1 |
|
1 |
|
1 |
|
γΔH0 |
, |
(44) |
|
T |
T |
T |
T |
2 |
|||||||
|
|
|
|
|
|
||||||
2 |
|
2 |
|
2 |
|
2 |
|
|
|
|
если присутствуют только две причины уширения линии: спин-спиновая релакса-
ция и неоднородность статического магнитного поля H0.
7. Спектрометр для регистрации ядерного магнитного резонанса
Вывод о существовании прецессии вектора ядерной макроскопической намаг-
ниченности M, подсказывает схему устройства для регистрации ядерного магнит-
ного резонанса. Очевидно, что если охватить образец, помещённый в постоянное
магнитное поле, катушкой индуктивности, ось которой перпендикулярна вектору
H0 , то прецессирующая намагниченность наведёт на концах катушки ЭДС, равную
E n Φ ,t
где n — число витков катушки; Ф — магнитный поток через площадь витков катуш-
ки, который связан с осциллирующей намагниченностью. По порядку величины сиг-
налы ядерного магнитного резонанса близки для всех способов регистрации и могут быть фиксированы в благоприятных случаях радиоприёмником среднего класса.
Рис. 4. Блок-схема экспериментальной установки
для регистрации ядерного магнитного резонанса
19
На рис. 4 показана блок-схема экспериментальной установки для регистрации ядерного магнитного резонанса. Для увеличения сигнала ядерного магнитного резонанса на входе приёмного устройства (усилителя) катушку индуктивности, охватывающую образец, включают в состав радиочастотного контура, настроенного на частоту резонанса 0. Кроме упомянутых элементов и источника магнитного поля,
прибор должен включать какое-то устройство для выведения вектора M из равно-
весного состояния вдоль вектора поля H0 . Обычно эту функцию выполняет ра-
диочастотный генератор (Г), работающий на частоте 0 в стационарном или импульсном режиме. Выход приёмника соединяется с блоком регистрации (Р), которым может быть осциллограф, самопишущий потенциометр, аналого-цифровой преобразователь и т. п. Специальное устройство (СУ) управляет режимом работы генератора и регистрирующими цепями.
8.Импульсный способ регистрации ядерного магнитного резонанса
Взависимости от метода приложения переменного магнитного поля различают стационарный и импульсный способы регистрации магнитного резонанса.
Внастоящее время основной экспериментальной методикой наблюдения ядерного магнитного резонанса является импульсный способ, который состоит в том, что на образец воздействуют переменным магнитным полем в виде коротких радиоимпульсов и после них (или между ними) регистрируют сигналы ЯМР. Первоначально импульсная методика использовалась лишь для измерения времен релаксации и коэффициентов диффузии ядер, но в середине 70-х годов широко распространился способ получения спектров ЯМР из импульсных сигналов с помощью Фурье-преобразования. В настоящее время импульсная методика применяется и для регистрации сигнала электронного парамагнитного резонанса.
Рассмотрим подробнее способ регистрации сигнала ЯМР после одиночного импульса.
Пусть в начальный момент ядерная намагниченность находится в равновес-
ном состоянии, то есть ориентирована вдоль вектора поля H0 и имеет значение
М0. Воздействуем на образец радиочастотным полем H1 в течение короткого
20
