Skip to main content

ORIGINAL RESEARCH article

Front. Phys., 15 April 2020
Sec. Statistical and Computational Physics
This article is part of the Research Topic Contact Interactions in Quantum Mechanics: Theory, Mathematical Aspects and Applications View all 16 articles

The Propagators for δ and δ′ Potentials With Time-Dependent Strengths

  • 1Department of Mathematics, İzmir Institute of Technology, İzmir, Turkey
  • 2Departamento de Física Teórica, Atómica y Óptica and IMUVA, Universidad de Valladolid, Valladolid, Spain
  • 3Department of Physics, Adnan Menderes University, Aydın, Turkey

We study the time-dependent Schrödinger equation with finite number of Dirac δ and δ′ potentials with time dependent strengths in one dimension. We obtain the formal solution for generic time dependent strengths and then we study the particular cases for single delta potential and limiting cases for finitely many delta potentials. Finally, we investigate the solution of time dependent Schrödinger equation for δ′ potential with particular forms of the strengths.

1. Introduction

Dirac delta potentials in quantum mechanics have been used to model different physical systems almost since the beginning of quantum mechanics. Kronig Penney model [1] is the well-known example of these models. These potentials are the particular cases of much more general class of potentials, namely point interactions. In one dimension, one rigorous way of defining the point interaction at the origin is based on the self-adjoint extension of the free symmetric Hamiltonian defined on ℝ \ {0}. In this approach, the initially ill-defined formal δ and δ′ function potentials appear naturally as two special cases of point interactions constructed from the self-adjoint extension theory. In general one has a 4 -parameter family of self-adjoint extension in one dimension. The monograph [2] includes a great deal of all the details and summarize the history of the literature about the δ interactions. The review article [3] and the book [4] are also good reference sources about the δ potentials from the physical point of view.

The δ′ perturbation of free Hamiltonian H0=-d2dx2 is defined as a limit of short range potentials in the distributional sense [57]. Although there are some controversial issues about δ′ interactions (see e.g., [811]), they are also getting considerable amount of interest. The ambiguities about δ′ interactions have been summarized in a very recent article [11], where the integral form of the Schrödinger equation for δ′ potential has been studied based on the work of Kurasov [12]. We also adopt the distributional approach developed by Kurasov [12] for the functions having a discontinuity at the point of δ′. It is possible to overcome these ambiguities by considering different choices, as different type of δ′ interactions [13]. Therefore, the different results on the spectrum of δ′ potential obtained in [2] and in [8] using the Kurasov's approach -as a special case of −(x) + ′(x) potential- can be interpreted consistently. In other words, the Kurasov's approach corresponds to different self-adjoint extension of the free Hamiltonian H0. These self-adjoint extensions are given by matching conditions at the origin (or at the point supporting the perturbation) and two of these matching conditions maybe identified as a δ′ interaction and receive the names of non-local and local δ′ interaction, respectively [13]. Both approaches are used to see how the spectrum of the V-shaped potential changes when it is perturbed by −(x) + ′(x) perturbation in [9, 13]. The scattering, and resonant tunneling for δ′ potential in one dimension is also a controversial issue because different results are obtained in the literature [8, 10], depending on whether the δ′ interaction is the non-local or the local one. The results obtained in [1416] for the non-local δ′ potential, show that it is opaque for all energies of an incoming beam. However, other authors [1720] claim that there are discrete energy values in the spectrum of δ′ potential which lead to resonant tunneling. As an another physical application, δ′ interactions are used to model Casimir effect in [21, 22].

The exact expression of the propagator for one dimensional single Dirac delta potentials have been found in different ways [2325]. The generalization to two center case have been studied in Cacciari and Moretti [26]. The propagator for general four parameter family of point interactions have been given in Albeverio et al. [27]. Propagators for systems involving δ potentials are also studied from various points of view in references [2834]. The propagator for derivatives of Dirac delta distribution for constant strengths has been recently studied in Lange [11]. The Cauchy problem for the non-local δ′ potential with a time dependent strength has also been studied rigorously in detail [35]. Moreover, time dependent one dimensional point interactions have been studied in Campbell [36] and the exact solution to the initial value problem for Schrödinger equation has been given for some particular form of strength λ(t) of the Dirac delta potential. Later on, the system has been investigated in Hmidi et al. [37] more rigorously and the regularity assumptions on λ(t) is determined for which the initial value problem defines a unitary strongly continuous dynamical system on L2(ℝ). Such time-dependent point interactions have been studied rigorously in order to model asymptotic complete ionization and suitable conditions on the function λ(t) has been determined for ionization problem [3841]. The higher dimensional generalizations of the problem have been studied in great detail and summarized in the thesis by Correggi [42]. Transmission properties of a monochromatic beam and wave packets by studying the scattering from the time-dependent δ potential are studied in Martinez and Reichl [43] and Kuhn et al. [32]. Utilizing δ potential with a time dependent coefficient in an infinite well, Baek et al. [44] showed that it is possible to split a wave function which may have applications in statistical mechanics and condensed matter physics. However, it may also lead to philosophical problems [44]. The time dependence of the Dirac delta potentials could also be expressed through the motion of its support [45]. As a physical application, a moving Dirac δ potential is used to describe particle displacement using a standard tunneling microscope [46].

