79 numerical,
oe_gllb,
oe_lb,
oe_none,
oe_saop,
oe_sic,
orb_dx2,
orb_dxy,
orb_dy2,
orb_dyz, &
166#include "./base/base_uses.f90"
171 CHARACTER(len=*),
PARAMETER,
PRIVATE :: moduleN =
'input_cp2k_print_dft'
186 TYPE(
section_type),
POINTER :: print_key, sub_print_key, subsection
188 cpassert(.NOT.
ASSOCIATED(section))
190 description=
"Section of possible print options in DFT code.", &
191 n_keywords=0, n_subsections=1, repeats=.false.)
193 NULLIFY (print_key, keyword, subsection)
196 description=
"Controls the printing of the banner of the MM program", &
202 description=
"Controls the printing of a file with all basis sets used.", &
208 description=
"Controls the printing of the kinetic energy", &
214 description=
"Print all derivatives after the DFT calculation", &
218 description=
"Specify the number of digits used to print derivatives", &
227 description=
"Controls the printing of the neighbor lists", &
231 description=
"Activates the printing of the orbital "// &
232 "orbital neighbor lists, "// &
233 "i.e. the overlap neighbor lists", &
234 default_l_val=.false., &
235 lone_keyword_l_val=.true.)
240 name=
"sab_aux_fit", &
241 description=
"Activates the printing of the orbital "// &
242 "orbital neighbor lists wavefunction fitting basis, "// &
243 "i.e. the overlap neighbor lists", &
244 default_l_val=.false., &
245 lone_keyword_l_val=.true.)
250 name=
"sab_aux_fit_vs_orb", &
251 description=
"Activates the printing of the orbital "// &
252 "orbital mixed neighbor lists of wavefunction fitting basis, "// &
253 "and the orbital basis, i.e. the overlap neighbor lists", &
254 default_l_val=.false., &
255 lone_keyword_l_val=.true.)
261 description=
"Activates the printing of the vdW SCP "// &
263 default_l_val=.false., &
264 lone_keyword_l_val=.true.)
270 description=
"Activates the printing of the vdW "// &
271 "neighbor lists (from DFT, DFTB, SE), "// &
272 "i.e. the dispersion neighbor lists", &
273 default_l_val=.false., &
274 lone_keyword_l_val=.true.)
280 description=
"Activates the printing of the "// &
281 "neighbor lists used for coordination numbers in vdW DFT-D3", &
282 default_l_val=.false., &
283 lone_keyword_l_val=.true.)
289 description=
"Activates the printing of the orbital "// &
290 "nuclear attraction neighbor lists (erfc potential)", &
291 default_l_val=.false., &
292 lone_keyword_l_val=.true.)
298 description=
"Activates the printing of the orbital "// &
299 "GTH-PPL neighbor lists (local part of the "// &
300 "Goedecker-Teter-Hutter pseudo potentials)", &
301 default_l_val=.false., &
302 lone_keyword_l_val=.true.)
308 description=
"Activates the printing of the orbital "// &
309 "GTH-PPNL neighbor lists (non-local part of the "// &
310 "Goedecker-Teter-Hutter pseudo potentials)", &
311 default_l_val=.false., &
312 lone_keyword_l_val=.true.)
318 description=
"Activates the printing of the orbital "// &
319 "PAW-projector neighbor lists (only GAPW)", &
320 default_l_val=.false., &
321 lone_keyword_l_val=.true.)
327 description=
"Activates the printing of the two-center "// &
328 "neighbor lists for Coulomb type interactions in NDDO ", &
329 default_l_val=.false., &
330 lone_keyword_l_val=.true.)
336 description=
"Activates the printing of the long-range SE correction "// &
337 "neighbor lists (only when doing long-range SE with integral scheme KDSO and KDSO-d)", &
338 default_l_val=.false., lone_keyword_l_val=.true.)
344 description=
"Activates the printing of the DFTB Ewald "// &
346 default_l_val=.false., &
347 lone_keyword_l_val=.true.)
353 description=
"Activates the printing of the xTB sr-Coulomb "// &
355 default_l_val=.false., &
356 lone_keyword_l_val=.true.)
362 description=
"Activates the printing of core interaction "// &
364 default_l_val=.false., &
365 lone_keyword_l_val=.true.)
371 description=
"Activates the printing of XB interaction from (xTB) "// &
373 default_l_val=.false., &
374 lone_keyword_l_val=.true.)
379 name=
"sab_xtb_nonbond", &
380 description=
"Activates the printing of nonbonded interaction from (xTB) "// &
382 default_l_val=.false., &
383 lone_keyword_l_val=.true.)
389 description=
"Activates the printing of RI orbital-orbital "// &
391 default_l_val=.false., &
392 lone_keyword_l_val=.true.)
398 description=
"Activates the printing of RI basis-projector interaction "// &
400 default_l_val=.false., &
401 lone_keyword_l_val=.true.)
407 description=
"Activates the printing of the nuclear orbital "// &
408 "nuclear repulsion neighbor lists (erfc potential)", &
409 default_l_val=.false., &
410 lone_keyword_l_val=.true.)
418 description=
"Activates the printing of the subcells used for the "// &
419 "generation of neighbor lists.", unit_str=
"angstrom", &
425 description=
"Controls the printing of the ao (i.e. contracted gaussian) matrices (debug).", &
428 description=
"Print only the matrix data, not the row and column headers", &
429 default_l_val=.false., lone_keyword_l_val=.true.)
433 description=
"Specify the number of digits used to print the AO matrices", &
437 CALL keyword_create(keyword, __location__, name=
"CORE_HAMILTONIAN", &
438 description=
"If the printkey is activated controls the printing of the hamiltonian matrix", &
439 default_l_val=.false., lone_keyword_l_val=.true.)
443 description=
"If the printkey is activated controls the printing of the density (P) matrix", &
444 default_l_val=.false., lone_keyword_l_val=.true.)
448 CALL keyword_create(keyword, __location__, name=
"KINETIC_ENERGY", &
449 description=
"If the printkey is activated controls the printing of the kinetic energy matrix", &
450 default_l_val=.false., lone_keyword_l_val=.true.)
453 CALL keyword_create(keyword, __location__, name=
"KOHN_SHAM_MATRIX", &
454 description=
"If the printkey is activated controls the printing of the kohn-sham matrix", &
455 default_l_val=.false., lone_keyword_l_val=.true.)
459 keyword, __location__, name=
"MATRIX_VXC", &
460 description=
"If the printkey is activated compute and print the matrix of the exchange and correlation potential. "// &
461 "Only the GGA part for GPW is printed", &
462 default_l_val=.false., lone_keyword_l_val=.true.)
466 description=
"If the printkey is activated controls the printing of the orthogonalization matrix", &
467 default_l_val=.false., lone_keyword_l_val=.true.)
471 description=
"If the printkey is activated controls the printing of the overlap matrix", &
472 default_l_val=.false., lone_keyword_l_val=.true.)
475 CALL keyword_create(keyword, __location__, name=
"COMMUTATOR_HR", &
476 description=
"Controls the printing of the [H,r] commutator matrix", &
477 default_l_val=.false., lone_keyword_l_val=.true.)
481 CALL keyword_create(keyword, __location__, name=
"FERMI_CONTACT", &
482 description=
"If the printkey is activated controls the printing of the Fermi contact matrix", &
483 default_l_val=.false., lone_keyword_l_val=.true.)
487 keyword, __location__, name=
"PSO", &
488 description=
"If the printkey is activated controls the printing of the paramagnetic spin-orbit matrices", &
489 default_l_val=.false., lone_keyword_l_val=.true.)
493 keyword, __location__, name=
"EFG", &
494 description=
"If the printkey is activated controls the printing of the electric field gradient matrices", &
495 default_l_val=.false., lone_keyword_l_val=.true.)
498 CALL keyword_create(keyword, __location__, name=
"POTENTIAL_ENERGY", &
499 description=
"If the printkey is activated controls the printing of the potential energy matrix", &
500 default_l_val=.false., lone_keyword_l_val=.true.)
504 description=
"If the printkey is activated controls the printing of the OCE HARD matrix", &
505 default_l_val=.false., lone_keyword_l_val=.true.)
509 description=
"If the printkey is activated controls the printing of the OCE SOFT matrix", &
510 default_l_val=.false., lone_keyword_l_val=.true.)
514 description=
"If the printkey is activated controls the printing of the w matrix", &
515 default_l_val=.false., lone_keyword_l_val=.true.)
518 CALL keyword_create(keyword, __location__, name=
"W_MATRIX_AUX_FIT", &
519 description=
"If the printkey is activated controls the printing of the w matrix", &
520 default_l_val=.false., lone_keyword_l_val=.true.)
524 description=
"If the printkey is activated controls the printing "// &
525 "of derivatives (for the matrixes that support this)", &
526 default_l_val=.false., lone_keyword_l_val=.true.)
533 print_key, __location__,
"MO", &
534 description=
"Controls the printing of the molecular orbital (MO) information. The requested MO information "// &
535 "is printed for all occupied MOs by default. Use the MO_INDEX_RANGE keyword to restrict the number "// &
536 "of the MOs or to print the MO information for unoccupied MOs. With diagonalization, additional MOs "// &
537 "have to be made available for printout using the ADDED_MOS keyword in the SCF section. With OT, "// &
538 "it is sufficient to specify the desired MO_INDEX_RANGE. The OT eigensolver can be controlled with "// &
539 "the EPS_LUMO and MAX_ITER_LUMO keywords in the SCF section.", &
542 description=
"Print the MOs in the Cartesian basis instead of the default spherical basis.", &
543 default_l_val=.false., lone_keyword_l_val=.true.)
547 variants=
s2a(
"EIGENVALUES",
"EIGVALS"), &
548 description=
"Print the MO energies (eigenvalues).", &
549 default_l_val=.false., lone_keyword_l_val=.true.)
553 variants=
s2a(
"EIGENVECTORS",
"EIGVECS"), &
554 description=
"Print the MO coefficients (eigenvectors).", &
555 default_l_val=.false., lone_keyword_l_val=.true.)
558 CALL keyword_create(keyword, __location__, name=
"OCCUPATION_NUMBERS", &
559 variants=
s2a(
"OCCNUMS"), &
560 description=
"Print the MO occupation numbers.", &
561 default_l_val=.false., lone_keyword_l_val=.true.)
564 CALL keyword_create(keyword, __location__, name=
"OCCUPATION_NUMBERS_STATS", &
565 variants=
s2a(
"OCCNUMSTATS"), &
566 description=
"Print some stats (max number of occupied MOs, etc.) of the MO occupation numbers."// &
567 " First logical toggles stats printing, first real is the occupied threshold.", &
568 type_of_var=
char_t, n_var=-1, &
569 default_c_vals=[
".FALSE.",
"1.0E-6 "], &
570 lone_keyword_c_vals=[
".TRUE."], &
571 usage=
"OCCUPATION_NUMBERS_STATS {Logical} [{Real}]")
575 description=
"Specify the number of digits used to print the MO information.", &
580 name=
"MO_INDEX_RANGE", &
581 variants=
s2a(
"MO_RANGE",
"RANGE"), &
582 description=
"Print only the requested subset of MOs. The indices of the first and "// &
583 "the last MO have to be specified to define the range. -1 as the last MO index "// &
584 "prints all available orbitals with diagonalisation (ADDED_MOS) and all orbitals with OT.", &
588 default_i_vals=(/0, 0/), &
589 usage=
"MO_INDEX_RANGE 10 15")
596 description=
"Write the molecular orbitals in Molden file format, for visualisation.", &
599 description=
"Specifies the number of significant digits retained. 3 is OK for visualization.", &
600 usage=
"NDIGITS {int}", &
605 description=
"Representation of Gaussian-type orbitals", &
607 enum_c_vals=
s2a(
"CARTESIAN",
"SPHERICAL"), &
609 "Cartesian Gaussian orbitals. Use with caution", &
610 "Spherical Gaussian orbitals. Incompatible with VMD"), &
617 CALL create_mo_cubes_section(print_key)
621 CALL create_stm_section(print_key)
625 CALL create_wfn_mix_section(subsection)
630 description=
"Write a TREXIO file to disk.", &
631 n_keywords=1, n_subsections=0, repeats=.false.)
633 description=
"Body of Filename for the trexio file.", &
634 usage=
"FILENAME {name}", default_c_val=
"TREXIO", &
639 description=
"Store the MOs in the Cartesian basis instead of the default spherical basis.", &
640 default_l_val=.false., lone_keyword_l_val=.true.)
647 description=
"Controls the printing of some gapw related information (debug).", &
648 n_keywords=0, n_subsections=1, repeats=.false.)
650 description=
"If the printkey is activated controls if information on"// &
651 " the projectors is printed.", &
656 description=
"If the printkey is activated controls if information on rho0 is printed.", &
664 description=
"Controls the printing of dft control parameters.", &
670 description=
"Controls the printing of kpoint information.", &
676 CALL create_bandstructure_section(subsection)
681 description=
"Controls the checking and printing of an estimate "// &
682 "of the overlap matrix condition number", &
685 description=
"Calculate an estimate of the 1-norm condition number", &
686 default_l_val=.true., lone_keyword_l_val=.true.)
689 CALL keyword_create(keyword, __location__, name=
"DIAGONALIZATION", &
690 description=
"Calculate the 1- and 2-norm condition numbers using diagonalization", &
691 default_l_val=.false., lone_keyword_l_val=.true.)
695 description=
"Calculate the 2-norm condition number using the Arnoldi code (may not be reliable)", &
696 default_l_val=.false., lone_keyword_l_val=.true.)
703 description=
"Controls the printing of cube files with "// &
704 "the electronic density and, for LSD calculations, the spin density.", &
707 description=
"The stride (X,Y,Z) used to write the cube file "// &
708 "(larger values result in smaller cube files). You can provide 3 numbers (for X,Y,Z) or"// &
709 " 1 number valid for all components.", &
710 usage=
"STRIDE 2 2 2", n_var=-1, default_i_vals=(/2, 2, 2/), type_of_var=
integer_t)
714 CALL keyword_create(keyword, __location__, name=
"DENSITY_INCLUDE", &
715 description=
"Which parts of the density to include. In GAPW the electronic density "// &
716 "is divided into a hard and a soft component, and the default (TOTAL_HARD_APPROX) "// &
717 "is to approximate the hard density as a spherical gaussian and to print the smooth "// &
718 "density accurately. This avoids potential artefacts originating from the hard density. "// &
719 "If the TOTAL_DENSITY keyword is used the hard density will be computed more accurately "// &
720 "but may introduce non-physical features. The SOFT_DENSITY keyword will lead to only the "// &
721 "soft density being printed. In GPW these options have no effect and the cube file will "// &
722 "only contain the valence electron density.", &
723 usage=
"DENSITY_INCLUDE TOTAL_HARD_APPROX", &
724 enum_c_vals=
s2a(
"TOTAL_HARD_APPROX",
"TOTAL_DENSITY",
"SOFT_DENSITY"), &
725 enum_desc=
s2a(
"Print (hard+soft) density where the hard components shape is approximated", &
726 "Print (hard+soft) density. Only has an effect "// &
727 "if PAW atoms are present. NOTE: The total "// &
728 "in real space might exhibit unphysical features "// &
729 "like spikes due to the finite and thus "// &
730 "truncated g vector", &
731 "Print only the soft density"), &
740 description=
"append the cube files when they already exist", &
741 default_l_val=.false., lone_keyword_l_val=.true.)
745 CALL keyword_create(keyword, __location__, name=
"XRD_INTERFACE", &
746 description=
"It activates the print out of exponents and coefficients for the"// &
747 " Gaussian expansion of the core densities, based on atom calculations for each kind."// &
748 " The resulting core dansities are needed to compute the form factors."// &
749 " If GAPW the local densities are also given in terms of a Gaussian expansion,"// &
750 " by fitting the difference between local-fhard and local-soft density for each atom."// &
751 " In this case the keyword SOFT_DENSITY is enabled.", &
752 default_l_val=.false., lone_keyword_l_val=.true.)
757 description=
"Number of Gaussian functions used in the expansion of atomic (core) density", &
758 usage=
"NGAUSS 10", n_var=1, default_i_val=12, type_of_var=
integer_t)
766 description=
"Controls printing of cube files with "// &
767 "the total density (electrons+atomic core). Note that "// &
768 "the value of the total density is positive where the "// &
769 "electron density dominates and negative where the core is. "// &
770 "When GPW is enabled this will simply print the combined density "// &
771 "of the valence electrons and charge-balanced core. In GAPW the "// &
772 "electronic density (hard+soft plus a correction term) is printed "// &
773 "together with the charge-balanced core density to produce a complete "// &
774 "representation of the total density.", &
777 description=
"The stride (X,Y,Z) used to write the cube file "// &
778 "(larger values result in smaller cube files). You can provide 3 numbers (for X,Y,Z) or"// &
779 " 1 number valid for all components.", &
780 usage=
"STRIDE 2 2 2", n_var=-1, default_i_vals=(/2, 2, 2/), type_of_var=
integer_t)
785 description=
"append the cube files when they already exist", &
786 default_l_val=.false., lone_keyword_l_val=.true.)
794 description=
"Controls the printing of a cube file with eletrostatic"// &
795 " potential generated by the total density (electrons+ions). It is"// &
796 " valid only for QS with GPW formalism."// &
797 " Note that by convention the potential has opposite sign than the expected physical one.", &
800 description=
"The stride (X,Y,Z) used to write the cube file "// &
801 "(larger values result in smaller cube files). You can provide 3 numbers (for X,Y,Z) or"// &
802 " 1 number valid for all components.", &
803 usage=
"STRIDE 2 2 2", n_var=-1, default_i_vals=(/2, 2, 2/), type_of_var=
integer_t)
807 description=
"append the cube files when they already exist", &
808 default_l_val=.false., lone_keyword_l_val=.true.)
816 description=
"Controls the printing of a cube file with external"// &
817 " potential from the DFT%EXTERNAL_POTENTIAL section only.", &
820 description=
"The stride (X,Y,Z) used to write the cube file "// &
821 "(larger values result in smaller cube files). You can provide 3 numbers (for X,Y,Z) or"// &
822 " 1 number valid for all components.", &
823 usage=
"STRIDE 2 2 2", n_var=-1, default_i_vals=(/2, 2, 2/), type_of_var=
integer_t)
827 description=
"append the cube files when they already exist", &
828 default_l_val=.false., lone_keyword_l_val=.true.)
837 description=
"Controls the output of the electron density to the losslessly"// &
838 " compressed BQB file format, see [Brehm2018]"// &
839 " (via LibBQB see <https://brehm-research.de/bqb>)."// &
840 " Currently does not work with changing cell vector (NpT ensemble).", &
845 description=
"Skips the first step of a MD run (avoids duplicate step if restarted).", &
846 usage=
"SKIP_FIRST T", default_l_val=.false., lone_keyword_l_val=.true.)
850 CALL keyword_create(keyword, __location__, name=
"STORE_STEP_NUMBER", &
851 description=
"Stores the step number and simulation time in the comment line of each BQB"// &
852 " frame. Switch it off for binary compatibility with original CP2k CUBE files.", &
853 usage=
"STORE_STEP_NUMBER F", default_l_val=.true., lone_keyword_l_val=.true.)
858 description=
"Performs an on-the-fly decompression of each compressed BQB frame to check"// &
859 " whether the volumetric data exactly matches, and aborts the run if not so.", &
860 usage=
"CHECK T", default_l_val=.false., lone_keyword_l_val=.true.)
865 description=
"Specify this keyword to overwrite the output BQB file if"// &
866 " it already exists. By default, the data is appended to an existing file.", &
867 usage=
"OVERWRITE T", default_l_val=.false., lone_keyword_l_val=.true.)
872 description=
"Controls how many previous steps are taken into account for extrapolation in"// &
873 " compression. Use a value of 1 to compress the frames independently.", &
874 usage=
"HISTORY 10", n_var=1, default_i_val=10, type_of_var=
integer_t)
878 CALL keyword_create(keyword, __location__, name=
"PARAMETER_KEY", &
879 description=
"Allows to supply previously optimized compression parameters via a"// &
880 " parameter key (alphanumeric character sequence starting with 'at')."// &
881 " Just leave away the 'at' sign here, because CP2k will otherwise"// &
882 " assume it is a variable name in the input", &
883 usage=
"PARAMETER_KEY <KEY>", n_var=1, default_c_val=
"", type_of_var=
char_t)
888 description=
"Controls the time spent to optimize the parameters for compression efficiency.", &
889 usage=
"OPTIMIZE {OFF,QUICK,NORMAL,PATIENT,EXHAUSTIVE}", repeats=.false., n_var=1, &
891 enum_c_vals=
s2a(
"OFF",
"QUICK",
"NORMAL",
"PATIENT",
"EXHAUSTIVE"), &
892 enum_desc=
s2a(
"No optimization (use defaults)",
"Quick optimization", &
893 "Standard optimization",
"Precise optimization",
"Exhaustive optimization"), &
907 CALL create_implicit_psolver_section(subsection)
913 description=
"Controls the printing of a cube file with xc"// &
914 " potential generated by the ZMP method (for the moment). It is"// &
915 " valid only for QS with GPW formalism .", &
918 description=
"The stride (X,Y,Z) used to write the cube file "// &
919 "(larger values result in smaller cube files). You can provide 3 numbers (for X,Y,Z) or"// &
920 " 1 number valid for all components.", &
921 usage=
"STRIDE 2 2 2", n_var=-1, default_i_vals=(/2, 2, 2/), type_of_var=
integer_t)
925 description=
"append the cube files when they already exist", &
926 default_l_val=.false., lone_keyword_l_val=.true.)
934 description=
"Controls the printing of cube files with electric"// &
935 " field generated by the total density (electrons+ions). It is"// &
936 " valid only for QS with GPW formalism.", &
939 description=
"The stride (X,Y,Z) used to write the cube file "// &
940 "(larger values result in smaller cube files). You can provide 3 numbers (for X,Y,Z) or"// &
941 " 1 number valid for all components.", &
942 usage=
"STRIDE 2 2 2", n_var=-1, default_i_vals=(/2, 2, 2/), type_of_var=
integer_t)
946 description=
"append the cube files when they already exist", &
947 default_l_val=.false., lone_keyword_l_val=.true.)
955 description=
"Controls printing of cube files with"// &
956 " the electron localization function (ELF). Note that"// &
957 " the value of ELF is defined between 0 and 1: Pauli kinetic energy density normalized"// &
958 " by the kinetic energy density of a uniform el. gas of same density.", &
961 description=
"The stride (X,Y,Z) used to write the cube file "// &
962 "(larger values result in smaller cube files). You can provide 3 numbers (for X,Y,Z) or"// &
963 " 1 number valid for all components.", &
964 usage=
"STRIDE 2 2 2", n_var=-1, default_i_vals=(/2, 2, 2/), type_of_var=
integer_t)
968 description=
"append the cube files when they already exist", &
969 default_l_val=.false., lone_keyword_l_val=.true.)
973 CALL keyword_create(keyword, __location__, name=
"density_cutoff", &
975 usage=
"density_cutoff 0.0001", &
979 default_r_val=1.0e-10_dp)
987 description=
"Controls the printing of cube files with the local"// &
988 " energy. It is valid only for QS with GPW/GAPW formalism."// &
989 " Meta and hybrid functionals are not possible.", &
992 description=
"The stride (X,Y,Z) used to write the cube file "// &
993 "(larger values result in smaller cube files). You can provide 3 numbers (for X,Y,Z) or"// &
994 " 1 number valid for all components.", &
995 usage=
"STRIDE 2 2 2", n_var=-1, default_i_vals=(/2, 2, 2/), type_of_var=
integer_t)
999 description=
"append the cube files when they already exist", &
1000 default_l_val=.false., lone_keyword_l_val=.true.)
1008 description=
"Controls the printing of cube files with the local"// &
1009 " stress. It is valid only for QS with GPW/GAPW formalism."// &
1010 " Meta and hybrid functionals are not possible.", &
1013 description=
"The stride (X,Y,Z) used to write the cube file "// &
1014 "(larger values result in smaller cube files). You can provide 3 numbers (for X,Y,Z) or"// &
1015 " 1 number valid for all components.", &
1016 usage=
"STRIDE 2 2 2", n_var=-1, default_i_vals=(/2, 2, 2/), type_of_var=
integer_t)
1020 description=
"append the cube files when they already exist", &
1021 default_l_val=.false., lone_keyword_l_val=.true.)
1028 CALL create_dos_section(print_key)
1036 CALL create_wannier_section(print_key)
1043 keyword, __location__, &
1044 name=
"MAX_MOMENT", &
1045 description=
"Maximum moment to be calculated. Values higher than 1 not implemented under periodic boundaries.", &
1046 usage=
"MAX_MOMENT {integer}", &
1055 description=
"Calculate also magnetic moments, only implemented without periodic boundaries", &
1056 usage=
"MAGNETIC yes", &
1059 default_l_val=.false., &
1060 lone_keyword_l_val=.true.)
1065 description=
"Calculate expectation values of the el. multipole moments in their velocity "// &
1066 "representation during RTP. Implemented up to el. quadrupole moment.", &
1067 usage=
"VEL_REPRS yes", &
1070 default_l_val=.false., &
1071 lone_keyword_l_val=.true.)
1076 description=
"Include non local commutator for velocity representations. "// &
1077 "Necessary for origin independent results.", &
1078 usage=
"COM_NL yes", &
1081 default_l_val=.false., &
1082 lone_keyword_l_val=.true.)
1086 name=
"SECOND_REFERENCE_POINT", &
1087 description=
"Use second reference point", &
1088 usage=
"SECOND_REFERENCE_POINT .TRUE.", &
1091 default_l_val=.false., &
1092 lone_keyword_l_val=.true.)
1096 variants=
s2a(
"REF_2"), &
1097 description=
"Define a second reference point for the calculation of the electrostatic moment.", &
1098 usage=
"REFERENCE_2 COM", &
1099 enum_c_vals=
s2a(
"COM",
"COAC",
"USER_DEFINED",
"ZERO"), &
1100 enum_desc=
s2a(
"Use Center of Mass", &
1101 "Use Center of Atomic Charges", &
1102 "Use User Defined Point (Keyword:REF_POINT)", &
1103 "Use Origin of Coordinate System"), &
1111 CALL keyword_create(keyword, __location__, name=
"REFERENCE_POINT_2", &
1112 variants=
s2a(
"REF_POINT_2"), &
1113 description=
"Fixed second reference point for the calculations of the electrostatic moment.", &
1114 usage=
"REFERENCE_POINT_2 x y z", &
1116 n_var=3, default_r_vals=(/0._dp, 0._dp, 0._dp/), &
1126 description=
"Controls the printing of the Mulliken (spin) population analysis", &
1128 common_iter_levels=1)
1130 keyword, __location__, &
1132 description=
"Print the gross orbital populations (GOP) in addition to the gross atomic populations (GAP) "// &
1133 "and net charges", &
1134 usage=
"PRINT_GOP yes", &
1137 default_l_val=.false., &
1138 lone_keyword_l_val=.true.)
1142 keyword, __location__, &
1144 description=
"Print all information including the full net AO and overlap population matrix", &
1145 usage=
"PRINT_ALL yes", &
1148 default_l_val=.false., &
1149 lone_keyword_l_val=.true.)
1157 description=
"Controls the printing of the Lowdin (spin) population analysis", &
1159 common_iter_levels=1)
1161 keyword, __location__, &
1163 description=
"Print the orbital populations in addition to the atomic populations and net charges", &
1164 usage=
"PRINT_GOP yes", &
1167 default_l_val=.false., &
1168 lone_keyword_l_val=.true.)
1172 keyword, __location__, &
1174 description=
"Print all information including the full symmetrically orthogonalised density matrix", &
1175 usage=
"PRINT_ALL yes", &
1178 default_l_val=.false., &
1179 lone_keyword_l_val=.true.)
1187 description=
"Controls the printing of the Hirshfeld (spin) population analysis", &
1189 common_iter_levels=1)
1190 CALL keyword_create(keyword, __location__, name=
"SELF_CONSISTENT", &
1191 description=
"Calculate charges from the Hirscheld-I (self_consistent) method."// &
1192 " This scales only the full shape function, not the added charge as in the original scheme.", &
1193 usage=
"SELF_CONSISTENT yes", repeats=.false., n_var=1, &
1194 default_l_val=.false., lone_keyword_l_val=.true.)
1197 CALL keyword_create(keyword, __location__, name=
"SHAPE_FUNCTION", &
1198 description=
"Type of shape function used for Hirshfeld partitioning.", &
1199 usage=
"SHAPE_FUNCTION {Gaussian,Density}", repeats=.false., n_var=1, &
1201 enum_c_vals=
s2a(
"GAUSSIAN",
"DENSITY"), &
1202 enum_desc=
s2a(
"Single Gaussian with Colvalent radius", &
1203 "Atomic density expanded in multiple Gaussians"), &
1207 CALL keyword_create(keyword, __location__, name=
"REFERENCE_CHARGE", &
1208 description=
"Charge of atomic partitioning function for Hirshfeld method.", &
1209 usage=
"REFERENCE_CHARGE {Atomic,Mulliken}", repeats=.false., n_var=1, &
1211 enum_c_vals=
s2a(
"ATOMIC",
"MULLIKEN"), &
1212 enum_desc=
s2a(
"Use atomic core charges",
"Calculate Mulliken charges"), &
1217 description=
"Use user defined radii to generate Gaussians."// &
1218 " These radii are defined by the keyword ATOMIC_RADII", &
1219 usage=
"USER_RADIUS yes", repeats=.false., n_var=1, &
1220 default_l_val=.false., lone_keyword_l_val=.true.)
1223 CALL keyword_create(keyword, __location__, name=
"ATOMIC_RADII", &
1224 description=
"Defines custom radii to setup the spherical Gaussians.", &
1225 usage=
"ATOMIC_RADII {real} {real} {real}", repeats=.false., &
1226 unit_str=
"angstrom", &
1227 type_of_var=
real_t, n_var=-1)
1235 description=
"Controls the printing of the EEQ charges", &
1237 common_iter_levels=1, &
1244 description=
"Controls the printing of the MAO (modified atomic orbital) analysis", &
1246 common_iter_levels=1, &
1249 description=
"Threshold for matrix elements in MAO determination.", &
1250 usage=
"EPS_FILTER reps", repeats=.false., n_var=1, &
1251 default_r_val=1.e-8_dp, type_of_var=
real_t)
1254 CALL keyword_create(keyword, __location__, name=
"REFERENCE_BASIS", &
1255 description=
"Basis set used to construct MAO's.", &
1256 usage=
"REFERENCE_BASIS {ORBITAL,PRIMITIVE,EXTERNAL}", repeats=.false., n_var=1, &
1258 enum_c_vals=
s2a(
"ORBITAL",
"PRIMITIVE",
"EXTERNAL"), &
1259 enum_desc=
s2a(
"Use standard orbital basis set",
"Construct basis from primitives of the orbital basis", &
1260 "Read external basis (MAO)"), &
1265 description=
"Print out MAO reference basis.", &
1266 usage=
"PRINT_BASIS {logical}", repeats=.false., n_var=1, &
1267 default_l_val=.false., lone_keyword_l_val=.true., type_of_var=
logical_t)
1271 description=
"Threshold for gradient in MAO optimization.", &
1272 usage=
"EPS_GRAD reps", repeats=.false., n_var=1, &
1273 default_r_val=1.e-4_dp, type_of_var=
real_t)
1276 CALL keyword_create(keyword, __location__, name=
"EPS_FUNCTION", &
1277 description=
"Threshold for electron defect in MAO optimization.", &
1278 usage=
"EPS_FUNCTION feps", repeats=.false., n_var=1, &
1279 default_r_val=1.e-3_dp, type_of_var=
real_t)
1283 description=
"Maximum allowed iterations for MAO optimization.", &
1284 usage=
"MAX_ITER iter", repeats=.false., n_var=1, &
1289 description=
"Neglect 3 atom terms in MAO analysis.", &
1290 usage=
"NEGLECT_ABC {logical}", repeats=.false., n_var=1, &
1291 default_l_val=.true., lone_keyword_l_val=.true., type_of_var=
logical_t)
1294 CALL keyword_create(keyword, __location__, name=
"AB_THRESHOLD", &
1295 description=
"Threshold for printing of AB shared electron numbers.", &
1296 usage=
"AB_THRESHOLD thr", repeats=.false., n_var=1, &
1297 default_r_val=1.e-2_dp, type_of_var=
real_t)
1300 CALL keyword_create(keyword, __location__, name=
"ABC_THRESHOLD", &
1301 description=
"Threshold for printing of ABC shared electron numbers.", &
1302 usage=
"ABC_THRESHOLD thr", repeats=.false., n_var=1, &
1303 default_r_val=1.e-5_dp, type_of_var=
real_t)
1306 CALL keyword_create(keyword, __location__, name=
"ANALYZE_UNASSIGNED_CHARGE", &
1307 description=
"Calculate atomic contributions to the unassigned charge.", &
1308 usage=
"ANALYZE_UNASSIGNED_CHARGE {logical}", repeats=.false., n_var=1, &
1309 default_l_val=.false., lone_keyword_l_val=.true., type_of_var=
logical_t)
1317 description=
"Controls the printing of the minimal localized basis analysis", &
1319 common_iter_levels=1, &
1322 description=
"Threshold for matrix elements in basis determination.", &
1323 usage=
"EPS_FILTER reps", repeats=.false., n_var=1, &
1324 default_r_val=1.e-8_dp, type_of_var=
real_t)
1327 CALL keyword_create(keyword, __location__, name=
"FULL_ORTHOGONALIZATION", &
1328 description=
"Orthogonalize the localized minimal basis.", &
1329 usage=
"FULL_ORTHOGONALIZATION {logical}", repeats=.false., n_var=1, &
1330 default_l_val=.false., lone_keyword_l_val=.true., type_of_var=
logical_t)
1334 description=
"Calculate Mayer Bond Orders.", &
1335 usage=
"BOND_ORDER {logical}", repeats=.false., n_var=1, &
1336 default_l_val=.false., lone_keyword_l_val=.true., type_of_var=
logical_t)
1340 NULLIFY (sub_print_key)
1342 description=
"Write the minimal basis on Cube files.", &
1345 description=
"The stride (X,Y,Z) used to write the cube file "// &
1346 "(larger values result in smaller cube files). You can provide 3 numbers (for X,Y,Z) or"// &
1347 " 1 number valid for all components.", &
1348 usage=
"STRIDE 2 2 2", n_var=-1, default_i_vals=(/2, 2, 2/), type_of_var=
integer_t)
1352 description=
"Indexes of the atoms minimal basis to be printed as cube files "// &
1353 "This keyword can be repeated several times "// &
1354 "(useful if you have to specify many indexes).", &
1355 usage=
"ATOM_LIST 1 2", &
1356 n_var=-1, type_of_var=
integer_t, repeats=.true.)
1362 NULLIFY (sub_print_key)
1364 description=
"Write the minimal basis in Molden file format, for visualisation.", &
1367 description=
"Specifies the number of significant digits retained. 3 is OK for visualization.", &
1368 usage=
"NDIGITS {int}", &
1373 description=
"Representation of Gaussian-type orbitals", &
1375 enum_c_vals=
s2a(
"CARTESIAN",
"SPHERICAL"), &
1377 "Cartesian Gaussian orbitals. Use with caution", &
1378 "Spherical Gaussian orbitals. Incompatible with VMD"), &
1390 description=
"Controls energy decomposition analysis", &
1392 common_iter_levels=1, &
1394 CALL keyword_create(keyword, __location__, name=
"REFERENCE_ORB_CANONICAL", &
1395 description=
"Use reference orbitals in canonical form.", &
1396 usage=
"REFERENCE_ORB_CANONICAL {logical}", repeats=.false., n_var=1, &
1397 default_l_val=.true., lone_keyword_l_val=.true., type_of_var=
logical_t)
1400 CALL keyword_create(keyword, __location__, name=
"SKIP_LOCALIZATION", &
1401 description=
"Don't localize the MOs.", &
1402 usage=
"SKIP_LOCALIZATION {logical}", repeats=.false., n_var=1, &
1403 default_l_val=.false., lone_keyword_l_val=.true., type_of_var=
logical_t)
1406 CALL keyword_create(keyword, __location__, name=
"DETAILED_ENERGY", &
1407 description=
"Calculate detailed atomic decomposition energies.", &
1408 usage=
"DETAILED_ENERGY {logical}", repeats=.false., n_var=1, &
1409 default_l_val=.false., lone_keyword_l_val=.true., type_of_var=
logical_t)
1412 CALL keyword_create(keyword, __location__, name=
"EWALD_ALPHA_PARAMETER", &
1413 description=
"Calculate Energy Decomposition for a specific alpha value. "// &
1414 "alpha = 1/(2*rc**2), see GTH pseudopotentials.", &
1415 usage=
"EWALD_ALPHA_PARAMETER alpha", repeats=.false., n_var=1, &
1416 default_r_val=0.0_dp, type_of_var=
real_t)
1425 description=
"Controls the printing of the IAO (intrinsic atomic orbital) analysis", &
1427 common_iter_levels=1, &
1430 description=
"Threshold for matrix inversion eigenvalues.", &
1431 usage=
"EPS_SVD reps", repeats=.false., n_var=1, &
1432 default_r_val=0.0_dp, type_of_var=
real_t)
1436 description=
"Threshold in occupation for vectors included.", &
1437 usage=
"EPS_OCC reps", repeats=.false., n_var=1, &
1438 default_r_val=0.0_dp, type_of_var=
real_t)
1441 CALL keyword_create(keyword, __location__, name=
"ATOMIC_CHARGES", &
1442 description=
"Calculate atomic charges from IAO.", &
1443 usage=
"ATOMIC_CHARGES {logical}", repeats=.false., n_var=1, &
1444 default_l_val=.true., lone_keyword_l_val=.true., type_of_var=
logical_t)
1448 NULLIFY (sub_print_key)
1450 description=
"Write the IAO basis in Molden file format, for visualisation.", &
1453 description=
"Specifies the number of significant digits retained. 3 is OK for visualization.", &
1454 usage=
"NDIGITS {int}", &
1459 description=
"Representation of Gaussian-type orbitals", &
1461 enum_c_vals=
s2a(
"CARTESIAN",
"SPHERICAL"), &
1463 "Cartesian Gaussian orbitals. Use with caution", &
1464 "Spherical Gaussian orbitals. Incompatible with VMD"), &
1471 NULLIFY (sub_print_key)
1473 description=
"Controls the printing of the IAO basis "// &
1474 "as *.cube files.", &
1478 description=
"The stride (X,Y,Z) used to write the cube file "// &
1479 "(larger values result in smaller cube files). You can provide 3 numbers (for X,Y,Z) or"// &
1480 " 1 number valid for all components.", &
1481 usage=
"STRIDE 2 2 2", n_var=-1, default_i_vals=(/2, 2, 2/), type_of_var=
integer_t)
1485 description=
"append the cube files when they already exist", &
1486 default_l_val=.false., lone_keyword_l_val=.true.)
1490 description=
"Indices of the atoms to be included in basis CUBE file printing. ", &
1491 usage=
"ATOM_LIST {integer} {integer} .. {integer} ", &
1492 n_var=-1, type_of_var=
integer_t, repeats=.true.)
1498 NULLIFY (sub_print_key)
1500 description=
"Calculates single center expansion of IAOs ", &
1504 description=
"Maximum l quantum number used in the expansion.", &
1505 usage=
"LMAX 2", n_var=1, default_i_val=3, type_of_var=
integer_t)
1509 description=
"Max number of basis functions used in the expansion."// &
1510 " Default is determined by the orbital basis set.", &
1511 usage=
"NBAS 10", n_var=1, default_i_val=-1, type_of_var=
integer_t)
1515 description=
"Append the OCE basis files when it already exists", &
1516 default_l_val=.false., lone_keyword_l_val=.true.)
1522 NULLIFY (sub_print_key)
1524 description=
"Calculate intrinsic bond orbitals using "// &
1525 "localized MOs in IAO basis.", &
1529 CALL keyword_create(keyword, __location__, name=
"LOCALIZATION_OPERATOR", &
1530 description=
"Operator to be optimized for orbital localization", &
1531 enum_c_vals=
s2a(
"PIPEK_MEZEY",
"PIPEK_MEZEY_4",
"L1NORM"), &
1533 enum_desc=
s2a(
"Use Pipek-Mezey operator (order 2)", &
1534 "Use Pipek-Mezey operator (order 4)", &
1539 CALL keyword_create(keyword, __location__, name=
"ENERGY_LOCALIZATION_FUNCTION", &
1540 description=
"Function for energy localization: f(e_i), e_i orbital energy", &
1541 enum_c_vals=
s2a(
"NONE",
"ENERGY",
"OCCUPATION"), &
1543 enum_desc=
s2a(
"Don't use energy localization.", &
1544 "Use orbital energies for localization.", &
1545 "Use occupation numbers for localization."), &
1549 CALL keyword_create(keyword, __location__, name=
"ENERGY_LOCALIZATION_WEIGHT", &
1550 description=
"Weight given to energy localization, using f(e_i) function", &
1551 usage=
"ENERGY_LOCALIZATION_WEIGHT 0.1", n_var=1, &
1552 default_r_val=0.0_dp, type_of_var=
real_t)
1557 NULLIFY (subsection)
1559 description=
"Calculation and printing of centers and spreads "// &
1560 "of localized orbitals.", &
1563 CALL keyword_create(keyword, __location__, name=
"POSITION_OPERATOR_BERRY", &
1564 description=
"Use Berry phase position operator.", &
1565 usage=
"POSITION_OPERATOR_BERRY T", n_var=1, &
1566 default_l_val=.true., lone_keyword_l_val=.true.)
1572 NULLIFY (subsection)
1574 description=
"Write the IBO orbitals in Molden file format, for visualisation.", &
1577 description=
"Specifies the number of significant digits retained. 3 is OK for visualization.", &
1578 usage=
"NDIGITS {int}", &
1583 description=
"Representation of Gaussian-type orbitals", &
1585 enum_c_vals=
s2a(
"CARTESIAN",
"SPHERICAL"), &
1587 "Cartesian Gaussian orbitals. Use with caution", &
1588 "Spherical Gaussian orbitals. Incompatible with VMD"), &
1595 NULLIFY (subsection)
1597 description=
"Controls the printing of the IBO orbitals "// &
1598 "as *.cube files.", &
1602 description=
"The stride (X,Y,Z) used to write the cube file "// &
1603 "(larger values result in smaller cube files). You can provide 3 numbers (for X,Y,Z) or"// &
1604 " 1 number valid for all components.", &
1605 usage=
"STRIDE 2 2 2", n_var=-1, default_i_vals=(/2, 2, 2/), type_of_var=
integer_t)
1609 description=
"append the cube files when they already exist", &
1610 default_l_val=.false., lone_keyword_l_val=.true.)
1614 description=
"Indices of the orbitals to be included in IBO CUBE file printing. ", &
1615 usage=
"STATE_LIST {integer} {integer} .. {integer} ", &
1616 n_var=-1, type_of_var=
integer_t, repeats=.true.)
1630 description=
"Controls the printing of the DOS from the density matrix. "// &
1631 "This allows the calculation of the DOS even in density matrix based "// &
1632 "REAL_TIME_PROPAGATION and LS_SCF. "// &
1633 "However, it requires a cubically scaling diagonalization of the Hamiltonian. "// &
1634 "Hartree-Fock NYI, values will be wrong. "// &
1635 "Careful, the orbitals in rtp/emd are not actually eigenstates of the Hamiltonian. "// &
1636 "Assumes absence of spin polarization (so far).", &
1638 each_iter_names=
s2a(
"MD"), each_iter_values=(/100/), &
1641 description=
"The number of energy windows.", &
1642 usage=
"N_WINDOWS 200", &
1647 description=
"Filtering threshold for sparse matrix operations.", &
1648 usage=
"EPS_FILTER 1.0E-6", &
1649 default_r_val=1.0e-14_dp)
1652 CALL keyword_create(keyword, __location__, name=
"RESTRICT_RANGE", &
1653 description=
"Restricts the energy windows to states close to the fermi level", &
1654 usage=
"RESTRICT_RANGE .TRUE.", &
1655 default_l_val=.false., lone_keyword_l_val=.true.)
1659 description=
"If the RESTRICT_RANGE keyword is set, then all energy widnows will"// &
1660 " be placed in an interval from from the fermi level minus to the fermi level plus this keyword", &
1662 default_r_val=1.0_dp)
1666 description=
"Print the energy windows to cube files", &
1667 usage=
"PRINT_CUBES .TRUE.", &
1668 default_l_val=.false., lone_keyword_l_val=.true.)
1672 description=
"The stride (X,Y,Z) used to write the energy windows cube files (if enabled) "// &
1673 "(larger values result in smaller cube files). You can provide 3 numbers (for X,Y,Z) or"// &
1674 " 1 number valid for all components.", &
1675 usage=
"STRIDE 2 2 2", n_var=-1, default_i_vals=(/2, 2, 2/), type_of_var=
integer_t)
1683 description=
"Write the KS matrix in CSR format into a file.", &
1686 description=
"Threshold on the absolute value of the elements to be printed out. "// &
1687 "In CP2K all the elements in a (atomic) matrix block are considered non-zero, "// &
1688 "if the block contains at least one non-zero element.", &
1689 usage=
"THRESHOLD {real}", &
1691 default_r_val=0.0_dp)
1694 CALL keyword_create(keyword, __location__, name=
"Upper_triangular", &
1695 description=
"Print only the upper triangular part of the matrix. ", &
1696 usage=
"UPPER_TRIANGULAR {logical}", &
1698 default_l_val=.false., &
1699 lone_keyword_l_val=.true.)
1703 description=
"Whether or not to generate the file in a binary format. ", &
1704 usage=
"BINARY {logical}", &
1706 default_l_val=.false., &
1707 lone_keyword_l_val=.true.)
1711 description=
"Print the KS matrix in real-space instead of k-space.. ", &
1712 usage=
"REAL_SPACE {logical}", &
1714 default_l_val=.false., &
1715 lone_keyword_l_val=.true.)
1723 description=
"Write the overlap matrix in CSR format into a file.", &
1726 description=
"Threshold on the absolute value of the elements to be printed out. "// &
1727 "In CP2K all the elements in a (atomic) matrix block are considered non-zero, "// &
1728 "if the block contains at least one non-zero element.", &
1729 usage=
"THRESHOLD {real}", &
1731 default_r_val=0.0_dp)
1734 CALL keyword_create(keyword, __location__, name=
"Upper_triangular", &
1735 description=
"Print only the upper triangular part of the matrix. ", &
1736 usage=
"UPPER_TRIANGULAR {logical}", &
1738 default_l_val=.false., &
1739 lone_keyword_l_val=.true.)
1743 description=
"Whether or not to generate the file in a binary format. ", &
1744 usage=
"BINARY {logical}", &
1746 default_l_val=.false., &
1747 lone_keyword_l_val=.true.)
1751 description=
"Print the overlap matrix in real-space instead of k-space.. ", &
1752 usage=
"REAL_SPACE {logical}", &
1754 default_l_val=.false., &
1755 lone_keyword_l_val=.true.)
1763 description=
"Writes an (upper-triangular) adjacency matrix indicating the "// &
1764 "interaction between atoms (according to overlapping basis functions). The "// &
1765 "columns are: iatom, jatom, ikind, jkind; where iatom and jatom are the atom "// &
1766 "indices (based on the coordinate file), ikind and jkind are the atomic kinds "// &
1767 "(indeces as shown in the ATOMIC KIND INFORMATION section of a CP2K output). ", &
1774 print_key, __location__, name=
"XRAY_DIFFRACTION_SPECTRUM", &
1775 description=
"Calculate and print the coherent X-ray "// &
1776 "diffraction spectrum", &
1781 keyword, __location__, &
1783 variants=(/
"Q_MAXIMUM"/), &
1784 description=
"Maximum Q value calculated for the spectrum", &
1785 usage=
"Q_MAX {real}", &
1790 unit_str=
"angstrom^-1"), &
1791 unit_str=
"angstrom^-1")
1798 description=
"Calculate and print the electric field gradients "// &
1799 "at atomic positions", &
1801 filename=
"__STD_OUT__")
1804 name=
"INTERPOLATION", &
1805 description=
"Use interpolation method from real space grid", &
1806 usage=
"INTERPOLATION {logical}", &
1809 default_l_val=.false., lone_keyword_l_val=.true.)
1814 name=
"GSPACE_SMOOTHING", &
1815 description=
"Use a G-space smoothing function", &
1816 usage=
"GSPACE_SMOOTHING cutoff {real}, width {real}", &
1818 n_var=2, default_r_vals=(/-1._dp, -1._dp/), &
1825 description=
"Print additional debug output", &
1826 usage=
"DEBUG {logical}", &
1829 default_l_val=.false., lone_keyword_l_val=.true.)
1841 description=
"Print the two quantities needed in the basis molopt generation:"// &
1842 " total energy and condition number of the overlap matrix (S matrix)", &
1844 filename=
"__STD_OUT__")
1849 description=
"Calculate and print the EPR hyperfine coupling tensor"// &
1850 " at atomic positions", &
1852 filename=
"__STD_OUT__")
1855 name=
"INTERACTION_RADIUS", &
1856 description=
"Radius of interaction for EPR hyperfine tensor calculation", &
1857 usage=
"INTERACTION_RADIUS radius {real}", &
1859 n_var=1, default_r_val=10._dp, &
1868 description=
"Optimize the exponents of the LRI basis set", &
1870 filename=
"OPTIMIZED_LRI_BASIS")
1875 print_key, __location__, name=
"PLUS_U", &
1876 description=
"Controls the printing for the DFT+U methods", &
1878 filename=
"__STD_OUT__", &
1879 each_iter_names=
s2a(
"QS_SCF"), &
1880 each_iter_values=(/0/), &
1886 print_key, __location__, name=
"CHARGEMOL", &
1887 description=
"Write .wfx input file for Chargemol", &
1889 filename=
"CHARGEMOL", &
1891 CALL keyword_create(keyword, __location__, name=
"BACKUP_COPIES", &
1892 description=
"Specifies the maximum number of backup copies.", &
1893 usage=
"BACKUP_COPIES {int}", &
1898 description=
"Write information about cell periodicity.", &
1899 usage=
"PERIODIC {LOGICAL}", &
1900 default_l_val=.false., lone_keyword_l_val=.true.)
1907 print_key, __location__, name=
"SCCS", &
1908 description=
"Controls the printing for the SCCS models", &
1910 filename=
"__STD_OUT__", &
1911 each_iter_names=
s2a(
"QS_SCF"), &
1912 each_iter_values=(/0/), &
1915 NULLIFY (sub_print_key)
1918 sub_print_key, __location__, name=
"DENSITY_GRADIENT", &
1919 description=
"Controls the printing of the cube files with "// &
1920 "the norm of the density gradient |∇ρ| "// &
1921 "used by the SCCS model.", &
1924 each_iter_names=
s2a(
"QS_SCF"), &
1925 each_iter_values=(/0/))
1927 description=
"The stride (X,Y,Z) used to write the cube file "// &
1928 "(larger values result in smaller cube files). You can provide 3 "// &
1929 "numbers (for X,Y,Z) or 1 number valid for all components.", &
1931 default_i_vals=(/2, 2, 2/), &
1937 description=
"Append the cube files when they already exist", &
1938 default_l_val=.false., &
1939 lone_keyword_l_val=.true., &
1947 sub_print_key, __location__, name=
"DIELECTRIC_FUNCTION", &
1948 description=
"Controls the printing of the cube files with "// &
1949 "the dielectric function used by the SCCS model. "// &
1950 "This function determines the cavity formed by a solute in "// &
1951 "a solvent and thus it can be used for the visualisaton of the cavity.", &
1954 each_iter_names=
s2a(
"QS_SCF"), &
1955 each_iter_values=(/0/), &
1958 description=
"The stride (X,Y,Z) used to write the cube file "// &
1959 "(larger values result in smaller cube files). You can provide 3 "// &
1960 "numbers (for X,Y,Z) or 1 number valid for all components.", &
1962 default_i_vals=(/2, 2, 2/), &
1968 description=
"Append the cube files when they already exist", &
1969 default_l_val=.false., &
1970 lone_keyword_l_val=.true., &
1978 sub_print_key, __location__, name=
"TOTAL_CHARGE_DENSITY", &
1979 description=
"Controls the printing of the cube files with the "// &
1980 "total charge density $\rho^\text{tot}$ used by the SCCS model.", &
1983 each_iter_names=
s2a(
"QS_SCF"), &
1984 each_iter_values=(/0/), &
1987 description=
"The stride (X,Y,Z) used to write the cube file "// &
1988 "(larger values result in smaller cube files). You can provide 3 "// &
1989 "numbers (for X,Y,Z) or 1 number valid for all components.", &
1991 default_i_vals=(/2, 2, 2/), &
1997 description=
"Append the cube files when they already exist", &
1998 default_l_val=.false., &
1999 lone_keyword_l_val=.true., &
2007 sub_print_key, __location__, name=
"POLARISATION_CHARGE_DENSITY", &
2008 description=
"Controls the printing of the cube files with the "// &
2009 "polarisation charge density $\rho^\text{pol}$ used by the SCCS model with the "// &
2010 "total charge density $\rho^\text{tot} = \rho^\text{sol} + \rho^\text{pol}", &
2013 each_iter_names=
s2a(
"QS_SCF"), &
2014 each_iter_values=(/0/), &
2017 description=
"The stride (X,Y,Z) used to write the cube file "// &
2018 "(larger values result in smaller cube files). You can provide 3 "// &
2019 "numbers (for X,Y,Z) or 1 number valid for all components.", &
2021 default_i_vals=(/2, 2, 2/), &
2027 description=
"Append the cube files when they already exist", &
2028 default_l_val=.false., &
2029 lone_keyword_l_val=.true., &
2037 sub_print_key, __location__, name=
"POLARISATION_POTENTIAL", &
2038 description=
"Controls the printing of the cube files with the "// &
2039 "polarisation potential $\phi^\text{pol}$ used by the SCCS model with the "// &
2040 "total potential $\phi^\text{tot} = \phi^\text{sol} + \phi^\text{pol}$", &
2043 each_iter_names=
s2a(
"QS_SCF"), &
2044 each_iter_values=(/0/), &
2047 description=
"The stride (X,Y,Z) used to write the cube file "// &
2048 "(larger values result in smaller cube files). You can provide 3 "// &
2049 "numbers (for X,Y,Z) or 1 number valid for all components.", &
2051 default_i_vals=(/2, 2, 2/), &
2057 description=
"Append the cube files when they already exist", &
2058 default_l_val=.false., &
2059 lone_keyword_l_val=.true., &
2076 SUBROUTINE create_bandstructure_section(section)
2082 cpassert(.NOT.
ASSOCIATED(section))
2083 CALL section_create(section, __location__, name=
"BAND_STRUCTURE", &
2084 description=
"Specifies the k-points used in band structure calculation.", &
2085 n_keywords=0, n_subsections=0, repeats=.false.)
2089 description=
"File name used for band structure", &
2090 usage=
"FILE_NAME <filename>", default_c_val=
"", &
2091 type_of_var=
char_t, n_var=1)
2096 variants=(/
"ADDED_BANDS"/), &
2097 description=
"Number of MOs/Bands added to the Band Structure calculation.", &
2102 NULLIFY (subsection)
2107 END SUBROUTINE create_bandstructure_section
2113 SUBROUTINE create_mo_cubes_section(print_key)
2121 description=
"Controls the printing of the molecular orbitals (MOs) as cube files."// &
2122 " It can be used during a Real Time calculation to print the MOs."// &
2123 " In this case, the density corresponding to the time dependent MO is printed"// &
2124 " instead of the wave-function.", &
2127 description=
"The stride (X,Y,Z) used to write the cube file "// &
2128 "(larger values result in smaller cube files). You can provide 3 numbers (for X,Y,Z) or"// &
2129 " 1 number valid for all components.", &
2130 usage=
"STRIDE 2 2 2", n_var=-1, default_i_vals=(/2, 2, 2/), type_of_var=
integer_t)
2135 description=
"If the MO cube file should be written. If false, the eigenvalues are still computed."// &
2136 " Can also be useful in combination with STM calculations", &
2137 default_l_val=.true., lone_keyword_l_val=.true.)
2142 description=
"If the printkey is activated controls the number of lumos"// &
2143 " that are printed and dumped as a cube (-1=all)", &
2148 keyword, __location__, name=
"nhomo", &
2149 description=
"If the printkey is activated controls the number of homos that dumped as a cube (-1=all),"// &
2150 " eigenvalues are always all dumped", &
2156 keyword, __location__, name=
"homo_list", &
2157 description=
"If the printkey is activated controls the index of homos dumped as a cube,"// &
2158 " eigenvalues are always all dumped. It overrides nhomo.", &
2159 usage=
"HOMO_LIST {integer} {integer} .. {integer} ", type_of_var=
integer_t, &
2160 n_var=-1, repeats=.true.)
2165 description=
"append the cube files when they already exist", &
2166 default_l_val=.false., lone_keyword_l_val=.true.)
2170 END SUBROUTINE create_mo_cubes_section
2176 SUBROUTINE create_dos_section(print_key)
2185 description=
"Print Density of States (DOS) (only available states from SCF)", &
2189 description=
"Append the DOS obtained at different iterations to the output file. "// &
2190 "By default the file is overwritten", &
2191 usage=
"APPEND", default_l_val=.false., &
2192 lone_keyword_l_val=.true.)
2197 description=
"Histogramm energy spacing.", &
2198 usage=
"DELTA_E 0.0005", type_of_var=
real_t, default_r_val=0.001_dp)
2203 description=
"Specify the number of digits used to print density and occupation", &
2209 description=
"Specify a Monkhorst-Pack grid with which to compute the density of states. "// &
2210 "Works only for a k-point calculation", &
2211 usage=
"MP_GRID {integer} {integer} {integer}", default_i_vals=(/-1/), &
2216 END SUBROUTINE create_dos_section
2229 NULLIFY (subsection)
2233 description=
"Print out the DOS projected per kind and angular momentum ", &
2236 description=
"Print out PDOS distinguishing all angular momentum components.", &
2237 usage=
"COMPONENTS", default_l_val=.false., &
2238 lone_keyword_l_val=.true.)
2242 description=
"Append the PDOS obtained at different iterations to the PDOS output file. "// &
2243 "By default the file is overwritten", &
2244 usage=
"APPEND", default_l_val=.false., &
2245 lone_keyword_l_val=.true.)
2249 keyword, __location__, name=
"NLUMO", &
2250 description=
"Number of virtual orbitals to be added to the MO set (-1=all)."//
newline// &
2251 "CAUTION: Setting this value to be higher than the "// &
2252 "number of states present may cause a Cholesky error."//
newline// &
2253 "This keyword only applies to MO sets, not to XAS_TDP states.", &
2254 usage=
"NLUMO integer", default_i_val=0)
2257 CALL keyword_create(keyword, __location__, name=
"OUT_EACH_STATE", &
2258 variants=(/
"OUT_EACH_MO"/), &
2259 description=
"Output on the status of the calculation every OUT_EACH_MO states. If -1 no output", &
2260 usage=
"OUT_EACH_STATE integer", default_i_val=-1)
2265 description=
"Controls the printing of local PDOS, projected on subsets"// &
2266 " of atoms given through lists", &
2267 n_keywords=4, n_subsections=0, repeats=.true.)
2269 description=
"Print out PDOS distinguishing all angular momentum components.", &
2270 usage=
"COMPONENTS", default_l_val=.false., &
2271 lone_keyword_l_val=.true.)
2276 description=
"Specifies a list of indexes of atoms where to project the DOS ", &
2277 usage=
"LIST {integer} {integer} .. {integer} ", type_of_var=
integer_t, &
2278 n_var=-1, repeats=.true.)
2286 description=
"Controls the printing of local PDOS, projected on 3D volume in real space,"// &
2287 " the volume is defined in terms of position with respect to atoms in the lists", &
2288 n_keywords=4, n_subsections=0, repeats=.true.)
2291 description=
"Specifies a list of indexes of atoms used to define the real space volume ", &
2292 usage=
"LIST {integer} {integer} .. {integer} ", type_of_var=
integer_t, &
2293 n_var=-1, repeats=.true.)
2298 description=
"range of positions in Cartesian direction x: all grid points within "// &
2299 "this range from at least one atom of the list are considered", &
2300 usage=
"XRANGE -10.0 10.0", unit_str=
"angstrom", n_var=2, type_of_var=
real_t)
2304 description=
"range of positions in Cartesian direction y: all grid points within "// &
2305 "this range from at least one atom of the list are considered", &
2306 usage=
"YRANGE -10.0 10.0", unit_str=
"angstrom", n_var=2, type_of_var=
real_t)
2310 description=
"range of positions in Cartesian direction z: all grid points within "// &
2311 "this range from at least one atom of the list are considered", &
2312 usage=
"ZRANGE -10.0 10.0", unit_str=
"angstrom", n_var=2, type_of_var=
real_t)
2317 description=
"Only project states with the energy values in the given interval. "// &
2318 "Default is all states.", &
2319 usage=
"ERANGE -1.0 1.0", unit_str=
"hartree", n_var=2, type_of_var=
real_t)
2332 SUBROUTINE create_wannier_section(print_key)
2341 description=
"Interface to Wannier90 code. (EXPERIMENTAL)", &
2345 description=
"The seedname for the Wannier90 calculation (body of filenames).", &
2346 usage=
"SEED_NAME filename", default_c_val=
"wannier90", &
2347 n_var=1, type_of_var=
char_t)
2352 description=
"The dimensions of the Monkhorst-Pack k-point grid. ", &
2353 usage=
"MP_GRID 6 6 6", n_var=-1, default_i_vals=(/10, 10, 10/), type_of_var=
integer_t)
2358 variants=(/
"ADDED_BANDS"/), &
2359 description=
"Number of MOs/Bands added to the Band Structure calculation.", &
2364 CALL keyword_create(keyword, __location__, name=
"EXCLUDE_BANDS", &
2365 description=
"List of Bands excluded in the Wannier calculation.", &
2366 usage=
"EXCLUDE_BANDS b1 b2 ...", n_var=-1, repeats=.true., &
2371 CALL keyword_create(keyword, __location__, name=
"WANNIER_FUNCTIONS", &
2372 description=
"Number of Wannier functions to be calculated. ", &
2373 usage=
"WANNIER_FUNCTIONS 6", n_var=1, default_i_val=0, &
2378 END SUBROUTINE create_wannier_section
2384 SUBROUTINE create_stm_section(print_key)
2392 description=
"Controls the printing of cubes for the generation of STM images.", &
2395 description=
"The stride (X,Y,Z) used to write the cube file "// &
2396 "(larger values result in smaller cube files). You can provide 3 numbers (for X,Y,Z) or"// &
2397 " 1 number valid for all components.", &
2398 usage=
"STRIDE 2 2 2", n_var=-1, default_i_vals=(/2, 2, 2/), type_of_var=
integer_t)
2403 description=
"If the printkey is activated controls the number of additional lumos"// &
2404 " that are computed to be able to reproduce STM images obtained"// &
2405 " from positive bias (imaging unoccupied states)", &
2411 description=
"Bias energy for scanning tunneling microscopy (STM) image generation. "// &
2412 "Orbital densities are summed according to the bias energy. "// &
2413 "For negative values, states in the range ]EF+bias,EF] are summed, "// &
2414 "While positive values sum states in the range [EF,EF+bias[. "// &
2415 "If positive biases are used, sufficiently many unoccupied stated "// &
2416 "(see ADDED_MOS and NLUMO ) should be computed.", &
2417 n_var=-1, type_of_var=
real_t, default_r_vals=(/0.0_dp/), unit_str=
'eV')
2422 description=
"Tip orbital symmetry in Tersoff-Hamann approximation to compute STM images", &
2424 default_i_val=
orb_s, &
2425 usage=
"TH_TORB s dz2", &
2426 enum_c_vals=
s2a(
"S",
"PX",
"PY",
"PZ",
"DXY",
"DYZ",
"DZX",
"DX2",
"DY2",
"DZ2"), &
2427 enum_i_vals=(/
orb_s,
orb_px,
orb_py,
orb_pz,
orb_dxy,
orb_dyz,
orb_dzx,
orb_dx2,
orb_dy2,
orb_dz2/), &
2428 enum_desc=
s2a(
"s orbital",
"px orbital",
"py orbital",
"pz orbital", &
2429 "dxy orbital",
"dyz orbital",
"dzx orbital",
"x^2 orbital",
"y^2 orbital",
"z^2 orbital"))
2434 description=
"By default the reference energy is the Fermi energy. In order to compare"// &
2435 " with STS experiments, where specific energy ranges are addressed, here"// &
2436 " one can set a different reference energy."// &
2437 " The energy range is anyway controlled by the BIAS", &
2438 type_of_var=
real_t, default_r_val=0.0_dp, unit_str=
'eV')
2443 description=
"append the cube files when they already exist", &
2444 default_l_val=.false., lone_keyword_l_val=.true.)
2448 END SUBROUTINE create_stm_section
2454 SUBROUTINE create_wfn_mix_section(section)
2461 NULLIFY (subsection)
2464 cpassert(.NOT.
ASSOCIATED(section))
2467 description=
"A section that allows manipulation of the MO coeffs, "// &
2468 "e.g. for changing a ground state into an excited state. "// &
2469 "Starting from a copy of the original MOs, changes can be made "// &
2470 "by adding linear combinations of HOMO/LUMO of the original MOs to the result MOs. "// &
2471 "This method is called after an SCF optimization or before an RTP run if "// &
2472 "INITIAL_WFN=RESTART_WFN. Note that if called after an SCF optimization, a restart file "// &
2473 "with the mixed MOs is saved. This is not the case for an RTP with "// &
2474 "INITIAL_WFN=RESTART_WFN.", &
2475 n_keywords=1, n_subsections=0, repeats=.false.)
2477 CALL keyword_create(keyword, __location__, name=
"OVERWRITE_MOS", &
2478 description=
"If set to True, the active molecular orbitals in memory will be replaced by the mixed wfn "// &
2479 "at the end of the wfn mixing procedure. For instance, you can then use this new set of MOs to perform "// &
2480 "RTP or EMD directly. Note that in the case of an RTP run with INITIAL_WFN=RESTART_WFN, the OVERWRITE_MOS "// &
2481 "keyword is not used.", &
2482 default_l_val=.false., lone_keyword_l_val=.true.)
2487 description=
"Update a result MO with with a linear combination of original MOs."// &
2488 " This section can be repeated to build arbitrary linear combinations using repeatedly y=a*y+b*x. "// &
2489 "RESULT is (y), RESULT_SCALE is (a), ORIG is (x), ORIG_SCALE is (b)", &
2490 n_keywords=1, n_subsections=0, repeats=.true.)
2492 CALL keyword_create(keyword, __location__, name=
"RESULT_MO_INDEX", &
2493 description=
"Index of the MO (y) to be modified. Counting down in energy: "// &
2494 "set to 1 for the highest MO, to 3 for the highest MO-2.", &
2495 usage=
"RESULT_MO_INDEX 1", type_of_var=
integer_t, default_i_val=0)
2499 CALL keyword_create(keyword, __location__, name=
"RESULT_MARKED_STATE", &
2500 description=
"Specifies the MO according to "// &
2501 "the marks set in MOLECULAR_STATES. The value corresponds to the repetition "// &
2502 "of MARK_STATES in MOLECULAR_STATES", &
2503 usage=
"RESULT_MARKED_STATE 1", type_of_var=
integer_t, default_i_val=0)
2507 CALL keyword_create(keyword, __location__, name=
"REVERSE_MO_INDEX", &
2508 description=
"Reverses the index order of the OCCUPIED and EXTERNAL MOs. With this keyword "// &
2509 "ORIG_MO_INDEX/RESULT_MO_INDEX 1 point to the lowest energy MO (instead of the highest) "// &
2510 "and counts up in energy. The VIRTUAL MOs indexing is unchanged.", &
2511 default_l_val=.false., lone_keyword_l_val=.true.)
2515 CALL keyword_create(keyword, __location__, name=
"RESULT_SPIN_INDEX", &
2516 description=
"Spin of the MO (y) to be modified.", &
2517 enum_c_vals=
s2a(
"Alpha",
"Beta"), &
2518 enum_i_vals=(/1, 2/), &
2520 enum_desc=
s2a(
"Majority spin",
"Minority spin"))
2524 CALL keyword_create(keyword, __location__, name=
"RESULT_SCALE", &
2525 description=
"Scaling factor of the result variable (a).", &
2526 usage=
"RESULT_SCALE 0.0", type_of_var=
real_t)
2530 CALL keyword_create(keyword, __location__, name=
"ORIG_MO_INDEX", &
2531 description=
"Index of the original MO (x). "// &
2532 "If ORIG_TYPE is OCCUPIED, it counts down in energy: set to 1 to point to "// &
2533 "the highest MO and to 3 for the highest MO-2. "// &
2534 "If ORIG_TYPE is VIRTUAL, it counts up in energy: set to 1 to point to "// &
2535 "the lowest virtual MO and to 3 for the lowest MO+2. "// &
2536 "If ORIG_TYPE is EXTERNAL, it counts down in energy for the external "// &
2537 "set of MOs: set to 1 to point to the highest MO and to 3 for the highest MO-2. "// &
2538 "Do not set to zero or negative values.", &
2539 usage=
"ORIG_MO_INDEX 1", type_of_var=
integer_t, default_i_val=0)
2543 CALL keyword_create(keyword, __location__, name=
"ORIG_MARKED_STATE", &
2544 description=
"Specifies the MO according to "// &
2545 "the marks set in MOLECULAR_STATES. The value corresponds to the repetition "// &
2546 "of MARK_STATES in MOLECULAR_STATES", &
2547 usage=
"ORIG_MARKED_STATE 1", type_of_var=
integer_t, default_i_val=0)
2551 CALL keyword_create(keyword, __location__, name=
"ORIG_SPIN_INDEX", &
2552 description=
"Spin of the MO (x) to be modified.", &
2553 enum_c_vals=
s2a(
"Alpha",
"Beta"), &
2554 enum_i_vals=(/1, 2/), &
2556 enum_desc=
s2a(
"Majority spin",
"Minority spin"))
2561 description=
"Scaling factor of the original variable (b).", &
2562 usage=
"ORIG_SCALE 0.0", type_of_var=
real_t)
2567 description=
"Type of the original MO. Note that if ADDED_MOS was used in the "// &
2568 "SCF construction of the MO matrix, these extra MOs are also treated as OCCUPIED. ", &
2569 enum_c_vals=
s2a(
"OCCUPIED",
"VIRTUAL",
'EXTERNAL'), &
2570 usage=
"ORIG_TYPE OCCUPIED", &
2572 enum_desc=
s2a(
"The original MO is the result of the SCF procedure. This can also contain "// &
2573 "unoccupied MOs if the SCF%ADDED_MOS keyword was used.", &
2574 "The original MO is taken from the result of additional MOs calculated a "// &
2575 "posteriori of the SCF by request of the user. E.g. by specifying print%mo_cubes%nlumo. ", &
2576 "The orginal MO is from an external .wfn file. Use the keyword "// &
2577 "ORIG_EXT_FILE_NAME to define its name."), &
2582 CALL keyword_create(keyword, __location__, name=
"ORIG_EXT_FILE_NAME", &
2583 description=
"Name of the wavefunction file to read the original MO from. "// &
2584 "For instance, a restart wfn file from SCF calculation or an excited state from XAS_TDP calculation. "// &
2585 "If no file is specified, the run will crash. "// &
2586 "Currently, a RTP restart file (.rtpwfn) cannot be used as reference. "// &
2587 "Currently, this file SHALL have the basis set, number of MO and the same number of spin as the one "// &
2588 "from the SCF cycle.", &
2589 usage=
"ORIG_EXT_FILE_NAME <FILENAME>", &
2590 default_lc_val=
"EMPTY")
2597 END SUBROUTINE create_wfn_mix_section
2606 SUBROUTINE create_implicit_psolver_section(section)
2612 cpassert(.NOT.
ASSOCIATED(section))
2613 CALL section_create(section, __location__, name=
"IMPLICIT_PSOLVER", &
2614 description=
"Controls printing of cube files for data from the implicit "// &
2615 "(generalized) Poisson solver.", &
2617 n_keywords=0, n_subsections=3, repeats=.false.)
2619 NULLIFY (keyword, print_key)
2623 description=
"Controls the printing of a cube file with dielectric constant from "// &
2624 "the implicit (generalized) Poisson solver.", &
2628 description=
"The stride (X,Y,Z) used to write the cube file "// &
2629 "(larger values result in smaller cube files). You can provide 3 numbers (for X,Y,Z) or"// &
2630 " 1 number valid for all components.", &
2631 usage=
"STRIDE 2 2 2", n_var=-1, default_i_vals=(/2, 2, 2/), type_of_var=
integer_t)
2635 description=
"append the cube files when they already exist", &
2636 default_l_val=.false., lone_keyword_l_val=.true.)
2645 print_key, __location__,
"DIRICHLET_BC_CUBE", &
2646 description=
"Controls the printing of cube files with unit step functions (constraints)"// &
2647 " representing Dirichlet-type (boundary) regions defined in the implicit (generalized) Poisson"// &
2648 " solver section. The regions remain unchanged throughout the calculations. If the Dirichlet"// &
2649 " regions are relatively large and/or the number of partitions is quite high, in order to save memory,"// &
2650 " generate the cube files in early steps and perform the rest of the calculations with this keyword"// &
2655 description=
"Print tiles that tessellate the Dirichlet regions into cube files. If TRUE, "// &
2656 "generates cube files as many as the total number of tiles.", &
2657 usage=
"TILE_CUBES <logical>", &
2658 default_l_val=.false., lone_keyword_l_val=.true.)
2662 description=
"The stride (X,Y,Z) used to write the cube file "// &
2663 "(larger values result in smaller cube files). You can provide 3 numbers (for X,Y,Z) or"// &
2664 " 1 number valid for all components.", &
2665 usage=
"STRIDE 2 2 2", n_var=-1, default_i_vals=(/2, 2, 2/), type_of_var=
integer_t)
2669 description=
"append the cube files when they already exist", &
2670 default_l_val=.false., lone_keyword_l_val=.true.)
2679 description=
"Controls the printing of cube files with penalty charges induced to "// &
2680 "Dirichlet regions by Lagrange multipliers (implicit Poisson solver).", &
2684 description=
"The stride (X,Y,Z) used to write the cube file "// &
2685 "(larger values result in smaller cube files). You can provide 3 numbers (for X,Y,Z) or"// &
2686 " 1 number valid for all components.", &
2687 usage=
"STRIDE 2 2 2", n_var=-1, default_i_vals=(/2, 2, 2/), type_of_var=
integer_t)
2691 description=
"append the cube files when they already exist", &
2692 default_l_val=.false., lone_keyword_l_val=.true.)
2699 END SUBROUTINE create_implicit_psolver_section
2712 cpassert(.NOT.
ASSOCIATED(section))
2713 CALL section_create(section, __location__, name=
"interpolator", &
2714 description=
"controls the interpolation for the G-space term", &
2715 n_keywords=5, n_subsections=0, repeats=.false.)
2717 NULLIFY (keyword, print_key)
2719 CALL keyword_create(keyword, __location__, name=
"aint_precond", &
2720 description=
"the approximate inverse to use to get the starting point"// &
2721 " for the linear solver of the spline3 methods", &
2722 usage=
"aint_precond copy", &
2724 enum_c_vals=
s2a(
"copy",
"spl3_nopbc_aint1",
"spl3_nopbc_precond1", &
2725 "spl3_nopbc_aint2",
"spl3_nopbc_precond2",
"spl3_nopbc_precond3"), &
2732 description=
"The preconditioner used"// &
2733 " for the linear solver of the spline3 methods", &
2734 usage=
"precond copy", &
2736 enum_c_vals=
s2a(
"copy",
"spl3_nopbc_aint1",
"spl3_nopbc_precond1", &
2737 "spl3_nopbc_aint2",
"spl3_nopbc_precond2",
"spl3_nopbc_precond3"), &
2744 description=
"accuracy on the solution for spline3 the interpolators", &
2745 usage=
"eps_x 1.e-15", default_r_val=1.e-10_dp)
2750 description=
"accuracy on the residual for spline3 the interpolators", &
2751 usage=
"eps_r 1.e-15", default_r_val=1.e-10_dp)
2756 variants=(/
'maxiter'/), &
2757 description=
"the maximum number of iterations", &
2758 usage=
"max_iter 200", default_i_val=100)
2764 description=
"if convergence information about the linear solver"// &
2765 " of the spline methods should be printed", &
2767 each_iter_values=(/10/), filename=
"__STD_OUT__", &
integer, parameter, public basis_sort_zet
integer, parameter, public basis_sort_default
collects all references to literature in CP2K as new algorithms / method are included from literature...
integer, save, public vandevondele2005b
integer, save, public blochl1995
integer, save, public guidon2010
integer, save, public vandevondele2003
integer, save, public bengtsson1999
integer, save, public kunert2003
integer, save, public iannuzzi2007
integer, save, public yin2017
integer, save, public stewart2007
integer, save, public holmberg2017
integer, save, public vandevondele2006
integer, save, public golze2017b
integer, save, public caldeweyher2020
integer, save, public krack2002
integer, save, public vandevondele2005a
integer, save, public caldeweyher2019
integer, save, public kuhne2007
integer, save, public golze2017a
integer, save, public schiffmann2015
integer, save, public lippert1999
integer, save, public eriksen2020
integer, save, public dewar1977
integer, save, public avezac2005
integer, save, public vanvoorhis2015
integer, save, public repasky2002
integer, save, public weber2008
integer, save, public andreussi2012
integer, save, public iannuzzi2006
integer, save, public iannuzzi2005
integer, save, public rocha2006
integer, save, public lippert1997
integer, save, public holmberg2018
integer, save, public fattebert2002
integer, save, public andermatt2016
integer, save, public merlot2014
integer, save, public thiel1992
integer, save, public pracht2019
integer, save, public ehrhardt1985
integer, save, public becke1988b
integer, save, public perdew1981
integer, save, public knizia2013
integer, save, public shigeta2001
integer, save, public stewart1982
integer, save, public schenter2008
integer, save, public lu2004
integer, save, public dudarev1997
integer, save, public heinzmann1976
integer, save, public brelaz1979
integer, save, public krack2000
integer, save, public dewar1985
integer, save, public banihashemian2016
integer, save, public stewart1989
integer, save, public brehm2018
integer, save, public kolafa2004
integer, save, public dudarev1998
routines to handle the output, The idea is to remove the decision of wheter to output and what to out...
integer, parameter, public debug_print_level
integer, parameter, public low_print_level
integer, parameter, public medium_print_level
integer, parameter, public high_print_level
integer, parameter, public add_last_numeric
integer, parameter, public silent_print_level
subroutine, public cp_print_key_section_create(print_key_section, location, name, description, print_level, each_iter_names, each_iter_values, add_last, filename, common_iter_levels, citations, unit_str)
creates a print_key section
utils to manipulate splines on the regular grid of a pw
integer, parameter, public pw_interp
integer, parameter, public spline3_nopbc_interp
integer, parameter, public spline3_pbc_interp
real(kind=dp) function, public cp_unit_to_cp2k(value, unit_str, defaults, power)
converts to the internal cp2k units to the given unit
Defines the basic variable types.
integer, parameter, public dp
This module defines the grid data type and some basic operations on it.
integer, parameter, public do_pw_grid_blocked_false
integer, parameter, public do_pw_grid_blocked_true
integer, parameter, public do_pw_grid_blocked_free
different utils that are useful to manipulate splines on the regular grid of a pw
integer, parameter, public precond_spl3_3
integer, parameter, public precond_spl3_aint
integer, parameter, public no_precond
integer, parameter, public precond_spl3_2
integer, parameter, public precond_spl3_aint2
integer, parameter, public precond_spl3_1
module that contains the definitions of the scf types
subroutine, public create_mixing_section(section, ls_scf)
Create CP2K input section for the mixing of the density matrix to be used only with diagonalization m...
manage control variables for the maximum overlap method
subroutine, public create_mom_section(section)
Create CP2K input section for variable occupancy using the Maximum Overlap Method....
Utilities for string manipulations.
character(len=1), parameter, public newline
1.9.8