Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Практическая квантовая радиофизика. Учебное пособие

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
1 Мб
Скачать

Раздел 3. Форма линии излучения. Ширина линии излучения

 

2

 

 

3I (I 1)

2 2 m

1

 

B2

0

 

 

 

 

 

 

(3cos2 jk 1)2 .

(4 )

2

 

4m

 

6

 

 

 

 

j 1 k

rjk

 

В этом выражении индексы k и j относятся к ядрам, на которых наблюдается резонанс; m — число структурно неэквивалентных ядер, в данном случае m = 1; rjk — значение радиусвектора, соединяющего ядра; jk — угол между приложенным полем и вектором rjk.

Из формулы следует, что с увеличением расстояния вклад от удаленных ядер будет уменьшаться. Атомные ядра, расположенные в кристалле по соседству с исследуемым ядром, можно сгруппировать по слоям, каждый из которых характеризуется определенным значением r. Для каждого такого слоя будет существовать одно или несколько значений направляющих косинусов, причем каждый из них может отвечать нескольким ядрам.

Воспользуемся таблицей для проведения вычислений второго момента.

 

 

 

 

Таблица 2 [6]

 

 

 

 

 

 

Оболочка

(r/a)2

cos2

f (r, )

Число атомов

 

1

1

0

1,0

4

 

1

4,0

2

 

 

 

 

2

2

0

0,125

4

 

1/2

0,0313

8

 

 

 

 

3

3

1/3

0

8

 

4

4

0

0,0156

4

 

1

0,625

2

 

 

 

 

 

 

0

0,008

8

 

5

5

1/5

0,0013

8

 

 

 

4/5

0,0157

8

 

6

6

1/6

0,0012

16

 

2/3

0,0046

8

 

 

 

 

31

Практическая квантовая радиофизика

Здесь f (r, )

a6

(3cos2 jk 1)2 .

r6

 

 

Суммируя по выбранным ядрам, получаем:

 

2

 

9,9I (I 1)

2 2

 

B2

0

 

 

 

a 6 ,

4 2

4

 

 

 

 

 

подставляя значения I = 1/2, = 2,67

108 рад

с–1 Т–1. Окон-

чательно получаем для монокристалла:

 

< B2> = 13,9 10–8 Т2.

Тогда полная ширина линии для монокристалла

B 2,36 B2 8,8 10–4 Т.

В случае поликристаллических материалов следует учесть, что диполи ориентированы случайным образом. Усреднение по всем углам дает

 

 

1

2

 

4

 

(3cos2 1)2

d (3cos2 1)2 sin d

.

4

5

 

 

0

0

 

 

 

 

 

 

В этом случае формула для второго момента имеет вид

B2

2

 

 

3I (I 1) 2 2

 

1

 

0

 

 

 

 

.

(4 )

2

5

6

 

 

jk

rjk

Подставляя численные значения, находим для второго момента:

< B2> = 7,1 10–8 Т2.

Тогда полная ширина линии для поликристалла:

B 2,36 B2 6,3 10–4 Т.

Грубая оценка для ширины линии может быть получена, если предположить, что ширина линии обусловлена локальными магнитными полями, создаваемыми одним диполем в точке расположения второго диполя, B 2Bлок . Причем

32

Раздел 3. Форма линии излучения. Ширина линии излучения

 

Bлок

0

 

[I (I

1)]1/ 2

.

 

4

 

 

 

a3

Вычисления дают: B = 2,2 10–4 Т.

Задача

3.7. Ширина

линии

магнитного резонанса 23Na

(I = 3/2) в металлическом натрии при увеличении температуры

испытывает резкое изменение в

 

 

области 200 К (рис. 4.1). Рас-

 

 

считать коэффициент диффу-

 

 

зии натрия, полагая, что темпе-

 

 

ратурная

зависимость

дается

 

 

выражением

D D0 exp( ED / RT ) ,

где ED — энергия активации, связанная с диффузией; R — газовая постоянная. Принять, что время релаксации для диффу-

Рис. 4.1. Ширина линии 23Na ядерного магнитного резонанса в металлическом натрии

зии равно: C r2 , где r — расстояние между ближайшими

8D

атомами натрия, равное 3,7 10–10 м [4].

Решение.

Ширина линии в магнитном резонансе определяется временем поперечной или спин-спиновой релаксации Т2.

1

2[ B2

 

 

B2

 

 

 

0

 

 

],

T

 

1

 

2 2

Z

 

0

Y

 

 

2

 

 

 

 

 

 

 

0

0

 

