- •1.2 Акустические волны в скалярных средах
- •1.3 Интенсивность ультразвукового излучения
- •1.4 Биологическое действие ультразвука
- •1.5 Сферические волны
- •Функция Грина и формулы Кирхгофа
- •2.4 Поле сферического фокусирующего преобразователя
- •3.2 Уравнение Вестервельта и волны комбинационных частот
- •3.2. Коэффициент нелинейности биологических сред
- •Таблица 3.2
- •4. ГИДРОДИНАМИЧЕСКИЕ МЕТОДЫ
- •4.1 Уравнения течения вязкой несжимаемой жидкости
- •4.2 Элементы гемодинамики
- •4.3 Аппарат искусственного кровообращения
- •4.4 Вискозиметрия
- •4.5 Измерение артериального давления крови
- •5.1 Методы ультрацентрифугирования
- •5.2 Метод скорости седиментации
- •5.1. Электрические свойства биотканей
- •5.2. Источник внеклеточного поля
- •5.3 Воздействие постоянными и переменными токами и полями
- •5.4 Механизмы поглощения высокочастотных полей
- •Спектр ЭМИ
- •Вопросы для самостоятельного изучения
- •Литература
- •Излучение
- •Ионизирующее
- •6. ПОГЛОЩЕНИЕ ЭМИ
- •6.1 Поглощение СВЧ волн в биологических тканях
- •6.2 Обратные переходы и процессы релаксации
- •6.3 Флуктуационно-диссипационная теорема
- •6.4 Индуктотермия
- •6.5 Биологическое действие высокочастотного ЭМИ
- •Вопросы для самостоятельного изучения
- •Литература
- •7.1 История развития метода
- •7.2 Условие резонанса
- •7.3 Кинетика переходов
- •7.4 Уравнение Блоха
- •7.5 Спектр ЯМР
- •Квадратурное детектирование ЯМР сигналов
- •Чувствительность cw-ЯМР
- •РЧ импульсы и импульсная ЯМР-спектроскопия
- •Чувствительность импульсного ЯМР
- •Применение ЯМР-спектроскопии в биохимии
- •Зависящие от времени процессы в ЯМР
- •Химический анализ
- •Тип движения
- •Область значений частоты, Гц
- •Вибрационные и торсионные движения
- •Латеральная диффузия в мембранах
- •Диффузия в растворах
- •Вращательная диффузия в растворах
- •Конформационные изменения протеинов
НЕПРЕРЫВНЫЙ И ИМПУЛЬСНЫЙ ЯМР
Квадратурное детектирование ЯМР сигналов
До сих пор мы имели дело с откликом на линейно поляризованное внешнее РЧ поле, излучаемое непрерывно. Метод непрерывно излучаемого РЧ поля наиболее интенсивно использовался первые 20 лет применения ЯМРспектроскопии и называется cw-ЯМР. Как было показано, внешнее РЧ поле с линейной поляризации при cw-ЯМР демонстрирует все особенности магнитного возбуждения системы ядерных спинов и поглощения энергии поля. В экспериментах как с непрерывным, так и с импульсными РЧ полями используют также поле с круговой поляризацией
H1 = i H1 cosωt + jH1 sinωt . |
(8.1) |
С точки зрения результата резонансного возбуждения системы и физики протекающих при этом процессов, например, процессов релаксации, такая замена линейно поляризованного поля на поле (8.1) с круговой поляризацией в
плоскости, перпендикулярной направлению поля H0 , ничего не меняет. Однако
при рассмотрении таких полей некоторые экспериментальные особенности всех ЯМР-методов становятся более наглядными. Это относится, в частности, к
квадратурному детектированию сигналов ЯМР отклика.
При наличии дополнительного магнитного поля вида (8.1) алгоритм решения уравнения Блоха и вычисления магнитной восприимчивости в линейном приближении такой же, как и при линейной поляризации поля. Но в этом случае гармоническое решение можно представить в виде
M + (t) = M + (ω)eiωt |
≡{M x (ω) +iM y (ω)}eiωt = (χ+′ |
+iχ+′′)H1eiωt , |
(8.2) |
где введена комплексная |
вращательная компонента намагниченности |
||
M+ (t) = M x (t) +iM y (t). Зависимость (8.2) компоненты |
Фурье или |
спектра |
комплексного сигнала вращательной намагниченности от поля описывается очень простой магнитной восприимчивостью:
|
ω −ω0 |
|
|
|
1 T2 |
|
|
|
|||
χ+ (ω) = −γM0 |
|
|
|
|
|
+i |
|
|
|
. |
(8.3) |
(ω −ω0 ) |
2 |
+ |
1 |
2 |
(ω0 −ω) |
2 |
2 |
||||
|
|
T2 |
|
+1 T2 |
|
|
Как видим, действительная и мнимая части восприимчивости (8.3) описывают по отдельности дисперсию системы ядерных спинов и поглощение ею энергии внешнего РЧ поля. При заданной частоте поля РЧ это выражение устанавливает связь между амплитудой колебаний намагниченности и величиной поля.
Если вектор намагниченности прецессирует вокруг направления поля H0 ,
то изменение во времени поперечных компонент намагниченности и сигналов, которые они дают, всегда сдвинуто по фазе на π2, как показано на рис.1б для
двух условных сигналов. По физическому смыслу именно это обстоятельство делает возможным представление этих компонент в виде действительной и мнимой частей некоторой комплексной вращательной намагниченности,
Мнимая часть
Действительная часть
Рис.1. а) зависимость от времени регистрируемых поперечных компонент вектора намагниченности; б) спектры сигналов поглощения и дисперсии
поскольку в комплексной плоскости ее действительная и мнимая части также сдвинуты на π2 по фазе. Отсюда вытекают два весьма полезных для
практических приложений вывода.
Первый сводится к тому, что возможно применение квадратурного детектирования, при котором регистрируемый РЧ сигнал, обычно в диапазоне 20-100МГц, преобразуется с помощью фазовых гетеродинов в два сигнала звуковой частоты (2-10кГц). Опорные частоты этих гетеродинов сдвинуты на π2, поэтому и получающиеся сигналы имеют такую же фазовую сдвижку
(аналогично детектированию в ультразвуковых доплеровских методах). Спектр комплексной величины, составленной из таких произвольных сигналов, с точностью до постоянного фазового множителя всегда описывает магнитную восприимчивость (8.3) системы ядерных спинов. Если регистрируются сигналы от составляющих намагниченности M x (t) и M y (t) , то фазовый множитель равен
единице.
Второе преимущество квадратурного детектирования вытекает непосредственно из выражения (8.2). В соответствии с этим выражением
сигналу, у которого компонента Фурье действительна |
и равна |
||
Re M x (ω) + Im M y (ω) |
соответствует только |
действительная часть |
магнитной |
восприимчивости, а |
сигналу со спектром |
Im M x (ω) + Re M y (ω) – |
только ее |
мнимая часть, как показано на рис.1б. Поэтому первый из этих двух сигналов называется сигналом дисперсии, а второй – сигналом поглощения. Если в процессе измерений формируется только один из этих сигналов, то говорят либо
о режиме поглощения, либо о режиме дисперсии.
Чувствительность cw-ЯМР
Для измерения спектров ЯМР высокого разрешения при непрерывном режиме излучения образец помещается в однородное магнитное поле и подвергается непрерывному воздействию РЧ поля. Это поле варьируется по частоте в области существования линий поглощения ЯМР для данного образца. В ЯМР высокого разрешения регистрации подлежит обычно линия поглощения.
Чтобы детальнее разобраться с сигналами поглощения и насыщения, чувствительностью cw-ЯМР и резонансным влиянием РЧ поля, снова воспользуемся уравнением Блоха, которое запишем в виде
dM |
|
|
|
M |
|
M0 |
− M|| |
|
|
dt |
=γ[M ×(H0 |
+ H1 )] − |
T2 |
+ |
|
|
. |
(4) |
|
|
T1 |
||||||||
|
|
|
|
|
|
|
где M|| = kM z , а магнитное поле обладает круговой поляризацией (1). С учетом этого обстоятельства удобно перейти к новой системе координат, которая вращается вокруг направления поля H0 в том же направлении с частотой ω ,
близкой к ларморовской частоте прецессии вектора намагниченности. Для этого достаточно учесть, что скорость изменения вектора намагниченности в лабораторной системе координат и во вращающейся связаны соотношением
ddtM = ddtM ′ +[ω × M ′].
Подставляя это выражение в (4) и производя несложные преобразования, приходим к уравнению Блоха вида
dM ′ |
|
ω |
|
ω |
|
M |
′ |
M0 |
− M|| |
|
|
|
=γ[M ′×(− |
0 |
+ |
|
+ H1′)] − |
+ |
|
|
, |
(5) |
|
dt |
γ |
|
T1 |
||||||||
|
γ |
|
|
T2 |
|
|
|
|
где ω0 = −γH0 - вектор ларморовской угловой скорости.
Из уравнения (5) следует, что во вращающейся системе координат магнитный момент должен прецессировать вокруг направления некоторого
эффективного поля He , которое является векторной суммой полей γ −1 (−ω0 +ω)
и неподвижного во вращающейся системе отсчета поля H1′. Особенно простым является вид уравнения Блоха в случае полного совпадения частот, когда эффективное поле совпадает с полем H1′. Если считать для определенности это
поле направленным вдоль вращающейся оси x′, то H1′ = i′H1 , а угловая скорость прецессии вокруг него равна
ω′ = −γ −1H1′. |
(6) |
Знак минус в выражении (6) указывает на направление вращения, которое показано на рис.2а вместе с направлением поля.
В явном виде связь между компонентами намагниченности во вращающейся системе отсчета и в лабораторной дается выражениями:
M x′ = M x cosωt − M y sinωt , |
(7) |
M ′y = M x sinωt + M y cosωt . |
(8) |
С течением времени система ядерных спинов релаксирует к некоторому постоянному равновесному значению вектора намагниченности, который
направлен вдоль эффективного поля He точно так же, как он направлен вдоль поля H0 в отсутствие внешнего РЧ поля. После этого величина намагниченности во вращающейся системе отсчета не изменяется, поэтому в уравнении (5) можно положить dM ′dt = 0 . Получающаяся в результате система
алгебраических уравнений дает стационарное решение уравнений Блоха, которое имеет простой вид
M x′ |
|
|
|
|
|
γH (ω −ω)T 2 |
|
|
|
(9) |
||||||
= M0 1+ (ω −ω)2 T 2 |
+γ 2 H 2TT , |
|||||||||||||||
|
|
|
|
|
|
|
1 |
|
|
0 |
|
2 |
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
2 |
1 |
1 |
2 |
|
|
M ′y |
= M0 |
|
|
|
|
|
γH1T2 |
|
|
, |
(10) |
|||||
1 |
+ (ω |
0 |
−ω)2 T 2 |
+γ 2 H |
2TT |
|||||||||||
|
|
|
|
|
|
|
|
|
|
2 |
1 |
1 |
2 |
|
|
|
M x′ |
|
|
|
|
|
1+ (ω −ω)2 T 2 |
|
|
|
(11) |
||||||
= M0 1+ (ω −ω)2 T 2 |
+γ 2 H 2TT . |
|||||||||||||||
|
|
|
|
|
|
|
|
|
0 |
|
|
2 |
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
2 |
1 |
1 |
2 |
|
|
Теперь видно, что сигналом поглощения в данном случае является сигнал |
||||||||||||||||
M ′y , а сигналом дисперсии |
– |
|
M x′. Для очень слабых РЧ полей, когда |
|||||||||||||
выполняется неравенство |
γ |
2 H |
1 |
2TT <<1, |
|
|
сигнал |
поглощения |
описывается |
|||||||
|
|
|
1 |
2 |
|
|
|
|
|
|
|
|
|
|
|
лоренцевской кривой, как это было в линейном приближении при решении уравнений Блоха. В общем случае выражения (9)-(11) дают более правильное представление о временах релаксации и о влиянии на них постоянных T1 и T2 .
Кроме того, из этих выражений видно, что при очень большой величине внешнего РЧ поля амплитуда всех сигналов уменьшается. Это связано с тем, что при очень большой скорости поглощения энергии поля по сравнению с временами релаксации T1 и T2 заселенность уровней выравнивается, а их
разность стремится к нулю. Этот эффект, который приводит к уширению линий поглощения, называется насыщением.
Условие максимальности амплитуды сигнала поглощения можно, очевидно, получить из условия равенства нулю производной амплитуды по величине поля: ∂M ′y ∂H1 = 0 . Решение этого уравнения приводит к оптимальной
величине внешнего РЧ поля
H opt = (γ |
|
)−1 . |
(12) |
|
TT |
||||
1 |
1 |
2 |
|
|
Вследствие насыщения амплитуда наблюдаемого сигнала поглощения убывает при значениях поля выше и ниже оптимального значения, причем в последнем случае малая частота переходов между состояниями связана со слабостью поля и недостаточно интенсивным поглощением его энергии. Условие (12)
ограничивает область линейной зависимости амплитуды сигнала от поля H1 . С ростом поля возникающее уширение линий требует введения коэффициента s = 1+ (H1 H1opt )2 . Иными словами, полуширина линии равна ∆ = sT2 .