Consider a box of width $L$ and the composed of the following potential $$V(x)=\frac{V_0x(x-L)}{L^2}, x\in[0,L]$$ and $V(x)=\infty$ elsewhere. Using perturbation theory - with a square box as the similar potential - one finds that the perturbation potential for $x\in[0,L]$ is simply the potential describing the sagging box (which makes sense because $V(x)=0$ for normal box).
When calculating the nth state energy shift to first order one also finds that $$\Delta E_n^{(1)}=-V_0\Big({1\over6}+{1\over2n^2\pi^2}\Big).$$
There are two interesting things about this result:
- It is negative,
- It is independent of the width of the box.
The first is easier to come up with an explanation for. My attempt is that the potential of the sagging box is such that the box appears to sag into the negative territory of the potential axis, thus making the change in energy negative.
Edit: Another potential explanation (pardon the pun) that I just thought of is similar: Because the box sags down this creates more space for lower energy eigenstates within the box, thus creating a net energy shift of negative value.
The second I am struggling to come up with an explanation for. Perhaps a sagging potential is one of ratios rather than sizes. One thought is that, visually, if you were to stretch the box then the 'sag' would move upwards (i.e. be less sagged), such as one would expect when pulling either side of an elastic band or something. Thus counteracting any dependancy on the width of the box.
Any suggestions?
Edit
Here is the calculation of $\Delta E_n^{(1)}$:
We consider the perturbation potential $\Delta V=V(x)-V_0(x)$, where $V_0(x)$ is some potential that is visually similar to the potential given. A very similar potential (you can check this by plotting both of them as boxes) is the simple situation of a 'box' (i.e. $V(x)=0$ for $x\in [0,L]$). Thus $$\Delta V(x)=V(x)-V_0(x)=\frac{V_0x(x-L)}{L^2}-0=\frac{V_0x(x-L)}{L^2}.$$ The energy shift is then given by $$\Delta E=\left< \Delta V\right>=\int_0^Ldx\Phi^*\Delta v\Phi.$$ Using the stationary state solutions to the Schrodinger equation, $\Phi_n(x)=\sqrt{L/2}\sin{\frac{n\pi x}{L}}$, we get $$\Delta E=\int_0^Ldx\frac{L}{2}\sin^2\Big(\frac{n\pi x}{L}\Big)\frac{V_0x(x-L)}{L^2}=...=-V_0\Big(\frac{1}{6}+\frac{1}{2n^2\pi^2}\Big).$$