где ВY,Z — компоненты флюктуирующего магнитного поля, вызывающего уширение линии, 0 — время корреляции молекулярного движения. Можно выделить два вклада в ширину линии магнитного резонанса — секулярное и несекулярное уширение. Секулярное уширение вызвано неоднородностью внешнего магнитного поля:

33

Практическая квантовая радиофизика

1

2 B2

 

0

.

 

T2

Z

 

 

 

 

 

 

Несекулярное уширение обусловлено конечностью времени жизни. На основании соотношения неопределенности можно записать:

E 1 ,

T1

где

1

2[ B2

B2

]

0

 

.

T

1 2 2

 

X

Y

 

 

 

1

 

 

 

0

0

 

Тогда, считая флюктуирующее поле изотропным, уширение, обусловленное обоими факторами, будет определяться соотношением

1 1' 1 .

T2 T2 2T1

Из приведенных соотношений видно, что при уменьшении времени корреляции молекулярного движения 0 время релак-

сации Т2 возрастает и линия сужается (случай быстрого движения в невязких жидкостях). С другой стороны, при возрас-

тании 0 до значений 0 T2 возникает ситуация, при которой

расфазировка спинов наступает раньше, чем происходит скачкообразное изменение поля. Поэтому расфазировку уже нельзя считать случайным процессом. В этом случае ширина линии не будет зависеть от частоты скачков и, следовательно, от температуры. Это случай жесткой дипольной решетки в твердом теле, который реализуется при низких температурах. На графике можно видеть, что в интервале температур от 0 до 170 К ширина линии постоянна. Оценим значение ширины линии, предполагая, что она полностью определяется локальным полем ближайших соседей:

Bлок

0

 

[I (I 1)]1/ 2

.

4

 

 

 

a3

34

Раздел 3. Форма линии излучения. Ширина линии излучения

Тогда 2 Bлок = 3 103 рад/с. Из графика видно, что это значение несколько меньше экспериментального. При температуре 170 К < T < 220 К диффузия атомов натрия сквозь решетку вызывает флуктуирующие локальные магнитные поля, что приводит к сужению линии. Критическое значение дости-

гается при 0 T2 , т. е. когда время корреляции молекулярного

движения становится сравнимо с обратной шириной диполь-

ной линии. С уменьшением времени корреляции 0 локальные дипольные поля усредняются и ширина линии еще больше сужается, что можно наблюдать при температурах T > 220 К.

Принимая, что

0

T

 

r2

и D D exp( E

D

/ RT ) , со-

 

 

2

 

8D

0

 

гласно условию задачи, и предполагая, что механизм сужения

линии начинает действовать с некоторой температуры Tn, по-

лучим:

 

 

 

 

 

r0

 

 

D

 

exp(E

D

/ RT ) .

 

0

8T2

n

 

 

 

Кроме того, считая, что линия имеет гауссову форму, имеем:

T

 

2

. Здесь

d d

— ширина линии, обусловленная

 

 

2

d d

 

 

 

 

 

 

 

диполь-дипольными взаимодействиями.

Для дальнейших вычислений необходимо определить Tn. Продолжая кривую в сторону низких температур и считая ее ход линейным, найдем точку пересечения с горизонтальным участком кривой. Это дает следующее значение для температуры Tn = 175 К, и при этом ширина линии (по графику)

d d = 17 103 рад/с. Вычисления дают следующие значения

для коэффициента диффузии D0 = 3,0 10–3 м2/с. Экспериментальное значение, полученное с помощью меченных атомов,

D0 = 2,42 10–5 м2/с. Это свидетельствует об очень грубой модели для определения температуры Tn. Кроме того, исследования показывают, что форма линии изменяется от гауссовой до лоренцевой при высоких температурах.

35

Практическая квантовая радиофизика

Задача 3.8. Линия магнитного резонанса ядер со спином I > 1/2 уширена вследствие взаимодействия квадрупольного момента ядра с градиентом электрического поля. Для одноквантовых переходов величина уширения описывается выражением

Q 3eQ(2m 1) VZZ , 4J (2J 1)

где Q — электрический квадрупольный момент ядра; <VZZ> — усредненный по ядрам градиент электрического поля. Такой эффект наблюдался для резонанса на ядрах 115In (I = 9/2) в InSb в присутствии примеси Te. Увеличение концентрации Te контролировалось путем измерения P максимума производной от амплитуды линии поглощения по частоте.

Оценить величину <VZZ>, если предположить, что при концентрации Te, равной n = 1024 м–3, максимум сигнала производной уменьшается в три раза для перехода 3/2 — 1/2 в образце без примеси. Считать, что форма линии — гауссова. Квадрупольный момент Q = 1,2 10–28 Кл м2, ширина линии в беспримесном образце = 4,6 103 рад/с, диэлектрическая проницаемость InSb ε = 15,9 [6].

