Hello everyone,
I am trying to simulate protein in a solvated environment with parallel strands of protein molecules. I am trying to compute the contribution of non-bonded interaction energy between proteins to the total potential energy of the system. Initially, I tried with "pe/atom" using the keywords pair and kspace, but according to the documentation provided it says "energy of each atom is due to its interaction with all other atoms in the simulation, not just with other atoms in the group." So the energy output of the syntax will be a combination of protein-protein and protein-water interactions.
Then I tried with group/group syntax which gives the total energy of interaction. So to obtain potential energy contribution, I subtracted the Kinetic energy of the protein atoms from the output. The output was coming quite absurdly high.
I would like to know can group/group command give the required output using the two groups as "protein" subtracting kinetic energy of the protein atoms. Or is there any other way to compute the non-bonded interactions between protein atoms?