# Matrix elements evaluation in Hartree-Fock method

In Roothaan's matrix formulation of the Hartree-Fock method, we eventually must evaluate the elements in the Fock matrix in the form:

$$F_{\mu\nu} = \langle \mu | {-\frac{1}{2}}\nabla^2 | \nu \rangle - \sum_m^{\text{nuclei}} Z_m \langle \mu | \frac{1}{r_m} | \nu \rangle + \sum_{\lambda\sigma} P_{\lambda\sigma}[\langle \mu\nu|\lambda\sigma\rangle - \frac{1}{2}\langle \mu\lambda | \nu\sigma\rangle]$$

I'm having difficulty "descending" from the realm of abstract symbols to the actual number crunching that a computer must do to evaluate each term. At the end of the day, each term in the equation above is a number. My question is in particular about the second term, which represents the coulomb attraction between electrons and nuclei.

How is it computed, exactly? Is $$1/r_m$$ an operator? I know it represents the distance from the nucleus (for each nucleus, since it's inside a sum over all nuclei), but what is it? How do I start with a basis function and end up with a number that depends on the distance from the nucleus? The distance of what, by the way? The basis function is not a point in space . . .

I realize I'm not asking a clear question, but that's because I'm confused, so please bear with me and try to address some of the stuff I asked. Thanks guys.

• Yes, it's an operator and the integral - which is commonly is referred to as $\hat J$ - is then solved in polar coordinates or in Fourier space. Dec 24 '21 at 17:41

Electronic structure theory at the end of the day is "advanced quantum mechanics". IMO, one should understand the symbols from a QM perspective before applying it here.

Is $$1/r_m$$ an operator?

Recall that in QM, observables are represented by operators. Even the humble $$x$$ in classical mechanics is represented by an operator $$\hat{x}$$ in quantum mechanics. When acting on a wavefunction, the effect of the $$\hat{x}$$ operator is to multiply the wavefunction by $$x$$:

$$\hat{x}[\psi(x)] = x\cdot\psi(x)$$

Now, that's kind of boring, but it doesn't stop it from being an operator. Formally, an operator is just something which takes a wavefunction (aka a state, or a ket) and transforms it into another wavefunction. The action of $$\hat{x}$$ did exactly that: it transforms $$\psi$$ into $$x\psi$$.

Generally speaking, the same is true of $$1/r$$. The only difference is that your wavefunction now is expressed as a function of three spherical coordinates $$(r, \theta, \phi)$$. In the equation you wrote, there aren't any hats involved; this is mostly for notational convenience, and also because it barely makes any difference (notice how similar the LHS and RHS are on the above equation). You can imagine that there's a hat there, though, if you're more comfortable with that.

How do I start with a basis function and end up with a number that depends on the distance from the nucleus?

You integrate over all space, which reduces the function into a number. Presumably you know what $$\mu$$ and $$\nu$$ are, because this is your basis set which you have chosen.* These are both functions of $$(r, \theta, \phi)$$. Then:

$$\langle \mu | \frac{1}{r} | \nu\rangle = \int \left(\mu^* \cdot \frac{1}{r} \cdot \nu\right) \,\mathrm{d}V = \int_0^\pi \int_0^{2\pi} \int_0^\infty \left(\mu^* \cdot \frac{1}{r} \cdot \nu\right) r^2 \cos\theta \,\mathrm{d}r\,\mathrm{d}\theta\,\mathrm{d}\phi$$

Actually evaluating this integral with a computer is a totally different challenge, and not one that I'm particularly familiar with, but fundamentally speaking, the above integral tells you exactly what to do.

The distance of what, by the way?

The off-diagonal elements of QM operators ($$\mu \neq \nu$$) are quite hard to interpret classically, so I'd suggest to not think too much about it. At best one may view it as some kind of "interaction" term between the two states $$\mu$$ and $$\nu$$, but that generally does not yield much insight.

For the diagonal elements (i.e. $$\mu = \nu$$), and for a generic operator $$\hat{A}$$, it's possible to interpret the integral $$\langle \mu | \hat{A} | \mu \rangle$$ as representing the average value of the observable $$A$$ when the system is in the state $$\mu$$. (This is a fundamental postulate of QM, often one of the first things that is taught, so it's worth revising that if you don't remember it.) In this case the system is an electron, the state is $$\mu$$, and the observable is $$1/r_m$$ (distance to the $$m$$-th nucleus). So the term above may be interpreted as one over the average distance of an electron in the state $$\mu$$ from the $$m$$-th nucleus.

The basis function is not a point in space . . .

Indeed, and that's the whole point of performing the integration over all space (i.e., all parts of the functions $$\mu$$ and/or $$\nu$$ contribute to the integral).

* It's been a long time since I studied this; please correct me if this is wrong. It doesn't change the point, though. $$\mu$$ and $$\nu$$ should be functions that are fully known.

• It's actually quite easy to think what the off-diagonal terms is once you realise the product of two basis functions describes a charge distribution, so the term is simply the electrostatic interaction between the (typically point charge) nucleus m and the charge distribution described by that product. Dec 24 '21 at 19:35
• @IanBush Hmm, I never thought of it that way before. To me, it still makes more sense for the diagonal term, because then $\mu^*\mu$ may be associated with the charge distribution of orbital $\mu$; but for $\mu^*\nu$, there isn't an object (or an orbital) which has that charge distribution...? Not too sure myself. Dec 24 '21 at 22:10
• Writing the total charge density as expanded in the basis and you get mu nu overlaps. Then write down the line interaction energy between that and the nuclei. Then you get the appropriate Hamiltonian element by taking the derivative of the energy wrt the corresponding density matrix element Dec 24 '21 at 22:21