Double-well potential

The so-called double-well potential is one of a number of quartic potentials of considerable interest in quantum mechanics, in quantum field theory and elsewhere for the exploration of various physical phenomena or mathematical properties since it permits in many cases explicit calculation without over-simplification.

Thus the "symmetric double-well potential" served for many years as a model to illustrate the concept of instanton as a pseudo-classical configuration in a Euclideanised field theory.^[1] In the simpler quantum mechanical context this potential served as a model for the evaluation of Feynman path integrals.^[2] or the solution of the Schrödinger equation by various methods for the purpose of obtaining explicitly the energy eigenvalues.

The "inverted symmetric double-well potential", on the other hand, served as a nontrivial potential in the Schrödinger equation for the calculation of decay rates^[3] and the exploration of the large-order behavior of asymptotic expansions.^[4]

The third form of the quartic potential is that of a "perturbed simple harmonic oscillator" or ″pure anharmonic oscillator″ having a purely discrete energy spectrum.

The fourth type of possible quartic potential is that of "asymmetric shape" of one of the first two named above.

The double-well and other quartic potentials can be treated by a variety of methods—the main methods being (a) a perturbation method (that of Dingle and Müller^[5]) which requires the imposition of boundary conditions, (b) the WKB method and (c) the path integral method. All cases are treated in detail in the book of Müller-Kirsten.^[6]

The symmetric double-well

The main interest in the literature has (for reasons related to field theory) focussed on the symmetric double-well (potential), and there on the quantum mechanical ground state. Since tunneling through the central hump of the potential is involved, the calculation of the eigenenergies of the Schrödinger equation for this potential is nontrivial. The case of the ground state is mediated by pseudoclassical configurations known as instanton and anti-instanton. In explicit form these are hyperbolic functions. As pseudoclassical configurations these naturally appear in semi-classical considerations—the summation of (widely separated) instanton-anti-instanton pairs being known as the dilute gas approximation. The ground state eigenenergy finally obtained is an expression containing the exponential of the Euclidean action of the instanton. This is an expression containing the factor $1/\hbar$ and is therefore described as a (classically) nonperturbative effect.

The stability of the instanton configuration in the path integral theory of a scalar field theory with symmetric double-well self-interaction is investigated by investigating the equation of small oscillations about he instanton. One finds that this equation is a Pöschl-Teller equation (i.e. a second order differential equation like the Schrödinger equation with Pöschl-Teller potential) with nonnegative eigenvalues. The nonnegativity of the eigenvalues is indicative of the stability of the instanton.^[7]

As stated above, the instanton is the pseudoparticle configuration defined on an infinite line of Euclidean time that communicates between the two wells of the potential and is responsible for the ground state of the system. The configurations correspondingly responsible for higher, i.e. excited, states are periodic instantons defined on a circle of Euclidean time which in explicit form are expressed in terms of Jacobian elliptic functions (the generalization of trigonometric functions). The evaluation of the path integral in these cases involves correspondingly elliptic integrals. The equation of small fluctuations about these periodic instantons is a Lamé equation whose solutions are Lamé functions. In cases of instability (as for the inverted double-well potential) this equation possesses negative eigenvalues indicative of this instability, i.e. decay (as shown in the paper of Liang, Müller-Kirsten and Tchrakian referred to above).

Application of the perturbation method of Dingle and Müller (applied originally to the Mathieu equation, i.e. a Schrödinger equation with cosine potential) requires exploitation of parameter symmetries of the Schrödinger equation for the quartic potential. One expands around one of the two minima of the potential. In addition this method requires matching of different branches of solutions in domains of overlap. The application of boundary conditions finally yields (as in the case of the periodic potential) the nonperturbative effect.

In terms of parameters as in the Schrödinger equation for the symmetric double-well potential in the following form

{\frac {d^{2}y(z)}{dz^{2}}}+[E-V(z)]y(z)=0,\;\;\;V(z)=-{\frac {1}{4}}z^{2}h^{4}+{\frac {1}{2}}c^{2}z^{4},\;\;c^{2}>0,h^{4}>0,

