计算化学公社

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

[Gaussian/gview] G09 势能面扫描问题

[复制链接 Copy URL]

156

帖子

0

威望

3743

eV
积分
3899

Level 5 (御坂)

跳转到指定楼层 Go to specific reply
楼主
我做relaxed势能面扫描,关键词Opt=modredundant。我想把某原子沿着X轴方向移动, X 33 S 15 0.3 0 0.0 0 0.0。

可是输出文件里读出来的却是 X 33 S 15 0.3 15 0.3 15 0.3。
如果沿着Z轴,我写 X 33 S 0 0.0 0 0.0 15 0.3
输出文件就是只有 X 33 S
但我看手册似乎可以只沿着某个方向,请问这是怎么回事?
Cartesian coordinates can also include scan specifications:
X atom S x-steps x-size y-steps y-size z-steps z-size
For example, the following ModRedundant input performs a relaxed potential energy surface scan. Starting at the initial position of atom 1 and moving in five 0.2 Angstrom steps in the X direction and three 0.1 Angstrom steps in the Z direction:
X 1 S 5 0.2 0 0.0 3 0.1

6万

帖子

99

威望

6万

eV
积分
125136

管理员

公社社长

2#
发表于 Post on 2015-3-16 23:49:45 | 只看该作者 Only view this author
你用的哪个版本?
若是G09 D.01,在读冗余内坐标的时候有bug,换老版本即可。
X 33 S 15 0.3 0 0.0 0 0.0这样的写法至少在G09 B.01上没问题。
北京科音自然科学研究中心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

156

帖子

0

威望

3743

eV
积分
3899

Level 5 (御坂)

3#
 楼主 Author| 发表于 Post on 2015-3-17 00:20:46 | 只看该作者 Only view this author
sobereva 发表于 2015-3-16 23:49
你用的哪个版本?
若是G09 D.01,在读冗余内坐标的时候有bug,换老版本即可。
X 33 S 15 0.3 0 0.0 0 0.0 ...

这么巧,还真是D01,换个版本。。

156

帖子

0

威望

3743

eV
积分
3899

Level 5 (御坂)

4#
 楼主 Author| 发表于 Post on 2015-3-20 04:44:13 | 只看该作者 Only view this author
sobereva 发表于 2015-3-16 23:49
你用的哪个版本?
若是G09 D.01,在读冗余内坐标的时候有bug,换老版本即可。
X 33 S 15 0.3 0 0.0 0 0.0 ...

Can't type Chinese here. Sorry for the inconvenience .
I want to do the relaxed energy surface scan along the vertical distance of two ion pairs (planar molecues). Each pair is constructed by one cation and one counterion. So I try to scan the distance between two atoms that are in the center of the two ion pairs, respectively. However, the distance never increase. Can you please check the attached input and output files. If possilbe, explain the reason. Anything can help.
Thank you!

C7PHNX_H1.log

1.63 MB, 下载次数 Times of downloads: 295

C7PHNX_H1.com

5.73 KB, 下载次数 Times of downloads: 21

6万

帖子

99

威望

6万

eV
积分
125136

管理员

公社社长

5#
发表于 Post on 2015-3-20 05:01:05 | 只看该作者 Only view this author
柔性扫描的第一步扫描的优化过程一直没收敛。所以显然也看不到距离在动。

对于你的目的,柔性扫描未必合适。只是不断变化中间那两个原子的距离的话,周围的分子片段在优化过程中会有很大变形,没法表现出两个平面分子的垂直距离对能量的影响。而且柔性扫描也太耗时。

就用刚性扫描就可以了。用笛卡尔坐标,把其中一个平面分子的所有原子的Z坐标都用相同的变量表示,然后扫描这个Z坐标。
北京科音自然科学研究中心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

156

帖子

0

威望

3743

eV
积分
3899

Level 5 (御坂)

6#
 楼主 Author| 发表于 Post on 2015-3-24 02:24:37 | 只看该作者 Only view this author
sobereva 发表于 2015-3-20 05:01
柔性扫描的第一步扫描的优化过程一直没收敛。所以显然也看不到距离在动。

对于你的目的,柔性扫描未必合 ...

谢谢Sob老师,我用了刚性扫描,感觉良好。现在分子式垂直移动的,请问那我想两个分子平行移动,应该怎么设置scan,昨天试了一下扫描角度+距离感觉不太对。应该距离+角度+二面角吗,感觉二面角不太好控制。又试了一下虚原子,好像还是没用。请赐教!

6万

帖子

99

威望

6万

eV
积分
125136

管理员

公社社长

