<div dir="ltr">Dear Dr.Nicola Marzari,<div><br></div><div>Thank You very much for your kind reply. I am changed my system and electron parameters according to your kind suggestions. I have few questions regarding the convergence of the grain boundary. When the grain boundary converged the magnetic moment per site of each atom has become zero and in some cases negative.Does this mean the system has lost magnetism? .Also</div><div>I have checked the  converged structure, there is the separation at the grain boundary.How can I reduce the forces  at grain boundary and how can I keep the GB stable without separation?</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On 6 October 2014 12:22, Nicola Marzari <span dir="ltr"><<a href="mailto:nicola.marzari@epfl.ch" target="_blank">nicola.marzari@epfl.ch</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
<br>
also, read this:<br>
<a href="http://theossrv1.epfl.ch/Main/ElectronicTemperature" target="_blank">http://theossrv1.epfl.ch/Main/<u></u>ElectronicTemperature</a><br>
<br>
nicola<br>
<br>
-------- Original Message --------<br>
Subject: Re: [Pw_forum] convergence of sigma 5 GB<br>
Date: Mon, 06 Oct 2014 20:12:51 +0200<br>
From: Nicola Marzari <<a href="mailto:nicola.marzari@epfl.ch" target="_blank">nicola.marzari@epfl.ch</a>><br>
Organization: École Polytechnique Fédérale de Lausanne<br>
To: PWSCF Forum <<a href="mailto:pw_forum@pwscf.org" target="_blank">pw_forum@pwscf.org</a>><br>
<br>
<br>
<br>
Provided your cell looks good in xcrysden (i.e. no errors in using an<br>
odd unphysical structure) I would start reducing the beta mixing to 0.2,<br>
or 0.1, or 0.05 - see if any is sufficient to converge, and if so, which<br>
is faster. Try using some of the more sophisticated mixing schemes.<br>
<br>
Before doing that, note that your pseudopotential is ultrasoft, so<br>
you need to have an ecutrho between 6 and 12 times larger than<br>
ecutwfc (ultrasoft pseudopotentials allow for smoother wavefunctions,<br>
but the charge density is what it is - it can't be made smoother.<br>
so the usual relation ecutrho=4*ecutwfc, that is correct for a<br>
rho=|wfc|^2, breaks if the wfc is smoother, and the real rho is<br>
reproduced by |wfv|^2 + augmentation charges, as it happens in ultrasoft<br>
pseudos).<br>
<br>
Last, given that your unit cell is much longer around z, you should<br>
sample with a 4 4 2 mesh, not a 4 4 4. Maybe a 2 2 1 mesh would<br>
be good enough to get things started, and then you increase it at the<br>
end for the last iterations around the converged geometry.<br>
<br>
Note also that typically a 1 1 1 shift (i.e. 2 2 1 1 1 1, or<br>
4 4 2 1 1 1) gives better sampling earlier - but something you could<br>
test in a unit cell of 4 atoms only, with one atom displaced by 5-10%:<br>
check how the forces converges to the true value as you increase the<br>
sampling using n n n 0 0 0 or n n n 1 1 1 meshes.<br>
<br>
                                nicola<br>
<br>
<br>
On 06/10/2014 19:52, Ravi Kiran wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Dear Quantum espresso users,<br>
<br>
I am looking to study the properties of sigma 5 grain boundary with 80<br>
atoms.I am using scf calculations for initial convergence. I am using<br>
the following script. The problem is that the values dont converge even<br>
after 100 iterations .Can anyone please tell me how can I converge the<br>
total energy in order to obtain relaxed atomic coordinates.<br>
&control<br>
     calculation='relax',<br>
     restart_mode='restart',<br>
     prefix='Ni',<br>
     tprnfor =.t.,<br>
     pseudo_dir = '/psuedo/',<br>
     outdir='/psuedo/tmp',<br>
/<br>
  &SYSTEM<br>
      nat =80,<br>
     ntyp = 1,<br>
     degauss = 0.02,<br>
     nspin = 2,<br>
     starting_magnetization = 0.62,<br>
     ecutwfc = 24,<br>
     ibrav = 0,<br>
     occupations = 'smearing',<br>
     smearing = 'mv',<br>
  /<br>
  &ELECTRONS<br>
        conv_thr = 1e-6,<br>
     mixing_beta = 0.3,<br>
     /<br>
&IONS<br>
  /<br>
&CELL<br>
/<br>
<br>
ATOMIC_SPECIES<br>
  Ni 58.69d0 Ni.pz-nd-rrkjus.UPF<br>
K_POINTS automatic<br>
    4 4 4   0 0 0<br>
