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

Ядерный магнитный резонанс и ЯМР-томография (110

..pdf
Скачиваний:
1
Добавлен:
15.11.2022
Размер:
838.71 Кб
Скачать

МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РФ ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ «ВОРОНЕЖСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ»

ЯДЕРНЫЙ МАГНИТНЫЙ РЕЗОНАНС И ЯМР-ТОМОГРАФИЯ

Учебное пособие для вузов

Составитель С. Г. Кадменский

Издательско-полиграфический центр Воронежского государственного университета

2012

Утверждено научно-методическим советом физического факультета 10 января 2012 г., протокол № 1

Рецензент д-р физ.-мат. наук, проф. В. Д. Овсянников

Учебное пособие подготовлено на кафедре ядерной физики физического факультета Воронежского государственного университета.

Рекомендуется для студентов 5-го курса физического факультета.

Для направления 010700 – Физика

2

ВВЕДЕНИЕ

Основные литературные источники, которые обычно используются в качестве учебных пособий при изучении явления ядерного магнитного резонанса (ЯМР) и его применения в исследованиях структуры различных конденсированных сред и в ядерной медицине, весьма ограниченны по числу и труднодоступны. Поэтому данное учебное пособие, ориентированное на анализ фундаментальных свойств ядерного магнетизма и ядерного магнитного резонанса, а также описание наиболее распространенных способов их применения, включая современную медицинскую ЯМР-томографию, может оказаться весьма полезным для широкого круга студентов и специалистов.

1. ИСТОРИЯ РАЗВИТИЯ МЕТОДОВ ЯМР

Развитие методов ядерного магнитного резонанса генетически связано с тем фактом, что ларморовские частоты ω0 ядерных магнитных подуров-

ней, возникающих при помещении атомных ядер во внешнее однородное и стационарное магнитное поле, оказываются лежащими в диапазоне частот, соответствующих частотам ω высокочастотных магнитных полей, которые могут быть созданы при использовании стандартных радиотехнических

средств. Заметим, что соответствующие ларморовские частоты ω0эл для

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

проблематичным экспериментальное получение высокочастотных переменных магнитных полей с частотами ω ω0эл .

Использование методов ЯМР началось достаточно давно. Еще до войны в 1939 году Раби продемонстрировал явление ЯМР при исследовании отклонения пучков нейтронов во внешних магнитных полях для определения дипольного магнитного момента нейтрона. Это явление связано с опрокидыванием вектора спина нейтрона, ориентированного в однородном и

постоянном магнитном поле с индукцией B0 , при резонансном поглощении нейтроном кванта высокочастотного магнитного поля с индукцией BG1, перпендикулярной вектору BG0 . В дальнейшем метод Раби [1] был применен и

для молекулярных пучков, что позволило измерять не только магнитные дипольные моменты, но и спины различных ядер. Наконец, позже метод ЯМР был использован для прецизионного определения напряженностей как внешних, так и внутренних магнитных полей широкого круга конденсированных сред.

3

В 1946 году независимо друг от друга двое ученых из США – Ф. Блох из Стэнфордского университета и Р. Парселл из Гарварда – описали явление ЯМР в конденсированных средах [2, 3], основанное на существенном изменении направления вектора намагниченности системы ядер, находящихся в однородном и постоянном магнитном поле B0 , при действии на эту

систему резонансного высокочастотного магнитного поля сG индукцией B1, ориентированной перпендикулярно направлению вектора B0 . Авторы опи-

сали спин-решеточные и спин-спиновые релаксационные явления, позволившие включить в метод ЯМР измерения соответствующих этим явлениям времен релаксации к равновесному состоянию вектора намагниченности ядерных систем для различных конденсированных сред. Хотя ЯМР-методы высокого разрешения разрабатывались как универсальные средства изучения химического состава и структуры твердых тел и жидкостей, особый интерес к этим методам с медицинской и биологической точки зрения был обусловлен возможностью проведения неинвазивных (неразрушающих) измерений внутренних характеристик биотканей. В 1971 году Дамадьян обнаружил разницу в магнитных свойствах нормальных и опухолевых тканей. В 1973 году Лаутербург, а также Мэнсфилд и Гранелл предложили использовать для исследования объемных характеристик ЯМР сдвиг резонансной частоты ω0 при наложении градиента магнитного поля. В 1974 году Лау-

тербург продемонстрировал явление ЯМР на живой мыши, а в 1977 году Дамадьян впервые получил ЯМР-томографические изображения человеческого тела, опираясь на свойства ядер водорода. Вслед за этим последовали технические и промышленные разработки с использованием компьютеров, что привело к возникновению разнообразных методов, позволяющих получать изображения по распределениям в биотканях не только ядер водорода, но и ядер других химических элементов, дающих информацию о плотностях спинов и временах релаксации. Были разработаны быстрые методы визуализации, включая визуализацию в масштабе реального времени, а также методы по разделению вкладов воды и жиров в формируемом изображении и по измерению кровотока в организме. В настоящее время в мире насчитывается более десяти тысяч ЯМР-томографов, работающих на базе постоянных магнитов как обычного, так и сверхпроводящего типов.

2.ЯДЕРНЫЙ МАГНЕТИЗМ

2.1.Атомное ядро во внешнем магнитном поле

Атомное ядро сGотличным от нуля полным спином I имеет магнитный дипольный момент μ [4]

4

 

 

 

μG =

 

e

gI IG =γI IG,

(1)

 

 

 

 

2mpc

 

 

 

 

 

 

 

где

e

– ядерный магнетон,

gI

– безразмерное гиромагнитное отноше-

2mpc

 

 

 

 

 

 

 

ние для указанного ядра, а γI

– гиромагнитное отношение, имеющее раз-

мерность, совпадающую с размерностью ядерного магнетона. Принято использовать величину МГц/Тл, где 1 мегагерц отвечает частоте колебательного процесса, связанного с 106 колебаний в секунду, а 1 тесла отвечает величине индукции магнитного поля, равной 104 гаусс, где 1 гаусс соответствует средней индукции магнитного поля Земли на ее поверхности. Из-за существования в атомных ядрах протонных и нейтронных куперовских пар, обладающих нулевым спином, спины основных состояний ядер равны нулю для четно-четных ядер, в которых числа протонов Z и нейтронов N являются четными, равны спину неспаренного нуклона для нечетных по атомному весу A = (Z+N) ядер и, наконец, равны векторной сумме спинов неспаренного протона и неспаренного нейтрона для нечетно-нечетных ядер. Отсюда следует, что при нормальных условиях дипольные магнитные моменты отличаются от нуля только для нечетных и нечетно-нечетных ядер. Как видно из табл. 1, величины гиромагнитных отношений γI для различных атомных

ядер, участвующих в формировании биологических тканей, меняются в диапазоне от 10 до 93 МГц/Тл и принимают особенно большие значения для нечетных изотопов водорода 1Н1, фтора 19F9 и калия 39K19.

 

 

 

 

 

 

Таблица 1

Ядра

1H1

13C6

19F9

23Na11

31P15

39K19

 

Гиромагнитное

42,58

10,71

40,05

11,26

17,12

93,1

 

отношение (МГц/Тл)

 

В теории электричества [5] известны следующие основные виды магнитных сред: диамагнетики, парамагнетики и ферромагнетики. Рассмотрим, к какому виду магнетиков относится система тождественных атомных ядер, обладающих отличным от нуля магнитным дипольным моментом. Для этого поместим данную систему во внешнее магнитное поле с вектором маг-

нитной индукции B0 , направленным по оси Z лабораторной системы коор-

динат (л.с.к.). Действие этого поля на атомное ядро с дипольным моментом μG приводит к трем основным эффектам [5].

Первый эффект связан с появлением дополнительной энергии E

рас-

сматриваемого ядра во внешнем магнитном поле

B [5]:

 

G

GG

G0

 

E = −(μGB0 ) = −γI (IB0 ) = −γI B0Iz = −γI B0=Iz ,

(2)

 

5

 

 

где проекция Iz на ось Z спина ядра I принимает при квантовомеханическом рассмотрении [6] значения =Iz ; здесь = – постоянная Планка, а Iz – безразмерная величина, которая меняется через единицу на интервале I Iz I , где положительно определенная величина I , называмая спином ядра и имеющая целые и полуцелые значения, определяет наблюдаемые значения вектора во второй степени спина ядра I 2 = =2I (I +1) . Если до включения магнитного поля B0 все состояния ядер со спином I и его различными проекциями Iz были вырождены (т.е. их энергии не зависели от Iz ), то после включения этого поля, как видно из формулы (2), происхо-

дит снятие вырождения и энергетический спектр ядра принимает вид, представленный на рис. 1 для gI >0.

Рис. 1

При этом расстояние E между соседними магнитными подуровнями,

отличающимися по Iz на ±1, определяется формулой

 

E = =ω0 = = γI B0 ,

(3)

где величина ω0 совпадает с модулем ларморовской частоты, которая вво-

дится при описании [7] классической прецессии магнитного дипольного

момента ядра μG вокруг вектора магнитной индукции B .

 

 

 

0

 

Второй эффект связан с появлением действующей на ядро во внешнем

магнитном поле силы FG , проекции которой на оси α (α = Х, Y. Z) системы

координат имеют вид

 

dB0 .

 

F

= − E = μ

(4)

α

α

z dx

 

 

 

α

 

Из формулы (4) следует, что сила F отличается от нуля только в случае неоднородногоG магнитного поля, когда величина его магнитной индук-

ции B0 зависит от координат Х, У и Z.

Наконец, третий эффект связан с появлением во внешнем магнитном

поле действующего на ядро момента сил P :

 

 

G

G

 

(5)

μ

 

P =

, B0 .

6

 

 

 

 

Тогда изменение полного момента количества движения ядраG, совпадающего с его полным спином I , под действием момента сил P может быть описано в рамках классической механики с помощью обобщенного уравнения Ньютона [7]:

dI

G

 

G

G

 

(6)

 

 

μ

 

dt

= P =

, B0 .

 

 

 

 

 

 

2.2. Ларморовская прецессия

Умножая уравнение (6) на γI , можно получить уравнение для магнит-

ного дипольного момента μG :

 

 

G

 

 

 

dμ

 

G

 

 

 

 

μ

 

(7)

 

 

 

 

dt

=γI

, B0 .

Если вектор μ(t ) представить как μ(t ) = eGxμx (t ) + eGy μy (t ) + eGz μz (t ), где

eGx , eGy , eGz – единичные орты осей координат Х,

У, Z соответственно, то

уравнение (7) можно представить в виде системы трех уравнений для проекций μx (t ), μy (t ) , μz (t ) :

dμx =γ

I

μ

B ,

 

 

dt

 

 

y 0

 

 

 

 

 

 

 

 

dμy

= −γI μx B0 ,

(8)

 

 

dt

 

 

 

 

 

 

 

dμ

z = 0.

 

 

 

dt

 

 

 

 

 

 

 

 

Умножая второе уравнение в формуле (8) на мнимую единицу i и складывая его с первым уравнением, можно получить уравнение для величины

A(t ) = μx (t ) +iμy (t ):

 

 

 

 

dA(t )

=iω0 A(t ) ,

 

 

 

 

 

 

(9)

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

ω0 = γI B0 .

 

 

 

 

 

 

(10)

 

 

 

 

 

 

 

 

 

 

При использовании начального условия μx (0) = μ0 ,

μy (0) = 0 решение

уравнения (9) представляется как

 

 

 

 

 

 

 

 

 

A(t ) = μ

 

eiωt ; μ

x

(t) = μ

 

cosω t ,

μ

y

(t ) = μ

 

sinω

t ,

(11)

 

 

 

 

0

 

 

0

 

 

где μ – начальное значение поперечной по отношению к оси Z компоненты вектора дипольного момента μ . В то же время решение третьего уравнения в формуле (8) имеет вид

7

μz (t) = μ& = const ,

(12)

где μ& – начальное значение продольной по отношению к оси Z компонен-

ты вектора μG .

Тогда временная зависимость вектора магнитного дипольного момента μG(t ), определяемая формулами (11) и (12), соответствует ларморовской

прецессии вектора μG вокруг вектора B0 с угловой скоростью ω0 , называемой ларморовской частотой. Как было показано выше, модуль этой частоты ω0 фигурирует в формуле (4) для расстояний между соседними магнитны-

ми подуровнями ядра во внешнем магнитном поле B0 .

 

 

 

 

 

 

 

 

 

Как видно из табл. 2, ларморовские частоты

 

ω0

 

при изменении векто-

 

 

 

 

 

 

ра магнитной индукции B0

внешнего однородного и постоянного во време-

ни магнитного поля от 0,15 Тл до 2 Тл для ядра 1H1 и от 1 Тл до 3 Тл для

ядра

 

 

31P

 

лежат в радиотехнических диапазонах 6,4

 

 

ω

 

85,2 МГц и 17,2

 

 

 

 

 

 

 

15

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

ω0

 

 

51,7 МГц соответственно.

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ядра

 

 

 

 

1H1

1H1

1H1

 

1H1

1H1

1H1

 

 

31P15

 

31P15

31P15

31P15

 

 

Внешнее

0,15

0,35

0,5

 

1

1,5

2

 

 

 

1

 

1,5

2

3

 

 

поле B

(Тл)

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Частота

 

ω0

 

 

6,4

14,9

21,3

 

42,6

63,9

85,2

 

17,2

 

25,9

34,5

51,7

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(МГц)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.3. Ядерный парамагнетизм

Рассмотрим равновесное состояние системы тождественных атомных ядер со спиномG I ≠ 0 в однородном и постоянном во времени внешнем маг-

нитном поле B0 при температуре T0 . Функция распределения W (Iz ) проек-

ций Iz спинов ядер для этого состояния носит характер распределения Гиб-

бса, зависящего от энергии E (2) ядра в поле B [5]:

 

 

 

 

0

 

 

 

W (Iz ) = Cexp

E

 

= Cexp(a Iz ),

(13)

kT0

 

 

 

 

где a = =γI B0 . Нормировочная константа С находится из условия норми- kT0

ровки распределения (13) и имеет значение

8

 

 

I

 

1

 

C =

exp(aIz )

.

(14)

 

z

=−I

 

 

 

I

 

 

 

 

Поскольку распределение (13) зависит только от значений проекции Iz спинаGядра на ось Z, направление которой совпадает с направлениями вектора B0 , то средние значения μx и μy проекций дипольного магнитного момента ядра μx и μy на оси X и Y будут равны нулю, а среднее значение μz проекции μz дипольного момента ядра на ось Z при использовании формулы (1) определится как

 

 

 

I

 

 

 

 

 

γI Iz exp(aIz )

 

 

μ

z

=

Iz =−I

.

(15)

I

 

 

 

 

 

 

 

exp(aIz )

 

 

Iz =−I

Рассмотрим случай, когда температура системы ядер T0 достаточно

велика, так что величина aIz = =ω0Iz << 1. Тогда можно разложить функ- kT0

цию exp(aIz )в ряд по aIz и ограничиться линейным приближением. В этом случае формула (15) может быть представлена как

 

 

=γI Iz (1+ aIz )

 

 

μz

=

Iz

 

.

(16)

 

(1+ aIz )

 

 

 

 

 

Поскольку суммы по Iz

 

 

Iz

 

 

в формуле (16) содержат одинаковое число

 

 

 

I

 

 

членов с Iz > 0 и Iz < 0 , то сумма

Iz = 0 и формула (16) приводится к

Iz =−I

виду

 

=aγI

I

μz =

Iz2 .

 

 

(2I +1)

 

 

Iz =−I

Учитывая хорошо известную в теории простых чисел формулу

 

 

 

I

(2I +1)I (I +1)

 

 

 

 

 

Iz2 =

,

 

 

 

Iz =−I

 

 

3

 

 

 

 

 

 

 

 

 

величину μz (17) можно представить в форме

 

 

μz

=

 

aγI I (I +1)

=

=2γI2I (I +1)

B0 .

3

 

3kT0

 

 

 

 

 

 

 

 

 

9

 

 

 

 

 

(17)

(18)

(19)

Если ввести плотность ядер n0 , совпадающую с числом ядер в единице

объема, то можно определить вектор намагничивания M системы ядер, определяемый суммарным магнитным дипольным моментом всех ядер, нахо-

дящихся в единице объема [5]:

 

 

eG

= χBG

 

 

M = n

μ

z

,

(20)

0

 

z

0

 

 

где коэффициент намагничивания χ оказывается положительным и имеет значение

 

n

=2γ 2I (I +1)

 

 

χ =

0

I

.

(21)

 

3kT0

 

 

 

 

Формула (21) соответствует формуле, полученной для парамагнетиков Пьером Кюри еще в конце XIX века. Таким образом, при не слишком низких температурах магнетизм системы атомных ядер со спинами, не равными нулю, имеет парамагнитный характер.

Следует подчеркнуть: как среднее значение проекции магнитногоG дипольного момента ядра μz (15) во внешнем магнитном поле B0 при тем-

пературе T0 , так и связанная с этой проекцией величина вектора намагни-

ченности системы ядер MG (20) могут меняться непрерывным образом и имеют классический характер. Это позволяет использовать для указанных величин классические уравнения движения.

2.4.Спин-спиновое взаимодействие ядер

иядерный ферромагнетизм

Заметим, что система атомных ядер может, в принципе, обладать и

ферромагнитными свойствами. Действительно, атомное ядро с отличным от нуля магнитным дипольным моментом μ на достаточно большихGрасстоя-

ниях R от себя создает магнитное поле с векторным потенциалом A1 [5, 8]

 

 

G

 

 

μG

, R

 

 

 

 

 

 

 

A

=

 

1

 

 

 

 

 

 

(22)

 

 

 

 

 

 

 

 

 

и напряженностью HG

 

1

 

 

R3

 

 

 

 

 

1 этого поля

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H1 = rotA1.

 

 

 

 

(23)

Это поле будет действовать на магнитный дипольный момент μG2 дру-

гого ядра и приводить к появлению энергии взаимодействия E12

двух рас-

сматриваемых ядер [5, 8]

 

 

 

 

 

G

G

 

 

 

G

G

 

 

G

 

 

 

 

 

 

 

μ1, R

 

 

 

 

E12 = −μ2H1

= −μ2rot

 

 

 

,

(24)

 

R3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10

 

 

 

 

 

 

 

 

 

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