Hello, Xiangxiong:
Please keep the list cc’d on future responses.
The other three components are the off-diagonal components of the pressure tensor. You can see the specific definitions for all six components of virial[] in the ev_tally routines that are found in pair.cpp.
–AEI