ATOMIC_POSITIONS<br>
Ni  0.000000  0.000000  0.500000<br>
Ni  0.000000  0.000000  0.250000<br>
Ni  0.400000  0.200000  0.500000<br>
Ni  0.400000  0.200000  0.250000<br>
Ni  0.200000  0.600000  0.500000<br>
Ni  0.200000  0.600000  0.250000<br>
Ni  0.800000  0.400000  0.500000<br>
Ni  0.800000  0.400000  0.250000<br>
Ni  0.600000  0.800000  0.500000<br>
Ni  0.600000  0.800000  0.250000<br>
Ni  0.100000  0.300000  0.500000<br>
Ni  0.100000  0.300000  0.250000<br>
Ni  0.700000  0.100000  0.500000<br>
Ni  0.700000  0.100000  0.250000<br>
Ni  0.500000  0.500000  0.500000<br>
Ni  0.500000  0.500000  0.250000<br>
Ni  0.300000  0.900000  0.500000<br>
Ni  0.300000  0.900000  0.250000<br>
Ni  0.900000  0.700000  0.500000<br>
Ni  0.900000  0.700000  0.250000<br>
Ni  0.500000  0.000000  0.375000<br>
Ni  0.500000  0.000000  0.125000<br>
Ni  0.300000  0.400000  0.375000<br>
Ni  0.300000  0.400000  0.125000<br>
Ni  0.100000  0.800000  0.375000<br>
Ni  0.100000  0.800000  0.125000<br>
Ni  0.900000  0.200000  0.375000<br>
Ni  0.900000  0.200000  0.125000<br>
Ni  0.700000  0.600000  0.375000<br>
Ni  0.700000  0.600000  0.125000<br>
Ni  0.200000  0.100000  0.375000<br>
Ni  0.200000  0.100000  0.125000<br>
Ni  0.000000  0.500000  0.375000<br>
Ni  0.000000  0.500000  0.125000<br>
Ni  0.600000  0.300000  0.375000<br>
Ni  0.600000  0.300000  0.125000<br>
Ni  0.400000  0.700000  0.375000<br>
Ni  0.400000  0.700000  0.125000<br>
Ni  0.800000  0.900000  0.375000<br>
Ni  0.800000  0.900000  0.125000<br>
Ni  0.000000  0.000000  0.750000<br>
Ni  0.000000  0.000000  1.000000<br>
Ni  0.600000  0.200000  0.750000<br>
Ni  0.600000  0.200000  1.000000<br>
Ni  0.200000  0.400000  0.750000<br>
Ni  0.200000  0.400000  1.000000<br>
Ni  0.800000  0.600000  0.750000<br>
Ni  0.800000  0.600000  1.000000<br>
Ni  0.400000  0.800000  0.750000<br>
Ni  0.400000  0.800000  1.000000<br>
Ni  0.300000  0.100000  0.750000<br>
Ni  0.300000  0.100000  1.000000<br>
Ni  0.900000  0.300000  0.750000<br>
Ni  0.900000  0.300000  1.000000<br>
Ni  0.500000  0.500000  0.750000<br>
Ni  0.500000  0.500000  1.000000<br>
Ni  0.100000  0.700000  0.750000<br>
Ni  0.100000  0.700000  1.000000<br>
Ni  0.700000  0.900000  0.750000<br>
Ni  0.700000  0.900000  1.000000<br>
Ni  0.800000  0.100000  0.625000<br>
Ni  0.800000  0.100000  0.875000<br>
Ni  0.400000  0.300000  0.625000<br>
Ni  0.400000  0.300000  0.875000<br>
Ni  0.000000  0.500000  0.625000<br>
Ni  0.000000  0.500000  0.875000<br>
Ni  0.600000  0.700000  0.625000<br>
Ni  0.600000  0.700000  0.875000<br>
Ni  0.200000  0.900000  0.625000<br>
Ni  0.200000  0.900000  0.875000<br>
Ni  0.500000  0.000000  0.625000<br>
Ni  0.500000  0.000000  0.875000<br>
Ni  0.100000  0.200000  0.625000<br>
Ni  0.100000  0.200000  0.875000<br>
Ni  0.700000  0.400000  0.625000<br>
Ni  0.700000  0.400000  0.875000<br>
Ni  0.300000  0.600000  0.625000<br>
Ni  0.300000  0.600000  0.875000<br>
Ni  0.900000  0.800000  0.625000<br>
Ni  0.900000  0.800000  0.875000<br>
CELL_PARAMETERS<br>
14.87  0     0<br>
0    14.87   0<br>
0     0    26.60<br>
<br>
<br>
<br>
--<br>
Thanks<br>
Ravi<br>
<br>
<br>
______________________________<u></u>_________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank">http://pwscf.org/mailman/<u></u>listinfo/pw_forum</a><br>
<br>
</blockquote>
<br>
-- <br>
<br>
------------------------------<u></u>------------------------------<u></u>----------<br>
Prof Nicola Marzari, Chair of Theory and Simulation of Materials, EPFL<span class="HOEnZb"><font color="#888888"><br>
<br>
-- <br>
------------------------------<u></u>------------------------------<u></u>---------<br>
Prof Nicola Marzari   Department of Materials Science and Engineering<br>
13-5066   MIT   77 Massachusetts Avenue   Cambridge MA 02139-4307 USA<br>
tel <a href="tel:617.4522758" value="+16174522758" target="_blank">617.4522758</a> fax 2586534 <a href="mailto:marzari@mit.edu" target="_blank">marzari@mit.edu</a> <a href="http://quasiamore.mit.edu" target="_blank">http://quasiamore.mit.edu</a><br>
<br>
<br>
<br>
-- <br>
<br>
------------------------------<u></u>------------------------------<u></u>----------<br>
Prof Nicola Marzari, Chair of Theory and Simulation of Materials, EPFL<br>
</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <br>Thanks<div>Ravi</div>
</div>