(git:d18deda)
Loading...
Searching...
No Matches
qs_force.F
Go to the documentation of this file.
1!--------------------------------------------------------------------------------------------------!
2! CP2K: A general program to perform molecular dynamics simulations !
3! Copyright 2000-2025 CP2K developers group <https://cp2k.org> !
4! !
5! SPDX-License-Identifier: GPL-2.0-or-later !
6!--------------------------------------------------------------------------------------------------!
7
8! **************************************************************************************************
9!> \brief Quickstep force driver routine
10!> \author MK (12.06.2002)
11! **************************************************************************************************
12MODULE qs_force
13 USE atomic_kind_types, ONLY: atomic_kind_type,&
14 get_atomic_kind_set
15 USE cp_control_types, ONLY: dft_control_type
16 USE cp_dbcsr_api, ONLY: dbcsr_copy,&
17 dbcsr_p_type,&
18 dbcsr_set
19 USE cp_dbcsr_operations, ONLY: dbcsr_allocate_matrix_set,&
20 dbcsr_deallocate_matrix_set
21 USE cp_dbcsr_output, ONLY: cp_dbcsr_write_sparse_matrix
22 USE cp_log_handling, ONLY: cp_get_default_logger,&
23 cp_logger_get_default_io_unit,&
24 cp_logger_type
25 USE cp_output_handling, ONLY: cp_p_file,&
26 cp_print_key_finished_output,&
27 cp_print_key_should_output,&
28 cp_print_key_unit_nr
29 USE dft_plus_u, ONLY: plus_u
30 USE ec_env_types, ONLY: energy_correction_type
31 USE efield_utils, ONLY: calculate_ecore_efield,&
32 efield_potential_lengh_gauge
33 USE energy_corrections, ONLY: energy_correction
34 USE excited_states, ONLY: excited_state_energy
35 USE hfx_exx, ONLY: calculate_exx
36 USE input_constants, ONLY: ri_mp2_laplace,&
37 ri_mp2_method_gpw,&
38 ri_rpa_method_gpw
39 USE input_section_types, ONLY: section_vals_get,&
40 section_vals_get_subs_vals,&
41 section_vals_type,&
42 section_vals_val_get
43 USE kinds, ONLY: dp
44 USE lri_environment_types, ONLY: lri_environment_type
45 USE message_passing, ONLY: mp_para_env_type
46 USE mp2_cphf, ONLY: update_mp2_forces
47 USE mulliken, ONLY: mulliken_restraint
48 USE particle_types, ONLY: particle_type
49 USE qs_core_energies, ONLY: calculate_ecore_overlap,&
50 calculate_ecore_self
51 USE qs_core_hamiltonian, ONLY: build_core_hamiltonian_matrix
52 USE qs_dftb_dispersion, ONLY: calculate_dftb_dispersion
53 USE qs_dftb_matrices, ONLY: build_dftb_matrices
54 USE qs_energy, ONLY: qs_energies
55 USE qs_energy_types, ONLY: qs_energy_type
56 USE qs_environment_methods, ONLY: qs_env_rebuild_pw_env
57 USE qs_environment_types, ONLY: get_qs_env,&
58 qs_environment_type
59 USE qs_external_potential, ONLY: external_c_potential,&
60 external_e_potential
61 USE qs_force_types, ONLY: allocate_qs_force,&
62 qs_force_type,&
63 replicate_qs_force,&
64 zero_qs_force
65 USE qs_ks_methods, ONLY: qs_ks_update_qs_env
66 USE qs_ks_types, ONLY: qs_ks_env_type,&
67 set_ks_env
68 USE qs_rho_types, ONLY: qs_rho_get,&
69 qs_rho_type
70 USE qs_scf_post_scf, ONLY: qs_scf_compute_properties
71 USE qs_subsys_types, ONLY: qs_subsys_set,&
72 qs_subsys_type
73 USE ri_environment_methods, ONLY: build_ri_matrices
74 USE rt_propagation_forces, ONLY: calc_c_mat_force,&
75 rt_admm_force
76 USE rt_propagation_velocity_gauge, ONLY: velocity_gauge_ks_matrix,&
77 velocity_gauge_nl_force
78 USE se_core_core, ONLY: se_core_core_interaction
79 USE se_core_matrix, ONLY: build_se_core_matrix
80 USE virial_types, ONLY: symmetrize_virial,&
81 virial_type
82 USE xtb_matrices, ONLY: build_xtb_matrices
83#include "./base/base_uses.f90"
84
85 IMPLICIT NONE
86
87 PRIVATE
88
89! *** Global parameters ***
90
91 CHARACTER(len=*), PARAMETER, PRIVATE :: moduleN = 'qs_force'
92
93! *** Public subroutines ***
94
95 PUBLIC :: qs_calc_energy_force
96
97CONTAINS
98
99! **************************************************************************************************
100!> \brief ...
101!> \param qs_env ...
102!> \param calc_force ...
103!> \param consistent_energies ...
104!> \param linres ...
105! **************************************************************************************************
106 SUBROUTINE qs_calc_energy_force(qs_env, calc_force, consistent_energies, linres)
107 TYPE(qs_environment_type), POINTER :: qs_env
108 LOGICAL :: calc_force, consistent_energies, linres
109
110 qs_env%linres_run = linres
111 IF (calc_force) THEN
112 CALL qs_forces(qs_env)
113 ELSE
114 CALL qs_energies(qs_env, calc_forces=.false., &
115 consistent_energies=consistent_energies)
116 END IF
117
118 END SUBROUTINE qs_calc_energy_force
119
120! **************************************************************************************************
121!> \brief Calculate the Quickstep forces.
122!> \param qs_env ...
123!> \date 29.10.2002
124!> \author MK
125!> \version 1.0
126! **************************************************************************************************
127 SUBROUTINE qs_forces(qs_env)
128
129 TYPE(qs_environment_type), POINTER :: qs_env
130
131 CHARACTER(len=*), PARAMETER :: routinen = 'qs_forces'
132
133 INTEGER :: after, handle, i, iatom, ic, ikind, &
134 ispin, iw, natom, nkind, nspin, &
135 output_unit
136 INTEGER, ALLOCATABLE, DIMENSION(:) :: atom_of_kind, kind_of, natom_of_kind
137 LOGICAL :: do_admm, do_exx, do_gw, do_im_time, &
138 has_unit_metric, omit_headers, &
139 perform_ec, reuse_hfx
140 REAL(dp) :: dummy_real, dummy_real2(2)
141 TYPE(atomic_kind_type), DIMENSION(:), POINTER :: atomic_kind_set
142 TYPE(cp_logger_type), POINTER :: logger
143 TYPE(dbcsr_p_type), DIMENSION(:), POINTER :: matrix_s, matrix_w, rho_ao
144 TYPE(dbcsr_p_type), DIMENSION(:, :), POINTER :: matrix_w_kp
145 TYPE(dft_control_type), POINTER :: dft_control
146 TYPE(energy_correction_type), POINTER :: ec_env
147 TYPE(lri_environment_type), POINTER :: lri_env
148 TYPE(mp_para_env_type), POINTER :: para_env
149 TYPE(particle_type), DIMENSION(:), POINTER :: particle_set
150 TYPE(qs_energy_type), POINTER :: energy
151 TYPE(qs_force_type), DIMENSION(:), POINTER :: force
152 TYPE(qs_ks_env_type), POINTER :: ks_env
153 TYPE(qs_rho_type), POINTER :: rho
154 TYPE(qs_subsys_type), POINTER :: subsys
155 TYPE(section_vals_type), POINTER :: hfx_sections, print_section
156 TYPE(virial_type), POINTER :: virial
157
158 CALL timeset(routinen, handle)
159 NULLIFY (logger)
160 logger => cp_get_default_logger()
161
162 ! rebuild plane wave environment
163 CALL qs_env_rebuild_pw_env(qs_env)
164
165 ! zero out the forces in particle set
166 CALL get_qs_env(qs_env, particle_set=particle_set)
167 natom = SIZE(particle_set)
168 DO iatom = 1, natom
169 particle_set(iatom)%f = 0.0_dp
170 END DO
171
172 ! get atom mapping
173 NULLIFY (atomic_kind_set)
174 CALL get_qs_env(qs_env, atomic_kind_set=atomic_kind_set)
175 CALL get_atomic_kind_set(atomic_kind_set=atomic_kind_set, &
176 atom_of_kind=atom_of_kind, &
177 kind_of=kind_of)
178
179 NULLIFY (force, subsys, dft_control)
180 CALL get_qs_env(qs_env, &
181 force=force, &
182 subsys=subsys, &
183 dft_control=dft_control)
184 IF (.NOT. ASSOCIATED(force)) THEN
185 ! *** Allocate the force data structure ***
186 nkind = SIZE(atomic_kind_set)
187 CALL get_atomic_kind_set(atomic_kind_set=atomic_kind_set, natom_of_kind=natom_of_kind)
188 CALL allocate_qs_force(force, natom_of_kind)
189 DEALLOCATE (natom_of_kind)
190 CALL qs_subsys_set(subsys, force=force)
191 END IF
192 CALL zero_qs_force(force)
193
194 ! Check if CDFT potential is needed and save it until forces have been calculated
195 IF (dft_control%qs_control%cdft) THEN
196 dft_control%qs_control%cdft_control%save_pot = .true.
197 END IF
198
199 ! recalculate energy with forces
200 CALL qs_energies(qs_env, calc_forces=.true.)
201
202 NULLIFY (para_env)
203 CALL get_qs_env(qs_env, &
204 para_env=para_env)
205
206 ! Now we handle some special cases
207 ! Maybe some of these would be better dealt with in qs_energies?
208 IF (qs_env%run_rtp) THEN
209 NULLIFY (matrix_w, matrix_s, ks_env)
210 CALL get_qs_env(qs_env, &
211 ks_env=ks_env, &
212 matrix_w=matrix_w, &
213 matrix_s=matrix_s)
214 CALL dbcsr_allocate_matrix_set(matrix_w, dft_control%nspins)
215 DO ispin = 1, dft_control%nspins
216 ALLOCATE (matrix_w(ispin)%matrix)
217 CALL dbcsr_copy(matrix_w(ispin)%matrix, matrix_s(1)%matrix, &
218 name="W MATRIX")
219 CALL dbcsr_set(matrix_w(ispin)%matrix, 0.0_dp)
220 END DO
221 CALL set_ks_env(ks_env, matrix_w=matrix_w)
222
223 CALL calc_c_mat_force(qs_env)
224 IF (dft_control%do_admm) CALL rt_admm_force(qs_env)
225 IF (dft_control%rtp_control%velocity_gauge .AND. dft_control%rtp_control%nl_gauge_transform) &
226 CALL velocity_gauge_nl_force(qs_env, particle_set)
227 END IF
228 ! from an eventual Mulliken restraint
229 IF (dft_control%qs_control%mulliken_restraint) THEN
230 NULLIFY (matrix_w, matrix_s, rho)
231 CALL get_qs_env(qs_env, &
232 matrix_w=matrix_w, &
233 matrix_s=matrix_s, &
234 rho=rho)
235 NULLIFY (rho_ao)
236 CALL qs_rho_get(rho, rho_ao=rho_ao)
237 CALL mulliken_restraint(dft_control%qs_control%mulliken_restraint_control, &
238 para_env, matrix_s(1)%matrix, rho_ao, w_matrix=matrix_w)
239 END IF
240 ! Add non-Pulay contribution of DFT+U to W matrix, since it has also to be
241 ! digested with overlap matrix derivatives
242 IF (dft_control%dft_plus_u) THEN
243 NULLIFY (matrix_w)
244 CALL get_qs_env(qs_env, matrix_w=matrix_w)
245 CALL plus_u(qs_env=qs_env, matrix_w=matrix_w)
246 END IF
247
248 ! Write W Matrix to output (if requested)
249 CALL get_qs_env(qs_env, has_unit_metric=has_unit_metric)
250 IF (.NOT. has_unit_metric) THEN
251 NULLIFY (matrix_w_kp)
252 CALL get_qs_env(qs_env, matrix_w_kp=matrix_w_kp)
253 nspin = SIZE(matrix_w_kp, 1)
254 DO ispin = 1, nspin
255 IF (btest(cp_print_key_should_output(logger%iter_info, &
256 qs_env%input, "DFT%PRINT%AO_MATRICES/W_MATRIX"), cp_p_file)) THEN
257 iw = cp_print_key_unit_nr(logger, qs_env%input, "DFT%PRINT%AO_MATRICES/W_MATRIX", &
258 extension=".Log")
259 CALL section_vals_val_get(qs_env%input, "DFT%PRINT%AO_MATRICES%NDIGITS", i_val=after)
260 CALL section_vals_val_get(qs_env%input, "DFT%PRINT%AO_MATRICES%OMIT_HEADERS", l_val=omit_headers)
261 after = min(max(after, 1), 16)
262 DO ic = 1, SIZE(matrix_w_kp, 2)
263 CALL cp_dbcsr_write_sparse_matrix(matrix_w_kp(ispin, ic)%matrix, 4, after, qs_env, &
264 para_env, output_unit=iw, omit_headers=omit_headers)
265 END DO
266 CALL cp_print_key_finished_output(iw, logger, qs_env%input, &
267 "DFT%PRINT%AO_MATRICES/W_MATRIX")
268 END IF
269 END DO
270 END IF
271
272 ! Check if energy correction should be skipped
273 perform_ec = .false.
274 IF (qs_env%energy_correction) THEN
275 CALL get_qs_env(qs_env, ec_env=ec_env)
276 IF (.NOT. ec_env%do_skip) perform_ec = .true.
277 END IF
278
279 ! Compute core forces (also overwrites matrix_w)
280 IF (dft_control%qs_control%semi_empirical) THEN
281 CALL build_se_core_matrix(qs_env=qs_env, para_env=para_env, &
282 calculate_forces=.true.)
283 CALL se_core_core_interaction(qs_env, para_env, calculate_forces=.true.)
284 ELSEIF (dft_control%qs_control%dftb) THEN
285 CALL build_dftb_matrices(qs_env=qs_env, para_env=para_env, &
286 calculate_forces=.true.)
287 CALL calculate_dftb_dispersion(qs_env=qs_env, para_env=para_env, &
288 calculate_forces=.true.)
289 ELSEIF (dft_control%qs_control%xtb) THEN
290 CALL build_xtb_matrices(qs_env=qs_env, calculate_forces=.true.)
291 ELSEIF (perform_ec) THEN
292 ! Calculates core and grid based forces
293 CALL energy_correction(qs_env, ec_init=.false., calculate_forces=.true.)
294 ELSE
295 ! Dispersion energy and forces are calculated in qs_energy?
296 CALL build_core_hamiltonian_matrix(qs_env=qs_env, calculate_forces=.true.)
297 ! The above line reset the core H, which should be re-updated in case a TD field is applied:
298 IF (qs_env%run_rtp) THEN
299 IF (dft_control%apply_efield_field) &
300 CALL efield_potential_lengh_gauge(qs_env)
301 IF (dft_control%rtp_control%velocity_gauge) &
302 CALL velocity_gauge_ks_matrix(qs_env, subtract_nl_term=.false.)
303
304 END IF
305 CALL calculate_ecore_self(qs_env)
306 CALL calculate_ecore_overlap(qs_env, para_env, calculate_forces=.true.)
307 CALL calculate_ecore_efield(qs_env, calculate_forces=.true.)
308 !swap external_e_potential before external_c_potential, to ensure
309 !that external potential on grid is loaded before calculating energy of cores
310 CALL external_e_potential(qs_env)
311 IF (.NOT. dft_control%qs_control%gapw) THEN
312 CALL external_c_potential(qs_env, calculate_forces=.true.)
313 END IF
314 ! RIGPW matrices
315 IF (dft_control%qs_control%rigpw) THEN
316 CALL get_qs_env(qs_env=qs_env, lri_env=lri_env)
317 CALL build_ri_matrices(lri_env, qs_env, calculate_forces=.true.)
318 END IF
319 END IF
320
321 ! MP2 Code
322 IF (ASSOCIATED(qs_env%mp2_env)) THEN
323 NULLIFY (energy)
324 CALL get_qs_env(qs_env, energy=energy)
325 CALL qs_scf_compute_properties(qs_env, wf_type='MP2 ', do_mp2=.true.)
326 CALL qs_ks_update_qs_env(qs_env, just_energy=.true.)
327 energy%total = energy%total + energy%mp2
328
329 IF ((qs_env%mp2_env%method == ri_mp2_method_gpw .OR. qs_env%mp2_env%method == ri_mp2_laplace .OR. &
330 qs_env%mp2_env%method == ri_rpa_method_gpw) &
331 .AND. .NOT. qs_env%mp2_env%do_im_time) THEN
332 CALL update_mp2_forces(qs_env)
333 END IF
334
335 !RPA EXX energy and forces
336 IF (qs_env%mp2_env%method == ri_rpa_method_gpw) THEN
337 do_exx = .false.
338 hfx_sections => section_vals_get_subs_vals(qs_env%input, "DFT%XC%WF_CORRELATION%RI_RPA%HF")
339 CALL section_vals_get(hfx_sections, explicit=do_exx)
340 IF (do_exx) THEN
341 do_gw = qs_env%mp2_env%ri_rpa%do_ri_g0w0
342 do_admm = qs_env%mp2_env%ri_rpa%do_admm
343 reuse_hfx = qs_env%mp2_env%ri_rpa%reuse_hfx
344 do_im_time = qs_env%mp2_env%do_im_time
345 output_unit = cp_logger_get_default_io_unit()
346 dummy_real = 0.0_dp
347
348 CALL calculate_exx(qs_env=qs_env, &
349 unit_nr=output_unit, &
350 hfx_sections=hfx_sections, &
351 x_data=qs_env%mp2_env%ri_rpa%x_data, &
352 do_gw=do_gw, &
353 do_admm=do_admm, &
354 calc_forces=.true., &
355 reuse_hfx=reuse_hfx, &
356 do_im_time=do_im_time, &
357 e_ex_from_gw=dummy_real, &
358 e_admm_from_gw=dummy_real2, &
359 t3=dummy_real)
360 END IF
361 END IF
362 ELSEIF (perform_ec) THEN
363 ! energy correction forces postponed
364 ELSEIF (qs_env%harris_method) THEN
365 ! Harris method forces already done in harris_energy_correction
366 ELSE
367 ! Compute grid-based forces
368 CALL qs_ks_update_qs_env(qs_env, calculate_forces=.true.)
369 END IF
370
371 ! Excited state forces
372 CALL excited_state_energy(qs_env, calculate_forces=.true.)
373
374 ! replicate forces (get current pointer)
375 NULLIFY (force)
376 CALL get_qs_env(qs_env=qs_env, force=force)
377 CALL replicate_qs_force(force, para_env)
378
379 DO iatom = 1, natom
380 ikind = kind_of(iatom)
381 i = atom_of_kind(iatom)
382 ! XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
383 ! the force is - dE/dR, what is called force is actually the gradient
384 ! Things should have the right name
385 ! The minus sign below is a hack
386 ! XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
387 force(ikind)%other(1:3, i) = -particle_set(iatom)%f(1:3) + force(ikind)%ch_pulay(1:3, i)
388 force(ikind)%total(1:3, i) = force(ikind)%total(1:3, i) + force(ikind)%other(1:3, i)
389 particle_set(iatom)%f = -force(ikind)%total(1:3, i)
390 END DO
391
392 NULLIFY (virial, energy)
393 CALL get_qs_env(qs_env=qs_env, virial=virial, energy=energy)
394 IF (virial%pv_availability) THEN
395 CALL para_env%sum(virial%pv_overlap)
396 CALL para_env%sum(virial%pv_ekinetic)
397 CALL para_env%sum(virial%pv_ppl)
398 CALL para_env%sum(virial%pv_ppnl)
399 CALL para_env%sum(virial%pv_ecore_overlap)
400 CALL para_env%sum(virial%pv_ehartree)
401 CALL para_env%sum(virial%pv_exc)
402 CALL para_env%sum(virial%pv_exx)
403 CALL para_env%sum(virial%pv_vdw)
404 CALL para_env%sum(virial%pv_mp2)
405 CALL para_env%sum(virial%pv_nlcc)
406 CALL para_env%sum(virial%pv_gapw)
407 CALL para_env%sum(virial%pv_lrigpw)
408 CALL para_env%sum(virial%pv_virial)
409 CALL symmetrize_virial(virial)
410 ! Add the volume terms of the virial
411 IF ((.NOT. virial%pv_numer) .AND. &
412 (.NOT. (dft_control%qs_control%dftb .OR. &
413 dft_control%qs_control%xtb .OR. &
414 dft_control%qs_control%semi_empirical))) THEN
415
416 ! Harris energy correction requires volume terms from
417 ! 1) Harris functional contribution, and
418 ! 2) Linear Response solver
419 IF (perform_ec) THEN
420 CALL get_qs_env(qs_env, ec_env=ec_env)
421 energy%hartree = ec_env%ehartree
422 energy%exc = ec_env%exc
423 IF (dft_control%do_admm) THEN
424 energy%exc_aux_fit = ec_env%exc_aux_fit
425 END IF
426 END IF
427 DO i = 1, 3
428 virial%pv_ehartree(i, i) = virial%pv_ehartree(i, i) &
429 - 2.0_dp*(energy%hartree + energy%sccs_pol)
430 virial%pv_virial(i, i) = virial%pv_virial(i, i) - energy%exc &
431 - 2.0_dp*(energy%hartree + energy%sccs_pol)
432 virial%pv_exc(i, i) = virial%pv_exc(i, i) - energy%exc
433 IF (dft_control%do_admm) THEN
434 virial%pv_exc(i, i) = virial%pv_exc(i, i) - energy%exc_aux_fit
435 virial%pv_virial(i, i) = virial%pv_virial(i, i) - energy%exc_aux_fit
436 END IF
437 ! The factor 2 is a hack. It compensates the plus sign in h_stress/pw_poisson_solve.
438 ! The sign in pw_poisson_solve is correct for FIST, but not for QS.
439 ! There should be a more elegant solution to that ...
440 END DO
441 END IF
442 END IF
443
444 output_unit = cp_print_key_unit_nr(logger, qs_env%input, "DFT%PRINT%DERIVATIVES", &
445 extension=".Log")
446 print_section => section_vals_get_subs_vals(qs_env%input, "DFT%PRINT%DERIVATIVES")
447 IF (dft_control%qs_control%semi_empirical) THEN
448 CALL write_forces(force, atomic_kind_set, 2, output_unit=output_unit, &
449 print_section=print_section)
450 ELSE IF (dft_control%qs_control%dftb) THEN
451 CALL write_forces(force, atomic_kind_set, 4, output_unit=output_unit, &
452 print_section=print_section)
453 ELSE IF (dft_control%qs_control%xtb) THEN
454 CALL write_forces(force, atomic_kind_set, 4, output_unit=output_unit, &
455 print_section=print_section)
456 ELSE IF (dft_control%qs_control%gapw .OR. dft_control%qs_control%gapw_xc) THEN
457 CALL write_forces(force, atomic_kind_set, 1, output_unit=output_unit, &
458 print_section=print_section)
459 ELSE
460 CALL write_forces(force, atomic_kind_set, 0, output_unit=output_unit, &
461 print_section=print_section)
462 END IF
463 CALL cp_print_key_finished_output(output_unit, logger, qs_env%input, &
464 "DFT%PRINT%DERIVATIVES")
465
466 ! deallocate W Matrix:
467 NULLIFY (ks_env, matrix_w_kp)
468 CALL get_qs_env(qs_env=qs_env, &
469 matrix_w_kp=matrix_w_kp, &
470 ks_env=ks_env)
471 CALL dbcsr_deallocate_matrix_set(matrix_w_kp)
472 NULLIFY (matrix_w_kp)
473 CALL set_ks_env(ks_env, matrix_w_kp=matrix_w_kp)
474
475 DEALLOCATE (atom_of_kind, kind_of)
476
477 CALL timestop(handle)
478
479 END SUBROUTINE qs_forces
480
481! **************************************************************************************************
482!> \brief Write a Quickstep force data structure to output unit
483!> \param qs_force ...
484!> \param atomic_kind_set ...
485!> \param ftype ...
486!> \param output_unit ...
487!> \param print_section ...
488!> \date 05.06.2002
489!> \author MK
490!> \version 1.0
491! **************************************************************************************************
492 SUBROUTINE write_forces(qs_force, atomic_kind_set, ftype, output_unit, &
493 print_section)
494
495 TYPE(qs_force_type), DIMENSION(:), POINTER :: qs_force
496 TYPE(atomic_kind_type), DIMENSION(:), POINTER :: atomic_kind_set
497 INTEGER, INTENT(IN) :: ftype, output_unit
498 TYPE(section_vals_type), POINTER :: print_section
499
500 CHARACTER(LEN=13) :: fmtstr5
501 CHARACTER(LEN=15) :: fmtstr4
502 CHARACTER(LEN=20) :: fmtstr3
503 CHARACTER(LEN=35) :: fmtstr2
504 CHARACTER(LEN=48) :: fmtstr1
505 INTEGER :: i, iatom, ikind, my_ftype, natom, ndigits
506 INTEGER, ALLOCATABLE, DIMENSION(:) :: atom_of_kind, kind_of
507 REAL(kind=dp), DIMENSION(3) :: grand_total
508
509 IF (output_unit > 0) THEN
510
511 IF (.NOT. ASSOCIATED(qs_force)) THEN
512 CALL cp_abort(__location__, &
513 "The qs_force pointer is not associated "// &
514 "and cannot be printed")
515 END IF
516
517 CALL get_atomic_kind_set(atomic_kind_set=atomic_kind_set, atom_of_kind=atom_of_kind, &
518 kind_of=kind_of, natom=natom)
519
520 ! Variable precision output of the forces
521 CALL section_vals_val_get(print_section, "NDIGITS", &
522 i_val=ndigits)
523
524 fmtstr1 = "(/,/,T2,A,/,/,T3,A,T11,A,T23,A,T40,A1,2( X,A1))"
525 WRITE (unit=fmtstr1(41:42), fmt="(I2)") ndigits + 5
526
527 fmtstr2 = "(/,(T2,I5,4X,I4,T18,A,T34,3F . ))"
528 WRITE (unit=fmtstr2(32:33), fmt="(I2)") ndigits
529 WRITE (unit=fmtstr2(29:30), fmt="(I2)") ndigits + 6
530
531 fmtstr3 = "(/,T3,A,T34,3F . )"
532 WRITE (unit=fmtstr3(18:19), fmt="(I2)") ndigits
533 WRITE (unit=fmtstr3(15:16), fmt="(I2)") ndigits + 6
534
535 fmtstr4 = "((T34,3F . ))"
536 WRITE (unit=fmtstr4(12:13), fmt="(I2)") ndigits
537 WRITE (unit=fmtstr4(9:10), fmt="(I2)") ndigits + 6
538
539 fmtstr5 = "(/T2,A//T3,A)"
540
541 WRITE (unit=output_unit, fmt=fmtstr1) &
542 "FORCES [a.u.]", "Atom", "Kind", "Component", "X", "Y", "Z"
543
544 grand_total(:) = 0.0_dp
545
546 my_ftype = ftype
547
548 SELECT CASE (my_ftype)
549 CASE DEFAULT
550 DO iatom = 1, natom
551 ikind = kind_of(iatom)
552 i = atom_of_kind(iatom)
553 WRITE (unit=output_unit, fmt=fmtstr2) &
554 iatom, ikind, " total", qs_force(ikind)%total(1:3, i)
555 grand_total(1:3) = grand_total(1:3) + qs_force(ikind)%total(1:3, i)
556 END DO
557 CASE (0)
558 DO iatom = 1, natom
559 ikind = kind_of(iatom)
560 i = atom_of_kind(iatom)
561 WRITE (unit=output_unit, fmt=fmtstr2) &
562 iatom, ikind, " overlap", qs_force(ikind)%overlap(1:3, i), &
563 iatom, ikind, " overlap_admm", qs_force(ikind)%overlap_admm(1:3, i), &
564 iatom, ikind, " kinetic", qs_force(ikind)%kinetic(1:3, i), &
565 iatom, ikind, " gth_ppl", qs_force(ikind)%gth_ppl(1:3, i), &
566 iatom, ikind, " gth_nlcc", qs_force(ikind)%gth_nlcc(1:3, i), &
567 iatom, ikind, " gth_ppnl", qs_force(ikind)%gth_ppnl(1:3, i), &
568 iatom, ikind, " core_overlap", qs_force(ikind)%core_overlap(1:3, i), &
569 iatom, ikind, " rho_core", qs_force(ikind)%rho_core(1:3, i), &
570 iatom, ikind, " rho_elec", qs_force(ikind)%rho_elec(1:3, i), &
571 iatom, ikind, " rho_lri_elec", qs_force(ikind)%rho_lri_elec(1:3, i), &
572 iatom, ikind, " ch_pulay", qs_force(ikind)%ch_pulay(1:3, i), &
573 iatom, ikind, " dispersion", qs_force(ikind)%dispersion(1:3, i), &
574 iatom, ikind, " gCP", qs_force(ikind)%gcp(1:3, i), &
575 iatom, ikind, " other", qs_force(ikind)%other(1:3, i), &
576 iatom, ikind, " fock_4c", qs_force(ikind)%fock_4c(1:3, i), &
577 iatom, ikind, " ehrenfest", qs_force(ikind)%ehrenfest(1:3, i), &
578 iatom, ikind, " efield", qs_force(ikind)%efield(1:3, i), &
579 iatom, ikind, " eev", qs_force(ikind)%eev(1:3, i), &
580 iatom, ikind, " mp2_non_sep", qs_force(ikind)%mp2_non_sep(1:3, i), &
581 iatom, ikind, " total", qs_force(ikind)%total(1:3, i)
582 grand_total(1:3) = grand_total(1:3) + qs_force(ikind)%total(1:3, i)
583 END DO
584 CASE (1)
585 DO iatom = 1, natom
586 ikind = kind_of(iatom)
587 i = atom_of_kind(iatom)
588 WRITE (unit=output_unit, fmt=fmtstr2) &
589 iatom, ikind, " overlap", qs_force(ikind)%overlap(1:3, i), &
590 iatom, ikind, " overlap_admm", qs_force(ikind)%overlap_admm(1:3, i), &
591 iatom, ikind, " kinetic", qs_force(ikind)%kinetic(1:3, i), &
592 iatom, ikind, " gth_ppl", qs_force(ikind)%gth_ppl(1:3, i), &
593 iatom, ikind, " gth_nlcc", qs_force(ikind)%gth_nlcc(1:3, i), &
594 iatom, ikind, " gth_ppnl", qs_force(ikind)%gth_ppnl(1:3, i), &
595 iatom, ikind, " all_potential", qs_force(ikind)%all_potential(1:3, i), &
596 iatom, ikind, " core_overlap", qs_force(ikind)%core_overlap(1:3, i), &
597 iatom, ikind, " rho_core", qs_force(ikind)%rho_core(1:3, i), &
598 iatom, ikind, " rho_elec", qs_force(ikind)%rho_elec(1:3, i), &
599 iatom, ikind, " rho_lri_elec", qs_force(ikind)%rho_lri_elec(1:3, i), &
600 iatom, ikind, " vhxc_atom", qs_force(ikind)%vhxc_atom(1:3, i), &
601 iatom, ikind, " g0s_Vh_elec", qs_force(ikind)%g0s_Vh_elec(1:3, i), &
602 iatom, ikind, " ch_pulay", qs_force(ikind)%ch_pulay(1:3, i), &
603 iatom, ikind, " dispersion", qs_force(ikind)%dispersion(1:3, i), &
604 iatom, ikind, " gCP", qs_force(ikind)%gcp(1:3, i), &
605 iatom, ikind, " fock_4c", qs_force(ikind)%fock_4c(1:3, i), &
606 iatom, ikind, " ehrenfest", qs_force(ikind)%ehrenfest(1:3, i), &
607 iatom, ikind, " efield", qs_force(ikind)%efield(1:3, i), &
608 iatom, ikind, " eev", qs_force(ikind)%eev(1:3, i), &
609 iatom, ikind, " mp2_non_sep", qs_force(ikind)%mp2_non_sep(1:3, i), &
610 iatom, ikind, " total", qs_force(ikind)%total(1:3, i)
611 grand_total(1:3) = grand_total(1:3) + qs_force(ikind)%total(1:3, i)
612 END DO
613 CASE (2)
614 DO iatom = 1, natom
615 ikind = kind_of(iatom)
616 i = atom_of_kind(iatom)
617 WRITE (unit=output_unit, fmt=fmtstr2) &
618 iatom, ikind, " all_potential", qs_force(ikind)%all_potential(1:3, i), &
619 iatom, ikind, " rho_elec", qs_force(ikind)%rho_elec(1:3, i), &
620 iatom, ikind, " total", qs_force(ikind)%total(1:3, i)
621 grand_total(1:3) = grand_total(1:3) + qs_force(ikind)%total(1:3, i)
622 END DO
623 CASE (3)
624 DO iatom = 1, natom
625 ikind = kind_of(iatom)
626 i = atom_of_kind(iatom)
627 WRITE (unit=output_unit, fmt=fmtstr2) &
628 iatom, ikind, " overlap", qs_force(ikind)%overlap(1:3, i), &
629 iatom, ikind, "overlap_admm", qs_force(ikind)%overlap_admm(1:3, i), &
630 iatom, ikind, " kinetic", qs_force(ikind)%kinetic(1:3, i), &
631 iatom, ikind, " gth_ppl", qs_force(ikind)%gth_ppl(1:3, i), &
632 iatom, ikind, " gth_nlcc", qs_force(ikind)%gth_nlcc(1:3, i), &
633 iatom, ikind, " gth_ppnl", qs_force(ikind)%gth_ppnl(1:3, i), &
634 iatom, ikind, " core_overlap", qs_force(ikind)%core_overlap(1:3, i), &
635 iatom, ikind, " rho_core", qs_force(ikind)%rho_core(1:3, i), &
636 iatom, ikind, " rho_elec", qs_force(ikind)%rho_elec(1:3, i), &
637 iatom, ikind, " ch_pulay", qs_force(ikind)%ch_pulay(1:3, i), &
638 iatom, ikind, " fock_4c", qs_force(ikind)%fock_4c(1:3, i), &
639 iatom, ikind, " mp2_non_sep", qs_force(ikind)%mp2_non_sep(1:3, i), &
640 iatom, ikind, " total", qs_force(ikind)%total(1:3, i)
641 grand_total(1:3) = grand_total(1:3) + qs_force(ikind)%total(1:3, i)
642 END DO
643 CASE (4)
644 DO iatom = 1, natom
645 ikind = kind_of(iatom)
646 i = atom_of_kind(iatom)
647 WRITE (unit=output_unit, fmt=fmtstr2) &
648 iatom, ikind, " all_potential", qs_force(ikind)%all_potential(1:3, i), &
649 iatom, ikind, " overlap", qs_force(ikind)%overlap(1:3, i), &
650 iatom, ikind, " rho_elec", qs_force(ikind)%rho_elec(1:3, i), &
651 iatom, ikind, " repulsive", qs_force(ikind)%repulsive(1:3, i), &
652 iatom, ikind, " dispersion", qs_force(ikind)%dispersion(1:3, i), &
653 iatom, ikind, " efield", qs_force(ikind)%efield(1:3, i), &
654 iatom, ikind, " ehrenfest", qs_force(ikind)%ehrenfest(1:3, i), &
655 iatom, ikind, " total", qs_force(ikind)%total(1:3, i)
656 grand_total(1:3) = grand_total(1:3) + qs_force(ikind)%total(1:3, i)
657 END DO
658 CASE (5)
659 DO iatom = 1, natom
660 ikind = kind_of(iatom)
661 i = atom_of_kind(iatom)
662 WRITE (unit=output_unit, fmt=fmtstr2) &
663 iatom, ikind, " overlap", qs_force(ikind)%overlap(1:3, i), &
664 iatom, ikind, " kinetic", qs_force(ikind)%kinetic(1:3, i), &
665 iatom, ikind, " rho_elec", qs_force(ikind)%rho_elec(1:3, i), &
666 iatom, ikind, " dispersion", qs_force(ikind)%dispersion(1:3, i), &
667 iatom, ikind, " all potential", qs_force(ikind)%all_potential(1:3, i), &
668 iatom, ikind, " other", qs_force(ikind)%other(1:3, i), &
669 iatom, ikind, " total", qs_force(ikind)%total(1:3, i)
670 grand_total(1:3) = grand_total(1:3) + qs_force(ikind)%total(1:3, i)
671 END DO
672 END SELECT
673
674 WRITE (unit=output_unit, fmt=fmtstr3) "Sum of total", grand_total(1:3)
675
676 DEALLOCATE (atom_of_kind)
677 DEALLOCATE (kind_of)
678
679 END IF
680
681 END SUBROUTINE write_forces
682
683END MODULE qs_force
cp_dbcsr_operations::dbcsr_allocate_matrix_set
Definition cp_dbcsr_operations.F:77
cp_dbcsr_operations::dbcsr_deallocate_matrix_set
Definition cp_dbcsr_operations.F:85
atomic_kind_types
Define the atomic kind types and their sub types.
Definition atomic_kind_types.F:23
atomic_kind_types::get_atomic_kind_set
subroutine, public get_atomic_kind_set(atomic_kind_set, atom_of_kind, kind_of, natom_of_kind, maxatom, natom, nshell, fist_potential_present, shell_present, shell_adiabatic, shell_check_distance, damping_present)
Get attributes of an atomic kind set.
Definition atomic_kind_types.F:223
cp_control_types
Defines control structures, which contain the parameters and the settings for the DFT-based calculati...
Definition cp_control_types.F:12
cp_dbcsr_api
Definition cp_dbcsr_api.F:8
cp_dbcsr_api::dbcsr_copy
subroutine, public dbcsr_copy(matrix_b, matrix_a, name, keep_sparsity, keep_imaginary)
...
Definition cp_dbcsr_api.F:368
cp_dbcsr_api::dbcsr_set
subroutine, public dbcsr_set(matrix, alpha)
...
Definition cp_dbcsr_api.F:1181
cp_dbcsr_operations
DBCSR operations in CP2K.
Definition cp_dbcsr_operations.F:18
cp_dbcsr_output
DBCSR output in CP2K.
Definition cp_dbcsr_output.F:17
cp_dbcsr_output::cp_dbcsr_write_sparse_matrix
subroutine, public cp_dbcsr_write_sparse_matrix(sparse_matrix, before, after, qs_env, para_env, first_row, last_row, first_col, last_col, scale, output_unit, omit_headers)
...
Definition cp_dbcsr_output.F:163
cp_log_handling
various routines to log and control the output. The idea is that decisions about where to log should ...
Definition cp_log_handling.F:41
cp_log_handling::cp_logger_get_default_io_unit
integer function, public cp_logger_get_default_io_unit(logger)
returns the unit nr for the ionode (-1 on all other processors) skips as well checks if the procs cal...
Definition cp_log_handling.F:515
cp_log_handling::cp_get_default_logger
type(cp_logger_type) function, pointer, public cp_get_default_logger()
returns the default logger
Definition cp_log_handling.F:234
cp_output_handling
routines to handle the output, The idea is to remove the decision of wheter to output and what to out...
Definition cp_output_handling.F:25
cp_output_handling::cp_print_key_unit_nr
integer function, public cp_print_key_unit_nr(logger, basis_section, print_key_path, extension, middle_name, local, log_filename, ignore_should_output, file_form, file_position, file_action, file_status, do_backup, on_file, is_new_file, mpi_io, fout)
...
Definition cp_output_handling.F:853
cp_output_handling::cp_print_key_finished_output
subroutine, public cp_print_key_finished_output(unit_nr, logger, basis_section, print_key_path, local, ignore_should_output, on_file, mpi_io)
should be called after you finish working with a unit obtained with cp_print_key_unit_nr,...
Definition cp_output_handling.F:1063
cp_output_handling::cp_p_file
integer, parameter, public cp_p_file
Definition cp_output_handling.F:103
cp_output_handling::cp_print_key_should_output
integer function, public cp_print_key_should_output(iteration_info, basis_section, print_key_path, used_print_key, first_time)
returns what should be done with the given property if btest(res,cp_p_store) then the property should...
Definition cp_output_handling.F:345
dft_plus_u
Add the DFT+U contribution to the Hamiltonian matrix.
Definition dft_plus_u.F:18
dft_plus_u::plus_u
subroutine, public plus_u(qs_env, matrix_h, matrix_w)
Add the DFT+U contribution to the Hamiltonian matrix. Wrapper routine for all "+U" methods.
Definition dft_plus_u.F:99
ec_env_types
Types needed for a for a Energy Correction.
Definition ec_env_types.F:14
efield_utils
all routins needed for a nonperiodic electric field
Definition efield_utils.F:12
efield_utils::calculate_ecore_efield
subroutine, public calculate_ecore_efield(qs_env, calculate_forces)
Computes the force and the energy due to a efield on the cores Note: In the velocity gauge,...
Definition efield_utils.F:187
efield_utils::efield_potential_lengh_gauge
subroutine, public efield_potential_lengh_gauge(qs_env)
Replace the original implementation of the electric-electronic interaction in the length gauge....
Definition efield_utils.F:65
energy_corrections
Routines for an energy correction on top of a Kohn-Sham calculation.
Definition energy_corrections.F:18
energy_corrections::energy_correction
subroutine, public energy_correction(qs_env, ec_init, calculate_forces)
Energy Correction to a Kohn-Sham simulation Available energy corrections: (1) Harris energy functiona...
Definition energy_corrections.F:223
excited_states
Routines for total energy and forces of excited states.
Definition excited_states.F:14
excited_states::excited_state_energy
subroutine, public excited_state_energy(qs_env, calculate_forces)
Excited state energy and forces.
Definition excited_states.F:63
hfx_exx
Routines to calculate EXX in RPA and energy correction methods.
Definition hfx_exx.F:16
hfx_exx::calculate_exx
subroutine, public calculate_exx(qs_env, unit_nr, hfx_sections, x_data, do_gw, do_admm, calc_forces, reuse_hfx, do_im_time, e_ex_from_gw, e_admm_from_gw, t3)
...
Definition hfx_exx.F:106
input_constants
collects all constants needed in input so that they can be used without circular dependencies
Definition input_constants.F:17
input_constants::ri_rpa_method_gpw
integer, parameter, public ri_rpa_method_gpw
Definition input_constants.F:1059
input_constants::ri_mp2_method_gpw
integer, parameter, public ri_mp2_method_gpw
Definition input_constants.F:1059
input_constants::ri_mp2_laplace
integer, parameter, public ri_mp2_laplace
Definition input_constants.F:1059
input_section_types
objects that represent the structure of input sections and the data contained in an input section
Definition input_section_types.F:15
input_section_types::section_vals_get_subs_vals
recursive type(section_vals_type) function, pointer, public section_vals_get_subs_vals(section_vals, subsection_name, i_rep_section, can_return_null)
returns the values of the requested subsection
Definition input_section_types.F:731
input_section_types::section_vals_get
subroutine, public section_vals_get(section_vals, ref_count, n_repetition, n_subs_vals_rep, section, explicit)
returns various attributes about the section_vals
Definition input_section_types.F:704
input_section_types::section_vals_val_get
subroutine, public section_vals_val_get(section_vals, keyword_name, i_rep_section, i_rep_val, n_rep_val, val, l_val, i_val, r_val, c_val, l_vals, i_vals, r_vals, c_vals, explicit)
returns the requested value
Definition input_section_types.F:1047
kinds
Defines the basic variable types.
Definition kinds.F:23
kinds::dp
integer, parameter, public dp
Definition kinds.F:34
lri_environment_types
contains the types and subroutines for dealing with the lri_env lri : local resolution of the identit...
Definition lri_environment_types.F:18
message_passing
Interface to the message passing library MPI.
Definition message_passing.F:23
mp2_cphf
Routines to calculate CPHF like update and solve Z-vector equation for MP2 gradients (only GPW)
Definition mp2_cphf.F:14
mp2_cphf::update_mp2_forces
subroutine, public update_mp2_forces(qs_env)
...
Definition mp2_cphf.F:1305
mulliken
compute mulliken charges we (currently) define them as c_i = 1/2 [ (PS)_{ii} + (SP)_{ii} ]
Definition mulliken.F:13
mulliken::mulliken_restraint
subroutine, public mulliken_restraint(mulliken_restraint_control, para_env, s_matrix, p_matrix, energy, order_p, ks_matrix, w_matrix)
computes the energy and density matrix derivate of a constraint on the mulliken charges
Definition mulliken.F:73
particle_types
Define the data structure for the particle information.
Definition particle_types.F:19
qs_core_energies
Calculation of the energies concerning the core charge distribution.
Definition qs_core_energies.F:14
qs_core_energies::calculate_ecore_overlap
subroutine, public calculate_ecore_overlap(qs_env, para_env, calculate_forces, molecular, e_overlap_core, atecc)
Calculate the overlap energy of the core charge distribution.
Definition qs_core_energies.F:199
qs_core_energies::calculate_ecore_self
subroutine, public calculate_ecore_self(qs_env, e_self_core, atecc)
Calculate the self energy of the core charge distribution.
Definition qs_core_energies.F:357
qs_core_hamiltonian
Calculation of the core Hamiltonian integral matrix <a|H|b> over Cartesian Gaussian-type functions.
Definition qs_core_hamiltonian.F:77
qs_core_hamiltonian::build_core_hamiltonian_matrix
subroutine, public build_core_hamiltonian_matrix(qs_env, calculate_forces)
Cosntruction of the QS Core Hamiltonian Matrix.
Definition qs_core_hamiltonian.F:157
qs_dftb_dispersion
Calculation of dispersion in DFTB.
Definition qs_dftb_dispersion.F:12
qs_dftb_dispersion::calculate_dftb_dispersion
subroutine, public calculate_dftb_dispersion(qs_env, para_env, calculate_forces)
...
Definition qs_dftb_dispersion.F:64
qs_dftb_matrices
Calculation of Overlap and Hamiltonian matrices in DFTB.
Definition qs_dftb_matrices.F:12
qs_dftb_matrices::build_dftb_matrices
subroutine, public build_dftb_matrices(qs_env, para_env, calculate_forces)
...
Definition qs_dftb_matrices.F:98
qs_energy_types
Definition qs_energy_types.F:14
qs_energy
Perform a QUICKSTEP wavefunction optimization (single point)
Definition qs_energy.F:14
qs_energy::qs_energies
subroutine, public qs_energies(qs_env, consistent_energies, calc_forces)
Driver routine for QUICKSTEP single point wavefunction optimization.
Definition qs_energy.F:70
qs_environment_methods
qs_environment methods that use many other modules
Definition qs_environment_methods.F:15
qs_environment_methods::qs_env_rebuild_pw_env
subroutine, public qs_env_rebuild_pw_env(qs_env)
rebuilds the pw_env in the given qs_env, allocating it if necessary
Definition qs_environment_methods.F:193
qs_environment_types
Definition qs_environment_types.F:14
qs_environment_types::get_qs_env
subroutine, public get_qs_env(qs_env, atomic_kind_set, qs_kind_set, cell, super_cell, cell_ref, use_ref_cell, kpoints, dft_control, mos, sab_orb, sab_all, qmmm, qmmm_periodic, sac_ae, sac_ppl, sac_lri, sap_ppnl, sab_vdw, sab_scp, sap_oce, sab_lrc, sab_se, sab_xtbe, sab_tbe, sab_core, sab_xb, sab_xtb_pp, sab_xtb_nonbond, sab_almo, sab_kp, sab_kp_nosym, particle_set, energy, force, matrix_h, matrix_h_im, matrix_ks, matrix_ks_im, matrix_vxc, run_rtp, rtp, matrix_h_kp, matrix_h_im_kp, matrix_ks_kp, matrix_ks_im_kp, matrix_vxc_kp, kinetic_kp, matrix_s_kp, matrix_w_kp, matrix_s_ri_aux_kp, matrix_s, matrix_s_ri_aux, matrix_w, matrix_p_mp2, matrix_p_mp2_admm, rho, rho_xc, pw_env, ewald_env, ewald_pw, active_space, mpools, input, para_env, blacs_env, scf_control, rel_control, kinetic, qs_charges, vppl, rho_core, rho_nlcc, rho_nlcc_g, ks_env, ks_qmmm_env, wf_history, scf_env, local_particles, local_molecules, distribution_2d, dbcsr_dist, molecule_kind_set, molecule_set, subsys, cp_subsys, oce, local_rho_set, rho_atom_set, task_list, task_list_soft, rho0_atom_set, rho0_mpole, rhoz_set, ecoul_1c, rho0_s_rs, rho0_s_gs, do_kpoints, has_unit_metric, requires_mo_derivs, mo_derivs, mo_loc_history, nkind, natom, nelectron_total, nelectron_spin, efield, neighbor_list_id, linres_control, xas_env, virial, cp_ddapc_env, cp_ddapc_ewald, outer_scf_history, outer_scf_ihistory, x_data, et_coupling, dftb_potential, results, se_taper, se_store_int_env, se_nddo_mpole, se_nonbond_env, admm_env, lri_env, lri_density, exstate_env, ec_env, harris_env, dispersion_env, gcp_env, vee, rho_external, external_vxc, mask, mp2_env, bs_env, kg_env, wanniercentres, atprop, ls_scf_env, do_transport, transport_env, v_hartree_rspace, s_mstruct_changed, rho_changed, potential_changed, forces_up_to_date, mscfg_env, almo_scf_env, gradient_history, variable_history, embed_pot, spin_embed_pot, polar_env, mos_last_converged, eeq, rhs)
Get the QUICKSTEP environment.
Definition qs_environment_types.F:514
qs_external_potential
Routines to handle an external electrostatic field The external field can be generic and is provided ...
Definition qs_external_potential.F:12
qs_external_potential::external_c_potential
subroutine, public external_c_potential(qs_env, calculate_forces)
Computes the force and the energy due to the external potential on the cores.
Definition qs_external_potential.F:171
qs_external_potential::external_e_potential
subroutine, public external_e_potential(qs_env)
Computes the external potential on the grid.
Definition qs_external_potential.F:67
qs_force_types
Definition qs_force_types.F:13
qs_force_types::replicate_qs_force
subroutine, public replicate_qs_force(qs_force, para_env)
Replicate and sum up the force.
Definition qs_force_types.F:360
qs_force_types::zero_qs_force
subroutine, public zero_qs_force(qs_force)
Initialize a Quickstep force data structure.
Definition qs_force_types.F:250
qs_force_types::allocate_qs_force
subroutine, public allocate_qs_force(qs_force, natom_of_kind)
Allocate a Quickstep force data structure.
Definition qs_force_types.F:79
qs_force
Quickstep force driver routine.
Definition qs_force.F:12
qs_force::qs_calc_energy_force
subroutine, public qs_calc_energy_force(qs_env, calc_force, consistent_energies, linres)
...
Definition qs_force.F:107
qs_ks_methods
routines that build the Kohn-Sham matrix (i.e calculate the coulomb and xc parts
Definition qs_ks_methods.F:22
qs_ks_methods::qs_ks_update_qs_env
subroutine, public qs_ks_update_qs_env(qs_env, calculate_forces, just_energy, print_active)
updates the Kohn Sham matrix of the given qs_env (facility method)
Definition qs_ks_methods.F:1080
qs_ks_types
Definition qs_ks_types.F:15
qs_ks_types::set_ks_env
subroutine, public set_ks_env(ks_env, v_hartree_rspace, s_mstruct_changed, rho_changed, potential_changed, forces_up_to_date, complex_ks, matrix_h, matrix_h_im, matrix_ks, matrix_ks_im, matrix_vxc, kinetic, matrix_s, matrix_s_ri_aux, matrix_w, matrix_p_mp2, matrix_p_mp2_admm, matrix_h_kp, matrix_h_im_kp, matrix_ks_kp, matrix_vxc_kp, kinetic_kp, matrix_s_kp, matrix_w_kp, matrix_s_ri_aux_kp, matrix_ks_im_kp, vppl, rho_core, rho_nlcc, rho_nlcc_g, vee, neighbor_list_id, kpoints, sab_orb, sab_all, sac_ae, sac_ppl, sac_lri, sap_ppnl, sap_oce, sab_lrc, sab_se, sab_xtbe, sab_tbe, sab_core, sab_xb, sab_xtb_pp, sab_xtb_nonbond, sab_vdw, sab_scp, sab_almo, sab_kp, sab_kp_nosym, task_list, task_list_soft, subsys, dft_control, dbcsr_dist, distribution_2d, pw_env, para_env, blacs_env)
...
Definition qs_ks_types.F:572
qs_rho_types
superstucture that hold various representations of the density and keeps track of which ones are vali...
Definition qs_rho_types.F:18
qs_rho_types::qs_rho_get
subroutine, public qs_rho_get(rho_struct, rho_ao, rho_ao_im, rho_ao_kp, rho_ao_im_kp, rho_r, drho_r, rho_g, drho_g, tau_r, tau_g, rho_r_valid, drho_r_valid, rho_g_valid, drho_g_valid, tau_r_valid, tau_g_valid, tot_rho_r, tot_rho_g, rho_r_sccs, soft_valid, complex_rho_ao)
returns info about the density described by this object. If some representation is not available an e...
Definition qs_rho_types.F:229
qs_scf_post_scf
Utility routines for qs_scf.
Definition qs_scf_post_scf.F:11
qs_scf_post_scf::qs_scf_compute_properties
subroutine, public qs_scf_compute_properties(qs_env, wf_type, do_mp2)
computes properties for a given hamilonian using the current wfn
Definition qs_scf_post_scf.F:45
qs_subsys_types
types that represent a quickstep subsys
Definition qs_subsys_types.F:12
qs_subsys_types::qs_subsys_set
subroutine, public qs_subsys_set(subsys, cp_subsys, local_particles, local_molecules, cell, cell_ref, use_ref_cell, energy, force, qs_kind_set, nelectron_total, nelectron_spin)
...
Definition qs_subsys_types.F:226
ri_environment_methods
Calculates integral matrices for RIGPW method.
Definition ri_environment_methods.F:18
ri_environment_methods::build_ri_matrices
subroutine, public build_ri_matrices(lri_env, qs_env, calculate_forces)
creates and initializes an lri_env
Definition ri_environment_methods.F:90
rt_propagation_forces
Routines needed for EMD.
Definition rt_propagation_forces.F:13
rt_propagation_forces::rt_admm_force
subroutine, public rt_admm_force(qs_env)
...
Definition rt_propagation_forces.F:218
rt_propagation_forces::calc_c_mat_force
subroutine, public calc_c_mat_force(qs_env)
calculates the three additional force contributions needed in EMD P_imag*C , P_imag*B*S^-1*S_der ,...
Definition rt_propagation_forces.F:71
rt_propagation_velocity_gauge
Routines to perform the RTP in the velocity gauge.
Definition rt_propagation_velocity_gauge.F:12
rt_propagation_velocity_gauge::velocity_gauge_nl_force
subroutine, public velocity_gauge_nl_force(qs_env, particle_set)
Calculate the force associated to non-local pseudo potential in the velocity gauge.
Definition rt_propagation_velocity_gauge.F:893
rt_propagation_velocity_gauge::velocity_gauge_ks_matrix
subroutine, public velocity_gauge_ks_matrix(qs_env, subtract_nl_term)
...
Definition rt_propagation_velocity_gauge.F:93
se_core_core
Split and build its own idependent core_core SE interaction module.
Definition se_core_core.F:14
se_core_core::se_core_core_interaction
subroutine, public se_core_core_interaction(qs_env, para_env, calculate_forces)
Evaluates the core-core interactions for NDDO methods.
Definition se_core_core.F:79
se_core_matrix
Calculation of the Hamiltonian integral matrix <a|H|b> for semi-empirical methods.
Definition se_core_matrix.F:13
se_core_matrix::build_se_core_matrix
subroutine, public build_se_core_matrix(qs_env, para_env, calculate_forces)
...
Definition se_core_matrix.F:80
virial_types
Definition virial_types.F:13
virial_types::symmetrize_virial
subroutine, public symmetrize_virial(virial)
Symmetrize the virial components.
Definition virial_types.F:67
xtb_matrices
Calculation of Overlap and Hamiltonian matrices in xTB Reference: Stefan Grimme, Christoph Bannwarth,...
Definition xtb_matrices.F:15
xtb_matrices::build_xtb_matrices
subroutine, public build_xtb_matrices(qs_env, calculate_forces)
...
Definition xtb_matrices.F:111
atomic_kind_types::atomic_kind_type
Provides all information about an atomic kind.
Definition atomic_kind_types.F:45
cp_control_types::dft_control_type
Definition cp_control_types.F:570
cp_dbcsr_api::dbcsr_p_type
Definition cp_dbcsr_api.F:170
cp_log_handling::cp_logger_type
type of a logger, at the moment it contains just a print level starting at which level it should be l...
Definition cp_log_handling.F:140
ec_env_types::energy_correction_type
Contains information on the energy correction functional for KG.
Definition ec_env_types.F:55
input_section_types::section_vals_type
stores the values of a section
Definition input_section_types.F:127
lri_environment_types::lri_environment_type
Definition lri_environment_types.F:271
message_passing::mp_para_env_type
stores all the informations relevant to an mpi environment
Definition message_passing.F:763
particle_types::particle_type
Definition particle_types.F:35
qs_energy_types::qs_energy_type
Definition qs_energy_types.F:25
qs_environment_types::qs_environment_type
Definition qs_environment_types.F:217
qs_force_types::qs_force_type
Definition qs_force_types.F:28
qs_ks_types::qs_ks_env_type
calculation environment to calculate the ks matrix, holds all the needed vars. assumes that the core ...
Definition qs_ks_types.F:136
qs_rho_types::qs_rho_type
keeps the density in various representations, keeping track of which ones are valid.
Definition qs_rho_types.F:62
qs_subsys_types::qs_subsys_type
Definition qs_subsys_types.F:56
virial_types::virial_type
Definition virial_types.F:26