1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
|
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2 Final//EN">
<!--Converted with LaTeX2HTML 2019.2 (Released June 5, 2019) -->
<HTML lang="EN">
<HEAD>
<TITLE>3.4 Optimization and dynamics</TITLE>
<META NAME="description" CONTENT="3.4 Optimization and dynamics">
<META NAME="keywords" CONTENT="user_guide">
<META NAME="resource-type" CONTENT="document">
<META NAME="distribution" CONTENT="global">
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=utf-8">
<META NAME="viewport" CONTENT="width=device-width, initial-scale=1.0">
<META NAME="Generator" CONTENT="LaTeX2HTML v2019.2">
<LINK REL="STYLESHEET" HREF="user_guide.css">
<LINK REL="next" HREF="node12.html">
<LINK REL="previous" HREF="node10.html">
<LINK REL="next" HREF="node12.html">
</HEAD>
<BODY >
<!--Navigation Panel-->
<A
HREF="node12.html">
<IMG WIDTH="37" HEIGHT="24" ALT="next" SRC="next.png"></A>
<A
HREF="node7.html">
<IMG WIDTH="26" HEIGHT="24" ALT="up" SRC="up.png"></A>
<A
HREF="node10.html">
<IMG WIDTH="63" HEIGHT="24" ALT="previous" SRC="prev.png"></A>
<A ID="tex2html164"
HREF="node1.html">
<IMG WIDTH="65" HEIGHT="24" ALT="contents" SRC="contents.png"></A>
<BR>
<B> Next:</B> <A
HREF="node12.html">3.5 Direct interface with</A>
<B> Up:</B> <A
HREF="node7.html">3 Using PWscf</A>
<B> Previous:</B> <A
HREF="node10.html">3.3 Electronic structure calculations</A>
<B> <A ID="tex2html165"
HREF="node1.html">Contents</A></B>
<BR>
<BR>
<!--End of Navigation Panel-->
<!--Table of Child-Links-->
<A ID="CHILD_LINKS"><STRONG>Subsections</STRONG></A>
<UL>
<LI><UL>
<LI><A ID="tex2html166"
HREF="node11.html#SECTION00044010000000000000">3.4.0.1 Structural optimization</A>
<LI><A ID="tex2html167"
HREF="node11.html#SECTION00044020000000000000">3.4.0.2 Molecular Dynamics</A>
<LI><A ID="tex2html168"
HREF="node11.html#SECTION00044030000000000000">3.4.0.3 Free-energy surface calculations</A>
<LI><A ID="tex2html169"
HREF="node11.html#SECTION00044040000000000000">3.4.0.4 Variable-cell optimization</A>
<LI><A ID="tex2html170"
HREF="node11.html#SECTION00044050000000000000">3.4.0.5 Variable-cell molecular dynamics</A>
</UL></UL>
<!--End of Table of Child-Links-->
<HR>
<H2><A ID="SECTION00044000000000000000">
3.4 Optimization and dynamics</A>
</H2>
<P>
<H4><A ID="SECTION00044010000000000000">
3.4.0.1 Structural optimization</A>
</H4>
For fixed-cell optimization, specify <TT>calculation='relax'</TT> and
add namelist &IONS. All options for a single SCF calculation apply,
plus a few others. You may follow a structural optimization with a
non-SCF band-structure calculation. See example 2.
<P>
<H4><A ID="SECTION00044020000000000000">
3.4.0.2 Molecular Dynamics</A>
</H4>
Specify <TT>calculation='md'</TT>, the time step <TT>dt</TT>,
and possibly the number of MD stops <TT>nstep</TT>.
Use variable <TT>ion_dynamics</TT> in namelist &IONS
for a fine-grained control of the kind of dynamics.
Other options for setting the initial temperature
and for thermalization using velocity rescaling are available.
Remember: this is MD on the electronic ground state, not
Car-Parrinello MD.
See example 3.
<P>
<H4><A ID="SECTION00044030000000000000">
3.4.0.3 Free-energy surface calculations</A>
</H4>
<P>
Once <TT>PWscf</TT> is patched with the <TT>PLUMED</TT> plug-in
(only v.1.3), it is possible to
use most PLUMED functionalities by running <TT>PWscf</TT> as:
<TT>./pw.x -plumed</TT> plus the other usual <TT>PWscf</TT> arguments.
The input file for <TT>PLUMED</TT> must be found in the specified
<TT>outdir</TT> with fixed name <TT>plumed.dat</TT>.
<P>
<H4><A ID="SECTION00044040000000000000">
3.4.0.4 Variable-cell optimization</A>
</H4>
<P>
Variable-cell calculations (both optimization and dynamics) are performed
with plane waves and G-vectors <EM>calculated for the starting cell</EM>.
Inly the last step, after convergence has been achieved, is performed
for the converged structure, with plane waves and G-vectors
<EM>calculated for the final cell</EM>. Small differences between the
two last steps are thus to be expected and give an estimate of the
convergence of the variable-cell optimization with respect to the
plane-wave basis. A large difference means that you are far from
convergence in the plane-wave basis set and you need to increase the
cutoff(s) <TT>ecutwfc</TT> and/or (if applicable) <TT>ecutrho</TT>.
<P>
<H4><A ID="SECTION00044050000000000000">
3.4.0.5 Variable-cell molecular dynamics</A>
</H4>
<P>
"A common mistake many new users make is to set the time step <TT>dt</TT>
improperly to the same order of magnitude as for CP algorithm, or
not setting <TT>dt</TT> at all. This will produce a ``not evolving dynamics''.
Good values for the original RMW (Wentzcovitch) dynamics are
<TT>dt</TT> <!-- MATH
$= 50 \div 70$
-->
= 50÷70. The choice of the cell mass is a
delicate matter. An
off-optimal mass will make convergence slower. Too small masses, as
well as too long time steps, can make the algorithm unstable. A good
cell mass will make the oscillation times for internal degrees of
freedom comparable to cell degrees of freedom in non-damped
Variable-Cell MD. Test calculations are advisable before extensive
calculation. I have tested the damping algorithm that I have developed
and it has worked well so far. It allows for a much longer time step
(dt=<!-- MATH
$100 \div 150$
-->
100÷150) than the RMW one and is much more stable with very
small cell masses, which is useful when the cell shape, not the
internal degrees of freedom, is far out of equilibrium. It also
converges in a smaller number of steps than RMW." (Info from Cesar Da
Silva: the new damping algorithm is the default since v. 3.1).
<P>
<HR>
<!--Navigation Panel-->
<A
HREF="node12.html">
<IMG WIDTH="37" HEIGHT="24" ALT="next" SRC="next.png"></A>
<A
HREF="node7.html">
<IMG WIDTH="26" HEIGHT="24" ALT="up" SRC="up.png"></A>
<A
HREF="node10.html">
<IMG WIDTH="63" HEIGHT="24" ALT="previous" SRC="prev.png"></A>
<A ID="tex2html164"
HREF="node1.html">
<IMG WIDTH="65" HEIGHT="24" ALT="contents" SRC="contents.png"></A>
<BR>
<B> Next:</B> <A
HREF="node12.html">3.5 Direct interface with</A>
<B> Up:</B> <A
HREF="node7.html">3 Using PWscf</A>
<B> Previous:</B> <A
HREF="node10.html">3.3 Electronic structure calculations</A>
<B> <A ID="tex2html165"
HREF="node1.html">Contents</A></B>
<!--End of Navigation Panel-->
</BODY>
</HTML>
|
