计算化学公社

标题: 如何用cp2k算周期性体系偶极矩 [打印本页]

作者
Author: Priscilla 时间: 2024-1-18 02:20
标题: 如何用cp2k算周期性体系偶极矩
本帖最后由 Priscilla 于 2024-1-19 10:58 编辑

师兄给了我一个任务，叫我计算一个稀土Eu的配合物分子晶体在不同温度条件下测得的不同构型的能量和里面一个片段的偶极矩，想考察能量与偶极矩之间的关系。公社群里问老师说用簇模型算偶极矩，那么我是不是也应该用簇模型算能量？然后我的结构是一个链状交替，大概就是Eu~O~C~O~Eu~O……这样的重复结构，如此该怎么用簇模型？另外是不是周期性结构都不能算偶极矩呢？翻了一下博文都是关于量化偶极矩的计算，如果能算，是不是也要用弥散函数的基组？完全不了解这方面，希望老师们指教

1.师兄要求我删除中间的抗衡分子和配体只计算C-C的键偶极矩和整体构型的能量，这是否合理？如果不删除会有大影响吗？我的策略是能量使用cp2k来算，偶极矩用Gaussian来算。
2.师兄说两层是一样的，只需要算一层即可，所以我采用这样的簇模型，边缘加H
3.算偶极矩时带正电的原子是否可以不加弥散函数？比如这里的Eu

作者
Author: 喝杯热水 时间: 2024-1-18 15:14
CP2K的print有一个moment关键字。你去查一下。

作者
Author: Priscilla 时间: 2024-1-18 23:54

喝杯热水发表于 2024-1-18 15:14
CP2K的print有一个moment关键字。你去查一下。

那个应该只能输出整体的，我需要片段的

作者
Author: sobereva 时间: 2024-1-19 02:56
没具体体系结构没法告诉你怎么用簇模型
算偶极矩和算能量没直接关系，毕竟你又不是通过能量对电场的有限差分算偶极矩
周期性体系没有分子那样的偶极矩的概念，看Journal of Solid State Chemistry, 195 (2012) 2–10

作者
Author: Priscilla 时间: 2024-1-19 10:52

sobereva 发表于 2024-1-19 02:56
没具体体系结构没法告诉你怎么用簇模型
算偶极矩和算能量没直接关系，毕竟你又不是通过能量对电场的有限差 ...

老师好，我已经补充了我要算的体系结构和自己弄得簇模型在一楼，希望老师能指教一下

作者
Author: sobereva 时间: 2024-1-20 06:43
没看出算这种片段偶极矩有什么实际意义

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