<div dir="ltr">Hi <span style="font-size:12.8px">S. Arash Sheikholeslam,</span><br style="font-size:12.8px"><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px">I'm not sure this is the problem but from what you described my understanding is that when you scan that angle manually you don't optimize the other degrees of freedom. On the other hand, when you relax the molecule, all degrees of freedom are relaxed. This discrepancy between the two approaches could explain the behavior you are seeing. If you are interested in doing the manual optimization, you could try to relax the molecule with a constraint on the angle with card CONSTRAINTS.</span></div><div><br></div><div><span style="font-size:12.8px">HTH</span></div><div><span style="font-size:12.8px">Marton Voros</span></div><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px">--</span></div><div><span style="font-size:12.8px">Materials Science Division</span></div><div><span style="font-size:12.8px">Argonne National Laboratory</span></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Apr 17, 2017 at 8:55 PM, sarashs <span dir="ltr"><<a href="mailto:sarashs@ece.ubc.ca" target="_blank">sarashs@ece.ubc.ca</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi,<br>
<br>
I've posted this question before but wasn't successful in getting an<br>
answer. I have a molecule ((OH)3ZrOSi(OH)3) and I want to find the<br>
equilibrium angle of ZrOSi. The equilibrium angle that I get with a<br>
"relax" calculation differs a lot from the result I get by increasing<br>
the angle and performing an "SCF" calculation. I tried various things<br>
including increasing the ecutwfc and ecutrho and fixing the number of<br>
bands for both calculations but had no success. Is there any chance<br>
someone can tell me what I'm missing?<br>
<br>
Here's a sample code my code:<br>
<br>
&CONTROL<br>
    calculation = "scf",<br>
    outdir       = "outputs",<br>
    prefix     ="ZrSi3",<br>
    restart_mode="from_scratch"<br>
    pseudo_dir = "pseudos",<br>
/<br>
&SYSTEM<br>
    nosym       = .TRUE.,<br>
    ibrav       = 0,<br>
    nat         = 15,<br>
    ntyp        = 4,<br>
    ecutwfc     = 600,<br>
    ecutrho = 6000,<br>
    vdw_corr    = "grimme-d2",<br>
/<br>
&ELECTRONS<br>
    conv_thr    = 1.D-11,<br>
    mixing_beta = 0.5D0,<br>
/<br>
ATOMIC_SPECIES<br>
O   15.9994 O.pz-mt.UPF<br>
H   1.00794 H.pz-vbc.UPF<br>
Si  28.0855 Si.pz-vbc.UPF<br>
Zr  91.224 Zr.pz-hgh.UPF<br>
ATOMIC_POSITIONS { angstrom }<br>
H       -0.217652040   2.144165004  -0.967842103<br>
O       -0.644926428   1.299751587  -1.173517450<br>
H       -0.763682000  -2.102234442  -0.637701677<br>
O       -0.631971073  -1.266428593  -1.106959840<br>
H       -2.244617658   0.059203828   1.041286276<br>
O       -1.287708239   0.028876445   1.183017241<br>
Si      -0.341076411   0.043743634  -0.163179170<br>
O        1.183131599   0.097692502   0.401792709<br>
Zr       2.959296806   0.743741575   0.268340451<br>
H        3.161927401   3.483669731   0.077078720<br>
O        2.777805281   2.601338899   0.023924414<br>
H        4.228044154   0.110551898   2.670358608<br>
O        4.003457303   0.315517992   1.757071324<br>
H        3.597999032  -0.225640653  -2.168546302<br>
O        3.766434111  -0.029092443  -1.238766538<br>
K_POINTS { gamma }<br>
CELL_PARAMETERS { angstrom }<br>
13 0 0<br>
0 13 0<br>
0 0 13<br>
<br>
Regards,<br>
<br>
S. Arash Sheikholeslam<br>
PhD candidate at electrical engineering department, University of<br>
British columbia<br>
______________________________<wbr>_________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/<wbr>listinfo/pw_forum</a><br>
</blockquote></div><br></div>