[QE-users] wrong celldm(3) with ibrav -13 in phonon
Pietro Delugas
pdelugas at sissa.it
Tue Aug 4 11:40:23 CEST 2020
Hi it should be fixed in the develop version and in the new release
which should eventually be ready.
in the meantime you can try to download the current release candidate
from this link.
https://gitlab.com/QEF/q-e/-/archive/qe-6.6rc/q-e-qe-6.6rc.tar.gz
Pietro
On 04/08/20 07:12, Zeeshan Ahmad wrote:
> Hi,
>
> I ran a ph.x and pw.x calculation on a ibrav=-13 system. When I run
> the q2r.x in the directory, I got the following error:
>
> Error in routine latgen (13):
> wrong celldm(3)
>
> After a little printing stuff inside the file Modules/latgen.f90, I
> found that the control is going to ibrav=13 (with celldm(3)<0.d0) and
> not ibrav=-13 as I specified in the scf.in file. I wonder if this has
> to do with the recent change of lattice vectors for ibrav=-13 in the
> 6.5 version since the pw.x and ph.x do not give any error.
>
> Here are the input scf, ph and q2r input files for reproducing.
>
> scf.in
>
> &CONTROL
> title = 'Li2CO3'
> calculation = 'scf'
> restart_mode = 'from_scratch'
> outdir = './scf1'
> pseudo_dir = '/home/azeeshan/pseudopot/all_pbe_UPF_v1.5'
> prefix = 'Li2CO3'
> disk_io = 'low'
> verbosity = 'default'
> etot_conv_thr = 0.00001
> forc_conv_thr = 0.0001
> nstep = 400
> tstress = .true.
> tprnfor = .true.
> /
> &SYSTEM
> ibrav = -13
> celldm(1) = 15.795917182
> celldm(2) = 0.595028736
> celldm(3) = 0.7409843951
> celldm(5) = 0.4192777936
> nat = 12
> ntyp = 3
> ecutwfc = 60
> ecutrho = 600
> /
> &ELECTRONS
> electron_maxstep = 200
> conv_thr = 1.0D-11
> diago_thr_init = 1e-4
> startingpot = 'atomic'
> startingwfc = 'atomic'
> mixing_mode = 'plain'
> mixing_beta = 0.5
> mixing_ndim = 8
> diagonalization = 'david'
> /
> &IONS
> ion_dynamics = 'bfgs'
> /
>
> ATOMIC_SPECIES
> Li 7.016003437 li_pbe_v1.4.uspp.F.UPF
> C 12.0107000000 c_pbe_v1.2.uspp.F.UPF
> O 15.9994000000 o_pbe_v1.2.uspp.F.UPF
>
> ATOMIC_POSITIONS crystal
> Li 0.75080411 0.35557873 0.16417001
> Li 0.64442127 0.24919589 0.66417001
> Li 0.35557873 0.75080411 0.33582999
> Li 0.24919589 0.64442127 0.83582999
> C 0.93252757 0.93252757 0.75
> C 0.06747243 0.06747243 0.25
> O 0.67497784 0.67497784 0.75
> O 0.32502216 0.32502216 0.25
> O 0.21140108 0.91630278 0.687171
> O 0.08369722 0.78859892 0.187171
> O 0.91630278 0.21140108 0.812829
> O 0.78859892 0.08369722 0.312829
>
> K_POINTS automatic
> 8 8 6 0 0 0
> ——————————————————————————————————
> ph.in
>
> Li2CO3 phonon dispersion
> &INPUTPH
> prefix = 'Li2CO3'
> tr2_ph = 1.0d-14
> ldisp = .true.
> nq1 = 2
> nq2 = 2
> nq3 = 2
> outdir = './scf1'
> fildyn = 'Li2CO31.dyn'
> amass(1) = 7.016003437
> amass(2) = 12.0107000000
> amass(3) = 15.9994000000
> !recover = .true.
> /
> ——————————————————————————————————
> q2r.in
>
> &input
> fildyn = 'Li2CO31.dyn'
> zasr = 'crystal'
> flfrc = 'Li2CO31.fc'
> /
> ——————————————————————————————————
>
>
> Thanks,
> -Zeeshan
> --
> Zeeshan Ahmad
> PhD candidate, Mechanical Engineering
> Carnegie Mellon University
> https://www.andrew.cmu.edu/~azeeshan/
> <http://www.andrew.cmu.edu/user/azeeshan/>
>
>
> _______________________________________________
> Quantum ESPRESSO is supported by MaX (www.max-centre.eu/quantum-espresso)
> users mailing list users at lists.quantum-espresso.org
> https://lists.quantum-espresso.org/mailman/listinfo/users
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20200804/edc85d1c/attachment.html>
More information about the users
mailing list