LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Mailing List Archives
Re: [lammps-users] Non-numeric pressure - simulation unstable
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [lammps-users] Non-numeric pressure - simulation unstable

From: Axel Kohlmeyer <akohlmey@...24...>
Date: Thu, 31 Aug 2017 16:34:29 -0400

On Thu, Aug 31, 2017 at 10:59 AM, Jian-Hui Zhai <jianhui.zhai@...4839...> wrote:
Dear All,

I run a simulation and encounter following error: Non-numeric pressure -
simulation unstable. I use FeNiCr.eam.alloy potential and read the model
from a input file.

​this is a frequently discussed error message. please consult the LAMMPS manual and the (many) previous discussions in the LAMMPS mailing​ archives for suggestions of how to address this and report back only in case none of the suggested remedies apply.



The detail information as following:

# ------------------------ INITIALIZATION ----------------------------
units         metal
dimension    3
boundary    p    p    p
atom_style    atomic

# ----------------------- ATOM DEFINITION ----------------------------

read_data "" group all

# ------------------------ FORCE FIELDS ------------------------------

pair_style    eam/alloy
pair_coeff  * * FeNiCr_2013.eam.alloy Fe Ni Cr

# ------------------------- SETTINGS ---------------------------------

compute peratom all pe/atom
compute virial all stress/atom NULL

#---------------------- EQUILIBRATION ----------------------------

# ---------------------set timestep
reset_timestep    0
timestep 0.001

# ---------------set temperature

velocity all create 300 12345

# Assign fixes

fix 1 all npt temp 300 300 0.5 iso 0 0 5 drag 1

# Set thermo output

thermo 400
thermo_style custom step lx ly lz press pxx pyy pzz pe temp

# ------------Use custom dump for Ovito output

# Run for at least 20 picosecond (assuming 1 fs timestep)

run 20000

Jian-Hui ZHAI

Institute of Materials Simulation (WW8)
Department of Materials Science and Engineering
University of Erlangen-Nürnberg
Dr.-Mack-Str. 77, 90762 Fürth, Germany

Check out the vibrant tech community on one of the world's most
engaging tech sites,!
lammps-users mailing list

Dr. Axel Kohlmeyer  akohlmey@...92......
College of Science & Technology, Temple University, Philadelphia PA, USA
International Centre for Theoretical Physics, Trieste. Italy.