Hi,
May I ask if it is possible to use a smearing type for metals that are not in the 'occopt' options, e.g. Marzari-Vanderbilt smearing? If I would have to change the code myself to use a different smearing type, what files do I need to change?
Thank you!
Best,
Ying Shi
Smearing types [SOLVED]
Moderator: bguster
Re: Smearing types
Dear Ying Shi,
Isn't it one of the occopt=4, 5 or 6?
Do you have a reference paper/thesis for the smearing you want (just to certify that it is not one already implemented)?
Best wishes,
Eric
Isn't it one of the occopt=4, 5 or 6?
Do you have a reference paper/thesis for the smearing you want (just to certify that it is not one already implemented)?
Best wishes,
Eric
Re: Smearing types
Dear Eric,
Thanks for your reply. I am referring to this paper:
https://doi.org/10.1103/PhysRevLett.82.3296
This particular form seems to be recommended compared to the forms in occopt=4 or occopt=5.
Best,
Ying Shi
Thanks for your reply. I am referring to this paper:
https://doi.org/10.1103/PhysRevLett.82.3296
This particular form seems to be recommended compared to the forms in occopt=4 or occopt=5.
Best,
Ying Shi
Re: Smearing types
Dear Ying Shi,
In the limit of what I understand, the PRL paper you mention would correspond to the cold smearing implemented within occopt 4/5, it is not the case?
Best wishes,
Eric
In the limit of what I understand, the PRL paper you mention would correspond to the cold smearing implemented within occopt 4/5, it is not the case?
Best wishes,
Eric
Re: Smearing types
Dear Eric,
My understanding is that the smearing types in the PRL paper, as well as in occopt=4 or 5 are all cold smearing but of different flavors. The cases of occopt=4 and occopt=5 correspond to the form in Marzari's thesis: http://theossrv1.epfl.ch/Main/Theses?ac ... s_1996.pdf (See eqn 4.45 in pg 134 for the broadened delta function). However there is a revised and hence more elegant version of cold smearing in the PRL paper shortly after (See eqn 1 in the PRL paper). The ones that are already implemented in Abinit may work as good as this in practice but I thought I would ask beforehand to confirm.
The other reason for my question is that I need to know the explicit expression of generalized entropy for my work. However I don't know it for the ones in occopt=4 or 5, while it is readily available in equation 2 for the form in the PRL paper. Therefore if you happen to know the generalized entropy expression, that'd be great since I can then simply switch to using occopt=4 or 5.
Sorry for the long message and thanks so much for your attention!
Best,
Ying Shi
My understanding is that the smearing types in the PRL paper, as well as in occopt=4 or 5 are all cold smearing but of different flavors. The cases of occopt=4 and occopt=5 correspond to the form in Marzari's thesis: http://theossrv1.epfl.ch/Main/Theses?ac ... s_1996.pdf (See eqn 4.45 in pg 134 for the broadened delta function). However there is a revised and hence more elegant version of cold smearing in the PRL paper shortly after (See eqn 1 in the PRL paper). The ones that are already implemented in Abinit may work as good as this in practice but I thought I would ask beforehand to confirm.
The other reason for my question is that I need to know the explicit expression of generalized entropy for my work. However I don't know it for the ones in occopt=4 or 5, while it is readily available in equation 2 for the form in the PRL paper. Therefore if you happen to know the generalized entropy expression, that'd be great since I can then simply switch to using occopt=4 or 5.
Sorry for the long message and thanks so much for your attention!
Best,
Ying Shi
Re: Smearing types [SOLVED]
Dear Yin Shi,
It sounds like it is not clear what is implemented exactly. In the following presentation (slide 53):
https://pitp.phas.ubc.ca/confs/sherbroo ... nanese.pdf
it is mentioned the PRL of Marzari for occopt=4/5 but if I grep this PRL reference in Abinit I find it in the bibliography but this reference is not used anywhere and in the documentation of occopt it is mentioned the Marzari's thesis...
For the management of occupancy, entropy, you can look the routines of the folder src/61_occeig.
In init_occ_ent.F90 you can find the following:
Does it corresponds to the thesis or the PRL (sorry, I'm not familiar with these)?
Best wishes,
Eric
It sounds like it is not clear what is implemented exactly. In the following presentation (slide 53):
https://pitp.phas.ubc.ca/confs/sherbroo ... nanese.pdf
it is mentioned the PRL of Marzari for occopt=4/5 but if I grep this PRL reference in Abinit I find it in the bibliography but this reference is not used anywhere and in the documentation of occopt it is mentioned the Marzari's thesis...
For the management of occupancy, entropy, you can look the routines of the folder src/61_occeig.
In init_occ_ent.F90 you can find the following:
Code: Select all
! Cold smearing of Marzari, two values of the "a" parameter being possible
! first value gives minimization of the bump
if(occopt==4)aa=-.5634
! second value gives monotonic occupation function
if(occopt==5)aa=-.8165
dsqrpi=1.0_dp/sqrt(pi)
do ii=0,nptsdiv2_def
xx=xgrid_prev(ii)
gauss=dsqrpi*exp(-xx**2)
smdfun_prev( ii,1)=gauss*(1.5_dp+xx*(-aa*1.5_dp+xx*(-1.0_dp+aa*xx)))
smdfun_prev(-ii,1)=gauss*(1.5_dp+xx*( aa*1.5_dp+xx*(-1.0_dp-aa*xx)))
end do
Best wishes,
Eric
Re: Smearing types
Dear Eric,
Thanks for all the information! The implementation corresponds to the one in the thesis, not the PRL paper. The expression of entropy is actually available in the thesis so that would work.
Thanks so much for all your replies.
Best,
Ying Shi
Thanks for all the information! The implementation corresponds to the one in the thesis, not the PRL paper. The expression of entropy is actually available in the thesis so that would work.
Thanks so much for all your replies.
Best,
Ying Shi