Thanks Justin, How about LJ and Coulomb energy? I want to see the trend of protein energy only. I defined 2 energy group: protein and SOL in .mdp file. Then, I saw LJ and Coulomb: protein-ptrotein option when using g_energy. Are they VDW and Coulomb force within the protein only?