# Derivation of Slater-Koster equations

I am trying to derive the Slater-Koster equations (Table. 1 of Ref. 1) for the two-centre approximation of hopping integrals between atomic orbitals. I understand that Slater-Koster approximates the two centre hopping integral as:

$$begin{equation} E_{n,m} = int{psi^*_n (textbf{r}-textbf{R})Hpsi_m (textbf{r}) dV},tag{1} end{equation}$$

assuming that potential, $$V$$ is spherically symmetric.

I tried substituting spherical harmonics, $$Y_{lm}$$ as the atomic orbitals $$psi (textbf{r})$$ but I am struggling to get to the result. First, I am not sure how to write the spherical harmonics $$Y_{lm}$$ not centred at the origin. Then, I can’t see how the integral would separate into two parts representing pi and sigma bonds.

Does anyone know how to derive one of the Slater-Koster equations? For example:

$$begin{equation} E_{p_z, d_{x^2-y^2}} = frac{1}{2}sqrt{3}n(l^2-m^2) V_{pdsigma} -n(l^2-m^2) V_{pdpi},tag{2} end{equation}$$

where $$l,m$$ and $$n$$ are directional cosines of connecting vector $$textbf{R}$$ (i.e. $$l=sin{beta}cos{alpha}=frac{R_x}{R}$$).

References

1. Slater, JC & Koster, GF, "Simplified LCAO Method for the Periodic Potential Problem" Physical Review (1954); link: https://journals.aps.org/pr/abstract/10.1103/PhysRev.94.1498

Matter Modeling Asked on August 19, 2021

The answer by Anyon does not address how to calculate the form of the Slater-Koster two-center integrals (2CI) but rather how to determine the direction cosines that appear in the 2CI's (the $$l$$,$$m$$,$$n$$'s in things like $$sqrt{3}lm*(sdsigma)$$, but not the $$sqrt{3}$$). I am currently trying to explicitly perform the same calculation to include in my comps so I can sympathize with your frustration in the lack of clarity in the literature...

Any who, the most explicit reference I have found that explains how to actually calculate the SK-matrix elements is in this paper (Takegahara et al 1980 J. Phys. C: Solid State Phys. 13 583, Slater-Koster tables for f electrons) where Takegahara extends the original SK table to f electrons.

The fundamental strategy he invokes, which is very tedious, is to calculate the form of hopping elements along the bond direction and then rotate it into a generic orientation to make use of the l,m,n- direction cosines we are familiar with.

Sorry that this is less of an answer and more of a comment, I don't have enough reputation to post any other way.

Answered by user3282375 on August 19, 2021

