Thank you Aidan. But what to do when I need to maintain zero pressure in other directions?? Although I tried with nvt ensemble and it is giving fine results. But in case of any pure mode (fracture) loading pressure in x direction has to be zero. I think nve can’t be used there.