BaseCorePotential_module.f90

Go to the documentation of this file.

! Copyright 2019
!
! For a comprehensive list of the developers that contributed to these codes
! see the UK-AMOR website.
!
! This file is part of UKRmol-in (UKRmol+ suite).
!
!     UKRmol-in is free software: you can redistribute it and/or modify
!     it under the terms of the GNU General Public License as published by
!     the Free Software Foundation, either version 3 of the License, or
!     (at your option) any later version.
!
!     UKRmol-in is distributed in the hope that it will be useful,
!     but WITHOUT ANY WARRANTY; without even the implied warranty of
!     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
!     GNU General Public License for more details.
!
!     You should have received a copy of the GNU General Public License
!     along with  UKRmol-in (in source/COPYING). Alternatively, you can also visit
!     <https://www.gnu.org/licenses/>.
 
module basecorepotential_module
 
    use precisn, only: wp
 
    implicit none
 
    type, abstract :: basecorepotential
        real(wp), allocatable :: target_energies(:,:)
        real(wp)              :: l2_ecp_deficit
        integer               :: num_target_sym
        integer, allocatable  :: num_target_per_symmetries(:)
        integer, allocatable  :: spatial_symmetry(:)
        integer, allocatable  :: spin_symmetry(:)
        real(wp)              :: ground_energy
    contains
        procedure, public                   :: construct => construct_core
        procedure(generic_ecp),    deferred :: parse_ecp
        procedure(generic_input),  deferred :: parse_input
        procedure(generic_energy), deferred :: modify_target_energies
        procedure, public                   :: modify_ecp_diagonal
        procedure, public                   :: modify_ecp_diagonal_l2
    end type basecorepotential
 
    abstract interface
 
        subroutine generic_ecp (this, filename)
            import :: basecorepotential
            class(basecorepotential) :: this
            character(len=*)         :: filename
        end subroutine generic_ecp
    end interface
 
    abstract interface
 
        subroutine generic_input(this)
            import :: basecorepotential
            class(basecorepotential) :: this
        end subroutine generic_input
    end interface
 
    abstract interface
        subroutine generic_energy (this, target_sym, target_energy)
            import :: basecorepotential
            class(basecorepotential) :: this
            integer, intent(in)      :: target_sym
            real(wp)                 :: target_energy(:)
        end subroutine generic_energy
    end interface
 
contains
 
    subroutine construct_core (this, num_target_sym, num_target_per_symmetries, spatial_symmetry, spin_symmetry)
        class(basecorepotential) :: this
        integer, intent(in) :: num_target_sym
        integer, intent(in) :: num_target_per_symmetries(num_target_sym), &
                               spatial_symmetry(num_target_sym), &
                               spin_symmetry(num_target_sym)
 
        this % num_target_sym = num_target_sym
        allocate(this % num_target_per_symmetries(num_target_sym))
 
        this % num_target_per_symmetries = num_target_per_symmetries
        allocate(this % target_energies(maxval(this % num_target_per_symmetries), this % num_target_sym))
 
        allocate(this % spatial_symmetry(num_target_sym))
        allocate(this % spin_symmetry(num_target_sym))
 
        this % spatial_symmetry = spatial_symmetry
        this % spin_symmetry = spin_symmetry
 
        this % target_energies = 0.0
        this % L2_ECP_deficit =  0.0
    end subroutine construct_core
 
 
    subroutine modify_ecp_diagonal (this, value)
        class(basecorepotential), intent(in) :: this
        real(wp), intent(inout) :: value
 
        value = value - this % ground_energy
        value = value + this % target_energies(1,1)
    end subroutine modify_ecp_diagonal
 
 
    subroutine modify_ecp_diagonal_l2 (this, value)
        class(basecorepotential),intent(in) :: this
        real(wp), intent(inout) :: value
 
        value = value + this % L2_ECP_deficit
    end subroutine modify_ecp_diagonal_l2
 
end module basecorepotential_module