| 查看: 1638 | 回复: 6 | ||||
[交流]
ETS-NOCV计算：聚合物的片段轨道相互作用、电子转移
|
||||
|
预备知识： ETS-NOCV计算：聚合物的片段轨道相互作用、电子转移ETS-NOCV方法，是将片段形成轨道过程中，电荷密度的形变，由所谓的化学价自然轨道NOCV（natural orbitals for chemical valence）来表征。也就是说，片段形成分子，带来的形变，是NOCV对应的密度的叠加：Δρ( r ) = Σ Δρk( r ) = Σ νk[-ψ2-k( r )+ψ2k( r )] 其中： 1）Δρ( r ) 表示总的形变密度； 2）这些NOCV是成对出现的，本征值（νk）符号相反，例如νk和-νk两个本征值，分别对应两个NOCV ψk( r )和ψ-k( r )，这二者对形变密度的贡献是νk[-ψ2-k( r )+ψ2k( r )]； 3）NOCV是SFO（片段的轨道）的线性组合； 4）NOCV也可以Lowdin基函数的线性组合，其中Lowdin基函数之间是互相正交归一化的，也是SFO的线性组合；直接的SFO基得到的NOCV本征值是一对对绝对值相等、符号相反的本征态，Lowdin基得到的，也是成对的，但本征值不相等。 5）ADF默认的片段分析（Bonding Energy Decomposition）中的Total Orbital Interactions也可以按化学键自然轨道拆分： ΔEorb = ΣΔEkorb = Σ νk[-FTS-k+FTSk]，在后面的out文件中，可以看到各个自然键轨道对ΔEorb的贡献。 ETS-NOCV分析是对ADF默认的片段分析功能的很好补充。 1，建模与结构优化： 聚合物的建模，参考：聚合物的建模 本例以聚乙烷为例，分析C＝C之间的自然价键轨道、能量分解。 结构优化过程，参考：以聚乙烯为例，演示周期性体系的结构优化 2，NOCV计算： 分区，参考：如何创建分区 提交任务，参考：如何提交作业 3，结果查看： 1）查看NOCV能量项： 在ADFinput窗口点击SCM > Output，点击Properties > PEDA-NOCV Energy Terms 各项意义，参考预备知识。 2）查看NOCV轨道： 在ADFinput或者Output窗口点击： SCM > View > Fields → Grid → Fine, Add → Isosurface (Double(+/-) 3）查看某NOCV引起的形变密度（电子从红色区域流向蓝色区域）： 设置NOCV轨道和形变密度透明度、背景色、导出图片，参考ETS-NOCV计算：表面吸附结构的片段轨道相互作用、电子转移 SCM>View>Add>Isosurface>NOCV densities，显示NOCV的密度空间分布； SCM>View>Add>Isosurface: Double(+/-)>NOCV Orbitals，显示NOCV轨道波函数的空间分布（包括两个相位都显示）； SCM>View>Add>Isosurface: Double(+/-)>NOCV Def Densities，显示NOCV之间的密度差（例如NOCV_1**2表示第一个NOCV的模方，也就是密度），电子从红色区域流向蓝色区域。 电子的转移 NOCV是成对出现的（alpha、beta），一般研究第i个NOCV带来的电子转移，要把第i个alpha的NOCV的NOCV Def Densities和第i个beta的NOCV的NOCV Def Densities加起来：SCM> View， Fields > Grid > Fine, Fields > calculate，之后窗口下方 第一个是编号，此处为C－1，后面会用到 第二个是选择要加的第1个量，选择例如NOCV Def Densities列表中第1个 第三个是选择符号，默认为➖，这里改为➕ 第四个是要加的第二个量，选择例如NOCV Def Densities列表中第2个(通常第2n个和第2n＋1个（n>=0）是同一个NOCV的两种自旋) 之后，Add>Isosurface: Double(+/-)，然后在下方Select Field区域选择Other，然后选中列表中的C－1。如此则显示了第1个NOCV的alpha和beta轨道，带来的电子的转移。 根据“NOCV对”（包含alpha和beta）对ΔEorb的贡献大小（例如前面看到，第一个NOCV对ΔEorb的贡献是－388.32211kcal/mol），找到形成该键，主要的是第几个“NOCV对”（包含alpha和beta）的贡献？从而根据该“NOCV对”的Def Densities加和，来判断出形成该键的时候，电子是从什么轨道转移到什么轨道？根据该“NOCV对”的Def Densities加和是sigma型还是Pi型、芳香型，来判断是不是从sigma轨道转移到Pi轨道？等等。 如果将所有“NOCV对”的Def Densities加和，则得到形成该键的总的电荷转移。但一般而言，实际上主要贡献的“NOCV对”只有1～2对。因此一般看前1对，2对的Def Densities加和就能看出规律来了。 |
回复此楼» 收录本帖的淘帖专辑推荐
 能量分解 |
» 猜你喜欢
关于碳量子点抗菌的研究
已经有15人回复
-
南大化学调剂
已经有11人回复
-
物理化学论文润色/翻译怎么收费?
已经有150人回复
-
材料与化工 323 英一+数二+物化，一志愿：哈工大 本人本科双一流
已经有5人回复
-
26申博
已经有2人回复
-
0703调剂
已经有2人回复
-
1
已经有20人回复
-
0703 调剂
已经有5人回复
-
西南某211招27届博士
已经有24人回复
-
大连工业杰青/长江团队-生物质材料-储能电池方向招收2026级博士生
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- Multiwfn支持的弱相互作用的分析方法概览
已经有10人回复
- 【求助】键能的分解
已经有5人回复
» 抢金币啦！回帖就可以得到:
-
持续创业者，A8离异男征婚，事业遇到瓶颈，寻找创业生活伴侣
+1/462
-
中国科学院杭州医学研究-常皓研究员课题组博士后招聘
+1/282
-
福建师范大学叶应祥课题组招收2026级博士研究生
+2/176
-
材料、化学类硕士研究生调剂机会，学硕专硕均可，有意请联系，非诚勿扰！
+1/91
-
浙江大学侯阳教授团队诚聘博士后---高性能固态电解质与固态电池器件设计
+1/90
-
愿遇见你，遇见不止擦肩而过的缘分
+1/74
-
招收2026年材料与化工专硕调剂，全国重点实验室，过国家线即可调剂，名额充足
+2/56
-
飘在美国，有点想家
+1/53
-
湘潭大学“过程强化与绿色化工”创新团队补招2026年秋入学博士生
+1/35
-
新加坡国立大学药学系化学生物学课题组招PhD
+1/11
-
南昌航空大学代威力教授招收光催化和光热催化研究方向博士生
+1/9
-
华中科技大学管理学院招聘社会用工（科研助理）1名
+1/5
-
南京林业大学化学工程学院柏惺峰课题组招收2026年申请考核制博士2名
+1/5
-
上海交通大学机械与动力工程学院博士后招聘
+3/5
-
东北林业大学（211，双一流）宋老师课题组招收高分子、化学、复合材料等方向研究生
+1/4
-
军科某研究院接收调剂研究生
+1/3
-
323分控制工程求调剂
+1/1
-
国家杰青低维材料与器件力学团队2026年招收博士研究生
+1/1
-
浙江大学高纯分离团队诚聘博士后
+1/1
-
南昌航空大学拟招收飞行器设计或计算流体力学研究方向博士生 1名
+1/1
简单回复
怕哦哦哦2楼
2016-06-06 09:58
回复
司徒傲然1981(金币+1): 谢谢参与
tzynew3楼
2016-06-11 15:05
回复
司徒傲然1981(金币+1): 谢谢参与
2016-06-11 15:11
回复
司徒傲然1981(金币+1): 谢谢参与
xhmaohan5楼
2016-06-12 01:25
回复
司徒傲然1981(金币+1): 谢谢参与
nono20096楼
2016-06-12 07:27
回复
司徒傲然1981(金币+1): 谢谢参与
2017-09-05 16:08
回复
