分子電子構造論（水曜日 3 コマ目）資料１
評価方法：授業の出欠はとらない。評価は試験、レポートによる。ただし、片方だけでもよい。(具体的にはレポートの点数を A、試験の点数を B とすると評価 S は S = A + B - A*B/100)で計算する。もし、他に評価項目 C がある場合には S = A + B + C - A*B/100 - B*C/100 - C*A/100 +A*B*C/1002 以下略）
試験：予想問題をだしますのでそれを参考に。
予備知識：線形代数、簡単な微積分、物理化学程度の量子論
進め方：プリント or 論文（英文 or 日本語）を読みながら理解していきます。適宜、当てます。また簡単な発表をしてもらうこともあります。予定は随時変わります。

H30 年度の講義予定
４／11 ガイダンス
４／18 １回目 ガイダンス＋数学的準備 1
４／25 ２回目 HF 理論＋α
5／09* ３回目 構造最適化、振動解析
5／16* ４回目 制限、非制限の HF 理論
5／23 ５回目 postHF（その１）
5／30 ６回目 波動関数の解説 （基底関数を中心に）
6／6 ７回目 postHF(その２)
6／13 8 回目 その他の話題 １
6／20 予備日(なにもなければレポートを宿題としてだします。)
6／27 9 回目

その他の話題 ２

6／29 10 回目

その他の話題 3

7／4 11 回目 DFT の種類、開発など
７／11 12 回目 計算相談
７／18 13 回目 試験

1

0.数学的準備
線形代数、微積、量子力学の簡単な復習をおこなう。
0-1. 基底と線形独立性
与えられたベクトル空間に対して任意のベクトルｘはｒ個の一次独立なベクトルの集合｛a1,a2,…,ar｝に対してそれらの一次結合として
x = c1a1 + c2a2 + … + crarと表され、その表現が一意的である。すなわち係数 ci (1 ≤ i ≤ r)は x によって一意的に定まる
このとき、ベクトルの集合｛a1,a2,…,ar｝を基底という。またこの数を次元という。
例、カルテシアンの座標系では x,y,z
ベクトル以外の関数にも拡張できて、これらも基底と呼ばれるテーラー展開における xm ただし m は自然数
フーリエ級数における sin mx, cos nx ただし m, n は自然数
0-2. 正規直交系
互いに直交してその内積が０であり、かつその大きさ(*ノルム)が１であるベクトルの集まり＜ai, aj＞= δij (i = j のとき１, i ≠ j のとき 0)*ノルム=ベクトルの長さを一般化した概念で距離を定義するために用いられる。距離空間を形成するための条件を満たす。たとえば 2 つのノルム||a||, ||b||は次の式を満たす。
||a|| ≥ 0,
|α|||a|| = ||αa||,||a + b|| ≤ ||a|| + ||b||

良く使われる例
距離２つの関数 f(x), g(x)の差、 max | f(x) - g(x) |0-3. 空間の完全性
任意のベクトル（関数も含む）を基底をつかって線形結合で表せるとき完全系をなす。
つまり任意のｘがつねに x = c1a1 + c2a2 + … + crar, で表せる。

2関数の場合の完全系のとき任意関数がある直交関数（基底が関数に相当する場合の直交系）で展開できるとき完全系という。
例、球面調和関数、ラゲール多項式、Hermite 多項式、フーリエ級数
0-4. Euler-Lagrange 方程式
ラグラジアン L は運動エネルギーmv2/2 とポテンシャル V(x)の差の形をしておりと書ける。
L を極小にするように物体が運動することが知られており、運動の軌跡は変分法で求めることができる。そのとき満たすべき条件式は作用が０になることでありこれを Euler-Lagrange 方程式という。L の式に代入するととなりの運動方程式に一致する。

変分法では Euler 方程式と呼ばれ汎関数に対してで定義される。この式が満たされるとき汎関数が極小値をとる。

（Euler 方程式は解くべき微分方程式を与えるもので解そのものではないことに注意）
エネルギーE が波動関数によって極小化される場合にも適用され、HF 法や DFT の定式であるロータン方程式やコーンシャム方程式の導出にも用いられる。
0-5. Lagrange 未定定数法
変分法において拘束条件に極値を与えるには

のもとで、

が

の極値曲線でない場合、3の極値曲線になるようなλが存在する。つまり、Euler 方程式はとなる。

