Syntax:
compute ID group-ID temp/region region-ID
Examples:
temperature mine flow region boundary
Description:
Define a computation that calculates the temperature of a group of atoms in a geometric region. This can be useful for thermostatting one portion of the simulation box. E.g. a McDLT simulation where one side is cooled, and the other side is heated. A compute of this style can be used by any command that computes a temperature, e.g. thermo_modify, fix temp/rescale, etc.
Note that a region-style temperature can be used to thermostat with fix temp/rescale or fix langevin, but should probably not be used with Nose/Hoover style fixes (>fix nvt, fix npt, or fix nph), if the degrees-of-freedom included in the computed T varies with time.
The temperature is calculated by the formula KE = dim/2 N k T, where KE = total kinetic energy of the group of atoms (sum of 1/2 m v^2), dim = 2 or 3 = dimensionality of the simulation, N = number of atoms in both the group and region, k = Boltzmann constant, and T = temperature.
A 6-component kinetic energy tensor is also calculated by this compute for use in the computation of a pressure tensor. The formula for the components of the tensor is the same as the above formula, except that v^2 is replaced by vx * vy for the xy component, etc.
The number of atoms contributing to the temperature is compute each time the temperature is evaluated since it is assumed atoms can enter/leave the region. Thus there is no need to use the dynamic option of the compute_modify command for this compute style.
Unlike other compute styles that calculate temperature, this compute does NOT currently subtract out degrees-of-freedom due to fixes that constrain molecular motion, such as fix shake and fix rigid. If needed the subtracted degrees-of-freedom can be altered using the extra option of the compute_modify command.
Restrictions: none
Related commands:
compute temp, compute pressure
Default: none