SciELO - Scientific Electronic Library Online

 
vol.62 número6Nanocrystals of Ni doped ZnO semiconductor by Sol-Gel combustion methodComparison of variational solutions of the Thomas-Fermi model in terms of the ionization energy índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • No hay artículos similaresSimilares en SciELO

Compartir


Revista mexicana de física

versión impresa ISSN 0035-001X

Rev. mex. fis. vol.62 no.6 México nov./dic. 2016

 

Review

Solutions of the Schrödinger equation given by solutions of the Hamilton-Jacobi equation

G.F. Torres del Castillo1 

C. Sosa Sánchez2 

1Departamento de Física Matemática, Instituto de Ciencias Universidad Autónoma de Puebla, 72570 Puebla, Pue., México

2Facultad de Ciencias Físico Matemáticas Universidad Autónoma de Puebla, 72570 Puebla, Pue., México.


Abstract

We find the form of the potential depending on the coordinates and the time such that a solution, S, of the Hamilton-Jacobi equation yields an exact solution, exp⁡(iS/ℏ), of the corresponding Schrödinger equation.

Keywords: Schrödinger equation; Hamilton-Jacobi equation

Resumen

Hallamos la forma del potencial dependiente de las coordenadas y el tiempo tal que una solución, S, de la ecuación de Hamilton-Jacobi produce una solución exacta, exp(iS=~), de la ecuación de Schrödinger correspondiente.

Descriptores: Ecuación de Schrödinger; ecuación de Hamilton-Jacobi

1. Introduction

As is well known, the Hamilton-Jacobi (HJ) equation can be regarded as an approximation to the Schrödinger equation. More precisely, these two equations are related in the following way. If ψ=exp⁡(i𝑆/ℏ) is substituted into the Schrödinger equation for a particle of mass 𝑚 and electric charge 𝑒, subject to a velocity-independent potential 𝑉 and a magnetic field with vector potential A.

-22m(-iecA)2ψ+Vψ=iψt,

imposing the condition ∇⋅A=0, one obtains

12m(S-ecA)2+V+St=i2m2S,

which is equivalent to ([sch]). Then, if one looks for a solution of this last equation of the form

S=S0+iS1+(i)2S2+,

one finds that 𝑆 0 satisfies the HJ equation

12m(S0-ecA)2+V+S0t=0

(see, e.g., Ref. , Sec. 2.8). However, Eq. ([equiv]) shows that if 𝑆 satisfies the Laplace equation, 𝛻 2 𝑆=0, and the HJ equation (3), then 𝜓=exp⁡(i𝑆/ℏ) is an exact solution of the Schrödinger equation (1).

The aim of this paper is to give some explicit examples of this relationship between solutions of the HJ equation and solutions of the corresponding Schrödinger equation. In Sec. 2 we show that from each solution, 𝑆, of the Laplace equation, and any given magnetic field, one can define a potential, 𝑉, in such a way that 𝑆 is a solution of the HJ equation and exp⁡(i𝑆/ℏ) is a solution of the corresponding Schrödinger equation. In Secs. 3-5 we present some examples, restricting ourselves to problems without magnetic field. An example of this class was given in Ref. , where the attention was focused on the HJ equation.

2. Basic results

According to Eq.(2), assuming that 𝑆 is a solution of the Laplace equation

2S=0,

satisfies the HJ equation (3) if and only if 𝜓=exp⁡(i𝑆/ℏ) is a solution of the Schrödinger equation ([sch]). Instead of specifying the potentials 𝑉 and 𝐀, and then try to find a simultaneous solution of the Laplace equation and the HJ equation, it is simpler to start from a solution, 𝑆, of the Laplace equation and define the potential 𝑉 so that 𝑆 automatically satisfies the HJ equation (3). That is, to take

V=-St-12m(S-ecA)2

or, if the magnetic field is absent, choosing A=0,

V=-St-12m(S)2.

Of course, the resulting potential may be of little interest, but if one makes use of the general solution of the Laplace equation, one can find all the potentials for which a solution 𝑆 of the HJ equation corresponds to a solution exp⁡(i𝑆/ℏ) of the Schrödinger equation.

