39#include "./base/base_uses.f90"
49 REAL(kind=
dp) :: alpha = 0.0_dp
50 REAL(kind=
dp) :: mu = 0.0_dp
51 REAL(kind=
dp) :: t = 0.0_dp
52 REAL(kind=
dp) :: eps_hp = 0.0_dp
53 INTEGER :: natoms = 0, &
56 INTEGER,
DIMENSION(:),
POINTER :: atom_ids => null()
63 LOGICAL :: spherical = .false.
64 LOGICAL :: fullspace = .false.
65 INTEGER,
DIMENSION(2) :: distribution_layout = 0
66 INTEGER :: blocked = 0
67 END TYPE pw_grid_option
73 INTEGER,
DIMENSION(:),
ALLOCATABLE :: ref_mo_index
74 INTEGER :: ref_mo_spin = 1
75 INTEGER :: ref_nlumo = 0
76 LOGICAL :: sum_on_all_ref = .false.
77 INTEGER,
DIMENSION(:),
ALLOCATABLE :: td_mo_index
78 REAL(
dp),
DIMENSION(:),
ALLOCATABLE :: td_mo_occ
79 INTEGER :: td_mo_spin = 1
80 LOGICAL :: sum_on_all_td = .false.
81 CHARACTER(LEN=default_path_length) :: ref_mo_file_name =
""
82 LOGICAL :: propagate_ref = .false.
89 END TYPE proj_mo_p_type
95 LOGICAL :: converged = .false.
96 REAL(kind=
dp) :: eps_ener = 0.0_dp
97 INTEGER :: max_iter = 0
98 INTEGER :: mat_exp = 0
99 INTEGER :: propagator = 0
100 LOGICAL :: fixed_ions = .false.
103 INTEGER :: initial_wfn = 0
104 REAL(
dp) :: eps_exp = 0.0_dp
105 LOGICAL :: initial_step = .false.
106 LOGICAL :: hfx_redistribute = .false.
107 INTEGER :: aspc_order = 0
108 INTEGER :: sc_check_start = 0
109 LOGICAL :: apply_wfn_mix_init_restart = .false.
110 LOGICAL :: apply_delta_pulse = .false.
111 LOGICAL :: apply_delta_pulse_mag = .false.
112 LOGICAL :: periodic = .false.
113 LOGICAL :: linear_scaling = .false.
114 LOGICAL :: write_restart = .false.
115 INTEGER :: mcweeny_max_iter = 0
116 INTEGER :: acc_ref = 0
117 REAL(
dp) :: mcweeny_eps = 0.0_dp
118 INTEGER,
DIMENSION(3) :: delta_pulse_direction = 0
119 REAL(kind=
dp) :: delta_pulse_scale = 0.0_dp
120 LOGICAL :: velocity_gauge = .false.
121 REAL(kind=
dp),
DIMENSION(3) :: field = 0.0_dp
122 REAL(kind=
dp),
DIMENSION(3) :: vec_pot = 0.0_dp
123 LOGICAL :: nl_gauge_transform = .false.
124 LOGICAL :: is_proj_mo = .false.
125 TYPE(proj_mo_p_type),
DIMENSION(:), &
126 POINTER :: proj_mo_list => null()
128 LOGICAL :: save_local_moments = .false.
130 REAL(
dp),
DIMENSION(:),
POINTER :: moment_trace_user_ref_point => null()
131 REAL(
dp) :: ft_damping = -1.0_dp
132 REAL(
dp) :: ft_t0 = 0.0_dp
139 INTEGER,
DIMENSION(:, :),
POINTER :: print_pol_elements => null()
140 LOGICAL :: pade_requested = .false.
141 REAL(
dp) :: pade_e_min = 0.0_dp
142 REAL(
dp) :: pade_e_step = 0.02_dp
143 REAL(
dp) :: pade_e_max = 100.0_dp
144 REAL(
dp) :: pade_fit_e_min = 0.0_dp
145 REAL(
dp) :: pade_fit_e_max = 300.0_dp
152 LOGICAL :: self_consistent = .false.
153 LOGICAL :: orthogonal_basis = .false.
154 LOGICAL :: dispersion = .false.
156 INTEGER :: dispersion_type = 0
160 LOGICAL :: dftb3_diagonal = .false.
161 LOGICAL :: hb_sr_damp = .false.
162 REAL(kind=
dp) :: hb_sr_para = 0.0_dp
168 REAL(kind=
dp) :: eps_disp = 0.0_dp
169 REAL(kind=
dp) :: epscn = 0.0_dp
170 REAL(kind=
dp) :: exp_pre = 0.0_dp
171 REAL(kind=
dp) :: scaling = 0.0_dp
172 REAL(kind=
dp) :: rcdisp = 0.0_dp
173 REAL(kind=
dp),
DIMENSION(3) :: sd3 = 0.0_dp
174 REAL(kind=
dp),
DIMENSION(4) :: sd3bj = 0.0_dp
175 LOGICAL :: do_ewald = .false.
176 CHARACTER(LEN=default_path_length) :: sk_file_path =
""
177 CHARACTER(LEN=default_path_length) :: sk_file_list =
""
178 CHARACTER(LEN=default_string_length), &
179 DIMENSION(:, :),
POINTER :: sk_pair_list => null()
180 CHARACTER(LEN=default_path_length) :: uff_force_field =
""
181 CHARACTER(LEN=default_path_length) :: dispersion_parameter_file =
""
187 TYPE xtb_reference_cli_guess_type
188 LOGICAL :: enabled = .false.
189 LOGICAL :: efield_active = .false.
190 LOGICAL :: grad = .false.
193 REAL(kind=
dp),
DIMENSION(3) :: efield = 0.0_dp
194 REAL(kind=
dp) :: electronic_temperature_guess = 0.0_dp
195 CHARACTER(LEN=default_path_length) :: input_file =
""
196 CHARACTER(LEN=default_path_length) :: json_file =
""
197 CHARACTER(LEN=default_string_length) :: input_format =
"gen"
198 END TYPE xtb_reference_cli_guess_type
200 TYPE xtb_reference_cli_param_type
201 LOGICAL :: enabled = .false.
202 LOGICAL :: method_explicit = .false.
203 INTEGER :: method = 0
204 CHARACTER(LEN=default_path_length) :: input_file =
""
205 CHARACTER(LEN=default_path_length) :: output_file =
""
206 END TYPE xtb_reference_cli_param_type
208 TYPE xtb_reference_cli_fit_type
209 LOGICAL :: enabled = .false.
210 LOGICAL :: dry_run = .false.
211 CHARACTER(LEN=default_path_length) :: copy_file =
""
212 CHARACTER(LEN=default_path_length) :: input_file =
""
213 CHARACTER(LEN=default_path_length) :: param_file =
""
214 END TYPE xtb_reference_cli_fit_type
216 TYPE xtb_reference_cli_tagdiff_type
217 LOGICAL :: enabled = .false.
218 LOGICAL :: fit = .false.
219 CHARACTER(LEN=default_path_length) :: actual_file =
""
220 CHARACTER(LEN=default_path_length) :: reference_file =
""
221 END TYPE xtb_reference_cli_tagdiff_type
224 LOGICAL :: enabled = .false.
225 LOGICAL :: check_energy = .true.
226 LOGICAL :: check_forces = .true.
227 LOGICAL :: check_virial = .false.
228 LOGICAL :: keep_files = .false.
229 LOGICAL :: stop_on_error = .false.
230 LOGICAL :: efield_active = .false.
231 LOGICAL :: solvation_active = .false.
236 REAL(kind=
dp),
DIMENSION(3) :: efield = 0.0_dp
237 REAL(kind=
dp) :: electronic_temperature_guess = 0.0_dp
238 REAL(kind=
dp) :: error_limit = 1.0e-8_dp
239 CHARACTER(LEN=default_path_length) :: program_name =
"tblite"
240 CHARACTER(LEN=default_path_length) :: grad_file =
""
241 CHARACTER(LEN=default_path_length) :: json_file =
""
242 CHARACTER(LEN=default_path_length) :: post_processing =
""
243 CHARACTER(LEN=default_path_length) :: post_processing_output_file =
""
244 CHARACTER(LEN=default_path_length) :: restart_file =
""
245 CHARACTER(LEN=default_path_length) :: solvation_solvent =
""
246 CHARACTER(LEN=default_path_length) :: work_directory =
"."
247 CHARACTER(LEN=default_path_length) :: prefix =
"tblite-reference"
248 CHARACTER(LEN=default_string_length) :: input_format =
"gen"
249 TYPE(xtb_reference_cli_guess_type) :: guess_cli
250 TYPE(xtb_reference_cli_param_type) :: param_cli
251 TYPE(xtb_reference_cli_fit_type) :: fit_cli
252 TYPE(xtb_reference_cli_tagdiff_type) :: tagdiff_cli
257 INTEGER :: gfn_type = 1
259 LOGICAL :: do_ewald = .false.
260 LOGICAL :: do_tblite = .false.
261 LOGICAL :: tblite_mixer_damping_explicit = .false.
264 INTEGER :: h_sto_ng = 0
265 INTEGER :: tblite_method = 0
270 REAL(kind=
dp) :: tblite_accuracy = 1.0_dp
276 CHARACTER(LEN=default_path_length) :: tblite_param_file =
""
278 INTEGER :: vdw_type = -1
279 CHARACTER(LEN=default_path_length) :: parameter_file_path =
""
280 CHARACTER(LEN=default_path_length) :: parameter_file_name =
""
282 CHARACTER(LEN=default_path_length) :: dispersion_parameter_file =
""
283 REAL(kind=
dp) :: epscn = 0.0_dp
284 REAL(kind=
dp) :: rcdisp = 0.0_dp
285 REAL(kind=
dp) :: s6 = 0.0_dp, s8 = 0.0_dp
286 REAL(kind=
dp) :: a1 = 0.0_dp, a2 = 0.0_dp
288 REAL(kind=
dp) :: ks = 0.0_dp, kp = 0.0_dp, kd = 0.0_dp, ksp = 0.0_dp, k2sh = 0.0_dp
289 REAL(kind=
dp) :: kg = 0.0_dp, kf = 0.0_dp
290 REAL(kind=
dp) :: kcns = 0.0_dp, kcnp = 0.0_dp, kcnd = 0.0_dp
291 REAL(kind=
dp) :: ken = 0.0_dp
292 REAL(kind=
dp) :: ksen = 0.0_dp, kpen = 0.0_dp, kden = 0.0_dp
293 REAL(kind=
dp) :: ben = 0.0_dp
294 REAL(kind=
dp) :: kxr = 0.0_dp, kx2 = 0.0_dp
295 REAL(kind=
dp) :: enscale = 0.0_dp
297 LOGICAL :: xb_interaction = .false.
298 LOGICAL :: do_nonbonded = .false.
299 LOGICAL :: coulomb_interaction = .false.
300 LOGICAL :: coulomb_lr = .false.
301 LOGICAL :: tb3_interaction = .false.
302 LOGICAL :: check_atomic_charges = .false.
303 LOGICAL :: var_dipole = .false.
305 REAL(kind=
dp) :: xb_radius = 0.0_dp
306 REAL(kind=
dp) :: coulomb_sr_cut = 0.0_dp
307 REAL(kind=
dp) :: coulomb_sr_eps = 0.0_dp
309 CHARACTER(LEN=default_string_length), &
310 DIMENSION(:, :),
POINTER :: kab_param => null()
311 INTEGER,
DIMENSION(:, :),
POINTER :: kab_types => null()
312 INTEGER :: kab_nval = 0
313 REAL,
DIMENSION(:),
POINTER :: kab_vals => null()
316 REAL(kind=
dp) :: eps_pair = 0.0_dp
317 REAL(kind=
dp),
DIMENSION(:, :), &
318 POINTER :: rcpair => null()
321 REAL(kind=
dp) :: ksrb = 0.0_dp, esrb = 0.0_dp, gscal = 0.0_dp
322 REAL(kind=
dp) :: c1srb = 0.0_dp, c2srb = 0.0_dp, shift = 0.0_dp
325 INTEGER :: enshift_type = 1
334 LOGICAL :: orthogonal_basis = .false.
335 LOGICAL :: analytical_gradients = .false.
336 LOGICAL :: force_kdsod_ex = .false.
337 LOGICAL :: do_ewald = .false., do_ewald_r3 = .false., do_ewald_gks = .false.
338 INTEGER :: integral_screening = 0, periodic_type = 0
339 INTEGER :: max_multipole = 0
340 INTEGER :: ga_ncells = 0
341 REAL(kind=
dp) :: delta = 0.0_dp
343 LOGICAL :: dispersion = .false.
344 REAL(kind=
dp) :: rcdisp = 0.0_dp
345 REAL(kind=
dp) :: epscn = 0.0_dp
346 REAL(kind=
dp),
DIMENSION(3) :: sd3 = 0.0_dp
347 CHARACTER(LEN=default_path_length) :: dispersion_parameter_file =
""
349 REAL(kind=
dp) :: cutoff_lrc = 0.0_dp, taper_lrc = 0.0_dp, range_lrc = 0.0_dp
350 REAL(kind=
dp) :: cutoff_cou = 0.0_dp, taper_cou = 0.0_dp, range_cou = 0.0_dp
351 REAL(kind=
dp) :: cutoff_exc = 0.0_dp, taper_exc = 0.0_dp, range_exc = 0.0_dp
352 REAL(kind=
dp) :: taper_scr = 0.0_dp, range_scr = 0.0_dp
359 INTEGER :: basis_1c = 0
360 REAL(kind=
dp) :: eps_fit = 0.0_dp, &
362 eps_vrho0 = 0.0_dp, &
365 INTEGER :: ladd_rho0 = 0, &
369 LOGICAL :: accurate_xcint = .false.
370 REAL(kind=
dp) :: aweights = 0.0_dp
371 REAL(kind=
dp),
DIMENSION(:),
POINTER :: aw => null()
372 LOGICAL :: alpha0_hard_from_input = .false., &
373 force_paw = .false., &
374 non_paw_atoms = .false., &
375 nopaw_as_gpw = .false.
376 REAL(kind=
dp) :: alpha0_hard = 0.0_dp
377 REAL(kind=
dp) :: max_rad_local = 0.0_dp
384 REAL(kind=
dp) :: actual_time = 0.0_dp
385 REAL(kind=
dp),
DIMENSION(:),
POINTER :: polarisation => null()
386 INTEGER :: envelop_id = 0
387 REAL(kind=
dp),
DIMENSION(:),
POINTER :: envelop_r_vars => null()
388 INTEGER,
DIMENSION(:),
POINTER :: envelop_i_vars => null()
389 REAL(kind=
dp) :: strength = 0.0_dp
390 REAL(kind=
dp) :: amplitude = 0.0_dp
391 REAL(kind=
dp) :: phase_offset = 0.0_dp
392 REAL(kind=
dp) :: wavelength = 0.0_dp
393 REAL(kind=
dp),
DIMENSION(3) :: vec_pot_initial = 0.0_dp
398 END TYPE efield_p_type
403 TYPE period_efield_type
404 LOGICAL :: displacement_field = .false.
405 REAL(kind=
dp),
DIMENSION(3) :: polarisation = 0.0_dp
406 REAL(kind=
dp),
DIMENSION(3) :: d_filter = 0.0_dp
407 REAL(kind=
dp) :: strength = 0.0_dp
408 REAL(kind=
dp),
ALLOCATABLE,
DIMENSION(:) :: strength_list
409 INTEGER :: start_frame = 0
410 INTEGER :: end_frame = -1
411 END TYPE period_efield_type
417 REAL(kind=
dp) :: strength = 0.0_dp
418 REAL(kind=
dp) ::
TARGET = 0.0_dp
419 INTEGER :: natoms = 0
420 INTEGER,
POINTER,
DIMENSION(:) :: atoms => null()
427 INTEGER :: ref_count = 0
428 REAL(kind=
dp) :: strength = 0.0_dp
429 REAL(kind=
dp) ::
TARGET = 0.0_dp
430 REAL(kind=
dp) :: ddapc_order_p = 0.0_dp
431 INTEGER :: functional_form = 0
432 INTEGER :: natoms = 0
433 INTEGER,
POINTER,
DIMENSION(:) :: atoms => null()
434 REAL(kind=
dp),
POINTER,
DIMENSION(:) :: coeff => null()
435 INTEGER :: density_type = 0
442 REAL(kind=
dp) :: strength = 0.0_dp
443 REAL(kind=
dp) ::
TARGET = 0.0_dp
444 REAL(kind=
dp) :: s2_order_p = 0.0_dp
445 INTEGER :: functional_form = 0
452 INTEGER,
DIMENSION(:),
ALLOCATABLE ::
list
453 END TYPE admm_block_type
456 REAL(kind=
dp) :: eps_filter = 0.0_dp
457 INTEGER :: admm_type = 0
458 INTEGER :: purification_method = 0
459 INTEGER :: method = 0
460 LOGICAL :: charge_constrain = .false.
461 INTEGER :: scaling_model = 0
462 INTEGER :: aux_exch_func = 0
463 LOGICAL :: aux_exch_func_param = .false.
464 REAL(kind=
dp),
DIMENSION(3) :: aux_x_param = 0.0_dp
465 TYPE(admm_block_type),
DIMENSION(:), &
466 ALLOCATABLE :: blocks
473 LOGICAL :: read_from_cube = .false.
474 LOGICAL :: maxwell_solver = .false.
475 LOGICAL :: static = .false.
476 REAL(kind=
dp) :: scaling_factor = 0.0_dp
483 LOGICAL :: log_test = .false.
484 INTEGER :: int_test = 0
485 REAL(kind=
dp) :: real_test = 0.0_dp
494 INTEGER :: method_id = 0
495 REAL(kind=
dp) :: eps_core_charge = 0.0_dp, &
496 eps_kg_orb = 0.0_dp, &
497 eps_pgf_orb = 0.0_dp, &
500 eps_rho_gspace = 0.0_dp, &
501 eps_rho_rspace = 0.0_dp, &
502 eps_filter_matrix = 0.0_dp, &
503 eps_gvg_rspace = 0.0_dp, &
504 progression_factor = 0.0_dp, &
505 relative_cutoff = 0.0_dp
506 LOGICAL :: do_almo_scf = .false.
507 LOGICAL :: do_ls_scf = .false.
508 LOGICAL :: do_kg = .false.
509 LOGICAL :: commensurate_mgrids = .false.
510 LOGICAL :: realspace_mgrids = .false.
511 LOGICAL :: gapw = .false., gapw_xc = .false., gpw = .false., pao = .false.
512 LOGICAL :: lrigpw = .false., rigpw = .false.
513 LOGICAL :: lri_optbas = .false.
514 LOGICAL :: ofgpw = .false.
515 LOGICAL :: dftb = .false.
516 LOGICAL :: xtb = .false.
517 LOGICAL :: semi_empirical = .false.
518 LOGICAL :: mulliken_restraint = .false.
519 LOGICAL :: ddapc_restraint = .false.
520 LOGICAL :: ddapc_restraint_is_spin = .false.
521 LOGICAL :: ddapc_explicit_potential = .false.
522 LOGICAL :: cdft = .false.
523 LOGICAL :: et_coupling_calc = .false.
524 LOGICAL :: s2_restraint = .false.
525 INTEGER :: do_ppl_method = 0
526 INTEGER :: wf_interpolation_method_nr = 0
527 INTEGER :: wf_extrapolation_order = 0
528 INTEGER :: periodicity = 0
529 REAL(kind=
dp) :: pairlist_radius = 0.0_dp
530 REAL(kind=
dp) :: cutoff = 0.0_dp
531 REAL(kind=
dp),
DIMENSION(:),
POINTER :: e_cutoff => null()
533 POINTER :: mulliken_restraint_control => null()
535 DIMENSION(:),
POINTER :: ddapc_restraint_control => null()
541 POINTER :: se_control => null()
543 TYPE(pw_grid_option) :: pw_grid_opt = pw_grid_option()
544 LOGICAL :: skip_load_balance_distributed = .false.
546 LOGICAL :: ref_embed_subsys = .false.
547 LOGICAL :: cluster_embed_subsys = .false.
548 LOGICAL :: high_level_embed_subsys = .false.
549 LOGICAL :: dfet_embedded = .false.
550 LOGICAL :: dmfet_embedded = .false.
557 LOGICAL :: sccs_activated = .false.
558 INTEGER :: derivative_method = 0, &
561 REAL(kind=
dp) :: alpha_solvent = 0.0_dp, &
563 beta_solvent = 0.0_dp, &
564 delta_rho = 0.0_dp, &
567 epsilon_solvent = 0.0_dp, &
568 gamma_solvent = 0.0_dp, &
581 LOGICAL :: do_ewald = .false.
582 LOGICAL :: do_exchange = .false.
583 REAL(kind=
dp) :: hfx_fraction = 0.0_dp
584 REAL(kind=
dp) :: eps_td_filter = 0.0_dp
585 REAL(kind=
dp) :: mn_alpha = 0.0_dp
586 REAL(kind=
dp) :: mn_beta = 0.0_dp
587 REAL(kind=
dp) :: coulomb_sr_cut = 0.0_dp
588 REAL(kind=
dp) :: coulomb_sr_eps = 0.0_dp
597 REAL(kind=
dp),
DIMENSION(:),
POINTER :: fermia => null()
598 REAL(kind=
dp),
DIMENSION(:, :),
POINTER :: fermib => null()
605 INTEGER :: surf_normal = 0
612 INTEGER :: surf_normal = 0
613 REAL(kind=
dp) :: dist_edge = 0.0_dp, &
623 LOGICAL :: enabled = .false.
625 INTEGER :: nstates = 0
627 INTEGER :: niters = 0
633 INTEGER :: nprocs = 0
635 INTEGER :: kernel = 0
637 INTEGER :: spinflip = 0
639 LOGICAL :: do_hfx = .false.
640 LOGICAL :: do_admm = .false.
642 LOGICAL :: do_hfxsr = .false.
643 LOGICAL :: hfxsr_re_int = .true.
644 INTEGER :: hfxsr_primbas = 0
645 LOGICAL :: do_hfxlr = .false.
646 REAL(kind=
dp) :: hfxlr_rcut = 0.0_dp, hfxlr_scale = 0.0_dp
647 LOGICAL :: do_exck = .false.
651 INTEGER :: oe_corr = 0
653 REAL(kind=
dp) :: ev_shift = 0.0_dp, eos_shift = 0.0_dp
655 INTEGER,
DIMENSION(2) :: nactive = -1
657 REAL(kind=
dp) :: conv = 0.0_dp
659 REAL(kind=
dp) :: min_excitation_amplitude = 0.0_dp
662 REAL(kind=
dp) :: orthogonal_eps = 0.0_dp
664 LOGICAL :: is_restart = .false.
666 LOGICAL :: rks_triplets = .false.
668 LOGICAL :: do_lrigpw = .false.
670 LOGICAL :: do_smearing = .false.
672 LOGICAL :: do_bse = .false.
673 LOGICAL :: do_bse_w_only = .false.
674 LOGICAL :: do_bse_gw_only = .false.
676 INTEGER :: auto_basis_p_lri_aux = 1
678 LOGICAL :: admm_symm = .false.
680 LOGICAL :: admm_xc_correction = .false.
682 LOGICAL :: do_exciton_descriptors = .false.
683 LOGICAL :: do_directional_exciton_descriptors = .false.
688 INTEGER :: dipole_form = 0
690 INTEGER :: dipole_reference = 0
692 REAL(kind=
dp),
DIMENSION(:),
POINTER :: dipole_ref_point => null()
695 LOGICAL :: do_soc = .false.
700 INTEGER :: mgrid_ngrids = 0
702 LOGICAL :: mgrid_commensurate_mgrids = .false.
706 LOGICAL :: mgrid_is_explicit = .false.
708 LOGICAL :: mgrid_realspace_mgrids = .false.
710 LOGICAL :: mgrid_skip_load_balance = .false.
712 REAL(kind=
dp) :: mgrid_cutoff = 0.0_dp
715 REAL(kind=
dp) :: mgrid_progression_factor = 0.0_dp
717 REAL(kind=
dp) :: mgrid_relative_cutoff = 0.0_dp
720 REAL(kind=
dp),
DIMENSION(:),
POINTER :: mgrid_e_cutoff => null()
730 LOGICAL :: enabled = .false.
731 INTEGER :: core_states = 0
732 INTEGER :: valence_states = 0
746 TYPE(period_efield_type),
POINTER :: period_efield => null()
758 TYPE(efield_p_type),
POINTER, &
759 DIMENSION(:) :: efield_fields => null()
761 DIMENSION(:) :: probe => null()
762 INTEGER :: nspins = 0, &
766 plus_u_method_id = 0, &
769 INTEGER :: sic_list_id = 0
770 INTEGER :: auto_basis_ri_aux = 1, &
771 auto_basis_aux_fit = 1, &
772 auto_basis_lri_aux = 1, &
773 auto_basis_p_lri_aux = 1, &
774 auto_basis_ri_hxc = 1, &
775 auto_basis_ri_xas = 1, &
776 auto_basis_ri_hfx = 1, &
778 REAL(kind=
dp) :: relax_multiplicity = 0.0_dp, &
779 sic_scaling_a = 0.0_dp, &
780 sic_scaling_b = 0.0_dp, &
781 pos_dir_surf_dip = 0.0_dp
782 LOGICAL :: do_xas_calculation = .false., &
783 do_xas_tdp_calculation = .false., &
784 drho_by_collocation = .false., &
785 use_kinetic_energy_density = .false., &
786 restricted = .false., &
791 apply_efield = .false., &
792 apply_efield_field = .false., &
793 apply_vector_potential = .false., &
794 apply_period_efield = .false., &
795 apply_external_potential = .false., &
796 eval_external_potential = .false., &
798 do_admm_dm = .false., &
799 do_admm_mo = .false., &
801 low_spin_roks = .false., &
802 apply_external_density = .false., &
803 read_external_density = .false., &
804 apply_external_vxc = .false., &
805 read_external_vxc = .false., &
806 correct_surf_dip = .false., &
807 surf_dip_correct_switch = .false., &
808 switch_surf_dip = .false., &
809 correct_el_density_dip = .false., &
811 apply_embed_pot = .false., &
812 apply_dmfet_pot = .false., &
820 CHARACTER(len=*),
PARAMETER,
PRIVATE :: modulen =
'cp_control_types'
869 SUBROUTINE mulliken_control_create(mulliken_restraint_control)
872 mulliken_restraint_control%strength = 0.1_dp
873 mulliken_restraint_control%target = 1.0_dp
874 mulliken_restraint_control%natoms = 0
875 NULLIFY (mulliken_restraint_control%atoms)
876 END SUBROUTINE mulliken_control_create
884 SUBROUTINE mulliken_control_release(mulliken_restraint_control)
887 IF (
ASSOCIATED(mulliken_restraint_control%atoms)) &
888 DEALLOCATE (mulliken_restraint_control%atoms)
889 mulliken_restraint_control%strength = 0.0_dp
890 mulliken_restraint_control%target = 0.0_dp
891 mulliken_restraint_control%natoms = 0
892 END SUBROUTINE mulliken_control_release
904 ddapc_restraint_control%strength = 0.1_dp
905 ddapc_restraint_control%ddapc_order_p = 0.0_dp
906 ddapc_restraint_control%functional_form = -1
907 ddapc_restraint_control%target = 1.0_dp
908 ddapc_restraint_control%natoms = 0
909 NULLIFY (ddapc_restraint_control%atoms)
910 NULLIFY (ddapc_restraint_control%coeff)
920 SUBROUTINE ddapc_control_release(ddapc_restraint_control)
923 IF (
ASSOCIATED(ddapc_restraint_control%atoms)) &
924 DEALLOCATE (ddapc_restraint_control%atoms)
925 IF (
ASSOCIATED(ddapc_restraint_control%coeff)) &
926 DEALLOCATE (ddapc_restraint_control%coeff)
927 ddapc_restraint_control%strength = 0.0_dp
928 ddapc_restraint_control%target = 0.0_dp
929 ddapc_restraint_control%natoms = 0
930 END SUBROUTINE ddapc_control_release
938 SUBROUTINE s2_control_create(s2_restraint_control)
941 s2_restraint_control%strength = 0.1_dp
942 s2_restraint_control%s2_order_p = 0.0_dp
943 s2_restraint_control%functional_form = -1
944 s2_restraint_control%target = 1.0_dp
945 END SUBROUTINE s2_control_create
953 SUBROUTINE s2_control_release(s2_restraint_control)
956 s2_restraint_control%strength = 0.0_dp
957 s2_restraint_control%target = 0.0_dp
958 END SUBROUTINE s2_control_release
970 NULLIFY (dft_control%xas_control)
971 NULLIFY (dft_control%qs_control)
972 NULLIFY (dft_control%tddfpt2_control)
973 NULLIFY (dft_control%rixs_control)
974 NULLIFY (dft_control%efield_fields)
975 NULLIFY (dft_control%period_efield)
976 NULLIFY (dft_control%admm_control)
977 NULLIFY (dft_control%expot_control)
978 NULLIFY (dft_control%maxwell_control)
979 NULLIFY (dft_control%smeagol_control)
980 NULLIFY (dft_control%rtp_control)
981 NULLIFY (dft_control%sccs_control)
982 NULLIFY (dft_control%probe)
983 NULLIFY (dft_control%paep_control)
984 NULLIFY (dft_control%pcc_control)
985 dft_control%do_sccs = .false.
986 dft_control%apply_embed_pot = .false.
987 dft_control%apply_dmfet_pot = .false.
988 dft_control%hairy_probes = .false.
989 dft_control%do_pcc = .false.
990 dft_control%do_paep = .false.
991 CALL qs_control_create(dft_control%qs_control)
992 CALL tddfpt2_control_create(dft_control%tddfpt2_control)
1009 CALL qs_control_release(dft_control%qs_control)
1010 CALL tddfpt2_control_release(dft_control%tddfpt2_control)
1012 IF (
ASSOCIATED(dft_control%xas_control))
THEN
1014 DEALLOCATE (dft_control%xas_control)
1017 CALL expot_control_release(dft_control%expot_control)
1018 CALL maxwell_control_release(dft_control%maxwell_control)
1020 CALL efield_fields_release(dft_control%efield_fields)
1021 IF (
ASSOCIATED(dft_control%probe))
THEN
1022 DO i = 1,
SIZE(dft_control%probe)
1023 DEALLOCATE (dft_control%probe(i)%atom_ids)
1025 DEALLOCATE (dft_control%probe)
1027 IF (
ASSOCIATED(dft_control%sccs_control))
DEALLOCATE (dft_control%sccs_control)
1028 IF (
ASSOCIATED(dft_control%period_efield))
THEN
1029 DEALLOCATE (dft_control%period_efield)
1031 IF (
ASSOCIATED(dft_control%rtp_control))
THEN
1032 CALL proj_mo_list_release(dft_control%rtp_control%proj_mo_list)
1033 DEALLOCATE (dft_control%rtp_control)
1035 IF (
ASSOCIATED(dft_control%pcc_control))
DEALLOCATE (dft_control%pcc_control)
1036 IF (
ASSOCIATED(dft_control%paep_control))
DEALLOCATE (dft_control%paep_control)
1044 SUBROUTINE qs_control_create(qs_control)
1047 cpassert(.NOT.
ASSOCIATED(qs_control))
1048 ALLOCATE (qs_control)
1050 NULLIFY (qs_control%e_cutoff)
1051 NULLIFY (qs_control%gapw_control)
1052 NULLIFY (qs_control%mulliken_restraint_control)
1053 NULLIFY (qs_control%ddapc_restraint_control)
1054 NULLIFY (qs_control%s2_restraint_control)
1055 NULLIFY (qs_control%se_control)
1056 NULLIFY (qs_control%dftb_control)
1057 NULLIFY (qs_control%xtb_control)
1058 NULLIFY (qs_control%cdft_control)
1059 NULLIFY (qs_control%ddapc_restraint_control)
1061 ALLOCATE (qs_control%mulliken_restraint_control)
1062 CALL mulliken_control_create(qs_control%mulliken_restraint_control)
1063 ALLOCATE (qs_control%s2_restraint_control)
1064 CALL s2_control_create(qs_control%s2_restraint_control)
1065 ALLOCATE (qs_control%gapw_control)
1066 CALL se_control_create(qs_control%se_control)
1067 CALL dftb_control_create(qs_control%dftb_control)
1068 CALL xtb_control_create(qs_control%xtb_control)
1069 ALLOCATE (qs_control%cdft_control)
1071 END SUBROUTINE qs_control_create
1077 SUBROUTINE qs_control_release(qs_control)
1082 IF (
ASSOCIATED(qs_control))
THEN
1083 CALL mulliken_control_release(qs_control%mulliken_restraint_control)
1084 DEALLOCATE (qs_control%mulliken_restraint_control)
1085 CALL s2_control_release(qs_control%s2_restraint_control)
1086 DEALLOCATE (qs_control%s2_restraint_control)
1087 CALL se_control_release(qs_control%se_control)
1088 CALL dftb_control_release(qs_control%dftb_control)
1089 CALL xtb_control_release(qs_control%xtb_control)
1090 IF (
ASSOCIATED(qs_control%cdft_control))
THEN
1092 DEALLOCATE (qs_control%cdft_control)
1095 IF (
ASSOCIATED(qs_control%e_cutoff))
THEN
1096 DEALLOCATE (qs_control%e_cutoff)
1098 IF (
ASSOCIATED(qs_control%gapw_control))
THEN
1099 IF (
ASSOCIATED(qs_control%gapw_control%aw))
THEN
1100 DEALLOCATE (qs_control%gapw_control%aw)
1102 DEALLOCATE (qs_control%gapw_control)
1104 IF (
ASSOCIATED(qs_control%ddapc_restraint_control))
THEN
1105 DO i = 1,
SIZE(qs_control%ddapc_restraint_control)
1106 CALL ddapc_control_release(qs_control%ddapc_restraint_control(i))
1108 DEALLOCATE (qs_control%ddapc_restraint_control)
1110 DEALLOCATE (qs_control)
1112 END SUBROUTINE qs_control_release
1120 SUBROUTINE tddfpt2_control_create(tddfpt_control)
1123 CHARACTER(len=*),
PARAMETER :: routinen =
'tddfpt2_control_create'
1127 cpassert(.NOT.
ASSOCIATED(tddfpt_control))
1128 CALL timeset(routinen, handle)
1130 ALLOCATE (tddfpt_control)
1131 tddfpt_control%do_soc = .false.
1133 CALL timestop(handle)
1134 END SUBROUTINE tddfpt2_control_create
1142 SUBROUTINE tddfpt2_control_release(tddfpt_control)
1145 CHARACTER(len=*),
PARAMETER :: routinen =
'tddfpt2_control_release'
1149 CALL timeset(routinen, handle)
1151 IF (
ASSOCIATED(tddfpt_control))
THEN
1152 DEALLOCATE (tddfpt_control)
1155 CALL timestop(handle)
1156 END SUBROUTINE tddfpt2_control_release
1165 cpassert(.NOT.
ASSOCIATED(rixs_control))
1166 ALLOCATE (rixs_control)
1168 NULLIFY (rixs_control%tddfpt2_control)
1169 CALL tddfpt2_control_create(rixs_control%tddfpt2_control)
1170 NULLIFY (rixs_control%xas_tdp_control)
1182 IF (
ASSOCIATED(rixs_control))
THEN
1183 CALL tddfpt2_control_release(rixs_control%tddfpt2_control)
1186 DEALLOCATE (rixs_control)
1195 SUBROUTINE proj_mo_list_release(proj_mo_list)
1196 TYPE(proj_mo_p_type),
DIMENSION(:),
POINTER :: proj_mo_list
1198 INTEGER :: i, mo_ref_nbr
1200 IF (
ASSOCIATED(proj_mo_list))
THEN
1201 DO i = 1,
SIZE(proj_mo_list)
1202 IF (
ASSOCIATED(proj_mo_list(i)%proj_mo))
THEN
1203 IF (
ALLOCATED(proj_mo_list(i)%proj_mo%ref_mo_index)) &
1204 DEALLOCATE (proj_mo_list(i)%proj_mo%ref_mo_index)
1205 IF (
ALLOCATED(proj_mo_list(i)%proj_mo%mo_ref))
THEN
1206 DO mo_ref_nbr = 1,
SIZE(proj_mo_list(i)%proj_mo%mo_ref)
1207 CALL cp_fm_release(proj_mo_list(i)%proj_mo%mo_ref(mo_ref_nbr))
1209 DEALLOCATE (proj_mo_list(i)%proj_mo%mo_ref)
1211 IF (
ALLOCATED(proj_mo_list(i)%proj_mo%td_mo_index)) &
1212 DEALLOCATE (proj_mo_list(i)%proj_mo%td_mo_index)
1213 IF (
ALLOCATED(proj_mo_list(i)%proj_mo%td_mo_occ)) &
1214 DEALLOCATE (proj_mo_list(i)%proj_mo%td_mo_occ)
1215 DEALLOCATE (proj_mo_list(i)%proj_mo)
1218 DEALLOCATE (proj_mo_list)
1220 END SUBROUTINE proj_mo_list_release
1226 SUBROUTINE efield_fields_release(efield_fields)
1227 TYPE(efield_p_type),
DIMENSION(:),
POINTER :: efield_fields
1231 IF (
ASSOCIATED(efield_fields))
THEN
1232 DO i = 1,
SIZE(efield_fields)
1233 IF (
ASSOCIATED(efield_fields(i)%efield))
THEN
1234 IF (
ASSOCIATED(efield_fields(i)%efield%envelop_r_vars))
THEN
1235 DEALLOCATE (efield_fields(i)%efield%envelop_r_vars)
1237 IF (
ASSOCIATED(efield_fields(i)%efield%envelop_i_vars))
THEN
1238 DEALLOCATE (efield_fields(i)%efield%envelop_i_vars)
1240 IF (
ASSOCIATED(efield_fields(i)%efield%polarisation)) &
1241 DEALLOCATE (efield_fields(i)%efield%polarisation)
1242 DEALLOCATE (efield_fields(i)%efield)
1245 DEALLOCATE (efield_fields)
1247 END SUBROUTINE efield_fields_release
1253 SUBROUTINE dftb_control_create(dftb_control)
1256 cpassert(.NOT.
ASSOCIATED(dftb_control))
1257 ALLOCATE (dftb_control)
1259 NULLIFY (dftb_control%sk_pair_list)
1260 END SUBROUTINE dftb_control_create
1266 SUBROUTINE dftb_control_release(dftb_control)
1269 IF (
ASSOCIATED(dftb_control))
THEN
1270 IF (
ASSOCIATED(dftb_control%sk_pair_list))
THEN
1271 DEALLOCATE (dftb_control%sk_pair_list)
1273 DEALLOCATE (dftb_control)
1275 END SUBROUTINE dftb_control_release
1281 SUBROUTINE xtb_control_create(xtb_control)
1284 cpassert(.NOT.
ASSOCIATED(xtb_control))
1285 ALLOCATE (xtb_control)
1287 NULLIFY (xtb_control%kab_param)
1288 NULLIFY (xtb_control%kab_vals)
1289 NULLIFY (xtb_control%kab_types)
1290 NULLIFY (xtb_control%nonbonded)
1291 NULLIFY (xtb_control%rcpair)
1293 END SUBROUTINE xtb_control_create
1299 SUBROUTINE xtb_control_release(xtb_control)
1302 IF (
ASSOCIATED(xtb_control))
THEN
1303 IF (
ASSOCIATED(xtb_control%kab_param))
THEN
1304 DEALLOCATE (xtb_control%kab_param)
1306 IF (
ASSOCIATED(xtb_control%kab_vals))
THEN
1307 DEALLOCATE (xtb_control%kab_vals)
1309 IF (
ASSOCIATED(xtb_control%kab_types))
THEN
1310 DEALLOCATE (xtb_control%kab_types)
1312 IF (
ASSOCIATED(xtb_control%rcpair))
THEN
1313 DEALLOCATE (xtb_control%rcpair)
1315 IF (
ASSOCIATED(xtb_control%nonbonded))
THEN
1318 DEALLOCATE (xtb_control)
1320 END SUBROUTINE xtb_control_release
1326 SUBROUTINE se_control_create(se_control)
1329 cpassert(.NOT.
ASSOCIATED(se_control))
1330 ALLOCATE (se_control)
1331 END SUBROUTINE se_control_create
1337 SUBROUTINE se_control_release(se_control)
1340 IF (
ASSOCIATED(se_control))
THEN
1341 DEALLOCATE (se_control)
1343 END SUBROUTINE se_control_release
1352 cpassert(.NOT.
ASSOCIATED(admm_control))
1353 ALLOCATE (admm_control)
1364 IF (
ASSOCIATED(admm_control))
THEN
1365 DEALLOCATE (admm_control)
1376 cpassert(.NOT.
ASSOCIATED(expot_control))
1377 ALLOCATE (expot_control)
1378 expot_control%read_from_cube = .false.
1379 expot_control%maxwell_solver = .false.
1380 expot_control%static = .true.
1381 expot_control%scaling_factor = 1.0_dp
1389 SUBROUTINE expot_control_release(expot_control)
1392 IF (
ASSOCIATED(expot_control))
THEN
1393 DEALLOCATE (expot_control)
1396 END SUBROUTINE expot_control_release
1405 cpassert(.NOT.
ASSOCIATED(maxwell_control))
1406 ALLOCATE (maxwell_control)
1414 SUBROUTINE maxwell_control_release(maxwell_control)
1417 IF (
ASSOCIATED(maxwell_control))
THEN
1418 DEALLOCATE (maxwell_control)
1421 END SUBROUTINE maxwell_control_release
Defines control structures, which contain the parameters and the settings for the DFT-based calculati...
subroutine, public rixs_control_release(rixs_control)
Releases the rixs_control_type.
subroutine, public admm_control_release(admm_control)
...
subroutine, public dft_control_create(dft_control)
allocates and perform a very basic initialization
subroutine, public expot_control_create(expot_control)
...
subroutine, public maxwell_control_create(maxwell_control)
...
subroutine, public ddapc_control_create(ddapc_restraint_control)
create the ddapc_restraint_type
subroutine, public admm_control_create(admm_control)
...
subroutine, public dft_control_release(dft_control)
...
subroutine, public rixs_control_create(rixs_control)
Creates and initializes the rixs_control_type.
represent a full matrix distributed on many processors
Add the DFT+U contribution to the Hamiltonian matrix.
Defines the basic variable types.
integer, parameter, public dp
integer, parameter, public default_string_length
integer, parameter, public default_path_length
An array-based list which grows on demand. When the internal array is full, a new array of twice the ...
subroutine, public pair_potential_p_release(potparm)
Release Data-structure that constains potential parameters.
Defines CDFT control structures.
subroutine, public cdft_control_release(cdft_control)
release the cdft_control_type
subroutine, public cdft_control_create(cdft_control)
create the cdft_control_type
Input control types for NEGF/SMEAGOL transport calculations.
subroutine, public smeagol_control_release(smeagol_control)
release SMEAGOL control object
subroutine, public smeagol_control_create(smeagol_control)
allocate control options for SMEAGOL calculation
CP2K-side tblite-compatible SCC Broyden mixer.
Defines control structures, which contain the parameters and the settings for the calculations.
subroutine, public xas_control_release(xas_control)
...
Define XAS TDP control type and associated create, release, etc subroutines, as well as XAS TDP envir...
subroutine, public xas_tdp_control_release(xas_tdp_control)
Releases the xas_tdp_control_type.
subroutine, public xas_tdp_control_create(xas_tdp_control)
Creates and initializes the xas_tdp_control_type.
SMEAGOL-related input parameters.
A type that holds controlling information for a xas calculation.
Type containing control information for TDP XAS calculations.