<div dir="ltr">Hi, Vega:<br><br>A recent result (Jpn. J. Appl. Phys. 39 (2000) L847) showed that different pseudopotential leaded to different conclusions. <br><br>Best,<br><br>Wei<br><br><div class="gmail_quote">On Thu, Aug 28, 2008 at 4:28 AM, vega <span dir="ltr">&lt;<a href="mailto:vegalew@hotmail.com">vegalew@hotmail.com</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">



<div style="padding-right: 10px; padding-left: 10px; padding-top: 15px;" bgcolor="#ffffff" name="Compose message area">
<div><font size="2" face="Arial">Dear all,</font></div>
<div><font size="2" face="Arial"></font>&nbsp;</div>
<div><font size="2" face="Arial">It is known that the rutile is more stable than 
anatase in thermodynamics. Surprisingly, when I did the vc-relax calculation by 
QE,&nbsp;I found the&nbsp;calculated energy&nbsp;of&nbsp;anatase lower than 
rutile.&nbsp;Both vc-relax calculation for rutile and anatase onverged in 6 scf 
cycles and 3 bfgs steps, giving the final enthalpy = -362.7585836890 Ry for 
rutile and -725.5447425835 Ry for antase. So the average energy of [TiO2] unit 
is&nbsp;-181.3792918445 Ry for rutile and -181.386185645875 Ry for anatase. It 
looks like anatase is more stable than rutile. Do you think it was quite 
strange?</font></div>
<div><font size="2" face="Arial"></font>&nbsp;</div>
<div><font size="2" face="Arial">To make the calculation both for anatase and rutile 
comparable, I used the same ecutwfc=40 Ry and ecutwfc=400 Ry, and a quite dense 
k-point mesh for both case, 4x4x7 for rutile and 5x5x2 for anatase.</font></div>
<div><font size="2" face="Arial"></font>&nbsp;</div>
<div><font size="2" face="Arial">for better understanding in my calculation, please 
read the input file below,</font></div>
<div><font size="2" face="Arial">for anatase case,</font></div>
<div><font size="2" face="Arial">&amp;CONTROL<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
title = &#39;Anatase lattice&#39; 
,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
calculation = &#39;vc-relax&#39; 
,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
restart_mode = &#39;from_scratch&#39; 
,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
outdir = &#39;/home/hjfeng/vega/TiO2/Anatase/lattice/vctest/tmp/&#39; 
,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
wfcdir = &#39;/tmp/&#39; 
,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
pseudo_dir = &#39;/home/hjfeng/vega/espresso-4.0/pseudo/&#39; 
,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
prefix = &#39;Anatase lattice&#39; 
,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
disk_io = &#39;none&#39; 
,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
nstep = 1000 
,<br>&nbsp;/<br>&nbsp;&amp;SYSTEM<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
ibrav = 
6,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
celldm(1) = 
7.1356,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
celldm(3) = 
2.5122,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
nat = 
12,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
ntyp = 
2,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
ecutwfc = 40 
,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
ecutrho = 400 ,<br>&nbsp;/<br>&nbsp;&amp;ELECTRONS</font></div>
<div>&nbsp;</div>
<div><font size="2" face="Arial">&nbsp;/<br>&nbsp;&amp;IONS<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
ion_dynamics = &#39;bfgs&#39; 
,<br>&nbsp;/<br>&nbsp;&amp;CELL<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
cell_dynamics = &#39;bfgs&#39; ,<br>&nbsp;/<br>ATOMIC_SPECIES<br>&nbsp;&nbsp; 
Ti&nbsp;&nbsp; 47.86700&nbsp; Ti.pw91-sp-van_ak.UPF <br>&nbsp;&nbsp;&nbsp; 
O&nbsp;&nbsp; 15.99940&nbsp; O.pw91-van_ak.UPF <br>ATOMIC_POSITIONS crystal 
<br>&nbsp;&nbsp; Ti&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; 
0.000000000&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp; 
Ti&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.500000000&nbsp;&nbsp;&nbsp; 
0.500000000&nbsp;&nbsp;&nbsp; 0.500000000&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp; 
Ti&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; 
0.500000000&nbsp;&nbsp;&nbsp; 0.250000000&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp; 
Ti&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.500000000&nbsp;&nbsp;&nbsp; 
0.000000000&nbsp;&nbsp;&nbsp; 0.750000000&nbsp;&nbsp;&nbsp; 
<br>&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
0.000000000&nbsp;&nbsp;&nbsp; 0.500000000&nbsp;&nbsp;&nbsp; 
0.042000000&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp;&nbsp; 
O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; 
0.000000000&nbsp;&nbsp;&nbsp; 0.208000000&nbsp;&nbsp;&nbsp; 
<br>&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
0.500000000&nbsp;&nbsp;&nbsp; 0.500000000&nbsp;&nbsp;&nbsp; 
0.292000000&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp;&nbsp; 
O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; 
0.500000000&nbsp;&nbsp;&nbsp; 0.458000000&nbsp;&nbsp;&nbsp; 
<br>&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
0.500000000&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; 
0.542000000&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp;&nbsp; 
O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.500000000&nbsp;&nbsp;&nbsp; 
0.500000000&nbsp;&nbsp;&nbsp; 0.708000000&nbsp;&nbsp;&nbsp; 
<br>&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
0.000000000&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; 
0.792000000&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp;&nbsp; 
O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.500000000&nbsp;&nbsp;&nbsp; 
1.000000000&nbsp;&nbsp;&nbsp; 0.958000000&nbsp;&nbsp;&nbsp; <br>K_POINTS 
automatic <br>&nbsp; 5 5 2&nbsp;&nbsp; 1 1 1 </font></div>
<div><font size="2" face="Arial"></font>&nbsp;</div>
<div><font size="2" face="Arial">for rutile case,</font></div>
<div><font size="2" face="Arial"></font>&nbsp;</div>
<div>&nbsp;&nbsp;&nbsp;&nbsp;&amp;CONTROL<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
title = Rutile_lattice 
,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
calculation = &#39;vc-relax&#39; 
,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
restart_mode = &#39;from_scratch&#39; 
,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
outdir = &#39;/home/vega32/TiO2/rutile/lattice/vcrelax/tmp&#39; 
,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
wfcdir = &#39;/tmp/&#39; 
,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
pseudo_dir = &#39;/home/vega32/espresso-4.0.1/pseudo/&#39; 
,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
prefix = &#39;Anatase lattice&#39; 
,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
disk_io = &#39;none&#39; 
,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
nstep = 1000 
,<br>&nbsp;/<br>&nbsp;&amp;SYSTEM<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
ibrav = 
6,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
celldm(1) = 
8.6814,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
celldm(3) = 
0.6441,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
nat = 
6,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
ntyp = 
2,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
ecutwfc = 40 
,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
ecutrho = 400 ,<br>&nbsp;/<br>&nbsp;&amp;ELECTRONS</div>
<div>&nbsp;</div>
<div>&nbsp;/<br>&nbsp;&amp;IONS<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
ion_dynamics = &#39;bfgs&#39; 
,<br>&nbsp;/<br>&nbsp;&amp;CELL<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
cell_dynamics = &#39;bfgs&#39; ,<br>&nbsp;/<br>ATOMIC_SPECIES<br>&nbsp;&nbsp;&nbsp; 
O&nbsp;&nbsp; 15.99940&nbsp; O.pw91-van_ak.UPF<br>&nbsp;&nbsp; Ti&nbsp;&nbsp; 
47.90000&nbsp; Ti.pw91-sp-van_ak.UPF<br>ATOMIC_POSITIONS 
crystal&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
0.304800000&nbsp;&nbsp;&nbsp; 0.304800000&nbsp;&nbsp;&nbsp; 
0.000000000&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp;&nbsp; 
O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.695200000&nbsp;&nbsp;&nbsp; 
0.695200000&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; 
<br>&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
0.195200000&nbsp;&nbsp;&nbsp; 0.804800000&nbsp;&nbsp;&nbsp; 
0.500000000&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp;&nbsp; 
O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.804800000&nbsp;&nbsp;&nbsp; 
0.195200000&nbsp;&nbsp;&nbsp; 0.500000000&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp; 
Ti&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.500000000&nbsp;&nbsp;&nbsp; 
0.500000000&nbsp;&nbsp;&nbsp; 0.500000000&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp; 
Ti&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; 
0.000000000&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; <br>K_POINTS 
automatic <br>&nbsp; 4 4 7&nbsp;&nbsp; 1 1 1 </div>
<div><font size="2" face="Arial"></font>&nbsp;</div>
<div><font size="2" face="Arial">vega</font></div>
<div><font size="2" face="Arial"></font>&nbsp;</div>
<div><font size="2" face="Arial">==============================================================<br><font color="#888888">Vega 
Lew ( weijia liu)<br>PH.D Candidate in Chemical Engineering<br>State Key 
Laboratory of Materials-oriented Chemical Engineering<br>College of Chemistry 
and Chemical Engineering<br>Nanjing University of Technology, 210009, Nanjing, 
Jiangsu, China</font></font></div></div>
<br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">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><br></div>