As pointed out after Eq. (3), if 𝑆 is any solution of the HJ equation (3) and, simultaneously, of the Laplace equation, then exp⁡(i𝑆/ℏ) is a solution of the corresponding Schrödinger equation. What would produce only a semiclassical approximation leads, in fact, to an exact solution of the Schrödinger equation. Since the modulus of exp⁡(i𝑆/ℏ) is equal to 1, all these wavefunctions cannot be normalized; however, just like the eigenfunctions of the linear momentum, which are also non-normalizable, in some cases they can be used to express any normalizable wavefunction (see also the comment at the end of Sec. 3.1). In fact, if one finds a complete solution of the HJ equation, 𝑆, then exp⁡(i𝑆/ℏ) possesses the appropriate number of parameters to form a complete set of functions in terms of which any wavefunction can be expanded (see the examples below). (It is interesting to note that the word “complete” appears here with two different, but related, meanings.) The solutions of the Schrödinger equation obtained here differ from the solutions usually found because the latter are commonly obtained by separation of variables (for instance, the wavefunction (11) is not separable and does not represent a state with a well-defined energy).

In the following sections we give some explicit examples, without magnetic field, starting from the highly illustrative one-dimensional case.

3. Examples in one dimension

In the one-dimensional case, the general solution of Eq. (4) is a linear function of the Cartesian coordinate

S(x,t)=α(t)x+β(t),

where 𝛼(𝑡) and ??(𝑡) are real-valued functions of 𝑡 only. Substituting (7) into (6) we get

V=-α'(t)x-β'(t)-12mα2(t),

which corresponds to a possibly time-dependent uniform field of force. Apart from the trivial case of a free particle, two specific examples contained in (8) are a time-independent uniform field, and a uniform field with an intensity that changes linearly with time.

3.1 Constant uniform field

From Eq.(8) we see that in order to have a potential

V=-Fx,

where 𝐹 is a constant, corresponding to a constant uniform force of magnitude 𝐹, 𝛼′(𝑡) must be equal to 𝐹 and 𝛽′(𝑡)=− 𝛼 2 (𝑡)/2𝑚, hence

α(t)=Ft+P,β(t)=-(Ft+P)36mF,

where 𝑃 is an arbitrary constant. Thus, substituting these expressions into (7), we obtain

S(x,t)=(Ft+P)x-(Ft+P)36mF.

The solution (10) of the HJ equation corresponding to the potential (9) is complete (by virtue of the presence of the arbitrary parameter P and it is not separable (because of the term 𝐹𝑡𝑥) and, for each value of 𝑃, the wavefunction

ψ(x,t)=expi[(Ft+P)x-(Ft+P)36mF]

is a solution of the Schrödinger equation with the potential (9), which is not (multiplicatively) separable (and is not a stationary state). The wavefunctions (11) are non-normalizable, just as the “plane waves” exp⁡(i𝑃𝑥/ℏ), and, as the latter, for each value of 𝑡, form a complete set, in terms of which any wavefunction can be expanded (by means of an integral on the parameter 𝑃). In fact, in the same manner as the plane waves are eigenfunctions of the momentum operator, the wavefunctions (11) are eigenfunctions of the (time-dependent) operator 𝑝−𝐹𝑡 (which is a conserved quantity), with eigenvalue 𝑃 .

3.2 Field with constant growth

We now consider the time-dependent potential

V(x)=-ktx,

where 𝑘 is a constant. This potential coincides with (8) if 𝛼′(𝑡)=𝑘𝑡 and 𝛽′(??)=− 𝛼 2 (𝑡)/2𝑚, hence

α(t)=kt22+P,β(t)=-12mk2t520+kPt33+P2t,

P is an arbitrary constant. Thus,

S(x,t)=kt22+Px-12mk2t520+kPt33+P2t

is a complete solution of the HJ equation with potential (12) and

ψ(x,t)=expi[(kt22+P)x-12m(k2t520+kPt33+P2t)]

is a solution of the corresponding Schrödinger equation.

The (non-normalizable) wavefunctions (14) form a complete set. In fact, (14) is eigenfunction of the conserved operator 𝑝−(1/2)𝑘 𝑡 2 , with eigenvalue 𝑃.

4. Examples in two dimensions

The solutions of the Laplace equation in two dimensions are the form

S(x,y,t)=Re[f(z,t)],

where 𝑧≡𝑥+i𝑦, and 𝑓(𝑧,𝑡) is an analytic function of 𝑧 that may also depend on 𝑡 and on other parameters. Substituting (15) into Eq. (6) we obtain the expression for the potential

V=-12(f+f¯)t-12mfz2.

In the following subsections we give two explicit examples corresponding to specific choices of 𝑓.

4.1 The repulsive isotropic harmonic oscillator

The function

f(z,t)=mω2z2+(P1e-ωt-iP2eωt)z+14mω(P1e-ωt-iP2eωt)2,

where 𝑚, 𝜔, 𝑃 1 , and 𝑃 2 are constants, is an analytic function of 𝑧 and, from Eq. ([pot2d]), we obtain the potential

V=-mω22(x2+y2),

which corresponds to a “repulsive harmonic oscillator” (see, e.g., Ref. 4). Then, making use of Eq. (15), we have the complete solution of the HJ equation

S(x,y,t)=mω2(x2-y2)+P1e-ωtx+P2eωty+14mω(P12e-2ωt-P22e2ωt)

and the solutions of the Schrödinger equation

ψ(x,y,t)=expi2mω2(x2-y2)+P1e-ωtx+P2eωty+14mω(P12e-2ωt-P22e2ωt)1.

Each wavefunction (19) is a common eigenfunction of the conserved operators e 𝜔𝑡 ( 𝑝 𝑥 −𝑚𝜔𝑥) and e −𝜔𝑡 ( 𝑝 𝑦 +𝑚𝜔𝑦) with eigenvalues 𝑃 1 and 𝑃 2 , respectively.

4.2 The central potential 𝑽=−𝒌/ 𝒓 𝟐

The function

f(z,t)=2mklnz,

where 𝑘 is a constant, is an analytic function of 𝑧 (except at 𝑧=0) and leads to the central potential

V(x,y)=-kx2+y2.

Thus,

S=mk2ln(x2+y2)

is a solution of the HJ equation with the potential (20). By contrast with the previous examples, this solution does not contain arbitrary parameters and does not depend on the time.

5. Examples in higher dimensions

Apart from the potentials that can be defined directly with the aid of Eq. (6), making use of a solution of the Laplace equation, as in the foregoing sections, we can also combine the potentials already obtained in low dimensions to construct potentials in two or more dimensions. For instance, from Eqs. (10) and (18) we can form the function

Sx,y,z,t)=mω2(x2-y2)+P1e-ωtx+P2eωty+14mω(P12e-2ωt-P22e2ωt)+(Ft+P3)z-(Ft+P3)36mF

