24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

111222000

金虫 (正式写手)

应助: 11 (小学生)
金币: 1350.7
红花: 4
帖子: 983
在线: 432.2小时
虫号: 1782940
注册: 2012-04-27
性别: GG
专业: 爆炸与冲击动力学

[交流] 计算断键？？已有3人参与

最近计算一种30多原子的CHON化合物，用hf、3-21g结构优化，结果高斯能正常结束，但是用GV打开一看，我的化合物被优化的四分五裂，本来一种化合物，结果被弄成了四部分碎片，虽然GV中的那根棒棒只是一种形式，单这也太离谱了，求教这是什么原因？如何修改？谢谢

[ Last edited by 111222000 on 2012-5-5 at 22:51 ]

回复此楼

» 猜你喜欢

东方理-中科大联合博士生招聘已经有0人回复
钴酸锂半电池小倍率容量上不去已经有1人回复
物理化学论文润色/翻译怎么收费? 已经有131人回复
吉林大学材料物理本科生求问调剂信息已经有24人回复
2026第二届光电子与半导体器件前沿技术研讨会——光电子赋能·半导体创芯！已经有1人回复
光电子赋能·半导体创芯！已经有1人回复
光电子赋能·半导体创芯！已经有0人回复
Chemical Bonding at Surfaces and Interfaces,最经典的一本表面上化学相互作用教材已经有0人回复
QE计算电声耦合的时候报错Error in routine lambda (100)wrong or too many modes 已经有2人回复
cif求助下载已经有1人回复
高效微纳米储能材料与产业化推进实验室苏州招收硕士研究生（含调剂）已经有2人回复

» 本主题相关价值贴推荐，对您同样有帮助:

【求助】自由基耦合反应过渡态有过渡态吗？怎样计算其过渡态？已经有4人回复
castep计算类金刚石结构的应力应变曲线已经有18人回复
VMD 动态成键如何定义成断键已经有11人回复
NBO计算对于键型判断已经有8人回复
模拟键的断裂是不是不需要做过渡态的计算？已经有10人回复
【500金币本版活动】“经验力场分子模拟+第一性原理计算问题讨论互助及资源共享” 已经有82人回复
【求助】用gaussian view打开优化好的结构后,将断裂键连上,是否影响后面计算已经有9人回复
【求助】如何计算化学键的离子性？已经有4人回复
【讨论】大家讨论下计算弱键的方法-氢键，卤键已经有9人回复
【求助】ms怎么计算键能已经有14人回复
【求助】高斯计算过渡态不收敛，请问如何解决？已经有8人回复
【求助】化学键断了已经有10人回复

1楼 2012-05-05 22:47:29

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

QC强帖: 38
应助: 165 (高中生)
贵宾: 16.999
金币: 554.6
散金: 39880
红花: 395
沙发: 2
帖子: 10406
在线: 1528.1小时
虫号: 448368
注册: 2007-11-01
专业: 理论和计算化学
管辖: 科研工具资源

★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
gmy1990: 金币+2 2012-05-06 22:13:31

引用回帖:

2楼: Originally posted by zhangmt at 2012-05-05 23:58:35:
没有足够的信息，初步猜测是初始结构太离谱导致的。。解决方案，给一个合理的初始结构，重新优化。。

张老师，小卒想就楼主的分子随便说两句
小卒觉得楼主的问题未必是初始结构（至少不是现在需要考虑的问题），而在于计算方法
楼主用的是hf/3-21g（如果仅仅是为了计算速度，为什么不用DMol3这种高速软件，甚至pm6之类的半经验方法呢？）
楼主使用的hf/3-21g，既缺乏电子相关、又缺乏角动量函数、又没有弥散函数。楼主的计算实在是太粗糙了。在我印象中，成熟的力场方法，怕是都比hf/3-21g准确罢。
hf的计算很容易优化“散掉”（换句话说是键长很容易被高估）
既然这个分子是全由短周期主族元素构成的那么不必动用hGGA了看起来最低要求GGA即可，随便挑一个传统GGA即可。比如blyp就挺优秀了
本来觉得极化函数和弥散函数可加可不加，但是考虑到F这种强烈的负电中心，还是加上好。不过高角动量函数就不必了。
得到该分子优化计算的最低要求blyp/6-31+g(d)
张老师您给评判一下小卒的言论？
下面的话说给楼主：
我在GaussView里随便画了一下，发现这个分子居然还有大量异构体的存在。六个氮原子，每个氮原子都有孤对电子的方向问题（这个问题带来的是N锥形结构的差异），而且这个六元环还有船式/椅式互变异构哦（虽然正常情况下椅式占优势，但是在周边被拥挤的取代基占满之后，某些特定的船式构象未必不稳定）。楼主看来有得忙了。
我看楼主的专业是应用化学。如果楼主看了我上面的话有点犯晕，那么我倒是觉得：与其楼主自己去计算，真不如找个专业的量子化学研究生给你去搞计算，你给人家第二作者。没有基础却要搞量子化学计算，得到的结论多半不靠谱。
下图是其中一种异构体

[ Last edited by yjcmwgk on 2012-5-6 at 22:20 ]