qs_dispersion_types.F

Go to the documentation of this file.

!--------------------------------------------------------------------------------------------------!
!   CP2K: A general program to perform molecular dynamics simulations                              !
!   Copyright 2000-2026 CP2K developers group <https://cp2k.org>                                   !
!                                                                                                  !
!   SPDX-License-Identifier: GPL-2.0-or-later                                                      !
!--------------------------------------------------------------------------------------------------!
 
! **************************************************************************************************
!> \brief Definition of disperson types for DFT calculations
!> \author JGH (20.10.2008)
! **************************************************************************************************
MODULE qs_dispersion_types
 
   USE eeq_input,                       ONLY: eeq_solver_type
   USE input_section_types,             ONLY: section_vals_type
   USE kinds,                           ONLY: default_path_length,&
                                              default_string_length,&
                                              dp
   USE qs_neighbor_list_types,          ONLY: neighbor_list_set_p_type,&
                                              release_neighbor_list_sets
#include "./base/base_uses.f90"
 
   IMPLICIT NONE
 
   PRIVATE
 
   CHARACTER(len=*), PARAMETER, PRIVATE :: moduleN = 'qs_dispersion_types'
 
! **************************************************************************************************
   INTEGER, PARAMETER                       :: dftd2_pp = 1
   INTEGER, PARAMETER                       :: dftd3_pp = 2
   INTEGER, PARAMETER                       :: dftd4_pp = 3
 
   TYPE qs_dispersion_type
      INTEGER                                :: TYPE = -1
      INTEGER                                :: pp_type = -1
      INTEGER                                :: nl_type = -1
      CHARACTER(LEN=default_string_length)   :: ref_functional = ""
      REAL(kind=dp)                          :: scaling = -1.0_dp
      REAL(kind=dp)                          :: rc_disp = -1.0_dp
      REAL(kind=dp)                          :: rc_d4 = -1.0_dp
      REAL(kind=dp)                          :: rc_cn = -1.0_dp
      REAL(kind=dp)                          :: d4_cutoff_width = 0.05_dp
      REAL(kind=dp)                          :: d4_3b_cutoff_width = 0.0_dp
      REAL(kind=dp)                          :: exp_pre = -1.0_dp
      TYPE(section_vals_type), POINTER       :: dftd_section => null()
      LOGICAL                                :: verbose = .false. !extended output
      CHARACTER(LEN=default_path_length)     :: parameter_file_name = ""
      CHARACTER(LEN=default_path_length)     :: kernel_file_name = ""
      !charges
      LOGICAL                                :: ext_charges = .false.
      REAL(kind=dp), DIMENSION(:), POINTER   :: charges => null() !charges for D4
      REAL(kind=dp), DIMENSION(:), POINTER   :: dcharges => null() !derivatives of D4 energy wrt charges
      TYPE(eeq_solver_type)                  :: eeq_sparam ! parameters for EEQ solver
      !DFT-D3 global parameters
      INTEGER                                :: max_elem = -1 !elements parametrized
      INTEGER                                :: maxc = -1 !max coordination number references per element
      REAL(kind=dp)                          :: k1 = -1.0_dp, k2 = -1.0_dp, k3 = -1.0_dp !ad hoc parameters
      REAL(kind=dp)                          :: alp = -1.0_dp !ad hoc parameters
      REAL(kind=dp)                          :: s6 = -1.0_dp, s8 = -1.0_dp, sr6 = -1.0_dp !scaling parameters
      REAL(kind=dp)                          :: sr8 = 1.0_dp !D3(0) C8 damping parameter
      REAL(kind=dp)                          :: a1 = -1.0_dp, a2 = -1.0_dp !BJ scaling parameters
      REAL(kind=dp)                          :: eps_cn = -1.0_dp
      LOGICAL                                :: d4_reference_code = .false. !Use D4 calculation from ext. library
      LOGICAL                                :: d4_debug = .false. !Output debug information for D4
      LOGICAL                                :: d3_reference_code = .false. !Use D3 parameters from s-dftd3 library
      LOGICAL                                :: d3_scaling_explicit = .false. !Scaling parameters were explicitly set
      LOGICAL                                :: doabc = .false. !neglect C9 terms
      LOGICAL                                :: c9cnst = .false. !use constant c9 terms
      LOGICAL                                :: lrc = .false. !calculate a long range correction
      LOGICAL                                :: srb = .false. !calculate a short range bond correction
      REAL(kind=dp), DIMENSION(4)            :: srb_params = -1.0_dp ! parameters for SRB (s,g,t1,t2)
      REAL(kind=dp)                          :: s9 = -1.0_dp !scale the many-body dispersion energy (default=1.0), dftd4
      TYPE(neighbor_list_set_p_type), &
         DIMENSION(:), POINTER               :: sab_vdw => null(), sab_cn => null() ! neighborlists for pair interactions
      REAL(kind=dp), DIMENSION(:, :, :, :, :), POINTER &
         :: c6ab => null()
      INTEGER, DIMENSION(:), POINTER         :: maxci => null()
      REAL(kind=dp), DIMENSION(:, :), POINTER :: r0ab => null()
      REAL(kind=dp), DIMENSION(:), POINTER   :: rcov => null() !covalent radii
      REAL(kind=dp), DIMENSION(:), POINTER   :: eneg => null() !electronegativity
      REAL(kind=dp), DIMENSION(:), POINTER   :: r2r4 => null() !atomic <r^2>/<r^4> values
      INTEGER                                :: cnfun = 1 ! CN function to be used
      REAL(kind=dp), DIMENSION(:), POINTER   :: cn => null()
      TYPE(cn_kind_list), DIMENSION(:), POINTER  :: cnkind => null()
      TYPE(cn_atom_list), DIMENSION(:), POINTER  :: cnlist => null()
      ! KG molecular corrections
      LOGICAL                                :: domol = .false.
      REAL(kind=dp)                          :: kgc8 = -1.0_dp !s8 scaling parameter
      !vdW-DF variables
      REAL(kind=dp)                          :: pw_cutoff = -1.0_dp
      !parameters for the rVV10 functional
      REAL(kind=dp)                          :: b_value = -1.0_dp, c_value = -1.0_dp, scale_rvv10 = -1.0_dp
      INTEGER                                :: nqs = -1, nr_points = -1
      !! The number of q points and radial points
      !! used in generating the kernel phi(q1*r, q2*r)
      !! (see DION 14-16 and SOLER 3)
      REAL(kind=dp)                          :: r_max = -1.0_dp, q_cut = -1.0_dp, q_min = -1.0_dp, dk = -1.0_dp
      !! The maximum value of r, the maximum and minimum
      !! values of q and the k-space spacing of grid points.
      !! Note that, during a vdW run, values of q0 found
      !! larger than q_cut will be saturated (SOLER 6-7) to
      !! q_cut
      REAL(kind=dp), DIMENSION(:), POINTER    :: q_mesh => null() !! The values of all the q points used
      REAL(kind=dp), DIMENSION(:, :, :), POINTER &
         :: kernel_table => null() !! Fourier transformed kernel and its second radial derivatives
      !! kernel_table(q1*(q1 - 1)/2 + q2, k_point, itable), with q2 <= q1 and k_point starting at zero.
      !! itable = 1 for kernel values and 2 for second derivatives.
      REAL(kind=dp), DIMENSION(:, :), POINTER  :: d2y_dx2 => null() !! 2nd derivatives of q_mesh for interpolation
      INTEGER, DIMENSION(:, :), POINTER  :: d3_exclude_pair => null()
      INTEGER  :: nd3_exclude_pair = -1
   END TYPE qs_dispersion_type
 
   TYPE qs_atom_dispersion_type
      INTEGER                                :: TYPE = -1
      LOGICAL                                :: defined = .false.
      REAL(kind=dp)                          :: vdw_radii = -1.0_dp !van der Waals radii
      REAL(kind=dp)                          :: c6 = -1.0_dp !c6 coefficients
   END TYPE qs_atom_dispersion_type
 
   TYPE cn_kind_list
      REAL(kind=dp)                          :: cnum = -1.0_dp
      INTEGER                                :: kind = -1
   END TYPE cn_kind_list
   TYPE cn_atom_list
      REAL(kind=dp)                          :: cnum = -1.0_dp
      INTEGER                                :: natom = -1
      INTEGER, DIMENSION(:), POINTER          :: atom => null()
   END TYPE cn_atom_list
 
