Table of Contents
This is an old revision of the document!
GetSlaterIntegrals
GetSlaterIntegrals(OrbitalNames,RadialWavefunctions1,RadialWavefunctions2) calculates all Slater integrals for the given orbitals.
Input
-
OrbitalNames : table of strings
-
RadialWavefunctions1 : table of InterpolatingFunction objects which belong to the orbitals determined by OrbitalNames
-
RadialWavefunctions2 : optional - table of InterpolatingFunction objects - if given, the function calculates relativistic Slater integrals where RadialWavefunctions1 are the large parts and RadialWavefunctions2 the small ones
Output
-
Rk(tau1,tau2,tau3,tau4) : dictionary of Slater integrals indexed by a string corresponding to the orbitals as given by OrbitalNames (see example)
Example
Given a grid r for the radial coordinate and the radial wave function R3dVal evaluated at that grid we can create interpolating functions and use them to calculate slater integrals
Input
- Example.Quanty
R3d = InterpolatingFunction.Spline(r,R3dVal) -- calculate Slater integrals for 3d shell Fk = GetSlaterIntegrals({"3d"},{R3d}) F0dd = Fk["3d 3d 3d 3d"][0]*EnergyUnits.Ha.value F2dd = Fk["3d 3d 3d 3d"][2]*EnergyUnits.Ha.value F4dd = Fk["3d 3d 3d 3d"][4]*EnergyUnits.Ha.value