Реферат: Метод моментов в определении ширины линии магнитного резонанса

упомянутая модель может быть использована лишь, когда теоретическое отношение M₄ /( M₂ )² оказывается большим числом., В этом случае

(IV.25б)

Ширина на половине высоты значительно меньше, чем среднеквадратичная ширина. С другой стороны, предположение о гауссовой форме линии может быть разумным всякий раз, когда отношение M₄ /( M₂ )² порядка 3.

§ 5. МЕТОД ВЫЧИСЛЕНИЯ МОМЕНТОВ

Основной недостаток метода моментов состоит в том, что важный вклад в значение момента (вклад тем существеннее, чем выше момент) дают крылья кривой, которые на практике не наблюдаются. Необходимо из вычисленных моментов линии магнитного резонанса с центром на ларморовской частоте w =w₀ исключить вклады от сопутствующих линий на частотах w = 0, 2w₀ , 3w₀ о которых упоминалось ранее. Легко видеть, что, несмотря на их малую интенсивность (благодаря удаленности от центральной частоты w₀ ) вклад во второй момент сравним с вкладом от главной линии и тем больше, чем выше порядок момента. Для исключения вкладов от них следует рассматри­вать в гамильтониане возмущения ħH₁ ответственного за уширение, только его секулярную часть ħH ¢₀ , которая коммутирует с H₀ и, следова­тельно, не может отвечать перемешиванию состояний с различными пол­ными М; такое смешивание является причиной появления побочных линий. Таким образом, сокращение дипольного гамильтониана до его секулярной части

не только упрощает вычисление моментов, но и делает его более точным.

Прежде чем начать расчет, отметим, что линия магнитного резонанса симметрична относительно центральной частоты w₀ . Убедимся в правиль­ности этого утверждения. Если | а > и | b > — два собственных состояния ħ(H₀ +H ¢₁ ) с разностью энергии ħ(Е_а — Е_b ) = ħw₀ + d_ab , то два состоя­ния | а^~ > и | b^~ >, полученные из | а > и | b > соответственно путем пово­рота всех спинов в обратном направлении, будут также собственными состояниями ħ(H₀ +H ¢₁ ) с ħ(Е_b ~ – Е_a ~) = ħw₀ + d_ab . Таким образом, каждо­му переходу с частотой w₀ + u соответствует переход равной интенсивности с частотой w₀ – u. Если f(w) — функция формы, то h (u) = f(w₀ + u)— четная функция u. Поскольку моменты кривой пропорциональны про­изводным в начале координат от их фурье-преобразования, мы будем применять для их вычисления формулу (13). Вследствие узости линии ядерного магнитного резонанса можно пренебречь изменением величины w в пределах ширины линии и предположить, что форма линии описывается c¢¢(w)/w, так же как и c¢¢(w). Тогда, поскольку f(w) — нормированная функция формы, (13) может быть переписано в виде

f(w) = A∫ G(t) cos wt dt, (IV.26)

где постоянная A определяется из условия нормировки f(w), а опреде­ленная ранее четная функция G (t) равна Sp{M _x (t) M _x }. Обратно

G(t) = 2/(pA)∫ f(w) cos wt dw, (IV.27)

Согласно вышеизложенному, в выражении

M _x (t) = е ^{i H t} M _x е^{– i H t} .

следует вместо H = H₀ +H₁ подставить H = H₀ +H ¢₁ что значи­тельно упрощает вычисления. Поскольку H₀ и H ¢₁ коммутируют, можно записать

exp{i(H₀ +H ¢₁ )t} = exp(iH₀ t) exp(iH ¢₁ t).

Учитывая, что зеемановский гамильтониан ħH₀ равен ħw₀ I_z функцию G (t) можно переписать в виде

(IV.28)

Шпур произведения операторов инвариантен относительно циклической перестановки, поэтому

(IV.28a)

В этом выражении оператор exp(iw₀ I_z t) определяет поворот на угол w₀ t вокруг оси z, и, следовательно, можно записать

(29)

Легко видеть, что второй член в (29) равен нулю, так как поворот спинов на 180°, например вокруг оси ох, не изменяет H ¢₁ и M _x но преоб­разует M _у в – M _y .

Заменяя в (27) G (t) на G₁ (t) cosw₀ t, где

G₁ (t)=Sp{е xp ( i H ‘₁ t ) M _x е(– i H ‘₁ t )M _x }

называется сокращенной функцией автокорреляции, и вводя обозначение

h (u) = f(w₀ + u),

получаем

Заменяя нижний предел на – ¥, что допустимо для узких линий, найдем

Поскольку h (и) является четной функцией, второй интеграл равен нулю и

G₁ (t)=Sp{е xp ( i H ‘₁ t ) M _x е(– i H ‘₁ t )M _x }

(30)

Различные моменты кривой распределения h (и) относительно резонансной частоты w =w₀ определяются выражением

Нечетные моменты равны нулю, а четные определяются формулой

(31)

Таким образом, для вычисления моментов резонансной кривой достаточно разложить G₁ (t) в выражении (30) по степеням t. При этом коэффициенты разложения представляют собой шпуры от операторов, которые являются полиномами от H ¢₁ и M _x .

Сущность метода заключается в том, что значения упомянутых шпуров не зависят от выбора основных состояний и могут быть вычислены, напри­мер, в представлении, где значения m^j = I^j _z отдельных спинов (поэтому представление называется m^j -представлением) являются хорошими кван­товыми числами. Таким образом, нет необходимости решать проблему отыскания собственных состояний | n > полного гамильтониана. Из опре­деления (30) функции G₁ (t) вытекает, что значение ее р-й произ­водной в момент t = 0 определяется выражением

(IV.32)