Spherical harmonics are not themselves full atomic orbitals. Consider the Hydrogen wave function, which separates into a radial part and an angular part. The latter is a spherical harmonic, but the former is some other function (in the case of Hydrogen it's a Laguerre polynomial). In general, we can approximate the angular part for other atoms with the same spherical harmonics, but we usually don't know the radial part analytically. Hence, if you want to directly evaluate Eq. (1) by integration you also need to find the radial part somewhere. The strength of Slater's and Koster's work, however, is that we can avoid this problem entirely, by hiding all radial dependence in the Slater-Koster parameters (in your example they are $$V_{pdsigma}$$ and $$V_{pdpi}$$).

For simplicity, I'll focus on the case of an $$s$$-$$p$$ overlap. This simplifies the geometry, and will let me borrow pictures and notation from this thesis. Let's say we have an $$s$$ orbital at site $$i$$ with wave function $$psi_{is}$$, and a $$p_alpha$$ orbital at site $$j$$ with wave function $$psi_{jp_alpha}$$, where $$alpha in {x,y,z}$$. In Dirac's bra-ket notation, the overlap between them can be written $$E_{s,p_alpha}=langle psi_{is}|H_{2c}|psi_{jp_alpha}rangle = langle S|H_{2c}|P_alpharangle, tag{3}$$ where $$H_{2c}$$ is the two-center Hamiltonian, and we've introduced short-hand notation for the the two wave functions. Of course, behind the bra-ket notation you have exactly the kind of overlap integral shown in Eq. (1).

Next step is to work out the geometry. Let $$vec{r}$$ be the vector connecting sites $$i$$ and $$j$$, and let $$vec{d}$$ be a unit vector along the same direction. We decompose the $$p$$ orbital at site $$j$$ into components parallel to ($$sigma$$) and perpendicular to ($$pi$$) the vector $$vec{d}$$, as shown in this figure:

Figure from Anthony Carlson's 2006 MSc thesis at University of Minnesota.

In this notation, $$sigma$$ and $$pi$$ (and also $$delta$$) is used to denote the component of angular momentum about the axis $$vec{d}$$. $$sigma$$ means zero, $$pi$$ means $$1$$, etc. To proceed, we also define $$vec{a}$$ as pointing along the $$p$$ orbital, and $$vec{n}$$ as a vector perpendicular to $$vec{d}$$ in the plane spanned by $$vec{a}$$ and $$vec{d}$$. Then, the $$p$$ orbital (at site $$j$$) can be decomposed $$|P_alpharangle = vec{a}cdotvec{d}|P_sigmarangle + vec{a}cdotvec{n}|P_pirangle. tag{4}$$ Then, the overlap is simply $$langle S|H_{2c}|P_alpharangle = left( vec{a}cdot vec{d} right) langle S | H_{2c} |P_sigmarangle + left( vec{a}cdotvec{n} right) langle S | H_{2c} | P_pirangle, tag{5}$$ where the second term is zero by symmetry. Further, we can easily express this in terms of directional cosines since we can without loss of generality choose $$vec{a}$$ to be parallel to one of the coordinate axes. With $$d_x=frac{vec{r}cdothat{x}}{|vec{r}|},quad d_y=frac{vec{r}cdothat{y}}{|vec{r}|},quad d_z=frac{vec{r}cdothat{z}}{|vec{r}|}, tag{6}$$ we get $$langle S|H_{2c}|P_alpharangle = d_alpha langle S|H_{2c}|P_alpharangle = d_alpha V_{spsigma}, tag{7}$$ where all radial overlap is hidden in the parameter $$V_{spsigma}$$. In tight-binding models it's common to denote this overlap integral $$t_{spsigma}$$ if it occurs in a hopping term.

The same approach works for other orbital combinations. You just need to set up the geometry and coordinate systems properly, and know where the atomic orbitals point. (Admittedly this can become quite complicated in some systems, e.g. transition metal oxides.) Then the Slater-Koster parameters can be treated as tuning parameters - either tuned to explore possible phenomena in some system, or fit to reproduce some experiment or calculated band structure.

Answered by Anyon on August 19, 2021

## Related Questions

### I get this error in Gaussian09 even though I have enough memory. “Error termination in NtrErr: NtrErr called from FIOCnC.” Any suggestions?

1  Asked on January 5, 2022 by quantumx

### Why are single excitations ignored in the MP2 component of double hybrid functional calculations?

1  Asked on January 5, 2022

### Regarding oscillatory strength theoretical units to experimental ones

1  Asked on January 3, 2022

### How well can we model chemical synthesis?

3  Asked on January 1, 2022 by anyon

### Is it required to sample the full reaction coordinate in umbrella sampling?

1  Asked on January 1, 2022 by bnd

### How to find good universities for matter modeling?

1  Asked on December 25, 2021

### What does the Neighbor command do in LAMMPS?

1  Asked on December 22, 2021

### For software that does not support FCIDUMP format, what format is used and how can we get the software to interact with FCIDUMP integrals?

2  Asked on December 22, 2021

### How do the various programs read or write integrals in FCIDUMP format?

1  Asked on December 20, 2021

### Is there an anisotropy factor (g factor) for TDDFT absorption and circular dichroism calculations?

1  Asked on December 20, 2021 by c-alexander

### Why am I losing bonds when I use PACKMOL on two proteins?

1  Asked on December 20, 2021

### Can we “invert” Density Functional Theory through sufficiently accurate experiment?

1  Asked on December 18, 2021 by kfc

### How do you calculate the “true” chemical potential in classical density functional theory?

1  Asked on December 16, 2021

### Converging the potential of a dilute binary alloy system using the Coherent Potential Approximation. What troubleshooting options could I take?

1  Asked on December 16, 2021 by luphys

### The Hund’s J – Why can this be quantified?

1  Asked on December 16, 2021

### How accurate are the most accurate calculations?

1  Asked on December 13, 2021 by roman-korol

### Is there any software that can generate the electronic DOS of GaAs?

1  Asked on December 13, 2021

### Which Linux distribution is best for Matter Modeling?

9  Asked on December 11, 2021

### What is the role of a chemist in computational chemistry?

4  Asked on December 9, 2021 by georgios-karaiskakis

### How to deduce phase transitions from a phonon calculation?

2  Asked on December 9, 2021 by koroma