モースポテンシャル

モースポテンシャル（英: Morse potential）は、二原子分子の原子間相互作用を表現するのに便利なポテンシャルである。名称は物理学者フィリップ・M・モース（英語版）にちなむ。このポテンシャルによる原子間振動のモデルは、切断の効果（非結合状態）を明示的に含むため調和振動子（英語版）（quantum harmonic oscillator、QHO）とみなすよりも良い近似となる。さらに振動の非調和性や（赤外分光における）倍音（overtone）、結合音（英語版）によるバンド（hot band）の存在をも説明できる。 モースポテンシャルは他の相互作用（例えば原子・界面間）の記述にも用いられる。その単純性（3つのパラメータしか含まない）から最新の分光学では用いられなくなっているが、その数学的形式はモース長距離ポテンシャル（Morse/Long-range, MLR）へと発展させられた。これは分光データに当てはめるのに最もよく利用されるポテンシャルである。 モースポテンシャルの関数形は V ′ ( r ) = D e ( 1 − e − a ( r − r e ) ) 2 {\displaystyle V'(r)=D_{\mathrm {e} }\left(1-e^{-a(r-r_{\mathrm {e} })}\right)^{2}} である。ここで r {\displaystyle r} は原子間距離、 r e {\displaystyle r_{\mathrm {e} }} は平衡結合距離、 D e {\displaystyle D_{\mathrm {e} }} はポテンシャルの井戸の深さ（well depth, 解離状態の原子に基づいて定められる）、 a {\displaystyle a} はポテンシャルの幅を調整する定数である（ a {\displaystyle a} が小さいほど井戸は広くなる）。結合解離エネルギーは井戸の深さ D e {\displaystyle D_{\mathrm {e} }} から零点エネルギー E 0 {\displaystyle E_{0}} を引くことで計算できる。結合の強度を見るには V ′ ( r ) {\displaystyle V'(r)} を r = r e {\displaystyle r=r_{\mathrm {e} }} のまわりでのテイラー展開すればよい。2階導関数が D e a 2 {\displaystyle D_{\mathrm {e} }a^{2}} となるので、パラメータ a {\displaystyle a} は a = k e / 2 D e {\displaystyle a={\sqrt {k_{\mathrm {e} }/2D_{\mathrm {e} }}}} と表せる。ここで k e {\displaystyle k_{\mathrm {e} }} は r = r e {\displaystyle r=r_{\mathrm {e} }} 付近で調和振動とみなしたときのばね定数である。 ポテンシャルエネルギーの基準点のとり方には任意性があるから、定数の加減によりモースポテンシャル関数形は何通りにも書くことができる。原子・界面相互作用の場合は、基準点を調節して V ( r ) = V ′ ( r ) − D e = D e ( 1 − e − a ( r − r e ) ) 2 − D e {\displaystyle V(r)=V'(r)-D_{\mathrm {e} }=D_{\mathrm {e} }\left(1-e^{-a(r-r_{\mathrm {e} })}\right)^{2}-D_{\mathrm {e} }} のように改良できる。これは V ( r ) = D e ( e − 2 a ( r − r e ) − 2 e − a ( r − r e ) ) {\displaystyle V(r)=D_{\mathrm {e} }\left(e^{-2a(r-r_{\mathrm {e} })}-2e^{-a(r-r_{\mathrm {e} })}\right)} と書かれることが多い。この場合 r {\displaystyle r} は界面から垂直に測った位置を表し、ポテンシャルは r {\displaystyle r} が無限大のときゼロ、極小点（つまり r = r e {\displaystyle r=r_{\mathrm {e} }} ）のとき − D e {\displaystyle -D_{\mathrm {e} }} をとる。この形に書くと、モースポテンシャルによって与えられる力が短距離での斥力（第1項）と長距離での引力（第2項）との組み合わせであることが明瞭になる。この点はレナード-ジョーンズ・ポテンシャルと類似している。 量子的調和振動子のときと同様に、モースポテンシャルの固有エネルギーおよびエネルギー固有状態は演算子の代数的処理により求まる。その一法はハミルトニアンを因子分解するものである。 モースポテンシャルに対する定常状態を得るため、次のシュレーディンガー方程式： ( − ℏ 2 2 m ∂ 2 ∂ r 2 + V ( r ) ) Ψ n ( r ) = E n Ψ n ( r ) {\displaystyle \left(-{\frac {\hbar ^{2}}{2m}}{\frac {\partial ^{2}}{\partial r^{2}}}+V(r)\right)\Psi _{n}(r)=E_{n}\Psi _{n}(r)} を満たす Ψ n ( r ) {\displaystyle \Psi _{n}(r)} と E n {\displaystyle E_{n}} を求めたい。次のように新しい変数を導入すると： x = a r ; x e = a r e ; λ = 2 m D e a ℏ ; ε n = 2 m a 2 ℏ 2 E n {\displaystyle x=ar;x_{\mathrm {e} }=ar_{\mathrm {e} };\lambda ={\frac {\sqrt {2mD_{\mathrm {e} }}}{a\hbar }};\varepsilon _{n}={\frac {2m}{a^{2}\hbar ^{2}}}E_{n}} シュレーディンガー方程式は次の簡単な形になる： ( − ∂ 2 ∂ x 2 + V ( x ) ) Ψ n ( x ) = ε n Ψ n ( x ) {\displaystyle \left(-{\frac {\partial ^{2}}{\partial x^{2}}}+V(x)\right)\Psi _{n}(x)=\varepsilon _{n}\Psi _{n}(x)} V ( x ) = λ 2 ( e − 2 ( x − x e ) − 2 e − ( x − x e ) ) {\displaystyle V(x)=\lambda ^{2}\left(e^{-2\left(x-x_{\mathrm {e} }\right)}-2e^{-\left(x-x_{\mathrm {e} }\right)}\right)} この固有エネルギーおよびエネルギー固有状態は次のように書ける。 ε n = λ 2 − ( λ − n − 1 2 ) 2 = 2 λ ( n + 1 2 ) − ( n + 1 2 ) 2 {\displaystyle \varepsilon _{n}=\lambda ^{2}-\left(\lambda -n-{\frac {1}{2}}\right)^{2}=2\lambda \left(n+{\frac {1}{2}}\right)-\left(n+{\frac {1}{2}}\right)^{2}} Ψ n ( z ) = N n z λ − n − 1 / 2 e − z / 2 L n ( 2 λ − 2 n − 1 ) ( z ) {\displaystyle \Psi _{n}(z)=N_{n}z^{\lambda -n-1/2}e^{-z/2}L_{n}^{(2\lambda -2n-1)}(z)} ここで z = 2 λ e − ( x − x e ) ;  N n = [ n ! ( 2 λ − 2 n − 1 ) Γ ( 2 λ − n ) ] 1 2 {\displaystyle z=2\lambda e^{-\left(x-x_{\mathrm {e} }\right)}{\text{; }}N_{n}=\left[{\frac {n!\left(2\lambda -2n-1\right)}{\Gamma (2\lambda -n)}}\right]^{\frac {1}{2}}} 、また L n ( α ) ( z ) {\displaystyle L_{n}^{(\alpha )}(z)} は一般化ラゲール多項式である。 L n ( α ) ( z ) = z − α e z n ! d n d z n ( z n + α e − z ) = Γ ( α + n + 1 ) / Γ ( α + 1 ) n ! 1 F 1 ( − n , α + 1 , z ) {\displaystyle L_{n}^{(\alpha )}(z)={\frac {z^{-\alpha }e^{z}}{n!}}{\frac {d^{n}}{dz^{n}}}\left(z^{n+\alpha }e^{-z}\right)={\frac {\Gamma (\alpha +n+1)/\Gamma (\alpha +1)}{n!}}\,_{1}F_{1}(-n,\alpha +1,z)} さらに、以下のような位置演算子の行列要素の解析的表現も重要である（ここで l > n , N = λ − 1 / 2 {\displaystyle l>n,N=\lambda -1/2} とする）。 ⟨ Ψ l | x | Ψ n ⟩ = 2 ( − 1 ) l − n + 1 ( l − n ) ( 2 N − n − l ) ( N − n ) ( N − l ) Γ ( 2 N − l + 1 ) l ! Γ ( 2 N − n + 1 ) n ! {\displaystyle \left\langle \Psi _{l}|x|\Psi _{n}\right\rangle ={\frac {2(-1)^{l-n+1}}{(l-n)(2N-n-l)}}{\sqrt {\frac {(N-n)(N-l)\Gamma (2N-l+1)l!}{\Gamma (2N-n+1)n!}}}} 固有エネルギーを元々の変数で書くと： E n = h ν 0 ( n + 1 / 2 ) − [ h ν 0 ( n + 1 / 2 ) ] 2 4 D e {\displaystyle E_{n}=h\nu _{0}(n+1/2)-{\frac {\left[h\nu _{0}(n+1/2)\right]^{2}}{4D_{\mathrm {e} }}}} となる。ここで h {\displaystyle h} はプランク定数、 n {\displaystyle n} は振動状態を表す量子数、 ν 0 {\displaystyle \nu _{0}} は振動数の単位量であり、質量 m {\displaystyle m} とポテンシャルのパラメータを使って次のように書ける。 ν 0 = a 2 π 2 D e / m {\displaystyle \nu _{0}={\frac {a}{2\pi }}{\sqrt {2D_{\mathrm {e} }/m}}} 量子的調和振動子の振動準位の間隔は一定値 h ν 0 {\displaystyle h\nu _{0}} であったが、モースポテンシャルによる振動子の固有エネルギーの間隔は n {\displaystyle n} の増加とともに減少してゆく。数学的には E n + 1 − E n = h ν 0 − ( n + 1 ) ( h ν 0 ) 2 / 2 D e {\displaystyle E_{n+1}-E_{n}=h\nu _{0}-(n+1)(h\nu _{0})^{2}/2D_{\mathrm {e} }} と書ける。この傾向は現実的に観測される分子の非調和性とも一致する。しかし、この式はある値 n m {\displaystyle n_{m}} を超えたところで E n + 1 − E n {\displaystyle E_{n+1}-E_{n}} がゼロまたは負になり破綻し、具体的には n m = ⌊ 2 D e − h ν 0 h ν 0 ⌋ {\displaystyle n_{m}=\left\lfloor {\frac {2D_{\mathrm {e} }-h\nu _{0}}{h\nu _{0}}}\right\rfloor } である。これはモースポテンシャルの下では有限個の束縛状態しかとれないことに由来する。 n m {\displaystyle n_{m}} を上回るところでは任意のエネルギーをとることが可能になり、上に掲げた E n {\displaystyle E_{n}} の式は成り立たなくなる。 n m {\displaystyle n_{m}} 以下であれば、 E n {\displaystyle E_{n}} は回転を考慮しない二原子分子の実際の振動構造をよく近似する。実際、現実の分子のスペクトルに対し次式を一般に当てはめ得る。 E n / h c = ω e ( n + 1 / 2 ) − ω e χ e ( n + 1 / 2 ) 2 {\displaystyle E_{n}/hc=\omega _{\mathrm {e} }(n+1/2)-\omega _{\mathrm {e} }\chi _{\mathrm {e} }(n+1/2)^{2}\,} ここで定数 ω e {\displaystyle \omega _{\mathrm {e} }} と χ e {\displaystyle \chi _{\mathrm {e} }} はモースポテンシャルのパラメータと直接的に結び付けることができる。 次元解析から明らかなように、歴史的な理由からこの等式では ω e {\displaystyle \omega _{\mathrm {e} }} は E = ℏ ω {\displaystyle E=\hbar \omega } となる角周波数ではなく、 E = h c ω {\displaystyle E=hc\omega } を満たすような波数を表している。 モースポテンシャルの重要な拡張がモース長距離ポテンシャルであり（これによってモース型の関数は現代的な分光学において非常に有用なものとなった）、二原子分子の分光データやビリアル係数等のデータを表現するのに標準的に用いられる。これまでに N2、Ca2、KLi、MgH、いくつかの電子状態に対してのLi2、Cs2、Sr2、ArXe、LiCa、LiNa、Br2、Mg2、HF、HCl、HBr、HI、MgD、Be2、BeH、NaHといった適用例がある。多原子分子に対してはさらに洗練されたポテンシャルが用いられる。 ^ F. Cooper, A. Khare, U. Sukhatme (2001). Supersymmetry in Quantum Mechanics. World Scientific. pp. Table 4.1. ISBN 978-9810246129  ^ Dahl, J.P.; Springborg, M. (1988). “The Morse Oscillator in Position Space, Momentum Space, and Phase Space”. J. Chem. Phys. 88: 4535. Bibcode: 1988JChPh..88.4535D. doi:10.1063/1.453761. http://orbit.dtu.dk/files/3620619/Dahl.pdf.  ^ Lima, Emanuel F de; Hornos, José E M (2005). “Matrix elements for the Morse potential under an external field”. J. Phys. B: At. Mol. Opt. Phys. 38 (7): 815–825. doi:10.1088/0953-4075/38/7/004.  ^ a b Le Roy, Robert J.; N. S. Dattani; J. A. Coxon; A. J. Ross; Patrick Crozet; C. Linton (25 November 2009). “Accurate analytic potentials for Li2(X) and Li2(A) from 2 to 90 Angstroms, and the radiative lifetime of Li(2p)”. Journal of Chemical Physics 131 (20): 204309. Bibcode: 2009JChPh.131t4309L. doi:10.1063/1.3264688.  ^ Le Roy, R. J.; Y. Huang; C. Jary (2006). “An accurate analytic potential function for ground-state N2 from a direct-potential-fit analysis of spectroscopic data”. Journal of Chemical Physics 125 (16): 164310. Bibcode: 2006JChPh.125p4310L. doi:10.1063/1.2354502.  ^ Le Roy, Robert J.; R. D. E. Henderson (2007). “A new potential function form incorporating extended long-range behaviour: application to ground-state Ca2”. Molecular Physics 105 (5–7): 663–677. Bibcode: 2007MolPh.105..663L. doi:10.1080/00268970701241656.  ^ Salami, H.; A. J. Ross; P. Crozet; W. Jastrzebski; P. Kowalczyk; R. J. Le Roy (2007). “A full analytic potential energy curve for the a3Σ+ state of KLi from a limited vibrational data set”. Journal of Chemical Physics 126 (19): 194313. Bibcode: 2007JChPh.126s4313S. doi:10.1063/1.2734973.  ^ a b Henderson, R. D. E.; A. Shayesteh; J. Tao; C. Haugen; P. F. Bernath; R. J. Le Roy (4 October 2013). “Accurate Analytic Potential and Born–Oppenheimer Breakdown Functions for MgH and MgD from a Direct-Potential-Fit Data Analysis”. The Journal of Physical Chemistry A 117 (50): 131028105904004. Bibcode: 2013JPCA..11713373H. doi:10.1021/jp406680r.  ^ a b c Le Roy, R. J.; C. C. Haugen; J. Tao; H. Li (February 2011). “Long-range damping functions improve the short-range behaviour of 'MLR' potential energy functions”. Molecular Physics 109 (3): 435–446. Bibcode: 2011MolPh.109..435L. doi:10.1080/00268976.2010.527304. http://scienide2.uwaterloo.ca/~rleroy/Pubn/11MolP_damping.pdf.  ^ Shayesteh, A.; R. D. E. Henderson; R. J. Le Roy; P. F. Bernath (2007). “Ground State Potential Energy Curve and Dissociation Energy of MgH”. The Journal of Physical Chemistry A 111 (49): 12495–12505. Bibcode: 2007JPCA..11112495S. doi:10.1021/jp075704a. PMID 18020428.  ^ Dattani, N. S.; R. J. Le Roy (8 May 2013). “A DPF data analysis yields accurate analytic potentials for Li2(a) and Li2(c) that incorporate 3-state mixing near the c-state asymptote”. Journal of Molecular Spectroscopy (Special Issue) 268: 199–210. arXiv:1101.1361. Bibcode: 2011JMoSp.268..199.. doi:10.1016/j.jms.2011.03.030.  ^ a b W. Gunton, M. Semczuk, N. S. Dattani, K. W. Madison, High resolution photoassociation spectroscopy of the 6Li2 A-state, https://arxiv.org/abs/1309.5870 ^ Semczuk, M.; Li, X.; Gunton, W.; Haw, M.; Dattani, N. S.; Witz, J.; Mills, A. K.; Jones, D. J. et al. (2013年). “High-resolution photoassociation spectroscopy of the 6Li2 c-state”. Phys. Rev. A 87: pp. 052505. Bibcode: 2013PhRvA..87e2505S. doi:10.1103/PhysRevA.87.052505  ^ Xie, F.; L. Li; D. Li; V. B. Sovkov; K. V. Minaev; V. S. Ivanov; A. M. Lyyra; S. Magnier (2011). “Joint analysis of the Cs2 a-state and 1 g (33Π1g ) states”. Journal of Chemical Physics 135 (2): 02403. Bibcode: 2011JChPh.135b4303X. doi:10.1063/1.3606397.  ^ Coxon, J. A.; P. G. Hajigeorgiou (2010). “The ground X 1Σ+g electronic state of the cesium dimer: Application of a direct potential fitting procedure”. Journal of Chemical Physics 132 (9): 094105. Bibcode: 2010JChPh.132i4105C. doi:10.1063/1.3319739.  ^ Stein, A.; H. Knockel; E. Tiemann (April 2010). “The 1S+1S asymptote of Sr2 studied by Fourier-transform spectroscopy”. The European Physical Journal D 57 (2): 171–177. arXiv:1001.2741. Bibcode: 2010EPJD...57..171S. doi:10.1140/epjd/e2010-00058-y.  ^ Piticco, Lorena; F. Merkt; A. A. Cholewinski; F. R. W. McCourt; R. J. Le Roy (December 2010). “Rovibrational structure and potential energy function of the ground electronic state of ArXe”. Journal of Molecular Spectroscopy 264 (2): 83–93. Bibcode: 2010JMoSp.264...83P. doi:10.1016/j.jms.2010.08.007.  ^ Ivanova, Milena; A. Stein; A. Pashov; A. V. Stolyarov; H. Knockel; E. Tiemann (2011). “The X2Σ+ state of LiCa studied by Fourier-transform spectroscopy”. Journal of Chemical Physics 135 (17): 174303. Bibcode: 2011JChPh.135q4303I. doi:10.1063/1.3652755.  ^ Steinke, M.; H. Knockel; E. Tiemann (27 April 2012). “X-state of LiNa studied by Fourier-transform spectroscopy”. Physical Review A 85 (4): 042720. Bibcode: 2012PhRvA..85d2720S. doi:10.1103/PhysRevA.85.042720.  ^ Yukiya, T.; N. Nishimiya; Y. Samejima; K. Yamaguchi; M. Suzuki; C. D. Boonec; I. Ozier; R. J. Le Roy (January 2013). “Direct-potential-fit analysis for the system of Br2”. Journal of Molecular Spectroscopy 283: 32–43. Bibcode: 2013JMoSp.283...32Y. doi:10.1016/j.jms.2012.12.006. http://www.sciencedirect.com/science/article/pii/S0022285212002354.  ^ Knockel, H.; S. Ruhmann; E. Tiemann (2013). “The X-state of Mg2 studied by Fourier-transform spectroscopy”. Journal of Chemical Physics 138 (9): 094303. Bibcode: 2013JChPh.138i4303K. doi:10.1063/1.4792725.  ^ a b c d Li, Gang; I. E. Gordon; P. G. Hajigeorgiou; J. A. Coxon; L. S. Rothman (July 2013). “Reference spectroscopic data for hydrogen halides, Part II:The line lists”. Journal of Quantitative Spectroscopy & Radiative Transfer 130: 284–295. Bibcode: 2013JQSRT.130..284L. doi:10.1016/j.jqsrt.2013.07.019.  ^ a b c d “Improved direct potential fit analyses for the ground electronic states of the hydrogen halides: HF/DF/TF, HCl/DCl/TCl, HBr/DBr/TBr and HI/DI/TI”. Journal of Quantitative Spectroscopy and Radiative Transfer 151: 133–154. Bibcode: 2015JQSRT.151..133C. doi:10.1016/j.jqsrt.2014.08.028.  ^ Meshkov, Vladimir V.; Stolyarov, Andrey V.; Heaven, Michael C.; Haugen, Carl; LeRoy, Robert J. (2014). Direct-potential-fit analyses yield improved empirical potentials for the ground X1Σg+ state of Be2. 140. pp. 064315. doi:10.1063/1.4864355.  ^ “Beryllium monohydride (BeH): Where we are now, after 86 years of spectroscopy”. J. Mol. Spectrosc. 311: 76–83. arXiv:1408.3301. Bibcode: 2015JMoSp.311...76D. doi:10.1016/j.jms.2014.09.005.  ^ Walji, Sadru-Dean; Sentjens, Katherine M.; Le Roy, Robert J. (2015). “Dissociation energies and potential energy functions for the ground X1Σ+ and “avoided-crossing” A1Σ+ states of NaH”. J. Chem. Phys. 142 (4): 044305. doi:10.1063/1.4906086.  CRC Handbook of chemistry and physics, Ed David R. Lide, 87th ed, Section 9, SPECTROSCOPIC CONSTANTS OF DIATOMIC MOLECULES pp. 9–82 Morse, P. M. (1929年). “Diatomic molecules according to the wave mechanics. II. Vibrational levels”. Phys. Rev. 34: pp. 57–64. Bibcode: 1929PhRv...34...57M. doi:10.1103/PhysRev.34.57  Girifalco, L. A.; Weizer, G. V. (1959年). “Application of the Morse Potential Function to cubic metals”. Phys. Rev. 114 (3): p. 687. Bibcode: 1959PhRv..114..687G. doi:10.1103/PhysRev.114.687  Shore, Bruce W. (1973年). “Comparison of matrix methods applied to the radial Schrödinger eigenvalue equation: The Morse potential”. J. Chem. Phys. 59 (12): p. 6450. Bibcode: 1973JChPh..59.6450S. doi:10.1063/1.1680025  Keyes, Robert W. (1975年). “Bonding and antibonding potentials in group-IV semiconductors”. Phys. Rev. Lett. 34 (21): pp. 1334–1337. Bibcode: 1975PhRvL..34.1334K. doi:10.1103/PhysRevLett.34.1334  Lincoln, R. C.; Kilowad, K. M.; Ghate, P. B. (1967年). “Morse-potential evaluation of second- and third-order elastic constants of some cubic metals”. Phys. Rev. 157 (3): pp. 463–466. Bibcode: 1967PhRv..157..463L. doi:10.1103/PhysRev.157.463  Dong, Shi-Hai; Lemus, R.; Frank, A. (2001年). “Ladder operators for the Morse potential”. Int. J. Quant. Chem. 86 (5): pp. 433–439. doi:10.1002/qua.10038  Zhou, Yaoqi; Karplus, Martin; Ball, Keith D.; Bery, R. Stephen (2002年). “The distance fluctuation criterion for melting: Comparison of square-well and Morse Potential models for clusters and homopolymers”. J. Chem. Phys. 116 (5): pp. 2323–2329. doi:10.1063/1.1426419  I.G. Kaplan, in Handbook of Molecular Physics and Quantum Chemistry, Wiley, 2003, p207. モース長距離ポテンシャル レナード-ジョーンズ・ポテンシャル 分子力学法