《结构化学基础》课程PPT教学课件（第三版）第三章 分子的结构和性质

①变分法原理 对于任意给定的一个标准(品优)波函数f,用体系的A算符求 得的能量近似值(期望值)ε(即能量平均值E),一定大于或 接近于体系基态的能量E,即: ∫f"r ≥E 0 f fat (f的平均能量E必是体系基态能量E的上限;若体系基态精确 波函数v0已知,则E=E0) 因此,可以任意选取一变分函数(试探函数),利用上式求 出能量的期望值,而且此值总是大于体系基态真实的能量。 能量的期望值越低,它就越接近体系基态真实的能量,相应 的试探函数也就越接近体系基态的真实波函数。 依据上式求体系近似解的方法称为变分法。 2021/1/21

2021/1/21 6 ① 变分法原理 因此，可以任意选取一变分函数（试探函数），利用上式求 出能量的期望值，而且此值总是大于体系基态真实的能量。 能量的期望值越低，它就越接近体系基态真实的能量，相应 的试探函数也就越接近体系基态的真实波函数。 依据上式求体系近似解的方法称为变分法。 E0 f fd f Hfd E =        ˆ 对于任意给定的一个标准（品优）波函数f，用体系的Ĥ算符求 得的能量近似值（期望值）ε（即能量平均值Ē），一定大于或 接近于体系基态的能量E0，即： （f 的平均能量Ē必是体系基态能量E0的上限；若体系基态精确 波函数0已知，则 E = E0 ）

②变分法解 Schrodinger方程的一般步骤 a:选择变分函数 常用的变分函数是选择已知标准函数的线性组合,即: f=c1φ1+c22+c33+……+n OE 然后求出e值最低时对应的c值,即 0 此时的E值已非常接近体系基态的能量E0,相应的也非 常接近体系基态的精确波函数v。如在H12中,选用两氢原子 a和b的基态波函数的线性组合作为H2+的变分函数,即: f=caat copE 式中的cn和C为待定参数。 由于分子轨道在一定程度上继承和反映原子轨道的规律,所 2021/1/21

2021/1/21 7 ② 变分法解Schrödinger方程的一般步骤 a:选择变分函数 常用的变分函数是选择已知标准函数的线性组合，即： f=c11+ c22+ c33+ ……+cnn = 0   i c E 然后求出ε值最低时对应的ci值，即 此时的E值已非常接近体系基态的能量E0，相应的f也非 常接近体系基态的精确波函数0。如在H2 +中，选用两氢原子 a和b的基态波函数的线性组合作为H2 +的变分函数，即： f=caa+ cbb 式中的ca和 cb为待定参数。 由于分子轨道在一定程度上继承和反映原子轨道的规律，所

以可用原子轨道的线性组合作为组成该分子的变分函数是合理 省时的。 用已知函数的线性组合作为变分函数的变分法称为线性变分法 -LCAo-mo( Linear Combination of Atomic Orbitals Molecular orbital)法。一般认为采用LCAO作为试探函数,有 可能是最佳的试探函数。 b.解久期行列式确定能量 将代入变分法原理公式中,得: 「(cnn+cb)H(cnh+c)dz E(Ca, Cb) ∫(cn+c)dz 由于H2中两个核是等同的,而中和中又都是归一化函数,展 开上式,并令: H ∫dz=Hb=」dr 2021/1/21 8

2021/1/21 8 H  H d H  H d a a a a b b b b ˆ ˆ     = = = 由于H2 +中两个核是等同的，而a和b又都是归一化函数，展 开上式，并令：         c c d c c H c c d E c c a a b b a a b b a a b b a b 2 ( ) ( ) ˆ ( ) ( , )   + + + = 将f代入变分法原理公式中，得： b. 解久期行列式确定能量 用已知函数的线性组合作为变分函数的变分法称为线性变分法 —LCAO-MO(Linear Combination of Atomic Orbitals￾Molecular Orbital)法。一般认为采用LCAO作为试探函数，有 可能是最佳的试探函数。 以可用原子轨道的线性组合作为组成该分子的变分函数是合理 省时的

b=如:dz=Hm=丁dr Sa=addr=Sb=%Bdr= ∫q4=S=∫中 得E(ca2Cb)==2 C2H +2c C Hab+Cb bb=Y/Z s+2c c,s bab +CiS bbb 根据变分原理,参数Ca,Cb的选择应使E最小,因此可令: OE OE =0即对Ca,C偏微商求极值,得 OC b OE aYY aZ OE 1 aY aZ 0 0 ac zz ac ac, Zac Z acb 2021/1/21

2021/1/21 9 H  H d H  H d a b a b b a b a ˆ ˆ     = = = = = = =1     S   d S   d aa a a bb b b =    = =        Sa b a b d Sb a b a d 根据变分原理，参数ca，cb的选择应使E最小，因此可令： a a a a b a b b b b a a a a b a b b b b a b c S c c S c S c H c c H c H E c c 2 2 2 2 2 2 + + + + 得 ( , ) = ： ＝Y/Z 0 1 2 =   −   =   a a a c Z Z Y c Y c Z E 0 1 2 =   −   =   b b b c Z Z Y c Y c Z E = 0   =   a b c E c E 即对Ca，Cb偏微商求极值，得：

消去Z,由Y=E得:OYp22=00)on0z=02) ac F ocb 将Y、团的表达式代入(1)得 a(ca H+.c,h+Cbhbb-E aca saa t2c cnsn+Cb Sbb=0 对上式微分得:2cnH2+2cbH2b-2cnSn2E-2cbS2bE=0 同理对(2)式有:2cbHb+2cH12b2 C SLE-2c2SaE=0 ap: ca(Haa-ESaa)+C,(Hab-ESab)=0 Ca(Hab-ESab+C (Hbb-ESbb)=0 2021/1/21 10

2021/1/21 10 消去Z，由Y/Z=E, 得： (1) (2) = 0   −   a a c Z E c Y = 0   −   b b c Z E c Y 0 ( 2 ) ( 2 ) 2 2 2 2 =   + + −   + + a a a a a b a b b b b a a a a a b a b b b b c c S c c S c S E c c H c c H c H 将Y、Z的表达式代入(1) 得： 对上式微分得：2caHaa+2cbHab-2caSaaE-2cbSabE=0 同理对(2)式有: 2cbHbb+2caHab-2cbSbbE-2caSabE=0 即：ca (Haa-ESaa)+ cb (Hab-ESab)=0 ca (Hab-ESab)+ cb (Hbb-ESbb)=0