Check for positive MP2 correlation energies (instability warning) Correlation energies should be negative; positive values indicate instability
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| class(mp2_energy_t), | intent(in) | :: | this |
subroutine mp2_check_stability(this) !! Check for positive MP2 correlation energies (instability warning) !! Correlation energies should be negative; positive values indicate instability use pic_logger, only: logger => global_logger class(mp2_energy_t), intent(in) :: this if (this%ss > 0.0_dp) then call logger%warning("MP2 same-spin correlation energy is positive - possible instability!") end if if (this%os > 0.0_dp) then call logger%warning("MP2 opposite-spin correlation energy is positive - possible instability!") end if end subroutine mp2_check_stability