CSIR-NET Quantum Chemistry — Approximation Methods

By Aniket Bhardwaj · 2 October 2026 · CSIR-NET Chemistry

The Schrödinger equation can be solved exactly only for a handful of idealised systems — the particle in a box, the rigid rotor, the harmonic oscillator, and the hydrogen atom. Every real multi-electron molecule needs an approximation, and CSIR-NET Part C regularly tests whether you understand what each approximation actually does, and what it necessarily gets wrong. This article covers the two foundational approximate methods — the variation method and perturbation theory — and the logic behind the Hartree-Fock self-consistent field procedure built on top of them.

The variation theorem

E_trial = ⟨ψ_trial|Ĥ|ψ_trial⟩ ⁄ ⟨ψ_trial|ψ_trial⟩ ≥ E_{true ground state}

For any normalisable trial wavefunction, the calculated energy can never be lower than the true ground-state energy. This gives a built-in check for any variational calculation: the lower the E_trial you obtain, the closer your trial function is to reality — you can freely try different trial functions and simply keep the one giving the lowest energy.

Worked example 1 — variation method for a particle in a box

Q. A particle is confined to a 1-D box of length L (0 ≤ x ≤ L). Use the simple trial function ψ_trial = x(L − x) (which correctly vanishes at both boundaries) to estimate the ground-state energy, and compare it with the exact result.

Step 1 — the numerator, ⟨ψ|Ĥ|ψ⟩. With Ĥ = −(ℏ²/2m)(d²/dx²) and ψ″ = −2 (a constant, since ψ = Lx − x²):

⟨ψ|Ĥ|ψ⟩ = (ℏ²/m) ∫₀ᴸ (Lx − x²) dx = (ℏ²/m) × (L³/6) = ℏ²L³/(6m)

Step 2 — the normalisation integral, ⟨ψ|ψ⟩.

⟨ψ|ψ⟩ = ∫₀ᴸ (Lx − x²)² dx = L⁵/3 − L⁵/2 + L⁵/5 = L⁵/30

Step 3 — divide.

E_trial = [ℏ²L³/(6m)] ÷ [L⁵/30] = 30ℏ²/(6mL²) = 5ℏ²/(mL²)

Compare with the exact answer: E₁(exact) = π²ℏ²/(2mL²) ≈ 4.935ℏ²/(mL²). The trial estimate, 5.000ℏ²/(mL²), sits just above the true value — exactly as the variation theorem guarantees — and is within about 1.3% of it, which is a remarkably good result for such a simple polynomial guess.

First-order perturbation theory

Perturbation theory is used when the true Hamiltonian, Ĥ = Ĥ⁰ + Ĥ′, is a small modification (Ĥ′) of a Hamiltonian Ĥ⁰ you can already solve exactly. The first-order energy correction is calculated using the unperturbed wavefunctions — no new wavefunction needs to be found at this order.

E_n⁽¹⁾ = ⟨ψ_n⁽⁰⁾|Ĥ′|ψ_n⁽⁰⁾⟩

Worked example 2 — a tilted particle-in-a-box potential

Q. A particle in a box of length L (unperturbed ground state ψ₁⁽⁰⁾ = √(2/L) sin(πx/L)) experiences a small linear perturbation, Ĥ′ = λx. Find the first-order energy correction.

E₁⁽¹⁾ = λ ∫₀ᴸ x · (2/L) sin²(πx/L) dx

Because sin²(πx/L) is symmetric about the box's midpoint x = L/2, the average value of x weighted by this probability density is exactly L/2 — so the integral simplifies to λ(2/L) × (L/2) × ∫₀ᴸ sin²(πx/L) dx / L ... more directly: ∫₀ᴸ x sin²(πx/L) dx = L²/4, giving:

E₁⁽¹⁾ = λ(2/L)(L²/4) = λL/2

This is a clean, intuitive result: the energy shifts upward by exactly λ times the value the perturbing potential takes at the box's centre — sensible, because the ground-state probability density is symmetric about that point, so on average the particle "feels" the potential at x = L/2.

The Hartree-Fock (SCF) idea, qualitatively

For a real multi-electron atom or molecule, the electron-electron repulsion term in the true Hamiltonian couples every electron's coordinates to every other electron's, making an exact solution impossible. The Hartree-Fock method sidesteps this by replacing the true, instantaneous repulsion with an averaged field: each electron moves in the mean field created by all the others, described as a single Slater determinant of one-electron orbitals (built as linear combinations of basis functions — the LCAO approach). Because each orbital depends on the average field of the others, and that field depends on the orbitals, the equations must be solved iteratively: guess orbitals → compute the field → solve for new orbitals → repeat until the orbitals stop changing (self-consistency).

What Hartree-Fock necessarily misses: because electrons only feel an averaged repulsion, HF systematically ignores the instantaneous correlation between electrons dodging each other in real time. The energy difference between the true energy and the Hartree-Fock limit is called the correlation energy, and it is exactly the quantity that post-Hartree-Fock methods (configuration interaction, Møller-Plesset perturbation theory, coupled-cluster theory) are built to recover.

Common mistakes that cost marks

Using an unnormalised trial function directly in the variation theorem. If ψ_trial is not already normalised, you must divide by ⟨ψ|ψ⟩ as shown in worked example 1 — skipping this step gives a numerically meaningless energy.
Applying the variation theorem to claim a bound on an excited-state energy directly. The simple inequality E_trial ≥ E_true only guarantees a bound on the ground state unless the trial function is constructed orthogonal to all lower states.
Forgetting that first-order perturbation theory requires Ĥ′ to be genuinely small. The method is a series expansion; a "perturbation" comparable in size to Ĥ⁰ makes the first-order term unreliable.
Confusing Hartree-Fock's mean-field approximation with an exact solution. HF is variational (it obeys E_HF ≥ E_true) but it is not exact — the gap is the correlation energy, not zero.

Method comparison table

Method	What it approximates	Main limitation
Variation method	Ground-state energy from a chosen trial function	Result quality depends entirely on the trial function's flexibility
Perturbation theory	Energy correction from a small Hamiltonian modification	Series only converges usefully when the perturbation is genuinely small
Hartree-Fock / SCF	Multi-electron wavefunction as a single Slater determinant, mean-field repulsion	Misses instantaneous electron correlation entirely

Working with quantum-derived formulas such as de Broglie wavelength or photon energy alongside an approximation-methods question? The Scientific Calculator and the modern-physics tools handle the arithmetic so you can focus on the derivation.

Open the ABC Chemistry Calculator Suite →

Preparing for CSIR-NET, IIT-JAM, GATE or CUET-PG? ABC Chemistry runs dedicated competitive-exam batches at the coaching centre and online across India — details at abcchemistry.in.