On Fri, Nov 6, 2020 at 11:14 PM Ashis Kundu <ashiskundu...@gmail.com> wrote:
Why is there a mismatch between k-points from SCF and q-point from > system.dyn0 (phonon calculation)? > [...] Ideally they should be the same > they are the same, but the algorithms for generation and selection of the irreducible sector are different for the two grids so the ordering and sign may differ. Each of the points in the first set is either equal to or equivalent to a point of the second set. Please let me know if I am doing something wrong. > you aren't doing anything wrong. Or maybe you are: the scf input you sent cannot produce the k-point grid you listed Paolo > > The k-points generated from SCF > > number of k points= 8 Methfessel-Paxton smearing, width (Ry)= > 0.0200 > cart. coord. in units 2pi/alat > k( 1) = ( 0.0000000 0.0000000 0.0000000), wk = 0.0740741 > k( 2) = ( -0.3333333 0.0000000 0.0977788), wk = 0.2962963 > k( 3) = ( -0.3333333 0.3333333 0.1955577), wk = 0.5925926 > k( 4) = ( 0.3333333 0.3333333 -0.0000000), wk = 0.2962963 > k( 5) = ( 0.0000000 0.0000000 0.1955577), wk = 0.1481481 > k( 6) = ( 0.6666667 0.0000000 -0.0000000), wk = 0.1481481 > k( 7) = ( 0.0000000 -0.3333333 -0.0977788), wk = 0.2962963 > k( 8) = ( 0.0000000 0.6666667 0.0000000), wk = 0.1481481 > > q-point from phonon calculations (a mesh of 3x3x3 is used) > > 0.00000000 0.00000000 0.00000000 > -0.33333333 0.00000000 0.09777884 > 0.00000000 0.33333333 0.09777884 > -0.33333333 0.33333333 0.19555767 > 0.33333333 0.33333333 -0.00000000 > 0.00000000 0.00000000 0.19555767 > 0.66666667 0.00000000 -0.00000000 > 0.00000000 -0.66666667 -0.00000000 > > Below are the SCF and phonon input file > > scf.in > ----- > &control > calculation = 'scf' > prefix = 'TaP' > pseudo_dir = '../pp' > outdir = './' > wf_collect = .true. > tprnfor = .true. > tstress = .true. > / > &system > ibrav = 7 > celldm(1) = 6.173529611 > celldm(3) = 3.409054 > nat = 4 > ntyp = 2 > ecutwfc = 100 > occupations = 'smearing' > smearing = 'mp' > degauss = 0.02 > / > &electrons > diagonalization = 'david' > mixing_beta = 0.7 > conv_thr = 1.0d-10 > / > ATOMIC_SPECIES > Ta 180.948 Ta.pw-mt_fhi.UPF > P 30.974 P.pw-mt_fhi.UPF > ATOMIC_POSITIONS {alat} > Ta 0.0000000 0.0000000 -0.0044625 > Ta 0.0000000 0.5000000 0.8478074 > P 0.0000000 0.0000000 1.4178478 > P 0.5000000 0.0000000 0.5655748 > K_POINTS {automatic} > 15 15 15 0 0 0 > > ph.in > ---- > phonon > &inputph > prefix = 'TaP' > fildyn = 'TaP.dyn' > fildvscf = 'dvscf' > ldisp = .true. > epsil = .false. > nq1 = 3 > nq2 = 3 > nq3 = 3 > tr2_ph = 1.0d-14 > / > > > Best regards, > Ashis > _______________________________________________ > Quantum ESPRESSO is supported by MaX (www.max-centre.eu) > users mailing list users@lists.quantum-espresso.org > https://lists.quantum-espresso.org/mailman/listinfo/users -- Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche, Univ. Udine, via delle Scienze 208, 33100 Udine, Italy Phone +39-0432-558216, fax +39-0432-558222
_______________________________________________ Quantum ESPRESSO is supported by MaX (www.max-centre.eu) users mailing list users@lists.quantum-espresso.org https://lists.quantum-espresso.org/mailman/listinfo/users