|
[交流]
【转帖】MS中原子对前线轨道贡献的计算实例已有17人参与
MS作量子化学计算,和GS一样,可以计算所有原子的轨道系数
因为要用到编辑输入文件,使用关键字,所以好多人不知道
在你的输入文件里任意一行加入
print eigenvectors
然后就可以在输出文件里,就可以看到 all the MO coefficients
To do this you need first to "save files" from the DMol3 or Castep calculation dialog, then edit the input file and finally use "run files."
看看这些轨道系数,找到你的HOMO和LUMO,就知道每一个原子对前线轨道的贡献了
对Castep只适合计算分子,不适合计算晶体。因为晶体的计算是固体物理的能带理论。
举一个简单例子:
比如利用Dmol计算水分子,想看一看每个原子对前线对前线轨道计算
1、建立水分子模型
2、打开Dmol计算对话框,设置计算,你要计算的那些性能都选上
3、点击Run运行,这是就会产生一个input文件,然后停止运行即可
4、打开input文件,在里面# Print options这一行下面(实际上任意一行都可以)写上print eigenvectors,然后保存文件
5、让input文件至于当前对话框,打开Dmol计算对话框,点击Files,然后Run File
6、查看你的outmol文件,里面就有所有原子的轨道组合系数,你就可以看出每个原子对分子轨道的贡献。找到你的HOMO、LUMO,就可以看出你的每个原子对前线轨道的贡献。
对Dmol3分子轨道计算做点补充说明,最近中国科学B化学刊登了一篇关于PbCNN结构的文章,改论文通过QM计算表明在Pb-C之间存在triple bond。我们以这个例子为研究对象,采用Dmol3这个功能进行了计算,泛函采用BLYP形式,波函数机组采用DNP,即Double numerical basis+plarized basis,数值计算精度为fine。计算过程为了反应波函数具体组合分成,输入print eigenvectors 关键词
最后得到如下计算结果:
Symmetry group of the molecule: c*v
Carbon nbas= 1, z= 6, nrfn= 7, rcut= 9.45, e_ref= -0.042518 Ha
rcore= 0.00 zval= 6.00 6.00
n=1 L=0 occ= 2.00 e= -10.059732Ha -273.7394eV
n=2 L=0 occ= 2.00 e= -0.497867Ha -13.5477eV
n=2 L=1 occ= 2.00 e= -0.190572Ha -5.1857eV
n=2 L=0 occ= 0.00 e= -1.480347Ha -40.2823eV
n=2 L=1 occ= 0.00 e= -1.169346Ha -31.8195eV
n=3 L=2 occ= 0.00 e= -2.722222Ha -74.0755eV
n=3 L=2 occ= 0.00 e= -1.388876Ha -37.7932eV eliminated
Nitrogen nbas= 2, z= 7, nrfn= 7, rcut= 9.45, e_ref= -0.107073 Ha
rcore= 0.00 zval= 7.00 7.00
n=1 L=0 occ= 2.00 e= -14.152027Ha -385.0964eV
n=2 L=0 occ= 2.00 e= -0.676218Ha -18.4008eV
n=2 L=1 occ= 3.00 e= -0.258966Ha -7.0468eV
n=2 L=0 occ= 0.00 e= -1.801528Ha -49.0221eV
n=2 L=1 occ= 0.00 e= -1.380332Ha -37.5608eV
n=3 L=2 occ= 0.00 e= -2.722222Ha -74.0755eV
n=3 L=2 occ= 0.00 e= -1.388876Ha -37.7932eV eliminated
Lead nbas= 3, z= 82, nrfn= 18, rcut= 9.45, e_ref= -0.022553 Ha
rcore= 0.00 zval= 82.00 82.00
n=1 L=0 occ= 2.00 e= -2903.263585Ha -79001.8554eV
n=2 L=0 occ= 2.00 e= -489.120791Ha -13309.6596eV
n=2 L=1 occ= 6.00 e= -471.031618Ha -12817.4279eV
n=3 L=0 occ= 2.00 e= -116.566037Ha -3171.9246eV
n=3 L=1 occ= 6.00 e= -107.955115Ha -2937.6094eV
n=3 L=2 occ=10.00 e= -91.875290Ha -2500.0549eV
n=4 L=0 occ= 2.00 e= -25.764798Ha -701.0961eV
n=4 L=1 occ= 6.00 e= -21.985690Ha -598.2613eV
n=4 L=2 occ=10.00 e= -15.007279Ha -408.3690eV
n=4 L=3 occ=14.00 e= -5.554004Ha -151.1322eV
n=5 L=0 occ= 2.00 e= -4.180772Ha -113.7646eV
n=5 L=1 occ= 6.00 e= -2.914053Ha -79.2954eV
n=5 L=2 occ=10.00 e= -0.874176Ha -23.7876eV
n=6 L=0 occ= 2.00 e= -0.331833Ha -9.0296eV
n=6 L=1 occ= 2.00 e= -0.117000Ha -3.1837eV
n=6 L=0 occ= 0.00 e= -0.882962Ha -24.0266eV
n=6 L=1 occ= 0.00 e= -0.626853Ha -17.0575eV
n=6 L=2 occ= 0.00 e= -0.303391Ha -8.2557eV
Point group symmetry c4v symmetry orbital prototypes generated (SYMDEC)
Symmetry orbitals c4v ! nsym ihmi icax nfmt ipr
n norb representation
1 36 a1.1
2 8 b1.1
3 8 b2.1
4 21 e.1
5 21 e.2
total number of valence orbitals: 94
PbCNN分子结构为直线型,属于Cv类型,但DMOL3不能处理某些特殊点群,因此对于这个分子结构DMOL3按照C4v处理的。上面信息可以看到。
我们机组采用的是all electrons,即所有原子轨道群等价处理,机组信息在DMOL3输出文件里面已经写的很清楚了,如给出了主量子数n,角动量子数l以及轨道本征能量,由于是采用DN机组,因此2s轨道有两个,即n=2,l=0出现两次,同理其他轨道,其次是还有极化机组如2p增加d成分,3d增加4f的成分以改善轨道的各向异性。
通过C4v的点群分析,得到PBCNN分子分子轨道的不可约表示类型和数目,如36个A1,8个B1和21个E类型。 |
» 收录本帖的淘贴专辑推荐
» 猜你喜欢
» 本主题相关价值贴推荐,对您同样有帮助:
|