calculate_diagonal_centrifugal_element Function

private function calculate_diagonal_centrifugal_element(total_angular_momentum_, j_, omega_) result(diagonal_element_)

calculates diagonal element of the centrifgual matrix, see Eq. 4 in "Coupling Matrix" section

Arguments

Type IntentOptional Attributes Name
integer(kind=int32), intent(in) :: total_angular_momentum_

total angular momentum
integer(kind=int32), intent(in) :: j_

rotational angular momentum
integer(kind=int32), intent(in) :: omega_

\(\bar{\Omega}\)

Return Value real(kind=dp)

(output) value of the diagonal element of the centrifgual matrix


Called by

proc~~calculate_diagonal_centrifugal_element~~CalledByGraph proc~calculate_diagonal_centrifugal_element calculate_diagonal_centrifugal_element proc~calculate_centrifugal_matrix calculate_centrifugal_matrix proc~calculate_centrifugal_matrix->proc~calculate_diagonal_centrifugal_element proc~initial_setup initial_setup proc~initial_setup->proc~calculate_centrifugal_matrix proc~numerov numerov proc~numerov->proc~initial_setup program~scattering SCATTERING program~scattering->proc~numerov

Contents

Source Code

calculate_diagonal_centrifugal_element

Source Code

      function calculate_diagonal_centrifugal_element(total_angular_momentum_, &
         j_, omega_) result(diagonal_element_)
         !! calculates diagonal element of the centrifgual matrix, see
         !! Eq. 4 in "Coupling Matrix" section
         !---------------------------------------------------------------------!
         integer(int32), intent(in) :: total_angular_momentum_
            !! total angular momentum
         integer(int32), intent(in) :: j_
            !! rotational angular momentum
         integer(int32), intent(in) :: omega_
            !! \\(\bar{\Omega}\\)
         real(dp) :: diagonal_element_
            !! (output) value of the diagonal element of the centrifgual matrix
         !---------------------------------------------------------------------!
         diagonal_element_ = real(                                             &
            total_angular_momentum_ * (total_angular_momentum_ + 1)            &
            + j_ * (j_ + 1) - 2 * omega_ **2, dp)
         !---------------------------------------------------------------------!
      end function calculate_diagonal_centrifugal_element