计算化学公社

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

[综合交流] 小分子硝基甲烷去1个H之后HOMO计算

[复制链接 Copy URL]

62

帖子

0

威望

1063

eV
积分
1125

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
硝基甲烷失去1个活性H之后,算它的HOMO能级,同时算硝基乙烷失去1个H之后的HOMO能级,为什么硝基甲烷能级会高于硝基乙烷的能级,是因为算法有问题么?用的是B3LYP/6–31G*。感觉实际中硝基乙烷的能级应该更高吧。

6万

帖子

99

威望

5万

eV
积分
120110

管理员

公社社长

2#
发表于 Post on 2015-7-18 06:11:12 | 只看该作者 Only view this author
如果你是指硝基甲烷失去一个氢原子,那么之后是自由基。这种情况没法说HOMO,HOMO仅对于闭壳层才有定义。此时只能说SOMO(限制性开壳层计算),或者独立讨论alpha和beta的HOMO(非限制性计算)。

在ROB3LYP/6-31G*下得出的失去氢原子的硝基甲烷的SOMO能量为-0.14920,低于硝基乙烷失H情况的-0.13387。优化在同级别进行。
北京科音自然科学研究中心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

62

帖子

0

威望

1063

eV
积分
1125

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2015-7-18 08:30:36 | 只看该作者 Only view this author
sobereva 发表于 2015-7-18 06:11
如果你是指硝基甲烷失去一个氢原子,那么之后是自由基。这种情况没法说HOMO,HOMO仅对于闭壳层才有定义。此 ...

对不起,我说的不是很清楚,我想说的是分子失去一个H+,得到的是一个带负电的基团,这样的情况下算的HOMO能级,这种情况需要变成您说的算SOMO么,感觉失去一个H+的时候应该是符合闭壳层的情况吧

6万

帖子

99

威望

5万

eV
积分
120110

管理员

公社社长

4#
发表于 Post on 2015-7-18 08:54:28 | 只看该作者 Only view this author
ca0yan9 发表于 2015-7-18 08:30
对不起,我说的不是很清楚,我想说的是分子失去一个H+,得到的是一个带负电的基团,这样的情况下算的HOMO ...

作为-1离子来算体系仍是闭壳层的,可以说HOMO。
北京科音自然科学研究中心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

62

帖子

0

威望

1063

eV
积分
1125

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2015-7-18 13:09:44 | 只看该作者 Only view this author
sobereva 发表于 2015-7-18 08:54
作为-1离子来算体系仍是闭壳层的,可以说HOMO。

那您看我选取的B3LYP/6–31G*方法合适么,算出来硝基乙烷负离子能级会低于硝基甲烷,感觉正常条件下应该硝基乙烷负离子能级会高于硝基甲烷负离子能级吧,因为同样条件下硝基甲烷能级是低于硝基乙烷能级的。得到这样的结果是因为精度不够造成 的么

6万

帖子

99

威望

5万

eV
积分
120110

管理员

公社社长

6#
发表于 Post on 2015-7-18 13:19:42 | 只看该作者 Only view this author
ca0yan9 发表于 2015-7-18 13:09
那您看我选取的B3LYP/6–31G*方法合适么,算出来硝基乙烷负离子能级会低于硝基甲烷,感觉正常条件下应该 ...


这种计算级别没问题。得到这样的结果就说明是这样。
你也可以改用M062X/6-31+G*再验证一下看看
北京科音自然科学研究中心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

62

帖子

0

威望

1063

eV
积分
1125

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2015-7-18 14:02:06 | 只看该作者 Only view this author
sobereva 发表于 2015-7-18 13:19
这种计算级别没问题。得到这样的结果就说明是这样。
你也可以改用M062X/6-31+G*再验证一下看看

多谢您了

62

帖子

0

威望

1063

eV
积分
1125

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2015-7-18 14:05:10 | 只看该作者 Only view this author

那我想请教您一下,用什么方法能够说明硝基乙烷负离子中C上的孤对电子给电子活性能强于硝基甲烷负离子中C上 的给电子活性呢?

62

帖子

0

威望

1063

eV
积分
1125

Level 4 (黑子)

9#
 楼主 Author| 发表于 Post on 2015-7-18 17:00:26 | 只看该作者 Only view this author
sobereva 发表于 2015-7-18 13:19
这种计算级别没问题。得到这样的结果就说明是这样。
你也可以改用M062X/6-31+G*再验证一下看看

请问您是不是因为给的mol2模型用Dmol3程序优化过,或者我模型弄错了造成的结果又问题。

6万

帖子

99

威望

5万

eV
积分
120110

管理员

公社社长

10#
发表于 Post on 2015-7-19 07:02:46 | 只看该作者 Only view this author
ca0yan9 发表于 2015-7-18 14:05
那我想请教您一下,用什么方法能够说明硝基乙烷负离子中C上的孤对电子给电子活性能强于硝基甲烷负离子中C ...


计算C的f-型condensed Fukui function,然后乘上全局软度(全局硬度的定义是IP减EA再除2,全局软度就是其倒数)或者亲核指数(当前体系的HOMO能量减去tetracyanoethylene的HOMO能量)。对两个体系的碳都这么算然后对比。
北京科音自然科学研究中心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

6万

帖子

99

威望

5万

eV
积分
120110

管理员

公社社长

11#
发表于 Post on 2015-7-19 07:06:13 | 只看该作者 Only view this author
ca0yan9 发表于 2015-7-18 17:00
请问您是不是因为给的mol2模型用Dmol3程序优化过,或者我模型弄错了造成的结果又问题。


这种体系结构很简单,模型应该不会有什么错误。
体系肯定是得经过优化的。
北京科音自然科学研究中心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

62

帖子

0

威望

1063

eV
积分
1125

Level 4 (黑子)

12#
 楼主 Author| 发表于 Post on 2015-7-19 20:32:46 | 只看该作者 Only view this author
sobereva 发表于 2015-7-19 07:02
计算C的f-型condensed Fukui function,然后乘上全局软度(全局硬度的定义是IP减EA再除2,全局软度就是 ...

sob老师,那我需要算“亲核指数”的话,因为HOMO能级算出来的硝基乙烷低于硝基甲烷中间体能级,亲核指数是不是还是硝基乙烷低于硝基甲烷

6万

帖子

99

威望

5万

eV
积分
120110

管理员

公社社长

13#
发表于 Post on 2015-7-20 04:24:04 | 只看该作者 Only view this author
ca0yan9 发表于 2015-7-19 20:32
sob老师,那我需要算“亲核指数”的话,因为HOMO能级算出来的硝基乙烷低于硝基甲烷中间体能级,亲核指数 ...

HOMO越高亲核指数越高
北京科音自然科学研究中心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

62

帖子

0

威望

1063

eV
积分
1125

Level 4 (黑子)

14#
 楼主 Author| 发表于 Post on 2015-7-20 08:52:43 | 只看该作者 Only view this author
sobereva 发表于 2015-7-20 04:24
HOMO越高亲核指数越高

感觉这个结果还是有问题呢,实际上应该硝基乙烷上的C电负性更高吧

62

帖子

0

威望

1063

eV
积分
1125

Level 4 (黑子)

15#
 楼主 Author| 发表于 Post on 2015-7-20 15:43:05 | 只看该作者 Only view this author
老师,我感觉有可能我模型弄错了,硝基甲烷乙烷负离子模型是我自己做的,可能有问题。
您能帮助弄个模型,我再算算,可行?
多谢

本版积分规则 Credits rule

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

GMT+8, 2025-8-14 18:47 , Processed in 0.147466 second(s), 20 queries , Gzip On.

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