which satisfies the Laplace equation in three dimensions and leads to the potential [cf. Eqs. (9) and (17)]

V(x,y,z)=-mω22(x2+y2)-Fz

and, therefore, exp⁡(i𝑆/ℏ) is a solution of the Schrödinger equation with this last potential, which contains three quantum numbers 𝑃 1 , 𝑃 2 , and 𝑃 3 .

6. Concluding remarks

It may be remarked that, for a given potential, the existence of a solution of the HJ equation that also satisfies the Laplace equation, does not imply that all solutions of this HJ equation will satisfy the Laplace equation. In the examples of Section 3, the principal functions (10) and (13) (which contain an arbitrary parameter, 𝑃) are, up to an additive constant, the only common solutions of the corresponding HJ equation and the Laplace equation.

One open question is what are the physical or geometrical properties of a potential such that the corresponding HJ equation admits solutions that also satisfy the Laplace equation. Another question is the meaning of the quantum states obtained in this manner.

The eikonal equation of geometrical optics can be regarded as an approximation to the scalar Helmholtz equation, in a similar manner as the HJ equation is an approximation to the Schrödinger equation and, therefore, one can expect that, in some cases, a solution of the eikonal equation would produce an exact solution of the Helmholtz equation. This connection is currently under investigation5.

Acknowledgements

The authors acknowledge Dr. Gilberto Silva Ortigoza for useful comments and to the referee for calling their attention to some points not clearly presented in the original version. One of the authors (C.S.S.) wishes to thank the Consejo Nacional de Ciencia y Tecnología for financial support.

REFERENCES

1. K. Gottfried and T.-M. Yan, Quantum Mechanics: Fundamentals, 2nd ed. (Springer-Verlag, New York, 2003). [ Links ]

2. G.F. Torres del Castillo, Rev. Mex. Fís. 59 (2013) 478. [ Links ]

3. G.F. Torres del Castillo, Rev. Mex. Fís. 61 (2015) 376. [ Links ]

4. K.B. Wolf, Rev. Mex. Fís. E 56 (2010) 83. [ Links ]

5. Omar de Jesús Cabrera-Rosas, Ernesto Espíndola-Ramos, Salvador Juárez-Reyes, Israel Julián-Macías, Paula Ortega-Vidals, Gilberto Silva-Ortigoza, Ramón Silva-Ortigoza, and Citlalli Sosa-Sánchez 2014 (submitted for publication). [ Links ]

Received: May 16, 2016; Accepted: August 15, 2016

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License