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

книги из ГПНТБ / Сысоев, А. Н. Гидродинамика сжимаемой жидкости учеб. пособие

.pdf
Скачиваний:
25
Добавлен:
19.10.2023
Размер:
5.05 Mб
Скачать

 

 

Основные понятия ЯМ Р

37

 

+/

 

 

Рис. 1.6. Аппроксимация

у-(і) о

 

 

прямоугольной функции

'

 

 

рядом Фурье по синусам

 

 

 

с учетом 1,3 и 5 членов.

 

 

 

 

- /

 

 

 

о

т

мощью так называемого гармонического анализа, или раз­ ложения в ряд Фурье. Это математический метод разложе­ ния сложного колебания на его спектральные компоненты. Сложное колебание как функцию времени часто называют функцией во временной области, тогда как соответствую­ щий спектр называют функцией в частотной области. С помощью методов гармонического спектрального анализа можно преобразовывать данные из одной области в другую. В гл. 4 мы применим преобразование Фурье, чтобы опреде­ лить частоты молекулярных вращений, характерные для типичных хаотических молекулярных движений. В гл. 5 преобразование Фурье используется для определения ха­ рактеристических частот и интенсивностей спектра ЯМР по временной реакции (интерферограмме) системы ядерных спинов на ВЧ-импульс.

Как известно, функцию f(t) обычно можно представить

в виде ряда Фурье, т. е. бесконечного

ряда синусов и коси­

нусов:

 

 

СО

со

 

f(t)= 2

Д ,c° s (n~/T)t + 2

Basin(n*/T)t. (1.45)

п = 0

п=1

 

В учебниках математики [11] показывается, что это разло­

жение в ряд действительно в области —Т

t

Т\ там же

выводятся

формулы для определения коэффициентов А Л

и Вп. Если

f(t) — симметричная, или

четная,

функция

[т. е. /(—t) = f(t)], то все Вп равны нулю, и в формуле ос­ таются только косинусы. Если fff) антисимметричная, или нечетная, функция [т. е. /(—t)—/(/)], то остается толь­ ко ряд по синусам. Асимметричная функция общего вида

38 Глава 1

определяется полным выражением (1.45). На рис. 1.6 по­ казано представление простой прямоугольной функции рядом Фурье по синусам; видно, что включение в ряд все большего числа членов дает все лучшее приближение к заданной функции.

Возможность выражать функцию в виде ряда Фурье из синусов и косинусов с частотами 1/2T, 2I2T, 3/2Т и т. д. приводит к мысли описывать функцию просто набором коэффициентов при последовательных членах ряда. Тогда о функции, в разложении которой коэффициенты A it А 2, А3 велики, а А 30, А31, ... малы, говорят, что она содержит большой вклад от низкочастотных компонент Фурье и ма­ лый вклад от высокочастотных. Внимательное изучение рис. 1.6 показывает, что высокочастотные компоненты обес­ печивают быстрое нарастание, необходимое для аппрокси­ мации функции типа прямоугольной, так что для функции такого типа требуется широкий диапазон компонент Фурье. Поскольку истинные частоты компонент пропорциональны 1/27", уменьшение длительности (или периода) прямоуголь­ ных колебаний также ведет к увеличению вклада высокочастотных компонент.

Пример прямоугольного колебания имеет непосред­ ственное отношение к импульсным методам, так как импульс . генерируется путем быстрого включения ВЧ-генератора, так что огибающая мощного ВЧ-сигнала, поступающего на образец, очень близка к форме положительной полувол­ ны на рис. 1.6. Если генератор работает непрерывно, как в обычном ЯМР-эксперименте, то на образец воздействует только основная частота ѵ0. Благодаря быстрому включению генерируются компоненты Фурье, которые добавляются (или вычитаются) к основной частоте, и на образец дейст­ вует целая полоса частот, приблизительно лежащая в диа­ пазоне ѵ0± 1//р, где /р— длительность импульса.

Когда мы фактически проводим математические опера­ ции спектрального анализа, то нам удобнее иметь дело не *- с рядами Фурье, а с соответствующими интегралами, в ко­ торых снято ограничение на область определения входящих

в них

функций:

переменная t не

ограничена областью ^

— Т

Т, а

может изменяться

в бесконечных преде­

лах. В этом случае удобно определить преобразование Фурье от /(/) следующим образом:

Основные понятия ЯМ Р 39

F(u>) = J f(t)exp(— iui)dt.

(1-46)

—СО

 

^'Здесь со — угловая частота (со = 2яѵ), і — У —1. Можно показать [11, 12], что верно и обратное соотношение

/

°°

(1-47)

f(t) = (1/2іг)

