qs_kernel_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                                                      !
!--------------------------------------------------------------------------------------------------!
 
MODULE qs_kernel_types
   USE admm_types,                      ONLY: admm_env_release,&
                                              admm_type
   USE cp_control_types,                ONLY: admm_control_release,&
                                              admm_control_type
   USE hfx_types,                       ONLY: hfx_release,&
                                              hfx_type
   USE input_section_types,             ONLY: section_vals_type
   USE kinds,                           ONLY: dp
   USE lri_environment_types,           ONLY: lri_density_type,&
                                              lri_environment_type
   USE qs_rho_types,                    ONLY: qs_rho_release,&
                                              qs_rho_type
   USE qs_tddfpt2_stda_types,           ONLY: stda_env_type
   USE xc_derivative_set_types,         ONLY: xc_derivative_set_type,&
                                              xc_dset_release
   USE xc_rho_set_types,                ONLY: xc_rho_set_release,&
                                              xc_rho_set_type
#include "./base/base_uses.f90"
 
   IMPLICIT NONE
 
   PRIVATE
 
   CHARACTER(len=*), PARAMETER, PRIVATE :: moduleN = 'qs_kernel_types'
 
   LOGICAL, PARAMETER, PRIVATE          :: debug_this_module = .false.
   ! number of first derivative components (3: d/dx, d/dy, d/dz)
   INTEGER, PARAMETER, PRIVATE          :: nderivs = 3
   INTEGER, PARAMETER, PRIVATE          :: maxspins = 2
 
   PUBLIC :: full_kernel_env_type, kernel_env_type
   PUBLIC :: release_kernel_env
 
! **************************************************************************************************
!> \brief Collection of variables required to evaluate adiabatic TDDFPT kernel.
!> \par History
!>   * 12.2016 created [Sergey Chulkov]
! **************************************************************************************************
   TYPE full_kernel_env_type
      ! ground state electron density
      TYPE(xc_rho_set_type), POINTER                     :: xc_rho_set => null()
      TYPE(qs_rho_type), POINTER                         :: rho0_struct => null()
      ! response density
      TYPE(xc_rho_set_type), POINTER                     :: xc_rho1_set => null()
      !> first and second derivatives of exchange-correlation functional
      TYPE(xc_derivative_set_type)                       :: xc_deriv_set = xc_derivative_set_type()
      !> XC input section
      LOGICAL                                            :: do_exck = .false.
      TYPE(section_vals_type), POINTER                   :: xc_section => null()
      !> scaling coefficients in the linear combination:
      !> K = alpha * K_{\alpha,\alpha} + beta * K_{\alpha,\beta}
      REAL(kind=dp)                                      :: alpha = -1.0_dp, beta = -1.0_dp
      !> flags for finite differences/analytic XC kernels
      LOGICAL                                            :: deriv2_analytic = .false.
      LOGICAL                                            :: deriv3_analytic = .false.
      ! Local resolution of the identity for Coulomb
      TYPE(lri_environment_type), POINTER                :: lri_env => null()
      TYPE(lri_density_type), POINTER                    :: lri_density => null()
      ! Short range HFX integral environment
      TYPE(hfx_type), DIMENSION(:, :), POINTER           :: x_data => null()
      TYPE(section_vals_type), POINTER                   :: hfxsr_section => null()
      TYPE(admm_type), POINTER                           :: admm_env => null()
      TYPE(admm_control_type), POINTER                   :: admm_control => null()
   END TYPE full_kernel_env_type
! **************************************************************************************************
!> \brief Type to hold environments for the different kernels
!> \par History
!>   * 04.2019 created [JHU]
! **************************************************************************************************
   TYPE kernel_env_type
      TYPE(full_kernel_env_type), POINTER                :: full_kernel => null()
      TYPE(full_kernel_env_type), POINTER                :: admm_kernel => null()
      TYPE(stda_env_type), POINTER                       :: stda_kernel => null()
   END TYPE kernel_env_type
 
CONTAINS
 
! **************************************************************************************************
!> \brief Release kernel environment.
!> \param kernel_env  kernel environment (destroyed on exit)
!> \par History
!>    * 02.2017 created [Sergey Chulkov]
! **************************************************************************************************
   SUBROUTINE release_kernel_env(kernel_env)
      TYPE(full_kernel_env_type), POINTER                :: kernel_env
 
      IF (ASSOCIATED(kernel_env)) THEN
         IF (ASSOCIATED(kernel_env%xc_rho1_set)) THEN
            CALL xc_rho_set_release(kernel_env%xc_rho1_set)
            DEALLOCATE (kernel_env%xc_rho1_set)
         END IF
         IF (ASSOCIATED(kernel_env%rho0_struct)) THEN
            CALL qs_rho_release(kernel_env%rho0_struct)
            DEALLOCATE (kernel_env%rho0_struct)
         END IF
         CALL xc_dset_release(kernel_env%xc_deriv_set)
         IF (ASSOCIATED(kernel_env%xc_rho_set)) THEN
            CALL xc_rho_set_release(kernel_env%xc_rho_set)
            DEALLOCATE (kernel_env%xc_rho_set)
         END IF
         IF (ASSOCIATED(kernel_env%x_data)) THEN
            CALL hfx_release(kernel_env%x_data)
         END IF
         IF (ASSOCIATED(kernel_env%admm_env)) THEN
            CALL admm_env_release(kernel_env%admm_env)
         END IF
         IF (ASSOCIATED(kernel_env%admm_control)) THEN
            CALL admm_control_release(kernel_env%admm_control)
         END IF
      END IF
 
   END SUBROUTINE release_kernel_env
 
END MODULE qs_kernel_types