计算化学公社

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

[Gaussian/gview] 进行两次结构优化的目的

[复制链接 Copy URL]

55

帖子

0

威望

377

eV
积分
432

Level 3 能力者

本帖最后由 11111111m 于 2023-1-18 20:34 编辑

各位老师,我在复现一篇文献,但是文中进行两次结构优化的操作实在是没看懂,它只是为了定向到坐标轴然后旋转九十度的话为什么要进行第二次结构优化,这不是就做了下几何变换,用multiwfn应该就能实现吧。而且它两次优化输出的坐标还是一样的。还请各位老师帮忙解答一下,实在被困住了。两次优化的文件放附件了, 第一次输出文件太大了上不来

202301182003441869..png (140.99 KB, 下载次数 Times of downloads: 3)

202301182003441869..png

202301182003006718..png (160.61 KB, 下载次数 Times of downloads: 1)

202301182003006718..png

1_6_DAP_core.com

3.78 KB, 下载次数 Times of downloads: 0

ground1 input

ground2.rar

6.8 MB, 下载次数 Times of downloads: 0

55

帖子

0

威望

377

eV
积分
432

Level 3 能力者

16#
 楼主 Author| 发表于 Post on 2023-1-28 20:27:48 | 只看该作者 Only view this author
两次优化后的输出坐标

202301282026338246..png (90.02 KB, 下载次数 Times of downloads: 3)

1

1

202301282027227460..png (123.17 KB, 下载次数 Times of downloads: 4)

2

2

55

帖子

0

威望

377

eV
积分
432

Level 3 能力者

15#
 楼主 Author| 发表于 Post on 2023-1-28 20:26:14 | 只看该作者 Only view this author
都是在同一水平下进行的

202301282022044184..png (44.92 KB, 下载次数 Times of downloads: 1)

第一次优化收敛

第一次优化收敛

202301282022327480..png (62.58 KB, 下载次数 Times of downloads: 2)

同一坐标第二次不收敛

同一坐标第二次不收敛

202301282022533440..png (43.48 KB, 下载次数 Times of downloads: 2)

更新坐标后收敛

更新坐标后收敛

55

帖子

0

威望

377

eV
积分
432

Level 3 能力者

14#
 楼主 Author| 发表于 Post on 2023-1-28 20:17:02 | 只看该作者 Only view this author
sobereva 发表于 2023-1-28 16:22
搞清楚当前几何优化任务的位移收敛限是多少
不可能收敛到无穷精确,只能让误差小于收敛限以内

老师,我看了下输出文件,第一次优化的输出坐标在相同水平上再优化的时候它一开始位移不收敛,更新了一次坐标之后才收敛。 我又试着做了第三次优化还是这样,都是位移一开始不收敛。
可是它每次优化的收敛限是一样的,这是精度问题嘛?我把截图放下一楼

5万

帖子

99

威望

5万

eV
积分
112354

管理员

公社社长

13#
发表于 Post on 2023-1-28 16:22:24 | 只看该作者 Only view this author
11111111m 发表于 2023-1-26 11:46
那怎么判断是不是收敛精度的问题呢,关键词都一样最后得出来的坐标也应该是完全一样的吗

搞清楚当前几何优化任务的位移收敛限是多少
不可能收敛到无穷精确,只能让误差小于收敛限以内
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

55

帖子

0

威望

377

eV
积分
432

Level 3 能力者

12#
 楼主 Author| 发表于 Post on 2023-1-26 11:46:34 | 只看该作者 Only view this author
sobereva 发表于 2023-1-26 08:46
可以近似当成一样
如果不是朝向原因,就是几何优化收敛精度原因

那怎么判断是不是收敛精度的问题呢,关键词都一样最后得出来的坐标也应该是完全一样的吗

5万

帖子

99

威望

5万

eV
积分
112354

管理员

公社社长

11#
发表于 Post on 2023-1-26 08:46:16 | 只看该作者 Only view this author
11111111m 发表于 2023-1-25 14:45
sob老师,我看它两次优化的输出坐标小数部分的最后两位有些不同,比如-1.97102和-1.97086,所以想问下这 ...

可以近似当成一样
如果不是朝向原因,就是几何优化收敛精度原因
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

55

帖子

0

威望

377

eV
积分
432

Level 3 能力者

10#
 楼主 Author| 发表于 Post on 2023-1-25 14:45:23 | 只看该作者 Only view this author
