The hydrogen molecular ion and bonding
Use a nonorthogonal two-orbital variational model to explain one-electron bonding while comparing total molecular energy with a declared dissociation threshold.
On this page
The question for today
What you’ll understand
- Normalize symmetric and antisymmetric combinations with orbital overlap.
- Calculate their variational energies at a fixed geometry.
- Distinguish electronic stabilization from a stable total molecular bond.
The hydrogen molecular ion contains two protons and one electron, making it a clean setting for studying a chemical bond without electron-electron repulsion. Its electron can occupy a coherent state spread over both centers. A simple linear combination of atomic orbitals captures the interference behind bonding and antibonding states, but it is still a restricted variational approximation. To decide whether a geometry is bound, we must also include proton repulsion and compare the total energy with the same separated-fragment energy reference.
One electron feels two attractive centers
At fixed proton separation R, use the electronic Hamiltonian excluding proton-proton repulsion, consistent with the preceding lesson. It contains electron kinetic energy and attraction to each proton. The electronic eigenvalue is then supplemented by positive nuclear repulsion to obtain the nuclear potential curve. At very small R the nuclear repulsion grows strongly; at large R the appropriate dissociation channel is a hydrogen atom plus a proton. Overall center-of-mass and finite-nuclear-mass corrections are omitted at this level.
Choose normalized 1s-like functions and centered on the two nuclei. They overlap, so their inner product S is generally nonzero. Their sum has enhanced amplitude between the nuclei, while their difference has a node on the median plane for equivalent centers. The sign of an orbital amplitude is not an electric charge: it specifies the relative quantum phase needed when amplitudes are added. Density is obtained only after squaring the complete normalized wavefunction.
Overlap changes both normalization and energy
For equivalent centers, symmetric and antisymmetric trial states have normalization factors involving and , respectively. Evaluating the electronic expectation value gives energies when the matrix elements are real and equal-center symmetry applies. Dividing only by two would be incorrect for overlapping functions. The same calculation is a generalized eigenproblem with an overlap matrix, not an ordinary eigenproblem in an orthonormal atomic basis.
The symmetric state is lower than the diagonal electronic expectation when . The overlap contribution means that the sign of an off-diagonal element by itself is insufficient without its convention and normalization. A variational calculation provides an upper bound to the exact ground electronic energy for the chosen Hamiltonian and geometry. Enlarging the trial space can improve that bound, but a two-function result does not give an exact equilibrium bond length, dissociation energy, or complete spectrum.
A bond is an energy comparison across geometry
Enhanced internuclear electron density can favor bonding, but the complete energy balance also includes electronic kinetic energy, electron-nucleus attraction, and nuclear repulsion. A bound minimum requires the total surface to fall below an appropriate separated-fragment threshold over a range of geometries. Comparing one electronic eigenvalue with a bare atomic number while omitting nuclear repulsion can falsely suggest excessive binding. Even a total surface minimum lies below the lowest vibrational state by its zero-point energy.
The elementary molecular-orbital bond-order convention assigns one half to a single electron in a bonding orbital, but that count is a qualitative model label, not a directly measured bond energy. Spectroscopy tests actual energy differences and nuclear motion. For the example below, all matrix elements and the separated threshold are deliberately hypothetical to make the arithmetic transparent. They illustrate the calculation’s logic rather than quote computed or measured constants of the physical hydrogen molecular ion.
Work it through
Electronic lowering is not the full bond energy
At one hypothetical geometry use , , and S = 0.20. Nuclear repulsion is +1 eV, and the stated separated-fragment threshold is −10 eV. Find the two electronic energies and decide whether the lower total energy at this geometry is below threshold.
The symmetric electronic energy is . Its normalization includes constructive overlap.
The antisymmetric electronic energy is . It is higher in this restricted model.
Adding nuclear repulsion gives the lower total energy . This lies 1/6 eV above the specified −10 eV threshold, despite electronic lowering relative to H AA.
The symmetric electronic state is lower, but this geometry does not bind relative to the stated threshold. Determining a bond requires scanning total energy versus separation, not inspecting electronic interference alone.
The living notebook
Explore the model
Interactive lab · waves
Relative phase changes an amplitude sum
Use the two fixed-frequency waves as an interference analogy. Reset, set Phase difference = 0° and Amplitude ratio = 1, and predict doubled wave amplitude. For a separate run Reset, set Phase difference = 180° and Amplitude ratio = 1, and predict cancellation.
- Sum at x = 0
- 0.866
- Resultant amplitude
- 1.732
- Phase difference
- 60°
Two waves share the same frequency and wave number. Playback advances t in sin(2x − t) + a sin(2x − t + phase). Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.
Look forThe lab sums traveling sine waves, not overlapping atomic orbitals. It does not normalize a molecular wavefunction, integrate electronic energy, or predict a bond; the relative-phase idea is the only intended connection.
Pause & reason
Check your understanding
Your turn
Try it yourself
For a different hypothetical geometry take S = 0, , and . Nuclear repulsion is +0.5 eV and the separated threshold is −5 eV. Find total symmetric-state energy and its stabilization below threshold.
Take a moment to make a prediction or work it out on paper. Open a hint when you need one.
Give me a hint
Orthogonal trial functions remove the overlap denominator, but nuclear repulsion still belongs in total energy.
Compare with the solution
The symmetric electronic energy is . Including repulsion gives , which is 0.5 eV below the stated threshold. This geometry is energetically favorable in the hypothetical model. Establishing an equilibrium distance or a nuclear bound level still requires a potential curve and nuclear kinetic energy, not this one point alone.
Keep the idea
Take it with you
- Atomic-orbital overlap enters both normalization and the variational eigenproblem.
- Relative orbital phase controls interference, while probability density is the squared full state.
- Bonding requires a total-energy curve and a consistent dissociation reference, beyond an electronic energy lowering.
Read further
Original sources and deeper explanations to follow at your own pace.
Optional progress, saved in this browser when available. Revisit any lesson at any time.
Lesson not marked complete.