このときのλを未定乗数と呼ぶ。波動関数は積分すると電子の個数 N となることからそのような拘束条件をいれて上記の最適化を実行することでエネルギーの定式を得る。

0-6. 対角化、行列式
一般的に正方行列に適当な変換してもとの行列と相似な対角行列を得るような変形をいう。
Hermite 行列という転置自己共役にたいして不変な行列はユニタリ行列によって対角化（＝対角成分以外０に変形）可能である。
ｎ個の線形方程式で表される場合、基底が直交しかつ演算子に対してスカラー倍として表現できる場合を指す。良く用いられるのは Hermite 演算子で任意の基底 a, b に対して内積＜ax, b＞ = ＜a, xb＞をみたす演算子を Hermite 演算子（自己共役演算子）という。量子化学でとりあつかうエネ
ルギーに対応する Hamilton 演算子は Hermite 演算子であり、固有値が実数でかつ、ユニタリ演算子によって対角可能である。

0-7. 固有値、固有ベクトル関数（ベクトル）が演算子（行列）に対してスカラー倍（固有値）として表現できる場合、固有関数（固有ベクトル）という。量子力学では主に波動関数に対する固有値問題を対象としてエネ
ルギーを波動関数の Hamilton 演算子に対する固有値として求める。通常の方程式は解を求めるのに対して、固有値問題は固有関数である波動関数と固有値であるエネルギーを同時に求めるという点が特徴である。また、量子力学において固有値がとびとびの値をとる量子化という現象が見られることも重要である。そのような値に対応する固有関数は固有状態とよばれる状態に対応する。これらの系は無限個の基底からなるが実際には有限基底の中で近似解を求めるためいろいろな問題点に直面するがそれについては後で詳細に説明する。

4

1.基礎方程式
量子力学の中で古典的な軌跡の概念に取って代わる波のことを波動関数という。1926 年、オーストリアの物理学者 Schrödinger は任意の系の波動関数を求めるための方程式を提案した。この方程式を Schrödinger 方程式という。Schrödinger 方程式は量子力学の基本的方程式であり、波動関数がどんな系をも完全に記述できる。今回は、時間を変数として含まない型の Schrödinger 方程式を提示し説明する。時間に依存しない Schrödinger 方程式は時間とは無関係なので、その解は定常状態の波動関数という。化学者にとって興味ある原子、分子の電子状態などの問題が、定常状態の波動関数のみを用いて扱える。

1-1.Schrödinger 方程式
時間を含まない Schrödinger 方程式は、(1.1.1)で与える。E は電子のエネルギー、は波動関数、る。

からなっている。N 個の電子と M 個のは電子 i の座標であり、空間座標 とスピン座標

核に対するは Hamilton 演算子であは、(1.1.2)である。Mαは核の質量と電子の質量比であり、ラシアン演算子とはおのおの i 番目の電子と

および

は核

と核

の電荷である。ラプ

番目の核の座標に関する微分を含む。第１項は、電子の運動エネルギーに対応する演算子、第２項は核の運動エネルギーに対応する演算子
で、第３項が電子と核の間のクーロン引力を、第４と第５項はそれぞれ電子間、核間の斥力をあらわしている。
次に、Born-Oppenheimer 近似を導入する。この近似は、核は電子と比べると非常に重いので、電子が核に対して相対的に動いているとき核は静止している、という考え方である。この近似のもとでは、式 (1.1.2) の第２項の核の運動エネルギーは無視できて、最後の項の核間反発は定数とみなせる。これより(1.1.3)5

となる。原子核間の反発エネルギーは(1.1.4)である。全エネルギーW は電子エネルギーE と原子核間の反発エネルギーの和で(1.1.5)である。式 (1.1.1) を E について解いてからいてを加えても、あるいは

の定義に

を含めてお

の形の Schrödinger 方程式を解いても、どちらでもかまわない。

1-2.Born の解釈
波動関数は、それによって記述される系に関するあらゆる力学的な情報を含んでいる、というのが量子力学の主張の中心にある。ここで粒子の位置に関して波動関数がもっている情報に注目してみる。粒子の位置に関しての波動関数の解釈として、Born の解釈がある。この解釈は、波動関数の２乗 (または、もしが複素関数なら絶対値の２乗、) の値が、その粒子を見出す確率に比例する、というものである。具体的にいうと、「ある点 r において、ある粒子の波動関数がという値をもつなら、r と r＋dr の間にその粒子を見出す確率はに比例する」ということである。

