Physics Library
 An open source physics library
Encyclopedia | Forums | Docs | Random | Template Test |  
Login
create new user
Username:
Password:
forget your password?
Main Menu
Sections

Talkback

Downloads

Information
Rayleigh-Ritz method (Definition)

The Rayleigh-Ritz method is an algorithm for obtaining approximate solutions to eigenvalue ODEs. It can be neatly summarized as follows:

  1. Choose an approximate form for the eigenfunction with the lowest eigenvalue (the ground state wavefunction, in the language of quantum mechanics). Include one or more free parameters.
  2. Find the expectation value of the eigenvalue with respect to the trial eigenfunction.
  3. Minimize the resulting equation with respect to the free parameter(s), hence finding a value for the free parameter.
  4. Substitute this new eigenfunction back into the expectation value.
  5. The expectation value obtained is an upper bound for the actual eigenvalue of the true eigenfunction.

Example

Consider the time independent Schrƶdinger equation for a one-dimensional harmonic oscillator potential:

$\displaystyle \left(-\frac{\hbar^2}{2m} \frac{\partial^2}{\partial x^2} + \frac{1}{2} m^2\omega^2\right)\psi = E\psi $
where $m$ is the mass of the particle in the well, and $\omega$ is the angular velocity a classical particle would move with in the well. This equation can be solved exactly using Frobenius' method, and leads to eigenfunctions of the form of Hermite polynomials multiplied by Gaussians, and half-integer eigenvalues of the form $E_n = (n+1/2)\hbar \omega$. Since the solutions are known, it is a good test case. We choose the ground state wavefunction of the infinite potential well as our trial eigenfunction:

$\displaystyle \psi = \frac{\cos(\frac{\pi x}{2a})}{\sqrt{a}} $
with $a$ as our free parameter. We now find the expectation value:

$\displaystyle \langle E \rangle = \langle \psi \vert \hat{H} \vert \psi \rangle = \int^a_{-a} \psi^* \hat{H} \psi \, dx $
Evaluating the integral, we find

$\displaystyle \langle E \rangle = \frac{\hbar^2 \pi^2}{8ma^2} +m\omega^2a^2\left(\frac{1}{6} - \frac{1}{pi^2}\right) $
We now minimise this with respect to $a$ to obtain:

$\displaystyle 2m\omega^2a\left(\frac{1}{6} - \frac{1}{\pi^2}\right) = \frac{\hbar^2 \pi^2}{4ma^2} $
Hence:

$\displaystyle a = \pi\left(\frac{3}{4(\pi^2-6)}\right)^{\frac{1}{4}}\left(\frac{\hbar}{m\omega}\right)^{\frac{1}{2}} $
Substituting this into the expecation value $\langle E \rangle$ we obtain

$\displaystyle \langle E \rangle = \frac{1}{2} \left(\frac{\pi^2-6}{3}\right)^{\frac{1}{2}}\hbar\omega $

$\displaystyle \langle E \rangle \approx 0.568\hbar \omega $
The analytical value is of course $0.5\hbar \omega$. Considering the crudeness of the approximation used, the result is impressive.



"Rayleigh-Ritz method" is owned by invisiblerhino.

View style:


Cross-references: Hermite polynomials, Frobenius method, velocity, mass, time independent Schrödinger equation, parameters, quantum mechanics, algorithm
There is 1 reference to this object.

This is version 2 of Rayleigh-Ritz method, born on 2008-03-08, modified 2008-03-10.
Object id is 265, canonical name is RayleighRitzMethod.
Accessed 966 times total.

Classification:
Physics Classification: 02.60.Lj (Ordinary and partial differential equations; boundary value problems)
Pending Errata and Addenda
None.
Discussion
Style: Expand: Order:

No messages.

Testing some escape charachters for html category with a generator has an injective cogenerator" now escape ā€¯ with "