<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META http-equiv=Content-Type content=text/html;charset=iso-8859-1>
<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>>A recent result (Jpn. J. Appl. Phys. 39 (2000) L847) showed that
different pseudopotential leaded to different conclusions. </DIV>
<DIV><FONT face=Arial size=2>thank you for your information.</FONT></DIV>
<DIV><FONT face=Arial size=2>I have no subscription for this <FONT
face="Times New Roman" size=3>japanese journal. Please sent it to <A
title="mailto:vegalew@gmail.com
CTRL + 单击以下链接"
href="mailto:vegalew@gmail.com">vegalew@gmail.com</A> if you like me to read
it.</FONT></FONT></DIV>
<DIV><SPAN class=en><FONT face=Arial><FONT size=2>I am <FONT
color=#000000>rather <SPAN class=heighlight>curious</SPAN> to know <SPAN
class=heighlight>about</SPAN> it.</FONT></FONT></FONT></SPAN></DIV>
<DIV><SPAN class=en><FONT face=Arial size=2></FONT></SPAN> </DIV>
<DIV><SPAN class=en><FONT face=Arial size=2>vega</FONT></SPAN></DIV>
<DIV><FONT face=Arial size=2></FONT> </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>
<DIV style="FONT: 10pt Tahoma">
<DIV><BR></DIV>
<DIV style="BACKGROUND: #f5f5f5">
<DIV style="font-color: black"><B>From:</B> <A
title="mailto:willykohn@gmail.com
CTRL + 单击以下链接"
href="mailto:willykohn@gmail.com">willy kohn</A> </DIV>
<DIV><B>Sent:</B> Saturday, August 30, 2008 2:02 AM</DIV>
<DIV><B>To:</B> <A title="mailto:pw_forum@pwscf.org
CTRL + 单击以下链接"
href="mailto:pw_forum@pwscf.org">PWSCF Forum</A> </DIV>
<DIV><B>Subject:</B> Re: [Pw_forum] something strange in rutile and anatase
calculation</DIV></DIV></DIV>
<DIV><BR></DIV>
<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><<A
href="mailto:vegalew@hotmail.com">vegalew@hotmail.com</A>></SPAN> wrote:<BR>
<BLOCKQUOTE class=gmail_quote
style="PADDING-LEFT: 1ex; MARGIN: 0pt 0pt 0pt 0.8ex; BORDER-LEFT: rgb(204,204,204) 1px solid">
<DIV style="PADDING-RIGHT: 10px; PADDING-LEFT: 10px; PADDING-TOP: 15px"
name="Compose message area" bgcolor="#ffffff">
<DIV><FONT face=Arial size=2>Dear all,</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </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, I found the calculated energy of anatase lower than
rutile. 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 -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> </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> </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>&CONTROL<BR>
title = 'Anatase lattice'
,<BR>
calculation = 'vc-relax'
,<BR>
restart_mode = 'from_scratch'
,<BR>
outdir = '/home/hjfeng/vega/TiO2/Anatase/lattice/vctest/tmp/'
,<BR>
wfcdir = '/tmp/'
,<BR>
pseudo_dir = '/home/hjfeng/vega/espresso-4.0/pseudo/'
,<BR>
prefix = 'Anatase lattice'
,<BR>
disk_io = 'none'
,<BR>
nstep = 1000
,<BR> /<BR> &SYSTEM<BR>
ibrav =
6,<BR>
celldm(1) =
7.1356,<BR>
celldm(3) =
2.5122,<BR>
nat =
12,<BR>
ntyp =
2,<BR>
ecutwfc = 40
,<BR>
ecutrho = 400 ,<BR> /<BR> &ELECTRONS</FONT></DIV>
<DIV> </DIV>
<DIV><FONT face=Arial
size=2> /<BR> &IONS<BR>
ion_dynamics = 'bfgs'
,<BR> /<BR> &CELL<BR>
cell_dynamics = 'bfgs' ,<BR> /<BR>ATOMIC_SPECIES<BR>
Ti 47.86700 Ti.pw91-sp-van_ak.UPF <BR>
O 15.99940 O.pw91-van_ak.UPF <BR>ATOMIC_POSITIONS crystal
<BR> Ti
0.000000000 0.000000000
0.000000000 <BR>
Ti 0.500000000
0.500000000 0.500000000 <BR>
Ti 0.000000000
0.500000000 0.250000000 <BR>
Ti 0.500000000
0.000000000 0.750000000
<BR> O
0.000000000 0.500000000
0.042000000 <BR>
O 0.000000000
0.000000000 0.208000000
<BR> O
0.500000000 0.500000000
0.292000000 <BR>
O 0.000000000
0.500000000 0.458000000
<BR> O
0.500000000 0.000000000
0.542000000 <BR>
O 0.500000000
0.500000000 0.708000000
<BR> O
0.000000000 0.000000000
0.792000000 <BR>
O 0.500000000
1.000000000 0.958000000 <BR>K_POINTS
automatic <BR> 5 5 2 1 1 1 </FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>for rutile case,</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV> &CONTROL<BR>
title = Rutile_lattice
,<BR>
calculation = 'vc-relax'
,<BR>
restart_mode = 'from_scratch'
,<BR>
outdir = '/home/vega32/TiO2/rutile/lattice/vcrelax/tmp'
,<BR>
wfcdir = '/tmp/'
,<BR>
pseudo_dir = '/home/vega32/espresso-4.0.1/pseudo/'
,<BR>
prefix = 'Anatase lattice'
,<BR>
disk_io = 'none'
,<BR>
nstep = 1000
,<BR> /<BR> &SYSTEM<BR>
ibrav =
6,<BR>
celldm(1) =
8.6814,<BR>
celldm(3) =
0.6441,<BR>
nat =
6,<BR>
ntyp =
2,<BR>
ecutwfc = 40
,<BR>
ecutrho = 400 ,<BR> /<BR> &ELECTRONS</DIV>
<DIV> </DIV>
<DIV> /<BR> &IONS<BR>
ion_dynamics = 'bfgs'
,<BR> /<BR> &CELL<BR>
cell_dynamics = 'bfgs' ,<BR> /<BR>ATOMIC_SPECIES<BR>
O 15.99940 O.pw91-van_ak.UPF<BR> Ti
47.90000 Ti.pw91-sp-van_ak.UPF<BR>ATOMIC_POSITIONS
crystal <BR>
O 0.304800000
0.304800000 0.000000000
<BR> O
0.695200000 0.695200000
0.000000000 <BR>
O 0.195200000
0.804800000 0.500000000
<BR> O
0.804800000 0.195200000
0.500000000 <BR>
Ti 0.500000000
0.500000000 0.500000000 <BR>
Ti 0.000000000
0.000000000 0.000000000 <BR>K_POINTS
automatic <BR> 4 4 7 1 1 1 </DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial size=2>vega</FONT></DIV>
<DIV><FONT face=Arial size=2></FONT> </DIV>
<DIV><FONT face=Arial
size=2>==============================================================<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>
<P>
<HR>
<P></P>_______________________________________________<BR>Pw_forum mailing
list<BR>Pw_forum@pwscf.org<BR>http://www.democritos.it/mailman/listinfo/pw_forum<BR></BODY></HTML>