1-3.波動関数の制限

Schrödinger 方程式は固有値方程式であるから、一般に意の定数とするとき N

がその固有関数の解であれば、N を任もその方程式の解である。これより、Born の解釈の比例を等号にかえる

ような、規格化定数 N をいつでも見つけることができることになる。また、度として解釈されるので、は粒子の確率密を全空間で積分したとき、その解は１でなければならない。これらの要請から、次の等式が成り立たなければならない。
(1.3.1)
さらに、波動関数には Pauli の原理を満たすことが条件づけられている。これはつまり、フェルミ粒子である二つの電子は座標の交換に対して反対称でなければならないというものである。
(1.3.2)6

1-4.スピン軌道と空間軌道１個の粒子、つまり１電子に対する波動関数を軌道と定義する。分子の中の１電子の波動関数に対しては分子軌道を使っていくことにする。空間軌道
電子の空間分布を記述しており、は位置ベクトル r の関数で、そのが r にある微小体積要素 dr の中にその電子を見出す

確率を与える。空間分子軌道は、普通は規格直交系をなすと仮定されている。
(1.4.1)
ここで、 は Kronecker のデルタ記号であり、が完全系だとすれば、任意の関数

では１、では０となる。空間軌道の組は正確に次のように展開できる。

(1.4.2)
ここで、 は定数の係数である。完全系であるためには無限個からなる組でなければならないが、実際には無限個の集合は扱えないので K 個の軌道からなる有限個の組のみを考える。この有限の組は、完全系の張る空間の中のある部分空間を張るにすぎないが、得られた結果はこの有限個の軌道で得られる部分空間の中での“正確”な結果といういい方ができる。
１個の電子を完全に記述するためには、スピンの特定が不可欠である。電子のスピンを記述する完全系は、２つの規格直交関数と、すなわちスピン上向き（↑）とスピン下向き（↓）

からなっている。１個の電子の空間分布とスピンをともに記述する波動関数が、スピン軌道である。x は空間とスピン、両方の座標をあらわしている。各空間軌道、から２つのスピン軌道をつくることができる。一方はスピン上向き、もう一方はスピン下向きに対応していて、おのおのその空間軌道にαまたはβのスピン関数をかけたものである。

(1.4.3)

K 個の空間軌道の組が与えられれば、２K 個のスピン軌道の組をつくることができる。

(1.4.4)

空間軌道が規格直交ならば、スピン軌道も規格直交である。

7

(1.4.5)

2.Hartree-Fock 理論
電子の Schrödinger 方程式に対する近似解を見つけることは、量子力学の誕生以来、量子化学者の第一の目的であり、我々は常に多電子問題を取り扱わねばならない。この多電子問題を解く試
みの中心となるのが Hartree-Fock 法である。この方法では、１個の配置関数で波動関数を近似し、変分法によって最良のスピン軌道の組を求めることができる。

2-1.Slater 行列式を近似的に表現するのに、ここでは規格直交化された N 個のスピン軌道の反対称積を用いる。この反対称積が Slater 行列式(2.1.1)である。変分原理法によると、電子座標に関するいかなる反対称規格化波動関数に対するエネルギーの期待値

は、常に正確な波動関数の、基底状態のエネルギーE0 より高くなる。
(2.1.2)
等号はのときのみ成り立つ。ここで規格化された Slater 行列式に対する、簡略で便利な記法を導入する。この記法では、規格化定数をすでに含んでいるとし、その行列式の対角成分のみを書く。
(2.1.3)もし電子の座標がつねにの順に並べてあるとすると、さらに次のように略記できる。
(2.1.4)

この記法を用いると、式 (2.1.1) は次のように書ける。

(2.1.5)8

また、は通常ブラケット表示を用い次のように書かれる。
(2.1.6)2-2.Hartree-Fock 方程式
N 電子系の基底状態を記述するのに使うことのできる最も単純な反対称波動関数は１個の Slater
行列式である。
(2.2.1)変分原理によれば、この汎関数形をもつ最良の波動関数は、可能な限り低いエネルギー(2.2.2)を極小にするようなものである。ここで(2.2.3)(2.2.4)(2.2.5)である。これらの積分はすべて実数値をとり、である。

