Differences

This shows you the differences between two versions of the page.

--- the_variational_principle [2020/08/21 22:54] – [5.i.4 Dealing with A Discontinuous Ansatz] admin
+++ the_variational_principle [2020/08/24 07:02] (current) – [In Class Activities] admin
@@ Line 1: / Line 1: @@
 ====== 5.i.1 Introduction ======
-This section discusses two related concepts: the //variational principle// and the //variational method// (also called the //Rayleigh-Ritz method//).  The variational principle shows that solving the Time Independent Schrödinger Equation (TISE) is equivalent to the stationarity of the energy functional $E[\psi]$ in the calculus of variations.  From classical mechanics, you should be familiar with the idea that dynamical laws can be formulated as variational principles.  In that case, the trajectories of a classical system are those for which the variation of the action $\delta S$ is zero.  The //variational method// is a method of approximating the energy spectrum and stationary states, particularly the ground state, that is based on the variational principle.
+This section discusses two related concepts: the //variational principle// and the //variational method// (also called the //Rayleigh-Ritz method//).  The variational principle shows that solving the Time Independent Schrödinger Equation (TISE) is equivalent to the stationarity of the energy functional $E[\psi]$ in the calculus of variations.  From classical mechanics, you should be familiar with the idea that dynamical laws can be formulated as variational principles.  In that case, the trajectories of a classical system are those for which the variation of the action $\delta S$ is zero.  The //variational method// is a method of approximating the energy spectrum and stationary states of a Hamiltonian $\hat{H}$, particularly the ground state, that is based on the variational principle.
 There are many potentials for which it is difficult or impossible to calculate the energy spectrum and stationary states exactly.  For this reason, we either have to use numerical methods or employ suitable approximation schemes to solve most physically realistic potentials.  Several approximation methods are covered in PHYS 452, but the variational method, which we discuss here, is perhaps the simplest.  A better name for it would be "guessing", but physicists do not like to admit that they are just guessing in their research papers, hence the more sophisticated sounding name.
-We also often use the word //ansatz// instead of "guess".  When quantum mechanics was being developed in the 1920's and 1930's, the language that most physics papers were published in was German, so it was natural to use the German word "ansatz".  It has several different meanings and is not precisely translatable into English, but it has the same sense as words like "attempt", "first approach", or "beginnings".  As the center of world physics moved westwards around the time of World War II, and the lingua franca of physics became English, and it proved more appealing to retain the word "ansatz" than replace it by the accurate, but less impressive sounding, word "guess".
+We also often use the word //ansatz// instead of "guess".  When quantum mechanics was being developed in the 1920's and 1930's, the language that most physics papers were published in was German, so it was natural to use the German word "ansatz".  It has several different meanings and is not precisely translatable into English, but it has the same sense as words like "attempt", "first approach", or "beginnings".  As the center of world physics moved westwards around the time of World War II, and the lingua franca of physics became English, it proved more appealing to retain the word "ansatz" than replace it by the accurate, but less impressive sounding, word "guess".
 ====== 5.i.2 The Variational Principle ======
@@ Line 119: / Line 119: @@
 and for the potential term we have,
-\[\langle \hat{T} \rangle = \frac{1}{2}m \omega^2 \int_{-\infty}^{+\infty} \mathrm{d}x \, \psi(x)^* x^2 \psi(x).\]
+\[\langle \hat{V} \rangle = \frac{1}{2}m \omega^2 \int_{-\infty}^{+\infty} \mathrm{d}x \, \psi(x)^* x^2 \psi(x).\]
 In an in class activity, you will show that these integrals evaluate to $\langle \hat{T}\rangle = \frac{\hbar^2 b}{2m}$ and $\langle \hat{V} \rangle = \frac{m\omega^2}{8b}$, so
@@ Line 143: / Line 143: @@
 Let's think about how to use the variational method more generally.  If we have a Hamiltonian $\hat{H}$ for which it is too difficult to solve for the eigenstates directly, then can still get an upper bound on the ground state energy by simply guessing any quantum state we like and calculating the expectation value of its energy.  Unfortunately, this is an //uncontrolled// approximation, meaning that we do not have a way of estimating how far above the ground state energy the expectation value is.  However, we can often use what we know about the potential to make an intelligent <del>guess</del>ansatz.
