Skip to main content

Reconstruction of potential function for Sturm-Liouville operator with Coulomb potential

Abstract

In this paper, we are concerned with an inverse problem for the Sturm-Liouville operator with Coulomb potential using a new kind of spectral data that is known as nodal points. We give a reconstruction of q as a limit of a sequence of functions whose n th term is dependent only on eigenvalue and its associated nodal data. It is mentioned that this method is based on the works of Law and Yang, but we have applied the method to the singular Sturm-Liouville problem.

MSC:34L05, 45C05.

1 Introduction

Inverse problems of spectral analysis imply the reconstruction of a linear operator from some or other of its spectral characteristics. Such characteristics are spectra (for different boundary conditions), normalizing constants, spectral functions, scattering data, etc. An early important result in this direction, which gave vital impetus for further development of inverse problem theory, was obtained in [1]. At present, inverse problems are studied for certain special classes of ordinary differential operators. Inverse problems from two spectra are the most simple in their formulation and well studied in [2, 3]. An effective method of constructing a regular and singular Sturm-Liouville operator from a spectral function or from two spectra is given in [4–7].

We note that the details of the inverse problem for singular equations are given in the monographs [8–11] and references therein.

In some recent interesting works [12, 13], Hald and McLaughlin and Browne and Sleeman have taken a new approach to inverse spectral theory for the Sturm-Liouville problem. The novelty of these works lies in the use of nodal points as the given spectral data. In recent years, inverse nodal problems have been studied by several authors [14–21]etc.

In this paper, we deal with an inverse nodal problem for the Sturm-Liouville operator with Coulomb potential. We have reconstructed the potential function q from the nodal points of eigenfunctions, provided q is smooth enough. The method is based on a series of works by Law and Yang [14, 17].

Before giving the main results, we mention some physical properties of the Sturm-Liouville operator with Coulomb potential. Learning about the motion of electrons moving under the Coulomb potential is of significance in quantum theory. Solving these types of problems allows us to find energy levels not only for a hydrogen atom but also for single valence electron atoms such as sodium. For hydrogen atom, the Coulomb potential is given by U= − e 2 r , where r is the radius of the nucleus, e is electronic charge. According to this, we use the time-dependent Schrödinger equation

iħ ∂ Ψ ∂ t =− ħ 2 2 m ∂ 2 Ψ ∂ x 2 +U(x,y,z)Ψ, ∫ R 3 | Ψ | 2 dxdydz=1,

where Ψ is the wave function, ħ is Planck’s constant and m is the mass of electron. In this equation, if the Fourier transform is applied

it will convert to energy equation dependent on the situation as follows:

ħ 2 2 m ▽ 2 Ψ ˜ + U ˜ Ψ ˜ =E Ψ ˜ .

Therefore, energy equation in the field with the Coulomb potential becomes

− ħ 2 2 m ▽ 2 Ψ ˜ + ( E + e 2 r ) Ψ ˜ =0.

If this hydrogen atom is substituted to other potential area, then the energy equation becomes

− ħ 2 2 m ▽ 2 Ψ ˜ + ( E + e 2 r + q ( x , y , z ) ) Ψ ˜ =0.

If we make the necessary transformation, then we can get a Sturm-Liouville equation with Coulomb potential

− y ′ ′ + [ A x + q ( x ) ] y=λy,

where λ is a parameter which corresponds to the energy [22].

We consider the singular Sturm-Liouville problem

(1.1)
(1.2)
(1.3)

in which the function q(x)∈ L 1 [0,π], A, H are finite numbers and y ( x ) x ∈C[0,π]. Next, we denote by φ(x,s) the solution of (1.1) satisfying the initial condition

φ(0,s)=0, φ ′ (0,s)=s.
(1.4)