7#
发表于 Post on 2015-3-24 05:21:59 | 只看该作者 Only view this author
superrice 发表于 2015-3-24 02:24
谢谢Sob老师,我用了刚性扫描,感觉良好。现在分子式垂直移动的,请问那我想两个分子平行移动,应该怎么 ...

写个脚本吧,利用循环,每次给要平移的分子的x或y坐标增加一点,生成一批单点输入文件。然后批量执行、批量提取能量。
北京科音自然科学研究中心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

156

帖子

0

威望

3743

eV
积分
3899

Level 5 (御坂)

8#
 楼主 Author| 发表于 Post on 2015-3-25 22:02:49 | 只看该作者 Only view this author
sobereva 发表于 2015-3-24 05:21
写个脚本吧,利用循环,每次给要平移的分子的x或y坐标增加一点,生成一批单点输入文件。然后批量执行、批 ...

Thx Sob. It helps!

349

帖子

0

威望

905

eV
积分
1258

Level 4 (黑子)

9#
发表于 Post on 2020-7-24 22:38:16 | 只看该作者 Only view this author
sobereva 发表于 2015-3-16 23:49
你用的哪个版本?
若是G09 D.01,在读冗余内坐标的时候有bug,换老版本即可。
X 33 S 15 0.3 0 0.0 0 0.0 ...

老师,您好,请问进行主客体复合物的势能面扫描时,是沿着Z轴方向移动,有没有什么方法使其沿着X轴或者Y轴进行移动(我想确定二聚体间堆积的最佳距离,都是有机小分子,需要从X,Y,Z三个方向分别进行考虑,最终确定最佳位置),辛苦啦

6万

帖子

99

威望

6万

eV
积分
125136

管理员

公社社长

10#
发表于 Post on 2020-7-25 00:10:30 | 只看该作者 Only view this author
xxzj 发表于 2020-7-24 22:38
老师,您好,请问进行主客体复合物的势能面扫描时,是沿着Z轴方向移动,有没有什么方法使其沿着X轴或者Y ...

自己写脚本产生要算的坐标最省事,用其它乱七八糟做法反倒麻烦
北京科音自然科学研究中心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

349

帖子

0

威望

905

eV
积分
1258

Level 4 (黑子)

11#
发表于 Post on 2020-7-25 11:15:36 | 只看该作者 Only view this author
sobereva 发表于 2020-7-25 00:10
自己写脚本产生要算的坐标最省事,用其它乱七八糟做法反倒麻烦

老师,我现在能力不足,还不会写脚本。通过帖子内容,很快确定了z轴位置,现在就差其它2两个方向。想问一下,还有其它方法吗

6万

帖子

99

威望

6万

eV
积分
125136

管理员

公社社长

12#
发表于 Post on 2020-7-26 01:37:16 | 只看该作者 Only view this author
xxzj 发表于 2020-7-25 11:15
老师,我现在能力不足,还不会写脚本。通过帖子内容,很快确定了z轴位置,现在就差其它2两个方向。想问一 ...

没其它办法。如果你特有耐心的话,在excel里手动一个一个生成坐标,比如x统一加0.2,保存一次文件,加0.4,保存一次文件...

不会写点脚本、程序,碰上点非标准类型的计算任务就会很麻烦。
北京科音自然科学研究中心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

349

帖子

0

威望

905

eV
积分
1258

Level 4 (黑子)

13#
发表于 Post on 2020-7-26 07:39:48 | 只看该作者 Only view this author
sobereva 发表于 2020-7-26 01:37
没其它办法。如果你特有耐心的话,在excel里手动一个一个生成坐标,比如x统一加0.2,保存一次文件,加0.4 ...

收到,谢谢老师

349

帖子

0

威望

905

eV
积分
1258

Level 4 (黑子)

14#
发表于 Post on 2020-7-26 13:47:09 | 只看该作者 Only view this author
sobereva 发表于 2020-7-26 01:37
没其它办法。如果你特有耐心的话,在excel里手动一个一个生成坐标,比如x统一加0.2,保存一次文件,加0.4 ...

老师,对二聚体的势能面扫描就是计算单点能吗

1665

帖子

5

威望

4788

eV
积分
6553

Level 6 (一方通行)

喵星人

15#
发表于 Post on 2020-7-26 15:30:20 | 只看该作者 Only view this author
xxzj 发表于 2020-7-26 13:47
老师,对二聚体的势能面扫描就是计算单点能吗

刚性扫描就是一系列单点能

本版积分规则 Credits rule

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

GMT+8, 2026-2-20 10:58 , Processed in 0.183899 second(s), 23 queries , Gzip On.

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