こんにちはー！今回も化学のお話やっていきます。今回のテーマはこちら！
動画はこちら！
動画で使ったシートはこちら(anharmonic oscillator)
では早速行きましょう!
非調和振動子とは
まず、ばねで繋がった2原子分子を考えてみましょう。
このばねが振動するわけですが、原子は質点ではなく大きさを持っているので、ばねが縮んだときに立体的な反発が生まれます。
そのため、ばねが短いときには、調和振動子のポテンシャルより実際のポテンシャルのほうが大きいという事になります。
そして、ばねの長さが無限大のときにはそれはもはや結合が解離してしまっているため、ばねの復元力ははたらかず、ポテンシャルは\(0\)になります。
つまり、結合長が長い領域では実際のポテンシャルは調和振動子より小さくなります。
それを考慮したのがこちらの図になります。
点線が調和振動子で、実線が実際のポテンシャルです。
調和振動子であればこの平衡核間距離\(r_\rm{e}\)\(\)を中心に左右対称となるのですが、実際の系では左が上がって右が下がっているので左右非対称となります。
この調和振動子からずれた分をハミルトニアンの摂動としてエネルギーを求めてみようというのが今回のテーマになります。
ちなみに、このポテンシャル図はあくまでわかりやすくしたものであり、実際のポテンシャルはもっとシャープな谷になっているのが普通なので知っておいてください。
ハミルトニアンの摂動
ではここから摂動の形を考えます。
とりあえず、ポテンシャルをマクローリン展開の形で書いてみましょう。
マクローリン展開についてはこちらの記事で解説しています。
ここで、ばねの長さが平衡核間距離に近い場合は立体反発や結合の解離が関係しないので、ほとんど調和振動子であると考えることができます。
そのため0次、1次の項が\(0\)になるとしておきます。
2次の項は\(\frac{1}{2}kx^2\)としておいて、今回は4次の項までを考えることにします。
その係数は\(a\)と\(b\)にしておきます。
ではここからエネルギーの変化分を考えていきます。
エネルギーの補正項
一次摂動エネルギー
前々回の記事で、エネルギーの補正項は補正前の波動関数から出せるということをお話ししたので、とりあえず調和振動子の波動関数を思い出してみます。
すると、このような形になっていて、\(N_n\)が規格化定数が規格化定数です。
\(H_n\)はエルミート多項式といってこの形で与えられました。
そして、実はエルミート多項式にはある漸化式が成り立っていて、それが今回重要になっていきます。
詳しい導出は割愛するので、気になる方だけ「エルミート多項式 漸化式」で検索してほしいのですが、\(\xi H_n=nH_{n-1}+\frac{1}{2}H_{n+1}\)という式が成り立ちます。
この式を使うと計算が簡単になるので、最初に少し変形をしていきます。
まず、この式の\(H_{n-1}\)と\(H_{n+1}\)をこの漸化式から書き換えます。
そして\(\xi H_n\)を表すと、こんな式になります。
両辺に\(\xi\)をかけると\(\xi ^2H_n\)を\(H_{n-2}\)と\(H_n\)と\(H_{n＋2}\)で表した式になります。
そして、さらに右辺の\(H_{n-2}\)、\(H_n\)、\(H_{n+2}\)をすべて書き換えて、また両辺に\(\xi\)をかけたものがこの式になります。
今度は\(\xi ^3H_n\)を\(H_{n-3}\)、\(H_{n-1}\)、\(H_{n+1}\)、\(H_{n+3}\)で表すことができました。
これをもう一度行うと\(\xi ^4H_n\)を\(H_{n-4}\)、\(H_{n-2}\)、\(H_n\)、\(H_{n+2}\)、\(H_{n+4}\)で表せるということになります。
それで、この関係を使えばエネルギーの補正項を考えることができます。
まず、一次の補正項は以前の記事の導出よりこれで表せます。
ここで、\(\mu\)は換算質量です。
エルミート多項式が直交性を満たすことから、\(0\)になる項がわかります。
一次摂動エネルギーについてはこちらをご覧ください。
換算質量についてはこちらをご覧ください。
まず、\(\xi ^3H_n\)について、先ほどの漸化式から求めたものを見ると、\(H_n\)がないため、直交性よりすべて消えてしまうことになります。
よって、\(\xi ^4\)部分だけが残ることになり、エネルギーの補正項をこのように求めることができます。
二次摂動エネルギー
同様に2次の補正項も考えると、これも前回の記事よりこのように書けるのですが、この場合は\(\xi ^3\)だけが残ってきます。
最終的には、このようなエネルギーの値が得られます。
モースポテンシャル
では最後に非調和な振動ポテンシャルの例としてモースポテンシャルと呼ばれるこちらの式を考えてみましょう。
ここで\(D\)がポテンシャルの最小値から解離したところのポテンシャルまでの差になっていて、\(\alpha\)はカーブの曲がり具合を表すパラメータです。
このポテンシャルをマクローリン展開した際には0次と1次の項は\(0\)になっています。
そして、2次の項の係数の関係より\(\alpha\)をこのように表すことができます。
3次、4次の係数は、導関数から下のように表せるため、\(a\)と\(b\)を求めることができます。
これをエネルギーの補正項に入れることで、数値として表すことができるようになります。
そしてエネルギーを2次まで補正した際のエネルギーの変化分は\(-\frac{\hbar ^2\omega ^2}{4D}(n+\frac{1}{2})^2\)ということになります。
調和振動子のエネルギー準位は等間隔でしたが、摂動を考えると\(n\)が大きくなるにつれて、だんだん間隔が小さくなっていくということをこの式は示しています。
まとめ
今回練習問題はありませんので、最後おさらいをやって終わります。
今回は摂動法を実際に使った例として非調和な振動子の例を紹介しました。
細かい計算をもう一度おさらいすることはありませんが、エルミート多項式の漸化式をうまく使うとエネルギーの補正項をシンプルに考えることができました。
最後には、非調和な振動ポテンシャルの例として摂動によって起こることを考えました。
その結果、調和振動子では等間隔だったエネルギー準位が等間隔ではなくなり、nが大きくなるほどその間隔が広くなっていくという結果になりました。
この効果は絶対零度付近の最安定に近い状態であれば見えませんが、高エネルギーにしていくほど顕著に表れるという事になります。
次回以降の量子化学の動画では変分法に入っていく予定ですが、その準備として行列の説明をしてから変分法に入っていきます。
それではどうもありがとうございました！
