<!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>&gt;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&#10;CTRL + &#21333;&#20987;&#20197;&#19979;&#38142;&#25509;" 
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>&nbsp;</DIV>
<DIV><SPAN class=en><FONT face=Arial size=2>vega</FONT></SPAN></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>
<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&#10;CTRL + &#21333;&#20987;&#20197;&#19979;&#38142;&#25509;" 
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&#10;CTRL + &#21333;&#20987;&#20197;&#19979;&#38142;&#25509;" 
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>&lt;<A 
href="mailto:vegalew@hotmail.com">vegalew@hotmail.com</A>&gt;</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>&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><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>