The paper is organized as follows: In section 2, we obtain a formal expression of the propagator for a finite number of Dirac δ potentials with time dependent strengths and solve the time dependent Schrödinger equation for this system. In the following subsections, we investigate one δ potential with time dependent strength in more detail and find the propagator for N Dirac δ potential in the limit that centers are infinitely separated. In section 3, we get an expression of the propagator for a finite number of δ′ potentials with time dependent strengths and solve the time dependent Schrödinger equation for this potential. Finally, we elaborate on one δ′ case.

2. The Propagator for N Dirac Delta Potentials With Time Dependent Strengths

We begin with a one dimensional model in which a free Hamiltonian of the type H0=-22md2dx2 is perturbed with a time dependent potential

V(x,t)=j=1Nλj(t)δ(x-xj).    (1)

The initial value problem of the time-dependent Schrödinger equation for this potential is

itψ(x,t)=[-d2dx2+j=1Nλj(t)δ(x-xj)]ψ(x,t),    (2)

with the given sufficiently smooth function ψ(x, 0). Here we have used the units such that ℏ = 2m = 1 for simplicity. It is well-known that the Laplace transform is a very useful tool to solve initial value problems so we first take the Laplace transformation of Equation (2) with respect to time variable t and get

ψ̄xx(x,s)-iψ(x,0)+isψ̄(x,s)=j=1Nδ(x-xj)L{λj(t)ψ(xj,t)},    (3)

where ψ̄(x,s)=L{ψ(x,t)} and ψ̄xx(x,s)=L{2x2ψ(x,t)}.

After this, we take the Fourier transformation of both sides of Equation (3) with respect to the coordinate variable x and get:

-k2ψ̄^(k,s)-iψ^(k,0)+isψ̄^(k,s)=j=1Ne-ikxjL{λj(t)ψ(xj,t)},    (4)

where ψ̄^(k,s)=F(ψ̄(x,s)) denotes the Fourier transform of ψ̄(x,s) with respect to the variable x. This equation is an algebraic equation for the unknown wave function ψ̄^(k,s) and the solution is easily obtained as

ψ̄^(k,s)=-1k2-is[iψ^(k,0)+j=1Ne-ikxjL{λj(t)ψ(xj,t)}].    (5)

Now, we immediately find the inverse Fourier transform of the ψ̄^(k,s):

ψ̄(x,s)=-i-dk2πeikxk2-isψ^(k,0)                   -j=1NL{λj(t)ψ(xj,t)}-dk2πeik(x-xj)k2-is.    (6)

Using ψ^(k,0)=-dxe-ikxψ(x,0), the Equation (6) becomes

ψ̄(x,s)=-i-dxψ(x,0)-dk2πeik(x-x)k2-is             -j=1NL{λj(t)ψ(xj,t)}-dk2πeik(x-xj)k2-is.    (7)

The integrals in Equation (7) are easily taken using residue theorem

-dk2πeikxk2-is=ieiis|x|2is,    (8)

so we obtain

ψ̄(x,s)=12is-dxeiis|x-x|ψ(x,0)             -i2isj=1Neiis|x-xj|L{λj(t)ψ(xj,t)}.    (9)

Now we need to take inverse Laplace transform to obtain the solution of the time dependent equation for the potential λ(t)δ(x). Using the Bromwich contour [47], one can easily find the inversion

L-1{eiis|x|is}=1iπtexp(ix24t),    (10)

and using the convolution theorem we get the formal solution as

ψ(x,t)=-dx14πitexp[i(x-x)24t]ψ(x,0)             -ij=1N0tdtλj(t)ψ(xj,t)4πi(t-t)exp[i(x-xj)24(t-t)].    (11)

Although this is an explicit formal expression for ψ(x, t), it is not completely expressed in terms of the initial condition ψ(x, 0) and includes the unknown factors ψ(xj, t). These can be found by simply inserting x = xj in the formal solution and then solving the resulting coupled Volterra type integral equations

ψ(xj,t)=-dx14πitexp[i(xj-x)24t]ψ(x,0)-i0tdtλj(t)ψ(xj,t)4πi(t-t)-ik=1kjN0tdtλk(t)ψ(xk,t)4πi(t-t)exp[i(xj-xk)24(t-t)],    (12)

where we split the term k = j in the summation over k. Since this is not an easy problem for a generic function λ(t), we will first investigate for particular cases, where λ is constant and λ is inversely proportional to t. All these results we present in the next subsection is previously obtained by Campbell [36] using a slightly different method, where only one integral transformation with the boundary conditions at the position of δ potential was used.

2.1. Single δ Potential With a Time Dependent Strength

As a particular case of (1), we consider a single delta potential with time-dependent strength

V(x,t)=λ(t)δ(x),    (13)

where N = 1 and x1 = 0. We can formally obtain the solution of the time dependent Schrödinger equation for this case using Equation (11)

ψ(x,t)=-dx14πitexp[i(x-x)24t]ψ(x,0)                -i0tdtλ(t)ψ(0,t)4πi(t-t)exp[ix24(t-t)].    (14)

Actually this result can be directly obtained from the Duhamel's formula [48] for time-dependent Schrödinger equation associated with the Hamiltonian H = H0 + V, where H0 is self-adjoint free Hamiltonian and V is bounded (or relatively H0 -bounded with relative bound <1):

e-itH|ψ0=e-itH0|ψ0+                               (-i)0tdte-i(t-t)H0Ve-itH|ψ0    (15)

for every |ψ0〉 = |ψ(t = 0)〉. This shows that even if we formally take V = λ(t′)|0〉〈0| which corresponds to our Dirac delta potential, one immediately sees that the Duhamel's formula is still formally valid for such singular interactions.

