neb_types.F

Go to the documentation of this file.

!--------------------------------------------------------------------------------------------------!
!   CP2K: A general program to perform molecular dynamics simulations                              !
!   Copyright 2000-2024 CP2K developers group <https://cp2k.org>                                   !
!                                                                                                  !
!   SPDX-License-Identifier: GPL-2.0-or-later                                                      !
!--------------------------------------------------------------------------------------------------!
 
! **************************************************************************************************
!> \brief Typo for Nudged Elastic Band Calculation
!> \note
!>      Numerical accuracy for parallel runs:
!>       Each replica starts the SCF run from the one optimized
!>       in a previous run. It may happen then energies and derivatives
!>       of a serial run and a parallel run could be slightly different
!>       'cause of a different starting density matrix.
!>       Exact results are obtained using:
!>          EXTRAPOLATION USE_GUESS in QS section (Teo 09.2006)
!> \author Teodoro Laino 10.2006
! **************************************************************************************************
MODULE neb_types
 
   USE force_env_types,                 ONLY: force_env_type
   USE input_constants,                 ONLY: band_md_opt,&
                                              do_b_neb,&
                                              do_band_cartesian,&
                                              do_band_collective,&
                                              pot_neb_full
   USE input_section_types,             ONLY: section_vals_type
   USE kinds,                           ONLY: default_string_length,&
                                              dp
#include "../base/base_uses.f90"
 
   IMPLICIT NONE
   PRIVATE
   CHARACTER(len=*), PARAMETER, PRIVATE :: moduleN = 'neb_types'
 
   PUBLIC :: neb_type, &
             neb_var_type, &
             neb_var_create, &
             neb_var_release
 
! **************************************************************************************************
   TYPE neb_type
      ! NEB parameters
      INTEGER       :: id_type = do_b_neb
      INTEGER       :: opt_type = band_md_opt
      INTEGER       :: pot_type = pot_neb_full
      INTEGER       :: number_of_replica = 0, nsize_xyz = 0, nsize_int = 0
      INTEGER       :: nsteps_it = 0, istep = 0
      INTEGER       :: nr_HE_image = 0
      LOGICAL       :: rotate_frames = .false., align_frames = .false.
      LOGICAL       :: optimize_end_points = .false.
      LOGICAL       :: use_colvar = .false.
      LOGICAL       :: reparametrize_frames = .false.
      INTEGER       :: spline_order = 0
      REAL(KIND=dp) :: k = 0.0_dp, spring_energy = 0.0_dp, avg_distance = 0.0_dp
      REAL(KIND=dp) :: smoothing = 0.0_dp
      CHARACTER(LEN=default_string_length) :: opt_type_label = ""
      ! Section used for restart
      TYPE(force_env_type), POINTER        :: force_env => null()
      TYPE(section_vals_type), POINTER     :: root_section => null()
      TYPE(section_vals_type), POINTER     :: motion_print_section => null()
      TYPE(section_vals_type), POINTER     :: force_env_section => null()
      TYPE(section_vals_type), POINTER     :: neb_section => null()
   END TYPE neb_type
 
! **************************************************************************************************
   TYPE neb_var_type
      INTEGER :: in_use = do_band_collective, size_wrk(2) = 0
      REAL(KIND=dp), DIMENSION(:, :), POINTER  :: xyz => null(), int => null(), wrk => null()
      REAL(KIND=dp), DIMENSION(:, :, :), POINTER  :: mmatrix => null()
   END TYPE neb_var_type
 
CONTAINS
! **************************************************************************************************
!> \brief Creates a variable type for BAND calculation
!> \param neb_var ...
!> \param neb_env ...
!> \param full_allocation ...
!> \date   05.2007
!> \author Teodoro Laino [tlaino] - University of Zurich
! **************************************************************************************************
   SUBROUTINE neb_var_create(neb_var, neb_env, full_allocation)
      TYPE(neb_var_type), POINTER                        :: neb_var
      TYPE(neb_type), POINTER                            :: neb_env
      LOGICAL, INTENT(IN), OPTIONAL                      :: full_allocation
 
      INTEGER                                            :: neb_nr_replica
      LOGICAL                                            :: allocate_all
 
      cpassert(.NOT. ASSOCIATED(neb_var))
      allocate_all = .false.
      IF (PRESENT(full_allocation)) allocate_all = full_allocation
      neb_nr_replica = neb_env%number_of_replica
      ALLOCATE (neb_var)
      NULLIFY (neb_var%xyz, neb_var%int, neb_var%wrk, neb_var%Mmatrix)
      IF (allocate_all) THEN
         ALLOCATE (neb_var%xyz(neb_env%nsize_xyz, neb_nr_replica))
         neb_var%xyz = 0.0_dp
      END IF
      IF (neb_env%use_colvar) THEN
         neb_var%in_use = do_band_collective
         cpassert(neb_env%nsize_int > 0)
         ALLOCATE (neb_var%int(neb_env%nsize_int, neb_nr_replica))
         neb_var%int = 0.0_dp
         neb_var%wrk => neb_var%int
      ELSE
         neb_var%in_use = do_band_cartesian
         IF (.NOT. allocate_all) THEN
            ALLOCATE (neb_var%xyz(neb_env%nsize_xyz, neb_nr_replica))
            neb_var%xyz = 0.0_dp
         END IF
         neb_var%wrk => neb_var%xyz
      END IF
      neb_var%size_wrk(1) = SIZE(neb_var%wrk, 1)
      neb_var%size_wrk(2) = SIZE(neb_var%wrk, 2)
 
   END SUBROUTINE neb_var_create
 
! **************************************************************************************************
!> \brief Releases a variable type for BAND calculation
!> \param neb_var ...
!> \date   05.2007
!> \author Teodoro Laino [tlaino] - University of Zurich
! **************************************************************************************************
   SUBROUTINE neb_var_release(neb_var)
      TYPE(neb_var_type), POINTER                        :: neb_var
 
      cpassert(ASSOCIATED(neb_var))
      IF (ASSOCIATED(neb_var%xyz)) THEN
         DEALLOCATE (neb_var%xyz)
      END IF
      IF (neb_var%in_use == do_band_collective) THEN
         DEALLOCATE (neb_var%int)
      END IF
      NULLIFY (neb_var%wrk)
      DEALLOCATE (neb_var)
 
   END SUBROUTINE neb_var_release
 
END MODULE neb_types