删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

Representation of the quantum mechanical wavefunction by orthogonal polynomials in the energy and ph

本站小编 Free考研考试/2022-01-02

A D Alhaidari,Saudi Center for Theoretical Physics, PO Box 32741, Jeddah 21438, Saudi Arabia

Received:2019-06-27Revised:2019-09-23Accepted:2019-10-31Online:2019-12-30


Abstract
We present a formulation of quantum mechanics based on the theory of orthogonal polynomials. The wavefunction is expanded over a complete set of square integrable basis where the expansion coefficients are orthogonal polynomials in the energy and physical parameters. Information about the corresponding physical systems (both structural and dynamical) are derived from the properties of these polynomials. We demonstrate that an advantage of this formulation is that the class of exactly solvable quantum mechanical problems becomes larger than in the conventional formulation (see, for example, table 3 in the text). We limit our investigation in this work to the Askey classification scheme of hypergeometric orthogonal polynomials and focus on the Wilson polynomial and two of its limiting cases (the Meixner–Pollaczek and continuous dual Hahn polynomials). Nonetheless, the formulation is amenable to other classes of orthogonal polynomials.
Keywords: tridiagonal representation;orthogonal polynomials;recursion relation;asymptotics;energy spectrum;phase shift


PDF (453KB)MetadataMetricsRelated articlesExportEndNote|Ris|BibtexFavorite
Cite this article
A D Alhaidari. Representation of the quantum mechanical wavefunction by orthogonal polynomials in the energy and physical parameters. Communications in Theoretical Physics, 2020, 72(1): 015104- doi:10.1088/1572-9494/ab5d00

1. Introduction

In physics, we are accustomed to writing vector quantities (e.g. force, velocity, electric field, etc) in terms of their components in some conveniently chosen vector space. For example, the force $\vec{F}$ is written in three dimensional space with Cartesian coordinates as $\vec{F}={f}_{x}\hat{x}+{f}_{y}\hat{y}+{f}_{z}\hat{z},$ where $\left\{{f}_{x},{f}_{y},{f}_{z}\right\}$ are the components of the force along the unit vectors $\left\{\hat{x},\hat{y},\hat{z}\right\}.$ These components contain all physical information about the quantity whereas the unit vectors (basis) are dummy, but must form a complete set to allow for a faithful physical representation. In fact, we can as well write the same force in another coordinates, say the spherical coordinates with basis $\left\{\hat{r},\hat{\theta },\hat{\varphi }\right\},$ as $\vec{F}={f}_{r}\hat{r}+{f}_{\theta }\hat{\theta }+{f}_{\varphi }\hat{\varphi },$ where $\left\{{f}_{r},{f}_{\theta },{f}_{\varphi }\right\}$ are the new components that contain the same physical information. And so on, where in general we can write $\vec{F}=\displaystyle {\sum }_{n}{f}_{n}{\hat{x}}_{n}.$ The basis (unit vectors) $\left\{{\hat{x}}_{n}\right\}$ is chosen conveniently depending on the symmetry of the problem (e.g. rectangular, spherical, cylindrical, elliptical, etc). It is basic knowledge that if two unit vectors are independent then they can be, but not necessarily, orthogonal (i.e. $\left\langle {\hat{x}}_{n}| {\hat{x}}_{m}\right\rangle ={\delta }_{nm}$), however they cannot be collinear (i.e. $\left\langle {\hat{x}}_{n}| {\hat{x}}_{m}\right\rangle =\pm 1$ for $n\ne m$). Therefore, the chosen basis set must be complete and consisting of independent elements but they do not have to be orthogonal to each other.

In quantum mechanics we can also think of the wavefunction $\psi (x)$ as a local vector and write it in terms of its components $\left\{{f}_{n}\right\}$ along some local unit vectors (basis) $\left\{{\phi }_{n}(x)\right\}$ as $\psi (x)=\displaystyle {\sum }_{n}{f}_{n}\,{\phi }_{n}(x).$ All physical information about the system are contained in the components (expansion coefficients) $\left\{{f}_{n}\right\}.$ On the other hand, the basis set $\left\{{\phi }_{n}(x)\right\}$ is dummy but, in analogy to the unit vectors in the above example, must be normalizable (square integrable) and complete. If the physical system is also associated with a set of real parameters $\left\{\mu \right\},$ then the components of the wavefunction at an energy E could be written as parameterized functions of the energy, $\left\{{f}_{n}^{\mu }(E)\right\},$ and the state of the system becomes$ \begin{eqnarray}{\psi }_{E}^{\mu }(x)=\displaystyle \sum _{n}{f}_{n}^{\mu }(E)\,{\phi }_{n}(x).\end{eqnarray}$

We have shown elsewhere [1] that if we write ${f}_{n}^{\mu }(E)={f}_{0}^{\mu }(E){P}_{n}^{\mu }(\varepsilon ),$ where ϵ is some proper function of E and $\left\{\mu \right\},$ then completeness of the basis and energy normalization of the density of state make $\left\{{P}_{n}^{\mu }(\varepsilon )\right\}$ a complete set of orthogonal polynomials with an associated non-negative weight function ${\rho }^{\mu }(\varepsilon )={[{f}_{0}^{\mu }(E)]}^{2}.$ That is,$ \begin{eqnarray}\displaystyle \int {\rho }^{\mu }(\varepsilon ){P}_{n}^{\mu }(\varepsilon ){P}_{m}^{\mu }(\varepsilon )\,{\rm{d}}\varepsilon ={\delta }_{n,m}.\end{eqnarray}$

Therefore, we can rewrite the wavefunction expansion (1) as follows$ \begin{eqnarray}{\psi }_{E}^{\mu }(x)=\sqrt{{\rho }^{\mu }(\varepsilon )}\displaystyle \sum _{n}{P}_{n}^{\mu }(\varepsilon )\,{\phi }_{n}(x).\end{eqnarray}$

In cases where ${\phi }_{n}(x)\sim \sigma (x){x}^{n},$ this formula suggests that the wavefunction could also be considered as the generating function of the energy polynomials $\left\{{P}_{n}^{\mu }(\varepsilon )\right\}.$ In their treatment of quasi-exactly solvable systems, Bender and Dunne were among the first to realize this fact [2].

