声子谱计算过程，首先是体相优化，其次扩胞，然后呢

在完成体相优化和扩胞后，声子谱计算的主要步骤包括：

1. 计算力常数矩阵：使用扩胞后的结构计算各个原子的受力情况，然后计算出力常数矩阵。

2. 计算原子质量和基态振动频率：根据原子的质量和基态振动频率来计算各个原子的热运动，从而在接下来的声子计算中考虑温度效应。

3. 计算动力学矩阵：使用力常数矩阵和原子质量计算动力学矩阵，即描述原子之间相互作用的矩阵。

4. 计算声子本征振动模式：根据动力学矩阵求解声子本征值和本征向量，得到各个振动模式的频率和振幅。

5. 计算声子态密度和热容：使用声子本征值和本征向量计算声子态密度和热容，从而得到材料的热学性质。

这些步骤是声子谱计算的基本流程，不同的计算方法和工具可能会有一些细微的差别。