I want to calculate surafce stress in Si (001) surface after reconstruction by using software quantum espresso. As in the output it is displaying stress matrix ..my z - component (sigma zz) of stress is not zero, but I have taken a enough vaccum , no external force acts on the surafce and therfore it should be zero .So I am not understanding why it is giving z - component (sigma zz) of stress ?