Building on our past and recent experience in dealing with various physical systems [1, 39], we limit the current investigation of the energy polynomials to the Askey classification scheme of hypergeometric orthogonal polynomials. We focus on the Wilson polynomial and two of its limiting cases [10]: the Wilson→continuous dual Hahn (${}_{4}F_{3}\to {}_{3}F_{2}$) and the Wilson→Meixner–Pollaczek (${}_{4}F_{3}\to {}_{3}F_{2}\to {}_{2}F_{1}$). There is another chain in the Askey diagram that deals with the discrete version of hypergeometric orthogonal polynomials. At the top of this chain sits the Racah polynomial and it includes in the limit other discrete polynomials like the Meixner, Krawtchouk, Hahn, dual Hahn, etc. Now, the connection between scattering and the asymptotics ($n\to \infty $) of continuous orthogonal polynomials is well-stablished [1113]. Using this connection, scattering information about the system, whose wavefunction is written as shown by equation (3), is readily available from the asymptotics ${P}_{n\to \infty }^{\mu }(\varepsilon ).$ Table 1 shows a list of the asymptotics of most of the continuous polynomials in the Askey scheme of hypergeometric orthogonal polynomials. The asymptotics of the Wilson polynomial is obtained in section 2 below whereas the asymptotic of the continuous Hahn polynomial is derived in appendix B of [9] (see equation (37) therein). The rest are either well known or derived elsewhere (see, for example, the appendix in [1] where the asymptotics of the Meixner–Pollaczek and continuous dual Hahn polynomials were obtained). We should note that throughout this paper we deal with the orthonormal version of orthogonal polynomials where the right-hand side of the orthogonality relation (2) is ${\delta }_{n,m},$ the corresponding three-term recursion relation is symmetric, and ${P}_{0}^{\mu }(\varepsilon )=1.$ Due to the prime significance of the Wilson polynomial (being at the top of the Askey scheme), we present the derivation of its asymptotics in section 2. We should note that despite the fact that the Meixner–Pollaczek and continuous dual Hahn polynomials are limiting cases of the Wilson polynomial, we could not use the asymptotics of the Wilson polynomial because we cannot interchange this limit with the asymptotic limit. Table 1 suggests that all orthogonal polynomials relevant to our study are those with the following asymptotic behavior$ \begin{eqnarray}\begin{array}{l}{P}_{n}^{\mu }(\varepsilon )\approx {n}^{-\tau }{A}^{\mu }(\varepsilon )\left\{\cos \left[{n}^{\xi }\,\theta (\varepsilon )+\varphi (\varepsilon )\mathrm{log}\,n+{\delta }^{\mu }(\varepsilon )\right]\right.\\ \,\,\,\,+\,\left.O({n}^{-1})\right\},\end{array}\end{eqnarray}$where τ and ξ are real dimensionless constants with the value of τ being dependent on the type of normalization of the polynomial, which we consistently choose as ortho-normalization where the right-hand side of the orthogonality relation (2) is ${\delta }_{n,m}$ and ${P}_{0}^{\mu }(\varepsilon )=1.$ In the asymptotic formula (4), ${A}^{\mu }(\varepsilon )$ is the scattering amplitude and ${\delta }^{\mu }(\varepsilon )$ is the phase shift. For the polynomials listed in table 1, the values of the parameters τ and ξ belong to the sets $\tau \in \left\{0,\tfrac{1}{4},\tfrac{1}{2}\right\},$$\xi \in \left\{0,\tfrac{1}{2},1\right\}$ and we observe the following three alternative scenarios:(1)$\theta (\varepsilon )\ne 0,$$\varphi (\varepsilon )=0.$
(2)$\theta (\varepsilon )=0,$$\varphi (\varepsilon )\ne 0.$
(3)$\theta (\varepsilon )\ne 0,$$\varphi (\varepsilon )\ne 0.$



Table 1.
Table 1.Asymptotics ($n\to \infty $) of most of the continuous polynomials in the Askey scheme of hypergeometric orthogonal polynomials. The polynomials are shown in the second column in their orthonormal version. The asymptotics of the Wilson polynomial is obtained here in section 2 whereas the asymptotics of the continuous Hahn polynomial is derived in [9]. The asymptotics of the Meixner–Pollaczek and continuous dual Hahn polynomials are obtained in [1]. The rest are well known.
PolynomialOrthonormal versionxAsymptotics [1, 9]τξ$\theta (x)$$\varphi (x)$
Laguerre$\sqrt{\displaystyle \frac{n!}{{(\nu +1)}_{n}}}\,{L}_{n}^{\nu }(x)$$x\geqslant 0$${n}^{-1/4}{A}_{L}(x)\cos \left[2\sqrt{nx}-\tfrac{\pi }{2}\left(\nu +\tfrac{1}{2}\right)\right]$$1/4$$1/2$$2\sqrt{x}$0
Jacobi$\begin{array}{c}\sqrt{\tfrac{2n+\mu +\nu +1}{\mu +\nu +1}\tfrac{n\,!{(\mu +\nu +1)}_{n}}{{(\mu +1)}_{n}{(\nu +1)}_{n}}}\\ \times {P}_{n}^{(\mu ,\nu )}(\cos \,x)\end{array}$$\pi \geqslant x\geqslant 0$${A}_{J}(x)\cos \left[\left(n+\tfrac{\mu +\nu +1}{2}\right)x-\tfrac{\pi }{2}\left(\mu +\tfrac{1}{2}\right)\right]$01x0
Meixner–Pollaczek${P}_{n}^{\mu }(x;\theta )$$x\in {\mathbb{R}}$${n}^{-1/2}{A}_{MP}(x)\cos \left[n\theta -x\,\mathrm{log}\,n+{\delta }_{MP}(x)\right]$$1/2$1θx
Continuous Hahn${{\mathscr{P}}}_{n}^{\mu }(x;\nu ;a,b)$$x\in {\mathbb{R}}$${n}^{-1/2}{A}_{H}(x)\cos \left[n\tfrac{\pi }{2}-(2x+a-b)\mathrm{log}\,n+{\delta }_{H}(x)\right]$$1/2$1$\pi /2$$-(2x+a-b)$
Continuous dual Hahn${S}_{n}^{\mu }({x}^{2};a,b)$$x\geqslant 0$${n}^{-1/2}{A}_{dH}(x)\cos \left[x\,\mathrm{log}\,n+{\delta }_{dH}(x)\right]$$1/2$**0x
Wilson${W}_{n}^{\mu }({x}^{2};\nu ;a,b)$$x\geqslant 0$${n}^{-1/2}{A}_{W}(x)\cos \left[2x\,\mathrm{log}\,n+{\delta }_{W}(x)\right]$$1/2$**02x

New window|CSV

In the third scenario, the simultaneous presence of the logarithmic term ($\mathrm{log}\,n$) and power term (${n}^{\xi }$), like in the case of the Meixner–Pollaczek and continuous Hahn polynomials, is very interesting and a source of curiosity. The understanding of this behavior and its physical implication should be very fruitful. Unfortunately, we do not have the needed expertise to address this issue at present.

