Stuck deriving the Hylleraas variational method

In summary, the conversation is about a person trying to re-derive Hylleraas' calculation on the ground state of Helium from 1929. The process involves parametricizing the wavefunction in terms of electron-nuclear and electron-electron distances, which leads to a self-adjoint equation and a variational problem. The person is stuck at one point and is seeking help to understand how equation (5) in the original paper leads to equation (6). They have tried various sources but have not found a detailed derivation. Another person suggests using the Euler Lagrange equation from the calculus of variations to solve the problem.
  • #1
alxm
Science Advisor
1,848
9
Hi, first post here.
I've been trying (out of personal interest, not homework) to re-derive http://www.springerlink.com/content/mm61h49j78656107/" relatively famous calculation on the ground state of Helium from 1929. And I'm stuck at one point.

What Hylleraas did, was to parametricize the wavefunction in terms of [tex]r_1,r_2,r_{12}[/tex] - the scalar electron-nuclear and electron-electron distances, which works since it's spherically symmetrical in the ground state. You insert these into the electronic Hamiltonian, do a little work with the chain rule, etc, and arrive at:
[tex]\frac{\partial^2\psi}{\partial r_1^2}+\frac{2}{r_1}\frac{\partial\psi}{\partial r_1}+\frac{\partial^2\psi}{\partial r_2^2}+\frac{2}{r_1}\frac{\partial\psi}{\partial r_2}+2\frac{\partial^2\psi}{\partial r_{12}^2}+\frac{4}{r_{12}}\frac{\partial\psi}{\partial r_{12}}+\frac{r_1^2-r_2^2+r_{12}^2}{r_1r_{12}}\frac{\partial^2\psi}{\partial r_1\partial r_{12}}+\frac{r_2^2-r_1^2+r_{12}^2}{r_1r_{12}}\frac{\partial^2\psi}{\partial r_1\partial r_{12}}
+(\frac{1}{r_1}+\frac{1}{r_2}-\frac{1}{2r_{12}})\psi=-\frac{\lambda}{4}\psi
[/tex]
Equation (5) in the original paper. Where lambda is the energy (Hylleraas used units of R*h=half a Hartree). So far, I'm all good. It's the next step is where I run into trouble. To quote the text (translated):

"This equation is self-adjoint after multiplying with the density function [tex]r_1r_2r_{12}[/tex] and is the Euler equation of a variational problem. This variational problem is:"

[tex]
\int^\infty_0dr_1\int^\infty_0dr_2\int^{r_2+r_1}_{|r_2-r_1|}dr_{12}
\{r_1r_2r_{12}
[(\frac{\partial\psi}{\partial r_1})^2+(\frac{\partial\psi}{\partial r_2})^2+
2(\frac{\partial\psi}{\partial r_{12}})^2]
+r_2(r_1^2-r_2^2+r^2_{12})\frac{\partial\psi}{\partial r_1}\frac{\partial\psi}{\partial r_{12}}+
r_1(r_2^2-r_1^2+r^2_{12})\frac{\partial\psi}{\partial r_2}\frac{\partial\psi}{\partial r_{12}}\\
-[r_{12}(r_1+r_2)-\frac{r_1r_2}{2}]\psi^2\}=\lambda
[/tex]
With the normalization condition:
[tex]
\frac{1}{4}\int^\infty_0dr_1\int^\infty_0dr_2\int^{r_2+r_1}_{|r_2-r_1|}dr_{12}r_1r_2r_{12}\psi^2=1.
[/tex]

Equation (6) in the original paper. So, both sides are multiplied with psi and r1r2r12 and integrated so that the normalization condition can be applied to the right side, giving the energy. (The later problem to is to fit the basis-function parameters to minimize the left side; the integrals can be manually evaluated for the basis used)

My problem here is that I still don't quite follow how he got from (5) to (6). The potential terms are easy, but the first ones..? E.g. How did [tex]\frac{\partial^2\psi}{\partial r_1^2}[/tex] turn into [tex](\frac{\partial\psi}{\partial r_1})^2[/tex] etc? I've looked around quite a bit for a more detailed derivation, but I can't find one. (Closest thing was in Bethe's "QM of one and two-electron atoms" which after a hand-waving reference to using Green's theorem, skipped to the final result)

So. Can anyone give me a clue here? I've stared at this thing so long it's probably something really simple...
 
Last edited by a moderator:
Physics news on Phys.org
  • #3



Hello there, thank you for reaching out. The Hylleraas variational method is a powerful tool for solving the Schrodinger equation for multi-electron atoms. I can understand how you might be stuck in the derivation process, as it can be quite complex. Let me try to explain it in a simpler manner.

The goal of the Hylleraas method is to find a variational wavefunction that will minimize the energy of the system. This wavefunction is parametrized in terms of r1, r2, and r12, as you mentioned. Now, in order to find this wavefunction, we need to solve the Schrodinger equation, which is a second-order differential equation. When we insert the parameters into the Schrodinger equation, we get the equation that you mentioned in your post (equation 5).

Now, in order to make the equation simpler, we can multiply both sides by the density function r1r2r12. This makes the equation self-adjoint, which means that the left and right sides are equal when integrated over all space. This is the step you mentioned where you got confused.

After multiplying by the density function, we can rewrite the equation in a more convenient form, by using the chain rule and some algebraic manipulation. This results in the equation that you mentioned in your post (equation 6).

Now, the variational problem that Hylleraas set up is to find the wavefunction that minimizes the energy. This is done by finding the minimum of the integral on the left side of equation 6, subject to the normalization condition on the right side. This is a variational problem, and the solution to this problem will give us the optimal wavefunction.

I hope this helps to clarify the derivation process for you. It is a complex method, but with patience and persistence, I am sure you will be able to understand it fully. Good luck with your research!
 

FAQ: Stuck deriving the Hylleraas variational method

1. What is the Hylleraas variational method?

The Hylleraas variational method is a mathematical technique used to approximate the ground state energy of atoms or molecules. It is based on the variational principle, which states that the true ground state energy is always lower than any other possible energy state.

2. How does the Hylleraas variational method work?

The method involves choosing a trial wavefunction, which is a mathematical function that describes the possible energy states of the system. This trial wavefunction is then used to calculate an upper bound for the ground state energy using the variational principle. The closer the trial wavefunction is to the true ground state, the more accurate the calculated energy will be.

3. What makes the Hylleraas variational method different from other methods?

The Hylleraas variational method is unique in that it allows for the inclusion of electronic correlations, which are interactions between electrons that can significantly affect the ground state energy of a system. This makes it a more accurate method for calculating ground state energies compared to other methods that do not account for these correlations.

4. What are the limitations of the Hylleraas variational method?

One limitation of the method is that it can be computationally intensive, especially for larger systems, as it requires the evaluation of multiple integrals. Additionally, the accuracy of the calculated energy is highly dependent on the choice of trial wavefunction, so finding the most accurate one can be challenging.

5. How is the accuracy of the Hylleraas variational method determined?

The accuracy of the method is typically determined by comparing the calculated energy with experimental results or with more accurate theoretical calculations. The closer the calculated energy is to the true ground state energy, the more accurate the method is considered to be.

Back
Top