Homec4science
Diffusion lammps (lammm-master-deprecated) Browse

lammps/examples/HEAT/README
f876973bcccblammm-master-deprecated

Last Change
Raw File

README
View Options

This directory contains 4 input scripts for carrying out NEMD
simulations of thermal gradients for a Lennard-Jones fluid and SPC/E
water using the HEX/a (fix heat) and eHEX/a (fix ehex) algorithms.
All input scripts are part of the supplementary (open access) material
supporting the publication of Wirnsberger et al. [J. Chem. Phys. 143,
124104 (2015)] and allow one to reproduce the key results reported in
that paper. The full article is available for download under
http://dx.doi.org/10.1063/1.4931597 or http://arxiv.org/pdf/1507.07081
and the supplementary material is available under
https://www.repository.cam.ac.uk/handle/1810/250539.
Each input script in.X produces two output files:
==================================================
-> The file out.TX contains the temperature profile. Plotting coloumn
4 against column 2 yields the results shown in Fig. 5
(Lennard-Jones) or Fig. 6 (SPC/E).
-> The file out.EX contains the time evolution of the total
energy. The difference between the initial and final values of the
total energy in the second column corresponds to a single point in
Fig. 2 (Lennard-Jones) or Fig. 3 (SPC/E).
The input files data.lj and data.spce contain input configurations
which were previously prepared in a non-equilibrium steady-state
(NESS).
For the Lennard-Jones system, the box contains 2000 atoms at a density
of rho*=0.8444 (reduced units). The equilibrium temperature was
T*=0.72 in reduced units. After switching on the thermostat, a
temperature gradient of about -0.02 was generated between the
reservoirs in the reference box. The SPC/E configuration contains 1024
molecules at a density of 0.95 g/cm^3. The equilibrium temperature was
about 400 K, and after switching on the thermostat a temperature
gradient of about -4.9 K/‎Å was generated.
More details can be found in the publication.
c4science · Help