[Pw_forum] Vectors in rhombohedral cells

Marino Vetuschi Zuccolini zucco at dipteris.unige.it
Mon Nov 10 14:11:03 CET 2008


Dear all,
working with a rhombohedral cell I've found something that appear to  
me to going wrong, and I call you for a help. Probably is only my  
mistake but I ask you a confirmation.

The problem arises when I compare the results of a vc-relax output as  
below:

...
                         ibrav = 5,
                     celldm(1) = 10.9072,
                     celldm(4) = 0.67279,
....
      End of BFGS Geometry Optimization

      Final enthalpy =    -325.7624962652 Ry
      new unit-cell volume =    666.19980 a.u.^3 (    98.72062 Ang^3 )

CELL_PARAMETERS (alat)
    0.411520633  -0.237591548   0.875166762
    0.000000000   0.475183096   0.875166762
   -0.411520633  -0.237591548   0.875166762
...

and the few summary lines of an scf output of the optimized structure  
with c=0.661301537
...
      bravais-lattice index     =            5
      lattice parameter (a_0)   =      10.8619  a.u.
      unit-cell volume          =     661.4872 (a.u.)^3
...

      crystal axes: (cart. coord. in units of a_0)
                a(1) = (  0.411521 -0.237592  0.879887 )
                a(2) = (  0.000000  0.475183  0.879887 )
                a(3) = ( -0.411521 -0.237592  0.879887 )

You can see the difference in what concerns the volume.

To compute a0 and alpha of the optimized cell I used what is reported  
in the PW manual

       v1 = a(tx,-ty,tz),   v2 = a(0,2ty,tz),   v3 = a(-tx,-ty,tz).
where c=cos(alpha) is the cosine of the angle alpha between any pair
of crystallographic vectors, tc, ty, tz are defined as
      tx=sqrt((1-c)/2), ty=sqrt((1-c)/6), tz=sqrt((1+2c)/3)

but...if I calculate c (as cos(alpha)) from tx, ty and tz I noticed  
that tx and ty are good but tz is not.

 From the vc-relax output I obtain:
c from tx ==> 0.661301537
c from ty ==> 0.661301537
c from tz ==> 0.648875292 <=== ?

 From scf output I obtain, considering the truncation of the matrix'  
elements
c from tx ==> 0.661300933
c from ty ==> 0.661300249
c from tz ==> 0.661301699

Thanks for any help

m.

*******************************************************
Marino Vetuschi Zuccolini
zucco at dipteris.unige.it
Researcher / Geochemist
Laboratory of Geochemistry

DIPartimento per lo studio della TErra e delle sue RISorse -  
Università  di Genova
Tel. ++39 010 3538136 Fax. ++39 010 352169
Corso Europa 26, 16132 - Genova - Italy




More information about the Pw_forum mailing list