Considering the equilibrium expression for E k in Eq. ( ), where Rk = Rs k + Rf k , a corresponding nonequilibrium expression is obtained by instead using Rs 0 + Rf k as the reaction-field operator for excited state k, corresponding to a state-specific polarization charge density σ k = σ s 0 + σ f k . This substitution affords