计算化学公社

标题: 用lammps算热传导系数如何算virial部分的贡献？ [打印本页]

作者
Author: xunxun87 时间: 2016-3-27 00:03
标题: 用lammps算热传导系数如何算virial部分的贡献？
rt. 本人想把kappa matrix分解成KE,PE 和virial部分，已搞定了标量的KE和PE部分，但由于stress 是6个element的vector，不太清楚如何算virial的heat flux vector S*v.

compute myStress all stress/atom virial
variable mySvx atom c_myStress*vx
variable mySvy atom c_myStress*vy
variable mySvz atom c_myStress*vz
compute mySvx all reduce sum v_mySvx
compute mySvy all reduce sum v_mySvy
compute mySvz all reduce sum v_mySvz
fix JJS all ave/correlate $s $p $d c_mySvx c_mySvy c_mySvz type auto file J0Jt.dat ave running
这样应该有问题吧？那么stress的维度该如何修改呢？不胜感激！

欢迎光临计算化学公社 (http://bbs.keinsci.com/)