はクーロン積分、は交換積分と呼ばれる。次の等式(2.2.6)は重要である。式 (2.2.2) の 2 重和にの項を含めてある理由は、この等式にある。式 (2.2.2) を規格直交化の条件のもとで最小にすると、Hartree-Fock の微分方程式(2.2.7)(2.2.8)が得られる。ここで、はスピン軌道とを占めたときの軌道エネルギー、の電子が存在することによって電子 i が感じる平均ポテンシャルであり、9は他は Fock 演算子と呼ばれる。また

、

ある。

は

と波動関数

より、

は Hermite 演算子、は Hermite 行列 の要素での間の任意の Unitary 変換に不変である。また、

Hermite 行列は適当な Unitary 行列で対角化することができる。この任意性を利用すると、式 (2.2.7)は次のように書き換えることができる。
(2.2.9)こうして多電子問題が１電子問題に帰着されたわけである。この方程式の解は反復計算により求めなければならない。それは解として求めるべきが演算子の中に現れるからである。

2-3.Hartree-Fock 法の利点と問題点
Hartree-Fock 法の利点は、電子構造をみながら相互作用を数え上げることにより、エネルギーの表現が得られることである。これは閉殻構造ばかりでなく、開殻構造の電子状態であっても、波動関数が単一行列式で表現されていれば、いつでも相互作用を数えあげられることを意味する。
Hartree-Fock 法は、多くの系について初期的な予測を提供してくれるので非常に便利である。また、安定分子や遷移状態のいくつかの構造や振動数についても妥当な計算結果を与える。つまり、基本的なレベルの良い方法ということができる。
Hartree-Fock 法では、二個の電子を交換してもは変化しないという要請を、波動関数を反対称化することで満たしている。これは、同じスピンを持った一対の電子から来る主な相関を自動的に含む。この相関を交換相関という。しかし、反対スピンを持った電子の運動は相関のないま
まとなっている。このため、反応や結合解離のエネルギーを正確に調べるのには不十分である。
そのため、我々は次に述べる電子相関を含む密度汎関数法を用いる。

3.密度汎関数理論
密度汎関数理論(DFT)とは、波動関数のかわりに電子密度ρを変数として用いて、電子密度のみで電子状態のエネルギーが表せるという点である。つまり、電子密度の空間分布さえ知ることができれば、基底状態の電子状態を記述する理論である。
N 電子系では 3N 次元空間の Schrödinger
方程式を解かなくても、三次元空間における電子密度と、これに関連した相関交換汎関数か

ら、基底状態のエネルギーが求められる。そのため、電子密度が定まれば、基底状態に対応

する反対称化された波動関数が得られるはずである。このとき、電子密度

10は N 表示可能であるという。また、基底状態のエネルギーと、基底状態に対応した反対称化しされた波動関数の二つが定まれば、固有値問題から基底状態の Hamilton 演算子も完全に決定される。
そのため、
Hamilton
演算子の外部ポテンシャルも決定される。このとき、電子密度

は v 表示可能であるという。N 表示可能の条件は、が非負で連続、かつ全空間を積分すると N になることである。しかし、v 表示可能の成立条件は、いくつかの必要条件しか明らかでない。電子密度表示可能であるとき密度汎関数理論は意味を持つ。逆に電子密度が N、vが N 表示可能でないときには、最小解に対応した波動関数から得られた電子密度ではないし、v 表示可能でないときには、エネルギーを密度の汎関数として最小化させることが、基底状態のエネルギーを求めることに一致しない。

3-1.Hohenberg-Kohn の定理
オリジナルの密度汎関数理論はHohenberg-Kohn の定理に基づいている。まず、大きな箱中の時間に依存しない局所外部ポテンシャルの影響下で運動するN 電子系を考える（ここでは縮退していない基底状態のみを考慮する。
）。Hohenberg-Kohnの第一の定理は、電子密度によってポテンシャルが決定されることを示している（ただし、重要でない定数が付随する。）。つまり、電子密度と波動関数には一対一の関係があるということである。このことは以下で背理法によって証明する。同じ密度およびで関係づけられる別のポテンシャルが存在すると仮定する。この仮定により、二つの外部ポテンシャルに対応する二つの異なる基底状態波動関数
ネルギー

および

および

を持つ異なるHamilton演算子、

および

とそれ

が存在する。言い換えると、基底状態エとが存在するということも表している。Hamilton演算子は(3.1.1)で与えられる。Rayleigh-Ritzの変分原理により(3.1.2)を得るが、ととを入れ替えることにより、

11

