El-ph calculation on spin-polarized (LSDA) system
Moderators: myrta gruning, andrea marini, Daniele Varsano, Conor Hogan
-
- Posts: 9
- Joined: Sun Sep 17, 2017 2:09 pm
El-ph calculation on spin-polarized (LSDA) system
Dear Yambo Developer and Yambo users,
I have a problem in calculating electron-phonon coupling on spin-polarized material/LSDA (nspin = 2). Scf, nscf and ph.x calculation on the system are finished well, I faced difficulties in running ypp_ph (generating gkkp files), there is always an error message "PW(ELPH) databases ...[PHONON] ...incorrect K-point correspondance" , I have followed Dr. Alejandro's suggestion viewtopic.php?f=12&t=870 to change decimal precision on Kpoint, but the same problem still arose.
However, if I use non-polarized approach (nspin = 1), everything running well; unfortunately I got the wrong electronic structure.
Could you please to give me any suggestion how to perform electron-phonon calculation on the top of LSDA?
thank you
Best regards,
Faozan
I have a problem in calculating electron-phonon coupling on spin-polarized material/LSDA (nspin = 2). Scf, nscf and ph.x calculation on the system are finished well, I faced difficulties in running ypp_ph (generating gkkp files), there is always an error message "PW(ELPH) databases ...[PHONON] ...incorrect K-point correspondance" , I have followed Dr. Alejandro's suggestion viewtopic.php?f=12&t=870 to change decimal precision on Kpoint, but the same problem still arose.
However, if I use non-polarized approach (nspin = 1), everything running well; unfortunately I got the wrong electronic structure.
Could you please to give me any suggestion how to perform electron-phonon calculation on the top of LSDA?
thank you
Best regards,
Faozan
Faozan Ahmad
PhD Student, ITB, Indonesia
PhD Student, ITB, Indonesia
- amolina
- Posts: 135
- Joined: Fri Jul 15, 2011 11:23 am
- Location: Valencia, Spain
- Contact:
Re: El-ph calculation on spin-polarized (LSDA) system
Dear Faozan,
can you upload your scf, nscf and input files for the phonons. In principle, pw and yambo work for spin-polarized materials. Which version of pw are you using?
Best,
Alejandro
can you upload your scf, nscf and input files for the phonons. In principle, pw and yambo work for spin-polarized materials. Which version of pw are you using?
Best,
Alejandro
Alejandro Molina-Sánchez
Institute of Materials Science (ICMUV)
University of Valencia, Spain
Institute of Materials Science (ICMUV)
University of Valencia, Spain
-
- Posts: 9
- Joined: Sun Sep 17, 2017 2:09 pm
Re: El-ph calculation on spin-polarized (LSDA) system
Dear Dr. Alejandro,
Thank you for your prompt reply, Please find my attached input-file.
I am using QE-6.1 and yambo-4.1.4
Best regards,
Faozan
Thank you for your prompt reply, Please find my attached input-file.
I am using QE-6.1 and yambo-4.1.4
Best regards,
Faozan
You do not have the required permissions to view the files attached to this post.
Faozan Ahmad
PhD Student, ITB, Indonesia
PhD Student, ITB, Indonesia
- amolina
- Posts: 135
- Joined: Fri Jul 15, 2011 11:23 am
- Location: Valencia, Spain
- Contact:
Re: El-ph calculation on spin-polarized (LSDA) system
Dear Faozan,
your script seems OK for me. Maybe you can try to use a uniform grid as this one:
DVSCF input
&inputph
nq2 = 6
nq3 = 6
nq1 = 6
fildvscf = 'dvscf'
search_sym = .false.
qplot = .false.
fildyn = 'ferh.dyn'
iverbosity = 1
ldisp = .true.
alpha_mix(1) = 0.4
tr2_ph = 5e-11
prefix = 'ferh'
electron_phonon = 'dvscf'
trans = .true.
last_q = 10
start_q = 10
/
ELPH input file:
&inputph
nq2 = 6
nq3 = 6
nq1 = 6
fildvscf = 'dvscf'
search_sym = .false.
qplot = .false.
fildyn = 'ferh.dyn'
iverbosity = 1
ldisp = .true.
alpha_mix(1) = 0.4
tr2_ph = 5e-11
prefix = 'ferh'
electron_phonon = 'yambo'
trans = .false.
last_q = 10
start_q = 10
/
your script seems OK for me. Maybe you can try to use a uniform grid as this one:
DVSCF input
&inputph
nq2 = 6
nq3 = 6
nq1 = 6
fildvscf = 'dvscf'
search_sym = .false.
qplot = .false.
fildyn = 'ferh.dyn'
iverbosity = 1
ldisp = .true.
alpha_mix(1) = 0.4
tr2_ph = 5e-11
prefix = 'ferh'
electron_phonon = 'dvscf'
trans = .true.
last_q = 10
start_q = 10
/
ELPH input file:
&inputph
nq2 = 6
nq3 = 6
nq1 = 6
fildvscf = 'dvscf'
search_sym = .false.
qplot = .false.
fildyn = 'ferh.dyn'
iverbosity = 1
ldisp = .true.
alpha_mix(1) = 0.4
tr2_ph = 5e-11
prefix = 'ferh'
electron_phonon = 'yambo'
trans = .false.
last_q = 10
start_q = 10
/
Alejandro Molina-Sánchez
Institute of Materials Science (ICMUV)
University of Valencia, Spain
Institute of Materials Science (ICMUV)
University of Valencia, Spain
-
- Posts: 9
- Joined: Sun Sep 17, 2017 2:09 pm
Re: El-ph calculation on spin-polarized (LSDA) system
Dear Dr. Alejandro,
Thank you so much for checking my script, I'll try the input-file you suggest soon.
Best regards,
faozan
Thank you so much for checking my script, I'll try the input-file you suggest soon.
Best regards,
faozan
Faozan Ahmad
PhD Student, ITB, Indonesia
PhD Student, ITB, Indonesia
-
- Posts: 9
- Joined: Sun Sep 17, 2017 2:09 pm
Re: El-ph calculation on spin-polarized (LSDA) system
Dear Dr. Alejandro,
I have tried your suggested input, but still I got the same error. Here attached the input and output file.
Best regards,
Faozan
I have tried your suggested input, but still I got the same error. Here attached the input and output file.
Best regards,
Faozan
You do not have the required permissions to view the files attached to this post.
Faozan Ahmad
PhD Student, ITB, Indonesia
PhD Student, ITB, Indonesia
- amolina
- Posts: 135
- Joined: Fri Jul 15, 2011 11:23 am
- Location: Valencia, Spain
- Contact:
Re: El-ph calculation on spin-polarized (LSDA) system
Dear Faozan,
can you upload the scf and nscf files, together with the pseudos? I will try to reproduce the error.
Alejandro.
can you upload the scf and nscf files, together with the pseudos? I will try to reproduce the error.
Alejandro.
Alejandro Molina-Sánchez
Institute of Materials Science (ICMUV)
University of Valencia, Spain
Institute of Materials Science (ICMUV)
University of Valencia, Spain
-
- Posts: 9
- Joined: Sun Sep 17, 2017 2:09 pm
Re: El-ph calculation on spin-polarized (LSDA) system
Dear Dr. Alejandro,
Sorry for my delayed reply, I just read your post. Please find the attached files and thank you for your help.
regards,
Faozan
Sorry for my delayed reply, I just read your post. Please find the attached files and thank you for your help.
regards,
Faozan
You do not have the required permissions to view the files attached to this post.
Faozan Ahmad
PhD Student, ITB, Indonesia
PhD Student, ITB, Indonesia
- amolina
- Posts: 135
- Joined: Fri Jul 15, 2011 11:23 am
- Location: Valencia, Spain
- Contact:
Re: El-ph calculation on spin-polarized (LSDA) system
Hi Faozan,
have you tried without the flag search_sym = .false.?
Cheers,
Alejandro.
have you tried without the flag search_sym = .false.?
Cheers,
Alejandro.
Alejandro Molina-Sánchez
Institute of Materials Science (ICMUV)
University of Valencia, Spain
Institute of Materials Science (ICMUV)
University of Valencia, Spain
-
- Posts: 9
- Joined: Sun Sep 17, 2017 2:09 pm
Re: El-ph calculation on spin-polarized (LSDA) system
Hi Dr. Alejandro,
Yes, I have tried without the flag, but with random kpoint and still failed. do you successfully reproduce the same error?
Best,
Faozan
Yes, I have tried without the flag, but with random kpoint and still failed. do you successfully reproduce the same error?
Best,
Faozan
Faozan Ahmad
PhD Student, ITB, Indonesia
PhD Student, ITB, Indonesia