GULP采用经典势能模型输出shengBTE计算用力常数文件

仅需在option参数中使用supercell控制2阶力常数的扩胞倍数,以及output shengBTE:

但是注意,gulp这里存在一个bug,输出的2阶力常数文件中超胞中原子次序与phonopy以及shengBTE用不一致:shengBTE中FORCE_CONSTANT中原子次序为:ucatom*uz*uy*ux;但是gulp中为ucatom*ux*uy*uz,对于各向同性系统,则计算没有区别,但是对于各向异性系统,会产生计算错误,需要转换2阶力常数的中的原子次序。

你可能感兴趣的:(gulp,node.js,前端)