Жорсткий ротатор

Матеріал з testwiki
Версія від 14:45, 2 лютого 2025, створена imported>A.sav (clean up, typos fixed: фунцій → функцій за допомогою AWB)
(різн.) ← Попередня версія | Поточна версія (різн.) | Новіша версія → (різн.)
Перейти до навігації Перейти до пошуку

Жорсткий ротатор — механічна модель, що використовується при описі тіл, що обертаються. Довільний жорсткий ротатор є тривимірним жорстким тілом. Прикладом може бути дзиґа. Задання орієнтації такого тіла в просторі вимагає використання кутів Ейлера. Особливим жорстким ротатором є лінійний ротатор, для опису якого потрібно тільки два кути. Прикладом може бути лінійна молекула. Загалом молекули, наприклад, молекула води (асиметричний ротатор), аміаку (симетричний ротатор) чи метану (сферичний ротатор), є тривимірними.

Абстракція жорсткого ротатора використовується в молекулярній та ядерній фізиці для моделювання обертання молекул та ядер атома, тоді коли в якомусь наближенні деформацією можна знехтувати.

Лінійний ротатор

Модель лінійного ротатора складається з двох матеріальних точок, розділених фіксованою відстанню. Ця відстань та маси точок є єдиними параметрами моделі. Однак для багатьох двоатомних молекул така модель накладає занадто великі обмеження, оскільки відстані між атомами не зовсім фіксовані. Врахування малих видовжень при моделюванні дозволяє компенсувати цей недолік. Усе ж навіть для таких молекул модель жорсткого ротатора є корисним вихідним пунктом (моделлю першого порядку).

Класичний лінійний ротатор

Класичний лінійний ротатор складається з двох матеріальних точок із масами m1 та m2 (зведена маса μ=m1m2m1+m2) на відстані R одна від одної. Ротатор жорсткий, якщо R не залежить від часу. Кінематика лінійного жорсткого ротатора зазвичай описується в сферичній системі координат. Кути задають орієнтацію ротатора в просторі. Кінетична енергія T лінійного ротатора

2T=μR2[θ˙2+(φ˙sinθ)2]=μR2(θ˙φ˙)(100sin2θ)(θ˙φ˙)=μ(θ˙φ˙)(hθ200hφ2)(θ˙φ˙),

де hθ=R та hφ=Rsinθ - множники Ламе.

Множники Ламе входять у вираз для лапласіана. У разі сталої відстані R

2=1hθhφ[θhφhθθ+φhθhφφ]=1R2[1sinθθsinθθ+1sin2θ2φ2].

Класична функція Гамільтона лінійного жорсткого ротатора має вигляд

H=12μR2[pθ2+pφ2sin2θ].

Квантовомеханічний лінійний ротатор

Модель жорсткого лінійного ротатора можна також використати в квантовій механіці для опису обертання двоатомних молекул[1]. Енергія обертання залежить від моменту інерції системи I. У системі відліку центру мас момент інерції дорівнює

I=μR2.

Квантовані рівні енергії системи знаходяться з розв'язку рівняння Шредінгера

H^Ψ=EΨ,

де Ψ — хвильова функція, а H^  — гамільтоніан. Для жорсткого ротататора у вільному просторі, гамільтоніан відповідає кінетичній енергії[2] системи:

H^=22μ2

де  — зведена стала Планка, а 2 — лапласіан. Якщо оператор Лапласа записати у сферичній системі координат, гамільтоніан задається формулою:

H^=22I[1sinθθ(sinθθ)+1sin2θ2φ2]

Аналогічний оператор фігурує також у рівнянні Шредінгера для атома водню після розділення змінних. Його власні значення та власні функції:

H^Ym(θ,φ)=22I(+1)Ym(θ,φ).

Ym(θ,φ) позначає сферичні гармоніки. Енергія не залежить від квантового числа m. Рівні енергії

E=22I(+1)

2+1 разів вироджені: функції з фіксованим , в яких m=,+1,,, мають однакову енергію.

Запроваджуючи сталу обертання B, можна записати:

E=B(+1),B22I.

В одиницях оберненої довжини стала обертання дорівнює

B¯Bhc=h8π2cI=4πcμRe2,

