量子化学计算方法试验 1. 应用量子化学计算方法进行计算的意义 化学是一门基础学科，具有坚实的理论基础，化学已经发展为实验和理论并重的科学。理论化学和实验化学的主要区别在于，实验化学要求把各种具体的化学物质放在一起做试验，看会产生什么新的物质，而理论化学则是通过物理学的规律来预测、计算它可能产生的结果，这种计算和预测主要借助计算机的模拟。也就是说，理论化学可以更深刻地揭示实验结果的本质并阐述规律，还可以对物质的结构和性能预测从而促进科学的发展。特别是近几年来，随着分子电子结构、动力学理论研究的不断深入以及计算机的飞速发展，理论与计算化学已经发展成为化学、生物化学及相关领域中不可缺少的重要方向。目前，已有多种成熟的计算化学程序和商业软件可以方便地用于定量研究分子的各种物理化学性质，是对化学实验的重要的补充，不仅如此，理论计算与模拟还是药物、功能材料研发环境科学的领域的重要实用工具。 理论化学运用非实验的推算来解释或预测化合物的各种现象。理论化学主要包括量子化学，（quantum chemistry）是应用量子力学的基本原理和方法研究化学问题的一门基础科学。研究范围包括稳定和不稳定分子的结构、性能及其结构与性能之间的关系；分子与分子之间的相互作用；分子与分子之间的相互碰撞和相互反应等问题。量子化学可分基础研究和应用研究两大类，基础研究主要是寻求量子化学中的自身规律，建立量子化学的多体方法和计算方法等，多体方法包括化学键理论、密度矩阵理论和传播子理论，以及多级微扰理论、群论和图论在量子化学中的应用等。理论与计算化学的巨大进展，正使化学学科经历着革命性的变化。今天的理论与计算化学几乎渗透到现代一切科技领域，与材料、生物、能源、信息和环保尤为密切，理论化学的应用范围将越来越广。理论与计算化学逐步发展成为一门实用、高效、富有创造性的基础科学，在化学、生物学等领域的影响越来越显著，且与日剧增。 2. 应用量子化学计算方法进行计算的目的 （1）了解量子化学计算的用途。 （2）了解量子化学计算的原理、方法和步骤。 （3）通过一两个计算实例进行量子化学计算的上机操作试验。 （4）学会简单的分析和应用计算结果。 3. 量子化学计算试验的原理
第二章 一维势场中的粒子 §2.2 方 势 一、一维运动 当粒子在势场V (x ,y ,z )中运动时，其 Schrodinger 方程为： 22 [(,,)](,,)(,,)2V x y z x y z E x y z m ψψ-?+= 若势可写成： V (x ,y ,z ) = V 1(x ) + V 2(y ) + V 3(z ) 形式， 2212 [()]()()2x d V x X x E X x m dx -+= 2222 [()]()()2y d V y Y y E Y y m dy -+= 2232 [()]()()2z d V z Z z E Z z m dz -+= ψ(x ,y ,z ) = X (x ) Y (y ) Z (z ) ψ1(x ) x y z E E E E =++ 二、一维无限深势阱 0(0)()(0,) x a V x x x a ?<?? 这是定态问题 一维无限深势阱（0~a ）的求解 解：（1）列出各势域的 S — 方程 22 2 [()]()()2d V x x E x m dx ψψ-+= 20222 2 2202 22()0202()0I I II II III III d m V E dx d mE dx d m V E dx ψψψψψψ?--=???+=???--=?? 00E V << 0()V →∞ ，令k = )(0>k ，β=方程可简化为：22 2 222 222 000I I II II III III d dx d k dx d dx ψβψψψψβψ?-=????+=???-=??
清华大学大学物理习题库：量子物理 一、选择题 1．4185：已知一单色光照射在钠表面上，测得光电子的最大动能是1.2 eV ，而钠的红限波长是5400 ?，那么入射光的波长是 (A) 5350 ? (B) 5000 ? (C) 4350 ? (D) 3550 ? ［ ］ 2．4244：在均匀磁场B 内放置一极薄的金属片，其红限波长为??。今用单色光照射，发现有电子放出，有些放出的电子(质量为m ，电荷的绝对值为e )在垂直于磁场的平面内作半径为R 的圆周运动，那末此照射光光子的能量是： (A) 0λhc (B) 0λhc m eRB 2)(2+ (C) 0λhc m eRB + (D) 0λhc eRB 2+ ［ ］ 3．4383：用频率为??的单色光照射某种金属时，逸出光电子的最大动能为E K ；若改用 频率为2??的单色光照射此种金属时，则逸出光电子的最大动能为： (A) 2 E K (B) 2h ??- E K (C) h ??- E K (D) h ??+ E K ［ ］ 4．4737： 在康普顿效应实验中，若散射光波长是入射光波长的1.2倍，则散射光光子能量?与反冲电子动能E K 之比??/ E K 为 (A) 2 (B) 3 (C) 4 (D) 5 ［ ］ 5．4190：要使处于基态的氢原子受激发后能发射赖曼系(由激发态跃迁到基态发射的各谱线组成的谱线系)的最长波长的谱线，至少应向基态氢原子提供的能量是 (A) 1.5 eV (B) 3.4 eV (C) 10.2 eV (D) 13.6 eV ［ ］ 6．4197：由氢原子理论知，当大量氢原子处于n =3的激发态时，原子跃迁将发出： (A) 一种波长的光 (B) 两种波长的光 (C) 三种波长的光 (D) 连续光谱 ［ ］ 7．4748：已知氢原子从基态激发到某一定态所需能量为10.19 eV ，当氢原子从能量为－0.85 eV 的状态跃迁到上述定态时，所发射的光子的能量为 (A) 2.56 eV (B) 3.41 eV (C) 4.25 eV (D) 9.95 eV ［ ］ 8．4750：在气体放电管中，用能量为12.1 eV 的电子去轰击处于基态的氢原子，此时氢原子所能发射的光子的能量只能是 (A) 12.1 eV (B) 10.2 eV (C) 12.1 eV ，10.2 eV 和 1.9 eV (D) 12.1 eV ，10.2 eV 和 3.4 eV ［ ］ 9．4241： 若?粒子(电荷为2e )在磁感应强度为B 均匀磁场中沿半径为R 的圆形轨道运动，则?粒子的德布罗意波长是 (A) )2/(eRB h (B) )/(eRB h (C) )2/(1eRBh (D) )/(1eRBh ［ ］ 10．4770：如果两种不同质量的粒子，其德布罗意波长相同，则这两种粒子的 (A) 动量相同 (B) 能量相同 (C) 速度相同 (D) 动能相同 ［ ］
3. 系综与密度算符 1）纯系综和混合系综 相同的物理体系构成系综，例如由具有自旋的粒子构成的系综。 一个自旋为1/2的粒子的自旋态（方位角,αβ） /2/2(,)(,)(,)cos sin 22i i c c e e ααβ β χαβαβχαβχχχ-++--+-=+=+， 其中,χχ+-是?z s 的本征态， cos(/2)sin(/2) i c c e αββ+-=。 如果所有粒子的自旋都取相同方向，则称体系是极化系统，构成的系综是纯系综。 如果粒子的自旋不在同一方向，则构成的系综叫混合系综。例如自旋向上的粒子数占70％，自旋向下的粒子数占30％，体系是部分极化。一个自旋方向完全随机的系综，其自旋向上，向下的几率各有50％，整的表现是相互抵销，自旋为零，完全没极化。 2）系综平均与态密度算符 系统的力学量平均值 ?A A ααα=， 这里态α是固定的，是量子平均。进入任意表象B ， ,' ?''b b A b b A b b ααα=∑， 对表象的维数求和。 系综平均 [ ]A w A ααα=∑ ， 这里w α是体系处于态α的几率，显然满足归一化条件 1w αα =∑， 是统计平均，求和指标不是对表象的维数，而是对态。例如自旋1/2的粒子构成的系综，自旋表象的维数为2，但不同粒子的自旋态可以有很多取向，求和就是对不同的取向。
[],,','??''''b b b b A w b b A b w b b b A b αααααααα??== ??? ∑∑∑。 定义态密度算符 ?w αα ρ αα=∑， 它在表象B 的矩阵元 '?''bb b w b b αα ρρ αα==∑， []() ,'??????''b b b A b b b A b b A b tr A ρ ρρ==≡∑∑。 这是量子统计力学的基本公式。注意：表象变换不改变矩阵的求迹，上式不依赖于表象的选取。 在连续表象，例如坐标表象，密度算符的矩阵元 *'?''()(')xx x x w x x w x x αααααα ρρααψψ===∑∑ ， 系综平均 []() 3????A tr A d x x A x ρρ==? 。 密度矩阵满足归一化条件 ,,? 1 b b tr w b b w b b w w αααααααα ρ ααα α=====∑∑∑∑完备性条件 态的量子归一化条件 态的统计归一化条件 这里用到了归一化条件1α=和表象的完备性条件1b b b =∑。 设密度算符?ρ的本征态为θ， 22 ?,??ρ θθθρθρθθθθ==＝ 对于纯系综，所有系统都取同一个态n ，
10、量子力学 一、选择题 1．已知一单色光照射在钠表面上，测得光电子的最大动能是1.2 eV ，而钠的红限波长是5400 ?，那么入射光的波长是 (A) 5350 ? (B) 5000 ? (C) 4350 ? (D) 3550 ? 2．在均匀磁场B 内放置一极薄的金属片，其红限波长为λ0。今用单色光照射，发现有电子放出，有些放出的电子(质量为m ，电荷的绝对值为e )在垂直于磁场的平面内作半径为R 的圆周运动，那末此照射光光子的能量是： (A) 0λhc (B) 0 λhc m eRB 2)(2 + (C) 0λhc m eRB + (D) 0λhc eRB 2+ 3．用频率为ν 的单色光照射某种金属时，逸出光电子的最大动能为E K ；若改用频率为2ν 的单色光照射此种金属时，则逸出光电子的最大动能为： (A) 2 E K (B) 2h ν - E K (C) h ν - E K (D) h ν + E K 4．在康普顿效应实验中，若散射光波长是入射光波长的1.2倍，则散射光光子能量ε与反冲电子动能E K 之比ε / E K 为 (A) 2 (B) 3 (C) 4 (D) 5 5．要使处于基态的氢原子受激发后能发射赖曼系(由激发态跃迁到基态发射的各谱线组成的谱线系)的最长波长的谱线，至少应向基态氢原子提供的能量是 (A) 1.5 eV (B) 3.4 eV (C) 10.2 eV (D) 13.6 eV 6．由氢原子理论知，当大量氢原子处于n =3的激发态时，原子跃迁将发出： (A) 一种波长的光 (B) 两种波长的光 (C) 三种波长的光 (D) 连续光谱 7．已知氢原子从基态激发到某一定态所需能量为10.19 eV ，当氢原子从能量为－0.85 eV 的状态跃迁到上述定态时，所发射的光子的能量为 (A) 2.56 eV (B) 3.41 eV (C) 4.25 eV (D) 9.95 eV 8．在气体放电管中，用能量为12.1 eV 的电子去轰击处于基态的氢原子，此时氢原子所能发射的光子的能量只能是 (A) 12.1 eV (B) 10.2 eV (C) 12.1 eV ，10.2 eV 和1.9 eV (D) 12.1 eV ，10.2 eV 和3.4 eV 9．若α粒子(电荷为2e )在磁感应强度为B 均匀磁场中沿半径为R 的圆形轨道运动，则α粒子的德布罗意波长是 (A) )2/(eRB h (B) )/(eRB h (C) )2/(1eRBh (D) )/(1eRBh 10．如果两种不同质量的粒子，其德布罗意波长相同，则这两种粒子的 (A) 动量相同 (B) 能量相同 (C) 速度相同 (D) 动能相同 11．已知粒子在一维矩形无限深势阱中运动，其波函数为： a x a x 23cos 1)(π?= ψ ( - a ≤x ≤a )，那么粒子在x = 5a /6处出现的概率密度为 (A) 1/(2a ) (B) 1/a (C) a 2/1 (D) a /1 12．设粒子运动的波函数图线分别如图(A)、(B)、(C)、(D)所示，那么其中确定粒子动量的精确度最高的波函数是哪个图？ 13．波长λ =5000 ?的光沿x 轴正向传播，若光的波长的不确定量?λ =10-3 ?，则利用不 确定关系式h x p x ≥??可得光子的x 坐标的不确定量至少为： (A) 25 cm (B) 50 cm (C) 250 cm (D) 500 cm x (A) x (C) x (B) x (D)
第三章 力学量用算符表达 §3.1 算符的运算规则 一、算符的定义： 算符代表对波函数进行某种运算或变换的符号。 ?Au v = 表示?把函数u 变成 v ， ?就是这种变换的算符。 为强调算符的特点，常常在算符的符号上方加一个“^”号。但在不会引起误解的地方，也常把“^”略去。 二、算符的一般特性 1、线性算符 满足如下运算规律的算符?，称为线性算符 11221122 ???()A c c c A c A ψψψψ+=+ 其中c 1, c 2是任意复常数，ψ1, ψ2是任意两个波函数。 例如：动量算符?p i =-? ， 单位算符I 是线性算符。 2、算符相等 若两个算符?、?B 对体系的任何波函数ψ的运算结果都相同，即??A B ψψ=，则算符?和算符?B 相等记为??A B =。 3、算符之和 若两个算符?、?B 对体系的任何波函数ψ有：?????()A B A B C ψψψψ+=+=，则???A B C +=称为算符之和。 ????A B B A +=+，??????()()A B C A B C ++=++ 4、算符之积 算符?与?B 之积，记为??AB ，定义为 ????()()AB A B ψψ=?C ψ= ψ是任意波函数。一般来说算符之积不满足交换律，即????AB BA ≠。 5、对易关系 若????AB BA ≠，则称?与?B 不对易。 若A B B A ????=，则称?与?B 对易。 若算符满足????AB BA =-， 则称?A 和?B 反对易。 例如：算符x , ?x p i x ? =-? 不对易
证明：(1) ?()x xp x i x ψψ?=-? i x x ψ? =-? (2) ?()x p x i x x ψψ?=-? i i x x ψψ?=--? 显然二者结果不相等，所以： ??x x xp p x ≠ ??()x x xp p x i ψψ-= 因为ψ是体系的任意波函数，所以 ??x x xp p x i -= 对易关系 同理可证其它坐标算符与共轭动量满足 ??y y yp p y i -= ，??z z zp p z i -= 但是坐标算符与其非共轭动量对易，各动量之间相互对易。 ??0??0y y z z xp p x xp p x -=??-=?，??0??0x x z z yp p y yp p y -=??-=?，??0??0x x y y zp p z zp p z -=???-=?? ????0x y y x p p p p -=，????0y z z y p p p p -=，????0z x x z p p p p -= ????0xy yx -=，????0y z z y p p p p -=，????0z x x z p p p p -= 写成通式(概括起来)： ??x p p x i αββααβδ-= (1) ????0x x x x αββα-= ????0p p p p αββα-= 其中,,,x y z αβ=或1,2,3 量子力学中最基本的对易关系。 注意：当?与?B 对易，?B 与?对易，不能推知?与?对易与否。 6、对易括号(对易式) 为了表述简洁，运算便利和研究量子力学与经典力学的关系，人们定义了对易括号： ??????[,]A B AB BA ≡- 这样一来，坐标和动量的对易关系可改写成如下形式： ?[,]x p i αβαβδ= 不难证明对易括号满足下列代数恒等式： 1) ????[,][,]A B B A =- 2) ???????[,][,][,]A B C A B A C +=+ 3) ?????????[,][,][,]A BC B A C A B C =+ ，?????????[,][,][,]AB C A B C A C B =+，]?,?[]?,?[B A k B k A = 4) ?????????[,[,]][,[,]][,[,]]0A B C B C A C A B ++= ——称为 Jacobi 恒等式。
清华大学《大学物理》习题库试题及答案----10-量子力学习题解读
一、选择题 1．4185：已知一单色光照射在钠表面上， 测得光电子的最大动能是1.2 eV ，而钠的红限波 长是5400 ?，那么入射光的波长是 (A) 5350 ? (B) 5000 ? (C) 4350 ? (D) 3550 ? ［ ］ 2．4244：在均匀磁场B 内放置一极薄的金 属片，其红限波长为λ0。今用单色光照射，发现 有电子放出，有些放出的电子(质量为m ，电荷 的绝对值为e )在垂直于磁场的平面内作半径为 R 的圆周运动，那末此照射光光子的能量是： (A) (B) (C) (D) ［ ］ 3．4383：用频率为ν 的单色光照射某种金 属时，逸出光电子的最大动能为E K ；若改用频 率为2ν 的单色光照射此种金属时，则逸出光电 子的最大动能为： (A) 2 E K (B) 2h ν - E K (C) h ν - E K (D) h ν + E K ［ ］ 4．4737： 在康普顿效应实验中，若散射光 波长是入射光波长的1.2倍，则散射光光子能量 ε与反冲电子动能E K 之比ε / E K 为 (A) 2 (B) 3 (C) 4 (D) 5 ［ ］ 0λhc 0λhc m eRB 2)(2+0λhc m eRB +0λhc eRB 2+
5．4190：要使处于基态的氢原子受激发后能发射赖曼系(由激发态跃迁到基态发射的各谱线组成的谱线系)的最长波长的谱线，至少应向基态氢原子提供的能量是 (A) 1.5 eV (B) 3.4 eV (C) 10.2 eV (D) 13.6 eV ［］ 6．4197：由氢原子理论知，当大量氢原子处于n =3的激发态时，原子跃迁将发出： (A) 一种波长的光(B) 两种波长的光(C) 三种波长的光(D) 连续光谱［］ 7．4748：已知氢原子从基态激发到某一定态所需能量为10.19 eV，当氢原子从能量为－0.85 eV的状态跃迁到上述定态时，所发射的光子的能量为 (A) 2.56 eV (B) 3.41 eV (C) 4.25 eV (D) 9.95 eV ［］ 8．4750：在气体放电管中，用能量为12.1 eV 的电子去轰击处于基态的氢原子，此时氢原子所能发射的光子的能量只能是 (A) 12.1 eV (B) 10.2 eV (C) 12.1 eV，10.2 eV和1.9 eV (D) 12.1 eV，10.2 eV和 3.4 eV ［］ 9．4241：若 粒子(电荷为2e)在磁感应
高等半导体物理 课程内容(前置课程: 量子力学,固体物理) 第一章能带理论,半导体中得电子态 第二章半导体中得电输运 第三章半导体中得光学性质 第四章超晶格,量子阱 前言:半导体理论与器件发展史 1926 Bloch 定理 1931 Wilson 固体能带论(里程碑) 1948 Bardeen, Brattain and Shokley 发明晶体管,带来了现代电子技术得革命,同时也促进了半导体物理研究得蓬勃发展。从那以后得几十年间,无论在半导体物理研究方面,还就是半导体器件应用方面都有了飞速得发展。 1954半导体有效质量理论得提出,这就是半导体理论得一个重大发展,它定量地描述了半导体导带与价带边附近细致得能带结构,给出了研究浅能级、激子、磁能级等得理论方法,促进了当时得回旋共振、磁光吸收、自由载流子吸收、激子吸收等实验研究。 1958 集成电路问世 1959 赝势概念得提出,使得固体能带得计算大为简化。利用价电子态与原子核心态正交得性质,用一个赝势代替真实得原子势,得到了一个固体中价电子态满足得方程。用赝势方法得到了几乎所有半导体得比较精确得能带结构。1962 半导体激光器发明 1968 硅MOS器件发明及大规模集成电路实现产业化大生产 1970 * 超晶格概念提出,Esaki (江歧), Tsu (朱兆祥) * 超高真空表面能谱分析技术相继出现,开始了对半导体表面、界面物理得研究 1971 第一个超晶格Al x Ga1x As/GaAs 制备,标志着半导体材料得发展开始进入人工设计得新时代。 1980 德国得V on Klitzing发现了整数量子Hall 效应——标准电阻 1982 崔崎等人在电子迁移率极高得Al x Ga1x As/GaAs异质结中发现了分数量子Hall 效应 1984 Miller等人观察到量子阱中激子吸收峰能量随电场强度变化发生红移得量子限制斯塔克效应,以及由激子吸收系数或折射率变化引起得激子光学非线性效应,为设计新一代光双稳器件提供了重要得依据。 1990 英国得Canham首次在室温下观测到多孔硅得可见光光致发光,使人们瞧到了全硅光电子集成技术得新曙光。近年来,各国科学家将选择生成超薄层外延技术与精细束加工技术密切结合起来,研制量子线与量子点及其光电器件,预期能发现一些新得物理现象与得到更好得器件性能。在器件长度小于电子平均自由程得所谓介观系统中,电子输运不再遵循通常得欧姆定律,电子运动完全由它得波动性质决定。人们发现电子输运得AharonovBohm振荡,电子波得相干振荡以及量子点得库仑阻塞现象等。以上这些新材料、新物理现象得发现产生新得器件设计思想,促进新一代半导体器件得发展。 半导体材料分类: ?元素半导体, Si, Ge IV 族金刚石结构 Purity 10N9, Impurity concentration 1012/cm3 , Dislocation densities <103 /cm3 Size 20 inches (50 cm) in diameter P V 族 S, Te, Se VI 族 ?二元化合物, 1.IIIV族化合物: GaAS系列,闪锌矿结构, 电荷转移 GaAs, 1、47 eV InAs 0、36 eV GaP, 2、23 eV GaSb, 0、68 eV GaN, 3、3 eV BN 4、6 eV AlN 3、8 eV
物理化学专业博士研究生课程 教学大纲 课程名称：量子化学计算（Computational Quantum Chemistry） 课程编号：B07030411 学分：3 总学时数：72 开课学期：第2学期 考核方式：学习论文 课程说明：（课程性质、地位及要求的描述）。 《量子化学计算》是在学习了《结构化学》、《量子化学》之后，为物理化学专业博士研究生开设的一门方向课，在每学年第二学期讲授。 如果说《结构化学》、《量子化学》还有更多的抽象，那么《量子化学计算》则直接对各研究体系进行可与实验对比的计算机模拟。近二十年来，随着计算机硬件和软件水平的迅速发展，计算化学已成为理论化学的重要分支，主要通过量子化学方法、分子力学方法以及分子动力学模拟来解决与化学相关的问题。目前，计算化学已广泛应用于化学及相关交叉学科的各个领域，迅速成为定量预测分子的结构、性质以及反应性能的有力工具。 本课程计划安排72个学时。采用授课与上机演习相结合的教学方法，使学生在较短时间内掌握当今国际流行的常用计算软件的原理、使用方法及技巧，着重培养同学们解决化学实际问题的能力。要求同学们通过本课程的学习，能对计算化学的原理和方法有一个初步的了解，并能够在化学合成、反应机理、生物、材料等各个领域中得到应用。 教学内容、要求及学时分配： 第一章绪论 内容： 1.1量子力学历史背景 1.221世纪的理论化学计算机模拟
要求：了解量子化学的背景知识、国际国内发展现状及其未来方向学时：4 第二章从头计算法的基本原理和概念 内容： 2.1量子力学基本假设2.2定态近似 2.3从头计算法的“头” 2.4自洽场方法2.5变分法和LCAO-MO近似 2.6量子化学中的一些基本原理和 概念 2.7量子化学中的基本近似 要求：了解从头计算法的基础知识、计算化学中的一些基本原理、概念和近似。 学时：12 第三章布居分析和基组专题 内容： 3.1布居分析 3.2基组专题 要求：理解基组概念及选择的原则，掌握布居分析的计算方法和基组的计数，了解Mulliken布居分析的优缺点及改进的思路。 学时：6 第四章计算方法简介 内容： 4.1半经验方法 4.2HF方法 4.3Post-HF方法 4.4DFT方法 4.5SCF-X 方法 4.6精确模型化学理论方法——Gn 和CBS 4.7赝势价轨道从头计算法 4.8激发态的计算——CIS和CAS 4.9溶剂效应 4.10分子力学和分子动力学基础 要求：了解一些常用计算方法的基本原理及优缺点，重点掌握AM1、INDO、MNDO/PM3、HF、MP、CI、CC、DFT、CAS、溶剂效应等方法的原理，掌握选择计算方法的思路和原则。
量子化学计算方法及应用 吴景恒 实验目的： （1）掌握Gaussian03W的基本操作 （2）掌握 Gaussian03W进行小分子计算的方法，比较不同方法与基组对计算结果的影响，并比较同分异构体的稳定性（3）通过运用量子力学方法计算分子的总电子密度，自旋密度，分子轨道及静电势 实验注意： （1）穿实验服；实验记录用黑色，蓝色或蓝黑色钢笔或签字笔记录；实验数据记录不需要画表格 （2）实验前请先仔细阅读前面的软件使用介绍，然后逐步按照实验步骤所写内容进行操作 （3）截图方法：调整视角至分子大小适中，按下键盘上的PrintScreen按键截图，从“Windows开始菜单”打开“画图”工具，按Ctrl+v或“编辑-粘贴”，去掉四周多余部分只留下分子图形，保存图片 （4）所有保存的文件全部存在E盘或D盘根目录用自己学号命名的文件夹下，不要带中文命名，实验完毕全部删除，不得在计算用机上使用自己携带的U盘或其他便携存储设备！ （5）HyperChem里面截图时候可以用工具栏以下几个工具调整视图： Rotate out-of-plane：平面外旋转工具，转换视角用 Mgnify/Shrink：放大镜工具，转换视角用 Gaussian03W使用介绍：（注意，下面只是界面示意图，实验時切勿按下图设置） 输入文件：Gaussian输入文件，以GJF为文件后缀名 联系命令行：设定中间信息文件（以CHK为后缀名）存放的位置、计算所需的内存、CPU数量等 作业行：指定计算的方法，基组，工作类型，如：#P HF/6-31G(d) Scf=tight Opt Pop=full #作业行开始标记 P 计算结果显示方式为详细， 选择还有T（简单）和 N（常规，默认） HF/6-31G(d) 方法/基组 Opt对分子做几何优化 Pop=full进行轨道布居分析，详尽输出轨道信息和能量 电荷 多重态：分子总电荷及自旋多重态（2S+1, S=n/2, n为成单电子数） 分子结构的表示 1、直角坐标：元素符号X坐标Y坐标Z坐标（如上图所示） 2、Z矩阵（参考后附内容）：元素符号（原子一）原子二键长原子三键角原子四二面角
(1)开壳层,闭壳层 指电子的自旋状态,对于闭壳层,采用限制性计算方法,在方法关键词前面加R 对于开壳层,采用非限制性计算方法,在方法关键词前面加U.比如开壳层的HF就是UHF.对于不加的,程序默认为是闭壳层. 一般采用开壳层的可能性是 1. 存在奇数个电子,如自由基,一些离子 2. 激发态 3. 有多个单电子的体系 4. 描述键的分裂过程 (2) 核磁是单点能计算中另外一个可以提供的数据,在计算的工作设置部分,就是以#开头的一行里,加入NMR关键词就可以了,如 #T RHF/6-31G(d) NMR Test 在输出文件中,寻找如下信息 GIAO Magnetic shielding tensor (ppm) 1 C Isotropic = Anisotropy = 这是采用上面的设置计算的甲烷的核磁结果,所采用的甲烷构形是用B3LYP密度泛函方法优化得到的. 一般的,核磁数据是以TMS为零点的,下面是用同样的方法计算的TMS(四甲基硅烷)的结果1 C Isotropic = Anisotropy = 这样,计算所得的甲烷的核磁共振数据就是,与实验值相比,还是很接近的. (3) 标准几何坐标. 找到输出文件中Standard Orientation一行,下面的坐标值就是输入分子的标准几何坐标. (4) stable 本例中采用SCF方法分析分子的稳定性.对于未知的体系,SCF稳定性是必须要做的.当分子本身不稳定的时候,所得到的SCF结果以及波函数等信息就没有
化学意义. (5)势能面 分子几何构型的变化对能量有很大的影响.由于分子几何构型而产生的能量的变化,被称为势能面.势能面是连接几何构型和能量的数学关系.对于双原子分子,能量的变化与两原子间的距离相关,这样得到势能曲线,对于大的体系,势能面是多维的,其维数取决与分子的自由度. （6）opt Opt=ReadFC 从频率分析(往往是采用低等级的计算得到的)所得到的heckpoint文件中读取初始力矩阵,这一选项需要在设置行之前加入%Chk= filename 一句,说明文件的名称. Opt=CalCFC 采用优化方法同样的基组来计算力矩阵的初始值. Opt=CalcAll 在优化的每一步都计算力矩阵.这是非常昂贵的计算方法,只在非常极端的条件下使用. 有时候,优化往往只需要更多的次数就可以达到好的结果,这可以通过设置MaxCycle来实现.如果在优化中保存了Checkpoint文件,那么使用Opt=Restart可以继续所进行的优化.当优化没有达到效果的时候,不要盲目的加大优化次数.这是注意观察每一步优化的区别,寻找没有得到优化结果的原因,判断体系是否收敛,如果体系能量有越来越小的趋势,那么增加优化次数是可能得到结果的,如果体系能量变化没有什么规律,或者,离最小点越来越远,那么就要改变优化的方法. (7) 频率分析的计算要采用能量对原子位置的二阶导数.HF方法,密度泛函方法(如B3LYP),二阶Moller-Plesset方法(MP2)和CASSCF方法(CASSCF)都可以提供解析二阶导数.对于其他方法,可以提供数值二阶导数. 一般的,对于HF方法,采用计算的频率乘以矫正因子, 方法频率矫正因子零点能矫正因子 HF/3-21G HF/6-31G(d) MP2(Full)/6-31G(d) MP2(FC)/6-31G(d) SVWN/6-31G(d)
1 第二章半导体中的电子状态 要求：掌握能带、有效质量和空穴的概念、 常见半导体的能带结构、杂质能级。 纲要 ?半导体中的电子状态和能带的概念?半导体中电子在外力下的运动，有效质量?材料导电性的差别（半导体如何导电）?空穴的概念 ?常见半导体的能带结构（结论性知识）?杂质和缺陷能级
2 其中 布洛赫波为一个被周期函数u k (r )所调制的平面波。电子的空间分布几率|Ψ*Ψ|=|u *u |是晶格的周期函数。 u k (r )反映电子在原胞内的运动。e ik·r 反映电子在整个晶体中的共有化运动（如何理解？）。 §2.1 半导体中的电子状态和能带（复习） 布洛赫波 采用单电子近似，在晶体的周期场中，电子的波函数为布洛赫波函数： （a n 为任意晶格矢） r k i k k e r u r v v v v v v ?=)()(ψ)()(n k k a r u r u v v v v v +=
3 半导体能带的两种处理方法（半导体中的电子介于自由和强束缚之间） 1、近自由电子近似：从自由出发，加束缚条件 2、紧束缚近似：从束缚出发，减弱束缚 ?电子在原胞中不同位置上出现的几率不同(归一化)。 ?电子在不同原胞的对应位置出现的几率相同（几率能否叠加？），真正反映电子的共有化运动。?如何理解电子的宏观运动？ 电子的波函数几率和统计力学得到电子空间分布的变化，电子空间分布变化表现为宏观的电子运动。
4 例：Si(金刚石结构）能带的形成(SP 3杂化和电子占据能态情况) 1、能级数目； 2、能级间隔； 3、电子不再是局域的。
5 轨道杂化和电子占据情况：价带完全填满；导带没有电子。
实验四 应用量子化学计算方法进行分子结构优化 以及异构化反应研究 Experiment 4. Study on Molecular Structure Optimization and Isomerization Reaction by Using Quantum Chemistry Method 4.1 目的要求 Purpose （1）了解量子化学计算的原理和用途以及几种常用的量子化学计算方法。 （2）熟悉常用量子化学计算软件Gaussian 03的基本使用方法和操作步骤。 （3）掌握如何使用Gaussian 03软件进行分子结构优化和异构化反应过渡态计算。 （4）本实验4学时。 4.2 背景介绍 Background Information 量子化学（quantum chemistry ）以量子力学为理论基础，以计算机为工具，主要通过计算来阐述物质（化合物、晶体、离子、过渡态、反应中间体等）的结构、性质、反应性能及反应机理，研究物质的微观结构与宏观性质的关系，揭示物质和化学反应所具有的特性的内在本质及其规律性[1-4]。随着量子化学计算方法不断发展，计算量以及计算速度不断提高，所计算的体系越来越复杂，现在可以计算有机分子甚至较大分子量的生物分子。 目前常用的量子化学计算软件有Gaussian （https://www.doczj.com/doc/6516660396.html, ）、GAMESS （https://www.doczj.com/doc/6516660396.html,/GAMESS ）、Spartan （https://www.doczj.com/doc/6516660396.html, ）和Molpro （https://www.doczj.com/doc/6516660396.html, ）等。Gaussian 软件是使用最为广泛的量子化学计算软件，支持几乎所有的量子化学计算方法，可以计算得到分子的几乎一切性质，如稳定结构、能量、振动频率、红外和拉曼光谱、NMR 化学位移、轨道能级、静电势、极化率、电离能、电子亲和力、电子密度分布、过渡态和反应途径等。可以模拟在气相和溶液中的体系，模拟基态和激发态等问题。它最早的版本是1970年的Gaussian 70，最新的版本是Gaussian 09。本实验使用的版本为Gaussian 03。 4.3 实验原理 Experimental Principles 4.3.1 量子化学计算方法和特点 多体理论是量子化学的核心问题。n 个粒子构成的量子体系的性质原则上可通过求解n 粒子体系的薛定谔（Schr?dinger ）方程得到体系的波函数来描述。 22 ,111122p q p p i p pq j pi P i p q i j p i Z Z Z E m R ri r ψψ<
公司量化指标计算方法 办公室 1、特殊工种/工序资格一次评定合格率 对特殊工种/工序需要进行资格评定的人员，可以采取抽查几批的办法，对每批计算出一次评定合格率，再把几批一次评定合格率相加除以批次之和。 ２、办公设施完好率 办公设施完好率即对主要办公设施进行检（抽）查，完好的办公设施除以所检（抽）查办公设施总和。 ３、年度失窃次数 本年度被发现失窃次数之和。 工会 经有效评价合格的人数1、职工培训有效性评价合格率＝————————————— — 接受有效性评价的总人数 全厂职工年内接受培训的总时间（学时）２、年度人均接受培训时间＝——————————————————
全厂职工人数 内部培训完成项目＋外出培训完成项目３、培训计划完成率＝——————————————————— 内部培训项目＋外出培训项目 本年度提出合理化建议数 ４、合理化建议人均提出率＝————————————— 全厂职工人数 已采纳和实施数 ５、合理化建议落实率＝—————————————— 本年度提出合理化建议数 采购部 1、计划按时完成率 采购部计划完成率＝按时完成计划数 / 应完成计划数 个人计划完成率＝按时完成计划数 / 应完成计划数 影响因素：a.计划要求时间不合理；b.一份计划包含多项物资，其中一项没有完成，记为该计划没有完成。 2、按时交付率 统计方法：
原材料、摩擦副组件按年度签订意向性合同，每个月按采购计划向供方下达订单，通知供方本月交货的型号、数量、价格、交货期，要求供方审核无误后盖章确认，我公司按供方确认的交货时间接收。 按时交付的数量 / 总的订货数量＝按时交付率 3、质量合格率 统计方法： （1）、摩擦副： 质量不合格有：进货检验不合格、装配废品、市场退回产品。 将出现的不合格分项进行统计，气缸套的不合格分为：金相、尺寸、外观、标识等。 分项不合格数量 / 交付总量X1000000＝分项质量合格率（PPM） 各分项质量合格率相加得到总的质量合格率 （2）原材料 我公司对原材料主要监控其成分。 不合格数量 / 总进货量X100％＝质量合格率（％） 4、超额运费统计 记录每一笔超额运费，每月月底分析超额运费产品的原因，进行整改和预防。 5、产品价格 制定物资采购的最高限价，对于超出限价的采购进行审核，分析原因。对于采购价格的分析，限于目前的知识水平和计算机的应用水平，还不能确定明确的指标。
2.角动量 力学量的一般定义由对易关系确定。 例如坐标和动量的定义是 ??????,0, ,0, ,i j i j i j ij x x p p x p i δ??????==?????? ==， 不依赖于表象。 角动量的一般定义由其对易关系确定： ? J ???,i j ijk k J J i J ε??=?? 或者 ????????????,, ,, ,, x y z y z x z x y J J i J J J i J J J i J J J i J ??????===×?????? 。 满足该定义的力学量称为角动量。显然，轨道角动量? L 满足这个定义。 定义角动量平方 2222?????x y z i J J J J J J =++= ?i ， 描述角动量的大小。 以下用类似于求解谐振子本征值的代数方法来求解2?? ,z J J 的本征值。 1） ， 2????????????????,,,j i i j i j i ijk i k ijk k i ijk i k kji i k J J J J J J J J i J J i J J i J J i J J εεεε??????=+=+=+???????? 0=此处用到了ijk ε的反对称性质：ijk jik εε=?。 故与有共同本征函数。 2?J ?z J 因为不仅仅依赖于，共同本征态至少应有两个量子数，记为2?J ?z J ,m λ。量子数一般无量纲。考虑到和的量纲，记 2?J ?z J 22?,,J m m λλλ= ， ?,,z J m m m λλ= ， 问题：?,?m λ== 2）构造新的算符组 ??? ???????,,? x y x y z x y z J J iJ J J J J J iJ J J ++?+?=+??→=?=????，
数字是具备唯一值的，它能把一片的思维定在一个点上。用空间的123求4法，这次行情的下跌我所有的个股买入操作上都解决了精确价格的问题。公式描述： 1代表第一个波段的起点 2代表第一个波段的终点 3代表第二个波段的起点 买入价格 举例1：600091（明天科技），第一波下跌起点1月17日（波段最高价格）1=9.97元；第一波下跌的终点2月1日（波段最低价格）2=7.50元；第二波下跌起点3月6日（波段最高价格）3=9.80元；求该股的买入价格。 公式：4=2乘3除以1。4=7.50*9.8/9.97=7.38元。 结论：该股再价格下穿7.38元之后，马上受到强大支撑，反身向上，用数字化定量分析得到的精确买入价格，有效。如图。 举例2：600070（浙江富润），第一波下跌起点1月14日（波段最高价格）1=10.10元；第一波下跌的终点2月1日（波段最低价格）2=7.50元；第二波下跌起点3月6日（波段最高价格）3=9.63元；求该股的买入价格。 公式：4=2乘3除以1。4=7.50*9.63/10.1=7.15元。 结论：该股再价格下穿7.15元之后，最低下探至7.02元。马上受到强大支撑，反身向上，用数字化定量分析得到的精确买入价格，有效。如图。
举例3：000595（西北轴承），第一波下跌起点1月25日（波段最高价格）1=10.72元；第一波下跌的终点2月1日（波段最低价格）2=8.10元；第二波下跌起点3月6日（波段最高价格）3=10.17元；求该股的买入价格。 公式：4=2乘3除以1。4=8.1*10.17/10.72=7.69元。 结论：该股再价格下穿7.69元之后，最低下探至7.02元, 受到强大支撑，反身向上，这次用数字化定量分析得到的买入价格虽不够太精确，但仍有效。如图。 举例4：最后一只002219（独一味），由于日线的k线数量比较少，如果你在日线上是找不到取点，就把它切换到小周期。还是可以进行空间定量。我们该股看30分钟线。 第一波下跌起点1=29.25元；第一波下跌的终点2=21.67元；第二波下跌起点3
揭秘量子化学中的计算方法 材料的腐蚀与防护是现代科学技术研究的重要领域之一。添加缓蚀剂是抑制材料腐蚀最为简单、高效、经济的防护手段。因此，缓蚀剂的研究是科研工作者极为关注的重要课题。1971年，Vosta和Eliasek首次用量子化学方法研究了缓蚀剂的缓蚀机理，开创了量子腐蚀电化学。此后，众多学者采用量子化学方法研究缓蚀剂及其作用机理，并取得了丰硕的成果。 什么是量子化学法呢？让我们现在就来看看吧！ 彩虹分割线 量子化学计算方法是……？ 量子化学是以量子力学为基础，利用量子力学的基本原理和方法来研究化学问题的一门学科。 量子化学计算方法的发展 量子化学的基本理论形成于20世纪30年代，随着计算机水平的飞速发展，量子化学计算方法也高速发展起来，并涌现出多种理论算法，并出现大量量子化学计算软件。它们被广泛运用于分子设计、材料性能、化学过程等领域，其方法和结果都显示出了与其他研究手段相比无可比拟的优越性。 密度泛函理论 1、简介 密度泛函理论（DFT）是20世纪60年代在Thomas-Fermi近似理论的基础上发展起来的一种量子力学表达方式。 密度泛函理论指出，电子密度决定分子的一切性质，体系的能量是电子密度的泛函，并指出只要知道体系基态的电子密度，那么体系的一切性质就可以通过量子力学计算得到。 2、优势 密度泛函理论相比于其它量化算法最突出的优势是计算速度快、计算精度高，可用于大分子体系的量化计算，对于含过渡金属体系的计算则更具优势。 3、分类 密度泛函理论种类较多，各种理论的差异在于选择了不同的交换相关能量泛函近似形式。 常见的密度泛函理论算法有以BDW，BLYP为代表的纯密度泛函和以B3LYP，BHANDHLYP为代表的杂密度泛函，而其中最常用于缓蚀剂研究的是B3LYP。 4、应用 随着密度泛函理论体系的日益完善和计算精度的逐步提高，密度泛函理论越来越受到人们的