(3.1.3)を得る。二つの式を足すとにより符号が矛盾することが分かる。従って密度から一意に外部ポテンシャルが求められ、さらには、系のHamilton演算子や電子特性を決めることが可能であると結論づけられる。

を電子の運動エネルギー項、

ルエネルギー、を核-電子引力と核間の反発のポテンシャ

を電子-電子反発項としたとき、全エネルギーはは(3.1.4)で与えられる。ここで(3.1.5)である。

は全運動エネルギーと電子間反発エネルギーの和である。

一方、Hohenberg-Kohnの第二の定理は、いかなる試行密度に対して(3.1.6)が成り立つことを示している。ただし、、である。証明は試行波動関数のエネルギー(3.1.7)に対する変分原理から示される。ここで、等号は真の基底状態の時を表す。さらに電子数を一定にする拘束条件を課すことにより、変分(3.1.8)から、Euler方程式(3.1.9)を得る。ここでLagrange未定乗数

は化学ポテンシャル(あるいは負の電気陰性度) を表す。もしの厳密な形を知っていれば、式 (3.1.7) は基底状態の電子密度についての厳密な方程式となる。ただしは外部ポテンシャルない。これはが、の普遍な汎関数であることを意味する。

Hohenberg-Kohnの定理は、電子密度
しているが、電子密度とは独立に定義されることに注意しなければならから外部ポテンシャルが一意的に定まることを示

から基底状態の外部ポテンシャルを与えることまでは保証していない。

12

これを補う役割を果たしているのがv表示可能性である。そこで、Hohenberg-Kohnの第一の定理を言い換える。すなわち、基底状態の波動関数と、v表示可能である電子密度とは一対一の対応がある。この対応関係によってv表示可能な密度から、これに付随する基底状態の性質が決まる。したがって、基底状態の全ての性質がその電子密度の汎関数であるというときには、これらの汎関数が全てv表示可能な密度についてのみ定義されていることを理解する必要がある。また、汎関数および変分法に現れる密度に対する条件を弱めても、いわゆるN表示可能性を要求すれば、密度汎関
数理論の定式化が可能であることが知られている。

3-2.Kohn-Sham 方程式
Kohn と Sham は、運動エネルギー汎関数

に対する巧妙な間接的アプローチを創出した。

これにより密度汎関数理論は厳密な計算を行うための道具となった。Kohn と Sham が提案した内容は、運動エネルギーを精度よく、しかも単純な方法で計算できるように軌道を取り入れることにし、このとき生じる小さな補正を別に取り扱うことであった。Kohn と Sham は、Hohenberg とKohn の普遍的な汎関数の定義（式 (3.1.5) ）にならい、これに対応する形で相互作用のない参照系を導入した。その Hamilton 演算子は

(3.2.1)である。この系には電子-電子反発項が含まれておらず、基底状態の電子密度は近似なしにとなっている。この系では、基底状態の波動関数は近似なしに(3.2.2)である。ここでは１電子 Hamilton 演算子の固有状態で、エネルギーの低いほうから N 個を採用する。すなわち(3.2.3)である。運動エネルギーは(3.2.4)13

である。ここで式 (3.1.5) を書き換えて(3.2.5)とする。ただし(3.2.6)である。ここで定義された量は交換-相関エネルギーと呼ばれる。そこに含まれるものは、おそらく十分に小さいと考えられるとの差、および非古典的な部分のである。

こうすると式 (3.1.9) の Euler 方程式は(3.2.7)となる。ここで Kohn-Sham 有効ポテンシャル

は(3.2.8)で定義される。さらに(3.2.9)は交換-相関ポテンシャルである。
制限条件のもとで、外部ポテンシャルの中で互いに相互作用なく運動する電子の系に、通常の密度汎関数法を適用したときに得られる式と、式 (3.2.7) は完全に一致する。従って、与えられたに対して、単に N 個の１電子方程式(3.2.10)を解いて(3.2.11)とすれば、式 (3.2.7) を満たすを得る。この N 個の方程式が Kohn-Sham 方程式と呼ばれる。

Korn-Sham 方程式は、相互作用のある電子の運動を独立粒子系として扱うことができることを示14

している。つまり、電子は共通の局所ポテンシャルの中で運動している１個の電子として記述される。非常に複雑な多電子の関与する電子間相互作用は全て局所ポテンシャル

へ押し込めている。
また、Korn と Sham の考え方の狙いは、真の運動エネルギーの中心的な部分である

