@brief Get third derivative of the XC functional contracted with response densities, spin-polarized version @param[in] xce XC engine @param[in] d_r \delta_rho (alpha, beta) @param[in] d_s \delta_sigma (alpha-alpha, beta-beta, alpha-beta) @param[in] d_t \delta_tau (alpha, beta) @param[in] ss (\nabla\rho(T)\nabla\rho(T)) (alpha-alpha, beta-beta, alpha-beta) @param[out] g_r \sum_i,j d2E_xc / (d_rho * d_zeta_id_zeta_j) (alpha, beta) @param[out] g_s \sum_i,j d2E_xc / (d_sigma * d_zeta_id_zeta_j) (alpha-alpha, beta-beta, alpha-beta) @param[out] g_t \sum_i,j d2E_xc / (d_tau * d_zeta_id_zeta_j) (alpha, beta)
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| class(xc_engine_t) | :: | xce | ||||
| integer, | intent(in) | :: | ipt | |||
| real(kind=fp), | intent(in) | :: | dr(2) | |||
| real(kind=fp), | intent(in) | :: | ds(3) | |||
| real(kind=fp), | intent(in) | :: | dt(2) | |||
| real(kind=fp), | intent(in) | :: | ss(3) | |||
| real(kind=fp), | intent(out) | :: | f_s(3) | |||
| real(kind=fp), | intent(out) | :: | g_r(2) | |||
| real(kind=fp), | intent(out) | :: | g_s(3) | |||
| real(kind=fp), | intent(out) | :: | g_t(2) |