Решение.

Амплитуда линии поглощения гауссовой формы на некоторой частоте ω относительно резонансной ω0 задается соот-

ношением A aT2 exp( 2T22 ) . Тогда производная будет равна:

 

 

 

 

 

A

2 aT 3 exp( 2T 2 ) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Наибольшее

значение

производной

достигается, когда

2 A

0

при (2T 2 ) 1/ 2.

Тогда P (

A

)

 

 

2aT2 exp( 1/ 2) .

2

 

 

 

2

 

 

 

 

max

 

2

Следовательно, ширина гауссовой линии в присутствии нескольких механизмов уширения равна:

36

Раздел 3. Форма линии излучения. Ширина линии излучения

( )2 ( 0 )2 ( Q )2 ,

где 0 — ширина линии в нелегированном образце, Q — ширина линии в присутствии примеси. Согласно условию задачи, величина Р уменьшилась в три раза по сравнению со

значением в нелегированном образце. Тогда величина ( )2

должна увеличиться в три раза по сравнению с ( 0)2. Следо-

вательно, ( )2 = 3( 0)2. Теперь можно оценить вклад, обусловленный изменением градиента электрического поля

( Q)2

= 2( 0)2. По условию задачи

= 4,6 103 рад/с, то-

гда

 

Q

= 6,5 103рад/с. Подставляя

значения, получим

 

 

 

eQ

V

, что позволяет рассчитать градиент поля

 

48

 

Q

 

ZZ

 

 

<VZZ> = 1,7 1018

В/м2. Рассчитаем величину градиента поля,

используя модель точечного заряда. Потенциал, создаваемый

ионом Te на расстоянии r, равен: V

 

 

e

 

. Тогда градиент

4 0r

 

 

 

 

 

 

вдоль оси

Z можно оценить: V

 

2V

 

 

e

. При кон-

r2

2 0r3

 

ZZ

 

 

центрации Te 1024 м–3 можно принять, что <r3> = 10–24 м3. То-

гда вычисленное значение равно <VZZ> = 1,8 1014 В/м2. Это значительно расходится с полученным экспериментальным результатом. Следует усовершенствовать модель, включив в рассмотрение эффекты, связанные с поляризацией диэлектрика, с поляризацией электронов атомного остова в электрическом поле ядра (антиэкранирующий эффект Стернхеймера).

Задача 3.9. Оценить сдвиг частоты центральной линии магнитного резонанса ядра 7Li (S = 3/2) вследствие взаимодействия квадрупольного момента ядра с градиентом электриче-

ского поля. Принять, что e2Qqzz / h 582 103 Гц, угол между

направлением магнитного поля и осью z, связанной с одной из главных осей тензора градиента электрического поля, β = 450 и 900. Ларморовская частота прецессии лития 12 МГц.

37

Практическая квантовая радиофизика

Решение.

Линия магнитного резонанса ядер со спином S > 1/2 не только уширяется, но и испытывает частотный сдвиг из-за взаимодействия квадрупольного момента ядра с градиентом электрического поля. Рассмотрим ядро со спином S > 1/2 в сильном магнитном поле В0, таком, что взаимодействия квадрупольного момента ядра с градиентом электрического поля можно рассматривать как возмущение. Тогда динамический частотной сдвиг, обусловленный квадрупольными взаимодействиями, во втором порядке теории возмущения можно рассчитать по формуле

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2)

 

 

 

 

VnmVmn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

En

 

 

 

 

 

 

 

 

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2)

 

 

 

