如何利用cp2k计算熔体中掺杂的一个原子的化学势

nusiew

大家好，请问一下大家有没有基于cp2k从头算计算过掺杂原子化学势呢？我想要计算一个掺杂原子在熔体中的化学势，看到文献中曾经有报道（https://doi.org/10.1038/35012056），通过计算一个掺杂原子在体系中进行原子取代的吉布斯自由能变或者亥姆霍兹自由能变来计算化学势，根据化学势的定义

想请问，如果我用cp2k的NVT系综或者NPT系综分别来跑掺杂前后的体系，获取能量变化，然后按照上述化学势与A和G的关系，能够准确计算出掺杂原子的化学势呢？谢谢大家！

get-it

原理上没错，但是你需要算自由能而不是势能，而算液体的自由能需要用TI（你引的文章中也提到了），这个计算量很大。

nusiew

get-it 发表于 2022-3-13 16:27
原理上没错，但是你需要算自由能而不是势能，而算液体的自由能需要用TI（你引的文章中也提到了），这个计算 ...

谢谢回复，文章我看的不是太明白，请问TI是啥，用cp2k计算自由能好算不，可不可以详细的说一下，我不太懂这一部分，希望学习一下，非常感谢！

get-it

nusiew 发表于 2022-3-13 16:46
谢谢回复，文章我看的不是太明白，请问TI是啥，用cp2k计算自由能好算不，可不可以详细的说一下，我不太懂 ...

thermodynamic integration,要看懂文章还得靠自己.
我没多少cp2k使用经验,没办法提供建议.

喵星大佬

get-it 发表于 2022-3-13 17:29
thermodynamic integration,要看懂文章还得靠自己.
我没多少cp2k使用经验,没办法提供建议.

AIMD是怎么实现炼金术的

get-it

喵星大佬 发表于 2022-3-13 23:25
AIMD是怎么实现炼金术的

我的理解是，因为固相和液相都有好的reference system，它们在固液相中的采样分布和AIMD足够接近，所以aimd的耗时也可以接受。例如固相用声子谱的简谐哈密顿量，液相用一个拟合的LJ势。具体的可以看楼主发的链接。