• PUBLIC 公開
Dear Lammps users,
I want to simulate the stretching of a mwcnt along its longitudinal direction (z).
I have tried with the “change_box” and the “fix deform” commands, both within a wrapping for loop.
In both cases, before getting the box to stretch along the z direction, it squeezes in the xy plane, which makes my cnt look like a square because of periodic images come too close on +/- x/y directions.
How can I solve this? Do I have to say something about not conserving the volume?
The system has been previously equilibrated and relaxed.
I perform a npt dynamics with the following fix:
"fix mynpt all npt temp 300 300 300 iso 1 1 1 "
The full input.file is given below.
Thank you for your support.
BR
Marco
BOF
variable T equal 300
units metal
atom_style full
boundary p p p
read_restart relax_restart.mdrun
pair_style airebo/morse 3.0 1 1
pair_coeff * * CH.airebo-m C
reset_timestep 0
thermo 100
Elongate cell
dump dump_elongate all atom 100 dump_elongate.lammpstrj
variable var loop 30
label loop
fix mynpt all npt temp {T} {T} ${T} iso 1 1 1
run 1000
unfix mynpt
fix 1 all deform 1 z scale 1.1 # to be removed for change_box case
run 1 #
unfix 1 #
#to be inserted for change_box_case
#variable fac equal 1+0.01*{var} #change_box all z scale {fac} remap
next var
jump input.lammps loop
undump dump_elongate
EOF
This e-mail may contain confidential information. If you are not an addressee or otherwise authorised to receive this message, you should not use, copy, disclose or take any action based on this e-mail. If you have received this e-mail in error, please inform the sender promptly and delete this message and any attachments immediately.