де c - швидкість світла. Якщо використати значення h, c, and I в гаусових одиницях, B¯ записується в обернених сантиметрах, см−1. Ця одиниця часто використовується у ротаційно-вібраційній спектроскопії. Стала обертання B¯(R) залежить від відстані R. Часто записують Be=B¯(Re), де Re рівноважне значення R (значення, при якому енергія взаємодії між атомами мінімальна).

Типовий ротаційний спектр складається із серій піків, що відповідають переходам між рівнями з різними значеннями квантового числа кутового моменту (). Відповідно, обертальні піки мають енергію, що відповідає цілому числу 2B¯.

Правила відбору

Обертальні переходи в молекулі відбуваються, коли молекула поглинає фотон (квант електромагнітного поля). Залежно від енергії фотона, тобто від довжини хвилі, це може призвести до переходу між обертальними рівнями або викликати інші переходи в молекулі (електронні чи коливальні). Чисто обертальні переходи, в яких вібронні (коливально-електронні) хвильові функції не змінюються, відбуваються у мікрохвильовому діапазоні електромагнітного спектру.

Зазвичай обертальні переходи можуть спостерігатися лише тоді, коли квантове число кутового моменту змінюється на одиницю (Δl=±1). Це правило відбору є наслідком наближення теорії збурень першого порядку в залежному від часу рівнянні Шредінгера. У рамках цього підходу обертальні переходи можна спостерігати лише тоді, коли одна чи кілька складових дипольного моменту ненульова. Якщо z - напрямок складової електричного поля електромагнітної хвилі, дипольний момент переходу дорівнює

ψ2|μz|ψ1=(μz)21=ψ2*μzψ1dτ.

Перехід відбувається, коли цей інтеграл не дорівнює нулю. Відділяючи у хвильовій функції обертальну складову від вібронної, можна показати, що це означає наявність у молекули постійного дипольного моменту. Проінтегрувавши по вібронним координатам, залишається обертальна частина моменту переходу

(μz)l,m;l,m=μ02πdϕ0πYlm(θ,ϕ)*cosθYlm(θ,ϕ)dcosθ.

Тут μcosθ - z-ва компонента постійного дипольного моменту. Момент μ є усередниним по вібронних координатах дипольним оператором. Залишається тільки компонента постійного моменту вздовж осі різноатомних молекул. Використовуючи, ортогональність сферичних гармонік Ylm(θ,ϕ) можна визначити значення l, m, l, та m, що даватимуть ненульові значення моменту переходу. Це обмеження визначає правила відбору для жорсткого ротатора

Δm=0,Δl=±1.

Нежорсткий лінійний ротатор

Зазвичай для опису обертальної енергії двоатомних молекул використовують модель жорсткого ротатора, але вона не є абсолютно точним описом таких молекул. Це зумовлено тим, що молекулярний зв'язок, а з ним і міжатомна відстань R не є точно фіксованими. При швидшому обертанні (більших значеннях квантового числа l) зв'язок розтягається. Врахувати цей ефект за рахунок множника, відомого як стала відцентрового спотворення D¯ (риска над змінною означає, що її величина виражається у см−1):

E¯l=Elhc=B¯l(l+1)D¯l2(l+1)2

де

D¯=4B¯3ω¯2.

ω¯ є фундаментальною частотою коливань молекули (в см−1). Ця частота зв'язана зі зведеною масою та коефіцієнтом жорсткості зв'язку за формулою

ω¯=12πckμ.

Попри те, що модель нежорсткого ротататора задовільна для двоатомної молекули, вона все ж недосконала. Причина в тому, що вона не враховує розтягання зв'язку за рахунок енергії у ньому (ангармонічність потенціалу).

Жорсткий ротатор довільної форми

Жорсткий ротатор довільної форми є абсолютно твердим тілом з фіксованим центром маси (він може також рівномірно рухатися) у вільному від полів просторі R3, тож його енергія скаладається з кінетичної (при рівномірному русі є також стала складова, яку можна не враховувати). Жорстке тіло можна характеризувати трьома власними значеннями тензора моменту інерції, що є невід'ємними числами, які називають головними моментами інерції. У мікрохвильовій спектроскопії, яка в основному вивчає обертальні переходи, зазвичай молекули класифікують так:

  • сферичні ротатори
  • симетричні ротатори
    • сплюснуті симетричні ротатори
    • витягнуті симетричні ротатори
  • асиметричні ротатори

Ця класифікація залежить від віденосних значень головних моментів інерції.

Координати

