КВАНТОВОМЕХАНИЧЕСКИЙ ГАМИЛЬТОНИАН ЖЕСТКОГО ВОЛЧКА

Вращение молекул подчиняется законам квантовой механики. Вращательная энергия не может принимать произвольные величины, а только дискретный ряд значений. Очень большое количество молекул при своем вращении почти не изменяют своих размеров. Обычно это соответствует случаю, когда изменение длин связей между атомами за счет центробежных сил значительно меньше равновесных длин связей. Такие молекулы называются нормальными или квазижесткими и в первом приближении их можно считать абсолютно жесткими. В модели жесткой молекулы (жесткого волчка) предполагается, что взаимное расположение атомов друг относительно друга не меняется во время вращения молекулы как целого. Приближение жесткого волчка хорошо описывает качественную структуру вращательных уровней энергии квазижесткой молекулы, и во многих случаях дает близкие к реальным численные значения энергий. В качестве примера квазижестких молекул можно назвать такие молекулы, как SO2, HCN, O3, CO2, OCS, HCOOH, CH4, CH3Cl.

В данной лекции проводится простой вывод квантовомеханического гамильтониана жесткого волчка. Ход вывода следующий. Сначала записывается классический лагранжиан и с его помощью выводится гамильтониан жесткого волчка, то есть тела, движущегося по законам классической механики. При выводе в качестве обобщенных скоростей берутся проекции угловой скорости, через которые записывается кинетическая энергия. Далее определяются классические моменты количества движения как обобщенные импульсы, через которые выражается классический гамильтониан. После этого проводится квантование гамильтониана путем замены классических импульсов соответствующими квантовомеханическими операторами.

Функция Лагранжа жесткого волчка

Рассмотрим жесткую систему заданной структуры, состоящую из N материальных точек, которые соответствуют атомам молекулы с известными массами. Эта система представляет собой тело, часто называемое жестким волчком. Оно закреплено в центре масс и вокруг этой точки может вращаться. Классическая энергия складывается из кинетической и потенциальной энергий. Внутренняя потенциальная энергия частиц во время вращательного движения не меняется, так как она зависит только от взаимного расположения материальных точек, которые при вращении жесткого тела остаются неизменными. В силу этого, а так же в силу того, что за нулевой уровень энергии можно принять любой произвольный уровень, потенциальную энергию, как постоянную величину, можем положить равной 0, то есть U=0. Кинетическая энергия может быть записана, как

T=1/2 S(i) MiVi^2
(1)

Здесь Mi - массы частиц, Vi - их линейные скорости.

Для перехода к гамильтоновой форме энергии необходимо сначала выбрать обобщенные скорости, а затем определить сопряженные им обобщенные импульсы. Твердое тело с покоящимся центром масс имеет три степени свободы и столько же должно быть независимых обобщенных скоростей. Линейные скорости частиц не могут быть взяты в качестве обобщенных скоростей, так как в твердом теле их больше трех и они линейно зависимы. Можно было бы взять какие-либо три линейных скорости в качестве обобщенных координат, а остальные выразить через них, но при этом будет нарушен принцип равноправности. Для вращающейся молекулы в качестве обобщенных скоростей удобно взять декартовые компоненты угловой скорости omega . Их ровно три, они независимые, и через них могут быть выражены все Vi .

Vi = [omega X Ri]
(2)

где Ri радиус вектор, определяющий положение i-го атома.

Подставим выражение (2) в (1) в результате чего получим

T

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

T
(3)

Здесь введенные для сокращения записи величины Iab называются компонентами тензора инерции. Они зависят от масс частиц и их координат и определяются следующими формулами.

Для диагональных элементов

Iaa
(4)

Для недиагональных элементов

Iab
(5)

В этих формулах под alpha и beta понимаются обозначения осей X, Y, Z декартовой системы координат. Отметим, что тензор инерции симметричный, то есть

Iab=Iba
(6)

Выражение (3) часто записывается в более компактной матричной форме:

Н
(7)

где omega - вектор столбец, составленный из omegaX,omegaY,omegaz, I – матрица 3 x 3 тензора инерции, знак плюс над omega - операция транспонирования, так что omega+ это вектор строка. По сути, выражение (3) или (7) есть функция Лагранжа жесткого вращающегося волчка.

Гамильтониан жесткого волчка

Теперь перейдем к гамильтоновой форме. Сначала определим обобщенные импульсы Jalpha, сопряженные компонентам угловой скорости, используя известную формулу классической механики

