Compute electronic entropy contribution.
S_elec = R * ln(2S+1) where 2S+1 is the spin multiplicity. For singlet (mult=1): S_elec = 0
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| integer, | intent(in) | :: | spin_multiplicity |
Electronic spin multiplicity (2S+1) |
||
| real(kind=dp), | intent(out) | :: | S_elec |
Electronic entropy in cal/(mol*K) |
pure subroutine compute_electronic_entropy(spin_multiplicity, S_elec) !! Compute electronic entropy contribution. !! !! S_elec = R * ln(2S+1) where 2S+1 is the spin multiplicity. !! For singlet (mult=1): S_elec = 0 integer, intent(in) :: spin_multiplicity !! Electronic spin multiplicity (2S+1) real(dp), intent(out) :: S_elec !! Electronic entropy in cal/(mol*K) S_elec = R_CALMOLK*log(real(spin_multiplicity, dp)) end subroutine compute_electronic_entropy