计算化学公社

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

[综合交流] 高斯怎样计算分子间的距离?谢谢指点

[复制链接 Copy URL]

1294

帖子

0

威望

6240

eV
积分
7534

Level 6 (一方通行)

1:背景:想计算分子之间的电荷转移,之前咨询过,http://bbs.keinsci.com/forum.php?mod=viewthread&tid=2357&extra=  
    后来我看相关文章的supporting。我看到了如下的描述。计算电荷转移。
2:我的问题是这样的:
   2.1:这个分子基本都是平行的,这是高斯优化的结构结果吗?
   2.2:这个分子之间的距离和能量之间的这个扫描怎么在高斯中实现?这个我比较想知道。
   2.3:这个分子之间的距离是怎么找的?分子质心之间的距离的话,这个分子的质心怎么借助multiwfn来寻找?然后,两个质心的距离手动算出来?

3:谢谢指点。

0MD([_OSTG8`LPNPS[V4[8V.png (48.4 KB, 下载次数 Times of downloads: 53)

1

1

2.bmp (228.32 KB, 下载次数 Times of downloads: 46)

2

2

3.bmp (187.55 KB, 下载次数 Times of downloads: 58)

3

3
https://www.x-mol.com/groups/fan_jianzhong

strive for greatness

5万

帖子

99

威望

5万

eV
积分
112548

管理员

公社社长

2#
发表于 Post on 2015-12-18 12:14:54 | 只看该作者 Only view this author
2.1 是

2.3 图中没说清楚怎么定义的距离。但对于当前情况,定义为位于中间区域的环之间的平面的间距是合理的
调节平面分子间距的方法
http://sobereva.com/178

2.2 高斯里不好这么扫描,你可以自行写个脚本,按照如上规则批量产生一批单点输入文件,批量执行、批量提取
北京科音自然科学研究中心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!

1294

帖子

0

威望

6240

eV
积分
7534

Level 6 (一方通行)

3#
 楼主 Author| 发表于 Post on 2015-12-18 12:29:57 | 只看该作者 Only view this author
本帖最后由 小范范1989 于 2015-12-18 15:52 编辑

谢谢sob老师指点,也谢谢@lyfchem 的热心指点。
lyfchem老师给我的那个文章我看了,计算流程写的很详细,但是主要是我就只会用高斯,前面的晶体结构的搜索什么的我不太懂。还有那个material studio软件。我再想想计算流程什么的。我就想用高斯+multiwfn近似模拟。
谢谢两位老师的热心指点。ps,我把当时lyfchem老师给我的那个文章再放上来,计算流程很清楚。

2015_NProtocol_Quantitative prediction of charge mobilities of π-stacked system.pdf

1.22 MB, 下载次数 Times of downloads: 126

https://www.x-mol.com/groups/fan_jianzhong

strive for greatness

1294

帖子

0

威望

6240

eV
积分
7534

Level 6 (一方通行)

4#
 楼主 Author| 发表于 Post on 2015-12-21 15:06:10 | 只看该作者 Only view this author
本帖最后由 小范范1989 于 2015-12-21 15:07 编辑
sobereva 发表于 2015-12-18 12:14
2.1 是

2.3 图中没说清楚怎么定义的距离。但对于当前情况,定义为位于中间区域的环之间的平面的间距是合 ...

sob老师好,我想再咨询老师几个问题:
1:如上图所示,我搭建了两个平行的分子A,然后,我想计算一下,两个平行的A之间的距离是多大的时候,能量最低?
我才用了下面的方式:我把距离r,认为的设置3-5ai。然后每隔0.05计算一个单点能,但是这样计算的话计算量太大了,而且我计算出来的数据是:随着两个分子的距离r的增加,他的能量是一直减小的,并没有出现和途中的形式一样。这是什么原因?是不是就应该是随着距离越大,能量越低?而文献中不是采用的这种方式可能,因为文献没说怎么做的。
2:我想学习的一点就是:能不能这样做。固定其中一个分子A,,然后调节另外一个分子和该固定分子的距离,做一个是能量曲线试试,这样的扫描怎么实现?我知道,可以固定分子A,但是我不会使另外一个分子A沿着z轴扫描。也就是使另外一个A分子的z坐标同时增加相同的距离。老师可有什么妙招?谢谢老师。
https://www.x-mol.com/groups/fan_jianzhong

strive for greatness

5万

帖子

99

威望

5万

eV
积分
112548

管理员

公社社长

5#
发表于 Post on 2015-12-21 15:29:30 | 只看该作者 Only view this author
小范范1989 发表于 2015-12-21 15:06
sob老师好,我想再咨询老师几个问题:
1:如上图所示,我搭建了两个平行的分子A,然后,我想计算一下, ...


1 没必要0.05这么小。可以取大点,拟合一下就听平滑了。
一直减小说明一开始距离太近,有交换互斥,但随着继续增加,肯定能量还会上去。

2 你可以先把体系旋转到一个合适的朝向,让A平行于XY平面,然后不断给另一个分子的所有原子的z值增加一个值,起到扫描效果。如果分子本身就完全是平面的话,直接让其z坐标成为被scan的变量即可。
北京科音自然科学研究中心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!

本版积分规则 Credits rule

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

GMT+8, 2024-11-27 20:33 , Processed in 0.426986 second(s), 24 queries , Gzip On.

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