J f((o) exp (г со/) dco.

со

Вэкспоненциальную функцию от мнимого аргумента в формулах (1.46) и (1.47) входят и синус- и косинус-компо­ ненты. Можно ввести также раздельные синус-преобразо­ вание Фурье

<

S(u>) = 2

J

f(i) sin (о idt

(1.48)

 

 

о

 

 

и косинус-преобразование

 

Фурье

 

г

C(cu)=2

j

f(t) cos uitdt.

(1.49)

 

Ч

 

о

 

 

 

 

 

 

Г л а в а 2

Свободная индукция и спиновое эхо

В этой главе мы рассмотрим поведение спиновой систе­ мы, подвергнутой воздействию одного или нескольких ВЧ-импульсов. Мы будем считать, что ВЧ от передатчика включается и выключается так быстро, что огибающая ВЧколебаний во время импульса имеет прямоугольную форму. Примем также, что длительность импульса /р мала по срав­ нению с Ті и Т2, так что релаксация за время /р не проис­ ходит. (Как мы увидим в гл. 3, часто используются импульсы длительностью всего несколько микросекунд, так что эти требования удовлетворяются без труда.) С помощью фор­ мулы (1.39) можно найти угол, на который повернется М под действием заданного импульса, и мы будем обычно ха­ рактеризовать импульсы как 90°-ный (тс/2), 180°-ный (тс) и т. д. Здесь мы будем рассматривать только ВЧ-сигнал, возникающий поем выключения импульса. Поведение сис­ темы во время импульса мы исследуем в гл. 6.

2.1. • Спад индуцированного сигнала (СИС)

Пусть вдоль оси х' в системе координат, вращающейся с частотой ВЧ-поля, приложен 90°-ный импульс. По оконча­ нии импульса намагниченность М направлена точно по оси у', как показано на рис. 2.1, а. Построение спектромет­ ра обычно таково, что в нем регистрируется сигнал, наве­ денный (индуцированный) в катушке, расположенной вдоль неподвижной оси х или у (см. гл. 3); поэтому интенсивность наблюдаемого сигнала определяется величиной Мху. Этот наведенный сигнал называют-сигналом свободной индукции, поскольку ядра прецессируют «свободно» в отсутствие при­ ложенного ВЧ-поля. По мере развития поперечной релак-

Свободная индукция и спиновое эхо 41

сации сигнал уменьшается (спадает). В идеально однород­ ном поле постоянная времени спада была бы равна Т2; однако фактически -сигнал свободной индукции спадает с характеристическим временем Та*, которое часто опре­

деляется

прежде

всего неоднородностью поля, посколь­

ку ядра

в разных

частях поля прецессируют с несколько

разными

частотами

и,

следовательно, быстро расходят­

ся по фазе. Поэтому

сигнал спадает с характеристичес­

ким временем Т2*,

определенным с помощью соотноше­

ния (1.11).

На рис. 2.1, г показан чисто экспоненциальный спад сигнала, который получается, если частота ВЧ-импульса в точности равна резонансной частоте единственного типа ядер в образце. Этот спад дает непосредственную меру уменьшения Мху, поскольку обычно в спектрометре ис­ пользуется детектор, в котором фаза опорного сигнала свя­ зана с фазой ВЧ-сигнала (гл. 3). Таким образом, хотя реги­ страция сигнала производится в отсутствие непосредствен­ ного воздействия ВЧ-импульса (т. е. после импульса) на

Рис. 2.1.

а — 90°-ный импульс вдоль оси х ’ поворачивает М от равновесного положения до совпадения с осью у ' ; б — спад М при расфаэировании магнитных моментов; в— входной сигнал — 90°-ный импульс, соответствующий а; г — экспоненциальный спад

индуцированного сигнала, соответствующий б .

42 Глава 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

систему

ядерных

спинов,

 

 

 

 

опорный

ВЧ-сигнал

пода­

 

 

 

 

ется на детектор непрерыв­

I

 

 

 

но.

 

Поэтому

детектор

 

 

 

реагирует на ту часть на­

 

 

 

магниченности,

 

которая

 

 

 

 

находится в

определенном

 

 

 

 

фазовом соотношении с Hj

 

 

 

 

и,

следовательно,

 

лежит

 

 

 

 

вдоль неподвижной

оси во

 

 

 

 

вращающейся

системе

 

 

 

 

данном

случае

вдоль

по­

 

 

 

 

ложительного направления

 

 

 

 

оси

у').

 

 

 

теперь,

 

 

 

 

 

Предположим

 

 

 

 

 

что частота

ВЧ-импульса

 

 

 

 

слегка

отличается

от лар­

 

 

 

 

моровой

частоты

 

ядер.

 

 

 

10

Если

опять

рассмотреть

 

Время,(лс

систему координат, враща­

 

 

 

 

ющуюся с частотой ВЧ-им­

 

Рис. 2.2.

 

 

 

пульса, то сразу после 90°-

а — спад индуцированного сигнала при вы­

ного импульса М будет

ле­

сокой частоте, точно

равной ларморовой;

б — спад

индуцированного

сигнала при

жать вдоль

оси

у'1. Од­

расстройке

ВЧ относительно резонанса.

 

 

 

 

нако М теперь

вращает­

ющейся

 

 

 

ся

относительно

 

враща­

системы, и детектор выделяет

не только

 

экспо­

ненциальный

спад

величины

Мху, но и эффект интер­

ференции Мху с опорным сигналом вследствие того, что они попеременно то совпадают по фазе, то расходятся. Типич­ ный график для этого случая показан на рис. 2.2; он ана­ логичен биениям («виглям»), наблюдаемым в обычном спек­ тре ЯМР высокого разрешения.

Приведенные выше рассуждения можно распространить на случай системы, содержащей несколько сортов ядер од­ ного типа (например, протонов), различающихся по лармо-1

1 Здесь, вероятно, возникнет вопрос, каким образом вектор намагниченности может быть повернут до плоскости ху полем Нъ Имеющим частоту, не равную резонансной. Мы ответим на этот вопрос в гл. 5.

Свободная индукция и спиновое эхо

43

ровой частоте вследствие химического сдвига и (или) спинспинового взаимодействия. Очевидно, теперь некоторые яд­ ра будут прецессировать с частотами, отличающимися от частоты вращения системы координат, и могут наблюдаться эффекты интерференции, похожие на показанный на рис. 2.2, но часто значительно более сложные. Мы вернем­ ся к подробному рассмотрению этого случая в разд. 2.6

ив гл. 5, а в разд. 2.2 — 2.5 рассмотрим спин-систему с един­ ственной ларморовой частотой.

Регистрация спада индуцированного сигнала — основ­ ной способ определения величины и других характеристик

М.Спад индуцированного сигнала, следующий за 90°-ным импульсом, несет спектральную-информацию, выявляемую в фурье-спектроскопии ЯМР, а спад индуцированного сиг­ нала, наблюдаемый вслед за последовательностями из двух

иболее импульсов, применяется для определения времен

релаксации.

2.2. Измерение Тх

Импульсные методы дают наиболее универсальный спо­ соб измерения времен Ті в широком диапазоне значений. Наиболее широко применяемый способ — так называемая импульсная последовательность 180°, т, 90° (где т — про­ межуток времени между соседними импульсами) — пока­ зан на рис. 2.3. Сначала 180°-ный импульс инвертирует намагниченность вдоль оси г'. Далее происходит продоль­ ная релаксация, под действием которой Мг■изменяется от значения —М 0, проходит через нуль и стремится к своему равновесному значению УИ0Если через время % после 180°-ного импульса к системе приложить 90°-ный импульс, также направленный по оси х ', то вектор намагниченности М повернется и окажется направленным по оси у'. В ре­ зультате будет наблюдаться сигнал свободной индукции, начальная амплитуда которого пропорциональна величи­ не М и, следовательно, величине Мг<в момент времени т. Если теперь позволить системе вернуться к равновесию, для чего необходимо выждать по крайней мере время 57Y, и снова воздействовать на нее последовательностью1

1 При 5Ті Мг' = 0,993 М„.

44 Глава 2

импульсов 180°, т, 90°, но с другим значением т, то можно найти скорость спада M z, как показано на рис. 2.3. Диа­

пазон изменения % зависит, конечно, от

и, вообще гово­

ря, должен доходить примерно до 3—4 7V

 

 

Количественно спад M z описывается уравнениями Бло­

ха (1.28) при Мж= М у= 0:

 

 

 

dMJdt = — (Mz — М0)/Т{.

 

(2.1)

Интегрирование

уравнения (2.1) с начальным

условием

М 2= М 0 при

t = 0 дает

 

 

MZ = M0[ \ — 2ехр(— t/Ti)].

 

(2.2)

Это и есть функция, построенная на рис. 2.3. На прак­ тике выражение (2.2) используется в преобразованном виде

г'

в

о

'■ н у л ь

Рис. 2.3. Измерение Ті с помощью-импульсной последовательности

 

1 8 0 °, т ,

9 0 ° .

а — в момент < = 0

М инвертируется 180°-ным импульсом: б — через промежуток

времени т 90°-ный

импульс поворачивает

М до совпадения с осью у ' (или — у ’ ) ;

в — строится зависимость от т начального значения индуцированного сигнала после

90°-ного импульса, которое пропорционально М в момент времени т. Для построе­ ния каждой точки нужна своя последовательность 180°, г, 90°. Стрелкой отмечена точка, соответствующая случаю, показанному на схеме б .

Свободная индукция и спиновое эхо

45

\п(Аа - А х) = \п2Аю- х / Т и

(2.3)

где А х — начальная амплитуда индуцированного

сигнала

после 90°-ного импульса, приложенного в момент t, а А т— предельное значение Лт при очень длительном интервале между 180°-ным и 90°-ным импульсами. Тогда Tt опреде­

ляется

по

наклону графика зависимости

In (Лга — Лт)

от

т:.

 

 

 

 

=

Из

выражения

(2.3) можно видеть,

что при т 0 =

Тііп 2 =

0,697^

Аг— 0. Поэтому Тх можно определить

по временному промежутку х между импульсами, при ко­ тором 90°-ный импульс не сопровождается сигналом сво­ бодной индукции. Этот способ полезен для грубой оценки Ти но не годится для точных измерений. Причину этого мы увидим, когда в гл. 3 рассмотрим аппаратуру импуль­

сного метода

ЯМР.

Если Г, »

Т2, то можно определить также с помощью

1импульсной последовательности 90°, х , 90°. В этом случае сигнал свободной индукции после первого 90°-ного импуль­ са спадает до нуля значительно быстрее, чем M z достигает своего равновесного значения. Поэтому второй 90°-ный импульс позволяет определить значение M z в момент х, как и с помощью последовательности 180°, х, 90°. Если ряд 90°-ных импульсов повторяется с интервалами х , то М бы­ стро достигает стационарного состояния, величина М в котором зависит отт и Tt. После фиксации стационарной величины сигнала последовательность 90°, т, 90°, х, 90°, т , ... повторяется с другим значением т . Поскольку М по­ ворачивается только на 90°, а не на 180°, то начальное ус­ ловие для интегрирования уравнения (2.1) теперь имеет

вид М г— 0

при

t — 0.

Поэтому график зависимости

In(До,— Лт)

отт

снова

является прямой линией, по нак­

лону которой можно найти 7V Последовательность 90°, т, 90°,т , 90°, ... при больших 7\ часто позволяет сэкономить время, особенно если требуется усреднение по времени (на­ копление) повторяющихся сигналов для повышения отно­ шения сигнала к шуму. Описанный способ пригоден лишь при Т2<С Тс, в противном случае оставшаяся-компонента намагниченности, направленная вдоль у", в момент време­ ни т будет повернута в направлении оси г' и соотношение, используемое для определения Ти окажется неприменимым.-

46 Глава 2

Рис.

2.4.

Эксперимент

Хана

со

спин-эхо.

 

 

 

а — 90°-ный импульс,

приложенный в момент времени t

=

0 вдоль осн х ' ,

застав­

ляет Лі повернуться

до

совпадения

с положительным

направлением осн у ' ;

б —

макроскопические намагниченности

т . ядер,

находящихся в разчых частях образда,

расфазируются под влиянием неоднородности

в

напряженности Н0. Ядра,

прецес­

сирующие быстрее, чем со средней

скоростью вращающейся системы

координат,

во вращающейся системе движутся к наблюдателю.

Если смотреть вниз со сторо­

ны положительной осн г ',

то они движутся по часовой

стрелке, тогда

как ядра,

движущиеся медленнее,

чем вращающаяся

система, движутся

против

часовой

стрелки; в — 180°-ный импульс,

приложенный в момент

времени т вдоль осн де',

заставляет все т . повернуться на

180° вокруг осн х ' \

г — более

быстрые

ядра,

продолжая двигаться во вращающейся

системе

по

часовой стрелке,

теперь уда­

ляются от наблюдателя,

тогда

как медленные

ядра движутся

против

часовой

стрелки к наблюдателю; 3 — в момент времени 2 т все ш . оказываются

совпадаю­

щими по фазе и направленными вдоль

осн — у ' ;

е — при

і

> 2 т т .

снова расфазн-

 

 

 

 

руются.

 

 

 

 

 

 

 

 

 

В разд. 5.6 мы покажем,что существуют методики, в которых указанное ограничение заменяется на менее строгое, а именно Т2* < 7\. Для наиболее точных измерений, вообще говоря, используется методика 180°, т, 90°.

2.3.Измерение Г2-методом спинового эхо

Вслучаях, когда не выполняется условие Т2 < (2/7 Д#„) [см. выражение (1.11)1, вклад неоднородности поля Н0 в скорость спада индуцированного сигнала не позволяет ис­ пользовать это время спада Т2* как меру Т2. Остроумный способ преодоления влияния неоднородности впервые был предложен Ханом [2], который назвал его методом спино­ вого эхо. Метод состоит в том, что на спиновую систему воз­

Соседние файлы в папке книги из ГПНТБ