The ABINIT Users Mailing List ( CLOSED )

[Xavier Gonze <xavier.gonze@uclouvain.be>]

Dear Petr,

This error message in ABINIT is to be improved.

Anyhow, this problem is due to the fact that 1.15470052489312
is not accurate enough. Should be 1.15470053838 .

Regards,
Xavier

On 07 Jun 2009, at 22:09, Petr Sestak wrote:

> Hi Abinit users,
>
> I am doing some shear deformation in Molybdem. I created a standard  
> BCC basis
> with these RPRIM parameters:
> RPRIM:
> -0.5 0.5 0.5
> 0.5 -0.5 0.5
> 0.5 0.5 -0.5
> and everything is all right. However, in my case the following set  
> of RPRIM is
> more suitable:
> 0.408248290463863 0.7071067811865475244 1.15470052489312
> 0 1.414213562373095049 0
> 0 0 0.86602540378444
> These parameters are creating the same structure as I mentioned  
> above but when
> I run the Abinit with these parameters the program stops and print  
> following
> error in the log file:
> ==========================================
> symbrav : BUG -
> One of the coordinates of rprimd in axes is non-integer,
> or non-half-integer (if centering).
> coord= -5.000000E-01  1.000000E+00 -5.000000E-01
>          0.000000E+00  0.000000E+00 -1.000000E+00
>         -5.000000E-01  5.000000E-01  0.000000E+00
> Action : contact ABINIT group.
>
> .Delivered    0 WARNINGs and   0 COMMENTs to log file.
> -P-0000
> -P-0000  leave_new : decision taken to exit ...
> ==========================================
> Does anybody have an idea what is wrong? In attachment I am sending  
> my input
> (see below).
>
> Best regards
>
> --
> Petr Sestak
> Institute of Engineering Physics
> Faculty of Mechanical Engineering
> Brno University of Technology
> Technicka 2, Brno 616 69, Czech Republic
>
> Mail:          sestak@fme.vutbr.cz
>
> .Version 5.6.4  of ABINIT
>
> -----------------------
> Input file:
> acell 3*5.9901113615
> rprim
> 0.408248290463863 0.7071067811865475244 1.15470052489312
> 0 1.414213562373095049 0
> 0 0 0.86602540378444
> xred
> 0 0 0
> ntypat 1
> znucl 42
> natom 1
> typat 1
> #chkprim 0
> #nsym 1
>
> tsmear 0.01
> #iscf 2
> diemix 0.3
> ecut 700 eV
> #pawecutdg 520 eV
> kptopt 1
> ngkpt 15 15 13
> nstep 30
> toldfe 1.0d-6 eV
> #tolmxf 5.0d-9
> #tolvrs 1.0d-19
> nband 7
> #getwfk -1
> #irdwfk 1
> occopt 4
> ecutsm 0.5
> diemac 1.0d+6
>
> -------------------------
> *.files:
> abinit_config.in
> abinit_config.out
> abinit_configi
> abinit_configo
> abinit_config
> Mo_GGA_abinit