Let λ n be the n th eigenvalue and 0< x 1 n < x 2 n <⋯< x i n <π, i=1,2,…,n−1 be nodal points of the n th eigenfunction. Also, let I i n =[ x i n , x i + 1 n ] be the i th nodal domain of the n th eigenfunction and let l i n =| l i n |= x i + 1 n − x i n be the associated nodal length. We also define the function j n (x) by j n (x)=max{i: x i n <x}.

2 Main results

In this section, we try to obtain some asymptotic results and a reconstruction formula for the potential q, which has been obtained as a solution of an inverse nodal problem.

Lemma 2.1 The solution of problem (1.1)-(1.3) has the following form:

φ(x,s)=sinsx+ ∫ 0 x sin s ( x − t ) s { A t + q ( t ) } φ(t,s)dt,
(2.1)

where φ ( t , s ) t ∈C[0,π].

Proof Because φ(x,s) satisfies equation (1.1), we get

By integrating the first term twice on the right-hand side by parts and taking the conditions into account (1.2), we find that

φ(x,s)=sinsx+ ∫ 0 x sin s ( x − t ) s { A t + q ( t ) } φ(t,s)dt,

where φ ( t , s ) t ∈C[0,π]. □

Lemma 2.2 The eigenvalues of problem (1.1)-1.3) are the roots of (1.3). This spectral characteristic satisfies the following asymptotic expression [23]:

s n = λ n =n+ 1 2 + A 2 π ln ( n + 1 2 ) ( n + 1 2 ) + c 0 ( n + 1 2 ) +O ( ln n n 2 ) ,
(2.2)

where

Lemma 2.3 Assume that q∈ L 1 (0,π). Then, as n→∞,

(2.3)
(2.4)

Proof By using some iterations and trigonometric calculations in (2.1), we obtain

φ ( x , s ) = sin s x + sin s x 2 s ∫ 0 x sin 2 s t { A t + q ( t ) } d t − cos s x 2 s ∫ 0 x ( 1 − cos 2 s t ) { A t + q ( t ) } d t + o ( 1 s 2 ) .

If φ(x,s) is equal to zero and cosλx is not close to zero, then

Now, we take s= s n and x= x i n . Because Taylor’s expansion for the arctangent function is given by

arctanx=πi− ∑ k = 0 ∞ ( − 1 ) k + 1 x 2 k + 1 2 k + 1

for some integer i, then

s n x i n =πi+ 1 2 s n ∫ 0 x i n (1−cos2st) { A t + q ( t ) } dt+o ( 1 s n 2 ) .

Therefore

x i n = π i s n + 1 2 s n 2 ∫ 0 x i n (1−cos2st) { A t + q ( t ) } dt+o ( 1 s n 3 ) .

The nodal length is

l i n = x i + 1 n − x i n , l i n = π s n + 1 2 s n 2 ∫ x i n x i + 1 n (1−cos2st) { A t + q ( t ) } dt+o ( 1 s n 3 ) .

This completes the proof of Lemma 2.3. □

Lemma 2.4 Suppose f∈ L 1 (0,π). Then, for almost every x∈(0,π) with j= j n (x),

lim n → ∞ s n π ∫ x j n x j + 1 n f(t)dt=f(x).

Proof Since f∈ L 1 (0,π), d d x ∫ a x f(t)dt=f(x) almost everywhere. Thus, given any ζ>0, when n is sufficiently large and for almost every x∈(0,π),

This proves Lemma 2.4. □

Theorem 2.1 The potential function q(x)∈ L 1 (0,π) satisfies

q(x)= lim n → ∞ [ 2 s n 2 ( s n l j n π − 1 ) − s n A ln ( x j + 1 n x j n ) + s n A π ∫ x j n x j + 1 n cos 2 s n t t d t ]

for almost every x∈(0,π) with j= j n (x). We note that the asymptotic expression for s n in Theorem 2.1 implies that q(x)= lim n → ∞ F n (x).

Proof When we consider (2.4) in the form

