~ubuntu-branches/debian/sid/lammps/sid

« back to all changes in this revision

Viewing changes to examples/min/in.min.box

  • Committer: Package Import Robot
  • Author(s): Anton Gladky
  • Date: 2015-04-29 23:44:49 UTC
  • mfrom: (5.1.3 experimental)
  • Revision ID: package-import@ubuntu.com-20150429234449-mbhy9utku6hp6oq8
Tags: 0~20150313.gitfa668e1-1
Upload into unstable.

Show diffs side-by-side

added added

removed removed

Lines of Context:
1
 
# 2d Lennard-Jones melt and subsequent energy minimization
2
 
 
 
1
# 2d Lennard-Jones melt and subsequent energy minimization,
 
2
# followed by box relaxation to a target stress state
 
3
#
3
4
units           lj
4
5
dimension       2
5
6
atom_style      atomic
36
37
 
37
38
run             1000
38
39
 
 
40
# perform the box relaxation by progressively
 
41
# adding more degrees of freedom, for smoother
 
42
# convergence
 
43
 
39
44
neigh_modify    delay 0 every 1 check yes
40
45
 
41
46
velocity        all create 0.0 1
42
47
thermo          50
 
48
 
 
49
# atoms only
 
50
 
43
51
minimize        1.0e-6 0.001 1000 10000
44
52
 
45
 
fix             3 all box/relax x 1.0 y 2.0 vmax 1.0e-4 nreset 100
46
 
 
47
 
thermo_style    custom step temp pe pxx pyy pxy
48
 
 
49
 
min_modify      line quadratic
50
 
minimize        0.0 1.0e-6 10000 100000
 
53
# isotropic volume relaxation to hydrostatic target
 
54
 
 
55
fix             3 all box/relax iso 1.5 vmax 1.0e-3
 
56
 
 
57
variable        emin equal pe+f_3/atoms
 
58
thermo_style    custom step temp pe pxx pyy pxy f_3 v_emin
 
59
 
 
60
# anisotropic volume relaxation to hydrostatic target
 
61
 
 
62
minimize        0.0 10.0e0 10000 100000
 
63
 
 
64
fix             3 all box/relax aniso 1.5 vmax 1.0e-3
 
65
 
 
66
minimize        0.0 10.0 10000 100000
 
67
 
 
68
# anisotropic volume relaxation to non-hydrostatic target
 
69
 
 
70
fix             3 all box/relax x 1.0 y 2.0 vmax 1.0e-3
 
71
 
 
72
minimize        0.0 10.0 10000 100000
 
73
 
 
74
minimize        0.0 1.0 10000 100000
 
75
 
 
76
minimize        0.0 1.0e-2 10000 100000
 
77
 
 
78
# Final refinement uses nreset to eliminate stress offset
 
79
 
 
80
fix             3 all box/relax x 1.0 y 2.0 vmax 1.0e-3 nreset 100
 
81
 
 
82
minimize        0.0 1.0e-4 10000 100000
 
83