どのように決めれば、Kohn と Sham の補助的軌道を用い、明示的な形で、をを書けるようにできるかというところにある。しかし、縮退していない基底状態に対する相互作用のない、v 表示
可能な密度については、(3.3.1)
が成り立つ。これより、式 (3.2.6) で定義した交換-相関エネルギーには、正の運動エネルギー分が含まれている。

3-3.密度汎関数理論利点の問題点
Kohn と Sham の理論では、一様な電子ガスモデルに基づく局所交換相関ポテンシャルを用いているので、局所密度近似(LDA)と呼ばれる。このモデルでは、実際には一様でない電子分布を微小部分にわけ、個々の微小部分 dr に関して一様な電子ガスモデルの結果を適用し、それらを足し合わせることで、電子分布からエネルギーを求めている。そのため個々の微小部分 dr の電子密度が分かればよいという局所的性質を持っている。
しかし、不対電子を有する系に LDA を当てはめて系を記述すると、半整数の電子を縮退した軌道に詰めねばならず、系を適切に記述できない。そのため、全電子密度
電子スピン密度だけでなく、新たにを導入し、開殻系でもより正確に記述できるようになった。このスピン分極を考慮した方法は、全電子密度だけの場合と区別してスピン密度汎関数と呼ばれ、LDA に対してスピン分極がある場合は局所スピン密度近似(LSDA)と表記される。スピン分極を持ち込むこと
で、密度汎関数理論による分子の計算精度は格段に上がる。しかし、電子が自分自身と相互作用してしまうこと(自己相互作用)や、交換エネルギーの誤差が大きいことなどが欠点として依然として残っている。
また、密度汎関数理論は、相互作用のない v 表示可能な密度についてだけ成り立つという意味
において、制限されたものとなっている。すなわち、きちんと定義された相互作用のない系の基底状態密度にだけしか成り立たない。

15

さらに密度汎関数理論ではしばしば Korn-Sham 方程式を解く際に、HOMO-LUMO ギャップが小さく、SCF の収束が遅いことが問題となっている。そのために安定な波動関数を SCF の初期状態に与えることが重要になる。もし最初に不適切な波動関数が与えられてしまうと、SCF の収束が困難になってしまうか、基底状態ではない状態に収束してしまうおそれがある。

3-4.hybrid 法
Hartree-Fock (HF) 法は同一電子のクーロン反発項と交換項は完全にキャンセルする性質を持つ。
このように電子が自分自身と相互作用しないことを“正確な交換”という。しかし密度汎関数法の交換エネルギーは、前節で述べたように自己相互作用をもつために精度が悪くなっている。特に交換エネルギーは相関エネルギーに比べて絶対値が大きいので、交換エネルギーの誤差は全エ
ネルギーの精度に影響する。そのため HF 法の交換エネルギーの厳密な取り扱いを密度汎関数法に取り込み、自己相互作用をなくそうという試みが Becke により提案された。今までの密度汎関数法を純粋な密度汎関数というのに対して、Becke により提案された新しい汎関数は hybrid 汎関数と呼ばれる。ここでは hybrid 汎関数の一つである、B3LYP 法(Becke Three Parameter Hybrid
Functionals) を紹介する。

(3.3.1)
ここで、は交換-相関エネルギー、の交換エネルギー、
エネルギー、は HF の交換エネルギー、は Slaterは Becke の非局所的な部分の交換エネルギー、の相関の相関エネルギーである。

パラメータ A は HF と LSDA の交換エネルギーの混合比を決定する。汎関数を局所項と非局所項に分割したときの、非局所的な部分の交換エネルギー、相関エネルギーの割合がパラメータ B とC である。そしていろいろな分子や原子に DFT 計算を実行し、汎関数の混合比を最適化した。計算の対称は G1 データベースと呼ばれる 56 分子の原子化エネルギー、42 分子のイオン化エネルギ
ー、8 分子のプロトン親和力、及び周期表第一周期元素の原子エネルギーである。最終的に得られた A = 0.20、B = 0.72、C = 0.81 は実験値に合致するように調整された半経験的に係数である。その結果、G1 データベース 116 個の実験値からのずれを根二乗平均したものは 3.35 kcal/mol となり満足のいくものとなった。また、LYP 汎関数は局所項と非局所項に分けることができないので、LYP

16

汎関数の局所部分と VWN 汎関数を足して 1.00 になるように工夫されている。

17