Dear all,
I made a monoatomic gas simulation using NVE ensemble. But when I run the code, it outputs zero atoms until step 30 and then goes directly to 1000 atoms as I have defined on the input file.
Why does it happens?
Thanks
# Monoatomic gas - Argon
# PART A - ENERGY MINIMIZATION
# 1) Initialization
units lj
dimension 3
atom_style atomic
pair_style lj/cut 2.5
boundary p p p #no interactions with walls
#sigma = 3.405
#e/kB = 125.2
# 2) System definition
region simulation_box block -450 450 -450 450 -450 450 #dist�ncia entre �tomos = 6.67 sigma
create_box 1 simulation_box
create_atoms 1 random 1000 260702 NULL
group all type 1
# 3) Simulation settings
mass 1 1
pair_coeff 1 1 1.0 1.0 #energy; distance
velocity all create 0.998 196466 dist gaussian #T* = TkB/e
# 4) Visualization (thermodynamic props)
thermo 1
thermo_style custom step atoms temp press pe ke etotal vol
# 5) Energy minimization (stopping criteria)
minimize 1e-6 1e-8 1000 10000 #etol; ftol; maxiter; maxeval
# PART B - MOLECULAR DYNAMICS
# 4) Visualization
thermo 1
dump monogas all atom 100 dump.monogas.nve.lammpstrj
atom_modify sort 0 0.0
# 5) Run
fix simpad all nve #plain time integration to update position and velocity
timestep 0.001
run 10000