求助：DFPT计算时固定了部分原子，如何提取力常数矩阵

正文摘要:

体系是含金属离子的水溶液，目前想要计算该金属离子的平均力常数，因此计算时只考虑了金属离子及其周围的一圈水分子。在利用vasp优化结构和计算声子频率时固定了其它原子，之后再利用phonopy --fc vasprun.xml提取力 ...