分子力学法のソースを表示

[[ファイル:Bond stretching energy.svg|thumb|right|[[力場 (化学)|力場]]がこのエタン分子の結合伸縮エネルギーを最小化するために使われる。]] 
'''分子力学法'''（ぶんしりきがくほう、Molecular Mechanicsの頭文字より'''MM法'''と略される）あるいは'''分子力場計算'''（ぶんしりきばけいさん）は、分子の[[立体配座]]の安定性や配座間のエネルギー差を[[原子]]間に働く力による[[ポテンシャルエネルギー]]の総和によって計算する手法のことである。

[[分子]]の持つエネルギーは[[シュレーディンガー方程式]]を解くことによって計算することが可能であるが、これは分子を構成する原子および[[電子]]の数が多くなると計算量が急激に増加し困難になる。

しかしその一方で、分子の内部の原子同士に働く力はその原子の種類や結合様式が同じならば、別の種類の分子でもほぼ同じである。例えば[[混成軌道|sp<sup>3</sup>混成]]の[[炭素]]原子と[[水素]]原子の結合距離はどのような分子でもほぼ0.11 nm、結合エネルギーはほぼ4.1×10<sup>2</sup> kJ mol<sup>−1</sup>、[[赤外分光法|赤外吸収スペクトル]]でほぼ2950 cm<sup>−1</sup>付近に吸収を示す。このことはsp<sup>3</sup>混成の炭素原子と水素原子の結合距離の伸縮に伴って2つの原子間に働く力が分子によらず、ある一つの数式で表すことができることを示唆している。

そこで原子間に働くすべての力を、原子間の結合を表すパラメータ（結合距離、[[結合角]]など）を変数とし、原子の種類や結合様式によって決まる関数で表す。そしてそれらの力によるポテンシャルエネルギーの総和が分子の持つエネルギーとなっていると考える。このような考えのもとに、様々な実験値をうまく説明できるような[[原子間ポテンシャル|原子間のポテンシャルエネルギー]]を表す式を経験的に、あるいは量子化学的手法によって導き、それによって分子の[[立体配座]]の安定性や配座間のエネルギー差を計算する手法が分子力学法である。このポテンシャルエネルギーによる力の[[場]]を'''[[力場 (化学)|分子力場]]'''というため、分子力場計算ともいう。

分子力学法での分子のイメージは原子を球として、その球をその両端の原子の種類によって強さが決まった[[ばね]]で結びつけたような感じである。それぞれのばねが引きのばされたり押し縮められたりして、すべてのばねのポテンシャルエネルギーの合計が最も小さくなったところで、その分子の持つエネルギーと立体配座が決まることになる。

分子力学法は計算量が[[量子化学]]的手法に比べて少ないため、原子数の多い分子においても容易に計算結果が得られるという利点がある。しかし、原子の種類や結合様式が異なるごとに別のポテンシャルエネルギーの式を使わなければならないため、前もって準備しておくパラメータの数が非常に多くなってしまうという欠点がある。

単純な分子力学法においては分子の持つポテンシャルエネルギーEは
:<math>E = \sum E_d + \sum E_a + \sum E_t + \sum E_n</math>
*E<sub>d</sub>は直接結合している2原子間の結合距離によるポテンシャルエネルギー
*E<sub>a</sub>は直接結合している3原子間の結合角によるポテンシャルエネルギー
*E<sub>t</sub>は直接結合している4原子間の二面角によるポテンシャルエネルギー
*E<sub>n</sub>は直接結合していない2原子が空間的に接近することによるポテンシャルエネルギー
と表される。

E<sub>d</sub>は結合の伸縮によって生じるポテンシャルエネルギーを表す。E<sub>d</sub>の経験的で良い近似とされる式は[[モース・ポテンシャル]]が良く知られているが、この式は指数関数を含んでいるため計算量が多くなり、分子力学法のメリットを打ち消してしまうためあまり用いられない。また古典的な方法としては[[フックの法則]]に従う[[調和振動子]]として取り扱う方法がある。近似の程度はあまり高くないが、計算量が少ないため高分子の計算に使用される。より近似の程度を改善したものとしてはモース・ポテンシャルの平衡結合距離付近での振る舞いを近似した結合距離の3次式、あるいは4次式などが使用される。しかしこれらの式では、平衡結合距離から離れた部分ではエネルギー値がまったく当てにならなくなってしまうため、反応の[[遷移状態]]のエネルギー計算などは不可能である。

E<sub>a</sub>は結合角の圧縮や拡大によって生じるポテンシャルエネルギーを表す。E<sub>a</sub>を表す経験的な式には、調和振動子として取り扱う方法がまず挙げられる。さらに近似を高める場合には結合角の4次式や6次式が使用される。

E<sub>t</sub>は単結合の回りの回転に伴うポテンシャルエネルギーや、二重結合の平面から原子がずれることによって生じるポテンシャルエネルギーを表している。E<sub>t</sub>は二面角を360度回転させると元の形に戻る周期関数であるので、[[フーリエ級数]]型の式が用いられる。近似を高める場合には高次のフーリエ級数の項を含む式が使用される。

E<sub>n</sub>は立体的な反発や[[非共有電子対]]による[[クーロン力|静電的な斥力]]や[[ファン・デル・ワールス力]]や[[水素結合]]による引力が含まれるので様々な形の項を含む。立体的な反発とファン・デル・ワールス力を良く表すものとしては、引力が原子間距離の6乗に反比例し斥力が原子間距離の12乗に反比例する[[レナード-ジョーンズ・ポテンシャル]]があり、これがよく使用されている。

より精度を追求した分子力場では、例えば結合距離が伸びたときには二面角によるエネルギーを補正するというような、2種類以上の構造を表す変数を含む交差項とよばれる項が導入されている。

代表的な分子力場の種類としては以下のようなものが挙げられる。それぞれターゲットとしている用途や分子に特色がある。
*MM2、MM3、MM4 (Molecular Mechanics program 2 or 3 or 4)
*[[MMFF]]94 (Merck Molecular Force Field 94)
*[[AMBER (分子動力学)|AMBER]] (Assisted Model Building and Energy Refinement)
*[[CHARMM]] (Chemistry at HARvard Macromolecular Mechanics)

分子力学法においては分子力場から受ける力に従って原子の空間的な配置を移動させていくことで、よりエネルギーの低い配座に変化させていき、最終的には停留点として安定配座を求めることができる。単純なソフトウェアでは初期に入力した立体配座に近い安定配座の1つが求まるのみであるが、より高機能なソフトウェアでは複数の初期配座を自動的に発生させて、複数の安定配座を求められるものもある。

通常の分子力学法は静止した分子を扱うのに対し、[[分子動力学法]]は分子やその集団の運動を扱う。分子動力学法の中でも分子内あるいは分子間に働く力を記述するために分子力学法が利用されている。

{{sci-stub}}

{{DEFAULTSORT:ふんしりきかくほう}}
[[Category:計算化学]]
[[Category:分子モデリング]]
[[Category:分子物理学]]