the eigenvalues for $q_{0}=1,3,5,...$ are found to be (see book of Müller-Kirsten, formula (18.175b), p. 525)

E_{\pm }(q_{0},h^{2})=-{\frac {h^{8}}{2^{5}c^{2}}}+{\frac {1}{\sqrt {2}}}q_{0}h^{2}-{\frac {c^{2}(3q_{0}^{2}+1)}{2h^{4}}}-{\frac {{\sqrt {2}}c^{4}q_{0}}{8h^{10}}}(17q_{0}^{2}+19)+O({\frac {1}{h^{16}}})

\;\;\;\mp {\frac {2^{q_{0}+1}h^{2}(h^{6}/2c^{2})^{q_{0}/2}}{{\sqrt {\pi }}2^{q_{0}/4}[(q_{0}-1)/2]!}}e^{-h^{6}/6{\sqrt {2}}c^{2}}.

Clearly these eigenvalues are asymptotically ( $h^{2}\rightarrow \infty$ ) degenerate as expected from the harmonic part of the potential. Observe that terms of the perturbative part of the result are alternately even or odd in $q_{0}$ and $h^{2}$ (as in corresponding results for Mathieu functions, Lamé functions, prolate spheroidal wave functions, oblate spheroidal wave functions and others).

In field theory contexts the above symmetric double-well potential is often written ( $\phi$ being a scalar field)

V(\phi )={\frac {m^{4}}{2g^{2}}}\left(1-{\frac {g^{2}\phi ^{2}}{m^{2}}}\right)^{2},

and the instanton is the solution $\phi _{c}(\tau )$ of the Newton-like equation

{\frac {d^{2}\phi }{d\tau ^{2}}}=V'(\phi ),\;\;\;{\frac {1}{2}}\left({\frac {d\phi }{d\tau }}\right)^{2}-V(\phi )=-E_{cl}=0

( $\tau$ being the Euclidean time), i.e

\phi _{c}(\tau )={\frac {m}{g}}\tanh \left[m(\tau -\tau _{0})\right].

The equation of small fluctuations $\eta$ about $\phi _{c},\phi =\phi _{c}+\eta ,$ is the Pöschl-Teller equation (see Pöschl-Teller potential)

\left[-{\frac {d^{2}}{d\tau ^{2}}}+V''(\phi _{c})\right]\eta _{n}(\tau )=\omega _{n}^{2}\eta _{n}(\tau ),

with

V''(\phi _{c})=4m^{2}-{\frac {6m^{2}}{\cosh ^{2}m(\tau -\tau _{0})}}.

Since all eigenvalues $\omega _{n}^{2}$ are positive or zero, the instanton configuration is stable and there is no decay.

In the more general case of $E_{cl}\neq 0$ the classical solution is the periodic instanton

\phi _{c}(\tau )={\frac {kb(k)}{g}}sn[b(k)(\tau -\tau _{0})],\;\;\;b(k)=m\left({\frac {2}{1+k^{2}}}\right)^{1/2},

where $k$ is the elliptic modulus of the periodic Jacobian elliptic function $sn$ . The small fluctuation equation is in this general case a Lamé equation. In the limit $k=1,E_{cl}=0$ the solution $\phi _{c}$ becomes the vacuum instanton solution,

\phi _{c}={\frac {m}{g}}\tanh[m(\tau -\tau _{0})].

The inverted double-well potential

Perturbation theory along with matching of solutions in domains of overlap and imposition of boundary conditions (different from those for the double-well) can again be used to obtain the eigenvalues of the Schrödinger equation for this potential. In this case, however, one expands around the central trough of the potential. The results are therefore different from those above.

In terms of parameters as in the Schrödinger equation for the inverted double-well potential in the following form