l j n = π s n + 1 2 s n 2 ∫ x j n x j + 1 n (1−cos2st) { A t + q ( t ) } dt+o ( 1 s n 3 )

so that

By Lemma 2.4

lim n → ∞ s n π ∫ x j n x j + 1 n q(t)dt=q(x)

for almost every x∈(0,π).

It remains to show that for almost every x∈(0,π),

T n (x):= s n π ∫ x j n x j + 1 n cos2 s n tq(t)dt

tends to zero as n→∞. Take a sequence of continuous functions q k which converges to q in L 1 (0,π). Then q k has a subsequence converging to q almost everywhere in (0,π). We call this subsequence q k . Take any x such that q k (x) converges to q(x). Then for a given ε>0, we can fix a large k such that | q k (x)−q(x)|<ε. Hence

T n ( x ) = s n π ∫ x j n x j + 1 n cos 2 s n t [ q ( t ) − q k ( t ) ] d t + s n π ∫ x j n x j + 1 n cos 2 s n t [ q k ( t ) − q k ( x ) ] d t + s n π ∫ x j n x j + 1 n cos 2 s n t q k ( x ) d t = A n + B n + C n .

By Lemma 2.3,

C n = q k ( x ) 2 π [ sin ( 2 s n x j + 1 n ) − sin ( 2 s n x j n ) ] = q k (x)O ( 1 n )

and so it tends to zero as n→∞. By Lemma 2.4, the first term A n satisfies, when n is sufficiently large,

| A n | = | s n π ∫ x j n x j + 1 n cos 2 s n t [ q ( t ) − q k ( t ) ] d t | ≤ s n π ∫ x j n x j + 1 n | q ( t ) − q k ( t ) | d t < | q ( x ) − q k ( x ) | + ε < 2 ε .

On the other hand,

| B n |= | s n π ∫ x j n x j + 1 n cos 2 s n t [ q k ( t ) − q k ( x ) ] d t | ≤ s n π ∫ x j n x j + 1 n | q k ( t ) − q k ( x ) | dt.

Because q k is continuous, this term is arbitrarily every x∈(0,π). Hence we conclude that lim n → ∞ T n (x)=0. This proves Theorem 2.1. □

Lemma 2.5 We take a sequence f k ∈C[0,π] converges to f∈ L 1 , then, for any large enough n, with j= j n (x) as k→∞

∥ s n π ∫ x j n x j + 1 n [ f k ( t ) − f ( t ) ] d t ∥ 1 →0.

Proof By (2.4) and observation that the integral ∫ x j n x j + 1 n [ f k (t)−f(t)]dt is constant on any nodal interval I j n , we obtain

and for k→∞ this term converges to zero. □

Lemma 2.6 Suppose that q∈ L 1 (0,π), then as n→∞ with j= j n (x),

∥ s n π ∫ x j n x j + 1 n q ( t ) d t − q ( x ) ∥ 1 →0.

Proof Firstly, let us show that if q is continuous on [0,π], the result is satisfied. Let M= max x ∈ [ 0 , π ] |q(x)|. By using the intermediate value theorem, there exists ξ∈(a,x) such that

| 1 x − a ∫ a x q ( t ) d t − q ( x ) | = | q ( ξ ) − q ( x ) | .

If x is close enough to a, the difference can be arbitrarily small. Then, for all ε>0, when n is large enough, with j= j n (x) we get

| s n π ∫ x j n x j + 1 n q ( t ) d t − q ( x ) | ≤ | s n ( x − x j n ) π [ 1 x − x j n ∫ x j n x q ( t ) d t − q ( x ) ] | + | s n ( x j + 1 n − x ) π [ 1 x j + 1 n − x ∫ x x j + 1 n q ( t ) d t − q ( x ) ] | + | q ( x ) | | ( s n l j n π − 1 ) | ≤ 2 s n l j n ε π + M ε ≤ ( M + 2 + 2 ε ) ε .

