[Pw_forum] bug in data input?

Carlo Nervi carlo.nervi at unito.it
Wed Feb 24 18:52:43 CET 2010


Dear all,
I did some geometry optimizations with pw.x starting from a monoclininc 
cell, using Angstrom and the "A, B, C, cosAB, cosAC, cosBC" input 
methods, for example:

  &SYSTEM
                        ibrav = 12,
                            A = 8.0828 ,
                            B = 12.5828 ,
                            C = 9.7644 ,
                        cosAB = 0.0d0 ,
                        cosAC = -0.10713175431387366149628407325294d0 ,
                        cosBC = 0.0d0 ,
                          nat = 104,
                         ntyp = 4,
                      ecutwfc = 120 ,
  /

The input cif file contains the usual a, b, c, and alfa, beta=96.150, 
gamma. Therefore, I think, it should be correct to set cosAC as shown 
above. However, after open the input file with XCrysden, I realized that 
the cell was 90-90-90 (a parallelepiped).
Where I am wrong?

Playing a little bit with cosAB, cosAC and cosBC I understood that only 
cosAB makes the correct angle (96 degrees) to appears in XCrysden.
Therefore I tried another approach: switch b and c, and to use the 
celldm() approach. Of course I had to use crystallographic fractional 
coordinates and to exchange y with z in the Atomic coordinates.
I'm not sure this is the best approach, but apparently it works.
I would appreciate any suggestions from the community.

  &SYSTEM
                        ibrav = 12,
                    celldm(1) = 15.274655166 ,
                    celldm(2) = 1.20797 ,
                    celldm(3) = 1.5567 ,
                    celldm(4) = -0.107131784 ,
                          nat = 104,
                         ntyp = 4,
                      ecutwfc = 120 ,
  /

Is it normal that pw.x behaves in this way? In other words, are there 
any restrictions using the "A, B, C, cosAB, cosAC, cosBC" input method?

Thank you in advance,
	Carlo

-- 
------------------------------------------------------
Carlo Nervi carlo.nervi at unito.it Tel:+39 011 6707507/8
Fax: +39 011 6707855   -   Dipartimento di Chimica IFM
via P. Giuria 7, 10125 Torino, Italy
http://lem.ch.unito.it/



More information about the users mailing list