Gaussian variational bounds for a quartic oscillator ($E_0\le (3{\hslash }^2/8m)(6mk/{\hslash }^2{)}^{1/3}$)

For $H=p^2/(2m)+k{x}^4$ with $k>0$, the normalized even trial wavefunction $f_{\alpha }=(\alpha /\pi {)}^{1/4}{e}^{-\alpha x^2/2}$ has Rayleigh quotientThe Gaussian moment recurrence gives the potential term, and integration by parts gives the kinetic term. Its minimum is at ${\alpha }^3=6mk/{\hslash }^2$, proving the displayed upper bound by the Rayleigh-Ritz variational principle.

An odd trial family $g_{\beta }=(4{\beta }^3/\pi {)}^{1/4}x{e}^{-\beta x^2/2}$ is orthogonal to the exact even ground state. The Sturm oscillation theorem identifies the lowest odd state as the first excited state, so minimizing within this family bounds $E_1$:Orthogonality to an approximate ground trial alone would not justify this excited-state bound; odd parity ensures orthogonality to the exact ground state. Both bounds exhibit the characteristic quartic-oscillator energy scale $({\hslash }^{4}k/m^2{)}^{1/3}$.

For $H=-d^2/d{x}^2+x^2+ϵx^4$, use the normalized real wavefunction ${\psi }_{\alpha }=(\alpha /\pi {)}^{1/4}{e}^{-\alpha x^2/2}$, $\alpha >0$. The Gaussian integral gives $⟨x^2⟩=1/(2\alpha )$ and $⟨x^4⟩=3/(4{\alpha }^2)$, while the kinetic contribution is $\alpha /2$. Thus the Rayleigh quotient is the displayed expression. For positive $ϵ$, its stationary width satisfies ${\alpha }^3-\alpha -3ϵ=0$. At $ϵ=2$ the unique positive solution is $\alpha =2$, giving $E_0\le 13/8$ by the Rayleigh-Ritz variational principle. For $ϵ<0$, taking $\alpha \downarrow 0$ makes the negative quartic contribution dominate, so the energy is unbounded below even within this normalized family. The positive-coupling bound and negative-coupling instability concern different signs of the same variational expression.