HyperChem基本操作文档格式.docx
《HyperChem基本操作文档格式.docx》由会员分享,可在线阅读,更多相关《HyperChem基本操作文档格式.docx(19页珍藏版)》请在冰豆网上搜索。
![HyperChem基本操作文档格式.docx](https://file1.bdocx.com/fileroot1/2022-12/13/fe87a8aa-a721-4e1e-8bea-33f8ac7ab07a/fe87a8aa-a721-4e1e-8bea-33f8ac7ab07a1.gif)
选择的级别[原子(atoms),基(residues),或分子(molecules)]。
这里设置为原子(atoms)。
1.左键单击Select菜单。
2.左键单击选择Atoms。
接下来,关闭MultipleSelections:
2.如果MultipleSelections打开,左键单击来关闭它。
这样当你重新选择一个原子或者一组原子的时候,前面的选择将会取消。
选择原子:
1.左键单击Selection工具,把指针移到工作区。
2.左键单击前面创建的原子。
这将激活这个原子。
可以在Preferences对话
框中设置激活的样式。
Preferences在File菜单中。
3.左键单击前面创建的价键中部,这个价键就会激活,价键两端的原子就会
被选择。
4.右键单击工作区的空白处,撤销对所有的原子和价键的选择。
选择一组原子
1.在Select菜单,确保SelectSphere没有打开。
2.在工作区左上方挑选一个点。
3.同时按下鼠标左右键拖向工作区右下角。
HyperChem显示一个代表选择区
域的矩形。
4.继续拖左右键,直到想要的原子被包括在矩形里。
5.释放鼠标按钮。
选择区域所有原子就会被激活。
选择第二组原子:
1.在Select菜单选择MultipleSelections。
打开MultipleSelections
后,在当前的选择中将添加新的选择。
2.像选择第一组原子一样,左右键同时拖。
3.释放鼠标按钮,第二组原子被激活。
如果选择工作区中所有的原子,在选择模式中左键单击工作区的空白区域。
撤销对原子的选择:
1.右键单击未成键原子。
这个原子不再激活。
2.右键单击选择的价键中部。
价键两端的原子不再激活。
3.右键单击工作区的空白区域,撤销对所有原子的选择。
画二维草图1.双击Drawing工具,打开ElementTable对话框。
2.打开Allowions,关闭ExplicitHydrogens。
如果关闭Explicit
Hydrogens,在画分子的时候,不自动添加氢。
3.选择碳,接着关闭对话框。
碳设置为缺省元素。
4.接着按下面操作。
标记原子
1.在Display菜单选择Labels。
2.在Atoms选项框,选择Symbol。
靠近Symbol的单选按钮会被填充。
3.左键单击OK。
Labels对话框消失,所有的原子将标注元素符号。
删除原子:
1.左键单击Drawing工具进入绘图模式。
2.右键单击想要删除的原子,或者想要删除的价键中部。
原子和价键就
会消失。
删除多个原子或价键:
1.左键单击Selection工具进入选择模式。
2.用左右键在工作区拖出选择的目标。
3.在Edit菜单选择Clear。
4.在提示框中选择Yes。
相关回复:
2008-01-02
原子类型(AtomTypes)
在分子动力学计算中,分子中每一个原子都有原子类型。
原子类型代表
具有相似化学环境的一系列同类元素原子。
环境参数包括杂化及其它因
素。
例如,AMBER原子类型C代表一个sp2羰基碳,CT代表四面体碳,CH
代表一个包括一个氢原子的sp3碳。
分子动力学计算用相同的方法处理
同一类型的所有原子。
每一种力场用于一系列不同类型的原子。
原子类
型一般需要分子动力学计算。
但不是半经验计算。
每一种原子类型都有
一个不同于其他力场的名称。
力场需要原子类型来获得计算参数。
MM+
力场可以处理缺省的原子类型(用未知类型**表示)。
当在画一个新的
分子的时候,HyperChem给每个原子一个类型名称**,表示没有指定原子
类型。
当在Build菜单选择ModelBuild或者CalculateTypes,或者改变
力场(使用Setup菜单的MolecularMechanics),HyperChem使用
MolecularMechanicsForceField对话框中选择的力场指定恰当的原子
指定类型后,保留的**意味着没有适当的原子类型可用。
Setup设置
Charge和SpinMultiplicity
在Abinitio和Semi-empirical计算的对话框中出现。
Charge:
指定额外的净剩电荷。
额外电荷定义当前的分子系统是一个电中性
系统,正电系统(阳离子),还是一个负电系统(阴离子)。
SpinMultiplicity:
自旋多重度。
闭壳分子的自旋多重度为1(单重态)。
一个自由基,有一个未成对电子,自旋多重度为2(双重态)。
有两个未成
对电子的系统,自旋多重度一般为3(三重态)。
然而在某些情况下,例如
两个自由基,两个未成对电子也可能产生单态。
State
描述系统中价电子的状态。
包括指定分子处于第一激发单态(Nextlowest)
或者Lowest电子态。
Lowest
给定自旋多重度的最低电子态。
它不一定是基态。
Nextlowest
给定自旋多重度(单重、双重、三重或四重态)的第一电子激发态。
在HyperChem6.0中,semi-empirical方法通常需要给定多重度的最低能态
(lowest)或者次最低能态(nextlowest)的计算。
由于偶数个电子的分子
没有未成对电子,是闭壳层单态,所以只有最低三重态是有效的,而次最低
三重态是无效的。
例如,苯有偶数个电子,并且基态是闭壳层单重态。
我们可
以计算基态(最低单重态),第一激发单重态(次最低单重态),或者第一激
发三重态(最低三重态)。
也就是说,或者HOMO被两个电子占据,或者一个电
子在HOMO,另一个电子在LUMO,产生了激发单重态或者三重态。
对双重态和四
重态,只有给定多重度的最低态(lowest)可用。
UHF选项仅允许给定多重度的最低态(lowest)可用。
例如,可以用UHF选项
研究苯的最低三重激发态,但是不能用来计算单重激发态。
这是因为HyperChem
中的UHF选项不允许任意的轨道占据,也不允许CI计算。
对于RHF选项,可以计算CI波函数。
这个计算由一系列计算得到的RHF轨道开始,
或者从最低单态(或双重态),或者从halfelectron的单态和三态(或双态和
四重态)轨道。
从这些初始态可以得到参考态,从这些参考态,产生一系列微
观态。
通过这些微观态适当的线性组合,可以获得一系列更准确的态。
PM3
PM3方法是对AM1方法的改善,它建立在NDDO近似的基础上。
PM3不同于
AM1仅仅在于参数值不同。
PM3的参数来自对分子性质更大规模的计算。
典型地,未成键相互作用在PM3中比在AM1中有较小的排斥作用。
PM3最
早用于有机分子的计算,但现在对许多主族元素也实现了参数化。
单点计算
单点计算,仅仅执行势能曲面上的一个单个点的计算。
例如对一个双原子分子
来说,这可能是在点a原子间距离R=2.0埃的计算。
单点计算的结果给出系统在
当前几何构型的势能,以及那个点的梯度。
在点b,c,d,或e的单点计算可能
将给出更高的能量。
如果在关键的部分取足够多的点,利用Origin或者Matlab
等数学工具软件,就能描出势能曲线,从而精确算出离解能De和核平衡距离re,
利用公式就能得到震动参数ωe。
对多原子系统来说,状况更加复杂,但是本质是一样的。
Orbital中显示的轨道信息
Alpha&
Beta
显示选择的轨道是alpha自旋还是beta自旋。
LUMO+
显示选择轨道相对于LUMO的位置关系。
例如,选择能量在LUMO之上的
轨道,文本框依次显示+1,+2,+3......;
如果选择能量低于LUMO的
轨道,文本框显示-1,-2,-3......。
HOMO-
显示选择轨道相对于HOMO的位置关系。
例如,选择能量在HOMO之上的
轨道,文本框依次显示-1,-2,-3......;
如果选择能量低于HOMO的
轨道,文本框显示+1,+2,+3......。
Number
显示从能量最低轨道开始的选择轨道绝对数值。
对于UHF计算,轨道的
alpha和beta列编号分开显示,对HOMO-和LUMO+选项也是这样。
Energy
以eV为单位显示选择轨道的能量。
Symmetry
显示选择轨道的不可约表示。
Labels
在轨道显示窗口中,显示每个轨道的电子占据情况和每个轨道的能
量。
按照约定俗成,向上的箭头代表alpha自旋,向下的箭头代表
beta自旋。
对于RHF计算,最大轨道占据是2;
对于UHF计算,最大
轨道占据是1。
ZoomOut和轨道的放大
如果轨道间距太密,在使用Labels时,文字符号会重合在一起,看
起来很不方便。
通过鼠标左键圈出一个或相邻几个选择的轨道,这
些轨道就会被放大。
放大之后,ZoomOut选项可以在窗口中重新回到显示全部轨道范围。
OrbitalDisplay
设置显示选项,显示等高线或等效面。
2DContours
这个选项画选择轨道的二维等高线。
3DIsosurface
这个选项画选择轨道的三维曲面。
Orbitalsquared
选择这一项,图中显示轨道几率密度值的平方。
Options...
显示轨道等值面选项(OrbitalIsosurfaceOptions)的属性菜单。
Plot
按照指定选项显示选择轨道的等高线或等值面。
水中质子的从头计算
用explicithydrogens画H3O+:
1.在Build菜单选择ExplicitHydrogens。
2.选择AllowIons。
3.在Display菜单选择Labels打开Labels对话框。
4.选择Symbols单击OK。
5.在Build菜单选择DefaultElement打开周期表,设置氧元素为缺省元素。
6.左键单击Drawing工具。
画H3O+
1.左键单击工作区创建氧原子。
2.从氧原子画三个价键得到H3O。
在量子力学计算中添加正电荷:
1.在Setup菜单选择AbInitio。
2.单击Options按钮。
3.设置总电荷(Totalcharge)值为1,单击OK关闭Options对话框。
4.单击OK按钮关闭AbInitio对话框。
这样就得到了H3O+。
选择基组
2.为基组选择Other。
3.按下AssignOtherBasisSet按钮。
4.从列表中选择4-31G单击OK。
5.为基组选择Minimal(STO-3G)。
6.选择ApplyBasisSet接着单击OK,或者直接选择OK关闭AbInitio对话框。
观察每个原子应用的基组:
2.选择BasisSet,单击OK。
最小化能量结构
计算H3O+的几何优化:
1.从Setup菜单选择AbInitio。
2.单击Options按钮,保证Totalcharge=1,Spinmultiplicity=1,
Spinpairing=RHF,Conver-gencelimit=0.01,Iterationlimit=50,
Accelerateconver-gence=Yes,为几何优化选择合适的选项。
SinglePoint
only选项在几何优化中不起作用,因而可以是任何值。
3.左键单击OK关闭Options对话框和AbInitio对话框。
4.从Compute菜单中选择GeometryOptimization。
5.选择Polak-Ribiere作为几何优化方法,RMSgradient值为0.1,左键单击OK
开始几何优化。
H3O+优化结果为键角113.74度键长0.99埃。
计算优化后H3O+包括相关能的总能量:
1.从setup菜单选择AbInitio,选择Options按钮并选择计算MP2correlation
energy,依次单击OK关闭对话框。
2.从Compute菜单中选择SinglePoint。
计算结果:
用STO-3G基组计算H3O+总能量,在SCF级别为为-47270.57kcal/mol,
在MP2级别为-47300.98,包括了-30.50kcal/mol相关能。
质子能量
计算H2O+H+优化的STO-3G能量:
1.选择Drawing工具,右键单击H3O+中的一个氢,把它删除形成H2O。
H+的电能是零,所以计算H2O的质子能量从H2O能量中减去前面计算的H3O+能量。
2.从Setup菜单选择AbInitio,选择Options按钮,把Totalcharge从1改变成
0,接着依次单击OK关闭这些对话框。
3.从compute菜单选择GeometryOptimization优化H2O,就像前面优化H3O+一样。
优化结果为键角100.00度,键长0.99埃。
用STO-3G基组计算H2O总能量,在SCF级
别为为-47041.82kcal/mol,在MP2级别为-47066.21kcal/mol,包括了计算的
-24.45kcal/mol相关能。
得到质子能量在SCF级别为228.75kcal/mol,在MP2
级别为234.77kcal/mol。
氨振动态和转动态的从头计算
构造氨:
1.左键单击Build菜单的ExplicitHydrogens使其不被选择。
2.在Build菜单选择DefaultElement。
3.把缺省元素设置为N。
4.左键单击Drawing工具。
5.左键单击工作区,建立一个N原子。
6.在Build菜单选择AddH和ModelBuild,构造NH3。
2.选择原子标志为BasisSet并单击OK关闭对话框。
原子标志显示为当前
的基组。
目前的显示是"
None"
。
3.从Setup菜单选择AbInitio。
4.为基组选择Small(3-21G)接着单击ApplyBasisSet。
选择OK。
这是
每个原子的标志为"
3-21G"
5.用选择工具选择N原子。
6.从Setup菜单选择AbInitio。
7.按下AdvancedOptions按钮,并保证NumberofdOrbitals选择为Six。
单击OK回到AbInitio对话框。
8.按下ExtraBasisFunction按钮,选择DShellType,Exponent设置为
1.0。
依次选择OK关闭AbInitio对话框。
现在添加了有六个基元d形轨道的
额外壳层,每个轨道指数为1.0。
9.右键单击空白区域,取消对N原子的选择。
现在准备用这个基组对氨进
行从头计算。
对锥状氨进行振动分析
1.从Compute菜单选择GeometryOptimization,执行Polak-Ribiere优化,
RMSgradientof的值为0.1。
更好的SCF收敛需要设置为0.01。
这需要更多
的时间。
结果(H-N-Hangle=106.5度,N-Hbondlength=1.01埃)非常
接近实验值。
这是平衡态的基组。
2.从Compute菜单选择SinglePoint,计算包括MP2相关能的总能量。
SCF级
别总能量是-35086.49kcal/mol,包括MP2相关能的总能量是-35191.11
kcal/mol。
3.从Compute菜单选择Vibrations,执行优化几何的振动分析。
这需要几分
钟。
4.从Compute菜单选择VibrationalSpectrum显示振动光谱。
5.确保选择AnimateVibrations。
单击光谱右侧的最低振动能量,并按下
Apply按钮。
移走光谱窗口,在HyperChem就可以看到氨分子最低频振动了。
观察其他振动的方法一样。
6.选择最右侧的最低频跃迁。
按下Apply按钮并单击OK关闭振动光谱窗口。
用X-Y旋转工具旋转氨分子可以彻底观察最低频振动。
7.选择Cancel菜单停止振动。
从Display菜单选择Labels取消里面的任何显
示标志。
从Display菜单选择Renderings。
单击IRVectors和OK显示本征振动
模式。
用X-Y旋转工具旋转氨分子可以彻底观察IRvectors。
8.从Compute菜单选择VibrationalSpectrum再次显示振动光谱。
单击
breathing跃迁并按下Apply按钮观察IRVectors的breathing运动。
注意:
在振动分析计算以前进行执行几何优化是非常重要的,因为基本振动
过程假定初始能量是零。
乙烯最低电子激发态的从头计算
优化乙烯的基态
通过前面几个例子的学习,你已经对从头计算有了初步的了解,所以下面
的一些指令将省略。
用STO-3G基组构造乙烯:
1.在File菜单中选择New,刷新工作区。
2.确保ExplicitHydrogens没有选择。
3.从DefaultElement中选择碳,并画一条C-C单键。
单击碳键中部使它变
成双键。
4.从Build菜单选择AddH和ModelBuild构造乙烯。
5.从Setup菜单选择AbInitio,并选基组为Minimal(STO-3G)。
同时令
Totalcharge=0,Spinmultiplicity=1,Spinpairing=RHF,
Accelerateconvergence=Yes,SCFConvergencelimit=0.0001.按下
CI按钮选择CIMethod为None。
优化乙烯基态:
1.选择Compute菜单的GeometryOptimization。
选择Polak-Ribiere方法,
RMSgradient为0.01。
选择OK关闭对话框。
得到的结果为:
C-Cbondlength:
1.31埃,C-Hbondlength:
1.08埃,H-C-Hangle:
115.7度。
计算相关能:
1.在Setup菜单选择AbInitio,按下Options按钮并选择MP2correlation
energy。
单击OK回到工作区。
2.在Compute菜单选择SinglePoint。
得到结果:
SCF能量:
-48364.64kcal/mol,MP2总能量:
-48438.61kcal/mol,包
括-74.97kcal/mol相关能。
具体运算结果可能会与这个值有微小的差别。
乙烯基态轨道
观察乙烯的轨道和轨道能量图:
1.从Compute菜单选择Orbitals。
拖鼠标放大最高占据轨道(HOMO)和最低空轨道(LUMO)。
这两个轨道归属
为π(HOMO)和π*(LUMO)轨道。
乙烯分子的最低激发态相当于把一个π电子
移到π*轨道上。
这个π-π*跃迁导致激发单重态和三重态。
2.选择这两个轨道,画轨道图。
乙烯的CI和电子光谱
计算乙烯的光谱:
3.单击CI按钮。
4.选择CI方法为SinglyExcited。
5.选择EnergyCriterion而不是OrbitalCriterion。
键入最大值:
1000,
作为MaximumExcitationEnergy。
一直单击OK直到回到工作区。
当执行单重激发CI计算的时候,把大量潜在的组态混合在一起。
这些组态
涉及单电子激发,需要决定包括多少行列式。
这可以通过定义有多少占据
轨道和空轨道或者定义与两个轨道之间的能量差有关的energycriterion
来确定。
这里选择1000将包括CI中所有的单重激发组态。
6.从Compute菜单选择SinglePoint。
7.计算结束,选择Compute菜单中的ElectronicSpectrum。
最强峰位于52.7nm与π-π*跃迁无关。
最低能量跃迁(波长最长)位于最右
端的313.7nm,是三重π-π*跃迁。
当然,这是一个禁戒跃迁(一级近
似)。
最低能量允许跃迁位于106.4nm,是单重π-π*跃迁。
乙烯三重态的几何与能量
由于乙烯的最低激发态(三重)与基态(单重)自旋不同,所以可以用
HyperChem很容易的计算出优化的SCF波函数。
当激发态和基态有相同的对
称性的时候,这不是一个简单的计算,目前Hyper-Chem还不能够完成这个
计算,除非像上面完成的一样,基态提供了CI信息。
而三重态就可以直接
做计算。
计算三重态乙烯的几何:
1.在工作区构造一个乙烯。
我们前面已经得到了这样一个乙烯基态的优化
结构。
2.在平面