计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 第一性原理 (First-principles) › 优化金属Mg块体和表面结构奇怪，怎样解决

优化金属Mg块体和表面结构奇怪，怎样解决

查看数: 2166 | 评论数: 9 | 收藏 Add to favorites 3

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

h840473807

发布时间: 2022-11-29 20:32

正文摘要:

想要用CP2K计算金属Mg对一个有机物的吸附能，但在优化金属镁晶胞时就出现了问题。构建了晶胞大小为16.0*16.7*22.8的板块模型，含有4层共120个原子，固定最底层优化，初始结构用MS构建，如图1所示，但优化后除了固定 ...

回复 Reply

ultramanm87 发表于 Post on 2025-12-19 12:27:54

lisasever 发表于 2025-12-18 20:08
我今天优化金属Mg时也出现和楼主一样的问题，想问问楼主这个cutoff=2000和REL_cutoff=80是怎么确定的，这个 ...

根据官网上的一个建议“This is for the impatient CP2K user: As a rule of thumb one can multiply the largest exponent of the Gaussian basis sets employed by 40 to obtain a reasonable guess for the cutoff value in Rydberg. The first column of the CP2K basis set format contains the Gaussian exponents, e.g. for the oxygen DZVP-MOLOPT-SR-GTH basis set you may find a value of about 10.4 and the corresponding hydrogen DZVP-MOLOPT-SR-GTH basis set employs a largest exponent of about 10.0. Thus a cutoff value of 10*40 = 400 Rydberg would be a reasonable guess for a cutoff value to start with. ”https://www.cp2k.org/faq:cutoff。Mg DZVP-MOLOPT-SR-GTH DZVP-MOLOPT-SR-GTH-q10的相应数字是30.653047963189，因此至少1200+。进一步的就是自己做截断能测试，相关步骤论坛里一搜就有

lisasever 发表于 Post on 2025-12-18 20:08:38

我今天优化金属Mg时也出现和楼主一样的问题，想问问楼主这个cutoff=2000和REL_cutoff=80是怎么确定的，这个取值确实太大了

h840473807 发表于 Post on 2023-3-28 20:20:57

感谢sob老师的CP2K培训班，让我终于解决了这个问题

h840473807 发表于 Post on 2023-3-28 18:18:24

七尺贱发表于 2023-3-28 17:28
这个问题居然过了三个月才解决，是个大麻烦

因为当时用6-31G*基组了，绕开了这一问题，最近又捡起来才发现是CUTOFF的问题

七尺贱 发表于 Post on 2023-3-28 17:28:14

h840473807 发表于 2023-3-28 15:57
解决了，应该增大CUTOFF。对于金属镁，原来的CUTOFF=800不够，增加到CUTOFF=2000（CUTOFF=1300仍不够）；RE ...

这个问题居然过了三个月才解决，是个大麻烦

h840473807 发表于 Post on 2023-3-28 15:57:08

解决了，应该增大CUTOFF。对于金属镁，原来的CUTOFF=800不够，增加到CUTOFF=2000（CUTOFF=1300仍不够）；REL_CUTOFF=80后，m-DZVP也能正常优化Mg表面结构了。

sobereva 发表于 Post on 2022-11-30 10:18:06

考虑k点。当前大小晶胞，对于金属只考虑gamma点根本不够

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register