Dear EPW experts,
I'm relatively new to using EPW. I need to do superconductivity calculation using Eliashberg equations and I have following queries. Kindly forgive me if these questions are basic and already been answered. I tried but couldn't locate the answers for these in the forum.
I have already done electron and phonon bandstructures in QE and EPW. Now, when I need to solve Iso/anisotropic ME equations,
1. what values should be set for fsthick, degaussw, degaussq, wscut and muc tags ? (At present, I'm using values given in EPW school hands-on examples)
2. I see that, you say we must do convergence tests for above tags. However, I'm confused against what parameter should we do convergence. ( for ex: in
QE, it was against Total energy for ecutwfc, kmesh etc). In EPW, Is it the decay ( i.e. decay of H, D and g) ?
3. I have a server with single node, 128 cores, 256 GB of RAM. When running, "mpirun -np (...) epw.x -nk (..) < epw.in > epw.out", what number of processors and pools I must use? Also, How do I set value for max_memlt ?
I sincerely apologize for these lengthy questions. Hoping for some reply from you soon.
Regards
Ashwath
confusion regarding setting the values of eliashberg equation parameters and parallezation
Moderator: stiwari
Re: confusion regarding setting the values of eliashberg equation parameters and parallezation
Hi Ashwath,
Before jumping into the ME calculations, the most important step is the wannierization. Please make sure that your W90 bands fits well with the DFT and the spreads are reasonably small (< 5 Ang.), then perform the Block to Wann calculation. for dynamical matrices, etc in the coarse grid. Then it comes now to calculate el-ph matrices in the fine k/q meshes, and calculate ME equations.
For the ME equations calculations, if you see EPW/examples/mgb2/epw example, following inputs are set, I have explained them here.
fsthick = 0.4 ! eV (4 times the maximum phonon frequency)
degaussw = 0.10 ! eV (maximum phonon frequency)
degaussq = 0.5 ! meV (greater than (max. phonon freq)/nqstep, by default nqstep = 500, see https://docs.epw-code.org/doc/Inputs.html#nqstep)
wscut = 1.0 ! (10 times the max. ph. freq.)
And, yes of course, you need to run some convergence tests of these parameters and also for the k/q meshes with respect to the el-ph coupling parameter \lambda and superconducting gap \Delta.
For memory related question, please see https://forum.epw-code.org/viewtopic.php?t=5887
Happy EPWing,
Hari
Before jumping into the ME calculations, the most important step is the wannierization. Please make sure that your W90 bands fits well with the DFT and the spreads are reasonably small (< 5 Ang.), then perform the Block to Wann calculation. for dynamical matrices, etc in the coarse grid. Then it comes now to calculate el-ph matrices in the fine k/q meshes, and calculate ME equations.
For the ME equations calculations, if you see EPW/examples/mgb2/epw example, following inputs are set, I have explained them here.
fsthick = 0.4 ! eV (4 times the maximum phonon frequency)
degaussw = 0.10 ! eV (maximum phonon frequency)
degaussq = 0.5 ! meV (greater than (max. phonon freq)/nqstep, by default nqstep = 500, see https://docs.epw-code.org/doc/Inputs.html#nqstep)
wscut = 1.0 ! (10 times the max. ph. freq.)
And, yes of course, you need to run some convergence tests of these parameters and also for the k/q meshes with respect to the el-ph coupling parameter \lambda and superconducting gap \Delta.
For memory related question, please see https://forum.epw-code.org/viewtopic.php?t=5887
Happy EPWing,
Hari
Re: confusion regarding setting the values of eliashberg equation parameters and parallezation
Hi Ashwath,
I would like to add a few things related to your convergence parameters:
1. fsthick value needs to be converged for your system. For isotropic superconductivity calculations, you can vary the fsthick window and check the convergence of the Eliashberg spectral function and lambda by plotting the prefix.a2f file.
degaussw is the smearing parameter for electronic states, currently, set as degaussw = (1/4)* fsthick. This parameter is crucial for isotropic calculations and needs to be converged with respect to k/q grid sizes.
degaussq is the smearing for phonon frequencies, used for isotropic calculation (as Hari clarified).
muc is a semi-empirical Coulomb parameter, typically ranging between 0.1 - 0.2 values are used for most cases.
2. EPW performs two steps. In the first step, the electron-phonon matrix elements are interpolated from the coarse k/q grid (nscf k-grid and phonon q-grid) to the Wannier representation. Therefore, decay plots for electronic Hamiltonian (decay.H), phonon dynamical matrices (decay.D), el-ph matrix elements (decay .Re, decay.Rp) must be converged for accurate MLWF representation.
The second step involves interpolating from the Wannier representation (epwread=.true., epwrite=.false.) to fine k/q-grids for el-ph matrix elements. These el-ph matrix elements are used for superconductivity, transport, etc. calculations.
3. You can run an EPW calculation using 128 cores. Since the current EPW version supports only pool parallelization, use the following commnd. mpirun -np 128 epw.x -nk 128 < epw.in > epw.out
Please let us know if anything is not clear.
Regards,
Shashi
I would like to add a few things related to your convergence parameters:
1. fsthick value needs to be converged for your system. For isotropic superconductivity calculations, you can vary the fsthick window and check the convergence of the Eliashberg spectral function and lambda by plotting the prefix.a2f file.
degaussw is the smearing parameter for electronic states, currently, set as degaussw = (1/4)* fsthick. This parameter is crucial for isotropic calculations and needs to be converged with respect to k/q grid sizes.
degaussq is the smearing for phonon frequencies, used for isotropic calculation (as Hari clarified).
muc is a semi-empirical Coulomb parameter, typically ranging between 0.1 - 0.2 values are used for most cases.
2. EPW performs two steps. In the first step, the electron-phonon matrix elements are interpolated from the coarse k/q grid (nscf k-grid and phonon q-grid) to the Wannier representation. Therefore, decay plots for electronic Hamiltonian (decay.H), phonon dynamical matrices (decay.D), el-ph matrix elements (decay .Re, decay.Rp) must be converged for accurate MLWF representation.
The second step involves interpolating from the Wannier representation (epwread=.true., epwrite=.false.) to fine k/q-grids for el-ph matrix elements. These el-ph matrix elements are used for superconductivity, transport, etc. calculations.
3. You can run an EPW calculation using 128 cores. Since the current EPW version supports only pool parallelization, use the following commnd. mpirun -np 128 epw.x -nk 128 < epw.in > epw.out
Please let us know if anything is not clear.
Regards,
Shashi