Bound states, if they exist, occur at a set (infinite or finite) of energies that make the scattering amplitude vanish. That is, the kth bound state occurs at an energy ${E}_{k}=E({\varepsilon }_{k})$ such that ${A}^{\mu }({\varepsilon }_{k})=0$ and the corresponding bound state wavefunction is written as11There is an alternative, but equivalent, method for obtaining the energy spectrum from the phase shift by calculating the poles of the corresponding scattering matrix. That is, the phase shift angle at those energies becomes half-odd integer of π making the tangent of the phase angle diverge. Frequently, such condition occurs when the argument of the gamma function that appears in the phase shift becomes a negative integer (or zero) corresponding to the energy level in the spectrum. This alternative method is addressed in the book by Landau and Lifshitz [14] and recently by Chen et al in [15, 16], etc.
$ \begin{eqnarray}{\psi }_{k}^{\mu }(x)=\sqrt{{\omega }^{\mu }({\varepsilon }_{k})}\displaystyle \sum _{n}{Q}_{n}^{\mu }({\varepsilon }_{k})\,{\phi }_{n}(x),\end{eqnarray}$where $\left\{{Q}_{n}^{\mu }({\varepsilon }_{k})\right\}$ are the discrete version of the polynomials $\left\{{P}_{n}^{\mu }(\varepsilon )\right\}$ and ${\omega }^{\mu }({\varepsilon }_{k})$ is the associated discrete weight function. That is, $\displaystyle {\sum }_{k}{\omega }^{\mu }({\varepsilon }_{k}){Q}_{n}^{\mu }({\varepsilon }_{k}){Q}_{m}^{\mu }({\varepsilon }_{k})={\delta }_{n,m}.$ Sometimes, the quantum system consists of continuous as well as discrete energy spectra simultaneously. In that case, the total wavefunction is written as follows$ \begin{eqnarray}{{\rm{\Psi }}}^{\mu }(x,t)=\displaystyle \int {{\rm{e}}}^{-{\rm{i}}Et}{\psi }_{E}^{\mu }(x){\rm{d}}E+\displaystyle {\sum }_{k}{{\rm{e}}}^{-{\rm{i}}{E}_{k}t}{\psi }_{k}^{\mu }(x),\end{eqnarray}$and the appropriate polynomial orthogonality becomes$ \begin{eqnarray}\begin{array}{l}\displaystyle \int {\rho }^{\mu }(\varepsilon ){P}_{n}^{\mu }(\varepsilon ){P}_{m}^{\mu }(\varepsilon )\,{\rm{d}}\varepsilon +\displaystyle {\sum }_{k}{\omega }^{\mu }({\varepsilon }_{k}){P}_{n}^{\mu }({\varepsilon }_{k}){P}_{m}^{\mu }({\varepsilon }_{k})\\ \,=\,{\delta }_{n,m}.\end{array}\end{eqnarray}$

Now, the type of orthogonal polynomials associated with a given physical system depends on the structure of its energy spectrum: purely continuous, purely discrete or a mix of both and on whether the discrete energy spectrum is infinite or finite. Table 2 summarizes this association.


Table 2.
Table 2.The orthogonal polynomial(s) associated with a given physical system as a function of the structure of its energy spectrum. Note the matching of the spectra of the physical system and that of its corresponding orthogonal polynomial except when the physical system has a mix of a continuous and an infinite discrete energy spectrum then there is a need for two polynomials to represent the system: a continuous one and a discrete one with infinite spectrum.
Energy spectrumPolynomial spectrum
DiscreteDiscrete
ContinuousFiniteInfiniteContinuousFiniteInfinite

New window|CSV


Table 3.
Table 3.Partial list of exactly solvable potential functions in the polynomial class associated with ${H}_{n}^{(\mu ,\nu )}(z;\alpha ,\theta ).$ The coordinate transformation y(x) enters in the basis (44) that supports a tridiagonal matrix representation for the corresponding wave operator. The presence of the ${V}_{1}$ term in all of these potentials inhibits exact solvability of the Schrödinger wave equation in the standard formulation of quantum mechanics.
V(x)xy(x)
${V}_{0}+\displaystyle \frac{{V}_{+}-{V}_{-}\,\sin (\pi x/L)}{{\cos }^{2}(\pi x/L)}+{V}_{1}\,\sin \left(\pi x/L\right)$$-\tfrac{1}{2}L\leqslant x\leqslant +\tfrac{1}{2}L$$\sin \left(\pi x/L\right)$
$\displaystyle \frac{1}{1-{(x/L)}^{2}}\left\{{V}_{0}+\displaystyle \frac{{V}_{+}}{{(x/L)}^{2}}+\displaystyle \frac{{V}_{-}}{1-{(x/L)}^{2}}+{V}_{1}\left[2{(x/L)}^{2}-1\right]\right\}$$0\leqslant x\leqslant L$$2{\left(x/L\right)}^{2}-1$
$\displaystyle \frac{1}{{{\rm{e}}}^{\lambda x}-1}\left[{V}_{0}+{V}_{-}{{\rm{e}}}^{\lambda x}+\displaystyle \frac{{V}_{+}}{1-{{\rm{e}}}^{-\lambda x}}+{V}_{1}\left(1-2{{\rm{e}}}^{-\lambda x}\right)\right]$$x\geqslant 0$$1-2{{\rm{e}}}^{-\lambda x}$
${V}_{-}+\displaystyle \frac{{V}_{+}}{{\sinh }^{2}(\lambda x)}+\displaystyle \frac{{V}_{0}+{V}_{1}[2{\tanh }^{2}(\lambda x)-1]}{{\cosh }^{2}(\lambda x)}$$x\geqslant 0$$2{\tanh }^{2}(\lambda x)-1$
${V}_{+}-{V}_{-}\,\tanh (\lambda x)+\displaystyle \frac{{V}_{0}+{V}_{1}\,\tanh (\lambda x)}{{\cosh }^{2}(\lambda x)}$$-\infty \lt x\lt +\infty $$\tanh (\lambda x)$
${V}_{0}+\displaystyle \frac{{V}_{+}}{{\sin }^{2}(\pi x/L)}+\displaystyle \frac{{V}_{-}}{{\cos }^{2}(\pi x/L)}-{V}_{1}\,\cos (2\pi x/L)$$0\leqslant x\leqslant \tfrac{1}{2}L$$2{\sin }^{2}(\pi x/L)-1$
${V}_{0}+\displaystyle \frac{{V}_{+}-{V}_{-}\,\cosh (\lambda x)}{{\sinh }^{2}(\lambda x)}+{V}_{1}\,\cosh \left(\lambda x\right)$$x\geqslant 0$$\cosh \left(\lambda x\right)$

New window|CSV

In the conventional (textbook) formulation of quantum mechanics, the potential function plays a central role in providing physical information about the system. Hence, in the present formulation, the set of orthogonal polynomials replaces the potential function in this role. In fact, it is more than just that. As we shall see below, the orthogonal polynomials also carry kinematic information (e.g. the angular momentum) whereas the potential function does not.

Since the polynomials that are relevant to our work satisfy the general orthogonality relation (7) and have the asymptotic behavior (4) then Favard’s theorem [17] dictates that such polynomials must satisfy a three-term recursion relation of the form $\varepsilon {P}_{n}^{\mu }(\varepsilon )\,=$${a}_{n}^{\mu }{P}_{n}^{\mu }(\varepsilon )\,+{b}_{n-1}^{\mu }{P}_{n-1}^{\mu }(\varepsilon )\,+{c}_{n}^{\mu }{P}_{n+1}^{\mu }(\varepsilon ),$ where ${b}_{n}^{\mu }{c}_{n}^{\mu }\gt 0$ for all n. Of course, not all polynomials satisfy three-term recursion relations. In fact, some satisfy higher order recursions (e.g. four-term and five-term). However, such polynomials are not in the scope of our present study. Consequently, the three-term recursion relation dictates that the basis set $\left\{{\phi }_{n}(x)\right\}$ must produce a tridiagonal matrix representation for the corresponding wave operator. As such, the matrix wave equation becomes equivalent, and amenable to the said recursion relation. In this paper, we leave out technical details but include the most relevant information concerning orthogonal polynomials in the Appendices. Interested readers may consult cited references for the derivation of applicable results, especially [1, 3].

