Dear Lammps/MD Expert,
I am going to calculate the enthalp of deformation(H=U+PV) inside my 2 dimensional model with LAMMPS. for this purpose:
1 first I calculate the internal energy of the region (not all the system):
variable etotal_region equal “(c_pe_region+c_ke_region)”
2 Then I calculate the pressure inside the 2 dimensional region with this command:
 variable press_region equal (c_p[1]+c_p[2])/(2*Area of the region)
3 then I calculate the enthalpy in my NVE ensemble with :
variable enthalpy equal "(v_etotal+ Area of the region *v_press_region)
and I reached to some results.
Howeve, I have no idea what should be the order of the enthalpy values actually!!! I guess they should be with the same order with internal energy.
Best regards,
Hi Bahman,
What you are doing is conceptually straight forward, but it is not a typical LAMMPS calculation, and so there are lots of extra things you have to get right. I suggest eliminating some of the complicating factors first, reproducing a known correct result, and then adding back in the complicating factors one by one.

A related LAMMPS calculation would be to calculate total internal energy as function of compression for a 3D structure. This could be done statically at zero temperature, so no kinetic contribution and no configurational sampling, or it could be done at finite temperature, with a kinetic contribution, and timeaveraging of the configurational distribution. So, first figure out how to do that calculation correctly, comparing to prior results.

Restrict the calculation to a finite region. This is a nontrivial extension. For example, the c_p[1] that you use below must be defined for a region, but the LAMMPS compute pressure always operates on the group all (From the doc page 'The specified group must be “all”. '). You will have to use stress/atom or else make some kind of assumption about the relation between pressure in region and global pressure.

Similar considerations apply to defining kinetic energy on a region.

Restrict to 2D
Since you say that you have no idea what values to expect, I strongly suggest you spend some time reading the prior literature, to get some idea.
Aidan