Jalpha
(8)

или в матричной форме

Jalpha
(9)

Обращая (9) относительно omega и подставляя omega в (7) получаем выражение для энергии через импульсы J, что и является классическим гамильтонианом жесткого волчка.

omega=I^(-1)J

T=H=J(+)I(-1)J
(10)

Импульсы Jalpha на самом деле есть компоненты момента импульса, что легко проверить, если вспомнить формулу для этих величин и выразить их через угловые скорости:

(11)

Квантование гамильтониана жесткого волчка

Теперь можно перейти к квантовомеханическому гамильтониану, заменив классический момент импульса J на его квантовомеханический аналог, - оператор момента количества движения:

J=
(12)

Компоненты тензора инерции, входящие в выражение гамильтониана (10), имеют разные значения в разных координатных системах. Оказывается можно подобрать такую систему координат, что форма гамильтониана существенно упростится. В теории вращающейся молекулы используются два вида декартовых систем координат: лабораторная и молекулярная. Первый вид - это система координат, неподвижная в пространстве. В этой системе компоненты Iab будут меняться при повороте молекулы, так как будут меняться проекции координат вращающихся атомов. Второй вид системы координат, молекулярная система, прикреплена к волчку и вращается вместе с ним. В такой системе компоненты Iab есть постоянные величины, что является большим плюсом. Кроме того, молекулярную систему можно по-разному ориентировать относительно волчка и существуют такие положения, когда недиагональные элементы тензора инерции Iab, a.ne.b, становятся равными 0. В таком случае говорят, что координатные оси совпадают с главными осями тензора инерции (или направлены вдоль главных осей тензора инерции). Гамильтониан в таких молекулярных системах координат имеет вид:

H=A^2J_a+B^2J_b+C^2J_c
(13)

Здесь латинскими буквами, в отличие от ранее используемых греческих букв, обозначены оси молекулярной системы координат, A, B, C – это так называемые вращательные параметры, причем обозначения осей и параметров всегда выбираются таким образом, что всегда выполняются неравенства A >= B >= C. Эти постоянные связаны с компонентами тензора инерции следующими соотношениями:

A=, B=, C=
(14)

Компоненты оператора момента количества движения со штрихом есть безразмерные операторы момента импульса:

(15)

Следует отметить, пока без доказательства, что проекции оператора J на оси молекулярной системы координат обладают несколько отличными коммутационными соотношениями, так называемыми коммутационными соотношениями с аномальным знаком. Если для компонент в лабораторной системе справедливы следующие коммутационные соотношения

[J_alpha,J_beta]=
(16)

то в молекулярной системе их вид несколько другой:

[J_a,J_b]=
(17)

В выражениях (16) и (17) e_abc есть тензор Леви-Чивиты, значение компонент которого может быть определено как смешанное произведение соответствующих орт осей декартовой системы координат:

eabc=
(18)

Компоненты могут быть равны 1 и -1 если все три индекса a, b, c разные, и 0, если совпадают хотя бы два индекса.

Уровни энергии жесткого асимметричного волчка

В дальнейшем для упрощения написания формул безразмерные компоненты операторов момента импульса будем записывать без штриха:

J'=>J
(19)

Если вращательные постоянные A, B, C различны, то такой волчок называется асимметричным. Уравнение Шредингера с гамильтнианом (13) для асимметричного волчка в общем случае не имеет аналитического решения за исключением нескольких низших уровней, поэтому его приходится решать численно. Гамильтониан преобразуют в матричную форму. При этом в качестве базисных функций используют функции симметричного волчка, описанные в разделе Алгебраическая теория оператора момента количества движения. Все компоненты оператора момента количества движения в этом базисе диагональны по квантовому числу j, поэтому матрица гамильтониана представляется в виде блочно диагональной матрицы и каждый диагональный блок характризуется своим значением j. В таком случае говорят, что j хорошее (точное) квантовое число.

Для проведения расчетов и для анализа полученных результатов удобно гамильтониан преобразовать к следующей форме:

H(k),
(20)

где

H(k)
(21)

и

H(k)
(22)

есть, так называемый, параметр асимметрии Рея. Этот параметр лежит в пределах

H(k)
(23)

Удобство гамильтониана, записанного в форме (20), заключается в том, что собственные значения и собственные вектора нужно находить для матрицы (21), зависящей только от одного параметра, а не от трех, как это имеет место для гамильтониана (13), и что значение этого параметра лежит в конечных пределах.

Лекция не закончена