[Wien] odd failure in 'qtl' - triclinic

Cormac McGuinness cmcguin at tcd.ie
Fri Aug 29 00:44:29 CEST 2008


Dear mailinglist

I've encountered an odd error in the use of the qtl program.
When studying a supercell of a wurtzite compound such as ZnO or AlN where one
metal atom is substituted by another metal dopant, I have ended up with a
triclinic supercell where the wurtzite c-axis is the new a -axis and alpha is
120 degrees. As the local rotation matrices of all the atoms are the identity
matrix, I then attempted to use 'qtl' to rotate the z-axis to be parallel to
the wurtzite c-axis, e.g. :
=========== snippet from case.inq ======
  16  -2  1  1       iatom,qsplit,symmetrize,locrot
3   0  1  2          nL, l-values
1.0 0.0 0.0
============

However the qtl fails with the following error:
==========
Input/Output Error 156: Read after Write

   In Procedure: main program
        At Line: 497

      Statement: Formatted READ
           Unit: 33
   Connected To: fort.33
           Form: Formatted
         Access: Sequential
Records Read   : 0
Records Written: 14170

End of diagnostics

303.190u 11.348s 16:03.30 32.6%	0+0k 0+0io 0pf+0w
error: command   /usr/local/Wien2k/qtl upqtl.def   failed
===========

while in the output of  case.outputqup contains:
========== snippet of case.outputq
atom 16; type  16; qsplit=  2; for L=  0  1  2
 Symmetrization over eq. k-points is performed
 New z axis ||    1.0000   0.0000   0.0000
 LATTICE:P
  Triclinic implemented, but never tested
  New local rotation matrix in global orthogonal system
                       new x     new y     new z
LOCAL ROT MATRIX:    0.0000000 0.0000000 1.0000000
                     1.0000000 0.0000000 0.0000000
                     0.0000000 1.0000000 0.0000000
===============

Has anyone else tested the application of the new 'qtl'  program to a triclinic
case ?
Is there anyway for me to diagnose the error further ?

Yours
Cormac
--
Dr. Cormac McGuinness, Lecturer, School of Physics, Trinity College Dublin.
Phone: +353 1 8963547    Fax: +353 1 6711759    Office: SNIAM 3.17



More information about the Wien mailing list