sobereva 发表于 2023-1-21 19:18
用这个做法

Multiwfn中非常实用的几何操作和坐标变换功能介绍

sob老师,我看它两次优化的输出坐标小数部分的最后两位有些不同,比如-1.97102和-1.97086,所以想问下这种坐标的精度是什么。因为想用这个坐标做些计算,这两次的坐标可以近似看成一样的吗?

55

帖子

0

威望

377

eV
积分
432

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2023-1-25 14:17:03 | 只看该作者 Only view this author
wzkchem5 发表于 2023-1-19 16:05
结构优化任务里面,高斯先把分子旋转到标准朝向,再做优化。这篇文章只用到了前一个操作,后一个操作是在 ...

老师,第一次优化的输出不是已经在标准朝向下了吗?那第二次优化的前一个操作也没用到吧

5万

帖子

99

威望

5万

eV
积分
112354

管理员

公社社长

8#
发表于 Post on 2023-1-21 19:18:57 | 只看该作者 Only view this author
11111111m 发表于 2023-1-19 12:45
谢谢老师,我还想问一下,如果我想像文章一样让原子坐标差值最大的在x轴其次在yz,只要让分子主轴沿着x轴 ...

用这个做法

Multiwfn中非常实用的几何操作和坐标变换功能介绍
http://sobereva.com/610http://bbs.keinsci.com/thread-24674-1-1.html

(11)令一批原子对应的平面平行于某个笛卡尔平面
程序会让你输入一批原子的序号,然后会通过最小二乘法得到它们的拟合平面,然后体系会被旋转,使得这个拟合平面平行于XY或YZ或XZ平面。这个功能很有实际用处,比如《通过Multiwfn绘制等化学屏蔽表面(ICSS)研究芳香性》(http://sobereva.com/216)里介绍的ICSS_ZZ图对于考察共轭环的芳香性非常有用,当被考察的环是歪斜的时候,你就可以用这个功能先让这个环平行于XY平面,之后就能按文中的做法绘制ICSS_ZZ图了。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

1万

帖子

0

威望

7393

eV
积分
18146

Level 6 (一方通行)

7#
发表于 Post on 2023-1-19 16:05:50 | 只看该作者 Only view this author
11111111m 发表于 2023-1-19 04:27
那是怎么通过结构优化旋转坐标的呢,我看它两次输入文件的关键词也都是一样的

结构优化任务里面,高斯先把分子旋转到标准朝向,再做优化。这篇文章只用到了前一个操作,后一个操作是在白白浪费机时。
所以我说,这个方法纯从人力角度考虑倒是还比较方便,但是考虑了机时以后显然就是Multiwfn更方便了,或者提交一个高斯任务然后过几秒杀死也行。
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

55

帖子

0

威望

377

eV
积分
432

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2023-1-19 12:45:17 | 只看该作者 Only view this author
sobereva 发表于 2023-1-19 11:56
结构优化本身并没有按照特定要求旋转坐标的用处。对于描述或做法不当的文章就别较真了

谢谢老师,我还想问一下,如果我想像文章一样让原子坐标差值最大的在x轴其次在yz,只要让分子主轴沿着x轴就能实现吗

5万

帖子

99

威望

5万

eV
积分
112354

管理员

公社社长

5#
发表于 Post on 2023-1-19 11:56:11 | 只看该作者 Only view this author
11111111m 发表于 2023-1-19 11:27
那是怎么通过结构优化旋转坐标的呢,我看它两次输入文件的关键词也都是一样的

结构优化本身并没有按照特定要求旋转坐标的用处。对于描述或做法不当的文章就别较真了
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

55

帖子

0

威望

377

eV
积分
432

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2023-1-19 11:27:33 | 只看该作者 Only view this author
wzkchem5 发表于 2023-1-18 21:23
可能就是他们不会用Multiwfn或Multiwfn用得不熟,又不怕浪费一次梯度计算的机时。。。
毕竟如果只考虑人力 ...

那是怎么通过结构优化旋转坐标的呢,我看它两次输入文件的关键词也都是一样的

5万

帖子

99

威望

5万

eV
积分
112354

管理员

公社社长

3#
发表于 Post on 2023-1-18 21:27:58 | 只看该作者 Only view this author
没任何做第二次优化的意义。整体旋转不改变内坐标,也不影响能量
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

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

GMT+8, 2024-11-23 23:30 , Processed in 0.203667 second(s), 25 queries , Gzip On.

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