Post here questions linked with issue while running the EPW code
Moderator: stiwari
sevenel
Posts: 9 Joined: Mon Jul 19, 2021 8:12 am
Affiliation: postdcotor
Post
by sevenel » Fri Dec 17, 2021 1:38 pm
Dear all,
I want to calculate the mobility of 2D InSe in the paper Nano Lett. 2019, 19, 1774−1781,
This is my input file
Code: Select all
&inputepw
prefix='InSe'
outdir='./'
dvscf_dir='../save'
elph=.true.
epbwrite = .true.
epbread = .false.
epwwrite=.true.
epwread=.false.
lpolar = .true.
vme = 'wannier'
use_ws = .true.
wannierize=.true.
bands_skipped='exclude_bands = 1-12'
nbndsub=14
dis_froz_max = 4.0
dis_froz_min = -7.5
proj(1)='In:s,p'
proj(2)='Se:p'
efermi_read = .true.
fermi_energy= -1.4537
fsthick = 100
nstemp = 1
temps = 0.075
degaussw = 0.001
band_plot = .true.
filkf='KPOINTS1'
filqf='KPOINTS1'
nk1 = 12
nk2 = 12
nk3 = 1
nq1 = 12
nq2 = 12
nq3 = 1
/
and the output stop update after wannier90,
output stop update after wannier90 123.png (44.45 KiB) Viewed 6715 times
and when I go to the node of calculating, I found only part of the task is running (8/28),
only part of the task is running 456.png (59.09 KiB) Viewed 6715 times
anyone can give me some advice?
best wishes
zhenqing
hpaudya1
Posts: 194 Joined: Tue Mar 21, 2017 7:11 pm
Affiliation:
Post
by hpaudya1 » Sun Dec 19, 2021 6:29 pm
Hi Sevenel,
Please make sure your "dvscf_dir='../save'" directory contains *.dyn_q*, *.dvscf_q*, and *.phsave/patterns.*.
Best,
Hari
sevenel
Posts: 9 Joined: Mon Jul 19, 2021 8:12 am
Affiliation: postdcotor
Post
by sevenel » Mon Dec 20, 2021 1:14 am
Dear Hari,
Thanks for your advice, and I confirm that the ../save/ folder contains *.dyn_q*, *.dvscf_q*, and *.phsave/patterns.*.
Any other suggestions?
best wishes
zhenqing
hlee
Posts: 415 Joined: Thu Aug 03, 2017 12:24 pm
Affiliation: The University of Texas at Austin
Post
by hlee » Mon Dec 20, 2021 9:35 am
Dear zhenqing:
Symmetries of Bravais lattice: 24
Symmetries of crystal: 2
The information above doesn't make sense. Could you provide the inputs and full outputs of non-self-consistent and EPW calculations?
In the present and previous versions of EPW, some input flags of pw.x for the symmetry such as nosym, etc. are problematic.
Sincerely,
H. Lee
sevenel
Posts: 9 Joined: Mon Jul 19, 2021 8:12 am
Affiliation: postdcotor
Post
by sevenel » Tue Dec 21, 2021 1:35 am
Dear H. Lee,
Thanks for your advice, and I solved the problem with the following tips
In order to generate the q-point list from the symmetry of the system, we shouldn’t use nosym=.true. in nscf calculations. If nscf calculation generates more k points than you provided, you have to use calculation='bands' instead of calculation='nscf' without using any symmetry-related flag
such as nosym.
Best wishes
zhenqing
Nacy Cui
Posts: 10 Joined: Tue Dec 14, 2021 2:47 am
Affiliation: theoretical physics
Post
by Nacy Cui » Fri Feb 25, 2022 9:22 am
Dear zhenqing,
I encountered the same problem as you, but I didn't use the parameter- nosym in the process. And I also tried to change 'nscf' to 'bands' in non-self consistent computing, but the same problem still occurs. Do you have any good ways to share?
Best regards,
Xiangyue
hpaudya1
Posts: 194 Joined: Tue Mar 21, 2017 7:11 pm
Affiliation:
Post
by hpaudya1 » Fri Feb 25, 2022 2:51 pm
Hi xiangyue,
Could you provide the inputs and all outputs of nscf and EPW calculations?
Best,
Hari Paudyal
Nacy Cui
Posts: 10 Joined: Tue Dec 14, 2021 2:47 am
Affiliation: theoretical physics
Post
by Nacy Cui » Mon Feb 28, 2022 2:08 am
Dear Hari Paudyal,
The attached are my inputs and outputs files in the epw folder.
Best regards,
Xiangyue
Attachments
nscf.in.jpg (67.96 KiB) Viewed 6244 times
scf.in.jpg (64.26 KiB) Viewed 6244 times
&inputepw prefix = 'MoSi2N4', amass(1) = 95.9400024, amass(2) = 28.0849991, amass(3) = 14.0070000, outdir = './', dvscf_dir = './save', filukk = './MoSi2N4.ukk', elph = .true. ! calculate e-ph coefficients epwwrite = .true. ! write e-ph matrices in Wann representation epwread = .false. ! lpolar = .true. ! correct electron-phonon matrix for polar materials system_2d = .true. fsthick = 6 ! energy windows (eV) temps = 300 ! K degaussw = 0.01 ! Smearing in the energy-conserving delta functions (eV) degaussq = 0.02 ! Smearing for sum over q in the e-ph coupling (meV) elecselfen = .true. phonselfen = .false. elecselfen = .true. phonselfen = .false. wannierize = .true. ! calculate Wannier functions using W90 library nbndsub = 23, ! number of Wannier functions to utilize num_iter = 0, iprint = 2, dis_win_max = 17.0 dis_win_min = -18.0 dis_froz_min = -9.0 dis_froz_max = 6.8 proj(1) = 'Mo:l=2', proj(2) = 'Si:l=1' proj(3) = 'N:l=1' wdata(1) = 'bands_plot = .true.' wdata(2) = 'begin kpoint_path' wdata(3) = 'K -0.333 0.667 0.000 G 0.000 0.000 0.000' wdata(4) = 'G 0.000 0.000 0.000 M 0.000 0.500 0.000' wdata(5) = 'M 0.000 0.500 0.000 K -0.333 0.667 0.000' wdata(6) = 'end kpoint_path' wdata(7) = 'bands_plot_format = gnuplot' nkf1 = 128 nkf2 = 128 nkf3 = 128 nqf1 = 128 nqf2 = 128 nqf3 = 128 nk1 = 16 nk2 = 16 nk3 = 1 nq1 = 8 nq2 = 8 nq3 = 1 / 10 cartesian 0.000000000000000E+00 0.000000000000000E+00 0.000000000000000E+00 0.000000000000000E+00 0.144337567297408E+00 0.000000000000000E+00 0.000000000000000E+00 0.288675134594816E+00 0.000000000000000E+00 0.000000000000000E+00 0.433012701892223E+00 0.000000000000000E+00 0.000000000000000E+00 -0.577350269189631E+00 0.000000000000000E+00 0.125000000000010E+00 0.216506350946112E+00 0.000000000000000E+00 0.125000000000010E+00 0.360843918243519E+00 0.000000000000000E+00 0.125000000000010E+00 0.505181485540927E+00 0.000000000000000E+00 0.250000000000020E+00 0.433012701892223E+00 0.000000000000000E+00 0.250000000000020E+00 0.577350269189631E+00 0.000000000000000E+00 contents in the epw folder.jpg (176.34 KiB) Viewed 6244 times
hpaudya1
Posts: 194 Joined: Tue Mar 21, 2017 7:11 pm
Affiliation:
Post
by hpaudya1 » Mon Feb 28, 2022 2:07 pm
Hi Xiangyue,
Basically you can attach scf.in, scf.out, nscf.in, nscf.out, epw.in, and epw.out so that we can read them.
Best,
Hari Paudyal