Skip to content
Permalink
Branch: master
Find file Copy path
Find file Copy path
Fetching contributors…
Cannot retrieve contributors at this time
8958 lines (7757 sloc) 500 KB
!--------------------------------------------------------------------------------------------------!
! CP2K: A general program to perform molecular dynamics simulations !
! Copyright (C) 2000 - 2019 CP2K developers group !
!--------------------------------------------------------------------------------------------------!
! **************************************************************************************************
!> \brief function that build the dft section of the input
!> \par History
!> 10.2005 moved out of input_cp2k [fawzi]
!> \author fawzi
! **************************************************************************************************
MODULE input_cp2k_dft
USE bibliography, ONLY: &
Andermatt2016, Andreussi2012, Avezac2005, BaniHashemian2016, Becke1988b, Bengtsson1999, &
Blochl1995, Brelaz1979, Dewar1977, Dewar1985, Dudarev1997, Dudarev1998, Ehrhardt1985, &
Fattebert2002, Golze2017a, Golze2017b, Guidon2010, Heinzmann1976, Holmberg2017, &
Holmberg2018, Hunt2003, Iannuzzi2005, Iannuzzi2006, Iannuzzi2007, Kolafa2004, Krack2000, &
Krack2002, Kunert2003, Lippert1997, Lippert1999, Lu2004, Perdew1981, Repasky2002, &
Rocha2006, Schenter2008, Schiffmann2015, Stewart1982, Stewart1989, Stewart2007, Thiel1992, &
Tozer1996, VanVoorhis2015, VandeVondele2003, VandeVondele2005a, VandeVondele2005b, &
VandeVondele2006, Weber2008, Zhao1994
USE cp_output_handling, ONLY: add_last_numeric,&
cp_print_key_section_create,&
debug_print_level,&
high_print_level,&
low_print_level,&
medium_print_level,&
silent_print_level
USE cp_units, ONLY: cp_unit_to_cp2k
USE input_constants, ONLY: &
atomic_guess, becke_cutoff_element, becke_cutoff_global, broyden_type_1, &
broyden_type_1_explicit, broyden_type_1_explicit_ls, broyden_type_1_ls, broyden_type_2, &
broyden_type_2_explicit, broyden_type_2_explicit_ls, broyden_type_2_ls, casci_canonical, &
cdft_alpha_constraint, cdft_beta_constraint, cdft_charge_constraint, &
cdft_magnetization_constraint, cholesky_dbcsr, cholesky_inverse, cholesky_off, &
cholesky_reduce, cholesky_restore, core_guess, diag_block_davidson, diag_block_krylov, &
diag_filter_matrix, diag_ot, diag_standard, dmft_model, do_admm_aux_exch_func_bee, &
do_admm_aux_exch_func_default, do_admm_aux_exch_func_none, do_admm_aux_exch_func_opt, &
do_admm_aux_exch_func_pbex, do_admm_basis_projection, do_admm_blocked_projection, &
do_admm_blocking_purify_full, do_admm_charge_constrained_projection, &
do_admm_exch_scaling_merlot, do_admm_exch_scaling_none, do_admm_purify_cauchy, &
do_admm_purify_cauchy_subspace, do_admm_purify_mcweeny, do_admm_purify_mo_diag, &
do_admm_purify_mo_no_diag, do_admm_purify_none, do_admm_purify_none_dm, do_arnoldi, &
do_bch, do_cn, do_ddapc_constraint, do_ddapc_restraint, do_em, do_etrs, do_full_density, &
do_gapw_gcs, do_gapw_gct, do_gapw_log, do_loc_both, do_loc_crazy, do_loc_direct, &
do_loc_homo, do_loc_jacobi, do_loc_l1_norm_sd, do_loc_lumo, do_loc_max, do_loc_min, &
do_loc_none, do_lri_inv, do_lri_inv_auto, do_lri_opt_all, do_lri_opt_coeff, &
do_lri_opt_exps, do_lri_pseudoinv_diag, do_lri_pseudoinv_svd, do_method_am1, &
do_method_dftb, do_method_gapw, do_method_gapw_xc, do_method_gpw, do_method_lrigpw, &
do_method_mndo, do_method_mndod, do_method_ofgpw, do_method_pdg, do_method_pm3, &
do_method_pm6, do_method_pm6fm, do_method_pnnl, do_method_rigpw, do_method_rm1, &
do_method_xtb, do_pade, do_ppl_analytic, do_ppl_grid, do_pwgrid_ns_fullspace, &
do_pwgrid_ns_halfspace, do_pwgrid_spherical, do_s2_constraint, do_s2_restraint, &
do_se_is_kdso, do_se_is_kdso_d, do_se_is_slater, do_se_lr_ewald, do_se_lr_ewald_gks, &
do_se_lr_ewald_r3, do_se_lr_none, do_spin_density, do_taylor, ehrenfest, embed_diff, &
embed_fa, embed_grid_angstrom, embed_grid_bohr, embed_level_shift, embed_none, &
embed_quasi_newton, embed_resp, embed_steep_desc, eri_method_full_gpw, eri_method_gpw_ht, &
eri_operator_coulomb, eri_operator_erf, eri_operator_erfc, eri_operator_gaussian, &
eri_operator_yukawa, gaussian, general_roks, hf_model, high_spin_roks, history_guess, &
jacobian_fd1, jacobian_fd1_backward, jacobian_fd1_central, jacobian_fd2, &
jacobian_fd2_backward, kg_cholesky, kg_color_dsatur, kg_color_greedy, &
kg_ec_diagonalization, kg_ec_functional_harris, kg_tnadd_atomic, kg_tnadd_embed, &
kg_tnadd_embed_ri, kg_tnadd_none, ls_2pnt, ls_3pnt, ls_gold, ls_none, mao_basis_ext, &
mao_basis_orb, mao_basis_prim, mao_projection, mopac_guess, no_excitations, no_guess, &
numerical, oe_gllb, oe_lb, oe_none, oe_saop, oe_sic, op_loc_berry, op_loc_boys, &
op_loc_pipek, orb_dx2, orb_dxy, orb_dy2, orb_dyz, orb_dz2, orb_dzx, orb_px, orb_py, &
orb_pz, orb_s, ot_algo_irac, ot_algo_taylor_or_diag, ot_chol_irac, ot_lwdn_irac, &
ot_mini_broyden, ot_mini_cg, ot_mini_diis, ot_mini_sd, ot_poly_irac, ot_precond_full_all, &
ot_precond_full_kinetic, ot_precond_full_single, ot_precond_full_single_inverse, &
ot_precond_none, ot_precond_s_inverse, ot_precond_solver_default, &
ot_precond_solver_direct, ot_precond_solver_inv_chol, ot_precond_solver_update, &
outer_scf_basis_center_opt, outer_scf_becke_constraint, outer_scf_cdft_constraint, &
outer_scf_ddapc_constraint, outer_scf_hirshfeld_constraint, outer_scf_none, &
outer_scf_optimizer_bisect, outer_scf_optimizer_broyden, outer_scf_optimizer_diis, &
outer_scf_optimizer_newton, outer_scf_optimizer_newton_ls, outer_scf_optimizer_none, &
outer_scf_optimizer_sd, outer_scf_optimizer_secant, outer_scf_s2_constraint, &
plus_u_lowdin, plus_u_mulliken, plus_u_mulliken_charges, pw_interp, radius_covalent, &
radius_default, radius_single, radius_user, radius_vdw, random_guess, &
real_time_propagation, ref_charge_atomic, ref_charge_mulliken, rel_dkh, rel_none, &
rel_pot_erfc, rel_pot_full, rel_sczora_mp, rel_trans_atom, rel_trans_full, &
rel_trans_molecule, rel_zora, rel_zora_full, rel_zora_mp, restart_guess, rsdft_model, &
sccs_andreussi, sccs_derivative_cd3, sccs_derivative_cd5, sccs_derivative_cd7, &
sccs_derivative_fft, sccs_fattebert_gygi, shape_function_density, shape_function_gaussian, &
sic_ad, sic_eo, sic_list_all, sic_list_unpaired, sic_mauri_spz, sic_mauri_us, sic_none, &
slater, smear_energy_window, smear_fermi_dirac, smear_list, sparse_guess, &
spline3_nopbc_interp, spline3_pbc_interp, tddfpt_davidson, tddfpt_excitations, &
tddfpt_lanczos, tddfpt_singlet, tddfpt_triplet, use_coulomb, use_diff, use_no, &
use_restart_wfn, use_rt_restart, use_scf_wfn, wannier_projection, weight_type_mass, &
weight_type_unit, wfi_aspc_nr, wfi_frozen_method_nr, wfi_linear_p_method_nr, &
wfi_linear_ps_method_nr, wfi_linear_wf_method_nr, wfi_ps_method_nr, &
wfi_use_guess_method_nr, wfi_use_prev_p_method_nr, wfi_use_prev_rho_r_method_nr, &
wfi_use_prev_wf_method_nr, xas_1s_type, xas_2p_type, xas_2s_type, xas_3d_type, &
xas_3p_type, xas_3s_type, xas_4d_type, xas_4f_type, xas_4p_type, xas_4s_type, xas_dip_len, &
xas_dip_vel, xas_dscf, xas_none, xas_tp_fh, xas_tp_flex, xas_tp_hh, xas_tp_xfh, &
xas_tp_xhh, xes_tp_val
USE input_cp2k_almo, ONLY: create_almo_scf_section
USE input_cp2k_distribution, ONLY: create_distribution_section
USE input_cp2k_field, ONLY: create_efield_section,&
create_per_efield_section
USE input_cp2k_kpoints, ONLY: create_kpoints_section
USE input_cp2k_ls, ONLY: create_ls_scf_section
USE input_cp2k_mm, ONLY: create_dipoles_section,&
create_neighbor_lists_section
USE input_cp2k_motion_print, ONLY: add_format_keyword
USE input_cp2k_poisson, ONLY: create_poisson_section
USE input_cp2k_rsgrid, ONLY: create_rsgrid_section
USE input_cp2k_tb, ONLY: create_dftb_control_section,&
create_xtb_control_section
USE input_cp2k_transport, ONLY: create_transport_section
USE input_cp2k_xc, ONLY: create_xc_section
USE input_keyword_types, ONLY: keyword_create,&
keyword_release,&
keyword_type
USE input_section_types, ONLY: section_add_keyword,&
section_add_subsection,&
section_create,&
section_release,&
section_type
USE input_val_types, ONLY: char_t,&
integer_t,&
lchar_t,&
logical_t,&
real_t
USE kinds, ONLY: dp
USE pw_grids, ONLY: do_pw_grid_blocked_false,&
do_pw_grid_blocked_free,&
do_pw_grid_blocked_true
USE pw_spline_utils, ONLY: no_precond,&
precond_spl3_1,&
precond_spl3_2,&
precond_spl3_3,&
precond_spl3_aint,&
precond_spl3_aint2
USE qs_density_mixing_types, ONLY: create_mixing_section
USE qs_fb_input, ONLY: create_filtermatrix_section
USE qs_mom_types, ONLY: create_mom_section
USE string_utilities, ONLY: newline,&
s2a
#include "./base/base_uses.f90"
IMPLICIT NONE
PRIVATE
CHARACTER(len=*), PARAMETER, PRIVATE :: moduleN = 'input_cp2k_dft'
PUBLIC :: create_dft_section
PUBLIC :: create_bsse_section, create_qs_section
PUBLIC :: create_scf_section
PUBLIC :: create_interp_section, create_localize_section
PUBLIC :: create_ddapc_restraint_section
PUBLIC :: create_mgrid_section
CONTAINS
! **************************************************************************************************
!> \brief creates the dft section
!> \param section the section to be created
!> \author fawzi
! **************************************************************************************************
SUBROUTINE create_dft_section(section)
TYPE(section_type), POINTER :: section
CHARACTER(len=*), PARAMETER :: routineN = 'create_dft_section', &
routineP = moduleN//':'//routineN
TYPE(keyword_type), POINTER :: keyword
TYPE(section_type), POINTER :: subsection
CPASSERT(.NOT. ASSOCIATED(section))
CALL section_create(section, __LOCATION__, name="DFT", &
description="Parameter needed by LCAO DFT programs", &
n_keywords=3, n_subsections=4, repeats=.FALSE.)
NULLIFY (keyword)
CALL keyword_create(keyword, __LOCATION__, name="BASIS_SET_FILE_NAME", &
description="Name of the basis set file, may include a path", &
usage="BASIS_SET_FILE_NAME <FILENAME>", &
type_of_var=lchar_t, repeats=.TRUE., &
default_lc_val="BASIS_SET", n_var=1)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="POTENTIAL_FILE_NAME", &
description="Name of the pseudo potential file, may include a path", &
usage="POTENTIAL_FILE_NAME <FILENAME>", &
default_lc_val="POTENTIAL")
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="WFN_RESTART_FILE_NAME", &
variants=(/"RESTART_FILE_NAME"/), &
description="Name of the wavefunction restart file, may include a path."// &
" If no file is specified, the default is to open the file as generated by the wfn restart print key.", &
usage="WFN_RESTART_FILE_NAME <FILENAME>", &
type_of_var=lchar_t)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, &
name="UKS", &
variants=s2a("UNRESTRICTED_KOHN_SHAM", &
"LSD", &
"SPIN_POLARIZED"), &
description="Requests a spin-polarized calculation using alpha "// &
"and beta orbitals, i.e. no spin restriction is applied", &
usage="LSD", &
default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, &
name="ROKS", &
variants=(/"RESTRICTED_OPEN_KOHN_SHAM"/), &
description="Requests a restricted open Kohn-Sham calculation", &
usage="ROKS", &
default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, &
name="MULTIPLICITY", &
variants=(/"MULTIP"/), &
description="Two times the total spin plus one. "// &
"Specify 3 for a triplet, 4 for a quartet, "// &
"and so on. Default is 1 (singlet) for an "// &
"even number and 2 (doublet) for an odd number "// &
"of electrons.", &
usage="MULTIPLICITY 3", &
default_i_val=0) ! this default value is just a flag to get the above
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="CHARGE", &
description="The total charge of the system", &
usage="CHARGE -1", &
default_i_val=0)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="EXCITATIONS", &
description="If excitations should be calculated", &
usage="EXCITATIONS", &
enum_c_vals=s2a("NONE", "TDLR", "TDDFPT"), &
enum_i_vals=(/no_excitations, tddfpt_excitations, &
tddfpt_excitations/), &
default_i_val=no_excitations)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, &
name="PLUS_U_METHOD", &
description="Method employed for the calculation of the DFT+U contribution", &
repeats=.FALSE., &
enum_c_vals=s2a("LOWDIN", "MULLIKEN", "MULLIKEN_CHARGES"), &
enum_i_vals=(/plus_u_lowdin, plus_u_mulliken, plus_u_mulliken_charges/), &
enum_desc=s2a("Method based on Lowdin population analysis "// &
"(computationally expensive, since the diagonalization of the "// &
"overlap matrix is required, but possibly more robust than Mulliken)", &
"Method based on Mulliken population analysis using the net AO and "// &
"overlap populations (computationally cheap method)", &
"Method based on Mulliken gross orbital populations (GOP)"), &
n_var=1, &
default_i_val=plus_u_mulliken, &
usage="METHOD Lowdin")
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, &
name="RELAX_MULTIPLICITY", &
variants=(/"RELAX_MULTIP"/), &
description="Do not enforce the occupation of alpha and beta MOs due to the initially "// &
"defined multiplicity, but rather follow an Aufbau principle. "// &
"A threshold value greater than zero activates this option. "// &
"Larger threshold values increase the probability for a spin flip. "// &
"This option is only valid for unrestricted (i.e. spin polarised) "// &
"Kohn-Sham (UKS) calculations.", &
usage="RELAX_MULTIPLICITY 0.00001", &
repeats=.FALSE., &
default_r_val=0.0_dp)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="SUBCELLS", &
description="Read the grid size for subcell generation in the construction of "// &
"neighbor lists.", usage="SUBCELLS 1.5", &
n_var=1, default_r_val=2.0_dp)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="AUTO_BASIS", &
description="Specify size of automatically generated auxiliary basis sets: "// &
"Options={small,medium,large,huge}", &
usage="AUTO_BASIS {basis_type} {basis_size}", &
type_of_var=char_t, repeats=.TRUE., n_var=-1, default_c_vals=(/"X", "X"/))
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, &
name="SURFACE_DIPOLE_CORRECTION", &
variants=s2a("SURFACE_DIPOLE", &
"SURF_DIP"), &
description="For slab calculations with asymmetric geometries, activate the correction of "// &
" the electrostatic potential with "// &
"by compensating for the surface dipole. Implemented only for slabs with normal "// &
" parallel to one Cartesian axis. The normal direction is given by the keyword SURF_DIP_DIR ", &
usage="SURF_DIP", &
default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE., &
citations=(/Bengtsson1999/))
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, &
name="SURF_DIP_DIR", &
description="Cartesian axis parallel to surface normal.", &
enum_c_vals=s2a("X", "Y", "Z"), &
enum_i_vals=(/1, 2, 3/), &
enum_desc=s2a("Along x", "Along y", "Along z"), &
n_var=1, &
default_i_val=3, &
usage="SURF_DIP_DIR Z")
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
NULLIFY (subsection)
CALL create_scf_section(subsection)
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
CALL create_ls_scf_section(subsection)
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
CALL create_almo_scf_section(subsection)
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
CALL create_kg_section(subsection)
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
CALL create_admm_section(subsection)
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
CALL create_qs_section(subsection)
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
CALL create_tddfpt_section(subsection)
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
CALL create_mgrid_section(subsection, create_subsections=.TRUE.)
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
CALL create_xc_section(subsection)
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
CALL create_relativistic_section(subsection)
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
CALL create_sic_section(subsection)
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
CALL create_low_spin_roks_section(subsection)
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
CALL create_efield_section(subsection)
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
CALL create_per_efield_section(subsection)
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
CALL create_ext_pot_section(subsection)
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
CALL create_transport_section(subsection)
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
! ZMP sections to include the external density or v_xc potential
CALL create_ext_den_section(subsection)
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
CALL create_ext_vxc_section(subsection)
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
CALL create_poisson_section(subsection)
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
CALL create_kpoints_section(subsection)
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
CALL create_implicit_solv_section(subsection)
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
CALL create_density_fitting_section(subsection)
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
CALL create_xas_section(subsection)
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
CALL create_localize_section(subsection)
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
CALL create_rtp_section(subsection)
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
CALL create_print_dft_section(subsection)
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
CALL create_sccs_section(subsection)
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
END SUBROUTINE create_dft_section
! **************************************************************************************************
!> \brief Implicit Solvation Model
!> \param section ...
!> \author tlaino
! **************************************************************************************************
SUBROUTINE create_implicit_solv_section(section)
TYPE(section_type), POINTER :: section
CHARACTER(len=*), PARAMETER :: routineN = 'create_implicit_solv_section', &
routineP = moduleN//':'//routineN
TYPE(keyword_type), POINTER :: keyword
TYPE(section_type), POINTER :: print_key, subsection
NULLIFY (keyword, subsection, print_key)
CPASSERT(.NOT. ASSOCIATED(section))
CALL section_create(section, __LOCATION__, name="SCRF", &
description="Adds an implicit solvation model to the DFT calculation."// &
" Know also as Self Consistent Reaction Field.", &
n_keywords=0, n_subsections=0, repeats=.FALSE.)
CALL keyword_create(keyword, __LOCATION__, name="EPS_OUT", &
description="Value of the dielectric constant outside the sphere", &
usage="EPS_OUT <REAL>", &
default_r_val=1.0_dp)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="LMAX", &
description="Maximum value of L used in the multipole expansion", &
usage="LMAX <INTEGER>", &
default_i_val=3)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL create_sphere_section(subsection)
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
CALL cp_print_key_section_create(print_key, __LOCATION__, "program_run_info", &
description="Controls the printing basic info about the method", &
print_level=low_print_level, add_last=add_last_numeric, filename="__STD_OUT__")
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
END SUBROUTINE create_implicit_solv_section
! **************************************************************************************************
!> \brief Create Sphere cavity
!> \param section ...
!> \author tlaino
! **************************************************************************************************
SUBROUTINE create_sphere_section(section)
TYPE(section_type), POINTER :: section
CHARACTER(len=*), PARAMETER :: routineN = 'create_sphere_section', &
routineP = moduleN//':'//routineN
TYPE(keyword_type), POINTER :: keyword
TYPE(section_type), POINTER :: subsection
NULLIFY (keyword, subsection)
CPASSERT(.NOT. ASSOCIATED(section))
CALL section_create(section, __LOCATION__, name="SPHERE", &
description="Treats the implicit solvent environment like a sphere", &
n_keywords=0, n_subsections=0, repeats=.FALSE.)
CALL keyword_create(keyword, __LOCATION__, name="RADIUS", &
description="Value of the spherical cavity in the dielectric medium", &
usage="RADIUS <REAL>", &
unit_str="angstrom", &
type_of_var=real_t)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL create_center_section(subsection)
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
END SUBROUTINE create_sphere_section
! **************************************************************************************************
!> \brief ...
!> \param section ...
!> \author tlaino
! **************************************************************************************************
SUBROUTINE create_center_section(section)
TYPE(section_type), POINTER :: section
CHARACTER(len=*), PARAMETER :: routineN = 'create_center_section', &
routineP = moduleN//':'//routineN
TYPE(keyword_type), POINTER :: keyword
NULLIFY (keyword)
CPASSERT(.NOT. ASSOCIATED(section))
CALL section_create(section, __LOCATION__, name="CENTER", &
description="Defines the center of the sphere.", &
n_keywords=0, n_subsections=0, repeats=.FALSE.)
CALL keyword_create(keyword, __LOCATION__, name="XYZ", &
description="Coordinates of the center of the sphere", &
usage="XYZ <REAL> <REAL> <REAL>", &
unit_str="angstrom", &
type_of_var=real_t, n_var=3)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="ATOM_LIST", &
description="Defines a list of atoms to define the center of the sphere", &
usage="ATOM_LIST <INTEGER> .. <INTEGER>", &
type_of_var=integer_t, n_var=-1)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="WEIGHT_TYPE", &
description="Defines the weight used to define the center of the sphere"// &
" (if ATOM_LIST is provided)", &
usage="WEIGHT (UNIT|MASS)", &
enum_c_vals=(/"UNIT", "MASS"/), &
enum_i_vals=(/weight_type_unit, weight_type_mass/), &
default_i_val=weight_type_unit)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="FIXED", &
description="Specify if the center of the sphere should be fixed or"// &
" allowed to move", &
usage="FIXED <LOGICAL>", &
default_l_val=.TRUE.)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
END SUBROUTINE create_center_section
! **************************************************************************************************
!> \brief parameters fo the localization of wavefunctions
!> \param section ...
!> \par History
!> 03.2005 created [MI]
! **************************************************************************************************
SUBROUTINE create_localize_section(section)
TYPE(section_type), POINTER :: section
CHARACTER(len=*), PARAMETER :: routineN = 'create_localize_section', &
routineP = moduleN//':'//routineN
TYPE(keyword_type), POINTER :: keyword
TYPE(section_type), POINTER :: print_key, print_section, subsection
CPASSERT(.NOT. ASSOCIATED(section))
NULLIFY (keyword, print_key)
CALL section_create(section, __LOCATION__, name="LOCALIZE", &
description="Use one of the available methods to define the localization "// &
" and possibly to optimize it to a minimum or a maximum.", &
n_keywords=8, n_subsections=0, repeats=.FALSE.)
CALL keyword_create(keyword, __LOCATION__, name="_SECTION_PARAMETERS_", &
description="controls the activation of the MOS localization procedure", &
usage="&LOCALIZE T", default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="MAX_ITER", &
description="Maximum number of iterations used for localization methods", &
usage="MAX_ITER 2000", default_i_val=10000)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create( &
keyword, __LOCATION__, name="MAX_CRAZY_ANGLE", &
description="Largest allowed angle for the crazy rotations algorithm (smaller is slower but more stable).", &
usage="MAX_CRAZY_ANGLE 0.1", unit_str="rad", default_r_val=0.2_dp)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="CRAZY_SCALE", &
description="scale angles", &
usage="CRAZY_SCALE 0.9", default_r_val=1.0_dp)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="CRAZY_USE_DIAG", &
description="Use diagonalization (slow) or pade based calculation of matrix exponentials.", &
usage="CRAZY_USE_DIAG ", default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create( &
keyword, __LOCATION__, name="USE_HISTORY", &
description="Generate an improved initial guess based on a history of results, which is useful during MD."// &
"Will only work if the number of states to be localized remains constant.", &
usage="USE_HISTORY ", default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create( &
keyword, __LOCATION__, name="EPS_OCCUPATION", &
description="Tolerance in the occupation number to select only fully occupied orbitals for the rotation", &
usage="EPS_OCCUPATION 1.E-5", default_r_val=1.0E-8_dp)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="OUT_ITER_EACH", &
description="Every how many iterations of the localization algorithm"// &
"(Jacobi) the tolerance value is printed out", &
usage="OUT_ITER_EACH 100", default_i_val=100)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="EPS_LOCALIZATION", &
description="Tolerance used in the convergence criterion of the localization methods.", &
usage="EPS_LOCALIZATION 1.0E-2", default_r_val=1.0E-4_dp)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="MIN_OR_MAX", &
description="Requires the maximization of the spread of the wfn", &
usage="MIN_OR_MAX (SPREADMIN|SPREADMAX)", &
enum_c_vals=(/"SPREADMIN", "SPREADMAX"/), &
enum_i_vals=(/do_loc_min, do_loc_max/), &
default_i_val=do_loc_min)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create( &
keyword, __LOCATION__, name="METHOD", &
description="Method of optimization if any", &
usage="METHOD (JACOBI|CRAZY|DIRECT|L1SD|NONE)", &
enum_c_vals=s2a("NONE", "JACOBI", "CRAZY", "L1SD", "DIRECT"), &
enum_i_vals=(/do_loc_none, &
do_loc_jacobi, &
do_loc_crazy, &
do_loc_l1_norm_sd, &
do_loc_direct/), &
enum_desc=s2a("No localization is applied", &
"Using 2 x 2 rotations of the orbitals, slow but robust", &
"A new fast method is applied, might be slightly less robust than jacobi, but usually much faster", &
"Steepest descent minimization of an approximate l1 norm", &
"Using a direct minimisation approach"), &
default_i_val=do_loc_jacobi)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="JACOBI_FALLBACK", &
description="Use Jacobi method in case no convergence was achieved"// &
" by using the crazy rotations method.", &
usage="JACOBI_FALLBACK", default_l_val=.TRUE., &
lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="RESTART", &
description="Restart the localization from a set of orbitals"// &
" read from a localization restart file.", &
usage="RESTART", default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="LOCHOMO_RESTART_FILE_NAME", &
description="File name where to read the MOS from"// &
"which to restart the localization procedure for occupied states", &
usage="LOCHOMO_RESTART_FILE_NAME <FILENAME>", &
type_of_var=lchar_t)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="LOCLUMO_RESTART_FILE_NAME", &
description="File name where to read the MOS from"// &
"which to restart the localization procedure for unoccupied states", &
usage="LOCLUMO_RESTART_FILE_NAME <FILENAME>", &
type_of_var=lchar_t)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="OPERATOR", &
description="Type of opertator which defines the spread functional", &
usage="OPERATOR (BERRY|BOYS|PIPEK)", &
enum_c_vals=s2a("BERRY", "BOYS", "PIPEK"), &
enum_i_vals=(/op_loc_berry, op_loc_boys, op_loc_pipek/), &
default_i_val=op_loc_berry)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="LIST", &
description="Indexes of the occupied wfn to be localized"// &
"This keyword can be repeated several times"// &
"(useful if you have to specify many indexes).", &
usage="LIST 1 2", &
n_var=-1, type_of_var=integer_t, repeats=.TRUE.)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="LIST_UNOCCUPIED", &
description="Indexes of the unoccupied states to be localized, "// &
"up to now only valid in combination with GPW. "// &
"This keyword has to be present if unoccupied states should be localized. "// &
"This keyword can be repeated several times"// &
"(useful if you have to specify many indexes).", &
usage="LIST 1 2", &
n_var=-1, type_of_var=integer_t, repeats=.TRUE.)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="STATES", &
description="Which states to localize, LUMO up to now only available in GPW", &
usage="STATES (HOMO|LUMO|ALL)", &
enum_c_vals=s2a("OCCUPIED", "UNOCCUPIED", "ALL"), &
enum_i_vals=(/do_loc_homo, do_loc_lumo, do_loc_both/), &
default_i_val=do_loc_homo)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create( &
keyword, __LOCATION__, &
name="ENERGY_RANGE", &
description="Select the orbitals to be localized within the given energy range."// &
"This type of selection cannot be added on top of the selection through a LIST. It reads to reals that are"// &
" lower and higher boundaries of the energy range.", &
usage=" ENERGY_RANGE lower_bound {real}, higher_bound {real}", &
repeats=.FALSE., &
n_var=2, default_r_vals=(/0._dp, 0._dp/), unit_str='eV', &
type_of_var=real_t)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
NULLIFY (print_section)
CALL section_create(print_section, __LOCATION__, name="PRINT", &
description="Collects all printing options related to the Wannier centers and "// &
"properties computed with Wannier centers.", &
n_keywords=0, n_subsections=1, repeats=.FALSE.)
NULLIFY (print_key)
CALL cp_print_key_section_create(print_key, __LOCATION__, "program_run_info", &
description="Controls the printing basic info about the method", &
print_level=low_print_level, add_last=add_last_numeric, filename="__STD_OUT__")
CALL section_add_subsection(print_section, print_key)
CALL section_release(print_key)
! Add printing of wannier infos
CALL print_wanniers(print_section)
NULLIFY (subsection)
! Total Dipoles with wannier
CALL create_dipoles_section(subsection, "TOTAL_DIPOLE", debug_print_level+1)
CALL section_add_subsection(print_section, subsection)
CALL section_release(subsection)
! Molecular Dipoles with wannier
CALL create_dipoles_section(subsection, "MOLECULAR_DIPOLES", debug_print_level+1)
CALL section_add_subsection(print_section, subsection)
CALL section_release(subsection)
! Molecular States with wannier
CALL create_molecular_states_section(subsection)
CALL section_add_subsection(print_section, subsection)
CALL section_release(subsection)
! Wannier States with wannier
CALL create_wannier_states_section(subsection)
CALL section_add_subsection(print_section, subsection)
CALL section_release(subsection)
CALL section_add_subsection(section, print_section)
CALL section_release(print_section)
END SUBROUTINE create_localize_section
! **************************************************************************************************
!> \brief Controls the printing of the basic info coming from the LOCALIZE
!> section
!> \param section ...
!> \author teo
! **************************************************************************************************
SUBROUTINE print_wanniers(section)
TYPE(section_type), POINTER :: section
CHARACTER(len=*), PARAMETER :: routineN = 'print_wanniers', routineP = moduleN//':'//routineN
TYPE(keyword_type), POINTER :: keyword
TYPE(section_type), POINTER :: print_key
CPASSERT(ASSOCIATED(section))
NULLIFY (print_key, keyword)
CALL cp_print_key_section_create(print_key, __LOCATION__, "WANNIER_CUBES", &
description="Controls the printing of the wannier functions ", &
print_level=high_print_level, add_last=add_last_numeric, filename="")
CALL keyword_create(keyword, __LOCATION__, name="stride", &
description="The stride (X,Y,Z) used to write the cube file "// &
"(larger values result in smaller cube files). You can provide 3 numbers (for X,Y,Z) or"// &
" 1 number valid for all components.", &
usage="STRIDE 2 2 2", n_var=-1, default_i_vals=(/2, 2, 2/), type_of_var=integer_t)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="CUBES_LU_BOUNDS", &
variants=(/"CUBES_LU"/), &
description="The lower and upper index of the states to be printed as cube", &
usage="CUBES_LU_BOUNDS integer integer", &
n_var=2, default_i_vals=(/0, -2/), type_of_var=integer_t)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="CUBES_LIST", &
description="Indexes of the states to be printed as cube files"// &
"This keyword can be repeated several times"// &
"(useful if you have to specify many indexes).", &
usage="CUBES_LIST 1 2", &
n_var=-1, type_of_var=integer_t, repeats=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="APPEND", &
description="append the cube files when they already exist", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
NULLIFY (print_key)
CALL cp_print_key_section_create(print_key, __LOCATION__, "WANNIER_CENTERS", &
description="Controls the printing of the wannier functions", &
print_level=high_print_level, add_last=add_last_numeric, filename="", &
unit_str="angstrom")
CALL keyword_create(keyword, __LOCATION__, name="IONS+CENTERS", &
description="prints out the wannier centers together with the particles", &
usage="IONS+CENTERS", default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL add_format_keyword(keyword, print_key, pos=.TRUE., &
description="Specifies the format of the output file when IONS+CENTERS is enabled.")
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
NULLIFY (print_key)
CALL cp_print_key_section_create(print_key, __LOCATION__, "WANNIER_SPREADS", &
description="Controls the printing of the wannier functions", &
print_level=high_print_level, add_last=add_last_numeric, filename="")
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
NULLIFY (print_key)
CALL cp_print_key_section_create(print_key, __LOCATION__, "LOC_RESTART", &
description="Controls the printing of restart file for localized MOS", &
print_level=high_print_level, add_last=add_last_numeric, filename="")
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
END SUBROUTINE print_wanniers
! **************************************************************************************************
!> \brief Create the print dft section
!> \param section the section to create
!> \author teo
! **************************************************************************************************
SUBROUTINE create_print_dft_section(section)
TYPE(section_type), POINTER :: section
CHARACTER(len=*), PARAMETER :: routineN = 'create_print_dft_section', &
routineP = moduleN//':'//routineN
TYPE(keyword_type), POINTER :: keyword
TYPE(section_type), POINTER :: print_key, sub_print_key, subsection
CPASSERT(.NOT. ASSOCIATED(section))
CALL section_create(section, __LOCATION__, name="PRINT", &
description="Section of possible print options in DFT code.", &
n_keywords=0, n_subsections=1, repeats=.FALSE.)
NULLIFY (print_key, keyword, subsection)
CALL cp_print_key_section_create(print_key, __LOCATION__, "PROGRAM_BANNER", &
description="Controls the printing of the banner of the MM program", &
print_level=silent_print_level, filename="__STD_OUT__")
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
CALL cp_print_key_section_create(print_key, __LOCATION__, "BASIS_SET_FILE", &
description="Controls the printing of a file with all basis sets used.", &
print_level=high_print_level, filename="LOCAL_BASIS_SETS")
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
CALL cp_print_key_section_create(print_key, __LOCATION__, "KINETIC_ENERGY", &
description="Controls the printing of the kinetic energy", &
print_level=high_print_level, filename="__STD_OUT__")
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
CALL cp_print_key_section_create(print_key, __LOCATION__, "DERIVATIVES", &
description="Print all derivatives after the DFT calculation", &
print_level=high_print_level, filename="__STD_OUT__")
CALL keyword_create(keyword, __LOCATION__, &
name="ndigits", &
description="Specify the number of digits used to print derivatives", &
default_i_val=6)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
CALL cp_print_key_section_create(print_key, __LOCATION__, name="neighbor_lists", &
description="Controls the printing of the neighbor lists", &
print_level=debug_print_level, filename="", unit_str="angstrom")
CALL keyword_create(keyword, __LOCATION__, &
name="sab_orb", &
description="Activates the printing of the orbital "// &
"orbital neighbor lists, "// &
"i.e. the overlap neighbor lists", &
default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, &
name="sab_aux_fit", &
description="Activates the printing of the orbital "// &
"orbital neighbor lists wavefunction fitting basis, "// &
"i.e. the overlap neighbor lists", &
default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, &
name="sab_aux_fit_vs_orb", &
description="Activates the printing of the orbital "// &
"orbital mixed neighbor lists of wavefunction fitting basis, "// &
"and the orbital basis, i.e. the overlap neighbor lists", &
default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, &
name="sab_scp", &
description="Activates the printing of the vdW SCP "// &
"neighbor lists ", &
default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, &
name="sab_vdw", &
description="Activates the printing of the vdW "// &
"neighbor lists (from DFT, DFTB, SE), "// &
"i.e. the dispersion neighbor lists", &
default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, &
name="sab_cn", &
description="Activates the printing of the "// &
"neighbor lists used for coordination numbers in vdW DFT-D3", &
default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, &
name="sac_ae", &
description="Activates the printing of the orbital "// &
"nuclear attraction neighbor lists (erfc potential)", &
default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, &
name="sac_ppl", &
description="Activates the printing of the orbital "// &
"GTH-PPL neighbor lists (local part of the "// &
"Goedecker-Teter-Hutter pseudo potentials)", &
default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, &
name="sap_ppnl", &
description="Activates the printing of the orbital "// &
"GTH-PPNL neighbor lists (non-local part of the"// &
"Goedecker-Teter-Hutter pseudo potentials)", &
default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, &
name="sap_oce", &
description="Activates the printing of the orbital "// &
"PAW-projector neighbor lists (only GAPW)", &
default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, &
name="sab_se", &
description="Activates the printing of the two-center "// &
"neighbor lists for Coulomb type interactions in NDDO ", &
default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, &
name="sab_lrc", &
description="Activates the printing of the long-range SE correction "// &
"neighbor lists (only when doing long-range SE with integral scheme KDSO and KDSO-d)", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, &
name="sab_tbe", &
description="Activates the printing of the DFTB Ewald "// &
"neighbor lists ", &
default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, &
name="sab_core", &
description="Activates the printing of core interaction "// &
"neighbor lists ", &
default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, &
name="sab_xb", &
description="Activates the printing of XB interaction from (xTB) "// &
"neighbor lists ", &
default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, &
name="soo_list", &
description="Activates the printing of RI orbital-orbital "// &
"neighbor lists ", &
default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, &
name="sip_list", &
description="Activates the printing of RI basis-projector interaction "// &
"neighbor lists ", &
default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
CALL cp_print_key_section_create(print_key, __LOCATION__, "SUBCELL", &
description="Activates the printing of the subcells used for the"// &
"generation of neighbor lists.", unit_str="angstrom", &
print_level=high_print_level, filename="__STD_OUT__")
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
CALL cp_print_key_section_create(print_key, __LOCATION__, "AO_MATRICES", &
description="Controls the printing of the ao (i.e. contracted gaussian) matrices (debug).", &
print_level=debug_print_level, filename="__STD_OUT__")
CALL keyword_create(keyword, __LOCATION__, name="OMIT_HEADERS", &
description="Print only the matrix data, not the row and column headers", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="NDIGITS", &
description="Specify the number of digits used to print the AO matrices", &
default_i_val=6)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="CORE_HAMILTONIAN", &
description="If the printkey is activated controls the printing of the hamiltonian matrix", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="DENSITY", &
description="If the printkey is activated controls the printing of the density (P) matrix", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="KINETIC_ENERGY", &
description="If the printkey is activated controls the printing of the kinetic energy matrix", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="KOHN_SHAM_MATRIX", &
description="If the printkey is activated controls the printing of the kohn-sham matrix", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create( &
keyword, __LOCATION__, name="MATRIX_VXC", &
description="If the printkey is activated compute and print the matrix of the exchange and correlation potential."// &
"Only the GGA part for GPW is printed", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="ORTHO", &
description="If the printkey is activated controls the printing of the orthogonalization matrix", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="OVERLAP", &
description="If the printkey is activated controls the printing of the overlap matrix", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="COMMUTATOR_HR", &
description="Controls the printing of the [H,r] commutator matrix", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="FERMI_CONTACT", &
description="If the printkey is activated controls the printing of the Fermi contact matrix", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create( &
keyword, __LOCATION__, name="PSO", &
description="If the printkey is activated controls the printing of the paramagnetic spin-orbit matrices", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create( &
keyword, __LOCATION__, name="EFG", &
description="If the printkey is activated controls the printing of the electric field gradient matrices", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="POTENTIAL_ENERGY", &
description="If the printkey is activated controls the printing of the potential energy matrix", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="OCE_HARD", &
description="If the printkey is activated controls the printing of the OCE HARD matrix", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="OCE_SOFT", &
description="If the printkey is activated controls the printing of the OCE SOFT matrix", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="W_MATRIX", &
description="If the printkey is activated controls the printing of the w matrix", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="W_MATRIX_AUX_FIT", &
description="If the printkey is activated controls the printing of the w matrix", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="DERIVATIVES", &
description="If the printkey is activated controls the printing "// &
"of derivatives (for the matrixes that support this)", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
CALL cp_print_key_section_create( &
print_key, __LOCATION__, "MO", &
description="Controls the printing of the molecular orbitals. "// &
"Note that this is only functional with diagonalization based methods, in particular not with OT (see MO_CUBES)", &
print_level=high_print_level, filename="__STD_OUT__")
CALL keyword_create(keyword, __LOCATION__, name="Cartesian", &
description="If the printkey is activated controls the printing of the mo in the cartesian basis", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="EIGENVALUES", variants=s2a("EIGVALS"), &
description="If the printkey is activated controls the printing of the eigenvalues of the mos", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="EIGENVECTORS", variants=s2a("EIGVECS"), &
description="If the printkey is activated controls the printing of the eigenvectors of the mos", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="OCCUPATION_NUMBERS", variants=s2a("OCCNUMS"), &
description="If the printkey is activated controls the printing of the occupation numbers of the mos", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="NDIGITS", &
description="Specify the number of digits used to print the MO eigenvalues and occupation numbers", &
default_i_val=6)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, &
name="MO_INDEX_RANGE", &
variants=s2a("MO_RANGE", "RANGE"), &
description="Allows to print only a subset of the MO eigenvectors or eigenvalues. "// &
"The indices of the first and the last MO have to be specified", &
repeats=.FALSE., &
n_var=2, &
type_of_var=integer_t, &
default_i_vals=(/0, 0/), &
usage="MO_INDEX_RANGE 10 15")
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
CALL create_mo_cubes_section(print_key)
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
CALL create_stm_section(print_key)
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
CALL create_wfn_mix_section(subsection)
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
CALL create_active_space_section(subsection)
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
CALL section_create(subsection, __LOCATION__, name="GAPW", &
description="Controls the printing of some gapw related information (debug).", &
n_keywords=0, n_subsections=1, repeats=.FALSE.)
CALL cp_print_key_section_create(print_key, __LOCATION__, "projectors", &
description="If the printkey is activated controls if information on"// &
" the projectors is printed.", &
print_level=debug_print_level, filename="__STD_OUT__")
CALL section_add_subsection(subsection, print_key)
CALL section_release(print_key)
CALL cp_print_key_section_create(print_key, __LOCATION__, "rho0_information", &
description="If the printkey is activated controls if information on rho0 is printed.", &
print_level=debug_print_level, filename="__STD_OUT__", unit_str="angstrom")
CALL section_add_subsection(subsection, print_key)
CALL section_release(print_key)
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
CALL cp_print_key_section_create(print_key, __LOCATION__, "dft_control_parameters", &
description="Controls the printing of dft control parameters.", &
print_level=medium_print_level, filename="__STD_OUT__")
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
CALL cp_print_key_section_create(print_key, __LOCATION__, "KPOINTS", &
description="Controls the printing of kpoint information.", &
print_level=medium_print_level, filename="__STD_OUT__")
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
NULLIFY (subsection)
CALL create_bandstructure_section(subsection)
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
CALL cp_print_key_section_create(print_key, __LOCATION__, "OVERLAP_CONDITION", &
description="Controls the the checking and printing of and estimate "// &
"of the overlap matrix condition number", &
print_level=debug_print_level, filename="__STD_OUT__")
CALL keyword_create(keyword, __LOCATION__, name="1-NORM", &
description="Calculate an estimate of the 1-norm condition number", &
default_l_val=.TRUE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="DIAGONALIZATION", &
description="Calculate the 1- and 2-norm condition numbers using diagonalization", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="ARNOLDI", &
description="Calculate the 2-norm condition number using the Arnoldi code (may not be reliable)", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
CALL cp_print_key_section_create(print_key, __LOCATION__, name="E_DENSITY_CUBE", &
description="Controls the printing of cube files with "// &
"the electronic density and, for LSD "// &
"calculations, the spin density", &
print_level=high_print_level, filename="")
CALL keyword_create(keyword, __LOCATION__, name="stride", &
description="The stride (X,Y,Z) used to write the cube file "// &
"(larger values result in smaller cube files). You can provide 3 numbers (for X,Y,Z) or"// &
" 1 number valid for all components.", &
usage="STRIDE 2 2 2", n_var=-1, default_i_vals=(/2, 2, 2/), type_of_var=integer_t)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, &
name="TOTAL_DENSITY", &
description="Print the total electronic density in the case "// &
"of a GAPW run. This keyword has only an effect, "// &
"if PAW atoms are present. The default is to print "// &
"only the soft part of the electronic density and "// &
"to ignore the hard part. NOTE: The total density in "// &
"real space might exhibit unphysical features like "// &
"spikes due to the finite and thus truncated g vector "// &
"expansion", &
usage="TOTAL_DENSITY {logical}", &
repeats=.FALSE., &
n_var=1, &
default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="APPEND", &
description="append the cube files when they already exist", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="XRD_INTERFACE", &
description="It activates the print out of exponents and coefficients for the "// &
" Gaussian expansion of the core densities, based on atom calculations for each kind."// &
" The resulting core dansities are needed to compute the form factors."// &
" If GAPW the local densities are also given in terms of a Gaussian expansion,"// &
" by fitting the difference between local-fhard and local-soft density for each atom."// &
" In this case the keyword TOTAL_DENSITY is set to FALSE", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="NGAUSS", &
description="Number of Gaussian functions used in the expansion of atomic (core) density", &
usage="NGAUSS 10", n_var=1, default_i_val=12, type_of_var=integer_t)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
CALL cp_print_key_section_create(print_key, __LOCATION__, "tot_density_cube", &
description="Controls printing of cube files with "// &
"the total density (electrons+atomic core). Note that "// &
"the value of the total density is positive where the "// &
"electron density dominates and negative where the core is.", &
print_level=high_print_level, filename="")
CALL keyword_create(keyword, __LOCATION__, name="stride", &
description="The stride (X,Y,Z) used to write the cube file "// &
"(larger values result in smaller cube files). You can provide 3 numbers (for X,Y,Z) or"// &
" 1 number valid for all components.", &
usage="STRIDE 2 2 2", n_var=-1, default_i_vals=(/2, 2, 2/), type_of_var=integer_t)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="APPEND", &
description="append the cube files when they already exist", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
CALL cp_print_key_section_create(print_key, __LOCATION__, "v_hartree_cube", &
description="Controls the printing of a cube file with eletrostatic "// &
" potential generated by the total density (electrons+ions). It is "// &
" valid only for QS with GPW formalism."// &
" Note that by convention the potential has opposite sign than the expected physical one.", &
print_level=high_print_level, filename="")
CALL keyword_create(keyword, __LOCATION__, name="stride", &
description="The stride (X,Y,Z) used to write the cube file "// &
"(larger values result in smaller cube files). You can provide 3 numbers (for X,Y,Z) or"// &
" 1 number valid for all components.", &
usage="STRIDE 2 2 2", n_var=-1, default_i_vals=(/2, 2, 2/), type_of_var=integer_t)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="APPEND", &
description="append the cube files when they already exist", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
CALL cp_print_key_section_create(print_key, __LOCATION__, "external_potential_cube", &
description="Controls the printing of a cube file with external "// &
" potential from the DFT%EXTERNAL_POTENTIAL section only.", &
print_level=high_print_level, filename="")
CALL keyword_create(keyword, __LOCATION__, name="stride", &
description="The stride (X,Y,Z) used to write the cube file "// &
"(larger values result in smaller cube files). You can provide 3 numbers (for X,Y,Z) or"// &
" 1 number valid for all components.", &
usage="STRIDE 2 2 2", n_var=-1, default_i_vals=(/2, 2, 2/), type_of_var=integer_t)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="APPEND", &
description="append the cube files when they already exist", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
! cube files for data generated by the implicit (generalized) Poisson solver
CALL create_implicit_psolver_section(subsection)
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
! ZMP adding the print section for the v_xc cube
CALL cp_print_key_section_create(print_key, __LOCATION__, "v_xc_cube", &
description="Controls the printing of a cube file with xc "// &
" potential generated by the ZMP method (for the moment). It is "// &
" valid only for QS with GPW formalism .", &
print_level=high_print_level, filename="")
CALL keyword_create(keyword, __LOCATION__, name="stride", &
description="The stride (X,Y,Z) used to write the cube file "// &
"(larger values result in smaller cube files). You can provide 3 numbers (for X,Y,Z) or"// &
" 1 number valid for all components.", &
usage="STRIDE 2 2 2", n_var=-1, default_i_vals=(/2, 2, 2/), type_of_var=integer_t)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="APPEND", &
description="append the cube files when they already exist", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
CALL cp_print_key_section_create(print_key, __LOCATION__, "efield_cube", &
description="Controls the printing of cube files with electric "// &
" field generated by the total density (electrons+ions). It is "// &
" valid only for QS with GPW formalism .", &
print_level=high_print_level, filename="")
CALL keyword_create(keyword, __LOCATION__, name="stride", &
description="The stride (X,Y,Z) used to write the cube file "// &
"(larger values result in smaller cube files). You can provide 3 numbers (for X,Y,Z) or"// &
" 1 number valid for all components.", &
usage="STRIDE 2 2 2", n_var=-1, default_i_vals=(/2, 2, 2/), type_of_var=integer_t)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="APPEND", &
description="append the cube files when they already exist", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
CALL cp_print_key_section_create(print_key, __LOCATION__, "ELF_CUBE", &
description="Controls printing of cube files with "// &
"the electron localization function (ELF). Note that "// &
"the value of ELF is defined between 0 and 1: Pauli kinetic energy density normalized "// &
" by the kinetic energy density of a uniform el. gas of same density.", &
print_level=high_print_level, filename="")
CALL keyword_create(keyword, __LOCATION__, name="stride", &
description="The stride (X,Y,Z) used to write the cube file "// &
"(larger values result in smaller cube files). You can provide 3 numbers (for X,Y,Z) or"// &
" 1 number valid for all components.", &
usage="STRIDE 2 2 2", n_var=-1, default_i_vals=(/2, 2, 2/), type_of_var=integer_t)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="APPEND", &
description="append the cube files when they already exist", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="density_cutoff", &
description=" ", &
usage="density_cutoff 0.0001", &
repeats=.FALSE., &
n_var=1, &
type_of_var=real_t, &
default_r_val=1.0e-10_dp)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
CALL create_pdos_section(print_key)
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
CALL create_wannier_section(print_key)
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
!
!Printing of Moments
CALL create_dipoles_section(print_key, "MOMENTS", high_print_level)
CALL keyword_create( &
keyword, __LOCATION__, &
name="MAX_MOMENT", &
description="Maximum moment to be calculated. Values higher than 1 not implemented under periodic boundaries.", &
usage="MAX_MOMENT {integer}", &
repeats=.FALSE., &
n_var=1, &
type_of_var=integer_t, &
default_i_val=1)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, &
name="MAGNETIC", &
description="Calculate also magnetic moments, only implemented without periodic boundaries", &
usage="MAGNETIC yes", &
repeats=.FALSE., &
n_var=1, &
default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
! Mulliken population analysis
CALL cp_print_key_section_create(print_key, __LOCATION__, "MULLIKEN", &
description="Controls the printing of the Mulliken (spin) population analysis", &
print_level=medium_print_level, filename="__STD_OUT__", &
common_iter_levels=1)
CALL keyword_create( &
keyword, __LOCATION__, &
name="PRINT_GOP", &
description="Print the gross orbital populations (GOP) in addition to the gross atomic populations (GAP) "// &
"and net charges", &
usage="PRINT_GOP yes", &
repeats=.FALSE., &
n_var=1, &
default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create( &
keyword, __LOCATION__, &
name="PRINT_ALL", &
description="Print all information including the full net AO and overlap population matrix", &
usage="PRINT_ALL yes", &
repeats=.FALSE., &
n_var=1, &
default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
! Lowdin population analysis (fairly expensive to compute, so only at high)
CALL cp_print_key_section_create(print_key, __LOCATION__, "LOWDIN", &
description="Controls the printing of the Lowdin (spin) population analysis", &
print_level=high_print_level, filename="__STD_OUT__", &
common_iter_levels=1)
CALL keyword_create( &
keyword, __LOCATION__, &
name="PRINT_GOP", &
description="Print the orbital populations in addition to the atomic populations and net charges", &
usage="PRINT_GOP yes", &
repeats=.FALSE., &
n_var=1, &
default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create( &
keyword, __LOCATION__, &
name="PRINT_ALL", &
description="Print all information including the full symmetrically orthogonalised density matrix", &
usage="PRINT_ALL yes", &
repeats=.FALSE., &
n_var=1, &
default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
! Hirshfeld population analysis
CALL cp_print_key_section_create(print_key, __LOCATION__, "HIRSHFELD", &
description="Controls the printing of the Hirshfeld (spin) population analysis", &
print_level=medium_print_level, filename="__STD_OUT__", &
common_iter_levels=1)
CALL keyword_create(keyword, __LOCATION__, name="SELF_CONSISTENT", &
description="Calculate charges from the Hirscheld-I (self_consistent) method."// &
" This scales only the full shape function, not the added charge as in the original scheme.", &
usage="SELF_CONSISTENT yes", repeats=.FALSE., n_var=1, &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="SHAPE_FUNCTION", &
description="Type of shape function used for Hirshfeld partitioning.", &
usage="SHAPE_FUNCTION {Gaussian,Density}", repeats=.FALSE., n_var=1, &
default_i_val=shape_function_gaussian, &
enum_c_vals=s2a("GAUSSIAN", "DENSITY"), &
enum_desc=s2a("Single Gaussian with Colvalent radius", &
"Atomic density expanded in multiple Gaussians"), &
enum_i_vals=(/shape_function_gaussian, shape_function_density/))
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="REFERENCE_CHARGE", &
description="Charge of atomic partitioning function for Hirshfeld method.", &
usage="REFERENCE_CHARGE {Atomic,Mulliken}", repeats=.FALSE., n_var=1, &
default_i_val=ref_charge_atomic, &
enum_c_vals=s2a("ATOMIC", "MULLIKEN"), &
enum_desc=s2a("Use atomic core charges", "Calculate Mulliken charges"), &
enum_i_vals=(/ref_charge_atomic, ref_charge_mulliken/))
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="USER_RADIUS", &
description="Use user defined radii to generate Gaussians."// &
" These radii are defined by the keyword ATOMIC_RADII", &
usage="USER_RADIUS yes", repeats=.FALSE., n_var=1, &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="ATOMIC_RADII", &
description="Defines custom radii to setup the spherical Gaussians.", &
usage="ATOMIC_RADII {real} {real} {real}", repeats=.FALSE., &
unit_str="angstrom", &
type_of_var=real_t, n_var=-1)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
! MAO (modified atomic orbital) analysis
CALL cp_print_key_section_create(print_key, __LOCATION__, "MAO_ANALYSIS", &
description="Controls the printing of the MAO (modified atomic orbital) analysis", &
print_level=debug_print_level, filename="__STD_OUT__", &
common_iter_levels=1, &
citations=(/Heinzmann1976, Ehrhardt1985/))
CALL keyword_create(keyword, __LOCATION__, name="EPS_FILTER", &
description="Threshold for matrix elements in MAO determination.", &
usage="EPS_FILTER reps", repeats=.FALSE., n_var=1, &
default_r_val=1.e-8_dp, type_of_var=real_t)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="REFERENCE_BASIS", &
description="Basis set used to construct MAO's.", &
usage="REFERENCE_BASIS {ORBITAL,PRIMITIVE,EXTERNAL}", repeats=.FALSE., n_var=1, &
default_i_val=mao_basis_orb, &
enum_c_vals=s2a("ORBITAL", "PRIMITIVE", "EXTERNAL"), &
enum_desc=s2a("Use standard orbital basis set", "Construct basis from primitives of the orbital basis", &
"Read external basis (MAO)"), &
enum_i_vals=(/mao_basis_orb, mao_basis_prim, mao_basis_ext/))
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="PRINT_BASIS", &
description="Print out MAO reference basis.", &
usage="PRINT_BASIS {logical}", repeats=.FALSE., n_var=1, &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE., type_of_var=logical_t)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="EPS_GRAD", &
description="Threshold for gradient in MAO optimization.", &
usage="EPS_GRAD reps", repeats=.FALSE., n_var=1, &
default_r_val=1.e-4_dp, type_of_var=real_t)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="EPS_FUNCTION", &
description="Threshold for electron defect in MAO optimization.", &
usage="EPS_FUNCTION feps", repeats=.FALSE., n_var=1, &
default_r_val=1.e-3_dp, type_of_var=real_t)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="MAX_ITER", &
description="Maximum allowed iterations for MAO optimization.", &
usage="MAX_ITER iter", repeats=.FALSE., n_var=1, &
default_i_val=0, type_of_var=integer_t)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="NEGLECT_ABC", &
description="Neglect 3 atom terms in MAO analysis.", &
usage="NEGLECT_ABC {logical}", repeats=.FALSE., n_var=1, &
default_l_val=.TRUE., lone_keyword_l_val=.TRUE., type_of_var=logical_t)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="AB_THRESHOLD", &
description="Threshold for printing of AB shared electron numbers.", &
usage="AB_THRESHOLD thr", repeats=.FALSE., n_var=1, &
default_r_val=1.e-2_dp, type_of_var=real_t)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="ABC_THRESHOLD", &
description="Threshold for printing of ABC shared electron numbers.", &
usage="ABC_THRESHOLD thr", repeats=.FALSE., n_var=1, &
default_r_val=1.e-5_dp, type_of_var=real_t)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="ANALYZE_UNASSIGNED_CHARGE", &
description="Calculate atomic contributions to the unassigned charge.", &
usage="ANALYZE_UNASSIGNED_CHARGE {logical}", repeats=.FALSE., n_var=1, &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE., type_of_var=logical_t)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
!Minimal localized basis analysis
CALL cp_print_key_section_create(print_key, __LOCATION__, "MINBAS_ANALYSIS", &
description="Controls the printing of the minimal localized basis analysis", &
print_level=debug_print_level, filename="__STD_OUT__", &
common_iter_levels=1, &
citations=(/Lu2004/))
CALL keyword_create(keyword, __LOCATION__, name="EPS_FILTER", &
description="Threshold for matrix elements in basis determination.", &
usage="EPS_FILTER reps", repeats=.FALSE., n_var=1, &
default_r_val=1.e-8_dp, type_of_var=real_t)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="FULL_ORTHOGONALIZATION", &
description="Orthogonalize the localized minimal basis.", &
usage="FULL_ORTHOGONALIZATION {logical}", repeats=.FALSE., n_var=1, &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE., type_of_var=logical_t)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="BOND_ORDER", &
description="Calculate Mayer Bond Orders.", &
usage="BOND_ORDER {logical}", repeats=.FALSE., n_var=1, &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE., type_of_var=logical_t)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
NULLIFY (sub_print_key)
CALL cp_print_key_section_create(sub_print_key, __LOCATION__, "MINBAS_CUBE", &
description="Write the minimal basis on Cube files.", &
print_level=debug_print_level+1, add_last=add_last_numeric, filename="MINBAS")
CALL keyword_create(keyword, __LOCATION__, name="STRIDE", &
description="The stride (X,Y,Z) used to write the cube file "// &
"(larger values result in smaller cube files). You can provide 3 numbers (for X,Y,Z) or"// &
" 1 number valid for all components.", &
usage="STRIDE 2 2 2", n_var=-1, default_i_vals=(/2, 2, 2/), type_of_var=integer_t)
CALL section_add_keyword(sub_print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="ATOM_LIST", &
description="Indexes of the atoms minimal basis to be printed as cube files"// &
"This keyword can be repeated several times"// &
"(useful if you have to specify many indexes).", &
usage="CUBES_LIST 1 2", &
n_var=-1, type_of_var=integer_t, repeats=.TRUE.)
CALL section_add_keyword(sub_print_key, keyword)
CALL keyword_release(keyword)
CALL section_add_subsection(print_key, sub_print_key)
CALL section_release(sub_print_key)
NULLIFY (sub_print_key)
CALL cp_print_key_section_create(sub_print_key, __LOCATION__, "MOS_MOLDEN", &
description="Write the minimal basis in Molden file format, for visualisation.", &
print_level=debug_print_level+1, add_last=add_last_numeric, filename="MINBAS")
CALL keyword_create(keyword, __LOCATION__, name="NDIGITS", &
description="Specifies the number of signficiant digits retained. 3 is OK for visualization.", &
usage="NDIGITS {int}", &
default_i_val=3)
CALL section_add_keyword(sub_print_key, keyword)
CALL keyword_release(keyword)
CALL section_add_subsection(print_key, sub_print_key)
CALL section_release(sub_print_key)
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
!DOS from density matrix
CALL cp_print_key_section_create(print_key, __LOCATION__, "ENERGY_WINDOWS", &
description="Controls the printing of the DOS from the density matrix. "// &
"This allows the calculation of the DOS even in density matrix based"// &
"REAL_TIME_PROPAGATION and LS_SCF. "// &
"However, it requires a cubically scaling diagonalization of the Hamiltonian. "// &
"Hartree-Fock NYI, values will be wrong. "// &
"Careful, the orbitals in rtp/emd are not actually eigenstates of the Hamiltonian."// &
"Assumes absence of spin polarization (so far).", &
print_level=high_print_level, common_iter_levels=3, &
each_iter_names=s2a("MD"), each_iter_values=(/100/), &
add_last=add_last_numeric, filename="energy-windows")
CALL keyword_create(keyword, __LOCATION__, name="N_WINDOWS", &
description="The number of energy windows.", &
usage="N_WINDOWS 200", &
default_i_val=100)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="EPS_FILTER", &
description="Filtering threshold for sparse matrix operations.", &
usage="EPS_FILTER 1.0E-6", &
default_r_val=1.0E-14_dp)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="RESTRICT_RANGE", &
description="Restricts the energy windows to states close to the fermi level", &
usage="RESTRICT_RANGE .TRUE.", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="RANGE", &
description="If the RESTRICT_RANGE keyword is set, then all energy widnows will"// &
" be placed in an interval from from the fermi level minus to the fermi level plus this keyword", &
usage="RANGE 1", &
default_r_val=1.0_dp)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="PRINT_CUBES", &
description="Print the energy windows to cube files", &
usage="DENSITY_PROPAGATION .TRUE.", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="STRIDE", &
description="The stride (X,Y,Z) used to write the energy windows cube files (if enabled) "// &
"(larger values result in smaller cube files). You can provide 3 numbers (for X,Y,Z) or"// &
" 1 number valid for all components.", &
usage="STRIDE 2 2 2", n_var=-1, default_i_vals=(/2, 2, 2/), type_of_var=integer_t)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
! Hamiltonian in CSR format
CALL cp_print_key_section_create(print_key, __LOCATION__, "KS_CSR_WRITE", &
description="Write the KS matrix in CSR format into a file.", &
print_level=debug_print_level, filename="")
CALL keyword_create(keyword, __LOCATION__, name="Threshold", &
description="Threshold on the absolute value of the elements to be printed out. "// &
"In CP2K all the elements in a (atomic) matrix block are considered non-zero, "// &
"if the block contains at least one non-zero element.", &
usage="THRESHOLD {real}", &
repeats=.FALSE., &
default_r_val=0.0_dp)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="Upper_triangular", &
description="Print only the upper triangular part of the matrix. ", &
usage="UPPER_TRIANGULAR {logical}", &
repeats=.FALSE., &
default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="Binary", &
description="Whether or not to generate the file in a binary format. ", &
usage="BINARY {logical}", &
repeats=.FALSE., &
default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
! Overlap in CSR format
CALL cp_print_key_section_create(print_key, __LOCATION__, "S_CSR_WRITE", &
description="Write the overlap matrix in CSR format into a file.", &
print_level=debug_print_level, filename="")
CALL keyword_create(keyword, __LOCATION__, name="Threshold", &
description="Threshold on the absolute value of the elements to be printed out. "// &
"In CP2K all the elements in a (atomic) matrix block are considered non-zero, "// &
"if the block contains at least one non-zero element.", &
usage="THRESHOLD {real}", &
repeats=.FALSE., &
default_r_val=0.0_dp)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="Upper_triangular", &
description="Print only the upper triangular part of the matrix. ", &
usage="UPPER_TRIANGULAR {logical}", &
repeats=.FALSE., &
default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="Binary", &
description="Whether or not to generate the file in a binary format. ", &
usage="BINARY {logical}", &
repeats=.FALSE., &
default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
! interaction adjacency matrix
CALL cp_print_key_section_create(print_key, __LOCATION__, "ADJMAT_WRITE", &
description="Writes an (upper-triangular) adjacency matrix indicating the "// &
"interaction between atoms (according to overlapping basis functions). The "// &
"columns are: iatom, jatom, ikind, jkind; where iatom and jatom are the atom "// &
"indices (based on the coordinate file), ikind and jkind are the atomic kinds "// &
"(indeces as shown in the ATOMIC KIND INFORMATION section of a CP2K output). ", &
print_level=debug_print_level, filename="")
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
! Xray diffraction
CALL cp_print_key_section_create( &
print_key, __LOCATION__, name="XRAY_DIFFRACTION_SPECTRUM", &
description="Calculate and print the coherent X-ray "// &
"diffraction spectrum", &
print_level=debug_print_level, &
filename="", &
citations=(/Krack2000, Krack2002/))
CALL keyword_create( &
keyword, __LOCATION__, &
name="Q_MAX", &
variants=(/"Q_MAXIMUM"/), &
description="Maximum Q value calculated for the spectrum", &
usage="Q_MAX {real}", &
repeats=.FALSE., &
n_var=1, &
type_of_var=real_t, &
default_r_val=cp_unit_to_cp2k(value=20.0_dp, &
unit_str="angstrom^-1"), &
unit_str="angstrom^-1")
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
CALL cp_print_key_section_create(print_key, __LOCATION__, name="ELECTRIC_FIELD_GRADIENT", &
description="Calculate and print the electric field gradients "// &
"at atomic positions", &
print_level=debug_print_level, &
filename="__STD_OUT__")
CALL keyword_create(keyword, __LOCATION__, &
name="INTERPOLATION", &
description="Use interpolation method from real space grid", &
usage="INTERPOLATION {logical}", &
repeats=.FALSE., &
n_var=1, &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, &
name="GSPACE_SMOOTHING", &
description="Use a G-space smoothing function", &
usage="GSPACE_SMOOTHING cutoff {real}, width {real}", &
repeats=.FALSE., &
n_var=2, default_r_vals=(/-1._dp, -1._dp/), &
type_of_var=real_t)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, &
name="DEBUG", &
description="Print additional debug output", &
usage="DEBUG {logical}", &
repeats=.FALSE., &
n_var=1, &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL create_gspace_interp_section(subsection)
CALL section_add_subsection(print_key, subsection)
CALL section_release(subsection)
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
CALL cp_print_key_section_create(print_key, __LOCATION__, name="BASIS_MOLOPT_QUANTITIES", &
description="Print the two quantities needed in the basis molopt generation:"// &
" total energy and condition number of the overlap matrix (S matrix)", &
print_level=debug_print_level, &
filename="__STD_OUT__")
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
CALL cp_print_key_section_create(print_key, __LOCATION__, name="HYPERFINE_COUPLING_TENSOR", &
description="Calculate and print the EPR hyperfine coupling tensor"// &
" at atomic positions", &
print_level=debug_print_level, &
filename="__STD_OUT__")
CALL keyword_create(keyword, __LOCATION__, &
name="INTERACTION_RADIUS", &
description="Radius of interaction for EPR hyperfine tensor calculation", &
usage="INTERACTION_RADIUS radius {real}", &
repeats=.FALSE., &
n_var=1, default_r_val=10._dp, &
type_of_var=real_t)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
CALL cp_print_key_section_create(print_key, __LOCATION__, name="OPTIMIZE_LRI_BASIS", &
description="Optimize the exponents of the LRI basis set", &
print_level=low_print_level, &
filename="OPTIMIZED_LRI_BASIS")
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
CALL cp_print_key_section_create( &
print_key, __LOCATION__, name="PLUS_U", &
description="Controls the printing for the DFT+U methods", &
print_level=high_print_level, &
filename="__STD_OUT__", &
each_iter_names=s2a("QS_SCF"), &
each_iter_values=(/0/), &
citations=(/Dudarev1997, Dudarev1998/))
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
CALL cp_print_key_section_create( &
print_key, __LOCATION__, name="SCCS", &
description="Controls the printing for the SCCS models", &
print_level=high_print_level, &
filename="__STD_OUT__", &
each_iter_names=s2a("QS_SCF"), &
each_iter_values=(/0/), &
citations=(/Fattebert2002, Andreussi2012/))
NULLIFY (sub_print_key)
CALL cp_print_key_section_create( &
sub_print_key, __LOCATION__, name="DENSITY_GRADIENT", &
description="Controls the printing of the cube files with "// &
"the norm of the density gradient |&nabla;&rho;| "// &
"used by the SCCS model.", &
print_level=debug_print_level, &
filename="", &
each_iter_names=s2a("QS_SCF"), &
each_iter_values=(/0/))
CALL keyword_create(keyword, __LOCATION__, name="STRIDE", &
description="The stride (X,Y,Z) used to write the cube file "// &
"(larger values result in smaller cube files). You can provide 3 "// &
"numbers (for X,Y,Z) or 1 number valid for all components.", &
n_var=-1, &
default_i_vals=(/2, 2, 2/), &
type_of_var=integer_t, &
repeats=.FALSE.)
CALL section_add_keyword(sub_print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="APPEND", &
description="Append the cube files when they already exist", &
default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE., &
repeats=.FALSE.)
CALL section_add_keyword(sub_print_key, keyword)
CALL keyword_release(keyword)
CALL section_add_subsection(print_key, sub_print_key)
CALL section_release(sub_print_key)
CALL cp_print_key_section_create( &
sub_print_key, __LOCATION__, name="DIELECTRIC_FUNCTION", &
description="Controls the printing of the cube files with "// &
"the dielectric function used by the SCCS model. "// &
"This function determines the cavity formed by a solute in "// &
"a solvent and thus it can be used for the visualisaton of the cavity.", &
print_level=debug_print_level, &
filename="", &
each_iter_names=s2a("QS_SCF"), &
each_iter_values=(/0/), &
citations=(/Fattebert2002, Andreussi2012/))
CALL keyword_create(keyword, __LOCATION__, name="STRIDE", &
description="The stride (X,Y,Z) used to write the cube file "// &
"(larger values result in smaller cube files). You can provide 3 "// &
"numbers (for X,Y,Z) or 1 number valid for all components.", &
n_var=-1, &
default_i_vals=(/2, 2, 2/), &
type_of_var=integer_t, &
repeats=.FALSE.)
CALL section_add_keyword(sub_print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="APPEND", &
description="Append the cube files when they already exist", &
default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE., &
repeats=.FALSE.)
CALL section_add_keyword(sub_print_key, keyword)
CALL keyword_release(keyword)
CALL section_add_subsection(print_key, sub_print_key)
CALL section_release(sub_print_key)
CALL cp_print_key_section_create( &
sub_print_key, __LOCATION__, name="POLARISATION_POTENTIAL", &
description="Controls the printing of the cube files with the "// &
"polarisation potential &phi;<sup>pol</sup> used by the SCCS model with "// &
"total potential &phi;<sup>tot</sup> = &phi;<sup>sol</sup> + &phi;<sup>pol</sup>", &
print_level=debug_print_level, &
filename="", &
each_iter_names=s2a("QS_SCF"), &
each_iter_values=(/0/), &
citations=(/Fattebert2002, Andreussi2012/))
CALL keyword_create(keyword, __LOCATION__, name="STRIDE", &
description="The stride (X,Y,Z) used to write the cube file "// &
"(larger values result in smaller cube files). You can provide 3 "// &
"numbers (for X,Y,Z) or 1 number valid for all components.", &
n_var=-1, &
default_i_vals=(/2, 2, 2/), &
type_of_var=integer_t, &
repeats=.FALSE.)
CALL section_add_keyword(sub_print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="APPEND", &
description="Append the cube files when they already exist", &
default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE., &
repeats=.FALSE.)
CALL section_add_keyword(sub_print_key, keyword)
CALL keyword_release(keyword)
CALL section_add_subsection(print_key, sub_print_key)
CALL section_release(sub_print_key)
CALL section_add_subsection(section, print_key)
CALL section_release(print_key)
END SUBROUTINE create_print_dft_section
! **************************************************************************************************
!> \brief ...
!> \param section ...
!> \author JGH
! **************************************************************************************************
SUBROUTINE create_bandstructure_section(section)
TYPE(section_type), POINTER :: section
CHARACTER(len=*), PARAMETER :: routineN = 'create_bandstructure_section', &
routineP = moduleN//':'//routineN
TYPE(keyword_type), POINTER :: keyword
TYPE(section_type), POINTER :: subsection
CPASSERT(.NOT. ASSOCIATED(section))
CALL section_create(section, __LOCATION__, name="BAND_STRUCTURE", &
description="Specifies the k-points used in band structure calculation.", &
n_keywords=0, n_subsections=0, repeats=.FALSE.)
NULLIFY (keyword)
CALL keyword_create(keyword, __LOCATION__, name="FILE_NAME", &
description="File name used for band structure", &
usage="FILE_NAME <filename>", default_c_val="", &
type_of_var=char_t, n_var=1)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="ADDED_MOS", &
variants=(/"ADDED_BANDS"/), &
description="Number of MOs/Bands added to the Band Structure calculation.", &
default_i_val=0)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
NULLIFY (subsection)
CALL section_create(subsection, __LOCATION__, name="KPOINT_SET", &
description="Specifies a k-point line to be calculated.", &
n_keywords=0, n_subsections=0, repeats=.TRUE.)
! keywords
NULLIFY (keyword)
CALL keyword_create(keyword, __LOCATION__, name="SPECIAL_POINT", &
description="Name and coordinates of a special k-point", &
usage="SPECIAL_POINT GAMMA 0.0 0.0 0.0", n_var=-1, type_of_var=char_t, repeats=.TRUE.)
CALL section_add_keyword(subsection, keyword)
CALL keyword_release(keyword)
!
CALL keyword_create(keyword, __LOCATION__, name="NPOINTS", &
description="Number of k-points along the line.", &
default_i_val=0)
CALL section_add_keyword(subsection, keyword)
CALL keyword_release(keyword)
!
CALL keyword_create(keyword, __LOCATION__, name="UNITS", &
description="Special k-points are defined either in units"// &
" of reciprocal lattice vectors or in Cartesian coordinates in uints of 2Pi/len."// &
" B_VECTOR: in multiples of the reciprocal lattice vectors (b)."// &
" CART_ANGSTROM: In units of 2*Pi/Angstrom."// &
" CART_BOHR: In units of 2*Pi/Bohr.", &
usage="UNITS <value>", type_of_var=char_t, default_c_val="B_VECTOR")
CALL section_add_keyword(subsection, keyword)
CALL keyword_release(keyword)
!
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
END SUBROUTINE create_bandstructure_section
! **************************************************************************************************
!> \brief creates the input section for dealing with homo lumos, including dumping cubes
!> \param print_key ...
! **************************************************************************************************
SUBROUTINE create_mo_cubes_section(print_key)
TYPE(section_type), POINTER :: print_key
CHARACTER(len=*), PARAMETER :: routineN = 'create_mo_cubes_section', &
routineP = moduleN//':'//routineN
TYPE(keyword_type), POINTER :: keyword
NULLIFY (keyword)
CALL cp_print_key_section_create(print_key, __LOCATION__, "MO_CUBES", &
description="Controls the printing of cubes of the molecular orbitals.", &
print_level=high_print_level, filename="")
CALL keyword_create(keyword, __LOCATION__, name="stride", &
description="The stride (X,Y,Z) used to write the cube file "// &
"(larger values result in smaller cube files). You can provide 3 numbers (for X,Y,Z) or"// &
" 1 number valid for all components.", &
usage="STRIDE 2 2 2", n_var=-1, default_i_vals=(/2, 2, 2/), type_of_var=integer_t)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="write_cube", &
description="If the MO cube file should be written. If false, the eigenvalues are still computed."// &
" Can also be useful in combination with STM calculations", &
default_l_val=.TRUE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="nlumo", &
description="If the printkey is activated controls the number of lumos"// &
" that are printed and dumped as a cube (-1=all)", &
default_i_val=0)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create( &
keyword, __LOCATION__, name="nhomo", &
description="If the printkey is activated controls the number of homos that dumped as a cube (-1=all),"// &
" eigenvalues are always all dumped", &
default_i_val=1)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create( &
keyword, __LOCATION__, name="homo_list", &
description="If the printkey is activated controls the index of homos dumped as a cube,"// &
" eigenvalues are always all dumped. It overrides nhomo.", &
usage="HOMO_LIST {integer} {integer} .. {integer} ", type_of_var=integer_t, &
n_var=-1, repeats=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="APPEND", &
description="append the cube files when they already exist", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
END SUBROUTINE create_mo_cubes_section
! **************************************************************************************************
!> \brief ...
!> \param print_key ...
! **************************************************************************************************
SUBROUTINE create_pdos_section(print_key)
TYPE(section_type), POINTER :: print_key
CHARACTER(len=*), PARAMETER :: routineN = 'create_pdos_section', &
routineP = moduleN//':'//routineN
TYPE(keyword_type), POINTER :: keyword
TYPE(section_type), POINTER :: subsection
NULLIFY (subsection)
NULLIFY (keyword)
CALL cp_print_key_section_create(print_key, __LOCATION__, "PDOS", &
description="Print out the DOS projected per kind and per angular momentum ", &
print_level=debug_print_level, common_iter_levels=1, filename="")
CALL keyword_create(keyword, __LOCATION__, name="COMPONENTS", &
description="Print out pdos distinguishing all angular momentum components.", &
usage="COMPONENTS", default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="APPEND", &
description="Append the pdos obtained at different iterations to the pdos output file."// &
"By defaut the file is overwritten", &
usage="APPEND", default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create( &
keyword, __LOCATION__, name="NLUMO", &
description="Number of virtual orbitals to be added to the MO set (-1=all)."//newline// &
"CAUTION: Setting this value to be higher than the number of states present may cause a Cholesky error.", &
usage="NLUMO integer", default_i_val=0)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="OUT_EACH_MO", &
description="Output on the status of the calculation every OUT_EACH_MO states. If -1 no output", &
usage="OUT_EACH_MO integer", default_i_val=-1)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
!
CALL section_create(subsection, __LOCATION__, name="LDOS", &
description="Controls the printing of local PDOS, projected on subsets"// &
" of atoms given through lists", &
n_keywords=4, n_subsections=0, repeats=.TRUE.)
CALL keyword_create(keyword, __LOCATION__, name="COMPONENTS", &
description="Print out pdos distinguishing all angular momentum components.", &
usage="COMPONENTS", default_l_val=.FALSE., &
lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(subsection, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="LIST", &
description="Specifies a list of indexes of atoms where to project the DOS ", &
usage="LIST {integer} {integer} .. {integer} ", type_of_var=integer_t, &
n_var=-1, repeats=.TRUE.)
CALL section_add_keyword(subsection, keyword)
CALL keyword_release(keyword)
CALL section_add_subsection(print_key, subsection)
CALL section_release(subsection)
CALL section_create(subsection, __LOCATION__, name="R_LDOS", &
description="Controls the printing of local PDOS, projected on 3D volume in real space,"// &
" the volume is defined in terms of position with respect to atoms in the lists", &
n_keywords=4, n_subsections=0, repeats=.TRUE.)
CALL keyword_create(keyword, __LOCATION__, name="LIST", &
description="Specifies a list of indexes of atoms used to define the real space volume ", &
usage="LIST {integer} {integer} .. {integer} ", type_of_var=integer_t, &
n_var=-1, repeats=.TRUE.)
CALL section_add_keyword(subsection, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="XRANGE", &
description="range of positions in Cartesian direction x: all grid points within "// &
" this range from at least one atom of the list are considered", &
usage="XRANGE -10.0 10.0", unit_str="angstrom", n_var=2, type_of_var=real_t)
CALL section_add_keyword(subsection, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="YRANGE", &
description="range of positions in Cartesian direction y: all grid points within "// &
" this range from at least one atom of the list are considered", &
usage="YRANGE -10.0 10.0", unit_str="angstrom", n_var=2, type_of_var=real_t)
CALL section_add_keyword(subsection, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="ZRANGE", &
description="range of positions in Cartesian direction z: all grid points within "// &
" this range from at least one atom of the list are considered", &
usage="ZRANGE -10.0 10.0", unit_str="angstrom", n_var=2, type_of_var=real_t)
CALL section_add_keyword(subsection, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="ERANGE", &
description="only project states with the eigenvalues in the given interval. "// &
"Default is all states.", &
usage="ERANGE -1.0 1.0", unit_str="hartree", n_var=2, type_of_var=real_t)
CALL section_add_keyword(subsection, keyword)
CALL keyword_release(keyword)
CALL section_add_subsection(print_key, subsection)
CALL section_release(subsection)
END SUBROUTINE create_pdos_section
! **************************************************************************************************
!> \brief ...
!> \param print_key ...
! **************************************************************************************************
SUBROUTINE create_wannier_section(print_key)
TYPE(section_type), POINTER :: print_key
CHARACTER(len=*), PARAMETER :: routineN = 'create_wannier_section', &
routineP = moduleN//':'//routineN
TYPE(keyword_type), POINTER :: keyword
NULLIFY (keyword)
CALL cp_print_key_section_create(print_key, __LOCATION__, "WANNIER90", &
description="Interface to Wannier90 code. (EXPERIMENTAL)", &
print_level=debug_print_level, common_iter_levels=1, filename="")
CALL keyword_create(keyword, __LOCATION__, name="SEED_NAME", &
description="The seedname for the Wannier90 calculation (body of filenames).", &
usage="SEED_NAME filename", default_c_val="wannier90", &
n_var=1, type_of_var=char_t)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="MP_GRID", &
description="The dimensions of the Monkhorst-Pack k-point grid. ", &
usage="MP_GRID 6 6 6", n_var=-1, default_i_vals=(/10, 10, 10/), type_of_var=integer_t)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="ADDED_MOS", &
variants=(/"ADDED_BANDS"/), &
description="Number of MOs/Bands added to the Band Structure calculation.", &
default_i_val=0)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="EXCLUDE_BANDS", &
description="List of Bands excluded in the Wannier calculation.", &
usage="EXCLUDE_BANDS b1 b2 ...", n_var=-1, repeats=.TRUE., &
type_of_var=integer_t)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="WANNIER_FUNCTIONS", &
description="Number of Wannier functions to be calculated. ", &
usage="WANNIER_FUNCTIONS 6", n_var=1, default_i_val=0, &
repeats=.TRUE., type_of_var=integer_t)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
END SUBROUTINE create_wannier_section
! **************************************************************************************************
!> \brief ...
!> \param print_key ...
! **************************************************************************************************
SUBROUTINE create_stm_section(print_key)
TYPE(section_type), POINTER :: print_key
CHARACTER(len=*), PARAMETER :: routineN = 'create_stm_section', &
routineP = moduleN//':'//routineN
TYPE(keyword_type), POINTER :: keyword
NULLIFY (keyword)
CALL cp_print_key_section_create(print_key, __LOCATION__, "STM", &
description="Controls the printing of cubes for the generation of STM images.", &
print_level=debug_print_level, filename="")
CALL keyword_create(keyword, __LOCATION__, name="stride", &
description="The stride (X,Y,Z) used to write the cube file "// &
"(larger values result in smaller cube files). You can provide 3 numbers (for X,Y,Z) or"// &
" 1 number valid for all components.", &
usage="STRIDE 2 2 2", n_var=-1, default_i_vals=(/2, 2, 2/), type_of_var=integer_t)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="nlumo", &
description="If the printkey is activated controls the number of additional lumos"// &
" that are computed to be able to reproduce STM images obtained"// &
" from positive bias (imaging unoccupied states)", &
default_i_val=0)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="BIAS", &
description="Bias energy for scanning tunneling microscopy (STM) image generation."// &
"Orbital densities are summed according to the bias energy."// &
"For negative values, states in the range ]EF+bias,EF] are summed,"// &
"While positive values sum states in the range [EF,EF+bias[."// &
"If postive biases are used, sufficiently many unoccupied stated"// &
" (see ADDED_MOS and NLUMO ) should be computed.", &
n_var=-1, type_of_var=real_t, default_r_vals=(/0.0_dp/), unit_str='eV')
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="TH_TORB", &
description="Tip orbital symmetry in Tersoff-Hamann approximation to compute STM images", &
repeats=.TRUE., &
default_i_val=orb_s, &
usage="TH_TORB s dz2", &
enum_c_vals=s2a("S", "PX", "PY", "PZ", "DXY", "DYZ", "DZX", "DX2", "DY2", "DZ2"), &
enum_i_vals=(/orb_s, orb_px, orb_py, orb_pz, orb_dxy, orb_dyz, orb_dzx, orb_dx2, orb_dy2, orb_dz2/), &
enum_desc=s2a("s orbital", "px orbital", "py orbital", "pz orbital", &
"dxy orbital", "dyz orbital", "dzx orbital", "x^2 orbital", "y^2 orbital", "z^2 orbital"))
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="REF_ENERGY", &
description="By default the reference energy is the Fermi energy. In order to compare"// &
" with STS experiments, where specific energy ranges are addressed, here"// &
" one can set a different reference energy."// &
" The energy range is anyway controlled by the BIAS", &
type_of_var=real_t, default_r_val=0.0_dp, unit_str='eV')
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="APPEND", &
description="append the cube files when they already exist", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
END SUBROUTINE create_stm_section
! **************************************************************************************************
!> \brief ...
!> \param section ...
! **************************************************************************************************
SUBROUTINE create_wfn_mix_section(section)
TYPE(section_type), POINTER :: section
CHARACTER(len=*), PARAMETER :: routineN = 'create_wfn_mix_section', &
routineP = moduleN//':'//routineN
TYPE(keyword_type), POINTER :: keyword
TYPE(section_type), POINTER :: subsection
NULLIFY (subsection)
NULLIFY (keyword)
CPASSERT(.NOT. ASSOCIATED(section))
CALL section_create(section, __LOCATION__, name="WFN_MIX", &
description="A section that allows manipulation of the MO coeffs,"// &
" e.g. for changing a ground state into an excited state."// &
"Starting from a copy of the original MOs, changes can be made"// &
"by adding linear combinations of HOMO/LUMO of the original MOs to the result MOs", &
n_keywords=1, n_subsections=0, repeats=.FALSE.)
CALL keyword_create(keyword, __LOCATION__, name="OVERWRITE_MOS", &
description="If the keyword is active molecular orbitals in memory will be replaced by the mixed wfn."// &
" In combination with RTP or EMD no restart will be required to use the mixed wfn.", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL section_create(subsection, __LOCATION__, name="UPDATE", &
description="update a result MO with with a linear combination of of original MOs."// &
" This section can be repeated to build arbitrary linear combinations using repeatedly y=a*y+b*x.", &
n_keywords=1, n_subsections=0, repeats=.TRUE.)
CALL keyword_create(keyword, __LOCATION__, name="RESULT_MO_INDEX", &
description="Index of the MO (y) to be modified. Counting down in energy with HOMO=1", &
usage="RESULT_MO_INDEX 1", type_of_var=integer_t, default_i_val=0)
CALL section_add_keyword(subsection, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="RESULT_MARKED_STATE", &
description="Specifies the MO according to "// &
"the marks set in MOLECULAR_STATES. The value corresponds to the repetition "// &
" of MARK_STATES in MOLECULAR_STATES", &
usage="ORIG_MARKED_STATE 1", type_of_var=integer_t, default_i_val=0)
CALL section_add_keyword(subsection, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="RESULT_SPIN_INDEX", &
description="Spin of the MO (y) to be modified.", &
enum_c_vals=s2a("Alpha", "Beta"), &
enum_i_vals=(/1, 2/), & ! direct index in array
default_i_val=1, &
enum_desc=s2a("Majority spin", "Minority spin"))
CALL section_add_keyword(subsection, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="RESULT_SCALE", &
description="Scaling factor of the result variable (a).", &
usage="RESULT_SCALE 0.0", type_of_var=real_t)
CALL section_add_keyword(subsection, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="ORIG_MO_INDEX", &
description="Index of the original MO (x). "// &
"Counting down in energy with HOMO=1 or up from LUMO=1, depending on ORIG_IS_VIRTUAL.", &
usage="ORIG_MO_INDEX 1", type_of_var=integer_t, default_i_val=0)
CALL section_add_keyword(subsection, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="ORIG_MARKED_STATE", &
description="Specifies the MO according to "// &
"the marks set in MOLECULAR_STATES. The value corresponds to the repetition "// &
" of MARK_STATES in MOLECULAR_STATES", &
usage="ORIG_MARKED_STATE 1", type_of_var=integer_t, default_i_val=0)
CALL section_add_keyword(subsection, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="ORIG_SPIN_INDEX", &
description="Spin of the MO (x) to be modified.", &
enum_c_vals=s2a("Alpha", "Beta"), &
enum_i_vals=(/1, 2/), & ! direct index in array
default_i_val=1, &
enum_desc=s2a("Majority spin", "Minority spin"))
CALL section_add_keyword(subsection, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="ORIG_SCALE", &
description="Scaling factor of the result variable (b).", &
usage="ORIG_SCALE 0.0", type_of_var=real_t)
CALL section_add_keyword(subsection, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="ORIG_IS_VIRTUAL", &
description="The original MO (x) is a LUMO.", &
usage="ORIG_IS_VIRTUAL", type_of_var=logical_t)
CALL section_add_keyword(subsection, keyword)
CALL keyword_release(keyword)
CALL section_add_subsection(section, subsection)
CALL section_release(subsection)
END SUBROUTINE create_wfn_mix_section
! **************************************************************************************************
!> \brief creates the input section for the molecular states
!> \param print_key ...
!> \author teo
! **************************************************************************************************
SUBROUTINE create_molecular_states_section(print_key)
TYPE(section_type), POINTER :: print_key
CHARACTER(len=*), PARAMETER :: routineN = 'create_molecular_states_section', &
routineP = moduleN//':'//routineN
TYPE(keyword_type), POINTER :: keyword
TYPE(section_type), POINTER :: print_key2
CPASSERT(.NOT. ASSOCIATED(print_key))
NULLIFY (print_key2, keyword)
CALL cp_print_key_section_create(print_key, __LOCATION__, "MOLECULAR_STATES", &
description="Controls printing of molecular states ", &
print_level=high_print_level, filename=" ", citations=(/Hunt2003/))
CALL keyword_create( &
keyword, __LOCATION__, name="CUBE_EVAL_RANGE", &
description="only write cubes if the eigenvalues of the corresponding molecular states lie in the given interval. "// &
"Default is all states.", &
usage="CUBE_EVAL_RANGE -1.0 1.0", unit_str="hartree", n_var=2, type_of_var=real_t)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="MARK_STATES", &
description="Can be used to mark given molecular states."// &
" Sets a mark to both, occupied and unoccupied states. "// &
"Occupied states are counted beginning with HOMO=1, "// &
"unoccupied states are counted beginning with LUMO=1, "// &
"This is only meaningful in combination with WFN_MIX. "// &
"First integer specifies the molecule, second integer specifies the state.", &
usage="MARK_STATES integer integer", &
n_var=2, default_i_vals=(/-1, -1/), type_of_var=integer_t, repeats=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL cp_print_key_section_create(print_key2, __LOCATION__, "cubes", &
description="Controls the printing of cube files", &
print_level=high_print_level, filename="")
CALL keyword_create(keyword, __LOCATION__, name="stride", &
description="The stride (X,Y,Z) used to write the cube file "// &
"(larger values result in smaller cube files). You can provide 3 numbers (for X,Y,Z) or"// &
" 1 number valid for all components.", &
usage="STRIDE 2 2 2", n_var=-1, default_i_vals=(/2, 2, 2/), type_of_var=integer_t)
CALL section_add_keyword(print_key2, keyword)
CALL keyword_release(keyword)
CALL section_add_subsection(print_key, print_key2)
CALL section_release(print_key2)
END SUBROUTINE create_molecular_states_section
! **************************************************************************************************
!> \brief ...
!> \param print_key ...
! **************************************************************************************************
SUBROUTINE create_wannier_states_section(print_key)
TYPE(section_type), POINTER :: print_key
CHARACTER(len=*), PARAMETER :: routineN = 'create_wannier_states_section', &
routineP = moduleN//':'//routineN
TYPE(keyword_type), POINTER :: keyword
TYPE(section_type), POINTER :: print_key2
CPASSERT(.NOT. ASSOCIATED(print_key))
NULLIFY (print_key2, keyword)
CALL cp_print_key_section_create(print_key, __LOCATION__, "WANNIER_STATES", &
description="Controls printing of molecular states ", &
print_level=high_print_level, filename=" ")
CALL keyword_create( &
keyword, __LOCATION__, name="CUBE_EVAL_RANGE", &
description="only write cubes if the eigenvalues of the corresponding molecular states lie in the given interval. "// &
"Default is all states.", &
usage="CUBE_EVAL_RANGE -1.0 1.0", unit_str="hartree", n_var=2, type_of_var=real_t)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="MARK_STATES", &
description="Can be used to mark given molecular states."// &
" Sets a mark to both, occupied and unoccupied states. "// &
"Occupied states are counted beginning with HOMO=1, "// &
"unoccupied states are counted beginning with LUMO=1, "// &
"This is only meaningful in combination with WFN_MIX. "// &
"First integer specifies the molecule, second integer specifies the state.", &
usage="MARK_STATES integer integer", &
n_var=2, default_i_vals=(/-1, -1/), type_of_var=integer_t, repeats=.TRUE.)
CALL section_add_keyword(print_key, keyword)
CALL keyword_release(keyword)
CALL cp_print_key_section_create(print_key2, __LOCATION__, "cubes", &
description="Controls the printing of cube files", &
print_level=high_print_level, filename="")
CALL keyword_create(keyword, __LOCATION__, name="stride", &
description="The stride (X,Y,Z) used to write the cube file "// &
"(larger values result in smaller cube files). You can provide 3 numbers (for X,Y,Z) or"// &
" 1 number valid for all components.", &
usage="STRIDE 2 2 2", n_var=-1, default_i_vals=(/2, 2, 2/), type_of_var=integer_t)
CALL section_add_keyword(print_key2, keyword)
CALL keyword_release(keyword)
CALL section_add_subsection(print_key, print_key2)
CALL section_release(print_key2)
END SUBROUTINE create_wannier_states_section
! **************************************************************************************************
!> \brief creates the input section for the qs part
!> \param section the section to create
!> \author teo
! **************************************************************************************************
SUBROUTINE create_qs_section(section)
TYPE(section_type), POINTER :: section
CHARACTER(len=*), PARAMETER :: routineN = 'create_qs_section', &
routineP = moduleN//':'//routineN
TYPE(keyword_type), POINTER :: keyword
TYPE(section_type), POINTER :: subsection
CPASSERT(.NOT. ASSOCIATED(section))
CALL section_create(section, __LOCATION__, name="qs", &
description="parameters needed to set up the Quickstep framework", &
n_keywords=1, n_subsections=0, repeats=.FALSE.)
NULLIFY (keyword, subsection)
! Reals
CALL keyword_create(keyword, __LOCATION__, name="EPS_DEFAULT", &
description="Try setting all EPS_xxx to values leading to an energy correct up to EPS_DEFAULT", &
usage="EPS_DEFAULT real", default_r_val=1.0E-10_dp)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="EPS_CORE_CHARGE", &
description="Precision for mapping the core charges.Overrides EPS_DEFAULT/100.0 value", &
usage="EPS_CORE_CHARGE real", type_of_var=real_t)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create( &
keyword, __LOCATION__, name="EPS_GVG_RSPACE", &
variants=(/"EPS_GVG"/), &
description="Sets precision of the realspace KS matrix element integration. Overrides SQRT(EPS_DEFAULT) value", &
usage="EPS_GVG_RSPACE real", type_of_var=real_t)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="EPS_PGF_ORB", &
description="Sets precision of the overlap matrix elements. Overrides SQRT(EPS_DEFAULT) value", &
usage="EPS_PGF_ORB real", type_of_var=real_t)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create( &
keyword, __LOCATION__, name="EPS_KG_ORB", &
description="Sets precision used in coloring the subsets for the Kim-Gordon method. Overrides SQRT(EPS_DEFAULT) value", &
usage="EPS_KG_ORB 1.0E-8", &
type_of_var=real_t)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="EPS_PPL", &
description="Adjusts the precision for the local part of the pseudo potential. ", &
usage="EPS_PPL real", type_of_var=real_t, default_r_val=1.0E-2_dp)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create( &
keyword, __LOCATION__, name="EPS_PPNL", &
description="Sets precision of the non-local part of the pseudo potential. Overrides sqrt(EPS_DEFAULT) value", &
usage="EPS_PPNL real", type_of_var=real_t)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="EPS_CPC", &
description="Sets precision of the GAPW projection. Overrides EPS_DEFAULT value", &
usage="EPS_CPC real", type_of_var=real_t)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="EPS_RHO", &
description="Sets precision of the density mapping on the grids.Overrides EPS_DEFAULT value", &
usage="EPS_RHO real", type_of_var=real_t)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="EPS_RHO_RSPACE", &
description="Sets precision of the density mapping in rspace.Overrides EPS_DEFAULT value."// &
".Overrides EPS_RHO value", &
usage="EPS_RHO_RSPACE real", type_of_var=real_t)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="EPS_RHO_GSPACE", &
description="Sets precision of the density mapping in gspace.Overrides EPS_DEFAULT value."// &
".Overrides EPS_RHO value", &
usage="EPS_RHO_GSPACE real", type_of_var=real_t)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="EPS_FILTER_MATRIX", &
description="Sets the threshold for filtering matrix elements.", &
usage="EPS_FILTER_MATRIX 1.0E-6", type_of_var=real_t, default_r_val=0.0E0_dp)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="EPSFIT", &
variants=(/"EPS_FIT"/), &
description="GAPW: precision to give the extention of a hard gaussian ", &
usage="EPSFIT real", default_r_val=1.0E-4_dp)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="EPSISO", &
variants=(/"EPS_ISO"/), &
description="GAPW: precision to determine an isolated projector", &
usage="EPSISO real", default_r_val=1.0E-12_dp)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="EPSSVD", &
variants=(/"EPS_SVD"/), &
description="GAPW: tolerance used in the singular value decomposition of the projector matrix", &
usage="EPS_SVD real", default_r_val=1.0E-8_dp)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="EPSRHO0", &
variants=s2a("EPSVRHO0", "EPS_VRHO0"), &
description="GAPW : precision to determine the range of V(rho0-rho0soft)", &
usage="EPSRHO0 real", default_r_val=1.0E-6_dp)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="ALPHA0_HARD", &
variants=s2a("ALPHA0_H", "ALPHA0"), &
description="GAPW: Exponent for hard compensation charge", &
usage="ALPHA0_HARD real", default_r_val=0.0_dp)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create( &
keyword, __LOCATION__, name="FORCE_PAW", &
description="Use the GAPW scheme also for atoms with soft basis sets, i.e. "// &
" the local densities are computed even if hard and soft should be equal. "// &
"If this keyword is not set to true, those atoms with soft basis sets are treated by a GPW scheme, i.e. "// &
"the corresponding density contribution goes on the global grid and is expanded in PW. "// &
" This option nullifies the effect of the GPW_TYPE in the atomic KIND", &
usage="FORCE_PAW", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="MAX_RAD_LOCAL", &
description="GAPW : maximum radius of gaussian functions"// &
" included in the generation of projectors", &
usage="MAX_RAD_LOCAL real", default_r_val=25.0_dp)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
! Logicals
CALL keyword_create(keyword, __LOCATION__, name="LS_SCF", &
description="Perform a linear scaling SCF", &
usage="LS_SCF", lone_keyword_l_val=.TRUE., &
default_l_val=.FALSE.)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="ALMO_SCF", &
description="Perform ALMO SCF", &
usage="ALMO_SCF", lone_keyword_l_val=.TRUE., &
default_l_val=.FALSE.)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="TRANSPORT", &
description="Perform transport calculations (coupling CP2K and OMEN)", &
usage="TRANSPORT", lone_keyword_l_val=.TRUE., &
default_l_val=.FALSE.)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="KG_METHOD", &
description="Use a Kim-Gordon-like scheme.", &
usage="KG_METHOD", lone_keyword_l_val=.TRUE., &
default_l_val=.FALSE., citations=(/Iannuzzi2006, Brelaz1979, Andermatt2016/))
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="MAP_CONSISTENT", &
description="Compute the exact derivative (Hks) of the energy with respect to the density matrix. "// &
"This is slightly more expensive than using an approximate computation, "// &
"but consistent mapping can improve the stability of the SCF procedure, "// &
"especially for a tight EPS_SCF and a less tight EPS_DEFAULT.", &
usage="MAP_CONSISTENT FALSE", &
default_l_val=.TRUE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="REF_EMBED_SUBSYS", &
description="A total, reference, system in DFT embedding. ", &
usage="REF_EMBED_SUBSYS FALSE", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="CLUSTER_EMBED_SUBSYS", &
description="A cluster treated with DFT in DFT embedding. ", &
usage="CLUSTER_EMBED_SUBSYS FALSE", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="HIGH_LEVEL_EMBED_SUBSYS", &
description="A cluster treated with a high-level method in DFT embedding. ", &
usage="HIGH_LEVEL_EMBED_SUBSYS FALSE", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="DFET_EMBEDDED", &
description="Calculation with DFT-embedding potential. ", &
usage="DFET_EMBEDDED FALSE", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="DMFET_EMBEDDED", &
description="Calculation with DM embedding potential. ", &
usage="DMFET_EMBEDDED FALSE", &
default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
! Integers
CALL keyword_create(keyword, __LOCATION__, name="STO_NG", &
description="Order of Gaussian type expansion of Slater orbital basis sets.", &
usage="STO_NG", default_i_val=6)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="LMAXN1", &
variants=(/"LMAXRHO1"/), &
description="GAPW : max L number for expansion of the atomic densities in spherical gaussians", &
usage="LMAXN1 integer", &
default_i_val=-1)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="LMAXN0", &
variants=(/"LMAXRHO0"/), &
description="GAPW : max L number for the expansion compensation densities in spherical gaussians", &
usage="LMAXN0 integer", &
default_i_val=2)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="LADDN0", &
description="GAPW : integer added to the max L of the basis set, used to determine the "// &
"maximum value of L for the compensation charge density.", &
usage="LADDN0 integer", &
default_i_val=99)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
! Characters
CALL keyword_create(keyword, __LOCATION__, name="QUADRATURE", &
description="GAPW: algorithm to construct the atomic radial grids", &
usage="QUADRATURE GC_SIMPLE", &
enum_c_vals=s2a("GC_SIMPLE", "GC_TRANSFORMED", "GC_LOG"), &
enum_i_vals=(/do_gapw_gcs, do_gapw_gct, do_gapw_log/), &
enum_desc=s2a("Gauss-Chebyshev quadrature", &
"Transformed Gauss-Chebyshev quadrature", &
"Logarithmic transformed Gauss-Chebyshev quadrature"), &
default_i_val=do_gapw_log)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="PW_GRID", &
description="What kind of PW_GRID should be employed", &
usage="PW_GRID NS-FULLSPACE", &
enum_c_vals=s2a("SPHERICAL", "NS-FULLSPACE", "NS-HALFSPACE"), &
enum_desc=s2a("- not tested", " tested", " - not tested"), &
enum_i_vals=(/do_pwgrid_spherical, do_pwgrid_ns_fullspace, do_pwgrid_ns_halfspace/), &
default_i_val=do_pwgrid_ns_fullspace)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="PW_GRID_LAYOUT", &
description="Force a particular real-space layout for the plane waves grids. "// &
"Numbers &le; 0 mean that this dimension is free, incorrect layouts will be ignored. "// &
"The default (/-1,-1/) causes CP2K to select a good value, "// &
"i.e. plane distributed for large grids, more general distribution for small grids.", &
usage="PW_GRID_LAYOUT 4 16", &
repeats=.FALSE., n_var=2, &
default_i_vals=(/-1, -1/))
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="PW_GRID_BLOCKED", &
description="Can be used to set the distribution in g-space for the pw grids and their FFT.", &
usage="PW_GRID_BLOCKED FREE", &
enum_c_vals=s2a("FREE", "TRUE", "FALSE"), &
enum_desc=s2a("CP2K will select an appropriate value", "blocked", "not blocked"), &
enum_i_vals=(/do_pw_grid_blocked_free, do_pw_grid_blocked_true, do_pw_grid_blocked_false/), &
default_i_val=do_pw_grid_blocked_free)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create( &
keyword, __LOCATION__, name="EXTRAPOLATION", &
variants=s2a("INTERPOLATION", "WF_INTERPOLATION"), &
description="Extrapolation strategy for the wavefunction during e.g. MD. "// &
"Not all options are available for all simulation methods. "// &
"PS and ASPC are recommended, see also EXTRAPOLATION_ORDER.", &
citations=(/Kolafa2004, VandeVondele2005a/), &
usage="EXTRAPOLATION PS", &
enum_c_vals=s2a("USE_GUESS", "USE_PREV_P", "USE_PREV_RHO_R", "LINEAR_WF", &
"LINEAR_P", "LINEAR_PS", "USE_PREV_WF", "PS", "FROZEN", "ASPC"), &
enum_desc=s2a( &
"Use the method specified with SCF_GUESS, i.e. no extrapolation", &
"Use the previous density matrix", &
"Use the previous density in real space", &
"Linear extrapolation of the wavefunction (not available for K-points)", &
"Linear extrapolation of the density matrix", &
"Linear extrapolation of the density matrix times the overlap matrix (not available for K-points)", &
"Use the previous wavefunction (not available for K-points)", &
"Higher order extrapolation of the density matrix times the overlap matrix (not available for K-points)", &
"Frozen ...", &
"Always stable predictor corrector, similar to PS, but going for MD stability instead of initial guess accuracy. "// &
"(not available for K-points)"), &
enum_i_vals=(/ &
wfi_use_guess_method_nr, &
wfi_use_prev_p_method_nr, &
wfi_use_prev_rho_r_method_nr, &
wfi_linear_wf_method_nr, &
wfi_linear_p_method_nr, &
wfi_linear_ps_method_nr, &
wfi_use_prev_wf_method_nr, &
wfi_ps_method_nr, &
wfi_frozen_method_nr, &
wfi_aspc_nr/), &
default_i_val=wfi_aspc_nr)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="EXTRAPOLATION_ORDER", &
description="Order for the PS or ASPC extrapolation (typically 2-4). "// &
"Higher order might bring more accuracy, but comes, "// &
"for large systems, also at some cost. "// &
"In some cases, a high order extrapolation is not stable,"// &
" and the order needs to be reduced.", &
usage="EXTRAPOLATION_ORDER {integer}", default_i_val=3)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="METHOD", &
description="Specifies the electronic structure method that should be employed", &
usage="METHOD GAPW", &
enum_c_vals=s2a("GAPW", "GAPW_XC", "GPW", "LRIGPW", "RIGPW", &
"MNDO", "MNDOD", "AM1", "PM3", "PM6", "PM6-FM", "PDG", "RM1", "PNNL", "DFTB", "xTB", "OFGPW"), &
enum_desc=s2a("Gaussian and augmented plane waves method", &
"Gaussian and augmented plane waves method only for XC", &
"Gaussian and plane waves method", &
"Local resolution of identity method", &
"Resolution of identity method for HXC terms", &
"MNDO semiempirical", "MNDO-d semiempirical", "AM1 semiempirical", &
"PM3 semiempirical", "PM6 semiempirical", "PM6-FM semiempirical", "PDG semiempirical", &
"RM1 semiempirical", &
"PNNL semiempirical", &
"DFTB Density Functional based Tight-Binding", &
"GFN-xTB Extended Tight-Binding", &
"OFGPW Orbital-free GPW method"), &
enum_i_vals=(/do_method_gapw, do_method_gapw_xc, do_method_gpw, do_method_lrigpw, do_method_rigpw, &
do_method_mndo, do_method_mndod, do_method_am1, do_method_pm3, &
do_method_pm6, do_method_pm6fm, do_method_pdg, do_method_rm1, &
do_method_pnnl, do_method_dftb, do_method_xtb, do_method_ofgpw/), &
citations=(/Lippert1997, Lippert1999, Krack2000, VandeVondele2005a, &
VandeVondele2006, Dewar1977, Dewar1985, Rocha2006, Stewart1989, Thiel1992, &
Repasky2002, Stewart2007, VanVoorhis2015, Schenter2008/), &
default_i_val=do_method_gpw)
CALL section_add_keyword(section, keyword)
CALL keyword_release(keyword)
CALL keyword_create(keyword, __LOCATION__, name="CORE_PPL", &
description="Specifies the method used to calculate the local pseudopotential contribution.", &
usage="CORE_PPL ANALYTIC", &