Calculate Expected Value of Hamiltonian Without Integrals

Calculating the expected value of a Hamiltonian operator is a fundamental task in quantum mechanics. While traditional methods involve integrals over wavefunctions, we can use matrix methods to compute this without explicit integration.

What is a Hamiltonian?

The Hamiltonian (H) is an operator in quantum mechanics that represents the total energy of a system. It's a key component in the Schrödinger equation:

iℏ ∂/∂t |ψ⟩ = H |ψ⟩

Where:

ℏ is the reduced Planck constant
|ψ⟩ is the quantum state
H is the Hamiltonian operator

The Hamiltonian is typically expressed in terms of position and momentum operators, and its eigenvalues correspond to the possible energy levels of the system.

What is the Expected Value of a Hamiltonian?

The expected value of the Hamiltonian, denoted as <H>, represents the average energy of a quantum system in a given state. It's calculated as:

<H> = ∫ ψ*(x) H ψ(x) dx

Where:

ψ(x) is the wavefunction
ψ*(x) is the complex conjugate of the wavefunction
H is the Hamiltonian operator

This integral represents the average energy of the system, weighted by the probability density |ψ(x)|².

Calculation Method Without Integrals

Instead of performing the integral directly, we can use matrix methods when the Hamiltonian is expressed in a basis set. The expected value can be computed as:

<H> = Σ c_i* c_j H_ij

Where:

c_i are the expansion coefficients of the wavefunction in the basis set
H_ij are the matrix elements of the Hamiltonian in the basis set

This approach is particularly useful in computational quantum chemistry where basis sets are commonly used to represent wavefunctions.

Key Assumption

The Hamiltonian must be expressed in a complete basis set for this method to be valid. The basis set should be chosen to accurately represent the system's wavefunction.

Worked Example

Consider a simple two-state system with Hamiltonian matrix:

H = [ 2 1 ] [ 1 3 ]

And a wavefunction represented by coefficients c₁ = 0.6 and c₂ = 0.8 (normalized to 1).

The expected value is calculated as:

<H> = (0.6*0.6*2) + (0.6*0.8*1) + (0.8*0.6*1) + (0.8*0.8*3) = 0.72 + 0.48 + 0.48 + 1.92 = 3.6

This shows the average energy of the system in this state is 3.6 units.

FAQ

Why avoid integrals when calculating expected values?

Integrals can be computationally intensive, especially for complex systems. Matrix methods provide a more efficient approach, particularly when using basis sets in quantum chemistry calculations.

What basis sets are typically used?

Common basis sets include Slater-type orbitals, Gaussian-type orbitals, and plane waves. The choice depends on the system and computational requirements.

How accurate is this method compared to direct integration?

When implemented properly with a complete basis set, the matrix method yields equivalent results to direct integration. The accuracy depends on the quality of the basis set representation.