
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">

<HTML><HEAD>

<META http-equiv=Content-Type content=text/html;charset=gb2312>

<META content="MSHTML 6.00.6000.16705" name=GENERATOR></HEAD>

<BODY id=MailContainerBody 

style="PADDING-RIGHT: 10px; PADDING-LEFT: 10px; PADDING-TOP: 15px" 

bgColor=#ffffff leftMargin=0 topMargin=0 CanvasTabStop="true" 

name="Compose message area">

<DIV><FONT face=Arial size=2>Dear all,</FONT></DIV>

<DIV><FONT face=Arial size=2></FONT>&nbsp;</DIV>

<DIV><FONT face=Arial size=2>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 face=Arial size=2></FONT>&nbsp;</DIV>

<DIV><FONT face=Arial size=2>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 face=Arial size=2></FONT>&nbsp;</DIV>

<DIV><FONT face=Arial size=2>for better understanding in my calculation, please 

read the input file below,</FONT></DIV>

<DIV><FONT face=Arial size=2>for anatase case,</FONT></DIV>

<DIV><FONT face=Arial 

size=2>&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 = 'Anatase lattice' 

,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

calculation = 'vc-relax' 

,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

restart_mode = 'from_scratch' 

,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

outdir = '/home/hjfeng/vega/TiO2/Anatase/lattice/vctest/tmp/' 

,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

wfcdir = '/tmp/' 

,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

pseudo_dir = '/home/hjfeng/vega/espresso-4.0/pseudo/' 

,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

prefix = 'Anatase lattice' 

,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

disk_io = 'none' 

,<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 face=Arial 

size=2>&nbsp;/<BR>&nbsp;&amp;IONS<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

ion_dynamics = 'bfgs' 

,<BR>&nbsp;/<BR>&nbsp;&amp;CELL<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

cell_dynamics = 'bfgs' ,<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 face=Arial size=2></FONT>&nbsp;</DIV>

<DIV><FONT face=Arial size=2>for rutile case,</FONT></DIV>

<DIV><FONT face=Arial size=2></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 = 'vc-relax' 

,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

restart_mode = 'from_scratch' 

,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

outdir = '/home/vega32/TiO2/rutile/lattice/vcrelax/tmp' 

,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

wfcdir = '/tmp/' 

,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

pseudo_dir = '/home/vega32/espresso-4.0.1/pseudo/' 

,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

prefix = 'Anatase lattice' 

,<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

disk_io = 'none' 

,<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 = 'bfgs' 

,<BR>&nbsp;/<BR>&nbsp;&amp;CELL<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 

cell_dynamics = 'bfgs' ,<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 face=Arial size=2></FONT>&nbsp;</DIV>

<DIV><FONT face=Arial size=2>vega</FONT></DIV>

<DIV><FONT face=Arial size=2></FONT>&nbsp;</DIV>

<DIV><FONT face=Arial 

size=2>==============================================================<BR>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></DIV></BODY></HTML>