In the above process, we assume that x≠ x j n . The estimate also holds if x= x j n . Hence if q∈C[0,π], s n π ∫ x j n x j + 1 n q(t)dt converges to q(x) uniformly on (0,π). Thus ∥ s n π ∫ x j n x j + 1 n q(t)dt−q(x)∥ can be arbitrarily small. Because C[0,π] is dense in L 1 (0,π), for any q∈ L 1 (0,π), there exists a sequence q k ∈C[0,π] convergent to q in L 1 (0,π). Hence, fix n sufficiently large,

From the above process and Lemma 2.5, when k is large enough, the first two terms are arbitrarily small. Hence, as k→∞,

∥ s n π ∫ x j n x j + 1 n [ f k ( t ) − f ( t ) ] d t ∥ 1 →0.

 □

Theorem 2.2 F n converges to q in L 1 .

Proof When we consider the value of F n , we obtain that

It suffices to show that as n→∞

∥ s n ( 2 s n 2 l j n π − 2 s n − A ln ( x j + 1 n x j n ) + A π ∫ x j n x j + 1 n cos 2 s n t t d t ) − q ∥ 1 →0.

By using (2.4) we have

Hence, we only need to prove that for n→∞

∥ s n π ∫ x j n x j + 1 n q ( t ) d t − q ∥ 1 →0

and

∥ s n π ∫ x j n x j + 1 n cos 2 s n t { A t + q ( t ) } d t ∥ 1 →0.
(2.5)

From Lemma 2.6, the first limit holds and the second limit also holds. On the other hand, the sequence of functions

c n (x)= s n π ∫ x j n x j + 1 n cos2 s n t { A t + q ( t ) } dt

converges to 0 for almost every x∈(0,π). Furthermore,

| c n ( x ) | ≤ s n π ∫ x j n x i + 1 n | A t + q ( t ) | dt= g n (x)

and

∫ 0 π g n (x)dx= ∑ i = 0 n − 1 s n l j n π ∫ x j n x j + 1 n | A t + q ( t ) | dt= [ 1 + O ( ln n n ) ] ∥ q ∥ 1 .

Then, we may apply the Lebesque dominated convergence theorem to show that (2.5) is valid. The proof of Theorem 2.2 is completed. □