Різні підрозділи фізики та області техніки використовують різні засоби для опису кінематики жорсткого ротатора. Молекулярна фізика майже винятково користується кутами Ейлера. У квантовій механіці використання кутів Ейлера теж має переваги, оскільки вони є простим узагальненням сферичної системи координат.

Першим кроком є уявне закріплення на тілі ортогональної системи координат. Цю систему координат можна закріпити на тілі будь-як, але здебільшого для цього використовують систему, а якій тензор інерції діагональний, тобто осі кооридинат збігаються з головними осями тензора інерції. Така система завжди ортогональна, оскільки тензор інерції задається Ермітовою матрицею. Коли ротатор має вісь симетрії, вона зазвичай збігається з одною із головних осей. Зручно обрати вісь симетрії найвищого порядку за вісь z.

У початковий момент часу лабораторну систему суміщають із закріпленою на тілі, тож осі x, y та z закріпленої системи збігаються з осями X, Y та Z просторової лабораторної системи. Потім тіло повертають на певний додатній кут α навколо осі z (проти годинникової стрілки), що змусить вісь y зміститися в положення y. Третім кроком, тіло та закріплену на ньому систему координат повертають на кут β навколо осі y. Вісь z закріпленої на тілі системи після цих двох поворотів матиме стосовно непорушної просторової системи полярний кут α (зазвичай його позначають φ) та азимутальний кут β (зазвичай його позначають θ). Якщо ротатор циліндричний вздовж осі z, на кшталт лінійного жорсткого ротатора, його положення в просторі є однозначно визначеним.

Якщо циліндричної симетрії нема, то останній поворот проводиться навколо осі z-axis (з кутовими координатами β та α ), що необхідно для повного визначення орієнтації тіла. Традиційно останній кут повороту називають γ (або ψ).

Така конвенція визначення кутів Ейлера відома як конвенція zyz; можна показати, що вона аналогічна конвенції zyz, у якій порядок поворотів обернений.

Повна матриця трьох послідовних поворотів є добутком

𝐑(α,β,γ)=(cosαsinα0sinαcosα0001)(cosβ0sinβ010sinβ0cosβ)(cosγsinγ0sinγcosγ0001)

Нехай 𝐫(0) є радіус-вектором довільної точки тіла 𝒫 в закріпленій на тілі системі відліку. Початково 𝐫(0) є також радіус-вектром точки 𝒫 в просторовій системі координат. Повороти тіла не змінюють координати в системі, закріпленій на тілі, але змінюють просторові координати, тож вектор 𝒫 стає

𝐫(α,β,γ)=𝐑(α,β,γ)𝐫(0).

Зокрема, якщо точка 𝒫 лежить спочатку на осі Z, її координати стають

𝐑(α,β,γ)(00r)=(rcosαsinβrsinαsinβrcosβ),

що демонструє відповідність сферичній системі координат (у конвенції, що використовується в фізиці).

Знання кутів Ейлера як функції часу t та початкових координат 𝐫(0) визначає кінематику обертання абсолютно твердого тіла.

Класична кінетична енергія

Це є узагальненням добре відомого виразу для енергії обертання тіла навколо одної осі.

Тут припускається, що закріплена на тілі система координат є системою головних осей; вона діагоналізує миттєве значення тензора моменту інерції 𝐈(t) , тобто

𝐑(α,β,γ)1𝐈(t)𝐑(α,β,γ)=𝐈(0), де 𝐈(0)=(I1000I2000I3),

Кути Ейлера вважаються тут залежними від часу, що в свою чергу визначає залежність від часу 𝐈(t). Це позначення означає, що при t=0 кути Ейлера нульові, тобто при t=0 закріплена на тілі система відліку збігається з просторовою.

Класична кінетична енергія T жорсткого ротатора може бути записана по різному:

  • як функція кутових швидкостей
  • у лагранжевій формі
  • як функція кутового моменту
  • у гамільтоновій формі

Оскільки кожна з цих форм використовується і записана в підручниках, тут буде наведено усі.

Через кутові швидкості

T, виражена через кутові швидкості, має вигляд

T=12[I1ωx2+I2ωy2+I3ωz2],

де

(ωxωyωz)=(sinβcosγsinγ0sinβsinγcosγ0cosβ01)(α˙β˙γ˙).

