You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
# Modify LMAXMIX if you have d or f electrons present. Note that if the user
# explicitly sets LMAXMIX in settings it will override this logic. Previously, this
# was only set if Hubbard U was enabled as per the VASP manual but following an
# investigation it was determined that this would lead to a significant difference
# between SCF -> NonSCF even without Hubbard U enabled. Thanks to Andrew Rosen for
# investigating and reporting.
if"LMAXMIX"notinincar_settings.keys():
# contains f-electrons
ifany(el.Z>56forelinstructure.composition):
incar["LMAXMIX"] =6
# contains d-electrons
elifany(el.Z>20forelinstructure.composition):
incar["LMAXMIX"] =4
However, this is based on Z value and not based on the presence or lack of s, p, d, or f electrons (see VASP manual). For instance, Sr has Z > 20 but does not have d electrons. We should be checking this based on the "block" of the electrons in the Structure object.
The text was updated successfully, but these errors were encountered:
Currently,
LMAXMIX
is set as follows:atomate2/src/atomate2/vasp/sets/base.py
Lines 915 to 927 in c2d098a
However, this is based on Z value and not based on the presence or lack of s, p, d, or f electrons (see VASP manual). For instance, Sr has Z > 20 but does not have d electrons. We should be checking this based on the "block" of the electrons in the
Structure
object.The text was updated successfully, but these errors were encountered: