[QE-users] Elastic_constant calculations for supercell using thermo_pw.x

Shiferaw Gadisa shifgadisa at gmail.com
Wed Feb 19 07:33:44 CET 2020


When I calculate the elastic constants for supercells of my system using
thermo_pw.x. I found the compatibility problem of ibrav. My system is
tetragonal
Please how can I solve the issue and continue my calculations?
Thank you in advance for any assistance
The input / out files of calculations are shown as follows

&INPUT_THERMO
  what='mur_lc_elastic_constants',
/
&CONTROL
                 calculation = 'scf' ,
                     pseudo_dir = './',
                         outdir='temp/'
                      prefix = 'psto50' ,
 /
 &SYSTEM
                       ibrav = 6,
                   celldm(1) = 14.74,
                   celldm(3) = 1.06,
                         nat = 20,
                        ntyp = 4,
                     ecutwfc = 40.0 ,
                        nbnd = 62,
 /
 &ELECTRONS
                    conv_thr = 1.0d-8 ,
                 mixing_mode = 'plain' ,
                 mixing_beta = 0.7 ,
             diagonalization = 'davidson' ,
 /
ATOMIC_SPECIES
   Pb  207.20000  Pb.pbe-mt_fhi.UPF
   Sn  118.70000  Sn.pbe-mt_fhi.UPF
   Ti   47.86700  Ti.pbe-mt_fhi.UPF
    O   15.99940  O.pbe-mt_fhi.UPF
ATOMIC_POSITIONS (angstrom)
Sn       0.000000000   0.000000000   0.000000000
Ti       1.995794146   2.008568413   2.009579978
O        1.984411263   2.008567974  -0.000000000
O        0.000000000   2.008568082   2.009579978
O        1.985281957   0.000000000   2.009579978
Sn       0.000000000   4.017139736   0.000000000
Ti       1.995794146   6.025711058   2.009579978
O        1.984411263   6.025711497  -0.000000000
O        0.000000000   6.025711389   2.009579978
O        1.985283246   4.017139736   2.009579978
Pb       4.016267350   0.000000000   0.000000000
Ti       6.036740554   2.008568413   2.009579978
O        6.048123437   2.008567974  -0.000000000
O        4.016267350   2.008566335   2.009579978
O        6.047252743   0.000000000   2.009579978
Pb       4.016267350   4.017139736   0.000000000
Ti       6.036740554   6.025711058   2.009579978
O        6.048123437   6.025711497  -0.000000000
O        4.016267350   6.025713136   2.009579978
O        6.047251454   4.017139736   2.009579978
K_POINTS automatic
 4  4  4  1  1  1

*Output*

 Reading input from _temporary_1
               file Pb.pbe-mt_fhi.UPF: wavefunction(s)  5f renormalized
               file Sn.pbe-mt_fhi.UPF: wavefunction(s)  5d 4f renormalized
               file Ti.pbe-mt_fhi.UPF: wavefunction(s)  4f renormalized
               file O.pbe-mt_fhi.UPF: wavefunction(s)  4f renormalized

   Info: using nr1, nr2, nr3 values from input

   Info: using nr1, nr2, nr3 values from input

     C_1 (1) is incompatible with the tetragonal Bravais lattice
     It is compatible with the
     triclinic Bravais lattice; ibrav=   14
     You might want to change the Bravais lattice or to
     understand why the symmetries are wrong before continuing
     The point group or the Laue class are not used to reduce the number of
     computed tensor components
Info: using nr1, nr2, nr3 values from input
     Computing the frozen ions elastic constants at the minimum volume
     FFT mesh: (   60,   60,   64 )

     Bravais lattice:

     ibrav=  6: tetragonal
     Cell parameters:

     alat=  14.740000 a.u., c/a=   1.060000


     Starting primitive lattice vectors:
     crystal axes: (cart. coord. in units of alat)

               a(1) = (   1.000000   0.000000   0.000000 )
               a(2) = (   0.000000   1.000000   0.000000 )
               a(3) = (   0.000000   0.000000   1.060000 )

     Starting reciprocal lattice vectors:
     reciprocal axes: (cart. coord. in units 2 pi/alat)

               b(1) = (  1.000000  0.000000  0.000000 )
               b(2) = (  0.000000  1.000000  0.000000 )
               b(3) = (  0.000000  0.000000  0.943396 )
     ......


Shiferaw G.
PhD student
Department of physics
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20200219/67fb19af/attachment.html>


More information about the users mailing list