I'm dealing with a mixture of water and hydrocarbons and want to calculate pair-wise energy and forces between water molecules as below:
compute 1 water group/group water pair yes
however, the thermo output prints zero force components!! but energy has a finite value.
thermo_style step c_1 c_1 c_1 c_1
Why does go wrong in my work?