Given the initial condition ψ(x, 0) and the function λ(t), the function ψ(0, t) can be determined by solving the following integral equation:

ψ(0,t)=-dx14πitexp[ix24t]ψ(x,0)                  -i0tdtλ(t)ψ(0,t)4πi(t-t).    (16)

However, this is in general hard to solve and one usually applies some approximation techniques, e.g., Dyson series [48]. Nevertheless, as shown in Campbell [36], there are cases where one can calculate the Green's function explicitly. We will show two such cases here explicitly. Although Green's functions for these cases are derived in Campbell [36], we repeat these results here for the sake of completeness. Instead of directly solving the above integral equation for particular cases, it is convenient to start with the Laplace transformed wave function Equation (9) for a single δ potential centered at the origin:

ψ̄(x,s)=12is-dxeiis|x-x|ψ(x,0)             -i2iseiis|x|L{λ(t)ψ(0,t)}.    (17)

The cases for constant λ and λ1t are reviewed in Appendices A and B in detail.

2.2. The Propagator for N Dirac δ Potentials in a Limiting Case

When λj's are constant, the Laplace transformed wave function ψ̄(x,s) given by (9) becomes

ψ̄(x,s)=12is-dxeiis|x-x|ψ(x,0)             -i2isk=1Nλkeiis|x-xk|ψ̄(xk,s).    (18)

The unknown functions ψ̄(xk,s) can be found by evaluating the above expression at x = xj:

ψ̄(xk,s)=j=1N[Φ-1(s)]kjρ̄(xj,s),    (19)

where