Вектор ω=(ωx,ωy,ωz) складений із компонент кутових швидкостей ротатора, щодо осей, закріплених на тілі. Можна показати, що ω не є похідною будь-якого вектора, на відміну від звичного означння швидкості[3]. Крапки над ейлеровими кутами означають часові похідні в нотації Ньютона. Кутові швидкості задовольняють систему рівнянь, відому під назвою ейлерових (із нульовим моментом сили, оскільки вважається, що ротатор обертається у просторі, вільному від сил).

Форма Лагранжа

Якщо знову підставити вираз для ω у T, можна отримати кінетичну енергію в формі функції Лагранжа (як функцію похідних від ейлерових кутів). У матричному запису

2T=(α˙β˙γ˙)𝐠(α˙β˙γ˙),

де 𝐠 - метричний тензор, виражений через ейлерові кути (в криволінійних координатах);

𝐠=(I1sin2βcos2γ+I2sin2βsin2γ+I3cos2β(I2I1)sinβsinγcosγI3cosβ(I2I1)sinβsinγcosγI1sin2γ+I2cos2γ0I3cosβ0I3).

Через кутовий момент

Доволі часто кінетичну енергію записують через кутовий момент 𝐋. В обертовій системі відліку він має компоненти Li. Можна показати, що він зв'язаний з кутовою швидкістю,

𝐋=𝐈(0)ω, або Li=Tωi,i=x,y,z.

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

Вираз для кінетичної енергії через кутовий момент має вигляд

T=12[Lx2I1+Ly2I2+Lz2I3].

Гамільтонова форма

У гамільтоновій формі кінетична енергія записується через узагальнені імпульси, що визначаються як

(pαpβpγ) =def (T/α˙T/β˙T/γ˙)=𝐠(α˙β˙γ˙),

де використано симетричність матриці 𝐠.

Сам вираз для кінетичної енергії має вигляд

2T=(pαpβpγ)𝐠1(pαpβpγ),

де обернений метричний тензор задається як

sin2β𝐠1=
(cos2γI1+sin2γI2(1I21I1)sinβsinγcosγcosβcos2γI1cosβsin2γI2(1I21I1)sinβsinγcosγsin2βsin2γI1+sin2βcos2γI2(1I11I2)sinβcosβsinγcosγcosβcos2γI1cosβsin2γI2(1I11I2)sinβcosβsinγcosγcos2βcos2γI1+cos2βsin2γI2+sin2βI3).

Цей обернений тензор потрібен для отримання оператора Лапласа — Бельтрамі, який (помножений на 2) задає оператор енергії жорсткого ротатора в квантовій механіці.

Наведену класичну функцію Гамільтона можна переписати у вигляді, необхідному для інтегрування у фазовому просторі класичної статистичної механіки

T=12I1sin2β((pαpγcosβ)cosγpβsinβsinγ)2+12I2sin2β((pαpγcosβ)sinγ+pβsinβcosγ)2+pγ22I3.

Жорсткий ротатор у квантовій механіці

Шаблон:See also

Як зазвичай перехід до квантової механіки здійснюється заміною узагальнених імпульсів на оператори, в яких фігурують похідні щодо канонічно спряжених узагальнених координат. Так,

pαiα,

і аналогічно щодо pβ та pγ. На диво, це правило зводить доволі складну функцію кутів Ейлера, їхніх похідних та моментів інерції до простого диференціального оператора, який не залежить від часу чи від моментів інерції, і в якому фігурує похідна тільки від одного з кутів Ейлера.

Цього правила квантування достатньо, щоб отримати оператори, що відповідають класичним кутовим моментам. Існує два види кутових моментів: ті, що визначаються у фіксованій просторовій системі відліку, та ті, що зв'язані з тілом. Як і ті, так і іншу є векторними операторами, тобто мають три складові, що перетворюються при обертанні (фіксованої проторової системи та системи, зв'язаної з тілом, відповідно) як вектори. Точна форма кутових моментів жорсткого ротататора задається D-матрицею Вігнера (яку, втім, слід помножити на ). У системі тіла оператори кутового моменту записуються як 𝒫^i. Вони мають незвичні комутаційні співвідношення.

Правила квантування не вистачає, щоб записати оператор кінетичної енергії, виходячи з класичного гамільтоніана. Оскільки в класичній фізиці pβ комутує з cosβ та sinβ й оберненими до них функціями, положення цих тригонометричних функцій у класичній функції Гамільтона довільне. Після квантування комутації уже більше нема, і порядок операторів та функцій в гамільтоніані важко визначити. Подольський[2] у 1928 запропонував використовувати оператор Лапласа — Бельтрамі (помножений на 122), що має відповідну форму для того, щоб грати роль квантомеханічного оператора енергії. Форма цього оператора така (використовується конвенція підсумовування: повторення індексів означає суму— у цьому випадку це стосується трьох кутів Ейлера q1,q2,q3α,β,γ):

