vdw_lennard_jones_Rijgrad Subroutine

private subroutine vdw_lennard_jones_Rijgrad(Rij, Rij0, Eij, Rijgrad)

Arguments

Type IntentOptional Attributes Name
real(kind=rp), intent(in) :: Rij
real(kind=rp), intent(in) :: Rij0
real(kind=rp), intent(in) :: Eij
real(kind=rp), intent(out) :: Rijgrad

Contents

Source Code

vdw_lennard_jones_Rijgrad

Source Code

    subroutine vdw_lennard_jones_Rijgrad(Rij, Rij0, Eij, Rijgrad)
        implicit none

        real(rp), intent(in) :: Rij
        real(rp), intent(in) :: Rij0
        real(rp), intent(in) :: Eij
        real(rp), intent(out) :: Rijgrad

        real(rp) :: sigma_ov_r

        sigma_ov_r = Rij0 / Rij
        Rijgrad = -12.0 * Eij * (sigma_ov_r ** 12 - sigma_ov_r ** 6) / Rij

    end subroutine vdw_lennard_jones_Rijgrad