Orthogonality of Sturm-Liouville eigenfunctions Created 2026-10-06 Updated 2026-10-07
For a regular self-adjoint Sturm-Liouville problem with positive weight and self-adjoint boundary conditions, eigenfunctions with distinct eigenvalues are orthogonal in the weighted inner product . Integration by parts proves this by equating and .
For , separated normal modes decay at rates . Clamping at zero and vanishing moment and shear at give . A mode is
The fourth-order operator is self-adjoint: two integration by parts operations give under these boundary conditions. Its zero eigenspace is trivial and its eigenfunctions form an orthogonal complete basis. Thus a nonzero projection on the first mode gives long-time decay proportional to , where . If this projection vanishes, the lowest excited mode replaces it.
For an infinite square well of width , the normalized trial obeys both Dirichlet boundary conditions and has no interior nodes of a wavefunction. Its kinetic energy is
The ratio to the exact ground-state energy is . Its expansion in the normalized sine eigenfunctions has coefficients for odd and zero for even . This supplies a concrete check of the quantum variational principle and shows that a small energy error can coexist with excited-state admixtures.
Choose the growth direction as and the well as . In the envelope function approximation, write a subband state as , where is the band-edge cell-periodic Bloch factor. The infinite square well has
The envelope has the same shape for electron and hole states: changing their effective masses changes the eigenenergies, but not these infinite-well eigenfunctions. With in-plane polarization, the leading interband momentum matrix element acts on the Bloch factors, while the slow envelopes give
orthonormality therefore yields the infinite-well interband selection rule . A nonzero band-level matrix element is still required; the overlap rule does not replace polarization selection rules. Finite barriers can make electron and hole envelopes different, and interface asymmetry or band mixing can relax this ideal same-index rule.
Use a three-dimensional electron gas with a uniform positive background, the jellium model. The Hartree-Fock method minimizes the expectation of the interacting Hamiltonian over normalized Slater determinants. Antisymmetry is exact within that trial class, but correlations requiring a superposition of determinants are omitted. For orthonormal occupied spin orbitals , the energy is
where is the direct Coulomb integral and the exchange integral; the spin overlap makes exchange vanish between opposite-spin orbitals. The self-interaction terms and cancel. Variation gives a self-consistent Fock operator ; exchange is nonlocal. Its eigenfunctions define occupied orbitals of a stationary determinant, and the determinant must be compared with competing self-consistent solutions to select a variational ground state.
The eigenvalues are orbital Lagrange multipliers. Their sum counts each interaction twice:
Frozen-orbital removal energies may be estimated using Koopmans' theorem, but measured quasiparticle energies can contain relaxation and correlations absent from this approximation. Unoccupied orbital energies are not automatically an exact excitation spectrum. In particular, unscreened exchange in a uniform three-dimensional gas produces an unphysical logarithmic slope singularity at its Fermi surface; a one-determinant mean field does not reproduce the interacting Fermi liquid's screened response.
For a spatially uniform plane-wave determinant, electron-background and background-background terms cancel the uniform direct Coulomb contribution. Retain the background when making this cancellation: an isolated infinite electron charge distribution would have a divergent direct energy. The surviving contributions per particle are kinetic energy and negative same-spin exchange energy. Let define the Dimensionless Wigner-Seitz radius, with , and . Then
For the unpolarized state this is approximately . The kinetic contribution penalizes spin polarization because it enlarges an occupied Fermi sphere; exchange favors polarization because equal-spin occupation gains more exchange energy. As density falls, the exchange scale overtakes the kinetic scale. The uniform fully polarized and unpolarized energies cross at
This is a variational comparison within the homogeneous Hartree-Fock electron gas, not an exact physical ferromagnetic phase boundary. It demonstrates how Pauli kinetic energy and exchange compete and why the allowed trial states matter.
Allowing spatial symmetry breaking changes the result: the Overhauser instability lowers the paramagnetic determinant's Hartree-Fock energy through a spin density wave, even in the high-density limit. A calculation restricted to uniform Fermi spheres therefore does not give the unrestricted Hartree-Fock ground state.[ref] At still lower density, localization and a Wigner crystal compete with fluid states. electron-gas correlation energy beyond a single determinant is important in assessing their competition. Quantum Monte Carlo provides a concrete route to comparing correlated fluid and crystal energies; the original stochastic calculation demonstrates why exchange-only energy differences are insufficient.[ref]
Hartree-Fock is a useful variational baseline with exact exchange, not a controlled complete low-density phase diagram. At high density the kinetic scale dominates, but even a small symmetry-breaking energy gain can invalidate a restricted determinant as the exact Hartree-Fock minimizer. At low density neglected correlations, screening and localization can change both transition locations and the ordering of candidate phases. A numerical phase boundary should therefore state dimension, spin constraints, background model and whether correlation and spatially modulated states were included.
Figure 1.
Original homogeneous Hartree-Fock kinetic and exchange energy comparison; the crossing is a restricted variational result
. Details.
Let
The normalized orbital energy eigenstates vanish outside the well and have one-particle energies . Define the normalized symmetric and antisymmetric combinations for by
The normalized spin functions are the singlet state and the triplet states , , . Their exchange energies are respectively and .
For , the three lowest distinct energy levels of the exchange-coupled spin-one-half particles in an infinite square well are
The spatially doubly occupied orbital can only occur with the singlet state. The next orbital sum is , where either exchange symmetry is possible. The next unused sum is , well separated under the small-coupling assumption. The paper's inequality is equivalent, up to a numerical factor, to ; the factor in is essential because the spin operators carry units of angular momentum.
The sign of is not explicitly specified. The table assumes the usual intended positive coupling in the displayed minus-sign convention. For with , the same three sets of eigenfunctions occur, but the two rows at exchange order. At those rows merge into a single energy level of multiplicity four; if three distinct energies are requested exactly at that point, the third is with the single state .
For distinguishable particles, the ordered spatial products are independent states, and neither singlet states nor triplet states are excluded from any orbital occupation. The unperturbed energy eigenvalues are
For , the three lowest distinct energy levels are
where . The first two levels do have zero first-order correction, since their spatial expectation value factorizes as .
The last instruction in the original PDF is false for the third level. Although every diagonal product-state matrix element vanishes, the two orbital states have degenerate energy levels. Degenerate perturbation theory requires diagonalizing the perturbation within their span. For each , its matrix in the ordered basis is
The eigenfunctions are , with
Thus the sixfold unperturbed level splits into two triplets. Its average shift is zero, and the diagonal expectation in an individual ordered product state is zero, but neither fact makes all its energy eigenvalue corrections zero. This is the degenerate orbital exchange splitting from a product perturbation; distinguishability allows both combinations rather than removing their coupling.
For negative , the lowest two rows exchange order and the third unperturbed level is instead the twofold singlet orbital pair at . It splits into with the same shifts , each nondegenerate. At , the three lowest distinct unperturbed energies are , with multiplicities . The first and third have zero first-order shift; the middle splits by , each fourfold. The printed zero-shift assertion therefore fails also in the zero-exchange limit.
For the multiplication operator, . Taking gives equality, so .
An eigenfunction obeys at every . If , its level set is the single point , and a continuous function supported at a single point is zero. If , that level set is empty. For , every continuous function vanishing on is an eigenfunction, and nonzero examples exist, such as . Thus .
If , multiplication by is a bounded inverse of . If , choose with . Every function in the range of vanishes at , so the range omits the constant function one. Hence it is not onto. Therefore
The squared angular momentum operator and one selected component are commuting operators, so a complete set of common eigenfunctions can be chosen: for the component these are the usual angular momentum eigenstates. Distinct components do not have a common complete eigenbasis, because their commutator is another nonzero component. A state is an exception in which all components vanish; noncommutation does not forbid every individual common eigenfunction.
Separation of variables makes the angular factor a spherical harmonic. Applying the Quantum Hamiltonian gives
Consequently is an energy eigenfunction provided the radial factor satisfies
This is the radial Schrödinger equation for a central potential, here with the attractive Coulomb potential. An arbitrary does not give an eigenfunction; regularity, decay and wavefunction normalization select admissible bound quantum states.
Twice differentiation the sine gives
Both Dirichlet boundary conditions hold when is an integer. The zero function at is not an eigenfunction, and negative merely changes an overall sign. Therefore
The corresponding normalized eigenfunction is .
Up to a global phase, is real, positive throughout the interior, has no interior nodes of a wavefunction, and is symmetric about . It is smooth inside the box and has one central maximum. These are the shape and reflection symmetry properties of the lowest particle in a box eigenfunction, ; they make it a plausible variational method trial without making it exact.
The trial is even under reflection about the midpoint. The normalized box eigenfunctions satisfy , so orthogonality between even and odd functions gives for even . If the higher odd contributions are neglected, the weighted mean from the quantum variational principle becomes
Rearranging yields the two-level variational contamination estimate
The required assumption concerns the neglected contributions to the energy as well as their small probabilities. Numerically the two-level estimate is . It is not an exact identity: direct expansion of this parabolic trial gives for odd , hence the exact ratio is . The omitted higher odd states account for the difference.
A symmetry orbital is a molecular orbital when it is also an eigenfunction of the effective one-electron Hamiltonian operator. If an irreducible representation occurs once in a complete invariant model space, a symmetry-invariant Hamiltonian acts as a scalar on that copy by Schur's lemma. Its basis functions are then degenerate eigenfunctions. When the same symmetry occurs in several copies, those copies can mix; one must diagonalize the corresponding block. Merely having a one-dimensional symmetry label is insufficient if other orbitals have that same label.
The unperturbed Hamiltonian satisfies
The supplied eigenfunctions can be chosen as an orthonormal eigenbasis, with .
Linearity and the shared angular eigenvalue give
Using the relative phases in part (c),
Thus each real orbital remains an energy eigenstate with hartree and has definite orbital angular momentum magnitude , but neither is an eigenfunction of . A measurement of gives or , each with probability . Its mean is zero and its variance is ; a zero mean does not mean that a measurement gives zero.
For two eigenfunctions, apply the same Hermitian operator identity across the two states:
Hence . Distinct eigenvalues therefore imply orthogonal eigenfunctions. Within a degenerate eigenspace, orthogonality is a possible basis choice, rather than a consequence for every pair of vectors.
The decisive observation is that is itself an eigenfunction:
On the regular Legendre polynomial domain, integration by parts has endpoint factor , so the boundary form vanishes for regular functions at . The interval alone is not an endpoint condition; regularity is the usual condition selecting the Legendre polynomials. Their eigenvalues are , and the normalized function is
Therefore
Equivalently, the Hermitian operator identity gives , and nondegeneracy leaves only one possible nonzero coefficient. The unnormalized function has squared norm , consistent with this coefficient.