sarkas.tools.observables.RadialDistributionFunction.compute_sum_rule_integrals
sarkas.tools.observables.RadialDistributionFunction.compute_sum_rule_integrals#
- RadialDistributionFunction.compute_sum_rule_integrals(potential)[source]#
Compute integrals of the RDF used in sum rules.
The species dependent integrals are
\[I_{AB}^{\rm (Hartree, k)} = 2^{D - 2} \pi n_{A} n_{B} \int_0^{\infty} dr \, r^{D - 1 + k} \frac{d^k}{dr^k} \phi_{AB}(r),\]\[I_{AB}^{\rm (Corr, k)} = 2^{D - 2} \pi n_{A} n_{B} \int_0^{\infty} dr \, r^{D - 1 + k} h_{AB} (r) \frac{d^k}{dr^k} \phi_{AB}(r),\]where \(D\) is the number of dimensions, \(k = {0, 1, 2}\), and \(\phi_{AB}(r)\) is the potential between species \(A\) and \(B\).
Only Coulomb and Yukawa potentials are supported at the moment.
- Parameters
potential (
sarkas.potentials.core.Potential
) – Sarkas Potential object. Needed for all its attributes.- Returns
hartrees (numpy.ndarray) – Hartree integrals with \(k = {0, 1, 2}\).
Shape = (
sarkas.tools.observables.Observable.no_obs
, 3).corrs (numpy.ndarray) – Correlational integrals with \(k = {0, 1, 2}\).
Shape = (
sarkas.tools.observables.Observable.no_obs
, 3).