Anisotropic FBW ME equations using IR basis functions (EPW Summer School 2024 exercise)

[krishnendu]

Hello EPW Forum,

I am trying to run the Exercise 3 of the Hands-on session (Wed. 9) of the "Calculations of superconducting Tc" (EPW Summer School 2024)

I ran upto 4 th step and successfully ran the epw1.in

In the 5th step when I give the command
$ ln -s /epw-lambda/*fmt

I am getting the following:

ln: target '../epw1-lambda/selecq.fmt' is not a directory

So, I guess I am not able to create the symbolic link for all the fmt files.

Kindly let me know what should I do so that I can run the epw2.in

Thank you,
Regards,
Krishnendu

[Shashi]

Dear Krishnedu,

You might have missed the destination path while creating a soft link.

As shown in Exercise -3, the command is "ln -s ../epw1-lambda/*fmt ." You might have skipped . at the end.

Please check again.

Regards,
Shashi

[krishnendu]

Dear Shashi,

Thank you for your reply, It worked.
Please find the epw2.in script below which I made by following the steps given in the Exercise-3, 5th step. When I try to run the script I get the error:

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Error in routine epw_readin (1):
Bad line in namelist &inputepw: "" (error could be in the previous line)
%%%%%%%%%%%%%%%%%%%%%%%

Kindly let me know what I should do so that I can run the epw2.in and get the superconducting gap for Nb for the given temperatures.

Thank you,
Regards,
Krishnendu
--------------------------------------------
--
&inputepw
prefix = 'nb',
outdir = './'
dvscf_dir = '../phonons-krishnendu/save' ! directory where .dyn, .dvscf and prefix.phsave/patterns.xx.yy
! files obtained from phonon calculation are stored

ep_coupling = .false. ! run e-ph coupling calculation

elph = .false. ! calculate e-ph coefficients

epwwrite = .false. ! write e-ph matrices in the Wann representation

epwread = .true. ! read e-ph matrices from the 'prefix.epmatwp' file


etf_mem = 0 ! less I/O but more memory is required


vme = 'dipole'

nbndsub = 9

bands_skipped = 'exclude_bands = 1:4'

wannierize = .false.

num_iter = 0

dis_win_min = 10

dis_froz_min = 11

dis_froz_max= 28

proj(1) = 'Nb:s'

proj(2) = 'Nb:p'

proj(3) = 'Nb:d'

iverbosity = 0 ! if you want to get .frmsf files, change this to 2.


fsthick = 0.5 ! Fermi window thickness [eV]

degaussw = 0.1 ! smearing in the energy-conserving delta functions in [eV]

degaussq = 0.5 ! smearing for sum over q in the e-ph coupling in [meV]



ephwrite = .false. ! write ephmatXX, egnv, freq, and ikmap files in prefix.ephmat directory

eliashberg = .true. ! calculate Eliashberg spectral function

laniso = .true.

limag = .true.

fbw = .true.

muchem = .true. ! update the chemical potential at every temperature

gridsamp = 2 ! If gridsamp = 2, use sparse-ir sampling

nsiter = 500 ! convergence threshold for solving ME eqs. on imaginary axis

filirobj = '../../irobjs/ir_nlambda6_ndigit8.dat'

! the file containing the sparse samping points and the IR functions.

broyden_beta = -0.7 ! If this factor is negative, use the linear mixing instead Broyden mixing.

conv_thr_iaxis = 1.0d-3 ! convergence threshold for solving ME eqs. on imaginary axis

wscut = 1.0

temps = 0.2 6.0 10.0 12.0 13.0 13.5 ! temperature list

! nsiter = 500 ! number of self-consistent iterations when solving ME eqs.

muc = 0.24 ! effective Coulomb potential used in ME eqs.





nk1=4
nk2= 4 ! dimensions of the coarse electronic grid
nk3= 4


nq1 =4
nq2=4 ! dimensions of the coarse phonon grid
nq3=4

mp_mesh_k = .true. ! use irreduciable electronic fine mesh
nkf1 = 20
nkf2 = 20 ! dimensions of the fine electronic grid
nkf3 = 20


nqf1 = 20
nqf2 = 20 ! dimensions of the fine phonon grid
nqf3 = 20

/

[Shashi]

Dear Krishnendu,

I think you might not have the "ir_nlambda6_ndigit8.dat" file in the irobjs folder. I could see that only the file included in irobjs is "
ir_nlambda5_ndigit6.dat". If you want to generate "ir_nlambda6_ndigit8.dat", you need to generate it using the command

make irobjs

as mentioned in q-e/EPWirobjs/README.md.

Please let me know if it resolves your issue.

Regards,
Shashi

[krishnendu]

Dear Shashi,

Thank you for your reply. Yes I already had the "ir_nlambda6_ndigit8.dat" file in the irobjs folder. I obtained it from the docs.epw-code.org/doc/Inputs.html: below the gridsamp.

So, I guess the problem is somewhere else. Kindly help me.

Thank you,
Best regards,
Krishnendu

[Shashi]

Hi Krishnendu,

I could not reproduce the error. I tested the tutorial, and it ran successfully without any issues. I also tested your EPW input for the 5th step, and it worked fine. I am not sure what might be causing the problem on your end.

Are you using the latest version of EPW? If not, I would suggest updating, as the writing .ephmat files have changed slightly.

Regards,
Shashi

[krishnendu]

Dear Shashi,

Thank you for your reply and also for testing the inputs.
I am using EPW-v.5.8.1 which is compatible with my existing Quantum Espresso v.7.3.1

So, would you recommend installing EPW-6.0.0 ?

Thank you,
Best regards,
Krishnendu

[Shashi]

Hi Krishnedu,

Please use EPW-6.0.0

wget https://gitlab.com/epw/q-e/-/archive/EP ... 6.0.tar.gz

Regards,
Shashi