Energies for Particle in a Gaussian Potential Well

The Schrödinger equation for a particle in a one-dimensional Gaussian potential well , given by , has never been solved analytically. This Demonstration derives an approximation for the first few bound-state energies, , using the linear variational method. The wavefunction is approximated by a linear combination . It is convenient to take the basis functions as the corresponding orthonormalized eigenfunction of the linear harmonic oscillator: , where is the Hermite polynomial and is a scaling constant to be determined variationally. After evaluating the matrix elements over the selected set of basis functions, Mathematica can calculate the eigenvalues in a single step, from which we select only those with negative values. For convenience, we set , so that all distances are expressed in bohrs (Bohr radii) and energy quantities in hartrees.
The graphic shows the computed energy levels for selected values of , , , and , superposed on the potential energy function. By an estimation based on the WKB method, the number of bound states is approximated by .
The approximate eigenfunctions can also be calculated, applying the the built-in Mathematica function Eigenvectors. Since the resulting functions have very similar shapes and nodal structures as the corresponding harmonic oscillator eigenfunctions, we did not consider it worthwhile to include these in the Demonstration.


  • [Snapshot]
  • [Snapshot]
  • [Snapshot]


D.A. McQuarrie, Quantum Chemistry, Sausalito, CA: University Science Books, 1983, pp. 266-275. Or numerous other texts on quantum mechanics or quantum chemistry.
    • Share:

Embed Interactive Demonstration New!

Just copy and paste this snippet of JavaScript code into your website or blog to put the live Demonstration on your site. More details »

Files require Wolfram CDF Player or Mathematica.