<html>
<head>
<style>
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
FONT-SIZE: 10pt;
FONT-FAMILY:Tahoma
}
</style>
</head>
<body class='hmmessage'>
Dear all,<br><br>Paulatto gave me the new file for BFGS calculation. I have recompiled my QE and finished the same calculation<br>for Anatase cell parameters. The results is my calculation can be finished without stop at errors. But the parameters<br>is still smaller than literature both majority DFT calculation and experimental data. The parameters are, a=3.7559 <br>b=3.75561 c=9.3978, respectively.<br>My input file as follows, <br><br>&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/vega/espresso-4.0/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/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 default' ,<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; etot_conv_thr = 0.000000735 ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; forc_conv_thr = 0.0011668141375 ,<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.1469,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; celldm(3) = 2.5124,<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 = 25 ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ecutrho = 200 ,<br>&nbsp;/<br>&nbsp;&amp;ELECTRONS<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; conv_thr = 7.3D-8 ,<br>&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 angstrom <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; 1.888000000&nbsp;&nbsp;&nbsp; 1.888000000&nbsp;&nbsp;&nbsp; 4.743000000&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp; Ti&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; 1.888000000&nbsp;&nbsp;&nbsp; 2.372000000&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp; Ti&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.888000000&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; 7.115000000&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; 1.973000000&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.888000000&nbsp;&nbsp;&nbsp; 1.888000000&nbsp;&nbsp;&nbsp; 6.716000000&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; 1.888000000&nbsp;&nbsp;&nbsp; 4.345000000&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.888000000&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; 9.088000000&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.888000000&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; 5.141000000&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; 1.888000000&nbsp;&nbsp;&nbsp; 0.398000000&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.888000000&nbsp;&nbsp;&nbsp; 1.888000000&nbsp;&nbsp;&nbsp; 2.770000000&nbsp;&nbsp;&nbsp; <br>&nbsp;&nbsp;&nbsp; O&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; 0.000000000&nbsp;&nbsp;&nbsp; 7.513000000&nbsp;&nbsp;&nbsp; <br>K_POINTS automatic <br>&nbsp; 4 4 2&nbsp;&nbsp; 1 1 1 <br><br>I think there must be something wrong with my input file. Could you please do me a <br>favor tell me what's wrong and why the cell parameters so small I'm thinking about it <br>day and night.<br><br>By the way, I note that CASTEP has find 32 symmetry operations of my cell. My cell is <br>a standard anatase crystall. But Q-E only can find 4 at most since I adjust parameters of <br>~/espresso-4.0/PW/eqvect.f90. Do this matter effect my results significantly?<br><br>I'm look forwar for your advice.<br><br><font size="4"><span style="font-family: Lucida Handwriting,Cursive;">Vega Lew</span></font><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<br><br><a href="http://imagine-msn.com/messenger/launch80/default.aspx?locale=en-us&amp;source=wlmailtagline" target="_blank"></a>
<br /><hr />Connect to the next generation of MSN Messenger   <a href='http://imagine-msn.com/messenger/launch80/default.aspx?locale=en-us&source=wlmailtagline' target='_new'>Get it now! </a></body>
</html>