-For example, if the potential is $V(x) = k x^4$ for some positive constant $k$, then this potential looks similar to the Harmonic oscillator potential $V(x) = \frac{1}{2} m\omega x^2$ in the sense that they both have a minimum at $x=0$, they are both even functions, and they both curve upwards on either side of $x=0$.  It would therefore make sense to choose the ground state of the Harmonic oscillator as an ansatz for the potential $V(x) = k x^4$.  We also know that the ground state of any potential has no nodes and, if $V(x)$ is even, then the ground state is an even function.
+For example, if the potential is $V(x) = k x^4$ for some positive constant $k$, then this potential looks similar to the Harmonic oscillator potential $V(x) = \frac{1}{2} m\omega x^2$ in the sense that they both have a minimum at $x=0$, they are both even functions, and they both curve upwards on either side of $x=0$.  It would therefore make sense to choose the Gaussian ansatz $\psi(x) = Ae^{-bx^2}$ for the potential $V(x) = kx^4$ because we know that this is the form of the ground state of the harmonic oscillator.  We also know that the ground state of any potential has no nodes and, if $V(x)$ is even, then the ground state is an even function.
-The other thing we can do is to make our ansatz depend on some parameters $a,b,c,\cdots$, which are initially undetermined.  We did this in the harmonic oscillator example with the parameter $b$.  We can then calculate $\langle H \rangle = f(a,b,c,\cdots)$ as a function of those parameters and minimize $f(a,b,c,\cdots)$ over different choices of the parameters.  The more parameters we use, the closer we will get to the true ground state energy, but the harder it will be to minimize $f$.  In the extreme case where we leave $\psi(x)$ completely undetermined, we would have an infinite number of parameters and would get the exact ground state energy, but the minimization problem is then equivalent to solving $\hat{H}|\psi\rangle = E_{\mathrm{gs}}|\psi\rangle$ for the ground state, so it is of no help in simplifying the problem.  The art of using the variational method is therefore to use both our knowledge of what the ground state looks like for similar Hamiltonians //and// a small number of undetermined parameters so that we are confident that we can get close to the true ground state with a minimization problem that is still tractable.
+The other thing we can do is to make our ansatz depend on several (i.e. more than one) parameters $a,b,c,\cdots$, which are initially undetermined.  We can then calculate $\langle H \rangle = f(a,b,c,\cdots)$ as a function of those parameters and minimize $f(a,b,c,\cdots)$ over different choices of the parameters.  The more parameters we use, the closer we will get to the true ground state energy, but the harder it will be to minimize $f$.  In the extreme case where we leave $\psi(x)$ completely undetermined, we would have an infinite number of parameters and would get the exact ground state energy, but the minimization problem is then equivalent to solving $\hat{H}|\psi\rangle = E_{\mathrm{gs}}|\psi\rangle$ for the ground state, so it is of no help in simplifying the problem.  The art of using the variational method is therefore to use both our knowledge of what the ground state looks like for similar Hamiltonians //and// a small number of undetermined parameters so that we are confident that we can get close to the true ground state with a minimization problem that is still tractable.
 All this <del>guessing</del>ansatzing might seem not very rigorous, but this is actually how physical chemists calculate the ground states of atoms and molecules that are more complicated than the hydrogen atom.  By introducing a large number of undetermined parameters and performing the minimizations numerically, they are able to achieve remarkable agreement with the ground state energies observed in the lab.
@@ Line 273: / Line 273: @@
 The true ground state energy of the harmonic oscillator is $\frac{1}{2} \hbar \omega$.  It should not be a surprise that the upper bound obtained from our ansatz is a bit larger than this, since we know that the ground state is gaussian.
 ====== 5.i.4 Finding Excited States ======
-{{:question-mark.png?nolink&50 |}}====== In Class Activities ======
+Since all of the energy eigenstates are stationary points under variations of the expectation value of the Hamiltonian we can, in principle, find the excited states by the same method.  If you know the exact ground state $|E_1\rangle$ then you can choose an ansatz $|\psi_2\rangle$ that is orthogonal to the ground state $\langle \psi_2 | E_1 \rangle = 0$.  Then, we know that $|\psi_2 \rangle$ is a superposition of excited states $|\psi_2 \rangle = \sum_{j\geq 2}^n c_j |E_j \rangle$, with no ground state component in the superposition.  By the same argument we gave for the ground state, we will have $\langle \hat{H} \rangle \geq E_2$, so by varying the parameters of the ansatz we can obtain a bound on the energy of the first excited state.
+We can proceed for the higher excited states analagously.  If we want an upper bound on $E_n$ and we know the states $|E_1\rangle, |E_2\rangle, \cdots, |E_{n-1}\rangle$ then we can choose an ansatz $|\psi_n\rangle$ that is orthogonal to all of them, and then we will have $\langle \hat{H} \rangle \geq E_n$.
+However, there are two major problems with this method.  First, as we increase $n$, the number of orthogonality conditions we have to impose increases, so it gets harder to choose an ansatz that satisfies all of the orthogonality conditions.  Second, and more seriously, when we want to know the energy $E_n$, we typically do not know the eigenstates $|E_1\rangle, |E_2\rangle, \cdots, |E_{n-1}\rangle$ exactly, but only the approximations to them $|\psi_1\rangle, |\psi_2\rangle, \cdots, |\psi_{n-1}\rangle$ that we have obtained from the variational method.  The errors in these states will accumulate as we calculate the energies of higher and higher excited states.
+For this reason, the variational method is typically only used to estimate the ground state and perhaps the first few excited states, but after that, other approximation methods that we will learn about in PHYS 452 will be more useful.
+{{:question-mark.png?nolink&50 |}}
+====== In Class Activities ======
   - For the harmonic oscillator Hamiltonian
@@ Line 283: / Line 294: @@
   has
   \[\langle \hat{T}\rangle = \frac{\hbar^2 b}{2m}, \quad \mathrm{and} \quad \langle \hat{V} \rangle = \frac{m\omega^2}{8b}.\]
+  - For the same Hamiltonian, consider the ansatz $\psi(x) = \sqrt{b} e^{-b|x|}$.  Show that the potential energy term is
+  \[\langle \hat{V} \rangle = \frac{m\omega^2}{4b^2}.\]