After another long hiatus I’ll be back to more regular blogging: I’m teaching a course entitled *the variational principle in quantum mechanics* for the summer semester here at the ITP Hannover, and I’ll be posting the notes. I hope to post other things along the way as there are many exciting developments I would like to mention…

Please feel free to comment or make suggestions on the lecture notes, I would sincerely appreciate your feedback!

**1. Lecture 1 **

** 1.1. Organisation **

Lectures will take place on fridays at 12:00 until 14:00 (negotiable). The course comprises 16 lectures in english. Assessment will be via a combination of weekly homework exercises and a final examination.

** 1.2. Course outline **

We will cover the following topics.

- Reminder of quantum mechanics and the variational method.
- The helium atom.
- The many body problem.
- Mean field theory and Hartree-Fock theory.
- The Gross-Pitaevskii equation.
- Density functional theory.
- Density functional theory cont.
- Matrix product states.
- Matrix product states cont.
- The density matrix renormalisation group (DMRG).
- The DMRG cont.
- The time-dependent variational principle (TDVP).
- Applications of the TDVP to density functional theory and Hartree-Fock.
- Applications of the TDVP to density functional theory and Hartree-Fock cont.
- QFT and beyond.
- Review.

** 1.3. Lecture notes **

The notes for the lectures will be posted on my blog site: https://tjoresearchnotes.wordpress.com/

** 1.4. A reminder of quantum mechanics **

See handout.

** 1.5. The variational principle **

The variational principle is the basis of a tremendous number of highly successful calculational tools in physics, which is surprising because the basic method is simple enough to cover in an introductory quantum mechanics course. Of course life is never simple, especially for interacting systems of many particles. The trouble is that it is very hard to come up with a variational wavefunction that you can actually calculate with and which bears some resemblence to the actual ground-state. You need to be extremely clever to design reasonable variational wavefunctions, and it often requires deep new insights into the physics of a system to come up with a good one. Some classic examples include the BCS wavefunction in superconductivity, the Gutzwiller wavefunction, and the Laughlin wave function for the fractional quantum Hall effect (ok, the last one is a bit of a cheat, it doesn’t really have any free parameters).

Consider an arbitrary physical system whose hamiltonian is time independent. Let’s initially assume, for simplicity, that the spectrum of is discrete and non degenerate:

with . Although the hamiltonian is assumed known, this isn’t necessarily the case for and . Generally speaking, the variational method is most useful in the case where we cannot diagonalise exactly, although there are many situations where it can be helpful to exploit the method when exact diagonalisation is, in principle, available because it can be much less expensive computationally.

*Ground-state properties*

Choose an arbitrary state from the state space of the system (we don’t assume is normalised). The expectation value of the energy of the system in this state is given by

with equality occurring if and only if is proportional to the eigenvector .

*Proof:* Exercise.

This simple property is the basis for an approximate determination of the ground-state energy : choose a family of *trial states* which depend on a number of variational parameters . Calculate the expectation value

with respect to the parameters . The minimal value so obtained is an approximation to which, thanks to the above, always overestimates the ground-state energy: . The closer gets to the more closely the corresponding trial state should approximate the properties of .

This is the *variational method*.

*The Ritz theorem*

Theorem 1The expectation value is stationary in the vicinity of its discrete eigenvalues

*Proof:* Consider the energy expectation value

as a functional of . Calculate the increment ; stationarity implies that, to first order in the components of :

Set

Now (5) becomes

when we write

The norm of is therefore zero, and must be consequently zero. This means that

Thus the expectation value is stationary if and only if is an eigenvector of , and the stationary values of are the eigenvalues.

Thus the variational method may be generalised to the approximate determination of eigenvalues of . If has several extrema, they give the approximate values of some of its energies.

**2. Application to a simple example **

To illustrate the discussion and to give an idea of how valid the approximations made in the variational method are, we apply the method to the one-dimensional harmonic oscillator

Since the hamiltonian is even we know the ground-state must be represented by an even wavefunction.

The trial wavefunctions we explore here are the *rational wavefunctions*

The normalisation is easily found to be, after a simple integration (exercise)

so that (exercise)

The minimum value of is found to be at (exercise)

where

This result is within a relative factor

of the exact answer .

** 2.1. Caveats **

The variational method suffers from two important drawbacks:

- There is no
*a priori*guarantee whatsoever that the trial state corresponding to the minima of the energy expectation value at has any resemblance to the ground state. - For a generally chosen set of trial states there is typically no way to analytically find the minimum of .

By developing ways to overcome these two crucial drawbacks researchers have built a surprisingly rich and innovative set of methods to apply the variational principle to various physical systems. Roughly speaking, as we’ll see, to overcome the first problem one argues that the variational class is complete, in the sense that a member of the class could, *in principle*, capture the physics of a given system. To overcome the second problem one relaxes the minimisation via some minimisation heuristic such as Newton’s method. The combination of these two strategies has culminated in what must be the most sophisticated application of the variational method in quantum mechanics: *the density matrix renormalisation group*.

great post, love it. I heard variational problem when i was reading Feynman integral

sir how can occur eqation No. 12, 13, 14 &15 in above method

Dear Ganesh Chorage,

The only method I know to evaluate (12) and (13) is via the calculus of residues: see, e.g., http://en.wikipedia.org/wiki/Residue_theorem

The other equations follow by simply differentiating to find the minima and substituting for the optimal value found.

Sincerely,

Tobias

I like what you guys are up too. This type of clever work and exposure!

Keep up the superb works guys I’ve incorporated you guys to our blogroll.

How can I get the equation (5)? I am confused about it. Thank you.

You have to substitute in place of in Eq. (3) and expand to first order. The denominator gives the extra terms in Eq. (5). I hope this helps; sincerely,

Tobias