In the following section and due to the prime significance of the Wilson polynomial that sits at the top of the Askey tree of the hypergeometric class of orthogonal polynomials, we derive its asymptotics using the Darboux’s method and show that it agrees with the original work of Wilson [18] and with the second scenario of formula (4). In sections 3–5, we present several examples of orthogonal polynomials from the Askey scheme of the hypergeometric type and derive properties of the corresponding physical systems. Finally, we conclude in section 6 by making relevant comments and discussing related issues. Throughout the paper, we adopt the atomic units, $\hslash =M=1.$

2. Asymptotics of the Wilson polynomial

We write the Wilson polynomial in a different notation from that which is given by equation (1.1.1) on page 24 of [10] as follows$ \begin{eqnarray}\begin{array}{l}{\tilde{W}}_{n}^{\mu }({z}^{2};\nu ;a,b)=\displaystyle \frac{{(\mu +a)}_{n}{(\mu +b)}_{n}}{{(a+b)}_{n}n!}\\ \,\times \,{}_{4}F_{3}\left(\left.\begin{array}{l}-n,n+\mu +\nu +a+b-1,\mu +{\rm{i}}z,\mu -{\rm{i}}z\\ \,\,\,\mu +\nu ,\mu +a,\mu +b\end{array}\right|1\right),\end{array}\end{eqnarray}$where $z\geqslant 0.$ We can relate this notation to that in [10] as ${\tilde{W}}_{n}^{\mu }({z}^{2};\nu ;a,b)=$${W}_{n}({z}^{2};\mu ,\nu ,a,b)/n!\,{(\mu +\nu )}_{n}{(a+b)}_{n}.$ If $\mathrm{Re}(\mu ,\nu ,a,b)\gt 0$ and non-real parameters occur in conjugate pairs, then the generating function of these Wilson polynomials becomes (see equation (1.1.12) in [10])$ \begin{eqnarray}\begin{array}{l}\displaystyle \sum _{n=0}^{\infty }{\tilde{W}}_{n}^{\mu }({z}^{2};\nu ;a,b)\,{t}^{n}={}_{2}F_{1}\left(\left.\begin{array}{l}\mu +{\rm{i}}z,\nu +{\rm{i}}z\\ \,\mu +\nu \end{array}\right|t\right)\\ \,\,\,\,\,\,\,\,\,\,\times \,{}_{2}F_{1}\left(\left.\begin{array}{l}a-{\rm{i}}z,b-{\rm{i}}z\\ \,a+b\end{array}\right|t\right).\end{array}\end{eqnarray}$

Now, we apply the Darboux method to this generating function to obtain the asymptotics of ${\tilde{W}}_{n}^{\mu }({z}^{2};\nu ;a,b).$ The method is described in section 9 of chapter 8 in [19]. In addition, we also need the contiguous relation (equation (7) of Ex. 21 in [20])$ \begin{eqnarray}\begin{array}{l}\left(a+b-c\right){}_{2}F_{1}\left(\left.\begin{array}{l}a,b\\ \,\,c\end{array}\right|z\right)=a\left(1-z\right){}_{2}F_{1}\left(\left.\begin{array}{l}a+1,b\\ \,\,c\end{array}\right|z\right)\\ \,\,\,\,\,\,\,\,\,\,\,\,-\,\left(c-b\right){}_{2}F_{1}\left(\left.\begin{array}{l}a,b-1\\ \,\,c\end{array}\right|z\right),\end{array}\end{eqnarray}$and the Euler transformation [20]$ \begin{eqnarray}{}_{2}F_{1}\left(\left.\begin{array}{l}a,b\\ \,\,c\end{array}\right|z\right)={\left(1-z\right)}^{c-a-b}{}_{2}F_{1}\left(\left.\begin{array}{l}c-a,c-b\\ \,\,c\end{array}\right|z\right).\end{eqnarray}$

Moreover, we will also employ the Gauss sum (Theorem 18 in section 32 of [20])$ \begin{eqnarray}{}_{2}F_{1}\left(\left.\begin{array}{l}a,b\\ \,\,c\end{array}\right|1\right)=\tfrac{{\rm{\Gamma }}(c){\rm{\Gamma }}(c-a-b)}{{\rm{\Gamma }}(c-a){\rm{\Gamma }}(c-b)},\,\mathrm{Re}\left(c-a-b\right)\gt 0.\end{eqnarray}$

The contiguous relation (10) makes the first ${}_{2}F_{1}$ on the right side of the generating function (9) equal to$ \begin{eqnarray}\begin{array}{l}\displaystyle \frac{1}{2{\rm{i}}z}\left[\left(\mu +{\rm{i}}z\right)\left(1-t\right){}_{2}F_{1}\left(\left.\begin{array}{l}\mu +1+{\rm{i}}z,\nu +{\rm{i}}z\\ \,\,\,\mu +\nu \end{array}\right|t\right)\right.\\ \,-\,\left.\left(\mu -{\rm{i}}z\right){}_{2}F_{1}\left(\left.\begin{array}{l}\mu +{\rm{i}}z,\nu -1+{\rm{i}}z\\ \,\,\,\mu +\nu \end{array}\right|t\right)\right].\end{array}\end{eqnarray}$

We use the Euler transformation (11) to rewrite the first term inside the square bracket above as$ \begin{eqnarray}\left(\mu +{\rm{i}}z\right){\left(1-t\right)}^{-2{\rm{i}}z}{}_{2}F_{1}\left(\left.\begin{array}{l}\nu -1-{\rm{i}}z,\mu -{\rm{i}}z\\ \,\,\,\mu +\nu \end{array}\right|t\right).\end{eqnarray}$

Therefore, the first ${}_{2}F_{1}$ term on the right side of the generating function (9) becomes$ \begin{eqnarray}\begin{array}{l}\displaystyle \frac{{\left(1-t\right)}^{-{\rm{i}}z}}{2{\rm{i}}z}\left[\left(\mu +{\rm{i}}z\right){\left(1-t\right)}^{-{\rm{i}}z}{}_{2}F_{1}\left(\left.\begin{array}{l}\mu -{\rm{i}}z,\nu -1-{\rm{i}}z\\ \,\,\mu +\nu \end{array}\right|t\right)\right.\\ \,-\,\left.{\rm{complex}}\,{\rm{conjugate}}\Space{0ex}{3.2ex}{0ex}\right].\end{array}\end{eqnarray}$

