Dear Giuseppe,<div><br></div><div>I tried ecutrho as 8 or 12 times of ecutwfc=30 for ultrasoft PP. Also tried norm conserving PP with ecutwfc=70 and ecutrho=280. The c/a ratio is still far off from the expected 1.4142. Very confused.</div>

<div><span><br></span></div><div>Regards, Wei</div><div><div>Wei Lai<br>Michigan State University<br>East Lansing, MI 48824</div><br><div class="gmail_quote">On Wed, Mar 28, 2012 at 10:43 AM, Giuseppe Mattioli <span dir="ltr"><<a href="mailto:giuseppe.mattioli@ism.cnr.it" target="_blank">giuseppe.mattioli@ism.cnr.it</a>></span> wrote:<br>


<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
Dear Wei<br>
<div><br>
<br>
>  &system<br>
>     ibrav=  7, celldm(1) = 5.6298, celldm(3) = 1.4, nat=  2, ntyp= 2,<br>
>     ecutwfc =30, input_dft='pw91'<br>
>  /<br>
<br>
</div><div>>  ATOMIC_SPECIES<br>
>  Mg  24.305 Mg.pw91-np-van.UPF<br>
>  O   15.999 O.pw91-van_ak.UPF<br>
<br>
</div>You are using ultrasoft pseudopotentials. You should explicitly assign a value to the ecutrho<br>
parameter (&system list). The default value is 4*ecutwfc, and it is a good value for norm conserving<br>
PPs. A 6*ecutwfc~10*ecutwfc value is considered as a safe~good one in the case of ultrasoft PPs<br>
<br>
HTH<br>
<br>
Giuseppe<br>
<br>
--<br>
********************************************************<br>
- Article premier - Les hommes naissent et demeurent<br>
libres et ègaux en droits. Les distinctions sociales<br>
ne peuvent être fondèes que sur l'utilitè commune<br>
- Article 2 - Le but de toute association politique<br>
est la conservation des droits naturels et<br>
imprescriptibles de l'homme. Ces droits sont la libertè,<br>
la propriètè, la sùretè et la rèsistance à l'oppression.<br>
********************************************************<br>
<br>
   Giuseppe Mattioli<br>
   CNR - ISTITUTO DI STRUTTURA DELLA MATERIA<br>
   v. Salaria Km 29,300 - C.P. 10<br>
   I 00015 - Monterotondo Stazione (RM)<br>
   Tel <a href="tel:%2B%2039%2006%2090672836" value="+390690672836" target="_blank">+ 39 06 90672836</a> - Fax <a href="tel:%2B39%2006%2090672316" value="+390690672316" target="_blank">+39 06 90672316</a><br>
   E-mail: <<a href="mailto:giuseppe.mattioli@ism.cnr.it" target="_blank">giuseppe.mattioli@ism.cnr.it</a>><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>
</blockquote></div><br></div>