Content deleted Content added
→See also: alpha Tags: Mobile edit Mobile web edit Advanced mobile edit |
|||
| (34 intermediate revisions by 26 users not shown) | |||
Line 1:
{{short description|Solution method for linear differential equations}}
{{Redirect2|WKB|WKBJ|other uses|WKB (disambiguation)|the television station in Live Oak, Florida|WKBJ-LD}}
In [[mathematical physics]], the '''WKB approximation''' or '''WKB method''' is a
The name is an initialism for '''Wentzel–Kramers–Brillouin'''. It is also known as the '''LG''' or '''Liouville–Green method'''. Other often-used letter combinations include '''JWKB''' and '''WKBJ''', where the "J" stands for Jeffreys.
== Brief history ==
This method is named after physicists [[Gregor Wentzel]], [[Hendrik Anthony Kramers]], and [[Léon Brillouin]], who all developed it in 1926.<ref name=Wentzel-1926/><ref name=Kramers-1926/><ref name=Brillouin-1926/><ref>{{harvnb|Hall|2013}} Section 15.1 </ref> In 1923,<ref name=Jefferys-1924/> mathematician [[Harold Jeffreys]] had developed a general method of approximating solutions to linear, second-order differential equations, a class that includes the [[Schrödinger equation]]. The Schrödinger equation itself was not developed until two years later, and Wentzel, Kramers, and Brillouin were apparently unaware of this earlier work, so Jeffreys is often neglected credit. Early texts in quantum mechanics contain any number of combinations of their initials, including WBK, BWK, WKBJ, JWKB and BWKJ. An authoritative discussion and critical survey has been given by Robert B. Dingle.<ref>{{cite book |first=Robert Balson |last=Dingle |title=Asymptotic Expansions: Their Derivation and Interpretation |publisher=Academic Press |year=1973 |isbn=0-12-216550-0 }}</ref>
Earlier appearances of essentially equivalent methods are: [[Francesco Carlini]] in 1817,<ref name=Carlini-1817/> [[Joseph Liouville]] in 1837,<ref name=Liouville/> [[George Green (mathematician)|George Green]] in 1837,<ref name=Green-1837/> [[Lord Rayleigh]] in 1912<ref name=Rayleigh-1912/> and [[Richard Gans]] in 1915.<ref name=Gans-1915/> Liouville and Green may be said to have founded the method in 1837, and it is also commonly referred to as the Liouville–Green or LG method.<ref>{{cite book
| title = Atmosphere-ocean dynamics
| author = Adrian E. Gill
Line 55:
}}</ref><ref>
{{Cite book
| author1=
| title=Introduction to Perturbation Methods, 2nd Ed
▲ | year=1996
| publisher=Springer
| isbn=
}}</ref><ref name=":0">{{cite book
| first1=Carl M.
Line 79 ⟶ 78:
<math display="block"> \epsilon^2 \frac{d^2 y}{dx^2} = Q(x) y, </math>
where <math>Q(x) \neq 0</math>. Substituting
<math display="block">y(x) = \exp \left[\frac{1}{\delta} \sum_{n=0}^\infty \delta^n
results in the equation
<math display="block">\epsilon^2\left[\frac{1}{\delta^2} \left(\sum_{n=0}^\infty \delta^
To [[leading-order|leading order]] in ''ϵ'' (assuming, for the moment, the series will be asymptotically consistent), the above can be approximated as
<math display="block">\frac{\epsilon^2}{\delta^2}
In the limit {{math|''δ'' → 0}}, the [[Method of dominant balance|dominant balance]] is given by
<math display="block">\frac{\epsilon^2}{\delta^2}
So {{mvar|δ}} is proportional to ''ϵ''. Setting them equal and comparing powers yields
<math display="block">\epsilon^0: \quad
which can be recognized as the [[eikonal equation]], with solution
<math display="block">
Considering first-order powers of {{mvar|ϵ}} fixes
<math display="block">\epsilon^1: \quad 2
This has the solution
<math display="block">
where {{math|''k''<sub>1</sub>}} is an arbitrary constant.
Line 104 ⟶ 103:
Higher-order terms can be obtained by looking at equations for higher powers of {{mvar|δ}}. Explicitly,
<math display="block">
for {{math|''n'' ≥ 2}}.
Line 112 ⟶ 111:
For the equation
<math display="block"> \epsilon^2 \frac{d^2 y}{dx^2} = Q(x) y, </math>
with {{math|''Q''(''x'') <0}} an [[analytic function]], the value <math>n_\max</math> and the magnitude of the last term can be estimated as follows:<ref>{{cite journal| last=Winitzki |first=S. |year=2005 |arxiv=gr-qc/0510001 |title=Cosmological particle production and the precision of the WKB approximation |journal=Phys. Rev. D |volume=72 |issue=10 |pages=104011, 14 pp |doi=10.1103/PhysRevD.72.104011 |bibcode = 2005PhRvD..72j4011W |s2cid=119152049 }}</ref>
<math display="block">n_\max \approx 2\epsilon^{-1} \left| \int_{x_0}^{x_{\ast}} \sqrt{-Q(z)}\,dz \right| , </math>
<math display="block">\delta^{n_\max}S_{n_\max}(x_0) \approx \sqrt{\frac{2\pi}{n_\max}} \exp[-n_\max], </math>
Line 123 ⟶ 122:
the number {{math|''n''<sub>max</sub>}} will be large, and the minimum error of the asymptotic series will be exponentially small.
== Application in non
[[File:WKB approximation example.svg|thumb|WKB approximation to the indicated potential. Vertical lines show the turning points]]
[[File:WKB approximation to probability density.svg|thumb|Probability density for the approximate wave function. Vertical lines show the turning points]]
Line 134 ⟶ 133:
The wavefunction can be rewritten as the exponential of another function {{math|S}} (closely related to the [[Action (physics)|action]]), which could be complex,
<math display="block">\Psi(\mathbf x) = e^{i S(\mathbf{x}) \over \hbar}, </math>
so that its substitution in
<math display="block">i\hbar \nabla^2 S(\mathbf x) - (\nabla S(\mathbf x))^2 = 2m \left( V(\mathbf x) - E \right),</math>
Line 140 ⟶ 139:
Next, the semiclassical approximation is used. This means that each function is expanded as a [[power series]] in {{mvar|ħ}}.
<math display="block">S = S_0 + \hbar S_1 + \hbar^2 S_2 + \cdots </math>
Substituting in the equation, and only retaining terms
<math display="block">(\nabla S_0+\hbar \nabla S_1)^2-i\hbar(\nabla^2 S_0) = 2m(E-V(\mathbf x)) </math>
which gives the following two relations:
Line 147 ⟶ 146:
2\nabla S_0 \cdot \nabla S_1 - i \nabla^2 S_0 = 0
\end{align}</math>
which can be solved for 1D systems, first equation resulting in:<math display="block">S_0(x) = \pm \int \sqrt{
Line 163 ⟶ 162:
Hence, when <math>E > V(x)</math>, the wavefunction can be chosen to be expressed as:<math display="block">\Psi(x') \approx C \frac{\cos{(\frac 1 \hbar \int |p(x)|\,dx} + \alpha) }{\sqrt{|p(x)| }} + D \frac{ \sin{(- \frac 1 \hbar \int |p(x)|\,dx} +\alpha)}{\sqrt{|p(x)| }} </math>and for <math>V(x) > E</math>,<math display="block">\Psi(x') \approx \frac{ C_{+} e^{
=== Validity of WKB solutions ===
Line 186 ⟶ 185:
\end{align} </math>
where <math display="inline">|S_0'(x)|= |p(x)| </math> is used and <math display="inline">\lambda(x) </math> is the local [[De Broglie waves|de Broglie wavelength]] of the wavefunction. The inequality implies that the variation of potential is assumed to be slowly varying.<ref name=":1" /><ref name=":2">{{Cite web |last=Zwiebach |first=Barton |title=Semiclassical approximation |url=https://ocw.mit.edu/courses/8-06-quantum-physics-iii-spring-2018/bf207c35150e1f5d93ef05d4664f406d_MIT8_06S18ch3.pdf}}</ref> This condition can also be restated as the fractional change of <math display="inline">E-V(x) </math> or that of the momentum <math display="inline">p(x) </math>, over the wavelength <math display="inline">\lambda </math>, being much smaller than <math display="inline">1 </math>.<ref>{{Cite book |last1=Bransden |first1=B. H. |url=https://books.google.com/books?id=ST_DwIGZeTQC |title=Physics of Atoms and Molecules |last2=Joachain |first2=Charles Jean |date=2003 |publisher=Prentice Hall |isbn=978-0-582-35692-4 |pages=140–141 |language=en}}</ref>
Line 219:
==== First classical turning point ====
For <math>U_1 < 0</math> ie. decreasing potential condition or <math>x=x_1
</math> in the given example shown by the figure, we require the exponential function to decay for negative values of x so that wavefunction for it to go to zero. Considering Bairy functions to be the required connection formula, we get:<ref name=":3">{{Cite journal |
<math display="block">\begin{align}
Line 228:
We cannot use Airy function since it gives growing exponential behaviour for negative x. When compared to WKB solutions and matching their behaviours at <math>\pm \infty </math>, we conclude:
<math>
<math>
Thus, letting some normalization constant be <math>N </math>, the wavefunction is given for increasing potential (with x) as:<ref name=":1" />
Line 250:
We cannot use Bairy function since it gives growing exponential behaviour for positive x. When compared to WKB solutions and matching their behaviours at <math>\pm \infty </math>, we conclude:
<math>
<math>D=
Thus, letting some normalization constant be <math>N' </math>, the wavefunction is given for increasing potential (with x) as:<ref name=":1" />
Line 269:
Either way, the condition on the energy is a version of the [[Bohr–Sommerfeld quantization]] condition, with a "[[Lagrangian Grassmannian#Maslov index|Maslov correction]]" equal to 1/2.<ref>{{harvnb|Hall|2013}} Section 15.2</ref>
It is possible to show that after piecing together the approximations in the various regions, one obtains a good approximation to the actual [[eigenfunction]]. In particular, the Maslov-corrected Bohr–Sommerfeld energies are good approximations to the actual eigenvalues of the Schrödinger operator.<ref>{{harvnb|Hall|2013}} Theorem 15.8</ref> Specifically, the error in the energies is small compared to the typical spacing of the quantum energy levels. Thus, although the "old quantum theory" of Bohr and Sommerfeld was ultimately replaced by the Schrödinger equation, some vestige of that theory remains, as an approximation to the eigenvalues of the appropriate Schrödinger operator.
==== General connection conditions ====
Line 315:
Thus we conclude that, for <math display="inline">n = 1,2,3,\cdots </math><math display="block">\int_{x_1}^{x_2} \sqrt{2m \left( E-V(x)\right)}\,dx = \left(n-\frac 1 4\right)\pi \hbar </math>In 3 dimensions with spherically symmetry, the same condition holds where the position x is replaced by radial distance r, due to its similarity with this problem.<ref>{{Cite book |last=Weinberg |first=Steven |url=http://dx.doi.org/10.1017/cbo9781316276105 |title=Lectures on Quantum Mechanics |date=2015-09-10 |publisher=Cambridge University Press |isbn=978-1-107-11166-0 |edition=2nd |pages=204|doi=10.1017/cbo9781316276105 }}</ref>
=== Bound states within 2 rigid wall ===
Line 349:
\end{cases}</math>
The wavefunction solutions of the above can be solved
<math display="block">\int_{x_1}^{x_2} \sqrt{2m \left( E-V(x)\right)}\,dx = (n_{\text{odd}}+1/2)\pi \hbar</math>
Letting <math display="inline">n_{\text{odd}}=2n-1 </math> where <math display="inline">n = 1,2,3,\cdots </math>, solving for <math display="inline">E </math> with given <math>V(x) = mg|x|</math>, we get the quantum mechanical energy of a bouncing ball:<ref>{{Cite book |
<math display="block">E = {\left(3\left(n-\frac 1 4\right)\pi\right)^{\frac 2 3} \over 2}(mg^2\hbar^2)^{\frac 1 3}. </math>
This result is also consistent with the use of equation from [[bound state]] of one rigid wall without needing to consider an alternative potential.
=== Quantum Tunneling ===
Line 372:
<math display="block">
A \exp({ i p_0 x \over \hbar} ) + B \exp({- i p_0 x \over \hbar}) & \text{if } x < x_1 \\
\frac{C}{\sqrt{|p(x)|}}\exp{(-\frac 1 \hbar \int_{x_1}^{x} |p(x)| dx )} & \text{if } x_2 \geq x \geq x_1\\
Line 381:
\end{cases} </math>
By the requirement of continuity of wavefunction and its derivatives, the following relation can be shown:<math display="block">\frac {|
where <math>a_1 = |p(x_1)|</math> and <math>a_2 = |p(x_2)| </math>.
Line 395:
<math display="inline">J_{\text{ref.}} = \frac{\hbar}{2m}(\frac{2p_0}{\hbar}|B|^2) </math>
<math display="inline">J_{\text{trans.}} = \frac{\hbar}{2m}(\frac{2p_0}{\hbar}|
Thus, the [[transmission coefficient]] is found to be:
<math display="block">T = \frac {|
where <math display="inline">p(x) = \sqrt {2m( E - V(x))} </math>, <math>a_1 = |p(x_1)|</math> and <math>a_2 = |p(x_2)| </math>. The result can be stated as <math display="inline">T \sim ~ e^{-2\gamma} </math> where <math display="inline">\gamma = \int_{x_1}^{x_2} |p(x')| dx' </math>.<ref name=":1" />
Line 413:
* [[Langer correction]]
* [[Maslov index]]
* [[Method of dominant balance]]
* [[Method of matched asymptotic expansions]]
* [[Method of steepest descent]]
Line 423 ⟶ 424:
==References==
{{Reflist
}}
===
*{{cite book | author=Child, M. S. | title=Semiclassical mechanics with molecular applications | year=1991 | publisher = Clarendon Press | ___location=Oxford | isbn=0-19-855654-3}}
* {{cite book| last1=Fröman|first1=N. |last2= Fröman |first2= P.-O.| title=JWKB Approximation: Contributions to the Theory | publisher=North-Holland |___location = Amsterdam | year=1965}}
*{{cite book | author=Griffiths, David J. | title=Introduction to Quantum Mechanics | edition = 2nd | publisher=Prentice Hall |year=2004 |isbn=0-13-111892-7}}
*{{citation|first=Brian C.|last=Hall | title=Quantum Theory for Mathematicians | series=Graduate Texts in Mathematics| volume=267 | publisher=Springer|year=2013|bibcode=2013qtm..book.....H | isbn=978-1461471158}}
*{{cite book | author=Liboff, Richard L. | title=Introductory Quantum Mechanics | edition = 4th | publisher=Addison-Wesley | year=2003 |isbn=0-8053-8714-5| author-link=Liboff, Richard L }}
*{{cite book | author=Olver, Frank William John |author-link=Frank William John Olver | title=Asymptotics and Special Functions | url=https://archive.org/details/asymptoticsspeci0000olve | url-access=registration | publisher=Academic Press | year=1974 | isbn=0-12-525850-X}}
*{{cite book | author=Razavy, Mohsen | title=Quantum Theory of Tunneling | url=https://archive.org/details/quantumtheoryoft0000raza | url-access=registration | publisher=World Scientific | year=2003 | isbn=981-238-019-1}}
▲*{{cite book | author=Carlini, Francesco | year=1817 | title=Ricerche sulla convergenza della serie che serva alla soluzione del problema di Keplero | publisher=Milano | author-link=Francesco Carlini }}
▲*{{cite journal | author=Liouville, Joseph | year=1837 | title=Sur le développement des fonctions et séries..| journal=Journal de Mathématiques Pures et Appliquées | volume=1 | pages=16–35 | author-link=Joseph Liouville }}
▲*{{cite journal | author=Green, George | year=1837 | title=On the motion of waves in a variable canal of small depth and width | journal=Transactions of the Cambridge Philosophical Society | volume=6 | pages=457–462 | author-link=George Green (mathematician) }}
▲*{{cite journal | author=Rayleigh, Lord (John William Strutt) | year=1912 | title=On the propagation of waves through a stratified medium, with special reference to the question of reflection | journal=[[Proceedings of the Royal Society A]] | volume=86 | pages=207–226 | doi=10.1098/rspa.1912.0014 |bibcode = 1912RSPSA..86..207R | issue=586 | author-link=Lord Rayleigh | doi-access=free }}
▲*{{cite journal | author=Gans, Richard | year=1915 | title=Fortplantzung des Lichts durch ein inhomogenes Medium | journal=Annalen der Physik | volume=47 | issue=14 | pages=709–736 | doi = 10.1002/andp.19153521402 |bibcode = 1915AnP...352..709G | url=https://zenodo.org/record/1447303 | author-link=Richard Gans }}
▲*{{cite journal | author=Jeffreys, Harold | year=1924 | title=On certain approximate solutions of linear differential equations of the second order | journal=Proceedings of the London Mathematical Society | volume=23 | pages=428–436 | doi=10.1112/plms/s2-23.1.428 | author-link=Harold Jeffreys }}
▲*{{cite journal | author=Brillouin, Léon | year=1926 | title=La mécanique ondulatoire de Schrödinger: une méthode générale de resolution par approximations successives | journal=Comptes Rendus de l'Académie des Sciences | volume=183 | pages=24–26 | author-link=Léon Brillouin }}
▲*{{cite journal | author=Kramers, Hendrik A. | year=1926 | title=Wellenmechanik und halbzahlige Quantisierung | journal=Zeitschrift für Physik | volume=39 |pages=828–840 | doi=10.1007/BF01451751 |bibcode = 1926ZPhy...39..828K | issue=10–11 | s2cid=122955156 | author-link=Hendrik Anthony Kramers }}
▲*{{cite journal | author=Wentzel, Gregor | year=1926 | title=Eine Verallgemeinerung der Quantenbedingungen für die Zwecke der Wellenmechanik | journal=Zeitschrift für Physik | volume=38 | pages=518–529 | doi=10.1007/BF01397171 |bibcode = 1926ZPhy...38..518W | issue=6–7 | s2cid=120096571 | author-link=Gregor Wentzel }}
==External links==
Line 451 ⟶ 449:
[[Category:Approximations]]
[[Category:Asymptotic analysis]]
[[Category:Mathematical physics]]
| |||