|
|
| Line 1: |
Line 1: |
| The '''Ritz method''' is a direct method to find an approximate solution for [[boundary value problem]]s. The method is named after [[Walter Ritz]].
| | WoodBuyer Crosser from Montebello, spends time with interests which include dogs, property developers in singapore and antiques. Finds motivation through travel and just spent 5 months at Wet Tropics of Queensland.<br><br>Here is my web site ... [http://umangthefest.prayasastep.org/groups-2/new-launch-rental-2013/ http://umangthefest.prayasastep.org/] |
| | |
| In [[quantum mechanics]], a system of particles can be described in terms of an "energy functional" or [[Hamiltonian (quantum mechanics)|Hamiltonian]], which will measure the energy of any proposed configuration of said particles. It turns out that certain privileged configurations are more likely than other configurations, and this has to do with the [[eigenvalue|eigenanalysis]] ("analysis of characteristics") of this [[Hamiltonian system]]. Because it is often impossible to analyze all of the infinite configurations of particles to find the one with the least amount of energy, it becomes essential to be able to approximate this Hamiltonian in some way for the purpose of [[numerical analysis|numerical computations]].
| |
| | |
| The Ritz method can be used to achieve this goal. In the language of mathematics, it is exactly the [[finite element method]] used to compute the [[eigenvector]]s and [[eigenvalue]]s of a Hamiltonian system.
| |
| | |
| == Discussion ==
| |
| | |
| As with other [[variational method]]s, a [[ansatz|trial wave function]], <math>\Psi</math>, is tested on the system. This trial function is selected to meet boundary conditions (and any other physical constraints). The exact function is not known; the trial function contains one or more adjustable parameters, which are varied to find a lowest energy configuration.
| |
| | |
| It can be shown that the ground state energy, <math>E_0</math>, satisfies an inequality:
| |
| | |
| :<math> E_0 \le \frac{\langle \Psi | \hat{H}| \Psi \rangle}{\langle \Psi | \Psi \rangle}. </math>
| |
| | |
| That is, the ground-state energy is less than this value.
| |
| The trial wave-function will always give an expectation value larger than or equal to the ground-energy.
| |
| | |
| If the trial wave function is known to be [[orthogonality|orthogonal]] to the ground state, then it will provide a boundary for the energy of some excited state.
| |
| | |
| The Ritz ansatz function is a linear combination of ''N'' known basis functions <math>\left\lbrace\Psi_i\right\rbrace</math>, parametrized by unknown coefficients:
| |
| | |
| :<math> \Psi = \sum_{i=1}^N c_i \Psi_i. </math>
| |
| | |
| With a known Hamiltonian, we can write its expected value as
| |
| | |
| :<math> \varepsilon = \frac{\left\langle \displaystyle\sum_{i=1}^N c_i\Psi_i \right| \hat{H} \left| \displaystyle\sum_{i=1}^Nc_i\Psi_i \right\rangle}{\left\langle \left. \displaystyle\sum_{i=1}^N c_i\Psi_i \right| \displaystyle\sum_{i=1}^Nc_i\Psi_i \right\rangle} = \frac{\displaystyle\sum_{i=1}^N\displaystyle\sum_{j=1}^Nc_i^*c_jH_{ij}}{\displaystyle\sum_{i=1}^N\displaystyle\sum_{j=1}^Nc_i^*c_jS_{ij}} \equiv \frac{A}{B}. </math>
| |
| | |
| The basis functions are usually not orthogonal, so that the [[overlap matrix]] '''''S''''' has nonzero nondiagonal elements. Either <math>\left\lbrace c_i \right\rbrace</math> or <math>\left\lbrace c_i^* \right\rbrace</math> (the conjugation of the first) can be used to minimize the expectation value. For instance, by making the partial derivatives of <math>\varepsilon</math> over <math>\left\lbrace c_i^* \right\rbrace</math> zero, the following equality is obtained for every ''k'' = 1, 2, ..., ''N'':
| |
| | |
| :<math> \frac{\partial\varepsilon}{\partial c_k^*} = \frac{\displaystyle\sum_{j=1}^Nc_j(H_{kj}-\varepsilon S_{kj})}{B} = 0, </math>
| |
| | |
| which leads to a set of ''N'' [[secular equation]]s:
| |
| | |
| :<math>\sum_{j=1}^N c_j \left( H_{kj} - \varepsilon S_{kj} \right) = 0 \quad \text{for} \quad k = 1,2,\dots,N. </math>
| |
| | |
| In the above equations, energy <math>\varepsilon</math> and the coefficients <math>\left\lbrace c_j \right\rbrace</math> are unknown. With respect to '''''c''''', this is a homogeneous set of linear equations, which has a solution when the [[determinant]] of the coefficients to these unknowns is zero:
| |
| | |
| :<math>\det \left( H - \varepsilon S \right) = 0, </math>
| |
| | |
| which in turn is true only for ''N'' values of <math>\varepsilon</math>. Furthermore, since the Hamiltonian is a [[hermitian operator]], the '''''H''''' matrix is also [[hermitian matrix|hermitian]] and the values of <math>\varepsilon_i</math> will be real. The lowest value among <math>\varepsilon_i</math> (i=1,2,..,N), <math>\varepsilon_0</math>, will be the best approximation to the ground state for the basis functions used. The remaining ''N-1'' energies are estimates of excited state energies. An approximation for the wave function of state ''i'' can be obtained by finding the coefficients <math>\left\lbrace c_j \right\rbrace</math> from the corresponding secular equation.
| |
| | |
| == The relationship with the finite element method ==
| |
| | |
| In the language of the finite element method, the matrix <math>H_{kj}</math> is precisely the ''stiffness matrix'' of the Hamiltonian in the piecewise linear element space, and the matrix <math>S_{kj}</math> is the ''mass matrix''. In the language of linear algebra, the value <math>\epsilon</math> is an eigenvalue of the discretized Hamiltonian, and the vector <math>c</math> is a discretized eigenvector.
| |
| | |
| == Papers ==
| |
| *Walter Ritz (1909) "Über eine neue Methode zur Lösung gewisser Variationsprobleme der mathematischen Physik" ''Journal für die Reine und Angewandte Mathematik'', vol. '''135''', pages 1–61. Available on-line at: http://gdz.sub.uni-goettingen.de/no_cache/dms/load/img/?IDDOC=261182 .
| |
| *J.K. MacDonald, "Successive Approximations by the Rayleigh–Ritz Variation Method", ''Phys. Rev.'' '''43''' (1933) 830
| |
| | |
| == Books ==
| |
| * R.Courant and D.Hilbert, p.175
| |
| * G. Arfken, p.800
| |
| * E. Butkov, p.564
| |
| | |
| == External links ==
| |
| * [http://eom.springer.de/R/r082500.htm SpringerLink – Ritz method]
| |
| | |
| == See also ==
| |
| *[[Rayleigh–Ritz method]]
| |
| *[[Sturm–Liouville theory]]
| |
| *[[Hilbert space]]
| |
| | |
| [[Category:Perturbation theory]]
| |
| [[Category:Quantum chemistry]]
| |
WoodBuyer Crosser from Montebello, spends time with interests which include dogs, property developers in singapore and antiques. Finds motivation through travel and just spent 5 months at Wet Tropics of Queensland.
Here is my web site ... http://umangthefest.prayasastep.org/