H^=22|g|1/2qi|g|1/2gijqj,

де |g| позначає детермінант тензора g:

|g|=I1I2I3sin2β. gij=(𝐠1)ij.

Отримавши обернений метричний тензор, можна записати оператор кінетичної енергії через ейлерові кути за допомогою простої підстановки (зауваження: відповідне рівняння на власні значення є рівнянням Шредінгера для жорсткого ротатора у формі, для якої було вперше знайдено розв'язок Кронігом та Рабі[4] (для випадку симетричного ротатора). Це один із небагатьох випадків, коли рівняння Шредінгера можна розв'язати аналітично. Усі розв'язки було знайдено упродовж року після формулювання рівняння Шредінгера.)

Зараз зазвичай поступають так: можна показати, що H^ виражається в обертовій системіі відліку через оператори кутового моменту (при доведенні потрібно бути обережним щодо комутації диференційних операторів та тригонометричних функцій). Результат виглядає аналогічно класичному:

H^=12[𝒫x2I1+𝒫y2I2+𝒫z2I3].

Дія операторів 𝒫^i на D-матриці Вігнера проста. Зокрема

𝒫2Dmmj(α,β,γ)*=2j(j+1)Dmmj(α,β,γ)*, 𝒫2=𝒫x2+𝒫y2+𝒫z2,

тож рівняння Шредінгера для сферичного ротора (I=I1=I2=I3) розв'язується, визначаєть (2j+1)2 вироджені рівні енергії 2j(j+1)2I.

Симетричний ротатор має I1=I2. Він вигягнутий (як сигара), якщо I3<I1=I2. У цьому випадку гамільтоніан записується

H^=12[𝒫2I1+𝒫z2(1I31I1)],

і використовується

𝒫z2Dmkj(α,β,γ)*=2k2Dmkj(α,β,γ)*.

Отже,

H^Dmkj(α,β,γ)*=EjkDmkj(α,β,γ)*,, деEjk/2=j(j+1)2I1+k2(12I312I1).

Власне значення Ej0 вироджено 2j+1 разів, усі власні функції з m=j,j+1,,j мають однакове власне значення. Енергії з |k| > 0 вироджені 2(2j+1) разів. Цей точний розв'язок рівняння Шредінгера для симетричного ротатора було знайдено в 1927[4].

Задача про обертання несиметричного ротатора (I1I2I3) точно не розв'язується.

Пряме експериментальне спостереження обертань молекул

Обертання молекули впродовж довгого часу безпосередньо спостерігати не вдавалося. Тільки розробка технології вимірювання з атомною роздільністю зробила його можливим[5][6]. При низьких температурах обертання молекул можна заморозити, повністю або частково. Тоді його можна бачити в тунельних мікроскопах, стабілізацію при вищих температурах можна пояснити обертальною ентропією[6].

Пряме спостереження збудження обертання у виокремленій одиночній молекулі вдалося побачити нещодавно, використовуючи спектроскопію непружного тунелювання електронів у тунельному мікроскопі[7][8].


Виноски

Шаблон:Reflist

  1. Шаблон:Cite book
  2. 2,0 2,1 Шаблон:Cite journal
  3. Chapter 4.9 of Goldstein, H.; Poole, C. P.; Safko, J. L. (2001). Classical Mechanics (Third ed.). San Francisco: Addison Wesley Publishing Company. ISBN 0-201-65702-3.
  4. 4,0 4,1 Шаблон:Cite journal
  5. Шаблон:Citation/core
  6. 6,0 6,1 Шаблон:Citation/core
  7. S. Li, A. Yu, A, F. Toledo, Z. Han, H. Wang, H. Y. He, R. Wu, and W. Ho, Phys. Rev. Lett. 111, 146102 (2013).http://journals.aps.org/prl/abstract/10.1103/PhysRevLett.111.146102
  8. F. D. Natterer, F. Patthey, and H. Brune, Phys. Rev. Lett. 111, 175303 (2013).http://journals.aps.org/prl/abstract/10.1103/PhysRevLett.111.175303