Page 1 of 1
Mapping the superconducting energy gap on the Fermi surface in BZ
Posted: Mon Jun 28, 2021 3:03 pm
by YSW
Dear all,
I have a problem about fermi surface.
I see the tutorial vedio that explains how to use the scripts in conjunction with VESTA.But it is a fermi surface in real space.
How can I get fermi surface in BZ?
A superconducting energy gap fermi surface in BZ of LaBH8.
Roxana shows it in the youtube video: Superconductors and Eliashberg theory - Roxana Margine.(55:23)
but she have not teach me how to get it.
I need you help badly.
Regards
Re: Mapping the superconducting energy gap on the Fermi surface in BZ
Posted: Mon Jun 28, 2021 7:25 pm
by hlee
Dear YSW:
But it is a fermi surface in real space.
I don't think so; very naturally, according to the definition of Fermi surface it is in reciprocal space.
Sincerely,
H. Lee
Re: Mapping the superconducting energy gap on the Fermi surface in BZ
Posted: Mon Jun 28, 2021 10:33 pm
by YSW
Dear H. Lee：
Yes,you are right.But,if I want to limit the image in 1st BZ.
What I need to do ?
Just like this
Sincerely,
YSW
Re: Mapping the superconducting energy gap on the Fermi surface in BZ
Posted: Mon Jun 28, 2021 10:37 pm
by YSW
Dear H. Lee：
This image is invisible likely.
So,I provide the link of this image.
https://www.researchgate.net/publicatio ... -right.jpg
Sincerely,
YSW
Re: Mapping the superconducting energy gap on the Fermi surface in BZ
Posted: Tue Jun 29, 2021 4:07 am
by hlee
Dear YSW:
It seems that you want to plot the Fermi surface in the first Brillouin zone instead of primitive Brillouin zone.
With Vesta you can't do it and you need to use other programs with the suitable file format.
One of programs I know is the Fermisurfer (
https://fermisurfer.osdn.jp/), but you need to convert the output generated by EPW into the format Fermisurfer can support.
It is not difficult and please check the chapter 4 of its manual at
https://fermisurfer.osdn.jp/en/_build/l ... surfer.pdf .
Sincerely,
H. Lee
Re: Mapping the superconducting energy gap on the Fermi surface in BZ
Posted: Tue Jun 29, 2021 9:07 am
by YSW
Dear H. Lee：
Thank you very much for your help and advice.
Code: Select all
real(4) :: bvec1(3), bvec2(3), bvec3(3) ! Resiplocal lattice vector
integer :: nk1, nk2, nk3 ! k-grid of each direction
integer :: ishift ! 1 for shifted grid, 0 for unshifted grid.
integer :: nbnd ! The number of bands
real(4) :: eig(nk3,nk2,nk1,nbnd) ! energy
real(4) :: x(nk3,nk2,nk1,nbnd) ! matrix element
integer :: ik1, ik2, ik3, ibnd, fo
open(fo, file = "sample.fs")
write(fo,*) nk1, nk2, nk3
write(fo,*) ishift
write(fo,*) nbnd
write(fo,*) real(bvec1(1:3))
write(fo,*) real(bvec2(1:3))
write(fo,*) real(bvec3(1:3))
do ibnd = 1, nbnd
do ik1 = 1, nk1
do ik2 = 1, nk2
do ik3 = 1, nk3
write(fo,*) real(eig(ik3,ik2,ik1,ibnd))
end do
end do
end do
end do
do ibnd = 1, nbnd
do ik1 = 1, nk1
do ik2 = 1, nk2
do ik3 = 1, nk3
write(fo,*) real(x(ik3,ik2,ik1,ibnd))
end do
end do
end do
end do
close(fo)
In this code,how shouid I set the matrix element ? Which output generated by EPW is included it ?
If I want to plot the superconducting energy gap fermi surface.
I use Prefix.imag_aniso_gap_FS_010.00, and the energy part is set delta(0) [meV],is it right?
Thank you again.
Sincerely,
YSW
Re: Mapping the superconducting energy gap on the Fermi surface in BZ
Posted: Wed Jun 30, 2021 4:50 am
by hlee
Dear YSW:
delta(0) should go to the part of "real(x(ik3,ik2,ik1,ibnd))" and band energies should go to the part of real(eig(ik3,ik2,ik1,ibnd)) .
Sincerely,
H. Lee
Re: Mapping the superconducting energy gap on the Fermi surface in BZ
Posted: Wed Jun 30, 2021 8:41 am
by YSW
Dear H. Lee：
Thank you for your guidance.
But,I have to say,I am a green hand.and I have some problems about using this script.
How should I execute this script ?
Need I use Prefix.imag_aniso_gap_FS_010.00 and nscf.out ? If I need, How should I use the date with this script ?
Sincerely,
YSW
