Dear Sidra

If you want to optimize the lattice constant of cubic SrtiO3 it would be better to do it for the simple 5 atom cubic cell, using a 2X2X2 supercell and constraining the system to remain cubic (cell_dofree = 'ibrav') may cause troubles.

best wishes Pietro



Il 6/18/2021 8:58 AM, Sidrah Younus Khan ha scritto:
Dear QE users,

I have 3 questions.As i have just started the calculations so i have a lot of questions 1)I am using v-c relax calculations for my system of strontium titanate. I am using 16 processors on my personal computer PC, but the system does not converge even after 12 hours. I am unable to detect the problem in the input & output file as some errors are coming. I tried to change potentials and also the mixing -beta from 0.6-0.7. 2) Secondly, I need a pseudopotential -Ti.pbesol-spn-rrkjus_psl.1.0.0.upf , if somebody could lead me to that. as i am unable to find it 3) Is there any way to check the band and spin orbit coupling change with temperature , as far as I know the band structure calculations are working on 0K, i need to check the temperature variation effect on band structure and spin orbit coupling.

I am attaching my input file and some errors of output.
&CONTROL
calculation         = 'vc-relax',
restart_mode        = 'from_scratch',
wf_collect          = .true.,
nstep               = 500
pseudo_dir          = " . "
prefix              = 'srTio3relax'
outdir              = './tmp'
etot_conv_thr       = 1.0d-6
forc_conv_thr       = 1.0d-3
verbosity           = 'high'
/

&SYSTEM
ibrav               = 1,
celldm(1)           = 14.78052,
nat                 = 40,
ntyp                = 3,
ecutwfc             = 40.D0,
ecutrho             = 320.D0,
occupations         = 'smearing',
smearing            = 'gauss',
degauss             = 0.01
/

&ELECTRONS
    diagonalization     = 'david',
    conv_thr         =  1.0d-10
    mixing_beta      =  0.6
/
&ions
/
&cell
    cell_dofree='ibrav'
/
K_POINTS {automatic}
4  4  4  0 0 0

ATOMIC_SPECIES
Sr     87.62000  Sr.pbe-spn-rrkjus_psl.1.0.0.UPF
Ti     47.86700  Tl.pbe-dn-rrkjus_psl.1.0.0.UPF
O      15.99940  O.pbe-n-rrkjus_psl.1.0.0.UPF

ATOMIC_POSITIONS {angstrom}
Sr      0.000000   0.000000   0.000000
Sr      0.000000   0.000000   3.945130
Sr      0.000000   3.945130   0.000000
Sr      0.000000   3.945130   3.945130
Sr      3.945130   0.000000   0.000000
Sr      3.945130   0.000000   3.945130
Sr      3.945130   3.945130   0.000000
Sr      3.945130   3.945130   3.945130
Ti      1.972565   1.972565   1.972565
Ti      1.972565   1.972565   5.917695
Ti      1.972565   5.917695   1.972565
Ti      1.972565   5.917695   5.917695
Ti      5.917695   1.972565   1.972565
Ti      5.917695   1.972565   5.917695
Ti      5.917695   5.917695   1.972565
Ti      5.917695   5.917695   5.917695
O       1.972565   0.000000   1.972565
O       1.972565   0.000000   5.917695
O       1.972565   3.945130   1.972565
O       1.972565   3.945130   5.917695
O       5.917695   0.000000   1.972565
O       5.917695   0.000000   5.917695
O       5.917695   3.945130   1.972565
O       5.917695   3.945130   5.917695
O       1.972565   1.972565   0.000000
O       1.972565   1.972565   3.945130
O       1.972565   5.917695   0.000000
O       1.972565   5.917695   3.945130
O       5.917695   1.972565   0.000000
O       5.917695   1.972565   3.945130
O       5.917695   5.917695   0.000000
O       5.917695   5.917695   3.945130
O       0.000000   1.972565   1.972565
O       0.000000   1.972565   5.917695
O       0.000000   5.917695   1.972565
O       0.000000   5.917695   5.917695
O       3.945130   1.972565   1.972565
O       3.945130   1.972565   5.917695
O       3.945130   5.917695   1.972565
O       3.945130   5.917695   5.917695


_Output errors:_
_
_
_1)_CASE: enthalpy_new < enthalpy_old
2) Check: negative core charge=   -0.000006
_3)_WARNING: bfgs curvature condition failed, Theta= 0.977
4) Total force =     0.364085     Total SCF correction =   0.000016
5)
 Parallelization info
     --------------------
     sticks:   dense  smooth     PW     G-vecs:    dense smooth      PW
     Min         395     197     54                23717 8376    1222
     Max         397     199     55                23719 8381    1227
     Sum        6333    3157    877               379485 134059   19597

http://www.quantum-espresso.org/quote <http://www.quantum-espresso.org/quote>

     Parallel version (MPI), running on    16 processors

     MPI processes distributed on     1 nodes
     R & G space division:  proc/nbgrp/npool/nimage =  16
     Reading input from sr.vc_relax.in <http://sr.vc_relax.in>

     Current dimensions of program PWSCF are:
     Max number of different atomic species (ntypx) = 10
     Max number of k-points (npk) =  40000
     Max angular momentum in pseudopotentials (lmaxx) =  3
               file Tl.pbe-dn-rrkjus_psl.1.0.0.UPF: wavefunction(s)  6S 5D renormalized

     Subspace diagonalization in iterative solution of the eigenvalue problem:
     one sub-group per band group will be used
     ELPA distributed-memory algorithm (size of sub-group:  2*  2 procs)

_
_
_
_
_Thanks_
_Best Regards_
_Sidra Younus_
_American university of Sharjah-UAE_
_
_
_
_
_
_



_______________________________________________
Quantum ESPRESSO is supported by MaX (www.max-centre.eu)
users mailing list [email protected]
https://lists.quantum-espresso.org/mailman/listinfo/users
_______________________________________________
Quantum ESPRESSO is supported by MaX (www.max-centre.eu)
users mailing list [email protected]
https://lists.quantum-espresso.org/mailman/listinfo/users

Reply via email to