量子化学软件计算:含CHNO中性分子的密度预测
本文是上一篇帖子的进阶,仍然沿用了ChemBioOffice 2014、Gaussian 03w、Multiwfn等软件,并依然采用了Politzer拟合的参数和公式对软件计算结果进行处理,进行密度预测。不同的是本文将以Chem3D为操作平台利用Chem3D对接上的Gaussian 03w进行计算,大大节省了操作步骤,并避免了格式转换中可能出现的错误,如:单键变双键、化学键消失变孤立原子等情况。同时将B3LYP /6 –311G** 更改为B3PW91/6-31G(d,p),主要是因为Politzer处理数据拟合时采用了B3PW91/6-31G(d,p)基组,更适应拟合得到的参数。
一、Chem3D
如图1.1所示,在Chem3D右侧的空白处单击鼠标左键就会弹出和ChemDraw一样的工具栏,这时按照在ChemDraw上同样的操作即可绘制出化学结构式,本文以RDX为例示范后续操作。
1.jpg
图1.1 Chem3D绘图界面
在菜单Calculations(计算)中找到Gaussian Interface(高斯接口),选择Optimize to Transition State(过渡态的优化),如图1.2所示。
2.jpg
图1.2 Chem3D中高斯接口
从弹出界面的Method(任务种类)中选择DFT=B3PW91,设置为6-31G(d,p)基组,如图1.3。
3.jpg
图1.3 设置基组
接着将opt修改为maxcycle = 200,这是为了解决L9999错误,即加大循环数来应对优化不收敛的问题。我在计算稍微复杂的分子就会出现Error termination request processed by link 9999的错误,若是平时计算完成没有错误也可以跳过这一步。
不要忘记设置fchk文件的输出路径,不然无法开始计算,而且fchk文件可直接用Multiwfn打开,不需要使用Gaussian 03将chk转换为fch文件,简化操作步骤,如图1.4。
4.jpg
图1.4 输出路径的设置
采用以上方法后即可点击Run(运行)开始优化,不过以上的方法极为耗时我的联想Y480搭配i5处理器,cpu使用率50%的情况下,优化RDX耗时达3.7h。对于不想耗时那么久的朋友,可以选择Compute Properties(计算性能)任务种类,将得到的chk转换为fch文件用于计算即可,该方法一般仅耗时数十分钟,但精度较上述过渡态优化的方法稍差。
二、计算精度
由于前期计算中我没有注意到L9999错误,在优化不完善的情况下计算了十余种含能材料,浪费了大量的时间。在后期的计算中又引入了加大循环数来解决L9999错误,但随之而来的是计算时间的大大延长如出现错误时RDX加大循环前时长为0.7 h,加大循环后时长为3.7 h,在有限的时间和精力下,我只能将已完成的内容先行发表出来,后续的内容等完善后再开新帖发布。
5.jpg
图2.1 预测结果
文献中使用Politzer公式平均绝对误差为0.036 g/cm3,本文计算12种典型的含能材料后平均绝对误差为0.035 g/cm3,如图2.1所示。其中的硝酸甲酯结果并未计入误差计算中,因为文献中提供的硝酸甲酯密度均为标准状态下的液态密度,而本文所使用的密度预测方法预测的是晶体密度。硝化甘油的液态密度为1.596 g/cm3,固态时密度为1.735 g/cm3,本文预测密度为1.773 g/cm3,可见对于得到的结果更偏向晶体密度。
三、结论
1、本文仅计算得到了12个有效数据,而文献中则计算了36个含能材料。所以本文计算得到的平均绝对误差与文献中略有出入。
2、使用过渡态的优化的任务类型进行密度预测得到的结果精度较高,但极为耗时。在处理复杂结构或需大批量处理数据时可选择计算性能任务类型,该方法平均绝对误差为0.049 g/cm3,依然有着一定的精度。
3、出现L9999错误后的数据等效于计算性能任务类型得到的数据,二者预测的密度一致。
4、在计算了脂肪族硝基类、芳香族硝基类、脂肪族硝胺类、氮杂环硝胺类、氮杂环硝基类、无氧氮杂环类、笼型硝基类、硝酸酯类7个类型的含能材料后,对比实测密度,平均绝对误差仅0.035 g/cm3,符合文献中计算了36个含能材料后平均绝对误差0.036 g/cm3,精度较高可用于新型含能材料设计工作中。
5、数据上仍与文献有出入,在本文方法还有着提升空间。
+1  学术分    虎哥   2016-10-10   优秀的研究工具教程
来自:数理化 / 化学
dracula1429 作者
4年7个月前 修改于 4年7个月前
1楼
在计算ONC时我用chemdraw画的图,导入chem3d后变成了这个德行,按道理结构错误了计算得结果不能计入统计。但是我运行了1天零18个小时,删除数据我实在于心不忍啊。现在先放在上面,我在chem3d上直接画图,以正确结构计算完了,再把数据替换掉。
这短短一千余字真是一个痛苦的过程,总是算着算着有了新发现、新体会,得到更为精确的数据,于是推翻前面已经算好的数据重新开始。这篇帖子后面是我的笔记本没日没夜满负荷运行了半个月的结果。在这期间我真怕这老爷机翘辫子不干了。幸好它撑了下来。
让老爷机歇几天,再投入到离子化合物的计算中。然后摸索看看能不能再简化方法、提高效率。
s.jpg
回复
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
smith
4年7个月前
2楼
不错的软件!学习了
回复
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
dracula1429作者
4年7个月前 修改于 4年7个月前
3楼
在换成正确的结构后,计算得到的八硝基立方烷密度为1.985 g/cm3,与实测密度1.978 g/cm3,仅差0.007 g/cm3。平均绝对误差缩小到0.034 g/cm3,可见绘制结构图时一定要正确,这样得到的结果才更为精确。
onc.jpg
回复
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
dracula1429作者
4年7个月前
4楼
以下是我所使用的软件下载地址:
ChemBioOffice.Ultra.v14.0.0.117安装程序及注册机
attachment icon 49226F9A17BFE2A4040A067B9377DEC50AB3BDE1.torrent 16.58KB TORRENT 49次下载
attachment icon keygen.exe 101.00KB EXE 39次下载
Multiwfn,我是用的是Multiwfn_3.3.8_bin_Win,即Windows版本
http://multiwfn.codeplex.com/
Gaussian03W-E01版本
https://pan.baidu.com/share/link?shareid=2589033093&uk=1275933442&fid=671609856353766
GaussView 5.0.8
https://pan.baidu.com/share/link?shareid=2589033093&uk=1275933442&fid=1063208042313767
MOPAC2016,需发邮件索取激活码
http://openmopac.net/Download_MOPAC_Executable_Step2.html
回复
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
ne555
4年7个月前
5楼
非常感谢分享
回复
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
博丽灵梦
4年7个月前
6楼
直接导入的话,如果不注意键的前后顺序,确实容易挂。。。
回复
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
博丽灵梦
4年7个月前
7楼
另外有个比较233的方法,就是用chem3D自带的简单MM2算法初步优化结构,这样可能有利于收敛
回复
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论
dracula1429作者
4年7个月前
8楼
看文献有先用PM3之类半经验算法处理一下的,不过若是只针对一个化合物,尤其是简单化合物的话不预先优化也行。现在电脑性能也越来越强。不是像我这样非要死磕验证Politzer公式精度的话,睡个觉差不多就能出结果了。
回复
评论
加载评论中,请稍候...
200字以内,仅用于支线交流,主线讨论请采用回复功能。
折叠评论