(2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

En

Em

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

HQ

 

 

2

 

 

 

n

 

HQ

 

n 1

2

 

 

n 1

 

HQ

 

n 2

2

 

 

 

 

 

 

 

 

 

 

(2)

 

 

 

 

2 n

n 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n n 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q(0)

 

 

 

 

 

 

 

 

 

 

Q(0)

 

 

 

 

 

 

 

 

 

 

Q(0)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

H

Q

 

n

2

2

 

 

n 1

 

H

Q

 

n 1

2

 

 

 

 

n

2

 

H

Q

 

 

n

2

 

 

 

n 1

 

H

Q

 

n 3

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

2 (0)

 

 

 

2 (0)

 

 

 

 

 

 

2 (0)

 

 

 

 

 

 

 

 

2 (0)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q

 

 

 

 

 

 

 

 

 

 

Q

 

 

 

 

 

 

 

 

 

 

 

Q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 2

Здесь

 

n(2) n 1 (En(2) En(2)1) /

 

 

 

 

частота

между

 

зееманов-

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

возмущения; Q(0) (En(0) En(0)1) / — частота между зееманов-

скими уровнями энергии в отсутствии квадрупольного взаимодействия; HQ — квадрупольный гамильтониан в лабораторной системе координат. В ЯМР лабораторная система координат связана с направлением магнитного поля и, как правило, ось Z совпадает с направлением внешнего магнитного поля B0. Тогда квадрупольный гамильтониан в лабораторной системе координат имеет вид

38

Раздел 3. Форма линии излучения. Ширина линии излучения

H лQ.с.к.

e2qQ

{(3IZ2 I 2 )(3cos2 1 sin2 cos )

8S(2S 1)

 

 

(3sin cos sin cos cos2 i sin sin2 )e i (I IZ IZ I )

(3sin cos sin cos cos2 i sin sin2 )ei (I IZ IZ I )

12 (3sin2 (1 cos2 )cos 2 i2 cos sin 2 )e i2 I 2

12 (3sin2 (1 cos2 )cos 2 i2 cos sin 2 )ei2 I 2 .

Вычисление матричных элементов с учетом формул

J ,m | ˆz| J ,m m J J m m ,

J ,m | ˆ | J ,m J (J 1) m(m 1) J J m m 1,

J ,m| ˆ 2 | J,m 2 (J m)(J m 1)(J m 1)(J m 2) J J

приводит к следующему результату:

 

 

 

 

 

2

 

 

 

 

 

2 ( 4S(S 1) 24n2 24n 9)

(2)

 

 

Q

 

F ( , ,

 

 

 

 

 

 

 

n n 1

 

 

 

2 Q(0)

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

2

 

 

 

 

 

Q

 

 

F ( , ,

(4S(S 1) 12n2 12n 6).

 

 

 

 

 

 

 

4 Q(0)

 

 

2

 

 

 

 

 

 

 

 

 

 

Здесь введена постоянная квадрупольного взаимодействия в

виде Q.

e2qQ

. Первое слагаемое учитывает вклады от

4S(2S 1)

 

 

одноквантовых переходов, а второе — от двухквантовых пе-

реходов спина S. F1(β,α,η) и F2(β,α,η) — функции, зависящие от параметров решетки.

| F1( , , ) |2

(3sin cos sin cos cos2 )2 ( sin sin2 )2

9sin2 cos2 2 sin2 2 sin4 cos2(2 ) 6 sin2 cos2 cos(2 ).

39

Практическая квантовая радиофизика

Вычисления приводят к следующему результату:

| F2 ( , , ) |2 14 [(3sin2 (1 cos2 )cos 2 )2

4( cos sin 2 )2 ] 14 [9sin4 2 (1 cos4 )cos2 2

2 2 sin2 2 2 cos2 sin2 (2 ) 6 (1 cos4 )cos(2 )].

Выражения получены для общего случая. Применим их для частных случаев.

Рассмотрим переход из состояния n = –1/2 к n = 1/2, если спин полуцелый:

(2)

 

 

 

 

 

Q2

 

 

[

 

F ( , ,

 

2

1

 

 

F ( , ,

 

2 ](3 4S(S 1)) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1/ 2 1/ 2

 

2 Q(0)

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пусть спин равен S = 3/2. Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2)

 

 

 

6

Q2

[

 

F ( , ,

 

2

1

 

F ( , ,

 

2 ].

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1/ 2 1/ 2

 

 

 

 

 

 

Q(0)

 

 

 

 

1

 

 

 

 

 

 

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя выражения для функций F1 и F2, получаем:

 

(2)

 

 

6

 

 

Q2

 

[

9

(1

cos

2

)(9cos

2

1)

 

 

 

1/ 2 1/ 2

(0)

8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6 sin2 cos2 (7cos2 1)

 

 

 

 

 

 

 

3

 

 

 

8

9 cos2

 

 

 

 

 

 

 

9 cos4

cos2 2 9 cos2

 

 

 

2 (1

 

cos2

 

 

2

cos2 2 )].

2

 

 

 

 

 

 

8

 

 

 

 

 

 

 

 

 

 

 

 

8

 

 

 

 

 

 

 

 

 

4

 

 

 

 

Если параметр асимметрии равен нулю

(η = 0),

выражение

упрощается:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2)

 

 

 

 

 

 

 

 

 

 

Q2

 

 

 

 

9

 

 

 

 

 

2

 

 

 

 

 

2

 

 

 

 

 

 

1/ 2 1/ 2

6

 

 

 

 

[

 

 

 

(1 cos

 

 

)(9cos

 

1)] .

 

 

Q(0)

8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из полученной формулы следует, что сдвиг уменьшается с увеличением индукции магнитного поля.

40

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]