计算化学公社

标题: 请问高斯周期性计算定义的盒子大小来避免相互间的影响范围是多少？ [打印本页]

作者
Author: cqu_xiaoq 时间: 2019-4-20 22:01
标题: 请问高斯周期性计算定义的盒子大小来避免相互间的影响范围是多少？
老师您好，如果用高斯做周期性体系计算(尽管高斯计算周期性体系能力较弱，但是对于我研究体系能够满足)，那么定义那个盒子的尺寸使其避免周围的影响，那这个盒子大小的临界值是多少呢？换句话说好比是这个晶胞，我要将这个定义多大，就能将其视为一个孤立体系来计算，使其不受其他晶胞的影响，消除边界效应的影响。谢谢老师

作者
Author: sobereva 时间: 2019-4-20 23:41
最好能把体系特征描述充分
这跟盒子里体系特征有关系。比如之前有PBC算单分子的测试文章，盒子里是个无极性分子，那么真空区不需要很大就可以令误差充分小，而如果是极性分子，达到相同的误差就需要更大的盒子

作者
Author: cqu_xiaoq 时间: 2019-4-21 09:17

sobereva 发表于 2019-4-20 23:41
最好能把体系特征描述充分
这跟盒子里体系特征有关系。比如之前有PBC算单分子的测试文章，盒子里是个无极 ...

老师那这个盒子的大小有没有一个什么临界值，比如说当盒子的X,Y,Z方向都设为20埃，会不会就不受其他的影响了。

作者
Author: sobereva 时间: 2019-4-22 01:03

cqu_xiaoq 发表于 2019-4-21 09:17
老师那这个盒子的大小有没有一个什么临界值，比如说当盒子的X,Y,Z方向都设为20埃，会不会就不受其他的影 ...

20埃绝对够大了

作者
Author: beefly 时间: 2019-4-22 11:49
本帖最后由 beefly 于 2019-4-22 11:50 编辑

设置真空层计算低维体系是因为用了平面波基组造成的。Gaussian、Crystal、ADF-Band等AO基组程序一般不需要设置真空层，因为非周期方向本来就是空的（除非把低维体系当3D周期体系处理）。但是有的AO基组程序在算法上仅支持3D周期体系（siesta可能是这样），这种情况下低维体系还是要设置真空层。

作者
Author: cqu_xiaoq 时间: 2019-4-22 14:58

beefly 发表于 2019-4-22 11:49
设置真空层计算低维体系是因为用了平面波基组造成的。Gaussian、Crystal、ADF-Band等AO基组程序一般不需要 ...

感谢回复，好像gaussian不能设置真空层，综合各位老师的回复，我基本理解了，设定20埃便不会受到边界效应的影响了。

作者
Author: sobereva 时间: 2019-4-23 01:48

cqu_xiaoq 发表于 2019-4-22 14:58
感谢回复，好像gaussian不能设置真空层，综合各位老师的回复，我基本理解了，设定20埃便不会受到边界效应 ...

Gaussian做PBC计算依然有真空层的概念，但是一般没有任何必要设而已，不像那些基于平面波的程序是不得不设

作者
Author: cqu_xiaoq 时间: 2019-4-23 09:31

sobereva 发表于 2019-4-23 01:48
Gaussian做PBC计算依然有真空层的概念，但是一般没有任何必要设而已，不像那些基于平面波的程序是不得不 ...

好的，感谢老师给予纠正。

作者
Author: 洁然不同 时间: 2022-1-14 23:21

sobereva 发表于 2019-4-22 01:03
20埃绝对够大了

Sob老师，请问用CP2K软件计算单个甲醇分子，使用周期性体系方法，cubic cell length取20 Angstrom足够了吗？

作者
Author: sobereva 时间: 2022-1-15 06:51

洁然不同发表于 2022-1-14 23:21
Sob老师，请问用CP2K软件计算单个甲醇分子，使用周期性体系方法，cubic cell length取20 Angstrom足够了 ...

够。而且如果你不要求CP2K考虑周期性的话，盒子也没必要用这么大，用小点的盒子可以大为节约耗时

作者
Author: 洁然不同 时间: 2022-1-16 03:02

sobereva 发表于 2022-1-15 06:51
够。而且如果你不要求CP2K考虑周期性的话，盒子也没必要用这么大，用小点的盒子可以大为节约耗时

Sob老师，我试过小的盒子，能量有差别正常吗？如果不设置周期性计算，结果还是一样的吗？

作者
Author: sobereva 时间: 2022-1-16 05:36

洁然不同发表于 2022-1-16 03:02
Sob老师，我试过小的盒子，能量有差别正常吗？如果不设置周期性计算，结果还是一样的吗？

结果一定会受到盒子尺寸影响，只不过是程度以及是否能接受的问题
你算的本来就是一个孤立分子，当然不应当考虑周期性

作者
Author: 洁然不同 时间: 2022-1-16 14:36

sobereva 发表于 2022-1-16 05:36
结果一定会受到盒子尺寸影响，只不过是程度以及是否能接受的问题
你算的本来就是一个孤立分子，当然不应 ...

谢谢老师！明白了

欢迎光临计算化学公社 (http://bbs.keinsci.com/)