想参与大家的讨论?现在就 登录 或者 注册

所属专业
上级专业
同级专业
dracula1429
学者 机友 笔友
文章
109
回复
1070
学术分
25
2006/07/09注册,3 天前活动
暂无简介
%7B%22isDisplay%22%3Atrue%7D

仅供内部学术交流或培训使用,请先保存到本地。本内容不代表科创观点,未经原作者同意,请勿转载。

插入资源
全部
图片
视频
音频
附件
全部
未使用
已使用
正在上传
空空如也~
上传中..{{f.progress}}%
处理中..
上传失败,点击重试
等待中...
{{f.name}}
空空如也~
(视频){{r.oname}}
{{selectedResourcesId.indexOf(r.rid) + 1}}
处理中..
处理失败
插入表情
我的表情
共享表情
Emoji
上传
注意事项
最大尺寸100px,超过会被压缩。为保证效果,建议上传前自行处理。
建议上传自己DIY的表情,严禁上传侵权内容。
点击重试等待上传{{s.progress}}%处理中...已上传
空空如也~
草稿箱
加载中...
此处只插入正文,如果要使用草稿中的其余内容,请点击继续创作。
{{fromNow(d.toc)}}
{{getDraftInfo(d)}}
标题:{{d.t}}
内容:{{d.c}}
继续创作
删除插入插入
{{forum.displayName}}
{{forum.countThreads}}
篇文章,
{{forum.countPosts}}
条回复
{{forum.description || "暂无简介"}}
ID: {{user.uid}}
学术分隐藏
{{submitted?"":"投诉或举报"}}
请选择违规类型:
{{reason.description}}
支持的图片格式:jpg, jpeg, png
插入公式
分享回复:{{shareId}}
加载中...
评论控制
加载中...
文号:{{pid}}
加载中...
详情
详情
推送到专栏从专栏移除
设为匿名取消匿名
查看作者
回复
只看作者
加入收藏取消收藏
加入关注取消关注
折叠回复
置顶取消置顶
评学术分
鼓励
设为精选取消精选
建议修改
编辑
通过审核
评论控制
退修或删除
历史版本
违规记录
投诉或举报
加入黑名单移除黑名单
查看IP
{{format('YYYY/MM/DD HH:mm:ss', toc)}}
下载资料
{{fileName}}
大小:{{size}}
下载当前附件将花费 {{costMessage}}
{{description}}
你当前剩余 {{holdMessage}}
{{fileName}}
大小:{{size}}
当前附件免费。
你已购买过此附件,下载当前附件不需要花费积分。
加载中...
{{errorInfo}}
附件已丢失
当前账号的附件下载数量限制如下:
时段 个数
{{f.startingTime}}点 - {{f.endTime}}点 {{f.fileCount}}