Dear all I wanted to run some vc-relax calculations. To start with I took pbtio3 system and tried to relax it. for the following input file it predicted the c/a ratio as 1.23. The output structure (both coordinates and lattice parameters) is grossly in error. In the forum i have seen people reporting c/a ratio which is in good agreement with the experimental value. Later I realized for different input configurations it is predicting different c/a ratios. At one time for pbtio3 I got cubic structure after relaxing!!!. As i understood it is stuck up in local minimum. Since i know the structure is tetragonal i could figure out the error but for an unknown system how can i relax the system with confidence. <br>
(if there is some literature that i can read up on this please refere me the same i will read)<br clear="all"> <br>entering subroutine stress ...<br><br>          total   stress  (Ry/bohr**3)                   (kbar)     P=    0.38<br>
   0.00000259   0.00000000   0.00000000          0.38      0.00      0.00<br>   0.00000000   0.00000259   0.00000000          0.00      0.38      0.00<br>   0.00000000   0.00000000   0.00000248          0.00      0.00      0.36<br>
<br><br>     bfgs converged in  34 scf cycles and  27 bfgs steps<br><br>     End of BFGS Geometry Optimization<br><br>     Final enthalpy =    -333.0490579145 Ry<br><br>CELL_PARAMETERS (alat)<br>   0.983673643   0.000000000   0.000000000<br>
   0.000000000   0.983673643   0.000000000<br>   0.000000000   0.000000000   1.208160764<br><br>ATOMIC_POSITIONS (crystal)<br>Ti       0.500000000   0.500000000   0.562608151<br>O        0.500000000   0.500000000  -0.068487782<br>
O        0.500000000   0.000000000   0.453495153<br>O        0.000000000   0.500000000   0.453495153<br>Pb       0.000000000   0.000000000   0.118889324<br><br><br>&control<br>calculation = 'vc-relax',<br>restart_mode='from_scratch',<br>
outdir='/home/sai/pbtio3/'<br>prefix='pbtio3'<br>pseudo_dir='/home/sai/softwares/espresso-4.0.1/pseudo'<br>tstress=.true.,<br>/<br>&system<br>ibrav=0<br>nat=5,<br>ntyp=3,<br>ecutwfc=55,<br>ecutrho=550,<br>
nbnd=22,<br>occupations='fixed',<br>degauss=0.00<br>/<br>&electrons<br>conv_thr=1.D-6,<br>mixing_beta= 0.3<br>/<br>&ions<br>/<br>&cell<br>cell_dynamics='bfgs'<br>/<br>ATOMIC_SPECIES<br>Pb  207.21  Pb.pbe-d-van.UPF<br>
Ti 47.87  Ti.pbe-sp-van_ak.UPF<br>O  16  O.pbe-rrkjus.UPF<br>ATOMIC_POSITIONS {crystal}<br>Ti 0.5 0.5 0.52<br>O 0.5 0.5 0<br>O 0.5 0 0.5<br>O 0 0.5 0.5<br>Pb 0 0 0<br>K_POINTS {automatic}<br>8 8 7 0 0 0<br>CELL_PARAMETERS<br>
7.374 0 0<br>0 7.374 0<br>0 0 7.669<br><br>-- <br>U.Saibabu<br>PhD student,<br>Deformation mechanisms modeling group,<br>Materials engineering department,<br>IISc Bangalore,<br>India.<br><br>