Hi
I am trying to simulation rigid sphere (C180 Buckyball) impact on the
graphene sheet using NVE ensemble with AIREBO potential.
I have relaxed the graphene and the sphere for 10 ps before impact.
According to the impact mechanics, during impact kinetic energy (KE) of the
sphere will transfer into potential energy (PE) and kinetic energy
of the graphene and total energy (PE + KE) of the sphere and graphene should
be conserved. In the simulation,
I am observing before and after impact system total energy is conserved.
However, during impact (around 1.6 ps) total energy is not conserved.
Total energy is changing from -113000 eV to - 112564 eV during the impact.
I have also attached the Fig as ppt. I am sure whether you
can see it.
I have also tried initial relaxation with NVT ensemble. No significant
change.
In another simulation, I have checked the only the rigid body motion of the
C180 Buckyball at 5 km/s with "rigid/nve single" fix command.
This simulation result matches with the theoretical prediction of KE,
displacement etc.
I am using 2012 version of lammps. I have attached input script, and data
file.
Would you please make comments regarding why energy is not conserved during
impact? How can energy be conserved during impact?
Thanks,
Kind regards,
Sanjib
University of Delaware
# input script for Sphere Impact on Graphene
#-------10--------20--------30--------40--------50--------60--------70------
--80
# Variables
variable input index input_1-20-20-gs-1-C180.txt
# INITIALIZATION
units metal
boundary p p p
atom_style atomic
processors 6 6 1
# ATOM DEFINITION
read_data ${input}
# Define regions
region left block INF -90.0 INF INF INF INF units box
region right block 90.0 INF INF INF INF INF units box
region front block INF INF INF -90.5 INF INF units box
region rear block INF INF 90.5 INF INF INF units box
region section block INF INF 0.00 INF INF INF units box
group left region left
group right region right
group front region front
group rear region rear
group section region section
group gs type 1
group bball type 2
group center subtract gs left right front rear
# FORCE FIELD
pair_style airebo 3.0 1 1
pair_coeff * * CH.airebo C C
# SETTINGS
neighbor 1.0 bin
neigh_modify delay 5
# ------------------------------ RELAXATION
GS-Impact-Results.pptx (182 KB)