Hello
I am doing simulation of CNT. I am trying to do a simulation in which a single carbon atom with a certain velocity collides with an atom of CNT located at the center of the tube(length). However, the single atom does not move from its initial position irrespective of the velocity values.
The code I use:
units metal
boundary s p p
lattice hcp 1
atom_style atomic
read_data read10.cnt
pair_style airebo 2.5 1 1
pair_coeff * * CH.airebo C C C
group mobile type == 1
group fixed type == 2
group atom_pt type == 3
# Equilibrate
fix 1 mobile nvt 300.0 300.0 0.1
velocity fixed set 0.0 0.0 0.0
timestep 0.004
run 100
# Collision
velocity atom_pt set 0.0 -1.0 0.0
dump mydump all atom 100 dump.cnt
run 1000
Could anyone please suggest what could be wrong?
Regards
APC
|
read10.cnt (12.4 KB)