In accordance with the Darboux method, the dominant term in a comparison function for the t-expansion of this factor becomes its limiting value at $t=1.$ The Gauss sum (12) evaluates the above ${}_{2}F_{1}$ at $t=1$ as $\tfrac{{\rm{\Gamma }}(\mu +\nu ){\rm{\Gamma }}(1+2{\rm{i}}z)}{{\rm{\Gamma }}(\mu +1+{\rm{i}}z){\rm{\Gamma }}(\nu +{\rm{i}}z)}=\tfrac{2{\rm{i}}z}{\mu +{\rm{i}}z}\tfrac{{\rm{\Gamma }}(\mu +\nu ){\rm{\Gamma }}(2{\rm{i}}z)}{{\rm{\Gamma }}(\mu +{\rm{i}}z){\rm{\Gamma }}(\nu +{\rm{i}}z)}.$ Thus, expression (15) near $t=1$ becomes$ \begin{eqnarray}\begin{array}{l}{\left(1-t\right)}^{-{\rm{i}}z}\left[{\left(1-t\right)}^{-{\rm{i}}z}\displaystyle \frac{{\rm{\Gamma }}(\mu +\nu ){\rm{\Gamma }}(2{\rm{i}}z)}{{\rm{\Gamma }}(\mu +{\rm{i}}z){\rm{\Gamma }}(\nu +{\rm{i}}z)}\right.\\ \,+\,\left.{\rm{complex}}\,{\rm{conjugate}}\right]\\ \,=\,{\left(1-t\right)}^{-{\rm{i}}z}\displaystyle \frac{{\rm{\Gamma }}(\mu +\nu )\left|{\rm{\Gamma }}(2{\rm{i}}z)\right|}{\left|{\rm{\Gamma }}(\mu +{\rm{i}}z){\rm{\Gamma }}(\nu +{\rm{i}}z)\right|}\\ \,\times \,\left[{\left(1-t\right)}^{-{\rm{i}}z}{{\rm{e}}}^{{\rm{i}}\alpha }+{\rm{c}}{\rm{.}}\,{\rm{c}}{\rm{.}}\Space{0ex}{3.0ex}{0ex}\right],\end{array}\end{eqnarray}$where $\alpha =\text{arg}\left[{\rm{\Gamma }}(2{\rm{i}}z)/{\rm{\Gamma }}(\mu +{\rm{i}}z){\rm{\Gamma }}(\nu +{\rm{i}}z)\right].$ Repeating the same treatment on the second ${}_{2}F_{1}$ on the right side of the generating function (9), we obtain the same result as (16) but with the replacement $z\to -z,$$\mu \to a$ and $\nu \to b$ giving$ \begin{eqnarray}\begin{array}{l}{\left(1-t\right)}^{+{\rm{i}}z}\left[{\left(1-t\right)}^{+{\rm{i}}z}\displaystyle \frac{{\rm{\Gamma }}(a+b){\rm{\Gamma }}(-2{\rm{i}}z)}{{\rm{\Gamma }}(a-{\rm{i}}z){\rm{\Gamma }}(b-{\rm{i}}z)}\right.\\ \,\,+\,\left.{\rm{complex}}\,{\rm{conjugate}}\right]\\ \,=\,{\left(1-t\right)}^{+{\rm{i}}z}\displaystyle \frac{{\rm{\Gamma }}(a+b)\left|{\rm{\Gamma }}(2{\rm{i}}z)\right|}{\left|{\rm{\Gamma }}(a+{\rm{i}}z){\rm{\Gamma }}(b+{\rm{i}}z)\right|}\\ \,\,\times \,\left[{\left(1-t\right)}^{-{\rm{i}}z}{{\rm{e}}}^{{\rm{i}}\beta }+{\rm{c}}{\rm{.}}\,{\rm{c}}{\rm{.}}\Space{0ex}{3.2ex}{0ex}\right],\end{array}\end{eqnarray}$where $\beta =\text{arg}\left[{\rm{\Gamma }}(2{\rm{i}}z)/{\rm{\Gamma }}(a+{\rm{i}}z){\rm{\Gamma }}(b+{\rm{i}}z)\right].$ Multiplication of (16) and (17) gives$ \begin{eqnarray}\begin{array}{l}\displaystyle \frac{{\rm{\Gamma }}(\mu +\nu ){\rm{\Gamma }}(a+b){\left|{\rm{\Gamma }}(2{\rm{i}}z)\right|}^{2}}{\left|{\rm{\Gamma }}(\mu +{\rm{i}}z){\rm{\Gamma }}(\nu +{\rm{i}}z){\rm{\Gamma }}(a+{\rm{i}}z){\rm{\Gamma }}(b+{\rm{i}}z)\right|}\\ \,\times \,\left[{\left(1-t\right)}^{-2{\rm{i}}z}{{\rm{e}}}^{{\rm{i}}(\alpha +\beta )}+{{\rm{e}}}^{{\rm{i}}(\alpha -\beta )}+{\rm{c}}{\rm{.}}\,{\rm{c}}{\rm{.}}\right].\end{array}\end{eqnarray}$Aside from t-independent factors, the comparison function near $t=1$ of the above expression is ${\left(1-t\right)}^{-2{\rm{i}}z}.$ The expansion of this term is$ \begin{eqnarray}{\left(1-t\right)}^{-2{\rm{i}}z}=\displaystyle \sum _{n=0}^{\infty }\displaystyle \frac{{\left({\rm{2i}}z\right)}_{n}}{{\rm{\Gamma }}(n+1)}{t}^{n}.\end{eqnarray}$

Therefore, applying the Darboux method to the generating function (9) gives the following asymptotics for the Wilson polynomial$ \begin{eqnarray}\begin{array}{l}{\tilde{W}}_{n}^{\mu }({z}^{2};\nu ;a,b)\\ \,\approx \,\displaystyle \frac{{\rm{\Gamma }}(\mu +\nu ){\rm{\Gamma }}(a+b)\left|{\rm{\Gamma }}(2{\rm{i}}z)\right|}{\left|{\rm{\Gamma }}(\mu +{\rm{i}}z){\rm{\Gamma }}(\nu +{\rm{i}}z){\rm{\Gamma }}(a+{\rm{i}}z){\rm{\Gamma }}(b+{\rm{i}}z)\right|}\\ \,\,\times \,\left[{n}^{{\rm{2i}}z-1}{{\rm{e}}}^{-{\rm{i}}\gamma }{{\rm{e}}}^{{\rm{i}}(\alpha +\beta )}+{\rm{c}}.\,{\rm{c}}.\right],\end{array}\end{eqnarray}$where $\gamma =\text{arg}\left[{\rm{\Gamma }}(2{\rm{i}}z)\right]$ and we have used ${\left(z\right)}_{n}=\tfrac{{\rm{\Gamma }}(n+z)}{{\rm{\Gamma }}(z)}$ and $\tfrac{{\rm{\Gamma }}(n+a)}{{\rm{\Gamma }}(n+b)}\approx {n}^{a-b}.$ Now, using $\text{arg}(a)+\text{arg}(b)=\text{arg}(ab),$$\text{arg}(a)-\text{arg}(b)=\text{arg}(a/b)$ and ${a}^{{\rm{i}}b}={{\rm{e}}}^{{\rm{i}}b\mathrm{ln}a},$ we obtain$ \begin{eqnarray}\begin{array}{l}{\tilde{W}}_{n}^{\mu }({z}^{2};\nu ;a,b)\approx \displaystyle \frac{2}{n}{\rm{\Gamma }}(\mu +\nu ){\rm{\Gamma }}(a+b)\left|{\mathscr{A}}({\rm{i}}z)\right|\\ \,\times \,\cos \left\{2z\,\mathrm{ln}(n)+\text{arg}\left[{\mathscr{A}}({\rm{i}}z)\right]\right\}+O({n}^{-1}),\end{array}\end{eqnarray}$where ${\mathscr{A}}(z)={\rm{\Gamma }}(2z)/{\rm{\Gamma }}(\mu +z){\rm{\Gamma }}(\nu +z){\rm{\Gamma }}(a+z){\rm{\Gamma }}(b+z).$
This finding agrees with Wilson’s result [18]22There is a typo in [18] by which the+sign inside the argument of the cosine in equation (21) is replaced by a−sign (private communication with Wilson on July 2016)., which was obtained using a convexity argument which is especially well suited to estimating certain hypergeometric series and their integral analogs. The asymptotics of the orthonormal version of the polynomial, which is given in appendix C by equation (C1), could easily be obtained from (21) as$ \begin{eqnarray}\begin{array}{l}{W}_{n}^{\mu }({z}^{2};\nu ;a,b)\approx B(\mu ,\nu ,a,b)\sqrt{\displaystyle \frac{2}{n}}\left\{2\left|{\mathscr{A}}({\rm{i}}z)\right|\right.\\ \,\times \,\left.\cos \left[2z\,\mathrm{ln}\,n+\text{arg}{\mathscr{A}}({\rm{i}}z)\right]+O({n}^{-1})\right\},\end{array}\end{eqnarray}$where
$ \begin{eqnarray*}B(\mu ,\nu ,a,b)=\sqrt{\tfrac{{\rm{\Gamma }}(\mu +\nu ){\rm{\Gamma }}(a+b){\rm{\Gamma }}(\mu +a){\rm{\Gamma }}(\mu +b){\rm{\Gamma }}(\nu +a){\rm{\Gamma }}(\nu +b)}{{\rm{\Gamma }}(\mu +\nu +a+b)}}\end{eqnarray*}$
and we have again used the asymptotic identity $\tfrac{{\rm{\Gamma }}(n+a)}{{\rm{\Gamma }}(n+b)}\approx {n}^{a-b}.$ Comparing this asymptotics to formula (4) leads to the following scattering amplitude and phase shift$ \begin{eqnarray}{A}^{\mu }(\varepsilon )=\displaystyle \frac{2\sqrt{2}B(\mu ,\nu ,a,b)}{\left|{\rm{\Gamma }}(\mu +{\rm{i}}z){\rm{\Gamma }}(\nu +{\rm{i}}z){\rm{\Gamma }}(a+{\rm{i}}z){\rm{\Gamma }}(b+{\rm{i}}z)/{\rm{\Gamma }}(2{\rm{i}}z)\right|},\end{eqnarray}$$ \begin{eqnarray}\begin{array}{l}{\delta }^{\mu }(\varepsilon )=\text{arg}{\rm{\Gamma }}(2{\rm{i}}z)-\text{arg}\left[{\rm{\Gamma }}(\mu +{\rm{i}}z){\rm{\Gamma }}(\nu +{\rm{i}}z)\right.\\ \,\,\,\times \,\left.{\rm{\Gamma }}(a+{\rm{i}}z){\rm{\Gamma }}(b+{\rm{i}}z)\right].\end{array}\end{eqnarray}$

