计算化学公社

 找回密码 Forget password
 注册 Register
Views: 12321|回复 Reply: 5
打印 Print 上一主题 Last thread 下一主题 Next thread

[Gaussian/gview] 这样的双自由基该如何定义片段?

[复制链接 Copy URL]

41

帖子

0

威望

832

eV
积分
873

Level 4 (黑子)

想优化一个断键后(C1-C2)形成的双自由基体系,如图所示,其中C1和C2各带一个电子自旋相反。我用对称态破损方法(guess=mix)来初步优化了一下体系,发现C1、C2相邻的C-N键由单键变成了双键,这个不是想要的结果。在公社论坛中拜读了Sob版主的“谈谈片段组合波函数”的帖子(http://sobereva.com/82),请问大家这种情况是否能将体系划分成片段来计算?如果可以,该如何划分?请大家指教,谢谢!

diradical.jpg (80.16 KB, 下载次数 Times of downloads: 40)

diradical.jpg

3097

帖子

29

威望

1万

eV
积分
17098

Level 6 (一方通行)

2#
发表于 Post on 2015-11-18 16:46:33 | 只看该作者 Only view this author
GV里的键级是假的,不用管。做一下自旋密度图看单电子在哪儿。

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
sobereva + 1

查看全部评分 View all ratings

6万

帖子

99

威望

5万

eV
积分
120134

管理员

公社社长

3#
发表于 Post on 2015-11-18 23:37:42 | 只看该作者 Only view this author
2L说得很对。

非要用片段初猜,划分这个问题不用太敏感,随便觉得合适划分一下就完了。毕竟只是提供个初猜而已,没那么高的要求。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

41

帖子

0

威望

832

eV
积分
873

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2015-11-19 09:00:58 | 只看该作者 Only view this author
谢谢liyuanhe211和sob版主的建议。

41

帖子

0

威望

832

eV
积分
873

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2015-11-19 17:40:20 | 只看该作者 Only view this author
根据计算结果,电子自旋密度集中在与C1和C2相邻的两个N原子上(-0.6,0.6),而C1和C2上的自旋密度基本为零(-0.03,0.03)。这个不是我希望的结果啊,请问大家这种情况该如何处理?谢谢!

6万

帖子

99

威望

5万

eV
积分
120134

管理员

公社社长

6#
发表于 Post on 2015-11-19 23:59:36 | 只看该作者 Only view this author
happyrabbit 发表于 2015-11-19 17:40
根据计算结果,电子自旋密度集中在与C1和C2相邻的两个N原子上(-0.6,0.6),而C1和C2上的自旋密度基本为零 ...

先用stable关键词测一下是否当前是基态波函数。如果是,换其它泛函算算试试,如果都是这样,说明实际就是这样。
自旋密度分布和预想的不一样是常有的事,这点电脑比人脑准。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-15 06:54 , Processed in 0.193705 second(s), 26 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list