9
1
mirror of https://github.com/QuantumPackage/qp2.git synced 2024-12-30 15:15:38 +01:00

added the possibility to unset frozen core orbitals

This commit is contained in:
eginer 2022-03-07 16:36:35 +01:00
parent 173ea9eab1
commit 980b48d906

View File

@ -7,12 +7,13 @@ setting all MOs as Active, except the n/2 first ones which are set as Core.
If pseudo-potentials are used, all the MOs are set as Active. If pseudo-potentials are used, all the MOs are set as Active.
Usage: Usage:
qp_set_frozen_core [-q|--query] [(-l|-s|--large|--small)] EZFIO_DIR qp_set_frozen_core [-q|--query] [(-l|-s|-u|--large|--small|--unset)] EZFIO_DIR
Options: Options:
-q --query Prints in the standard output the number of frozen MOs -q --query Prints in the standard output the number of frozen MOs
-l --large Use a small core -l --large Use a small core
-s --small Use a large core -s --small Use a large core
-u --unset Unset frozen core
Default numbers of frozen electrons: Default numbers of frozen electrons:
@ -88,7 +89,9 @@ def main(arguments):
elif charge <= 54: n_frozen += 9 elif charge <= 54: n_frozen += 9
elif charge <= 86: n_frozen += 18 elif charge <= 86: n_frozen += 18
elif charge <= 118: n_frozen += 27 elif charge <= 118: n_frozen += 27
elif arguments["--unset"]:
n_frozen = 0
else: # default else: # default
for charge in ezfio.nuclei_nucl_charge: for charge in ezfio.nuclei_nucl_charge:
if charge <= 4: pass if charge <= 4: pass