The scattering amplitude shows that a discrete finite spectrum occur if $\mu +{\rm{i}}z=-k,$ where $k=0,1,2,\mathrm{..},N$ and N is the largest integer less than or equal to $-\mu .$ Thus, the spectrum formula associated with the Wilson polynomial is$ \begin{eqnarray}{z}_{k}^{2}=-{(k+\mu )}^{2}.\end{eqnarray}$

Despite the fact that the Meixner–Pollaczek and continuous dual Hahn polynomials are limiting cases of the Wilson polynomial, their asymptotics could not be obtained using the asymptotics of the Wilson polynomial calculated above simply because we cannot interchange this limit with the asymptotic limit. Therefore, those asymptotics have to be derived independently. This is done elsewhere (see, for example, the appendix in [1]).

3. The Meixner–Pollaczek polynomial class of problems

In this section, we consider the class of quantum systems whose continuum scattering states are described by the wavefunction (3) where the expansion coefficients are the Meixner–Pollaczek polynomials ${P}_{n}^{\mu }(z,\theta )$ shown in appendix A by equation (A1).

3.1. The Coulomb problem

We start by choosing the physical parameters in the polynomial as: $\mu ={\ell }+1,$$\cos \,\theta =$$\tfrac{2E-{(\lambda /2)}^{2}}{2E+{(\lambda /2)}^{2}}$ and $z=Z/\sqrt{2E},$ where ${\ell }$ is a non-negative integer, Z is a real number and λ is a positive length scale parameter. Depending on the range of values of the physical parameters, this system can have a continuous or discrete energy spectrum. For example, if E is negative then z becomes pure imaginary and $\left|\cos \,\theta \right|\gt 1.$ As explained in appendix A, this is equivalent to the replacement of $z\to {\rm{i}}z$ and $\theta \to {\rm{i}}\theta ,$ which turns the Meixner–Pollaczek polynomials into its discrete version, the Meixner polynomials. Substituting these parameters in equation (A7) gives the following scattering phase shift$ \begin{eqnarray}\delta (E)=\text{arg}{\rm{\Gamma }}\left({\ell }+1+{\rm{i}}Z/\sqrt{2E}\right).\end{eqnarray}$Whereas, substitution in equation (A8) gives the following energy spectrum formula$ \begin{eqnarray}{E}_{k}=-\displaystyle \frac{1}{2}\displaystyle \frac{{Z}^{2}}{{\left(k+{\ell }+1\right)}^{2}}.\end{eqnarray}$

These results are identical to those of the well-known Coulomb interaction in three dimensions, $V(r)=Z/r,$ where r is the radial coordinate, Z is the electric charge and ${\ell }$ is the angular momentum quantum number. The discrete polynomial in the bound states wavefunction expansion (5) is the Meixner polynomial shown in appendix A by equation (A9). The requirement on the basis to support a tridiagonal matrix representation of the Schrödinger wave operator, $-\tfrac{1}{2}\tfrac{{{\rm{d}}}^{2}}{{\rm{d}}{r}^{2}}+\tfrac{{\ell }({\ell }+1)}{2{r}^{2}}+\tfrac{Z}{r}-E,$ gives (see section II.A.2 in [3])$ \begin{eqnarray}{\phi }_{n}(r)=\sqrt{\tfrac{{\rm{\Gamma }}(n+1)}{{\rm{\Gamma }}(n+2{\ell }+2)}}\,{(\lambda r)}^{{\ell }+1}{{\rm{e}}}^{-\lambda r/2}{L}_{n}^{2{\ell }+1}(\lambda r),\end{eqnarray}$where ${L}_{n}^{\nu }(z)$ is the Laguerre polynomial. This is not the only exactly solvable problem in the Meixner–Pollaczek polynomial class. Next, we give two other examples; one with only an infinite bound states spectrum and another with a continuous as well as discrete finite spectrum (however, in this class it is exactly solvable only for the discrete bound states).

3.2. The oscillator problem

For the second example in this class, we take the polynomial parameters as $\mu =\tfrac{1}{2}\left({\ell }+\tfrac{3}{2}\right)$ and $z={\rm{i}}E/2\omega ,$ where ω is a real number. The choice of z as pure imaginary mandates the replacement $\theta \to {\rm{i}}\theta $ so that reality is maintained for the polynomial (A1) and its recursion relation (A4). As seen in appendix A, this leads only to bound states and the discrete Meixner polynomial. The infinite energy spectrum formula is obtained from equation (A8) as$ \begin{eqnarray}{E}_{k}=\omega \left(2k+{\ell }+\tfrac{3}{2}\right).\end{eqnarray}$

