| 查看: 6222 | 回复: 29 | |||||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||||
[交流]
【转帖】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类型。 |
回复此楼» 收录本帖的淘帖专辑推荐
 MaterialStudio学习使用 | MS资料 | 量子化学+实验研究 | 第一性原理 |
» 猜你喜欢
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有6人回复
-
孩子确诊有中度注意力缺陷
已经有14人回复
-
三甲基碘化亚砜的氧化反应
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
Gaussian计算得到的前线分子轨道图能说明什么？
已经有7人回复- 如何计算原子某轨道的能量？
已经有12人回复
- MS 计算 晶体 LaZnAsO 的 原子布居分析La的P轨道出现负值？
已经有4人回复
- 【求助】Ms计算磁性要不要设置每个原子的磁矩？？
已经有10人回复
- 【转帖】MS作量子化学计算，和GS一样，可以计算所有原子的轨道系数
已经有8人回复
- 【讨论】深问：原子轨道对分子轨道的贡献，函数如何如何归属到原子轨道里？
已经有3人回复
- 【求助】求：如何分析dmol计算出的氧原子原子轨道图形
已经有4人回复
a94259425
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2337.2
- 散金: 88
- 红花: 5
- 帖子: 195
- 在线: 194.1小时
- 虫号: 860878
- 注册: 2009-09-30
- 专业: 理论和计算化学
★
小木虫(金币+0.5):给个红包，谢谢回帖
小木虫(金币+0.5):给个红包，谢谢回帖
|
楼主还能再详细点么？ 我算的NO结果是这样的 eigenvalue -18.891 -14.125 -1.143 -0.617 -0.465 occupation 1.00 1.00 1.00 1.00 1.00 N 1 1s 0 -0.000009 1.000407 -0.034621 0.002353 0.000006 N 1 2s 0 -0.000148 0.002534 0.388294 -0.667129 0.000033 N 1 2s 0 -0.000061 -0.000305 0.042448 -0.015243 -0.000017 N 1 2p-1 0.000058 -0.000347 -0.033457 -0.006704 0.009159 N 1 2p-1 -0.000001 -0.000093 -0.009896 0.002238 0.001047 N 1 2p 0 0.000083 -0.000496 -0.048063 -0.009616 0.506251 N 1 2p 0 -0.000002 -0.000133 -0.014216 0.003217 0.057731 N 1 2p 1 -0.000364 0.002168 0.210171 0.042069 0.117267 N 1 2p 1 0.000009 0.000578 0.062166 -0.014066 0.013355 N 1 3d-2 0.000080 -0.000131 -0.008255 -0.002118 -0.000697 N 1 3d-1 -0.000019 0.000035 0.001896 0.000418 -0.006055 N 1 3d 0 0.000125 -0.000052 -0.013520 -0.005496 -0.014685 N 1 3d 1 0.000111 -0.000119 -0.011749 -0.003841 0.035104 N 1 3d 2 -0.000246 0.000408 0.025283 0.006347 0.008695 O 2 1s 0 1.000062 -0.000343 -0.027607 0.020070 -0.000006 O 2 2s 0 0.000588 -0.001736 0.630719 0.704510 -0.000046 O 2 2s 0 0.000184 0.000915 0.043490 -0.032370 0.000021 O 2 2p-1 0.000141 -0.000212 0.032394 -0.062756 0.012590 O 2 2p-1 0.000105 0.000125 0.008975 -0.005479 0.000668 O 2 2p 0 0.000205 -0.000305 0.046535 -0.090140 0.695436 - O 2 2p 0 0.000155 0.000180 0.012893 -0.007874 0.036869 - O 2 2p 1 -0.000891 0.001335 -0.203495 0.394180 0.161048 O 2 2p 1 -0.000671 -0.000785 -0.056382 0.034425 0.008520 O 2 3d-2 -0.000048 0.000164 -0.005896 0.007033 0.000654 O 2 3d-1 0.000013 -0.000039 0.001305 -0.001588 0.005695 O 2 3d 0 -0.000032 0.000246 -0.011113 0.012338 0.013809 O 2 3d 1 -0.000050 0.000225 -0.008985 0.010345 -0.033013 O 2 3d 2 0.000150 -0.000508 0.017956 -0.021487 -0.008175 2p（3d）后面的数字是怎么排列的？与2px,2py,2pz是如何对应的？ |
24楼2011-09-23 11:47:36