Φkj(s)={1+iλj2isif j=kiλk2iseiis|xj-xk|if jk,    (20)

and

ρ̄(xk,s)=12is-dxeiis|xk-x|ψ(x,0).    (21)

Substituting (19) into (18), we obtain

ψ̄(x,s)=ρ̄(x,s)-i2isj=1Nk=1Nλkeiis|x-xk|[Φ-1(s)]kjρ̄(xj,s).    (22)

Although we have obtained the Laplace transformed wave function ψ̄(x,s), it is not explicitly given since one has to invert the matrix Φ and find the inverse Laplace transform of the resulting expression to get the final solution ψ(x, t). In general, it is difficult to find the inverse Laplace transforms so one may apply some approximation schemes [49]. Moreover, one could use some numerical computations, but we will here simply show the limiting case, where the centers are infinitely far away from each other.

When all the centers are infinitely separated from each other, that is, |xjxk| → ∞, we expect that the off-diagonal elements of the matrix Φ given in (20) vanish, so that

[Φ-1(s)]kj=(11+iλk2is)δkj.    (23)

Then, ψ̄(x,s) can be explicitly found as

ψ̄(x,s)=ρ̄(x,s)-ij=1Neiis|x-xj|(λj2is+iλj)ρ̄(xj,s),    (24)

from which the propagator G(x, x′, s) reads

G(x,x,s)=eiis|x-x|2is-j=1Neiis|x-xj|(λj2s+iλj)eiis|xj-x|2s.    (25)

Hence,

G(x,x,t)=12iπtexp[i(x-x)24t]-j=1Nλj4exp[λj2(|x-xj|+|xj-x|)+iλj24t]erfc [(|x-xj|+|xj-x|)2it+itλj2].    (26)

This is actually the superposition of the individual propagators associated with single delta centers. This is expected since there is no correlation among the centers when they are far away from each other. Another limiting case is the case where all the centers coincide.

3. The Propagator for N Dirac δ′ Potentials With Time Dependent Strengths

In this section we first obtain a formal solution of the time-dependent Schrödinger equation, where the potential term is chosen formally as

V(x,t)=i=1Nλi(t)δ(x-xi).    (27)

The time dependent Schrödinger equation for this potential is

itψ(x,t)=[-d2dx2+j=1Nλj(t)δ(x-xj)]ψ(x,t).    (28)

As we mentioned in the introduction we adopt the distributional approach given by [12] for the definition of δ′(x), for functions having discontinuity at the point of δ′

δ(x)f(x)=f(0)δ(x)-fx(0)δ(x),    (29)

where we define f(y)=f(y+)+f(y-)2 and fx(y)=fx(y+)+fx(y-)2 and f(y±) denote the limits limxy±f(x). Note that this definition reduces to the well-known property of δ′(x) [50]

δ(x)f(x)=f(0)δ(x)-f(0)δ(x).    (30)

for continuous functions. As in the previous section we proceed by taking the Laplace transform of all the terms in the Equation (28) for time variable t and find

ψ̄xx(x,s)-iψ(x,0)+isψ̄(x,s)=j=1N[-L{λj(t)ψ(xj)}δ(x-xj)+L{λj(t)ψx(xj)}δ(x-xj)].    (31)

Now, we take the Fourier transform with respect to the variable x and solve ψ̄^(k,s) to get:

ψ̄^(k,s)=-1k2-is(iψ^(k,0)+j=1N-ike-ikxjL{λj(t)ψ(xj)}                                   +e-ikxjL{λj(t)ψ(xj)}).    (32)

Before taking the inverse Fourier transform of this equation we write ψ^(k,0)=-dxe-ikxψ(x,0) then take the inverse Fourier transform and get

ψ̄(x,s)=12is-dxeiis|x-x|ψ(x,0)-12j=1Nsgn(x-xj)eiis|x-xj|L{λj(t)ψ(xj,t)}-i12isj=1Neiis|x-xj|L{λj(t)ψx(xj,t)}.    (33)

where sgn(x) denotes sign function sgn(x)={1x>01x<0. We have also used -dk2πeikxk2-is=-i2iseiis|x| and -dk2πikeikxk2-is=-12iseiis|x| sgn(x). Now using L{exp[ix24t]2iπt}=eiis|x|2is and L{(-1)3/4exp[ix24t]x2πt3/2}=eiis|x| sgn(x) and convolution theorem for Laplace transform we get the formal expression of the wave function

ψ(x,t)=-dx14πitei(x-x)24tψ(x,0)-0tdt{j=1N(-1)3/44πλj(t)ψ(xj)exp[i(x-xj)24(t-t)](t-t)3/2(x-xj)+i2iπj=1Nλj(t)ψx(xj)exp[i(x-xj)24(t-t)](t-t)}.    (34)

Similar to the time dependent δ potential case, this formal equation contains unknown functions ψ(xj±)'s and ψx(xj±)'s. One can find equations for ψ(xj±)'s by inserting x=xj± to the Equation (34). In order to find equations for ψx(xj)±'s, the derivative of the Equation (34) with respect to the variable x has to be calculated. By taking this derivative we get

ψx(x,t)=14iπt3/2-dxi(x-x)ei(x-x)24tψ(x,0)-0tdt{-18iπj=1Nλj(t)ψ(xj)i(x-xj)2exp[i(x-xj)24(t-t)](t-t)5/2+(-1)3/44πj=1Nλj(t)ψ(xj)exp[i(x-xj)24(t-t)](t-t)3/2-14iπj=1Nλj(t)ψx(xj)(x-xj)exp[i(x-xj)24(t-t)](t-t)3/2}.    (35)

Now putting x=xj± in the Equation (35) one can get also integral equations for ψx(xj±)'s. Finally one has to solve the system of integral equations for ψ(xj±)'s and ψx(xj±)'s to get the complete solution.

3.1. Single δ′ Potential With a Time Dependent Strength

Now we will elaborate more on the single δ′ interaction which is described by the potential

V(x,t)=λ(t)δ(x)    (36)

The formal solution of the time dependent Schrödinger equation for single δ′ can be obtained from the above section for N = 1, x1 = 0 and λ1 = λ

ψ(x,t)=-dx14πitei(x-x)24tψ(x,0)-0tdt{(-1)3/44πλ(t)ψ(0)xexp[ix24(t-t)](t-t)3/2+i2iπλ(t)ψx(0)exp[ix24(t-t)](t-t)}.    (37)

Similarly the expression for the ψx(x, t) is obtained from the Equation (35)

ψx(x,t)=14iπt3/2-dxi(x-x)ei(x-x)24tψ(x,0)-0tdt{-18iπλ(t)ψ(0)x2exp[ix24(t-t)](t-t)5/2+(-1)3/44πλ(t)ψ(0)exp[ix24(t-t)](t-t)3/2-14iπλ(t)ψx(0)xexp[ix24(t-t)](t-t)3/2}.    (38)

For a given λ(t) one gets integral equations for ψ(x, 0) and ψx(x, 0) by inserting x = 0 in Equations (37) and (38), respectively:

ψ(0,t)=-dx14πiteix24tψ(x,0)-0tdt{(-1)3/44πλ(t)ψ(0)1(t-t)3/2+i2iπλ(t)ψx(0)1(t-t)}.    (39)
ψx(0,t)=14iπt3/2-dx(-x)eix24tψ(x,0)          -(-1)3/44π0tdtλ(t)ψ(0)1(t-t)3/2.    (40)

The Equations (39) and (40) constitute an equation system for ψ(0, t) and ψx(0, t). Solving this system one can determine ψ(0, t) and ψx(0, t) and insert them to Equation (37) to get the wave function for all times. When studying special cases the expression of the Laplace transform of the wave function is necessary. Therefore, utilizing Equation (33) we write the general formula of the Laplace transform of the wave function for a delta prime at x1 = 0 and λ1(t) = λ(t):

ψ̄(x,s)=12is-dxeiis|x-x|ψ(x,0)-12sgn(x)eiis|x| L{λ(t)ψ(0,t)}-i12iseiis|x| L{λ(t)ψx(0,t)}.    (41)

Now will investigate some special cases.

3.1.1. Case 1: λ Is Constant

When the strength of the single δ′ interaction is constant, we obtain from Equation (33)

ψ̄(x,s)=12is-dxeiis|x-x|ψ(x,0)-λ2sgn(x)eiis|x|ψ(0̄,s)-iλ2iseiis|x|ψx(0̄,s),    (42)

where ψ̄(0,s)=L{ψ(x,t)} and ψ̄x(0,s)=L{ψx(x,t)}. We need also derivative of the Laplace transformed wave function ψ̄(x,s) with respect to x.

ψ̄x(x,s)=i2-dx sgn(x-x)eiis|x-x|ψ(x,0)-iλis2sgn2(x)eiis|x|ψ(0̄,s)-λ2d sgn(x)dxeiis|x|ψ(0̄,s)+λ sgn(x)eiis|x|2ψx(0̄,s),    (43)

where we have used d|x|dx=sgn(x). Now we find wave function and its derivative at 0± by choosing x = 0± in Equations (42) and (43):

ψ̄(0±,s)=12is-dxeiis|x|ψ(x,0)λ2ψ(0̄,s)                    -iλ2isψx(0̄,s),    (44)

where we take sgn(0±) = ±1. In order to find ψ̄(x,s) we need to calculate ψ̄(0,s) and ψ̄x(0,s). From Equation (44) we get

ψ̄(0,s)+iλ2isψ̄x(0,s)=12is-dxeiis|x|ψ(x,0).    (45)

In order to find another equation for ψ̄(0,s) and ψ̄x(0,s) we calculate ψ̄x(0±,s)

ψ̄x(0±,s)=i2-dx sgn(-x)eiis|x|ψ(x,0)-iλis2ψ̄(0,s)λ2ψ̄x(0,s).    (46)

where we take d sgn(x)dx|x=0±=0. Using this equation we get

iλis2ψ̄(0,s)+ψ̄x(0,s)=i2-dx sgn(-x)eiis|x|ψ(x,0).    (47)

Solving Equations (45) and (47) we get

ψ̄(0,s)=12is(1+λ24)-1                         -dx(1-λ2sgn(x)) eiis|x| ψ(x,0),    (48)

and

ψ̄x(0,s)=-dx[-i2sgn(x)-iλ4(1+λ24)-1                              (1-λ2sgn(x))] eiis|x|ψ(x,0),    (49)

respectively. Substituting the expressions in Equations (48) and (49) into Equation (42) we get

ψ̄(x,s)=-dx{eiis|x-x|2is+(1+λ24)-1[-λ4sgn(x)-λ28-λ4sgn(x)+λ28sgn(x) sgn(x)]eiis(|x|+|x|)is}ψ(x,0).    (50)

The factor that multiplies ψ(x′, 0) in the integral of the equation above is the Green's function:

(x,x,s)={eiis|x-x|2is+(1+λ24)-1[-λ4sgn(x)-λ28-λ4sgn(x)+λ28sgn(x) sgn(x)]eiis(|x|+|x|)is}.    (51)

The Green's function in terms of s variable can be converted using Equation (10) to get

G(x,x,t)=14πit{ei(x-x)24t+(1+λ24)-1[-λ2sgn(x)-λ24-λ2sgn(x)+λ24sgn(x) sgn(x)]ei(|x|+|x|)24t}.    (52)

So, we can write the wave function for all times as:

ψ(x,t)=14πit-dx{ei(x-x)24t+(1+λ24)-1[-λ2sgn(x)-λ24-λ2sgn(x)+λ24sgn(x) sgn(x)]ei(|x|+|x|)24t}ψ(x,0).    (53)

This result is completely consistent with the one given in Lange [11] except the convention λ = −c.

3.1.2. Case 2: λ(t) = α/t

The general expression (41) for the Laplace transform of the wave function in this particular case λ(t)=αt, where α is a constant, becomes

ψ̄(x,s)=12is-dxeiis|x-x| ψ(x,0)-12sgn(x)eiis|x|L{αtψ(0,t)}-i12iseiis|x| L{αtψx(0,t)}.    (54)

Using the identity given in Equation (B.2) in Appendix B this equation becomes

ψ̄(x,s)=12is-dx eiis|x-x| ψ(x,0)-12sgn(x)eiis|x|αsdsψ̄(0,s)-i12iseiis|x|αsdsψ̄x(0,s).    (55)

In order to calculate this expression we need ψ(0±, s) and ψx(0±,s). The first of this is easily calculated by choosing x = 0± in Equation (55)

ψ̄(0±,s)=12is-dx eiis|x| ψ(x,0)α2sdsψ̄(0,s)-iα2issdsψ̄x(0,s).    (56)

In order to find ψx(0±,s) we take the derivative of the Equation (54) with respect to x variable and get

ψ̄x(x,s)=i2-dx sgn(x-x)eiis|x-x| ψ(x,0)-iisα2sgn2(x)eiis|x|sdsψ̄(0,s)-2eiis|x|αd sgn(x)dxsdsψ̄(0,s)+α2sgn(x)eiis|x|sdsψ̄x(0,s).    (57)

From this equation we get

ψ̄x(0±,s)=-i2-dx sgn(x)eiis|x| ψ(x,0)-iisα2sdsψ̄(0,s)±α2sdsψ̄x(0,s).    (58)

We denote u1(s)=sdsψ̄(0,s) and u2(s)=sdsψ̄x(0,s). We obtain from Equations (56) and (58)

ψ̄(0,s)=I0(s)2is-iα2isu1(s)    (59)

and

ψ̄x(0,s)=-iI1(s)2-iisα2u2(s).    (60)

Here I0(s)=-dxeiis|x|ψ(x,0) and I1(s)=-dx sgn xeiis|x|ψ(x,0). Note that ψ̄(0,s)=-du1(s)ds and ψ̄x(0,s)=-du2(s)ds. Inserting these equalities into Equations (59) and (60) we obtain two coupled differential equations:

du1(s)ds-iα2isu2(s)=-I0(s)2is    (61)

and

du2(s)ds-iisα2u1(s)=iI1(s)2is.    (62)

Although the coefficients of the unknown functions are not constants in these coupled equations, by taking the derivative of the Equations (61) and (62) with respect to the variable s one can uncouple these equations and get second order differential equations for u1(s) and u2(s):

d2u1(s)ds2+12sdu1(s)ds+α24u1(s)=α4isI1(s)-12isdI0(s)ds    (63)

and

d2u2(s)ds2-12sdu2(s)ds+α24u2(s)=14isI1(s)-iα4I0(s)+i2dI1(s)ds.    (64)

The solutions of these equations are elementary and easily obtained after some algebra as

u1(s)=-π42-dxψ(x,0) {0sdss1/4s3/4isg1(s,s)                  [α sgn(x)+|x|is]eiis|x|}    (65)

where

g1(s,s)=J14(αs2)J-14(αs2)-J14(αs2)J-14(αs2).    (66)

Here and in the following expressions Ja(x), stands for the first kind of Bessel's function. The solution of the differential equation in Equation (64) is

u2(s)=-π42-dxψ(x,0) {0sdss3/4s1/4g2(s,s)                  [sgn(x)is-iα-ixis]eiis|x|}    (67)

where

g2(s,s)=J34(αs2)J-34(αs2)-J34(αs2)J-34(αs2).    (68)

Note that

L{1tψ(0,t)}=u1(s)=sdsψ̄(0,s)    (69)

and

L{1tψx(0,t)}=u2(s)=sdsψ̄x(0,s).    (70)

Thus using Equations (69) and (70) in (41) we get

ψ̄(x,s)=-dxψ(x,0){eiis|x-x|2iseiis|x-x|+sgn(x)π α s1/4eiis|x|82isdss1/4g1(s,s)[α sgn(x)+|x|is]eiis|x|+iπ α s1/4eiis|x|82sdss1/4g2(s,s)[sgn(x)is-iα-ixis]eiis|x|}.    (71)

Thus the Green's function in Laplace transformed space for the Schrödinger equation with a potential V(x,t)=αtδ(x) is

G(x,x,s)={eiis|x-x|2iseiis|x-x|+sgn(x)π α s1/4eiis|x|82isdss1/4g1(s,s)[α sgn(x)+|x|is]eiis|x|+iπ α s1/4eiis|x|82sdss1/4g2(s,s)[sgn(x)is-iα-ixis]eiis|x|}.    (72)

This Green' s function cannot be converted in terms of elementary functions but it is possible to use numerical methods to obtain Green' s function in the position-time space.

4. Conclusion

In this work, we have studied some analytically solvable time-dependent point interactions. First, we have obtained a formal expression of the propagator for finite number of Dirac δ potentials with time dependent strengths and solved the time dependent Schrödinger equation for this system. Then we have investigated one δ potential with various time dependent strengths in more detail and found the propagator for N Dirac δ potential in the limit that centers are infinitely separated. Furthermore, we have found an expression of the propagator for finite number of δ′ potentials with time dependent strengths and solved the time dependent Schrödinger equation for this potential. We believe that these results obtained are useful in models of ionization problems, where the particle is initially bound to the time dependent δ or δ′ potentials. Such type of models have been studied (see e.g., [42, 51]) and the results obtained are compared with the experiment [52].

Data Availability Statement

All datasets analyzed for this study are included in the article/Supplementary Material.

Author Contributions

All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.

Funding

The authors thank TUBITAK for the 2221-Visiting Scientist Fellowship Programme (2016). MG acknowledged partial financial support from the Spanish MINECO (Grant no. MTM2014-57129-C2-1-P) and the Junta de Castilla y León (Grant nos. VA137G18 and BU229P18).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphy.2020.00065/full#supplementary-material

References

1. Kronig R de L, Penney WG. A quantum mechanics of electrons in crystal lattices. Proc R Soc. (1931) 130:499.

Google Scholar

2. Albeverio S, Gesztesy F, Hoegh-Krohn R, Holden H. Solvable Models in Quantum Mechanics. 2nd ed. Chelsea, RI: American Mathematical Society (2004).

Google Scholar

3. Belloni M, Robinett RW. The infinite well and Dirac delta function potentials as pedagogical, mathematical and physical models in quantum mechanics. Phys Rep. (2014) 540:25. doi: 10.1016/j.physrep.2014.02.005

CrossRef Full Text | Google Scholar

4. Demkov Yu N, Ostrovskii VN. Zero-range Potentials and Their Applications in Atomic Physics. New York, NY: Plenum Press (1988).

Google Scholar

5. Golovaty, YD, Hryniv, RO. Norm resolvent convergence of singularly scaled Schrödinger operators and δ′-potentials. Proc R Soc Edinb A Math. (2013) 143:791–816. doi: 10.1017/S0308210512000194

CrossRef Full Text | Google Scholar

6. Golovaty Y. Two-parametric δ′-interactions: approximation by Schrödinger operators with localized rank-two perturbations. J Phys A Math Theor. (2018) 51:255202. doi: 10.1088/1751-8121/aac110

CrossRef Full Text | Google Scholar

7. Golovaty Y. Schrödinger operators with singular rank-two perturbations and point interactions. Integr Equat Oper Theory. (2018) 90:57. doi: 10.1007/s0002

CrossRef Full Text | Google Scholar

8. Gadella M, Negro J, Nieto LM. Bound states and scattering coefficients of the −(x) + ′(x) potential. Phys Lett A. (2009) 373:1310–3. doi: 10.1016/j.physleta.2009.02.025

CrossRef Full Text | Google Scholar

9. Gadella M, Glasser ML, Nieto LM. One dimensional models with a singular potential of the type −αδ(x) + βδ′(x). Int J Theor Phys. (2010) 50:2144–52. doi: 10.1007/s10773-010-0641-6

CrossRef Full Text | Google Scholar

10. Zolataryuk AV. Boundary conditions for the states with resonant tunnelling across the δ′ potential. Phys Lett A. (2010) 374:1636–41. doi: 10.1016/j.physleta.2010.02.005

CrossRef Full Text | Google Scholar

11. Lange RJ. Distribution theory for Schrödinger's integral equation. J Math Phys. (2015) 56:122105. doi: 10.1063/1.4936302

CrossRef Full Text | Google Scholar

12. Kurasov P. Distribution theory for discontinuous test functions and differential operators with generalized coefficients. J Math Anal Appl. (1996) 201:297–323. doi: 10.1006/jmaa.1996.0256

CrossRef Full Text | Google Scholar

13. Fassari S, Gadella M, Nieto LM, Rinaldi F. Spectroscopy of a one-dimensional V-shaped quantum well with a point impurity. Ann Phys. (2018) 389:48–62. doi: 10.1016/j.aop.2017.12.006

CrossRef Full Text | Google Scholar

14. Seba P. Some remarks on the δ′-interaction in one dimension. Rep Math Phys. (1986) 24:111. doi: 10.1016/0034-4877(86)90045-5

CrossRef Full Text | Google Scholar

15. Patil SH. Schrödinger equation with δ′ and δ″ potentials. Phys Scr. (1994) 49:645. doi: 10.1088/0031-8949/49/6/002

CrossRef Full Text | Google Scholar

16. Coutinho FAB, Nogami Y, Perez JF. Generalized point interactions in one-dimensional quantum mechanics. J Phys A Math Gen. (1997) 30:3937. doi: 10.1088/0305-4470/30/11/021

CrossRef Full Text | Google Scholar

17. Christiansen PL, Arnbak NC, Zolotaryuk AV, Ermakov VN, Gaididei YB. On the existence of resonances in the transmission probability for interactions arising from derivatives of Dirac's delta function. J Phys A Math Gen. (2003) 36:7589. doi: 10.1088/0305-4470/36/27/311

CrossRef Full Text | Google Scholar

18. Zolotaryuk AV, Christiansen PL, Iermakova SV. Scattering properties of point dipole interactions. J Phys A Math Gen. (2006) 36:9329. doi: 10.1088/0305-4470/39/29/023

CrossRef Full Text | Google Scholar

19. Toyama FN, Nogami Y. Transmission-reflection problem with a potential of the form of the derivative of the delta function. J Phys A Math Theor. (2007) 40:F685. doi: 10.1088/1751-8113/40/29/F05

CrossRef Full Text | Google Scholar

20. Zolotaryuk AV. Two-parametric resonant tunneling across the δ′(x) potential. Adv Sci Lett. (2008) 1:187. doi: 10.1166/asl.2008.019

CrossRef Full Text | Google Scholar

21. Muñoz-Castañeda JM and Mateos-Guilarte J. Delta-delta' generalized Robin boundary conditions and quantum vacuum fluctuations. Phys Rev D. (2015) 91:025028. doi: 10.1103/PhysRevD.91.025028

CrossRef Full Text | Google Scholar

22. Silva, JDL, Braga, AN, Alves, DT. Dynamical Casimir effect with delta-delta' mirrors. Phys Rev D. (2016) 94:105009. doi: 10.1103/PhysRevD.94.105009

CrossRef Full Text | Google Scholar

23. Gaveau B, Schulman LS. Explicit time-dependent Schrödinger propagators. J Phys A Math Gen. (1986) 19:1833–46. doi: 10.1088/0305-4470/19/10/024

CrossRef Full Text | Google Scholar

24. Manoukian EB. Explicit derivation of the propagator for a Dirac delta potential. J Phys A Math Gen. (1989) 22:67. doi: 10.1088/0305-4470/22/1/013

CrossRef Full Text | Google Scholar

25. Blinder SM. Greens function and propagator for the one-dimensional delta function potential. Phys Rev A. (1988) 37:973. doi: 10.1103/PhysRevA.37.973

CrossRef Full Text | Google Scholar

26. Cacciari I, Moretti P. Propagator for the double delta potential. Phys Lett A. (2006) 359:396–401. doi: 10.1016/j.physleta.2006.06.061

CrossRef Full Text | Google Scholar

27. Albeverio S. Brzezniak Z, Dabrowski L. Time-dependent propagator with point interaction. J Phys A Math Gen. (1994) 27:4933–43.

Google Scholar

28. Elberfeld W, Kleber M. Time-dependent tunneling through thin barriers: a simple analytical solution. Am J Phys. (1988) 56:154. doi: 10.1119/1.15695

CrossRef Full Text | Google Scholar

29. Crandall RE. Combinatorial approach to Feynman path integration. J Phys A Math Gen. (1993) 26:3627. doi: 10.1088/0305-4470/26/14/024

CrossRef Full Text | Google Scholar

30. Cheng BK, da Luz MGE. Propagator for the δ-function potential moving with constant velocity. Phys Rev A. (1993) 47:4720–4. doi: 10.1103/PhysRevA.47.4720

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Scheitler G, Kleber M. Propagator for two dispersing δ-function potentials. Phys Rev A. (1990) 42:55–60. doi: 10.1103/PhysRevA.42.55

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Kuhn J, Zanetti FM, Azevedo AL, Schmidt AGM, Cheng BK, da Luz MGE. Time-dependent point interactions and infinite walls: some results for wavepacket scattering. J Opt B Quant Semiclass Opt. (2005) 7:S77–85. doi: 10.1088/1464-4266/7/3/011

CrossRef Full Text | Google Scholar

33. Laissaoui A, Chetouani L. Propagator for symmetric double delta potential using path decomposition method. J Math Phys. (2018) 59:022106. doi: 10.1063/1.5008516

CrossRef Full Text | Google Scholar

34. Refaei A, Kheirandish F. Quantum propagator and characteristic equation in the presence of a chain of δ potentials. Int J Mod Phys B. (2015) 29:1550099. doi: 10.1142/S021797921550099X

CrossRef Full Text | Google Scholar

35. Cacciapuoti C, Mantile A, Posilicano A. Time dependent delta-prime interactions in dimension one. Nanosyst Phys Chem Math. (2016) 2:303–14. doi: 10.17586/2220-8054-2016-7-2-303-314

CrossRef Full Text | Google Scholar

36. Campbell J. Some exact results for the Schrödinger wave equation with a time dependent potential. J Phys A Math Theor. (2009) 42:365212. doi: 10.1088/1751-8113/42/36/365212

CrossRef Full Text | Google Scholar

37. Hmidi T, Mantile A, Nier F. Time-dependent delta-interactions for 1D Schrödinger Hamiltonians. Math Phys Anal Geom. (2010) 13:83–103. doi: 10.1007/s11040-009-9070-2

CrossRef Full Text | Google Scholar

38. Costin O, Costin RD, Lebowitz JL, Rokhlenko A. Evolution of a model quantum system under time periodic forcing: conditions for complete ionization. Comm Math Phys. (2001) 221:1–26. doi: 10.1007/s002200100

CrossRef Full Text | Google Scholar

39. Costin O, Lebowitz JL, Rokhlenko A. Decay versus survival of a localized state subjected to harmonic forcing: exact results. J Phys A Math Gen. (2002) 35:8943–51. doi: 10.1088/0305-4470/35/42/305

CrossRef Full Text | Google Scholar

40. Costin O, Lebowitz JL, Rokhlenko A. Exact results for the ionization of a model quantum system. J Phys A Math Gen. (2000) 33:6311–9. doi: 10.1088/0305-4470/33/36/303

CrossRef Full Text | Google Scholar

41. Costin O, Costin RD, Lebowitz JL. Transition to the continuum of a particle in time-periodic potentials. In: Karpeshina Y, Stolz G, Weikard R, Zeng Y, editors. Advances in Differential Equations and Mathematical Physics, AMS Contemporary Mathematics Series. Providence, RI: AMS. (2003). p. 75–86.

Google Scholar

42. Correggi M. Time-dependent singular interactions (Ph.D. thesis), Mathematical Physics Sector, International School for Advanced Studies (2004).

Google Scholar

43. Martinez DF, Reichl LE. Transmission properties of the oscillating delta-function potential. Phys Rev B. (2001) 64:245315. doi: 10.1103/PhysRevB.64.245315

CrossRef Full Text | Google Scholar

44. Baek SK, Yi SD, Kim M. Particle in a box with a time dependent delta-function potential. Phys Rev A. Upper Saddle River: Prentice Hall (2016) 96:052124. doi: 10.1103/PhysRevA.94.052124

CrossRef Full Text | Google Scholar

45. Griffiths DJ. Introduction to Quantum Mechanics. 2nd ed. Printice Hall (2005).

46. Granot E, Marchewka A. Quantum particle displacement by a moving localized potential trap. Europhys Lett. (2009) textbf47:20007. doi: 10.1209/0295-5075/86/20007

CrossRef Full Text | Google Scholar

47. Appel W. Mathematics for Physics and Physicists. Princeton, NJ: Princeton University Press (2007).

Google Scholar

48. Reed M, Simon B. Methods of Modern Mathematical Physics. II. New York, NY: Academic Press (1975).

Google Scholar

49. Akin JE, Counts J. On rational approximation to the inverse laplace transform. SIAM J Appl Math. (1969) 17:1035–40.

Google Scholar

50. Hoskins RF. Delta Functions, Introduction to Generalized Functions. 2nd ed. Cambridge: Woodhead Publishing (2011).

Google Scholar

51. Qiu Z. Study of ionization of quantum systems with delta potentials in damped and undamped time periodic fields (Ph.D. thesis), The Ohio State University, Columbus, OH, United States (2009).

Google Scholar

52. Koch PM, Van Leeuwen KAH. The importance of resonances in microwave “ionization” of excited hydrogen atoms. Phys Rep. (1995) 255:289–403.

Google Scholar

Keywords: propagator, delta potentials, delta prime potentials, Green's function, time dependent Schrödinger equation

Citation: Erman F, Gadella M and Uncu H (2020) The Propagators for δ and δ′ Potentials With Time-Dependent Strengths. Front. Phys. 8:65. doi: 10.3389/fphy.2020.00065

Received: 28 February 2019; Accepted: 28 February 2020;
Published: 15 April 2020.

Edited by:

Horacio Sergio Wio, Institute of Interdisciplinary Physics and Complex Systems (IFISC), Spain

Reviewed by:

Abdullahi Yusuf, Federal University, Dutse, Nigeria
Yaroslav Zolotaryuk, Bogolyubov Institute for Theoretical Physics, Ukraine

Copyright © 2020 Erman, Gadella and Uncu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Haydar Uncu, huncu@adu.edu.tr

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.