This corresponds to the well-known energy spectrum of the isotropic oscillator $V(r)=\tfrac{1}{2}{\omega }^{2}{r}^{2}$ with oscillator frequency ω. The corresponding eigenstates are written as in equation (5) in terms of the discrete Meixner polynomial. Moreover, the requirement that the basis yield a tridiagonal matrix representation for the wave operator, $-\tfrac{1}{2}\tfrac{{{\rm{d}}}^{2}}{{\rm{d}}{r}^{2}}+$$\tfrac{{\ell }({\ell }+1)}{2{r}^{2}}\,+\tfrac{1}{2}{\omega }^{2}{r}^{2}-E,$ gives (see section II.A. 1 in [3])$ \begin{eqnarray}{\phi }_{n}(r)=\sqrt{\tfrac{{\rm{\Gamma }}(n+1)}{{\rm{\Gamma }}(n+{\ell }+3/2)}}\,{(\lambda r)}^{{\ell }+1}{{\rm{e}}}^{-{\lambda }^{2}{r}^{2}/2}{L}_{n}^{{\ell }+1/2}({\lambda }^{2}{r}^{2}),\end{eqnarray}$

where λ is a length scale parameter such that ${\lambda }^{2}\leqslant 4\omega .$ Additionally, the parameter β in the Meixner polynomial (A9) is obtained as $\beta ={{\rm{e}}}^{-2\theta }$ where $\cosh \,\theta =\tfrac{{\omega }^{2}+{(\lambda /2)}^{4}}{{\omega }^{2}-{(\lambda /2)}^{4}}.$

3.3. The Morse problem

The final problem in the Meixner–Pollaczek polynomial class corresponds to the parameter assignments: $\mu =\tfrac{1}{2}+\sqrt{-\varepsilon }$ and $z={\rm{i}}{u}_{1}/2\sqrt{{u}_{0}},$ where all parameters are real with $\varepsilon \lt 0,$${u}_{0}\gt 0$ and ${u}_{1}\lt 0.$ Thus, it is required that $\theta \to {\rm{i}}\theta $ in equations (A1) and (A4) turning the polynomial into one of its two discrete versions. Formula (A8) gives the energy spectrum as$ \begin{eqnarray}{\varepsilon }_{k}=-{\left(k+\tfrac{1}{2}+{u}_{1}/2\sqrt{{u}_{0}}\right)}^{2},\end{eqnarray}$where $k=0,1,\mathrm{..},N$ and N is the largest integer less than or equal to $-{u}_{1}/2\sqrt{{u}_{0}}-\tfrac{1}{2}.$ If we introduce an inverse length parameter α and write $E=\tfrac{1}{2}{\alpha }^{2}\varepsilon $ and ${V}_{i}=\tfrac{1}{2}{\alpha }^{2}{u}_{i},$ then we can rewrite the spectrum formula (31) as follows$ \begin{eqnarray}{E}_{k}=-\tfrac{1}{2}{\alpha }^{2}{\left(k+\tfrac{1}{2}+{V}_{1}/\alpha \sqrt{2{V}_{0}}\right)}^{2}.\end{eqnarray}$

This, in fact, is the energy spectrum formula of the one-dimensional Morse potential $V(x)={V}_{0}{{\rm{e}}}^{2\alpha x}+{V}_{1}{{\rm{e}}}^{\alpha x}$ where $-\infty \lt x\lt +\infty $ [21]. The corresponding bound states are written as in equation (5) in terms of the discrete version of the Meixner–Pollaczek polynomial with a finite spectrum, which is the Krawtchouk polynomial not the Meixner polynomial. The orthonormal version of this polynomial is given in appendix A by equation (A11). The requirement that the corresponding basis gives a tridiagonal matrix representation for the wave operator, $-\tfrac{1}{2}\tfrac{{{\rm{d}}}^{2}}{{\rm{d}}{x}^{2}}+{V}_{0}{{\rm{e}}}^{2\alpha x}+{V}_{1}{{\rm{e}}}^{\alpha x}-E,$ results in (see section II.A.3 in [3])$ \begin{eqnarray}{\phi }_{n}(x)=\sqrt{\tfrac{{\rm{\Gamma }}(n+1)}{{\rm{\Gamma }}(n+\nu +1)}}\,{y}^{\nu /2}{{\rm{e}}}^{-y/2}{L}_{n}^{\nu }(y),\end{eqnarray}$where $y(x)={{\rm{e}}}^{\alpha x}$ and $\nu =\tfrac{2}{\alpha }\sqrt{-2E}.$ Moreover, the parameter $\gamma ={{\rm{e}}}^{-2\theta }$ in the Krawtchouk polynomial is obtained from $\cosh \,\theta =\tfrac{2{V}_{0}+{(\alpha /2)}^{2}}{2{V}_{0}-{(\alpha /2)}^{2}}$ with ${V}_{0}\geqslant {\alpha }^{2}/8.$

In this class, we were able to obtain full solutions for two problems, the Coulomb and the isotropic oscillator. The latter has only discrete bound states whereas the former has both discrete bound states as well as continuum scattering states. Additionally, we were able to obtain only partial solution to the 1D Morse oscillator. We could find only the discrete bound states solution but not the continuum scattering states. In the following section, we will remedy that.

4. The continuous dual Hahn polynomial class of problems

In this section, we consider the class of problems whose continuum scattering states are described by the wavefunction (3) where the expansion coefficients are the continuous dual Hahn polynomial ${S}_{n}^{\mu }({z}^{2};a,b)$ shown in appendix B by equation (B1). Throughout this section, we restrict our investigation to the special case where the two polynomial parameters a and b are equal.

