计算化学公社

标题: 求助体心立方铁使用高斯优化的问题 [打印本页]

作者
Author: AdamN 时间: 2024-6-10 09:48
标题: 求助体心立方铁使用高斯优化的问题
求助各位，高斯可以用来优化体心立方Fe的结构嘛？之前我构建了一个体心立方的Fe，晶格常数设置为2.866，优化6步之后变成了2.442，但是在尝试扩展其结构时出现了难以收敛的问题是什么原因啊？
# opt freq sdd tpssh

Title Card Required

0 1
Fe                0.00000000 0.00000000 0.00000000
Fe                0.00000000 0.00000000 2.86600000
Fe                0.00000000 2.86600000 0.00000000
Fe                0.00000000 2.86600000 2.86600000
Fe                8.59900000 0.00000000 0.00000000
Fe                8.59900000 0.00000000 2.86600000
Fe                8.59900000 2.86600000 0.00000000
Fe                8.59900000 2.86600000 2.86600000
Fe                1.43300000 1.43300000 1.43300000
Fe                2.86600000 0.00000000 0.00000000
Fe                2.86600000 0.00000000 2.86600000
Fe                2.86600000 2.86600000 0.00000000
Fe                2.86600000 2.86600000 2.86600000
Fe                4.30000000 1.43300000 1.43300000
Fe                5.73300000 0.00000000 0.00000000
Fe                5.73300000 0.00000000 2.86600000
Fe                5.73300000 2.86600000 0.00000000
Fe                5.73300000 2.86600000 2.86600000
Fe                7.16600000 1.43300000 1.43300000

作者
Author: ys_song 时间: 2024-6-10 14:45
感觉你想算的是晶格，但你这样的的输入文件算的是Fe纳米颗粒

作者
Author: sobereva 时间: 2024-6-10 15:34
输入文件里得通过Tv（translation vector）才能定义晶胞，Gaussian才知道当周期性体系计算
然而Gaussian的三维周期性体系计算方面做得很不理想，效率非常低、可控性很差，如果你就是要算当前体系的晶体状态，甭用Gaussian，没希望。CP2K、QE等第一性原理程序是应当优先考虑的

作者
Author: wzkchem5 时间: 2024-6-10 16:14
铁是铁磁性的，自旋多重度显然不可能是1

作者
Author: AdamN 时间: 2024-6-11 10:17

sobereva 发表于 2024-6-10 15:34
输入文件里得通过Tv（translation vector）才能定义晶胞，Gaussian才知道当周期性体系计算
然而Gaussian的 ...

老师，请问这个平移矢量是什么意思啊，应该如何应用呢

作者
Author: sobereva 时间: 2024-6-12 02:49

AdamN 发表于 2024-6-11 10:17
老师，请问这个平移矢量是什么意思啊，应该如何应用呢

算周期性体系自然得定义晶胞，对三维周期性体系需要定义三个晶格矢

怎么应用几句话说不清楚，北京科音中级量子化学培训班（http://www.keinsci.com/workshop/KBQC_content.html）里讲周期性计算部分专门给了详细讲解和很多具体例子

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