forked from lijiext/lammps
119 lines
4.0 KiB
Groff
119 lines
4.0 KiB
Groff
LAMMPS (1 Feb 2014)
|
|
# 2d Lennard-Jones melt and subsequent energy minimization
|
|
|
|
units lj
|
|
dimension 2
|
|
atom_style atomic
|
|
|
|
lattice sq2 0.8442
|
|
Lattice spacing in x,y,z = 1.53919 1.53919 1.53919
|
|
region box block 0 20 0 20 -0.1 0.1
|
|
create_box 1 box
|
|
Created orthogonal box = (0 0 -0.153919) to (30.7838 30.7838 0.153919)
|
|
2 by 2 by 1 MPI processor grid
|
|
create_atoms 1 box
|
|
Created 800 atoms
|
|
mass 1 1.0
|
|
|
|
velocity all create 5.0 87287 loop geom
|
|
|
|
pair_style lj/cut 2.5
|
|
pair_coeff 1 1 1.0 1.0 2.5
|
|
pair_modify shift yes
|
|
|
|
neighbor 0.3 bin
|
|
neigh_modify delay 0 every 1 check yes
|
|
|
|
fix 1 all nve
|
|
fix 2 all enforce2d
|
|
|
|
#dump 1 all atom 100 dump.min
|
|
|
|
#dump 2 all image 50 image.*.jpg type type # zoom 1.6 adiam 1.5
|
|
#dump_modify 2 pad 4
|
|
|
|
#dump 3 all movie 50 movie.mpg type type # zoom 1.6 adiam 1.5
|
|
#dump_modify 3 pad 4
|
|
|
|
thermo 100
|
|
|
|
run 1000
|
|
Memory usage per processor = 2.06115 Mbytes
|
|
Step Temp E_pair E_mol TotEng Press
|
|
0 5 -2.461717 0 2.532033 5.0190509
|
|
100 3.2788864 -0.74311698 0 2.5316708 15.912832
|
|
200 3.3041082 -0.77000277 0 2.5299752 15.602653
|
|
300 3.3840228 -0.84859211 0 2.5312006 15.188203
|
|
400 3.3916063 -0.85694601 0 2.5304208 15.383853
|
|
500 3.3136052 -0.77935264 0 2.5301106 15.833296
|
|
600 3.3888915 -0.85213742 0 2.532518 15.162759
|
|
700 3.2123632 -0.67641807 0 2.5319297 16.310592
|
|
800 3.3016497 -0.76571656 0 2.531806 15.639198
|
|
900 3.4621454 -0.9260857 0 2.531732 14.772965
|
|
1000 3.2803075 -0.74465936 0 2.5315477 15.835597
|
|
Loop time of 0.066597 on 4 procs for 1000 steps with 800 atoms
|
|
|
|
Pair time (%) = 0.0324514 (48.728)
|
|
Neigh time (%) = 0.012701 (19.0714)
|
|
Comm time (%) = 0.0137331 (20.6212)
|
|
Outpt time (%) = 0.000190914 (0.28667)
|
|
Other time (%) = 0.00752062 (11.2927)
|
|
|
|
Nlocal: 200 ave 202 max 198 min
|
|
Histogram: 1 0 1 0 0 0 0 1 0 1
|
|
Nghost: 170.25 ave 176 max 167 min
|
|
Histogram: 1 0 2 0 0 0 0 0 0 1
|
|
Neighs: 1936.75 ave 1963 max 1896 min
|
|
Histogram: 1 0 0 0 0 1 0 0 1 1
|
|
|
|
Total # of neighbors = 7747
|
|
Ave neighs/atom = 9.68375
|
|
Neighbor list builds = 203
|
|
Dangerous builds = 0
|
|
|
|
neigh_modify delay 0 every 1 check yes
|
|
|
|
#dump_modify 1 every 25
|
|
|
|
thermo 50
|
|
minimize 1.0e-6 0.001 1000 10000
|
|
Memory usage per processor = 2.7478 Mbytes
|
|
Step Temp E_pair E_mol TotEng Press
|
|
1000 3.2803075 -0.74465936 0 2.5315477 15.835597
|
|
1050 3.2803075 -2.8198371 0 0.45636995 1.6553949
|
|
1100 3.2803075 -2.8554989 0 0.42070813 1.6571551
|
|
1150 3.2803075 -2.8646695 0 0.41153758 1.7045388
|
|
1200 3.2803075 -2.8688213 0 0.40738582 1.6922085
|
|
1250 3.2803075 -2.8764988 0 0.39970826 1.6948532
|
|
1300 3.2803075 -2.8776606 0 0.39854651 1.6821387
|
|
1350 3.2803075 -2.8789379 0 0.39726915 1.6715024
|
|
1365 3.2803075 -2.879076 0 0.3971311 1.66923
|
|
Loop time of 0.100178 on 4 procs for 365 steps with 800 atoms
|
|
|
|
Minimization stats:
|
|
Stopping criterion = energy tolerance
|
|
Energy initial, next-to-last, final =
|
|
-0.744659361981 -2.87907436106 -2.8790759758
|
|
Force two-norm initial, final = 2006.88 0.477401
|
|
Force max component initial, final = 395.193 0.0879763
|
|
Final line search alpha, max atom move = 0.00334856 0.000294593
|
|
Iterations, force evaluations = 365 1437
|
|
|
|
Pair time (%) = 0.0639488 (63.8349)
|
|
Neigh time (%) = 0.00472444 (4.71603)
|
|
Comm time (%) = 0.015418 (15.3905)
|
|
Outpt time (%) = 0.000124454 (0.124233)
|
|
Other time (%) = 0.0159628 (15.9343)
|
|
|
|
Nlocal: 200 ave 204 max 194 min
|
|
Histogram: 1 0 0 0 0 1 0 0 0 2
|
|
Nghost: 171.75 ave 173 max 170 min
|
|
Histogram: 1 0 0 1 0 0 0 0 0 2
|
|
Neighs: 1762.75 ave 1797 max 1711 min
|
|
Histogram: 1 0 0 0 0 0 1 1 0 1
|
|
|
|
Total # of neighbors = 7051
|
|
Ave neighs/atom = 8.81375
|
|
Neighbor list builds = 78
|
|
Dangerous builds = 0
|