Welcome to Quantum Mechanics! If you come from a biology background, the word "Perturbation" might sound intimidating. But think of it this way: In biology, you have a "Control Group" (perfect, untouched) and a "Test Group" (where you change one small thing).
Perturbation theory is exactly that. We take a system we know perfectly (the Control) and add a small disturbance (the Test) to see what happens to the energy.
Since you are comfortable with biology, the best way to understand Perturbation is not through complex physics, but through the concept of Homeostasis.
Imagine a perfectly healthy human body in a comfortable room. The heart rate is steady, the temperature is 37°C, and everything is balanced. In physics, we call this the Unperturbed System.
Now, imagine a small mosquito bites that person.The body doesn't fall apart. It doesn't completely reinvent how it works. Instead, it stays 99% the same, but makes a tiny adjustment (a small bump, a little itch, a slight immune response).
Key Insight: Perturbation theory is the art of saying: "I don't need to recalculate the whole human body just for one mosquito bite. I only need to calculate the itch."
In the real world (and in quantum mechanics), exact solutions are rare.We can solve the math for a Hydrogen atom (1 electron, 1 proton) perfectly. But as soon as we add a second electron (Helium), the math becomes impossible to solve exactly because the two electrons repel each other.
So, we cheat slightly:
We express this "Mosquito Bite" logic using a simple equation:
$$H = H^0 + \lambda H'$$
If the disturbance is huge (like getting hit by a truck instead of a mosquito bite), Perturbation Theory fails. It only works when the change is small enough that the original system is still recognizable.
Imagine a pendulum swinging in a vacuum. It is perfect. We can calculate its energy easily. This is our Unperturbed System.
Now, imagine a gentle breeze blowing on it. The motion changes slightly. The math for the "breeze" is hard to solve directly. So, instead of solving it from scratch, we say:
"The new energy is just the old energy plus a tiny correction."
Mathematically, we write the Total Hamiltonian (Energy Operator) H as:
$$H = H^0 + \lambda H'$$
We want to solve the Schrödinger Equation: Hψ=Eψ.
Since the disturbance is small, we assume the solution (Wavefunction ψ and Energy E) can be written as a series of corrections, just like adding small details to a sketch.
Energy Expansion:
$$E_n = E_n^{(0)} + \lambda E_n^{(1)} + \lambda^2 E_n^{(2)} + \dots$$
Wavefunction Expansion:
$$\psi_n = \psi_n^{(0)} + \lambda \psi_n^{(1)} + \lambda^2 \psi_n^{(2)} + \dots$$
Tutor Note: Don't panic! The superscript (0) means "original", (1) means "first correction", and (2) means "second correction". It is not an exponent!
We plug our expansions into the main equation Hψ=Eψ.
$$(H^0 + \lambda H') (\psi_n^{(0)} + \lambda \psi_n^{(1)} + \dots) = (E_n^{(0)} + \lambda E_n^{(1)} + \dots) (\psi_n^{(0)} + \lambda \psi_n^{(1)} + \dots)$$
Now, we act like we are doing algebra. We multiply everything out and group the terms based on the power of λ.
These are the terms with no λ.
$$H^0 \psi_n^{(0)} = E_n^{(0)} \psi_n^{(0)}$$
This is just the original equation. We already know the solution to this! We move on.
We collect all terms that have exactly one λ.
$$H^0 \psi_n^{(1)} + H' \psi_n^{(0)} = E_n^{(0)} \psi_n^{(1)} + E_n^{(1)} \psi_n^{(0)}$$
Goal: Find En(1).
1. Rearrange to put unknown terms on the left:
$$(H^0 - E_n^{(0)}) \psi_n^{(1)} = (E_n^{(1)} - H') \psi_n^{(0)}$$
2. Multiply by the conjugate wavefunction ⟨ψn(0)∣ and integrate. This is a standard quantum trick.
$$\langle \psi_n^{(0)} | H^0 - E_n^{(0)} | \psi_n^{(1)} \rangle = \langle \psi_n^{(0)} | E_n^{(1)} - H' | \psi_n^{(0)} \rangle$$
The Magic Trick: The Left Hand Side becomes ZERO. Why? Because H0 acting on its own eigenstate gives En(0), so we get (En(0) − En(0)) = 0.
3. Now we just look at the Right Hand Side:
$$0 = E_n^{(1)} \langle \psi_n^{(0)} | \psi_n^{(0)} \rangle - \langle \psi_n^{(0)} | H' | \psi_n^{(0)} \rangle$$
Since ⟨ψn(0)∣ψn(0)⟩ = 1 (Normalization rule), we get our final formula:
$$E_n^{(1)} = \langle \psi_n^{(0)} | H' | \psi_n^{(0)} \rangle$$
Meaning: The first-order energy correction is simply the average value of the perturbation.
We collect all terms with λ2. This is for when we need more precision.
$$(H^0 - E_n^{(0)}) \psi_n^{(2)} = (E_n^{(1)} - H') \psi_n^{(1)} + E_n^{(2)} \psi_n^{(0)}$$
Using similar integration tricks (multiplying by ⟨ψn(0)∣), the left side vanishes again. We eventually derive:
$$E_n^{(2)} = \sum_{m \neq n} \frac{|\langle \psi_m^{(0)} | H' | \psi_n^{(0)} \rangle|^2}{E_n^{(0)} - E_m^{(0)}}$$
Breakdown for Biology Students:
Here is the cheat sheet for your exam:
You asked for the full details, and that is a great way to learn. In the previous section, we used a shortcut ("similarly"). Now, we will walk through the Second Order derivation step-by-step, exactly like an algebra problem. This requires two main parts:
We start by collecting all the terms from our original expansion that had $\lambda^2$ attached to them.
$$H^0 \psi_n^{(2)} + H' \psi_n^{(1)} = E_n^{(0)} \psi_n^{(2)} + E_n^{(1)} \psi_n^{(1)} + E_n^{(2)} \psi_n^{(0)}$$
We want to find $E_n^{(2)}$. To do this, we use the "Projection" trick again. We multiply the entire equation by the conjugate of the original unperturbed wavefunction, $\langle \psi_n^{(0)} |$, and integrate.
Let's look at each term one by one after multiplying by $\langle \psi_n^{(0)} |$:
Putting it together:
$$E_n^{(0)} \langle \psi_n^{(0)} | \psi_n^{(2)} \rangle + \langle \psi_n^{(0)} | H' | \psi_n^{(1)} \rangle = E_n^{(0)} \langle \psi_n^{(0)} | \psi_n^{(2)} \rangle + 0 + E_n^{(2)}$$
Notice that the first term on the Left and the first term on the Right are identical. They cancel out!
$$E_n^{(2)} = \langle \psi_n^{(0)} | H' | \psi_n^{(1)} \rangle$$
The Problem: We have a formula for $E_n^{(2)}$, but it contains $\psi_n^{(1)}$ inside it. We don't know what $\psi_n^{(1)}$ is yet! We have to find it.
To find $\psi_n^{(1)}$, we go back to the First Order Equation (the $\lambda^1$ equation):
$$(H^0 - E_n^{(0)}) \psi_n^{(1)} = (E_n^{(1)} - H') \psi_n^{(0)}$$
The "LEGO" Concept: In Quantum Mechanics, we can build any unknown wavefunction by adding up (stacking) the known wavefunctions ($\psi_m^{(0)}$) like LEGO bricks. We write:
$$\psi_n^{(1)} = \sum_{m \neq n} c_m \psi_m^{(0)}$$
We just need to find the values of the coefficients $c_m$.
Substitute this sum into the Left Hand Side of our equation:
$$(H^0 - E_n^{(0)}) \sum_{m \neq n} c_m \psi_m^{(0)} = \sum_{m \neq n} c_m (E_m^{(0)} - E_n^{(0)}) \psi_m^{(0)}$$
Now, multiply this by a different state $\langle \psi_m^{(0)} |$. Because all states are orthogonal (independent), this filters out everything except the coefficient $c_m$.
Result for the coefficient $c_m$:
$$c_m = \frac{\langle \psi_m^{(0)} | H' | \psi_n^{(0)} \rangle}{E_n^{(0)} - E_m^{(0)}}$$
So, the First Order Wavefunction is:
$$\psi_n^{(1)} = \sum_{m \neq n} \frac{\langle \psi_m^{(0)} | H' | \psi_n^{(0)} \rangle}{E_n^{(0)} - E_m^{(0)}} \psi_m^{(0)}$$
Now we take the result from Part B and plug it back into the equation from Part A.
From Part A: $E_n^{(2)} = \langle \psi_n^{(0)} | H' | \psi_n^{(1)} \rangle$
Substitute $\psi_n^{(1)}$:
$$E_n^{(2)} = \langle \psi_n^{(0)} | H' \left( \sum_{m \neq n} \frac{\langle \psi_m^{(0)} | H' | \psi_n^{(0)} \rangle}{E_n^{(0)} - E_m^{(0)}} \psi_m^{(0)} \right)$$
We can pull the summation and the numbers (the denominator) outside the bracket. We are left with two matrix elements multiplying each other:
$$E_n^{(2)} = \sum_{m \neq n} \frac{\langle \psi_n^{(0)} | H' | \psi_m^{(0)} \rangle \langle \psi_m^{(0)} | H' | \psi_n^{(0)} \rangle}{E_n^{(0)} - E_m^{(0)}}$$
Since $\langle \psi_n | H' | \psi_m \rangle$ is just the complex conjugate of the other term, multiplying them gives the square of the absolute value.
$$E_n^{(2)} = \sum_{m \neq n} \frac{|\langle \psi_m^{(0)} | H' | \psi_n^{(0)} \rangle|^2}{E_n^{(0)} - E_m^{(0)}}$$
And there you have it! The second-order energy correction derived from scratch.
Welcome back! Let's dive deep into the Helium atom. You asked for the full Hamiltonian first, so we will start there. Think of the Hamiltonian as the "Total Energy Manager" of the atom. It tracks every bit of kinetic and potential energy in the system.
The Helium atom consists of a nucleus with charge $+Z$ (where $Z=2$) and two electrons (let's call them electron 1 and electron 2).
The full Hamiltonian operator ($\hat{H}$) is the sum of five distinct terms:
$$ \hat{H} = \underbrace{-\frac{1}{2}\nabla_1^2}_{\text{Kinetic E (e1)}} \underbrace{-\frac{Z}{r_1}}_{\text{Attraction (e1)}} \underbrace{-\frac{1}{2}\nabla_2^2}_{\text{Kinetic E (e2)}} \underbrace{-\frac{Z}{r_2}}_{\text{Attraction (e2)}} + \underbrace{\frac{1}{r_{12}}}_{\text{Repulsion}} $$
Tutor Note: In Physics, we use "Atomic Units" to make life easier. We set mass $m=1$, charge $e=1$, and $\hbar=1$. This removes all the messy constants so we can focus on the algebra!
We cannot solve the full Hamiltonian directly because of that last repulsion term. So, we split it:
$$ \hat{H} = \hat{H}^{(0)} + \hat{H}' $$
We want to find the first-order energy correction, which is the average value of the repulsion calculated using the unperturbed wavefunctions.
The ground state wavefunction for two independent electrons in the 1s orbital is:
$$ \psi^{(0)} = \psi_{1s}(r_1) \times \psi_{1s}(r_2) = \frac{Z^3}{\pi} e^{-Z(r_1 + r_2)} $$
The formula for the energy correction is:
$$ E^{(1)} = \langle \psi^{(0)} | \frac{1}{r_{12}} | \psi^{(0)} \rangle = \iint |\psi_{1s}(1)|^2 \frac{1}{r_{12}} |\psi_{1s}(2)|^2 d\tau_1 d\tau_2 $$
This looks scary, but physically it just means: "Sum up the repulsion between every speck of Cloud 1 and every speck of Cloud 2."
Instead of doing raw calculus on $1/r_{12}$, we use a trick from electrostatics. We treat Electron 1 as a "smear" of charge and calculate the potential it creates.
The Potential Integral ($V$): We fix electron 1 at radius $r_1$ and integrate over electron 2. The mathematical rule is:
This gives us the potential $V(r_1)$ seen by the first electron:
$$ V(r_1) = \frac{1}{r_1} - \left( \frac{1}{r_1} + Z \right) e^{-2Zr_1} $$
Now we calculate the energy of Electron 1 interacting with this potential $V(r_1)$.
$$ E^{(1)} = \int_0^\infty V(r_1) \cdot \underbrace{4\pi r_1^2 |\psi_{1s}(r_1)|^2}_{\text{Charge of Elec 1}} dr_1 $$
Substituting the density $|\psi_{1s}|^2 = \frac{Z^3}{\pi}e^{-2Zr_1}$:
$$ E^{(1)} = 4Z^3 \int_0^\infty \left[ \frac{1}{r_1} - \left(\frac{1}{r_1} + Z\right)e^{-2Zr_1} \right] e^{-2Zr_1} r_1^2 dr_1 $$
We split this into two integrals:
Now we simply subtract the two results and multiply by the constants in front:
$$ E^{(1)} = 4Z^3 \left[ \frac{1}{4Z^2} - \frac{3}{32Z^2} \right] $$
Find a common denominator (which is 32):
$$ E^{(1)} = 4Z^3 \left[ \frac{8}{32Z^2} - \frac{3}{32Z^2} \right] = 4Z^3 \left( \frac{5}{32Z^2} \right) $$
$$ E^{(1)} = \frac{20Z^3}{32Z^2} $$
$$ \mathbf{E^{(1)} = \frac{5}{8}Z} $$
The first-order perturbation theory predicts that the electron-electron repulsion adds an energy of $5/8 Z$ Hartrees to the system.
This document contains the rigorous algebraic steps for Time-Dependent Perturbation Theory, the Stark Effect, and the Zeeman Effect, designed for M.Sc. Chemistry students.
Goal: Calculate the probability of a system jumping from an initial state ($i$) to a final state ($f$) due to a time-varying force $H'(t)$.
We start with the Time-Dependent Schrödinger Equation (TDSE):
$$i\hbar \frac{\partial \Psi(t)}{\partial t} = \hat{H}(t) \Psi(t)$$Where $\hat{H}(t) = \hat{H}^0 + \hat{H}'(t)$. Since the unperturbed states $\psi_n$ form a complete set, the actual wavefunction $\Psi(t)$ can be written as a changing mixture:
$$\Psi(t) = \sum_n c_n(t) \cdot \psi_n \cdot e^{-i E_n t / \hbar}$$Substitute $\Psi(t)$ into the TDSE. On the Left Side, we use the product rule:
$$i\hbar \sum_n \left[ \dot{c}_n \psi_n e^{-i E_n t / \hbar} + c_n \psi_n \left( \frac{-i E_n}{\hbar} \right) e^{-i E_n t / \hbar} \right]$$On the Right Side:
$$(\hat{H}^0 + \hat{H}') \sum_n c_n \psi_n e^{-i E_n t / \hbar} = \sum_n c_n E_n \psi_n e^{-i E_n t / \hbar} + \sum_n c_n \hat{H}' \psi_n e^{-i E_n t / \hbar}$$Cancellation Trick: Notice that the term $\sum c_n E_n \psi_n \dots$ appears on both sides. They cancel out perfectly.
We want to find $\dot{c}_f$. We multiply by $\psi_f^*$ and integrate over all space. Using the property of orthonormality ($\int \psi_f^* \psi_n = \delta_{fn}$), the sum collapses to a single term:
$$\dot{c}_f(t) = -\frac{i}{\hbar} \sum_n c_n(t) H'_{fn} e^{i (E_f - E_n) t / \hbar}$$At $t=0$, the system is 100% in the initial state ($c_i \approx 1$, others $\approx 0$). Integrating with respect to time gives the final coefficient:
$$c_f(t) = -\frac{i}{\hbar} \int_0^t \langle f | \hat{H}'(t') | i \rangle e^{i \omega_{fi} t'} dt'$$Goal: Mathematically prove why the ground state shift is zero and how the excited state splits under an Electric Field ($\mathcal{E}$).
The first-order correction is given by:
$$E^{(1)} = e\mathcal{E} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} |\psi_{100}|^2 z \, dx dy dz$$For $n=2$, we have degenerate states ($2s, 2p_x, 2p_y, 2p_z$). The integral $\langle 2s | z | 2p_z \rangle$ involves Even × Odd × Odd = Even, so it does not vanish.
Solving the perturbation matrix determinant:
$$\begin{vmatrix} -E & -3e\mathcal{E}a_0 \\ -3e\mathcal{E}a_0 & -E \end{vmatrix} = 0$$This yields the energy splitting:
$$E = \pm 3e\mathcal{E}a_0$$Goal: Calculate the energy splitting in a magnetic field ($B$) due to angular momentum.
The interaction Hamiltonian includes Orbital ($\vec{L}$) and Spin ($\vec{S}$) angular momentum:
$$\hat{H}' = \frac{e B}{2m} (L_z + 2S_z)$$Since Hydrogen wavefunctions are already eigenfunctions of $L_z$ and $S_z$, we can simply apply the operators directly to the ket vector:
Substituting these values into the energy equation:
$$E^{(1)} = \frac{e B}{2m} (m_l \hbar + 2 m_s \hbar) \langle \psi | \psi \rangle$$Using the Bohr Magneton constant ($\mu_B = \frac{e \hbar}{2m}$):
$$E^{(1)} = \mu_B B (m_l + 2m_s)$$Physical Interpretation: If $B=0$, all $m_l$ states have the same energy. If $B \neq 0$, the $+m_l$ states go up in energy and $-m_l$ states go down, splitting the spectral lines.
When an atom is placed in an Electric Field, the field acts as the perturbation \( H' \).
This math justifies why certain bonds form and others do not.
This explains how MRI machines and UV-Vis spectrophotometers work. The perturbation is no longer static; it is a light wave oscillating with time: \( H'(t) \).
Instead of static coefficients, the mixing coefficients change with time \( c_b(t) \). We want to find the probability that an electron jumps from state \( a \) to state \( b \).
Result: Fermi's Golden Rule
After integrating the time-dependent Schrödinger equation, we get the rate of transition:
$$ Rate_{a \to b} = \frac{2\pi}{\hbar} \left| \int \psi_b^* H' \psi_a d\tau \right|^2 \times \delta(E_b - E_a - \hbar\omega) $$This derivation connects Quantum Mechanics to Spectroscopy. We will start from the Time-Dependent Schrödinger Equation and derive the formula that predicts the rate of absorption or emission of light.
We start with the fundamental equation of how quantum states change with time. Our Hamiltonian (Total Energy) is split into a static part and a time-dependent perturbation (like a light wave).
$$ i\hbar \frac{\partial \Psi}{\partial t} = (\hat{H}^0 + \hat{H}'(t)) \Psi $$Since the system is changing, we define the total wavefunction \( \Psi(t) \) as a mixture of all possible states, where the mixing coefficients \( c_n(t) \) change over time.
We substitute the sum into the Schrödinger equation. Using the product rule for derivatives on the left side, and cancelling out the static terms (\( H^0 \)) on both sides, we are left with:
Note: \( \dot{c}_n \) represents the rate of change of the coefficient over time.
We want to know the transition to a specific state b. We multiply by \( \psi_b^* \) and integrate. Due to orthogonality, the sum on the left collapses to a single term.
To make this solvable, we assume that at time \( t=0 \), the system is 100% in the ground state a. Therefore, \( c_a \approx 1 \) and all others are 0.
Where \( \omega_{ba} \) is the frequency gap between states.
We define the perturbation as an oscillating light wave: \( \hat{H}'(t) = \hat{V} \cos(\omega t) \). Using Euler's formula and keeping only the resonant term (Absorption):
We integrate from time 0 to \( t \) to find \( c_b(t) \). This yields a complex exponential function.
The probability is the square modulus \( P = |c_b|^2 \). Using trigonometric identities, this converts into the famous Sinc Function form:
We want the Rate (\( W = dP/dt \)). As time \( t \to \infty \), the Sinc function mathematically becomes a Dirac Delta Function, which represents strict Conservation of Energy.
A guide designed for students approaching Quantum Mechanics for the first time.
Imagine you are blindfolded in a valley and you want to find the absolute deepest point (this represents the Ground State Energy). You cannot see the bottom, but you can walk around.
We don't know the true Wavefunction ($\psi_{true}$) or the true Energy ($E_{true}$). So, we guess a "Trial Wavefunction" ($\phi_{trial}$). The energy calculated from this guess ($E_{trial}$) will always be greater than or equal to the true energy.
$$E_{trial} \ge E_{true}$$To calculate the energy of our guess, we use the following formula:
$$E = \frac{\int \phi^* \hat{H} \phi \, d\tau}{\int \phi^* \phi \, d\tau}$$Where:
Let's apply this mathematically. Don't worry about the calculus; we will use standard integration rules.
The System: A particle in a box of length $L$ (from $x=0$ to $x=L$).
The True Answer: We know the true ground state is a Sine wave.
The Guess (Trial): Let's guess a simple parabola because it is easier to integrate than trigonometry.
For a free particle inside the box, potential energy $V=0$.
$$\hat{H} = -\frac{h^2}{8m\pi^2} \frac{d^2}{dx^2}$$First, let's differentiate our guess $\phi = Lx - x^2$.
So:
$$\hat{H}\phi = -\frac{h^2}{8m\pi^2} (-2) = \frac{2h^2}{8m\pi^2} = \frac{h^2}{4m\pi^2}$$We integrate our guess multiplied by the result from Step 2.
$$Numerator = \int_{0}^{L} [x(L-x)] \cdot \left[ \frac{h^2}{4m\pi^2} \right] dx$$Pull the constants out and integrate:
$$= \frac{h^2}{4m\pi^2} \int_{0}^{L} (Lx - x^2) dx$$ $$= \frac{h^2}{4m\pi^2} \left[ \frac{Lx^2}{2} - \frac{x^3}{3} \right]_{0}^{L}$$Substitute limit $L$:
$$= \frac{h^2}{4m\pi^2} \left( \frac{L^3}{2} - \frac{L^3}{3} \right) = \frac{h^2 L^3}{24m\pi^2}$$Integrate term by term:
$$= \left[ \frac{L^2 x^3}{3} - \frac{2L x^4}{4} + \frac{x^5}{5} \right]_{0}^{L} = \frac{L^5}{30}$$Divide Numerator by Denominator:
$$E_{trial} = \frac{\frac{h^2 L^3}{24m\pi^2}}{\frac{L^5}{30}} = \frac{30 h^2 L^3}{24 m \pi^2 L^5} = \frac{1.25 h^2}{m \pi^2 L^2}$$Conclusion: Our trial energy is slightly higher than the true energy, which confirms the Variational Principle works!
This guide breaks down the complex quantum mechanics of the Helium atom into logical steps. We use Atomic Units (where mass of electron $m=1$, charge $e=1$) to simplify the math.
Figure 1: The Helium System ($Z=2$) with two electrons interacting.
We need to define the energy of the system. The total energy operator ($\hat{H}$) is the sum of Kinetic Energy (movement) and Potential Energy (attractions/repulsions).
If we ignore the difficult repulsion term ($1/r_{12}$), Helium looks like two independent Hydrogen atoms. We assume both electrons are in a 1s orbital.
Our Trial Wavefunction ($\phi$):
$$ \phi(r_1, r_2) = \psi_{1s}(r_1) \times \psi_{1s}(r_2) = \frac{\zeta^3}{\pi} e^{-\zeta(r_1 + r_2)} $$We calculate the total energy in three parts.
For one electron in a 1s orbital, the average Kinetic Energy is $\zeta^2/2$ and Potential Energy is $-Z\zeta$.
For two electrons combined:
$$ E_{two\_electrons} = 2 \times \left( \frac{\zeta^2}{2} - Z\zeta \right) = \zeta^2 - 2Z\zeta $$We calculate the average repulsion $\langle 1/r_{12} \rangle$. This involves a complex integral, but the standard result for 1s orbitals is:
$$ E_{repulsion} = \frac{5}{8}\zeta $$Add Step A and Step B together:
$$ E(\zeta) = (\zeta^2 - 2Z\zeta) + \frac{5}{8}\zeta $$Group the terms with $\zeta$:
According to the Variational Principle, we must find the value of $\zeta$ that gives the lowest possible energy. We differentiate $E$ with respect to $\zeta$ and set it to zero.
$$ \frac{dE}{d\zeta} = \frac{d}{d\zeta} \left( \zeta^2 - 2Z\zeta + \frac{5}{8}\zeta \right) = 0 $$Using the Power Rule ($dx^2/dx = 2x$):
$$ 2\zeta - \left(2Z - \frac{5}{8}\right) = 0 $$Solving for $\zeta$:
$$ \zeta = Z - \frac{5}{16} $$Substitute the optimal $\zeta$ back into the energy equation ($E = \zeta^2 - 2\zeta \cdot \zeta = -\zeta^2$):
$$ E_{min} = -\left( Z - \frac{5}{16} \right)^2 $$Plugging in $Z=2$:
$$ E_{min} = -(1.6875)^2 = -2.8476 \text{ Hartrees} $$Converting to Electron Volts (1 Hartree ≈ 27.211 eV):
$$ E = -2.8476 \times 27.211 = -77.48 \text{ eV} $$| Method | Energy (eV) | Accuracy |
|---|---|---|
| Experimental Value | -79.00 eV | Reference |
| First Order Perturbation | -74.80 eV | ~5.3% Error |
| Variational Method (Screening) | -77.48 eV | ~1.9% Error |
Conclusion: By optimizing the "screening" ($\zeta$), the Variational Method gets us much closer to the true experimental value.
Time Remaining: 900s