LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Mailing List Archives
[lammps-users] Movement of atom layer
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[lammps-users] Movement of atom layer


From: "Jeong, Jae Young" <JaeyoungJeong@...1598...>
Date: Mon, 30 Apr 2018 22:33:47 +0000

Hi, LAMMPS users

I have a question about fixing the initial position of atoms. I'm calculating the velocity autocorrelation of graphene on Cu with different distance between two materials. When I set 40angstrom for the distance and ran LAMMPS. Graphene layer moved to the Cu surface while calculating. In order word, graphene layer was at the initial position but it closed the Cu when the calculation started as attached files. How can I fix graphene layer above Cu surface with 40angstrong distance?

Here is my input script. I'd really appreciate it if you give me comments.

Thanks,

Jaeyoung

#  This is the control script for LAMMPS
echo both
#=======================================#
# Input file for phonon density of state
#=======================================#
units           metal
atom_style      atomic
boundary        p p p

read_data       g40-cu.data

group           C  type 1
group           Cu  type 2

#=======================================#
# Define interatomic potential
#=======================================#
pair_style       hybrid eam lj/cut 9.0 tersoff
pair_coeff       * * eam Cu_u3.eam
pair_coeff       * * tersoff tersoff.coeffs C NULL
pair_coeff       1 2 lj/cut 0.02578 3.6

neighbor        2.0 bin
neigh_modify    delay 0 every 1 check yes

#===========================================#
# Minimization
#============================================#

thermo 50
thermo_style custom step temp pe ke etotal press vol
min_style       cg
minimize        1.0e-8 1.0e-8 5000 10000

#==================================================#
# Relaxation of the lattice to the potential in NPT
#==================================================#
velocity        all create 300.00 156467 rot yes mom yes dist gaussian
timestep 0.0005
fix             1 all npt temp 300.00 300.00 0.02 iso 0.0 0.0 0.2
reset_timestep 0
run 50000
unfix           1

#================================#
# Equilibriate the system in NVE #
#================================#
reset_timestep  0
fix             2 all nve
fix             3 all momentum 1000 linear 1 1 1 
dump DUMP_NVE all xyz 100000 dump_NVE.xyz
dump_modify DUMP_NVE element C Cu
variable        tt0 equal step
run             500000
unfix           3  

#================================================================#
# Calculation of properties velocity autocorelation  in NVE #
#================================================================#
variable        tt1 equal ${tt0}
compute 1 C  vacf
compute 2 Cu vacf
variable        Cxxvacf equal c_1[1]
variable        Cyyvacf equal c_1[2]
variable        Czzvacf equal c_1[3]
variable        Callvacf equal c_1[4]
variable        Cuxxvacf equal c_2[1]
variable        Cuyyvacf equal c_2[2]
variable        Cuzzvacf equal c_2[3]
variable        Cuallvacf equal c_2[4]
variable        tt equal (step)*0.0005

reset_timestep 0
thermo          10
thermo_style    custom step c_1[1] c_1[2] c_1[3] c_1[4] c_2[1] c_2[2] c_2[3] c_2[4]
fix             vacfresult all print 10 "${tt} ${Cxxvacf} ${Cyyvacf} ${Czzvacf} ${Callvacf} ${Cuxxvacf} ${Cuyyvacf} ${Cuzzvacf} ${Cuallvacf}" screen no file vacf.txt title "# Velocity Autocorrelation"

run             1000000

write_restart   graphene.restart



Attachment: inintial.png
Description: inintial.png

Attachment: calculating.png
Description: calculating.png