| 查看: 1397 | 回复: 12 | |||
111222000版主
|
[交流]
计算断键？？ 已有3人参与
|
|
最近计算一种30多原子的CHON化合物，用hf、3-21g结构优化，结果高斯能正常结束，但是用GV打开一看，我的化合物被优化的四分五裂，本来一种化合物，结果被弄成了四部分碎片，虽然GV中的那根棒棒只是一种形式，单这也太离谱了，求教这是什么原因？如何修改？谢谢 [ Last edited by 111222000 on 2012-5-5 at 22:51 ] |
回复此楼» 猜你喜欢
七嗪类物质合成求助
已经有0人回复
-
[电子教材]天津大学物理化学（第七版）上、下册
已经有4人回复
-
物理化学论文润色/翻译怎么收费?
已经有169人回复
-
阴离子交换膜电解二氧化碳还原
已经有0人回复
-
九江学院2026年最新高层次人才招聘公告
已经有0人回复
-
深圳信息职业技术学院2026年教师招聘公告（最新）
已经有0人回复
-
研究生招生
已经有0人回复
-
山西大学资源与环境工程研究所化工分离/环境功能材料小组招收2026年博士生
已经有3人回复
-
一志愿南师大0703化学 275求调剂
已经有5人回复
-
A 区过线就可申！重庆双非公办高校 0856调剂，接受跨专业！纳米材料 / 润滑热门方向
已经有0人回复
-
欢迎调剂长江师范学院汤老师课题组的硕士研究生
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】自由基耦合反应过渡态有过渡态吗？怎样计算其过渡态？
已经有4人回复
- castep计算类金刚石结构的应力应变曲线
已经有18人回复
- VMD 动态成键 如何定义成断键
已经有11人回复
- NBO计算对于键型判断
已经有8人回复
- 模拟键的断裂是不是不需要做过渡态的计算？
已经有10人回复
- 【500金币本版活动】“经验力场分子模拟+第一性原理计算问题讨论互助及资源共享”
已经有82人回复
- 【求助】用gaussian view打开优化好的结构后,将断裂键连上,是否影响后面计算
已经有9人回复
- 【求助】如何计算化学键的离子性？
已经有4人回复
- 【讨论】大家讨论下计算弱键的方法-氢键，卤键
已经有9人回复
- 【求助】ms怎么计算键能
已经有14人回复
- 【求助】高斯计算过渡态不收敛，请问如何解决？
已经有8人回复
- 【求助】化学键断了
已经有10人回复
zhangmt
管理员
我叫MT
- QC强帖: 5
- 应助: 99 (初中生)
- 金币: 6961.8
- 散金: 10406
- 红花: 49
- 帖子: 1761
- 在线: 763.2小时
- 虫号: 880392
- 注册: 2009-10-22
- 性别: GG
- 专业: 理论和计算化学
2楼2012-05-05 23:58:35
111222000
主管区长
- 应助: 11 (小学生)
- 金币: 1350.7
- 红花: 4
- 帖子: 983
- 在线: 432.2小时
- 虫号: 1782940
- 注册: 2012-04-27
- 性别: GG
- 专业: 爆炸与冲击动力学
3楼2012-05-06 11:43:33
lishijunzong
版主
- 应助: 214 (大学生)
- 金币: 5034
- 红花: 10
- 帖子: 574
- 在线: 195.1小时
- 虫号: 1121020
- 注册: 2010-10-13
- 性别: GG
- 专业: 理论和计算化学
4楼2012-05-06 12:13:40
111222000
超级版主
- 应助: 11 (小学生)
- 金币: 1350.7
- 红花: 4
- 帖子: 983
- 在线: 432.2小时
- 虫号: 1782940
- 注册: 2012-04-27
- 性别: GG
- 专业: 爆炸与冲击动力学
5楼2012-05-06 21:04:03
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
小木虫: 金币+0.5, 给个红包，谢谢回帖
gmy1990: 金币+2 2012-05-06 22:13:31
|
张老师，小卒想就楼主的分子随便说两句 小卒觉得楼主的问题未必是初始结构（至少不是现在需要考虑的问题），而在于计算方法 楼主用的是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 ] |
6楼2012-05-06 22:04:16
111222000
主管区长
- 应助: 11 (小学生)
- 金币: 1350.7
- 红花: 4
- 帖子: 983
- 在线: 432.2小时
- 虫号: 1782940
- 注册: 2012-04-27
- 性别: GG
- 专业: 爆炸与冲击动力学
7楼2012-05-07 08:35:53
zhangmt
兑换贵宾
我叫MT
- QC强帖: 5
- 应助: 99 (初中生)
- 金币: 6961.8
- 散金: 10406
- 红花: 49
- 帖子: 1761
- 在线: 763.2小时
- 虫号: 880392
- 注册: 2009-10-22
- 性别: GG
- 专业: 理论和计算化学
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
根据2001年高压物理学报《利用分子轨道法计算黑索金的生成热和分子结构》文章给出的黑索金结构，属于扭曲的船式结构，但三个硝基的N都是sp2杂化，硝基链接的N原子也是sp2杂化（近似），所以，很难说是平伏键还是直立键。当然，它的结构是PM3优的，可信度难说。另外，NO2换成NF2，N的杂化似乎应该变成Sp3，这样结构的变化又不好预测。感觉上小虫画的两个直立键一个平伏键的椅式结构有一定道理。而三个平伏键的船式结构也未尝不可。这个就像小虫说的那样，楼主有的忙了，要多算几个构型对比一下了。DFT应该是要优于HF的，鉴于F的电负性，弥散应该是要的。变形估计不小，极化也是必须的。就这么几个小原子，计算起来应该是比较快的。 其实，不是说HF就算不好，但是如果结果太离谱，我感觉初始结构不好一般来说占的可能性更大，尤其不是搞计算化学的，做实验的把计算看得很简单，似乎就是个黑匣子，扔里就能出结果。呵呵。要真是这样，该有多好。 |
8楼2012-05-07 13:17:58
zhangmt
兑换贵宾
我叫MT
- QC强帖: 5
- 应助: 99 (初中生)
- 金币: 6961.8
- 散金: 10406
- 红花: 49
- 帖子: 1761
- 在线: 763.2小时
- 虫号: 880392
- 注册: 2009-10-22
- 性别: GG
- 专业: 理论和计算化学
9楼2012-05-07 13:29:18
yjcmwgk
版主
密度泛函·小卒
- QC强帖: 38
- 应助: 165 (高中生)
- 贵宾: 16.999
- 金币: 554.6
- 散金: 39880
- 红花: 395
- 沙发: 2
- 帖子: 10406
- 在线: 1528.1小时
- 虫号: 448368
- 注册: 2007-11-01
- 专业: 理论和计算化学
- 管辖: 科研工具资源
10楼2012-05-07 14:24:46
