Ядерный магнитный резонанс и ЯМР-томография (110
..pdfМИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РФ ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ «ВОРОНЕЖСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ»
ЯДЕРНЫЙ МАГНИТНЫЙ РЕЗОНАНС И ЯМР-ТОМОГРАФИЯ
Учебное пособие для вузов
Составитель С. Г. Кадменский
Издательско-полиграфический центр Воронежского государственного университета
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 |
|
|
|
|
|
|
|
|
|