Acta Polytechnica doi:10.14311/AP.2017.57.0385 Acta Polytechnica 57(6):385–390, 2017 © Czech Technical University in Prague, 2017 available online at http://ojs.cvut.cz/ojs/index.php/ap ON THE SPECTRUM OF THE ONE-DIMENSIONAL SCHRÖDINGER HAMILTONIAN PERTURBED BY AN ATTRACTIVE GAUSSIAN POTENTIAL Silvestro Fassaria, b, c, ∗, Manuel Gadellaa, Luis Miguel Nietoa, Fabio Rinaldib, c a Departamento de Física Teórica, Atómica y Óptica, and IMUVA, Universidad de Valladolid, 47011 Valladolid, Spain b CERFIM, PO Box 1132, Via F. Rusca 1, CH-6601 Locarno, Switzerland c Dipartimento di Fisica Nucleare, Subnucleare e delle Radiazioni, Universitá degli Studi Guglielmo Marconi, Via Plinio 44, I-00193 Rome, Italy ∗ corresponding author: silvestro.fassari@uva.es Abstract. We propose a new approach to the problem of finding the eigenvalues (energy levels) in the discrete spectrum of the one-dimensional Hamiltonian H = −∂2x −λe−x 2/2, by using essentially the well-known Birman-Schwinger technique. However, in place of the Birman-Schwinger integral operator we consider an isospectral operator in momentum space, taking advantage of the unique feature of this Gaussian potential, that is to say its invariance under Fourier transform. Given that such integral operators are trace class, it is possible to determine the energy levels in the discrete spectrum of the Hamiltonian as functions of λ with great accuracy by solving a finite number of transcendental equations. We also address the important issue of the coupling constant thresholds of the Hamiltonian, that is to say the critical values of λ for which we have the emergence of an additional bound state out of the absolutely continuous spectrum. Keywords: Schrödinger equation; Gaussian potential; Birman-Schwinger method; trace class operators; Fredholm determinant. 1. Introduction Given its increasing relevance in the field of Nanophysics, it is particularly interesting to investigate the Schrödinger Hamiltonian with an attractive Gaussian potential V = −λe−x 2/2, since the latter has the typical properties of short-range potentials, which implies the existence of bound states with negative energies below the absolutely continuous spectrum given by the semibounded interval [0, +∞), but also those of the harmonic oscillator near the bottom of the well. Hence, the original Schrödinger eigenvalue problem we are going to consider is ( − d2 dx2 −λe−x 2/2 ) ψ = −�2ψ, � > 0 (1) Although there have been quite a few papers on this model, the most rigorous of which being [1] by F. M. Fernández (see also [2, 3]), it seems that the implications of a remarkable feature of this potential, namely its invariance with respect to the Fourier transform, have been missed. Our method, whilst being essentially based on the (by now) classical Birman-Schwinger method (see, e.g., [4, 5]), considers instead of the integral operator associated to (1) λe−x 2/4 ( − d2 dx2 + �2 )−1 e−x 2/4, (2) another integral operator which is isospectral to the Birman-Schwinger one given in (2) (see [6–8]): λB� = λ ( − d2 dx2 + �2 )−1/2 e−x 2/2 ( − d2 dx2 + �2 )−1/2 . (3) Given that the function e−x 2/2 is invariant under the Fourier transform, the unitary equivalent of the above integral operator, still denoted by B�, is B� = 1 √ 2π ( p2 + �2 )−1/2 e−p 2/2 ∗ ( p2 + �2 )−1/2 , (4) 385 http://dx.doi.org/10.14311/AP.2017.57.0385 http://ojs.cvut.cz/ojs/index.php/ap S. Fassari, M. Gadella, L. M. Nieto, F. Rinaldi Acta Polytechnica where the star denotes the convolution. Therefore, the Schrödinger equation (1) can be reformulated in terms of the following integral equation: χ(p) = 1 √ 2π ( p2 + �2 )−1/2 ∫ +∞ −∞ e−(p−q) 2/2(q2 + �2)−1/2χ(q) dq, (5) with ψ̂(p) = (p2 + �2)−1/2χ(p), in the sense that if � is a value for which the above integral operator has an eigenvalue equal to one, then −�2 is an eigenvalue of the original Schrödinger equation. As a consequence, understanding in depth the properties of the integral operator is crucial in order to get a detailed description of the discrete spectrum of the original Hamiltonian. 2. The integral operator B� We wish to open this section by stating and proving the key property of the integral operator B�. Theorem 2.1. The integral operator B� is trace class. Proof. As a consequence of a well-known theorem (see [9]), given the positivity of the integral kernel and its smoothness, the trace is simply given by the integral of the kernel evaluated along the diagonal p = q, that is to say: 1 √ 2π ∫ +∞ −∞ dp p2 + �2 = √ π 2 1 � . (6) The fact that the trace diverges as � → 0+, guarantees that there is always at least one bound state even for very small values of the coupling constant (shallow wells). We take this opportunity to remind the reader that the latter property is typical of one-dimensional quantum Hamiltonians, as shown in [5, 10]. We also notice that, by expanding the square in the exponent of the Gaussian in (5), the integral kernel B�(p,q) of the operator can be recast as: B�(p,q) = 1 √ 2π e−p 2/2 (p2 + �2)1/2 epq e−q 2/2 (q2 + �2)1/2 , (7) As an immediate consequence of (7), the operator B� can be written as a direct sum of two operators, one acting onto the symmetric subspace (s) the other onto the antisymmetric one (a), whose kernels are: Bs� (p,q) = 1 √ 2π e−p 2/2 (p2 + �2)1/2 cosh(pq) e−q 2/2 (q2 + �2)1/2 , (8) Ba� (p,q) = 1 √ 2π e−p 2/2 (p2 + �2)1/2 sinh(pq) e−q 2/2 (q2 + �2)1/2 . (9) By using the Taylor expansion for each hyperbolic function, the two integral operators can be written as infinite sums of rank one operators, namely: Bs� = 1 √ 2π ∞∑ n=0 1 (2n)! ∣∣∣∣ p2ne−p 2/2 (p2 + �2)1/2 〉〈 p2ne−p 2/2 (p2 + �2)1/2 ∣∣∣∣, (10) Ba� = 1 √ 2π ∞∑ n=0 1 (2n + 1)! ∣∣∣∣p2n+1e−p 2/2 (p2 + �2)1/2 〉〈 p2n+1e−p 2/2 (p2 + �2)1/2 ∣∣∣∣. (11) Obviously, as a consequence of (6), both operators are trace class, which implies that both series are convergent in the trace class norm. Furthermore, although either operator has not yet been diagonalised at this stage (because the rank one operators are not mutually orthogonal), their diagonalisation can be achieved starting from the first rank one operator in each expansion (n = 0), that is to say: Bs,0� = 1 √ 2π ∣∣∣∣ e−p 2/2 (p2 + �2)1/2 〉〈 e−p 2/2 (p2 + �2)1/2 ∣∣∣∣, (12) Ba,0� = 1 √ 2π ∣∣∣∣ pe−p 2/2 (p2 + �2)1/2 〉〈 pe−p 2/2 (p2 + �2)1/2 ∣∣∣∣. (13) Their norms can be easily computed as the required improper integrals are well known: ‖Bs,0� ‖ = 1 √ 2π ∫ +∞ −∞ e−p 2 p2 + �2 dp = √ π 2�2 e� 2 erfc(�), (14) ‖Ba,0� ‖ = 1 √ 2π ∫ +∞ −∞ p2e−p 2 p2 + �2 dp = 1 √ 2 ( 1 − √ π�e� 2 erfc(�) ) . (15) 386 vol. 57 no. 6/2017 On the Spectrum of the One-Dimensional Schrödinger Hamiltonian As can be understood from (14), the divergence of the trace of the entire operator as � → 0+ is only due to the the divergence of Bs,0� . In fact, the trace class norm of the positive operator Bs,1� = Bs� −Bs,0� can be calculated as follows: ‖Bs,1� ‖1 = 1 √ 2π ∫ +∞ −∞ e−p 2 (cosh p2 − 1) p2 + �2 dp. (16) By writing the hyperbolic cosine as a combination of two exponentials, the right hand side of (16) is easily seen to be equal to ‖Bs,1� ‖1 = √ π 23�2 ( 1 + e2� 2 erfc(21/2�) − 2e� 2 erfc(�) ) , (17) which, given its removable singularity at the origin, is almost immediately seen to converge to √ 2−1 as � → 0+. Before determining the equation that will enable us to compute the ground state energy as a function of the coupling parameter λ, let us consider also the antisymmetric component Ba� . As a result of the explicit expression of the integral kernel of Ba� , it is not difficult to show that the latter operator converges weakly to Ba0 , where Ba0 (p,q) = 1 √ 2π e−p 2/2 |p| sinh(pq) e−q 2/2 |q| . (18) The latter is a positive trace class operator with trace equal to ‖Ba0‖1 = 1 √ 2π ∫ +∞ −∞ e−p 2 p2 sinh p2 dp = 1, (19) as follows easily by writing the hyperbolic sine as a combination of two exponentials and using integration by parts. It is worth pointing out that the latter quantity is nothing else but the limit, as � → 0+ of the trace class norm of Ba� , since ‖Ba�‖1 = 1 √ 2π ∫ +∞ −∞ e−p 2 sinh p2 p2 + �2 dp = √ π 23�2 ( 1 −e2� 2 erfc( √ 2�) ) , (20) which ensures the convergence in the norm topology of trace class operators, as a consequence of a well-known theorem on operators belonging not only to the trace class T1 but to any ideal Tp (see [11]). This result can be summarised in the following statement. Theorem 2.2. The positive trace class operator Ba� converges, as � → 0+, to Ba0 , the positive trace class operator defined by its kernel in (18), in the norm topology of trace class operators. It is crucial to point out that, whilst Bs,0� , the first rank one summand in the expansion of the symmetric component of our integral operator, diverges as � → 0+, Ba,0� obviously converges to the following rank one operator: Ba00 = 1 √ 2π ∣∣∣∣pe−p 2/2 |p| 〉〈 pe−p 2/2 |p| ∣∣∣∣ = 1√2π∣∣sgn(p)e−p2/2〉〈sgn(p)e−p2/2∣∣. (21) Although the antisymmetric function sgn(p)e−p 2/2 has a jump discontinuity at the origin, its square coincides with the one of the unnormalised ground state eigenfunction of the harmonic oscillator in momentum space. 3. The two lowest eigenvalues of −∂2x −λe−x 2/2 3.1. The ground state As pointed out earlier, the divergent behaviour of the first term (n = 0) of the expansion of the symmetric part (12) implies that, no matter how shallow the Gaussian well may be (small values of the coupling constant λ), there will always exist at least one bound state, the ground state, whose energy is −�0(λ)2 with �0(λ) given by the solution of a transcendental equation that can be derived from the application of well-known facts regarding the Fredholm determinant of a trace class operator (see, e.g., [12]), namely: det (1 −λBs� ) = 0. (22) By isolating the first divergent rank one operator in the expansion of λBs� and taking advantage of the boundedness of λBs,1� , the left hand side of (22) can be rewritten as det ( 1 −λBs,0� (1 −λB s,1 � ) −1) det(1 −λBs,1� ). (23) 387 S. Fassari, M. Gadella, L. M. Nieto, F. Rinaldi Acta Polytechnica 0 1 2 3 4 -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 λ E0 Figure 1. The ground state energy E0(λ) = −�0(λ)2, as a function of the coupling parameter λ. As the first factor is the only one involving the divergent rank one operator Bs,0� , the equation leading to the determination of the ground state energy reduces to det ( 1 −λBs,0� (1 −λB s,1 � ) −1) = 1 − tr(λBs,0� (1 −λBs,1� )−1) = 0, (24) given that the second term inside the deteminant is a rank one operator. By taking only the terms up to λ in the expansion of the inverse inside the trace in (24), we get the following quadratic equation in λ: g0(�) 2π λ2 + f0(�)√ 2π λ− 1 = 0, (25) with f0(�) = ∫ +∞ −∞ e−p 2 p2 + �2 dp = π � e� 2 erfc(�), (26) g0(�) = ∫ +∞ −∞ dp ∫ +∞ −∞ e−p 2 p2 + �2 ( cosh(pq) − 1 ) e−q2 q2 + �2 dq. (27) By using the standard Taylor expansion of the hyperbolic cosine, the latter double integral can be recast as the following convergent series whose coefficients are expressed in terms of the Gamma and the incomplete Gamma functions: g0(�) = ∞∑ n=1 1 (2n)! (∫ +∞ −∞ p2ne−p 2 p2 + �2 dp )2 = e2� 2 ∞∑ n=1 �2(2n−1) (Γ(n + 1/2)Γ(−n + 1/2,�2))2 Γ(2n + 1) . (28) Hence, the positive solution of (25) is given by the function from [0, +∞) to [0, +∞): λ(�) = 2 √ 2π (f20 (�) + 4g0(�))1/2 + f0(�) , (29) which can be inverted to get the required �0(λ), leading to the ground state energy E0(λ) = −�0(λ)2, the plot of which is shown in Figure 1. 3.2. The first excited state Let us now consider the antisymmetric component given by (9) or (11), so that we have to study the equation det (1 −λBa� ) = 0. As a result of the explicit expression of the integral kernel of Ba� , it is not difficult to show that λBa� converges weakly to λBa0 , where the trace class operator Ba0 is given by (18). As a consequence of Theorem 2.2, and in analogy to what has been done for the determination of the ground state energy, E1(λ) = −�1(λ)2, the energy of the first antisymmetric bound state, can be determined by means of the following equation: det ( 1 −λBa,0� (I −λB a,1 � ) −1) = 0, (30) with Ba,1� = Ba� − Ba,0� , which is obviously trace class. Taking account of the fact that Ba,0� is a rank one operator, the latter equation becomes 1 −λ tr ( Ba,0� (I −λB a,1 � ) −1) = 0 (31) 388 vol. 57 no. 6/2017 On the Spectrum of the One-Dimensional Schrödinger Hamiltonian 1.5 2.0 2.5 3.0 3.5 4.0 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0.0 λ E1 Figure 2. The energy of the first excited state E1(λ) = −�1(λ)2, as a function of the coupling parameter λ. By taking only the terms up to λ in the expansion of the inverse, we get the equation: g1(�) 2π λ2 + f1(�)√ 2π λ− 1 = 0, (32) with f1(�) = ∫ +∞ −∞ p2e−p 2 p2 + �2 dp = √ π −π�e� 2 erfc(�), (33) g1(�) = ∫ +∞ −∞ dp pe−p 2 p2 + �2 ∫ +∞ −∞ dq ( sinh(pq) −pq ) qe−q2 q2 + �2 . (34) By using now the standard Taylor expansion of the hyperbolic sine, the latter double integral can be recast as the following convergent series whose coefficients are expressed in terms of the Gamma and the incomplete Gamma function: g1(�) = ∞∑ n=1 1 (2n + 1)! (∫ ∞ −∞ p2n+2e−p 2 p2 + �2 dp )2 = e2� 2 ∞∑ n=1 �2(2n+1) (Γ(n + 3/2)Γ(−n− 1/2,�2))2 Γ(2n + 2) . (35) Hence, the positive solution of (32) is given by λ1(�) = 2 √ 2π (f21 (�) + 4g1(�))1/2 + f1(�) , (36) with domain [0, +∞) and codomain given by [λ1(0), +∞), λ1(0) being λ1(0) = 2 √ 6 √ 4π − 9 + √ 3 ≈ 1.35311. (37) Hence, the latter function can be inverted to get �1(λ), as well as E1(λ) = −�1(λ)2 whose plot is given in Figure 2. The plot of both eigenvalues E0 and E1 as functions of the coupling parameter λ is shown in Figure 3. The isolated points are those resulting from Table 1 in the aforementioned paper by Fernández [1]. 4. Conclusions By combining a variation of the renowned Birman-Schwinger principle and the use of Fredholm determinants, given that all the integral operators involved are positive trace class operators, we have been able to determine the two lowest eigenenergies, the one of the ground state and that of the lowest antisymmetric bound state, of the one-dimensional Hamiltonian −∂2x−λe−x 2/2 as functions of the coupling parameter. Whilst the ground state energy emerges out of the absolutely continuous spectrum at λ = λ0(0) = 0, the latter emerges at λ = λ1(0), approximately equal to 1.35311. The method can be further exploited to determine other eigenvalues of the Hamiltonian. For example, by starting with the function ψ̂2(p) = ( p2 − �(1 − �) √ πe� 2 erfc(�) ) e−p2/2 (p2 + �2)1/2 , (38) 389 S. Fassari, M. Gadella, L. M. Nieto, F. Rinaldi Acta Polytechnica 0 1 2 3 4 -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 λ E Figure 3. The curves of both energies, E0 (blue) and E1 (yellow), as functions of the coupling parameter λ, and the points resulting from the table provided by Fernández in [1]. orthogonal to the ground state eigenfunction ψ̂0(p) = e−p 2/2/(p2 + �2)1/2 (by using the Gram-Schmidt procedure), and the Fredholm determinant det(1 −λBs,1� ), one can determine the energy of the first excited symmetric bound state as a function of the coupling parameter emerging out of the absolutely continuous spectrum at the next coupling constant threshold. Work in this direction is in progress. 5. Acknowledgements Fabio Rinaldi would like to thank Prof. M. Znojil for kindly inviting him to present the main points of this article in the final session of the Conference “Analytic and Algebraic Methods in Physics” held in Prague, 11–14 September 2017. Partial financial support is acknowledged to the Spanish Junta de Castilla y León and FEDER (Project VA057U16), and MINECO (Project MTM2014-57129-C2-1-P). S. Fassari also wishes to thank the entire staff at Departamento de Física Teórica, Atómica y Óptica, Universidad de Valladolid, for their warm hospitality throughout his stay. References [1] Fernández F.M.: Quantum Gaussian Wells and Barriers, Am. J. Physics, 79 (7), 2011, 752-754. [2] Muchatibaya G., Fassari S., Rinaldi F., Mushanyu J.: A Note on the Discrete Spectrum of Gaussian Wells (I): The Ground State Energy in One Dimension, Advances in Mathematical Physics, 2016. [3] Nandi S.: The quantum Gaussian Well, Am. J. Physics, 78, 2010, 1341-1345. [4] Klaus M.: Some applications of the Birman-Schwinger principle, Helv. Physics Acta 55, 1982, 413-419. [5] Klaus M.: On the bound state of Schrödinger operators in one dimension, Annals of Physics 108, 1977, 288-300. [6] Klaus M.: A remark about weakly coupled one-dimensional Schrödinger operators Helv. Phys. Acta 52, 1979, 223. [7] Fassari S.: An estimate regarding one-dimensional point interactions Helv. Phys. Acta 68, 1995, 121-125. [8] Fassari S., Rinaldi F.: On the spectrum of the Schrödinger Hamiltonian with a particular configuration of three point interactions, Reports on Mathematical Physics, 64 (3) 2009, 367-393. [9] Reed M, Simon B.: Methods in Modern Mathematical Physics III - Scattering Theory, Academic Press, NY 1979 [10] Simon B.: The bound state of weakly coupled Schrödinger operators in one and two dimensions, Annals of Physics, 97, 1976, 279-288. [11] Simon B.: Trace Ideals and Their Applications, Cambridge University Press, Cambridge 1979 [12] Reed M, Simon B.: Methods in Modern Mathematical Physics IV - Analysis of Operators, Academic Press, NY 1978 390 Acta Polytechnica 57(6):385–390, 2017 1 Introduction 2 The integral operator B 3 The two lowest eigenvalues of -x2 -e-x2/2 3.1 The ground state 3.2 The first excited state 4 Conclusions 5 Acknowledgements References