! **************************************************************************************************
 
   PUBLIC :: qs_atom_dispersion_type, qs_dispersion_type
   PUBLIC :: dftd2_pp, dftd3_pp, dftd4_pp
   PUBLIC :: qs_dispersion_release
 
! **************************************************************************************************
CONTAINS
! **************************************************************************************************
!> \brief ...
!> \param dispersion_env ...
! **************************************************************************************************
   SUBROUTINE qs_dispersion_release(dispersion_env)
 
      TYPE(qs_dispersion_type), POINTER                  :: dispersion_env
 
      INTEGER                                            :: i
 
      IF (ASSOCIATED(dispersion_env)) THEN
         ! DFT-D3 arrays
         IF (ASSOCIATED(dispersion_env%maxci)) THEN
            DEALLOCATE (dispersion_env%maxci)
         END IF
         IF (ASSOCIATED(dispersion_env%c6ab)) THEN
            DEALLOCATE (dispersion_env%c6ab)
         END IF
         IF (ASSOCIATED(dispersion_env%r0ab)) THEN
            DEALLOCATE (dispersion_env%r0ab)
         END IF
         IF (ASSOCIATED(dispersion_env%rcov)) THEN
            DEALLOCATE (dispersion_env%rcov)
         END IF
         IF (ASSOCIATED(dispersion_env%eneg)) THEN
            DEALLOCATE (dispersion_env%eneg)
         END IF
         IF (ASSOCIATED(dispersion_env%r2r4)) THEN
            DEALLOCATE (dispersion_env%r2r4)
         END IF
         IF (ASSOCIATED(dispersion_env%cn)) THEN
            DEALLOCATE (dispersion_env%cn)
         END IF
         IF (ASSOCIATED(dispersion_env%cnkind)) THEN
            DEALLOCATE (dispersion_env%cnkind)
         END IF
         IF (ASSOCIATED(dispersion_env%cnlist)) THEN
            DO i = 1, SIZE(dispersion_env%cnlist)
               DEALLOCATE (dispersion_env%cnlist(i)%atom)
            END DO
            DEALLOCATE (dispersion_env%cnlist)
         END IF
         ! vdD-DF
         IF (ASSOCIATED(dispersion_env%q_mesh)) THEN
            DEALLOCATE (dispersion_env%q_mesh)
         END IF
         IF (ASSOCIATED(dispersion_env%kernel_table)) THEN
            DEALLOCATE (dispersion_env%kernel_table)
         END IF
         IF (ASSOCIATED(dispersion_env%d2y_dx2)) THEN
            DEALLOCATE (dispersion_env%d2y_dx2)
         END IF
         IF (ASSOCIATED(dispersion_env%d3_exclude_pair)) THEN
            DEALLOCATE (dispersion_env%d3_exclude_pair)
         END IF
         ! neighborlists
         CALL release_neighbor_list_sets(dispersion_env%sab_vdw)
         CALL release_neighbor_list_sets(dispersion_env%sab_cn)
         ! charges
         IF (ASSOCIATED(dispersion_env%charges)) THEN
            DEALLOCATE (dispersion_env%charges)
         END IF
         IF (ASSOCIATED(dispersion_env%dcharges)) THEN
            DEALLOCATE (dispersion_env%dcharges)
         END IF
 
         DEALLOCATE (dispersion_env)
 
      END IF
 
   END SUBROUTINE qs_dispersion_release
 
END MODULE qs_dispersion_types