闂備胶绮崝妤呭箠閹捐鍚规い鏂垮⒔閸楁岸鎮楅敐搴濈盎缂佷緡鍣i弻鐔煎礂閼测晝鐓傞梺绋跨焿閹凤拷2濠电偞鍨堕幐绋棵洪妸鈺嬬稏闁圭儤顨嗛崵鈧梺鍛婂姦娴滅偤宕洪敓鐘崇厽闁靛繈鍊栧☉褔鏌i埄鍐噰闁诡啫鍥ч唶闁挎繂娲㈤崑鎺楁⒑閸濆嫬鈧綊鎮锋潏鈺傤潟闁跨噦鎷�
濠电姰鍨归悥銏ゅ炊瑜嶆慨銈夋⒑閸涘﹤绗掓俊顐g洴椤㈡棃濮€閵堝棭妫勯柣搴秵閸犳牠宕㈤幘顔界厸闁告洟娼ч悘锝嗐亜閹存繃澶勭紒瀣樀閸┾偓妞ゆ巻鍋撻柍璇查叄濡鹃亶鏌嶈閸撴瑩宕导瀛樺亯婵炲樊浜濋弲顒勬倶閻愮數鎽傞柛銈囧Т闇夋繝濠傚暣椤庢銇勯埞顓炲婵挳鏌¢崶鈺佹灁闁告瑢鍋撻梻浣哥秺濞佳嗐亹閻愮數绠旈柟鎯ь嚟閳绘梹鎱ㄥΟ璇插闁搞倧绠撻弻鐔虹矙閹稿孩鎮欓梺浼欑秮缁犳牕顕i鈶╂瀻闁归偊鍘剧粙鍕⒑閹稿海鈽夐柡鍫墴瀹曞綊濡歌婵ジ鏌涘☉姗堟敾缂佺姵甯為埀顒€鐏氬姗€鎮ч崱娴板洭宕稿Δ浣镐痪闂佺鎻梽鍕晬閺嶎厽鐓忛柛鈩冩礀椤b暜ee濠电姰鍨圭紞濠囧焵椤掍胶鈯曢柕鍡楀暣閺屾盯骞掗幋鐑嗘濡炪倖甯為崰鎰矙婵犲洦鍋愰柣銏㈡暩鏁堥梻浣稿悑濠㈡﹢宕导瀛樺亯闁告繂濯辨惔銏$秶妞ゆ劗鍠庢禍楣冩煛閸ャ劍鐨戦柣鐔叉櫅閳藉骞樼紙鐘卞濡炪倖娲濆▍鏇炨缚韫囨稑宸濇い鎾楀啯顔�20婵°倗濮烽崑鐘诲箵椤忓棙顫曟繝闈涱儏缁犳垿鏌ゆ慨鎰偓妤€鈻旈姀鐘嗙懓饪伴崘鈺婃%缂備礁顦顓㈠焵椤掆偓濠€閬嶅磻閻旂厧鏋侀柕鍫濐槹閸庡秹鏌涢弴銊ュ闁伙箑缍婇幃妤冩喆閸曨収鏆¢梺鍝勬閸嬫捇姊洪崫鍕垫Ч闁告梹鐗犻幃锟犳晬閸曨剙鐝伴梺闈涚箚閸撴繈鎮″▎鎰濠㈣泛顑嗙粈鈧悗娈垮櫍閺€鍗烆嚗閸曨偒鍚嬮柛鏇ㄥ幘濡叉垿姊洪崫鍕偓浠嬶綖婢跺本鍏滈柛顐f礃閺咁剟鎮橀悙闈涗壕缂佺姵甯″濠氬炊閿濆懍澹曢梺鑽ゅ枑濞叉垿鎳楃捄琛℃灁闁硅揪闄勯崕鎴︽倵閿濆骸骞樼紒鐘崇墵閺屸剝寰勫☉娆忣伓
相关话题/Representation quantum mechanical

闁诲孩顔栭崰鎺楀磻閹剧粯鈷戞い鎰剁悼椤e弶绻濋埀顒勫箥椤旀儳宕ュ┑顔筋殘椤︾硽闂備焦瀵х粙鎺楁儗椤斿墽鍗氶悗娑欋缚閳绘棃鏌涢妷銏℃珦婵炲吋宀搁弻銈夊级閸喗娈堕梺缁樼墱閸庛倗绮欐径鎰闁肩⒈鍓涢幊婵嬫煟閻樺弶澶勭€规洘锕㈤敐鐐碘偓锝庡亾缁憋綁鏌熸潏楣冩妞ゆ挾鎳撻妴鎺戭潩閾忣偆銆婂┑鐘亾闁告稑鐡ㄩ弲顒勬倶閻愭彃鈷旀俊灞傚姂閺岋繝宕奸銏犲箰缂備焦顨呴ˇ闈涚暦濮橆叏绱eù锝勮娴煎洭姊洪崨濠傜瑨婵☆偅绻堥獮鎰板醇閺囩喐娅栨繛杈剧秬濞咃綁宕″⿰鍫熷€甸柣銏犵仛閸も偓闂佹悶鍔嶇换鍕垝鐠囧弬鏃傗偓锝庡墰閿涳拷
547闂備礁婀遍。浠嬪磻閹剧粯鈷掗柛鏇楁櫅閻忣亪鏌eΔ鈧柊锝夊箠閹捐绀冩い蹇撴閻撴盯姊洪崗鍏肩凡闁哥噥鍋勯悾鐑芥晸閿燂拷1130缂傚倷绀侀ˇ顖滅矓閻㈢鍋撻崹顐g殤闁逞屽墲椤鍠婂澶婃辈闁逞屽墴閺屸剝寰勭€n亜顫庡┑鐐茬墛閸ㄥ灝鐣烽敓鐘茬鐟滃繒绮欓崶鈺冪<濠㈣泛锕︽晥闂佸憡菧閸婃牜缂撻挊澹╂棃宕担瑙勭槣闂佸湱鍘ч悺銊╁箰閸洖鐒垫い鎴炲缁佺増銇勯銏╁剱闁挎稒鍔欓獮瀣敍濠婂拋妲锋繝鐢靛仦閸ㄥ綊寮粙妫电儤绻濋崶銊ユ闁哄鐗滈崑鎺楀吹閺冨牊鐓忛柛鈩冩礉閸忓瞼绱掗鍏夊亾鐡掍浇顫夐幆鏂库槈閹烘垳澹曟繛杈剧悼閺咁偄危閸儲鐓曢柟鐑樻尰閸嬬娀鏌嶈閸忔稓娆㈠璺洪棷濡わ絽鍟幊姘扁偓骞垮劚閸熺娀宕戦幘瀛樺闁绘垶锚閳ь剛鍋熼埀顒冾潐閹爼宕曢鐐茬劦妞ゆ垼鍎婚崗灞俱亜閹惧瓨鍊愰柟顔肩埣瀹曢亶骞囬妸銉ゅ婵炶揪绲炬禍鑺ョ閿曗偓闇夐柛蹇曞帶閹兼悂鏌嶈閸忔稑霉閸ヮ剙纾奸柕濠忕畱椤曡鲸鎱ㄥΟ绋垮姉闁稿鎸婚幏鍛喆閸曨剛鏆氶梻浣哄帶瀵儼銇愰崘顏嗙处濡わ絽鍟崑鐘绘煕閳╁啫濮€闁稿鎸婚幏鍛存偪椤栨艾绠戦梻浣告惈閸婄ǹ煤閵忋倕鐒垫い鎴炲缁佹澘顭跨憴鍕磳鐎殿喚鏁婚、娑樜熷畡棰佸婵炶揪缍€椤鎮¢埀顒勬⒒閸屾艾鈧粙顢欐繝鍕潟闁割偅娲栫粻缁樸亜閹炬潙顥氶柛瀣尰閹峰懘宕烽婧惧亾婵犲洦鍊垫繛鎴濈枃缁€瀣煃瑜滈崗娑氱矆娴h桨鐒婇柟娈垮枓閸嬫挸鈽夌€圭姷顦伴梺閫炲苯鍘告繛鏉戞喘椤㈡﹢宕妷褌绗夊┑掳鍊撻悞锔捐姳濮樿埖鐓忛柛鈩冩礈椤︼妇鈧湱枪椤嘲鐣烽敐鍥︽勃闁稿本顨呮禍鎯归敐鍛暈闁告洟绠栭弻锝夋倷閸欏妫戦梺閫炲苯鍘搁柣鎺炵畵瀵剟宕掑锝嗙參濠殿喚鎳撳ú鐘诲磻閹惧瓨濯撮柛娑橈攻閸f悂鏌f惔銏犲枙閻犳劗鍠栭崺鈧い鎴炲椤﹂绱撳鍜佸剶闁硅櫕鐗犻幊鐘活敆閸愮偓钑夌紓鍌欑劍閸愬骞忛敓锟�28缂傚倷绶¢崑澶愵敋瑜旈、妤呮偄閾忓湱鐓嬮梺瑙勬儗閸ㄥ磭澹曢敓锟