LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Mailing List Archives
[lammps-users] units metal vs real -Reaxff

# [lammps-users] units metal vs real -Reaxff

 From: Alexandra Davila Date: Wed, 18 Oct 2017 15:46:16 +0200

```Hi lammps-users,

```
I have been calculating with reax force field some gold systems, I have always used units real.
```
```
By error I did a calculation of the lattice with units metal, I was wondering why the output of the potential energy for both runs (with real and metal units) are the same. For real units energy should be printed in kcal/mol and for metal in eV, or with reax the default is real units (I think not)?
```
Probably I am making a very dummy mistake!

I attach the *in file!

thanks for the help!

lammps version:  Stable version (11 Aug 2017)

```
```# ---------- Initialize Simulation ---------------------
clear
units real
dimension 3
boundary p p p
#atom_style atomic
#atom_modify map array
atom_style charge

# ---------- Create Atoms ---------------------
lattice         fcc 4.080
region  box block 0 1 0 1 0 1 units lattice
create_box      1 box

lattice fcc 4.080 orient x 1 0 0 orient y 0 1 0 orient z 0 0 1

create_atoms 1 box
replicate 1 1 1

mass 1 196.966
set type 1 charge 0.0
# ---------- Define Interatomic Potential ---------------------
pair_style reax/c NULL
pair_coeff * * ffield.reax.AuO Au
#neighbor 2.0 bin
#neigh_modify delay 10 check yes

# ---------- Define Settings ---------------------
compute eng all pe/atom
compute eatoms all reduce sum c_eng

# ---------- Run Minimization ---------------------
#reset_timestep 0
fix 1 all qeq/reax 1 0.0 10.0 1e-6 reax/c
run 0

variable natoms equal "count(all)"
variable teng equal "pe"
variable length equal "lx"
variable ecoh equal "v_teng/v_natoms"

print "Total energy (eV) = \${teng};"
print "Number of atoms = \${natoms};"
print "Lattice constant (Angstoms) = \${length};"
print "Cohesive energy (eV) = \${ecoh};"

print "All done!"
```