Dear Eduardo,<br><br>thank you very much for expanded answer and sharing the practical tricks.<br>I've done some computational experiments on bulk Si (cubic conventional cell) vc-relaxation.<br><br>Here is the result (V is volume of initial unit cell, and V0 is equilibrium volume):<br>
<br>1) starting from V > V0, i.e. 1/V < 1/V0 -> more G-vectors for vc-relax:<br><br>G cutoff = 837.7995 ( 101505 G-vectors) FFT grid: ( 60, 60, 60) - vc-relax<br>G cutoff = 837.7995 ( 97137 G-vectors) FFT grid: ( 60, 60, 60) - post-scf<br>
<br>! total energy = -372.89634728 Ry - the last energy in the course of vc-relax<br>! total energy = -372.89587589 Ry - post-scf energy<br><br>NB1: # of G-vectors (vc-relax) > # G-vectors(post-scf), and E(the last point vc-relax) < E(post-scf). <br>
<br>1) starting from V < V0, i.e. 1/V > 1/V0 -> more G-vectors for post-scf:<br><br>G cutoff = 775.1830 ( 90447 G-vectors) FFT grid: ( 60, 60, 60) - vc-relax<br>G cutoff = 775.1830 ( 97137 G-vectors) FFT grid: ( 60, 60, 60) - post-scf<br>
<br>! total energy = -372.89498529 Ry - the last energy in the course of vc-relax<br>! total energy = -372.89587142 Ry - post-scf energy<br><br>NB2: # of G-vectors(vc-relax) < # G-vectors(post-scf), and E(the last point vc-relax) > E(post-scf). <br>
<br>Comparing these two experiments, one can make a preliminary conclusion: the more G-vectors, the lower<br>the total Energy, provided all other parameters to be fixed.<br>This is easy to understand: plane-wave basis set is complete, that means 2 things (when dealing with truncated bases):<br>
1) E(N+M) < E(N), where N,M - number of plane waves(G-vectors);<br>2) lim N->infinity of [ E(N+M)-E(N)] = 0.<br><br>Now it seems to be more clear for me :)<br>Correct me if I'm wrong somewhere.<br><br>-- <br>Best regards, Max Popov<br>
Ph.D. student<br>Materials center Leoben (MCL), Leoben, Austria.<br><br><div class="gmail_quote">2011/4/19 Eduardo Ariel Menendez Proupin <span dir="ltr"><<a href="mailto:eariel99@gmail.com" target="_blank">eariel99@gmail.com</a>></span><br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div><span style="font-family: arial,sans-serif; font-size: 13px; border-collapse: collapse;"><div>
>Dear Dr. Giannozzi,<br>>thank you for the answer! I could find it myself looking in the output file a bit >more carefully...<br>
>One thing, which is somehow contrary to my expectations, is that the final scf >energy is higher<br>>than the last one from vc-relax. Could you, please, elaborate a bit on the >matter?<br><br></div>Dear Maxim, <br>
</span></div><div><span style="font-family: arial,sans-serif; font-size: 13px; border-collapse: collapse;"><br></span></div><div><span style="font-family: arial,sans-serif; font-size: 13px; border-collapse: collapse;">I followed this discussion with interest, and thanks to that I learned about the new scf calculation with final G-vectors. Concerning your last question, </span></div>
<div><span style="font-family: arial,sans-serif; font-size: 13px; border-collapse: collapse;">the energy is higher because the vc-relaxed energy was optimized for a different basis set, than the final scf calculation (different G-vectors). Hence, the energy of the final scf calculation is is made for a structure that is slightly out of the minimum for the new basis set. Remember than the G-vectors used in a scf calculation are all the reciprocal lattice vectors contained in a sphere that has a radius determined by the cutoff. These vectors are selected at the first step of the vc-relaxation. When the unit cell gets deformed, the G-vectors vary accordingly, and the region that the G-vectors occupy is a deformation from the initial sphere, maybe an elipsoid. When the vc-relax stops, the final scf calculation takes the G-vectors contained inside a sphere. Hence, some of the old G-vectors that were in the border of the deformed sphere may be eliminated, and some that were absent are now included. </span></div>
<div><font face="arial, sans-serif"><span style="border-collapse: collapse;">If you had used an (impossible) infinite cutoff, the basis set would be complete in both cases (G-vectors contained in an infinite sphere or in an infinite elipsoid) and there would be no difference. Usually, I repeat the vc-relax procedure starting 'from_scratch' with the last structure (coordinates and lattice vectors) in the new input file, until the vc-relax procedure performs only one step. In this case there is no difference. If it never happens that vc-relax stops at the first step, then I increase the cutoffs. In your case, the energy difference of 0.5 mRy may be small enough and do not need to do that. It depends on the property that you want. E.g., if you are interested in elastic properties, you may need that the minimal energy structure also gives a stress tensor below 0.1 kbar or so. If you cannot get it, increase the cutoff. </span></font></div>
<div><font face="arial, sans-serif"><span style="border-collapse: collapse;"><br></span></font></div><div><font face="arial, sans-serif"><span style="border-collapse: collapse;">The following link may help</span></font></div>
<br clear="all"><a href="http://www.quantum-espresso.org/wiki/index.php/Methodological_Background#Stress" target="_blank">http://www.quantum-espresso.org/wiki/index.php/Methodological_Background#Stress</a><br><br><div>Best regards</div>
<div>
<br></div><div>-- <br><div><br></div>
<div><br></div>Eduardo Menendez<br>Departamento de Fisica<br>Facultad de Ciencias<br>Universidad de Chile<br>Phone: (56)(2)9787439<br>URL: <a href="http://fisica.ciencias.uchile.cl/%7Eemenendez" target="_blank">http://fisica.ciencias.uchile.cl/~emenendez</a><br>
</div>
<br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a><br>
<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
<br></blockquote></div>