计算化学公社

标题: 关于对聚合物链的力学性能的模拟程序的编写 [打印本页]

作者
Author: alexanderwdz 时间: 2019-2-19 15:46
标题: 关于对聚合物链的力学性能的模拟程序的编写
大家好，我是刚刚入坑cp2k的纯新手。最近在做PETG的拉伸模拟实验，准备用CP2K计算一下各个键的键能。可是很多地方不是很清楚：
1.两个cutoff的值应该怎么选择？
2.基组选的是GTH-BLYP，是不是交换关联势也要选择这个？
3.收敛算法是选择OT还是应该选择DIAG？我看有的帖子说OT适合于像半导体那样能隙大的材料，所以我应该选择DIAG，是这样吗？
目前能像想到的问题就是这些，提前谢过大家

作者
Author: alexanderwdz 时间: 2019-2-19 16:13
在自学CP2K的时候遇到一个词叫real-space integration grid，这个词对应的中文解释是什么？具体是说啥？

作者
Author: sobereva 时间: 2019-2-19 17:23

alexanderwdz 发表于 2019-2-19 16:13
在自学CP2K的时候遇到一个词叫real-space integration grid，这个词对应的中文解释是什么？具体是说啥？

实空间积分格点
很多积分没法解析积分，只能用基于数值格点的方法
比如
密度泛函计算中的格点积分方法
http://sobereva.com/69

单纯对于计算键能的问题，取一段聚合物模型，用Gaussian、ORCA等主流量化程序来算，比CP2K方便得多，而且算法选择上余地大得多，从而可以得到明显更精确的结果

作者
Author: alexanderwdz 时间: 2019-2-19 18:20

sobereva 发表于 2019-2-19 17:23
实空间积分格点
很多积分没法解析积分，只能用基于数值格点的方法
比如

谢谢，经你提醒我想起来这部分知识了！
另外，请问哪里可以找到Gaussian或者ORCA计算键能或者类似计算的tutorial？

作者
Author: sobereva 时间: 2019-2-19 18:30

alexanderwdz 发表于 2019-2-19 18:20
谢谢，经你提醒我想起来这部分知识了！
另外，请问哪里可以找到Gaussian或者ORCA计算键能或者类似计算的 ...

论坛首页google框搜索键能，足够找到能令你正确计算这个问题的信息

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