References

  1. Ambartsumyan VA: Ãœber eine frage der eigenwerttheorie. Z. Phys. 1929, 53: 690-695. 10.1007/BF01330827

    Article  Google Scholar 

  2. Levitan BM: On the determination of the Sturm-Liouville operator from one and two spectra. Math. USSR, Izv. 1978, 12: 179-193. 10.1070/IM1978v012n01ABEH001844

    Article  Google Scholar 

  3. Isaacson EL, Trubowitz E: The inverse Sturm-Liouville problem. I. Commun. Pure Appl. Math. 1983, 36: 767-783. 10.1002/cpa.3160360604

    Article  MathSciNet  Google Scholar 

  4. Gelfand IM, Levitan BM: On the determination of a differential equation by its spectral function. Izv. Akad. Nauk SSSR, Ser. Mat. 1951, 15: 309-360. Ams, 253-304 (1955)

    MathSciNet  Google Scholar 

  5. Hochstadt H: The inverse Sturm-Liouville problem. Commun. Pure Appl. Math. 1973, 26: 715-729. 10.1002/cpa.3160260514

    Article  MathSciNet  Google Scholar 

  6. Pöschel J, Trubowitz E: Inverse Spectral Theory. Academic Press, Boston; 1987.

    Google Scholar 

  7. Rundell W, Sack EP: Reconstruction of a radially symmetric potential from two spectral sequences. J. Math. Anal. Appl. 2001, 264: 354-381. 10.1006/jmaa.2001.7664

    Article  MathSciNet  Google Scholar 

  8. Carlson R: Borg-Levinson theorem for Bessel operator. Pac. J. Math. 1997, 177: 1-26. 10.2140/pjm.1997.177.1

    Article  Google Scholar 

  9. Chadan K, Colton D, Paivarinta L, Rundell W: An Introduction to Inverse Scattering and Inverse Spectral Problems. SIAM, Philadelphia; 1997.

    Book  Google Scholar 

  10. Panakhov ES, Sat M: On the determination of the singular Sturm-Liouville operator from two spectra. Comput. Model. Eng. Sci. 2012, 84: 1-11.

    MathSciNet  Google Scholar 

  11. Hald OH: Discontinuous inverse eigenvalue problem. Commun. Pure Appl. Math. 1984, 37: 539-577. 10.1002/cpa.3160370502

    Article  MathSciNet  Google Scholar 

  12. Browne PJ, Sleeman BD: Inverse nodal problems for Sturm-Liouville equation with eigenparameter dependent boundary conditions. Inverse Probl. 1996, 12: 377-381. 10.1088/0266-5611/12/4/002

    Article  MathSciNet  Google Scholar 

  13. Hald OH, McLaughlin JR: Solution of inverse nodal problems. Inverse Probl. 1989, 5: 307-347. 10.1088/0266-5611/5/3/008

    Article  MathSciNet  Google Scholar 

  14. Chen YT, Cheng YH, Law CK, Tsa J: Convergence of reconstruction formula for the potential function. Proc. Am. Math. Soc. 2002, 130: 2319-2324. 10.1090/S0002-9939-02-06297-4

    Article  Google Scholar 

  15. Yang FX: A solution of the inverse nodal problem. Inverse Probl. 1997, 13: 203-213. 10.1088/0266-5611/13/1/016

    Article  Google Scholar 

  16. McLaughlin JR: Inverse spectral theory using nodal points as a data - a uniqueness result. J. Differ. Equ. 1988, 73: 354-362. 10.1016/0022-0396(88)90111-8

    Article  Google Scholar 

  17. Law CK, Shen CL, Yang CF: The inverse nodal problem on the smoothness of the potential function. Inverse Probl. 1999, 15: 253-263. 10.1088/0266-5611/15/1/024

    Article  MathSciNet  Google Scholar 

  18. Yurko VA, Freiling G: Inverse nodal problems for differential operators on graphs with a cycle. Tamkang J. Math. 2010, 41: 15-24.

    MathSciNet  Google Scholar 

  19. Yang CF: Inverse nodal problems for the Sturm-Liouville operator with eigenparameter dependent boundary conditions. Oper. Matrices 2012, 6(1):63-77.

    Article  MathSciNet  Google Scholar 

  20. Koyunbakan H, Panakhov ES: A uniqueness theorem for inverse nodal problem. Inverse Probl. Sci. Eng. 2007, 15: 517-524. 10.1080/00423110500523143

    Article  MathSciNet  Google Scholar 

  21. Koyunbakan H: Reconstruction of potential function for diffusion operator. Numer. Funct. Anal. Optim. 2009, 30: 111-120. 10.1080/01630560802279256

    Article  MathSciNet  Google Scholar 

  22. Blohincev DI: Foundations of Quantum Mechanics. GITTL, Moscow; 1949.

    Google Scholar 

  23. Amirov RK, Çakmak Y, Gulyaz S: Boundary value problem for second order differential equations with Coulomb singularity on a finite interval. Indian J. Pure Appl. Math. 2006, 37: 125-140.

    MathSciNet  Google Scholar 

Download references

Acknowledgements

The authors would like to thank the editor and referees for their valuable comments and remarks which led to a great improvement of the article.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Murat Sat.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

MS wrote the first draft and ESP corrected and improved the final version. Both authors read and approved the final draft.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Panakhov, E.S., Sat, M. Reconstruction of potential function for Sturm-Liouville operator with Coulomb potential. Bound Value Probl 2013, 49 (2013). https://doi.org/10.1186/1687-2770-2013-49

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2013-49

Keywords