{\frac {d^{2}y(z)}{dz^{2}}}+[E-V(z)]y(z)=0,\;\;\;V(z)={\frac {1}{4}}h^{4}z^{2}-{\frac {1}{2}}c^{2}z^{4},\;\;h^{4}>0,\;c^{2}>0,

the eigenvalues for $q_{0}=1,3,5,...$ are found to be (see book of Müller-Kirsten, formula (18.86), p. 503)

E={\frac {1}{2}}q_{0}h^{2}-{\frac {3c^{2}}{4h^{4}}}(q_{0}^{2}+1)-{\frac {q_{0}c^{4}}{h^{10}}}(4q_{0}^{2}+29)+O({\frac {1}{h^{16}}})+i{\frac {2^{q_{0}}h^{2}(h^{6}/2c^{2})^{q_{0}/2}}{(2\pi )^{1/2}[(q_{0}-1)/2]!}}e^{-h^{6}/6c^{2}}.

The imaginary part of this expression agrees with the result of Bender and Wu^[8] (see their formula (3.36) and set $\hbar = 1$ , and in their notation $q_{0}=2K+1,h^{6}/2c^{2}=\epsilon )$ . This result plays an important role in the discussion and investigation of the large order behavior of perturbation theory.

The pure anharmonic oscillator

In terms of parameters as in the Schrödinger equation for the pure anharmonic oscillator in the following form

{\frac {d^{2}y(z)}{dz^{2}}}+[E-V(z)]y(z)=0,\;\;\;V(z)={\frac {1}{4}}h^{4}z^{2}+{\frac {1}{2}}c^{2}z^{4},\;\;h^{4}>0,\;c^{2}>0,

the eigenvalues for $q=q_{0}=1,3,5,...$ are found to be

E={\frac {1}{2}}qh^{2}+{\frac {3c^{2}}{4h^{4}}}(q^{2}+1)-{\frac {c^{4}}{h^{10}}}q(4q^{2}+29)+O({\frac {1}{h^{16}}}).

More terms can easily be calculated. Observe the coefficients of the expansion are alternately even or odd in $q$ and $h^{2}$ , as in all other cases. This is an important aspect of the solutions of the differential equation for quartic potentials.

General comments

The above results for the double-well and the inverted double-well can also be obtained by the path integral method (there via periodic instantons, cf. instantons), and the WKB method, though with the use of elliptic integrals and the Stirling approximation of the gamma function, all of which make the calculation more difficult.

References

↑ S. Coleman, The Whys of Subnuclear Physics, ed. A. Zichichi (Plenum Press, 1979), 805-916; S. Coleman, The Uses of Instantons, 1977 International School of Subnuclear Physics, Ettore Majorana.
↑ E. Gildener and A. Patrascioiu, Phys. Rev. D16 (1977) 423 (ground state calculation); J.-Q. Liang and H.J.W. Müller-Kirsten, Phys. Rev. D46 (1992) 4685 (all states)
↑ J.-Q. Liang and H.J.W. Müller-Kirsten, Phys. Rev. D50 (1994) 6519
↑ C.M. Bender and T.T. Wu, Phys. Rev. Lett. 21 (1968) 406, Phys. Rev. Lett. 27 (1971) 461, Phys. Rev. 184 (1969) 1231
↑ R.B.Dingle and H.J.W. Müller, J. reine angew. Math. 211 (1962) 11; in this reference the perturbation method is developed for the cosine potential, i.e. the Mathieu equation; see Mathieu function.
↑ Harald J.W. Müller-Kirsten, Introduction to Quantum Mechanics: Schrödinger Equation and Path Integral, 2nd ed., (World Scientific, 2012, ISBN 978-981-4397-73-5)
↑ J.-Q. Liang, H.J.W. Müller-Kirsten and D.H. Tchrakian, Phys. Lett. B282 (1992) 105
↑ C.M. Bender and T.T. Wu, Phys. Rev. D7 (1973) 1620.

This article is issued from Wikipedia. The text is licensed under Creative Commons - Attribution - Sharealike. Additional terms may apply for the media files.