Exodus 发表于 2021-7-2 13:31
首先谢谢老师您的回复！
1. 我通过nat_todo指定了一个双原子分子，也确实得到了6个实频，但同时也出现6 ...

你好，我用QE算频率也遇到相同的问题了，你的问题解决了吗？

我最近也要做类似的计算。如果写一点代码，应该可以用phonopy实现。
但是我在考虑对称性的问题。如果结构加上衬底，对称性会改变。到底如何处理才是正确的？

卡开发发 发表于 2021-7-2 13:19
1、按道理应该通过nat_todo指定参与计算原子数目，然后再后面指定这些原子的编号就行，理论上能实现固定的 ...

首先谢谢老师您的回复！
1. 我通过nat_todo指定了一个双原子分子，也确实得到了6个实频，但同时也出现6个虚频。其次，除了收敛精度，ecuwfc,ecutrho,kpoints，tr2_ph我也测试过了，并没有得到有效的改进。我再尝试更换asr试试。
2.好的，我尝试一下您的这个方法

1、按道理应该通过nat_todo指定参与计算原子数目，然后再后面指定这些原子的编号就行，理论上能实现固定的目的。你可以检查下频率的数目是否与未被固定的原子呈现3N的关系。如果没问题，那就得想办法消除虚频，除了scf和relax的精度，你的动能截断k点以及声子的tr2_ph也有影响。如果这些都严格测试了，也可以尝试使用声学支求和规则(asr)。

2、Phonopy做固定原子我不确定行不行，我这种处理不是很了解，看看其他人是否有相关心得。如果实在不行可以试试这样，比如不需要参与移动的自由度对应的image的结果可以拷贝未移动原子时得到的数据来骗过程序，这样这部分对Hessian的贡献就是0。