Arxiv:1802.04761V1 [Math.SP] 13 Feb 2018 Eesr N Ucetcniin O T Solvability

Home , Inverse problem

arXiv:1802.04761v1 [math.SP] 13 Feb 2018 eesr n ucetcniin o t solvability. inv th its the of for of part conditions solution a sufficient the using and for part, necessary algorithm kn remaining an is the provide kernel theorem, on convolution uniqueness it the recover that and interval, suppose the We considered. is volution e,nnoa operator. nonlocal lem, h olwn nes pcrlpolmhsbe netgtdi [13]. in investigated been has problem spectral inverse following The sy Dirac the pro for inverse and of an [10–12]) solution (see of delay for integral reduction applied an the with successfully ser operators on containing been equation, has based integral method proposed, nonlinear this a was to method operator new nonlocal obta a were [9], operators integro-differential differe paper for than investigation, problems for inverse difficult on rec more In ma therefore of and [1–5]. attention nonlocal monographs much attract the f operators in problems integro-differential found for spectral be inverse can on operators differential results Classical characteristics. 47E05 where h functions the au rbe o h ytm()wt h onayconditions boundary the with (1) system the for problem value nivrepolmfrteitgodffrnilDrcsyste Dirac integro-differential the for problem inverse An Abstract. Keywords: M ahmtc ujc lsicto (2010): Classification Subject Mathematics AMS nes rbe 1. Problem Inverse nti ae,w td nivreseta rbe o h integr the for problem spectral inverse an study we paper, this In Introduction 1. h anrslso 1]cnb omltda follows. as formulated be can [13] of results main The edsrb h eut fteppr[3 nmr eal Let detail. more in [13] paper the of results the describe We oper differential recovering in consist problems spectral Inverse λ steseta parameter, spectral the is ata nomto ie ntecnouinkernel convolution the on given information partial p , nitgodffrnilDrcsse iha nerltr ntefor the in term integral an with system Dirac integro-differential An B q ner-ieeta ia ytm nes pcrlpolm half problem, spectral inverse system, Dirac integro-differential = r ope-audadbln oteclass the to belong and complex-valued are − L 1 0 0 1 2 By ,π ie h pcrmof spectrum the Given ( ,π c, ′ + M , := ) Z 0 x aai .Bondarenko P. Natalia M { ( r x ( ( : x y = ) 1 − 0 = (0) π − t ) y − x p ( ) q t y ( 1 r ) ( x 1 ( x ( dt ) x D π ) ) 0 = ) = find , ∈ p q ( ( λy, L x x 2 ) ) . ( p ,π c, y , and 0 ) 45 40 40 41 34L40 34B10 34B09 34B05 34A55 π, < x < } hmtcas uhoeaosare operators Such thematicians. L c , ( ndi h ok 68.I the In [6–8]. works the in ined q 2 x tmi h om()(e [13]). (see (1) form the in stem ,π rSumLovleadDirac and Sturm-Liouville or . D = ) e fcnouin.Ltron, Later convolutions. of ies ta prtr.Frtresults First operators. ntial (0 n er,ivreproblems inverse years, ent = ≥ repolmtgte with together problem erse -ieeta ia system Dirac o-differential π , pcrm epoethe prove We spectrum. e lm o Sturm-Liouville for blems tr rmterspectral their from ators w roio atof part a on priori a own D 0 ), y y . ( 1 2 nes rbe o a for problem inverse ,q p, ( ( x x ) ) eteboundary the be ) , ivreprob- -inverse with m fcon- of m (2) (1) Theorem 1. The problem D has a countable set of eigenvalues, which can be numbered as {λk}k∈Z counting with their multiplicities, so that

λk = k + κk, {κk} ∈ l2. (3)

For any sequence of complex numbers {λk}k∈Z of the form (3), there exist unique functions p, q ∈ L2,π(0, π) such that {λk}k∈Z is the spectrum of the boundary value problem D(p, q) of the form (1)-(2). The proof of Theorem 1 is constructive, and gives an algorithm for the solution of Inverse Problem 1. In this paper, we suppose that the kernel of the integral term M(x) is known a priori on the π subinterval (0, a), 2 ≤ a < π. We show that in this case, only a fractional part of the spectrum is sufficient to reconstruct the matrix M(x) on the subinterval (a, π). Such problem statement is similar in some sense to the Hochstadt-Lieberman problem [14], which consists in recovering the potential of the Sturm-Liouville operator on a finite interval from a spectrum, while the potential on the half of the interval is known a priori. The Hochstadt-Lieberman problem, also called the half-inverse problem, and its generalizations for various classes of differential operators and pencils were studied in [15–29]. However, we know the only paper [30] on a half- inverse problem for an integro-differential operator, where the Sturm-Liouville operator with an integral delay in the form of convolution is recovered from a half of the spectrum and the convolution kernel is given a priori on the half of the interval. We also mention the paper [31], where an inverse problem for the integro-differential Sturm-Liouville operator on a graph was studied under the assumption, that the operator coefficients are known on a part of the graph. The paper is organized as follows. In Section 2, we provide the system of the main equations for Inverse Problem 1, derived in [13], and show the benefits of the knowledge of the functions p and q on a subinterval. In Section 3, the partial inverse problem is formulated, which consists in recovering the functions p and q on the remaining part of the interval from an abstract subset of the spectrum {λk}k∈I, I ⊂ Z. We prove the uniqueness theorem in terms of completeness of some system of vector functions, associated with the given subspectrum. We also provide a constructive algorithm for the solution of the partial inverse problem. In Section 4, we consider π N Z the case a = π − m , m ∈ , m ≥ 2, I = {sm: s ∈ }. In this case the partial inverse problem is uniquely solvable. Using the results of Section 3, we obtain the necessary and sufficient conditions for the solvability of the inverse problem in this special case. 2. Main equations In this section, we modify the system of the main equations from the paper [13], assuming that the functions p and q are known on the subinterval (0, a). Let us introduce the following notations for the convolutions: x (f ∗ g)(x)= f(t)g(x − t) dt, f ∗n(x)= f ∗ f ∗···∗ f(x), f ∗0(x)= δ(x), 0 Z n where δ is the Dirac delta function: δ ∗ f = f ∗ δ = |f. {z } In the paper [13], Inverse Problem 1 was reduced to the system of nonlinear integral equations with respect to the functions p and q:

∞ n a (π − t)n −w (π − t)=(π − t)q(t) − nj p∗j ∗ q∗(n−j) (t), 1 n!  n=2 j=0  X∞ Xn (4)  b (π − t)n −w (π − t)=(π − t)p(t)+ nj p∗j ∗ q∗(n−j) (t). 2 n! n=2 j=0  X X   2 Here anj and bnj are certain coeﬃcients, satisfying the estimates

n n n−1 n−1 |anj| ≤ 2 , |bnj| ≤ 2 , n ∈ N, j=0 j=0 X X and the functions w1, w2 can be found from the following representation of the characteristic function ∆(λ), whose zeros coincide with the eigenvalues of D:

π ∆(λ) = sin λπ + w1(t) sin λt + w2(t) cos λt dt, w1,w2 ∈ L2(0, π). (5) 0 Z Suppose that the matrix-function M(x) is known on the interval (0, a). Introduce the functions p(x), x ∈ (0, a), q(x), x ∈ (0, a), p (x)= q (x)= (6) 1 0, x ∈ (a, π), 1 0, x ∈ (a, π), 0, x ∈ (0, a), 0, x ∈ (0, a), p (x)= q (x)= 2 p(x), x ∈ (a, π), 2 q(x), x ∈ (a, π).

Thus, p(x)= p1(x)+ p2(x), q(x)= q1(x)+ q2(x). One can easily show that

∗j ∗(n−j) ∗j ∗(n−j) ∗(j−1) ∗(n−j) ∗j ∗(n−j−1) (p1 + p2) ∗ (q1 + q2) = p1 ∗ q1 + jp1 ∗ q1 ∗ p2 +(n − j)p1 ∗ q1 ∗ q2, so the system (4) takes the form

∞ (π − t)n −w (π − t)=(π − t)q (t)+(π − t)q (t)+ A (t)+ (A ∗ p + A ∗ q )(t), 1 1 2 1 n! n1 2 n2 2  n=2  X∞  (π − t)n −w (π − t)=(π − t)p (t)+(π − t)p (t)+ B (t)+ (B ∗ p + B ∗ q )(t), 2 1 2 1 n! n1 2 n2 2 n=2  X (7)  where

∞ n n a (π − t) ∗ − A (t)= − nj p∗j ∗ q (n j) (t), (8) 1 n! 1 1 n=2 j=0 X X n n−1 ∗(j−1) ∗(n−j) ∗j ∗(n−j−1) An1(t)= − janj p1 ∗ q1 (t), An2(t)= − (n − j)anj p1 ∗ q1 (t), j=1 j=0 X X (9) ∞ n n b (π − t) ∗ − B (t)= nj p∗j ∗ q (n j) (t), (10) 1 n! 1 1 n=2 j=0 X X n n−1 ∗(j−1) ∗(n−j) ∗j ∗(n−j−1) Bn1(t)= jbnj p1 ∗ q1 (t), Bn2(t)= (n − j)bnj p1 ∗ q1 (t). (11) j=1 j=0 X X

For arbitrary functions f,g ∈ L2(0, π) their convolution f ∗ g also belong to L2(0, π), and

kf ∗ gk ≤ πkfkkgk, kf ∗nk ≤ πn−1kfkn, n ∈ N, where k.k is the norm in L2(0, π). Consequently, the series A1 and B1 converge in L2(0, π), and the functions Ank, Bnk satisfy the estimate

n−2 n−1 n−1 kAnkk, kBnkk ≤ nπ 2 max{kp1k, kq1k} , n ≥ 2, k =1, 2. (12)

3 For t ∈ (0, a), the system (7) takes the form

−w (π − t)=(π − t)q (t)+ A (t), 1 1 1 (13) −w (π − t)=(π − t)p (t)+ B (t), 2 1 1

Thus, known the functions p, q ∈ L2(0, a), one can ﬁnd w1,w2 ∈ L2(b, π), b := π − a. For t ∈ (a, π), we have the system (7) with respect to the functions p2 and q2.

Lemma 1. For arbitrary functions w1,w2 ∈ L2(0, b), A1, B1 ∈ L2(a, π) and Ank, Bnk ∈ L2(0, π), n ≥ 2, k = 1, 2, satisfying the estimate (12), the system (7) has a unique solution p2, q2 ∈ L2,π(a, π), i.e. (π − t)p2(t), (π − t)q2(t) ∈ L2(a, π). Proof. Denote

q (t) −w (π − t) − A (t) y(t)=(π − t) 2 , µ(t)= 1 1 , p2(t) −w2(π − t) − B1(t) ∞ π − t (π − t)n−1 A (t − s) A (t − s) H(t, s)= n1 n2 . π − s n! Bn1(t − s) Bn2(t − s) n=2 X Then the system (7) for t ∈ (a, π) can be rewritten in the form

t y(t)+ H(t, s)y(s) ds = µ(t). (14) Za

According to the conditions of the theorem, µ ∈ L2(a, π)⊕L2(a, π) and the entries of the matrix function H(t, s) belong to L2 ((a, π) × (a, π)). Thus, the Volterra integral equation (14) of the second kind has a unique solution y ∈ L2(a, π) ⊕ L2(a, π).

3. Abstract partial inverse problem

Let I be a subset of Z. In this section, we study the following partial inverse problem. Inverse Problem 2. Given the functions p, q on the interval (0, a) and the subset of the spectrum {λk}k∈I, construct p and q on the interval (a, π). Of course, the uniqueness of the solution and the possibility of the constructive solution of Inverse Problem 2 depend on a and a subset I. Note that the subset {λk}k∈I can contain multiple eigenvalues. Denote by {λk}k∈J the subset of all the distinct eigenvalues among {λk}k∈I, and let M := {mk}k∈J be the corresponding multiplicities, mk ≥ 1. The analysis in the previous section shows, that we need to know w1 and w2 on the interval (0, b), in order to ﬁnd p2 and q2. Introduce the function

π E(λ) := ∆(λ) − sin λπ − (w1(t) sin λt + w2(t) cos λt) dt. (15) Zb Relation (5) implies b E(λ)= (w1(t) sin λt + w2(t) cos λt) dt. (16) Z0 It follows from (15) and (16), that

dj dj dj π E(λ )= − sin λπ − (w (t) sin λt + w (t) cos λt) dt , (17) dλj k dλj dλj 1 2 λ=λk Zb λ=λk

4 dj b dj dj E(λ )= w (t) sin λt + w (t) cos λt dt (18) dλj k 1 dλj 2 dλj Z0 λ=λk λ=λk ! dj for k ∈J , j = 0, m − 1, since ∆(λ ) = 0 for such k. k dλj k Consider the complex Hilbert space H := L2(0, b) ⊕ L2(0, b) with the scalar product

b g1 h1 (g, h)= (g1(t)h1(t)+ g2(t)h2(t)) dt, g, h ∈ H, g = , h = . g2 h2 Z0

In view of (18), the functions w1 and w2 on (0, b) are uniquely specified by the values of dj E(λ ), k ∈J , j = 0, m − 1, if the system of vector functions dλj k k dj sin λt dλj SI :=    (19) dj  λt   j cos   dλ λ=λk k∈J ,j=0,mk−1     dj   is complete in H. The values E(λ ) can be determined by (17). dλj k Along with D, we consider another boundary value problem D˜ = D(˜p, q˜) of the same form, but with different coefficients. We agree that if a certain symbol γ denotes an object related to D, the corresponding symbolγ ˜ denotes an analogous object related to D˜. Now we are ready to formulate a uniqueness theorem for Inverse Problem 2.

Theorem 2. Let p(x)=˜p(x) and q(x)=˜q(x) a.e. on (0, a), and let the eigenvalues {λk}k∈Z and {λ˜k}k∈Z be numbered in such a way, that the asymptotic formula (3) holds for the both sequences, and λk = λ˜k for k ∈ I, where the subset I is such that the system SI is complete in H. Then p(x)=˜p(x) and q(x)=˜q(x) a.e. on (a, π).

Proof. Since p(x)=˜p(x) and q(x)=˜q(x) a.e. on (0, a), we have A1(t)= A˜1(t), B1(t)= B˜1(t), Ank(t) = A˜nk(t), Bnk(t) = B˜nk(t) for n ≥ 2, k = 1, 2, t ∈ [0, π). The relations (13) yield w1(x)=w ˜1(x), w2(x)=w ˜2(x) a.e. on (b, π). By virtue of (17) and the equality λk = λ˜k, dj dj k ∈ I, we have E(λ ) = E˜(λ˜ ) for k ∈ J , j = 0, m − 1. Since the system SI is dλj k dλj k k complete in H, we obtain w1(t)=w ˜1(t) and w2(t)=w ˜2(t) a.e. on (0, b). Applying Lemma 1, we arrive at the assertion of the theorem.

If the system SI is a Riesz basis in H, one can solve Inverse Problem 2 by the following algorithm.

Algorithm 1. Let the functions p(x) and q(x) for x ∈ (0, a) and the eigenvalues {λk}k∈I be given.

1. Construct the functions p1, q1 by (6) and then A1, B1 by (8), (10).

2. Find w1(t) and w2(t) for t ∈ (b, π) by formulas (13). dj 3. Find E(λ ), k ∈J , j = 0, m − 1, using (17). dλj k k

4. Construct the system of vector functions SI , and determine the functions w1(t) and w2(t) dj from the coordinates E(λ ) with respect to the Riesz basis SI (see (18)). dλj k

5 5. Construct the functions Ank, Bnk, n ≥ 2, k =1, 2 by (9), (11).

6. Solving the system of Volterra integral equations (7), ﬁnd the functions p2(x)= p(x) and q2(x)= q(x) for x ∈ (a, π). 4. Speciﬁc partial inverse problem

The formulation of the results in the previous section contain abstract conditions of the completeness or the Riesz-basicity of the system SI . In this section, we show that these conditions hold for some certain subsets I with the certain a. N π Z Let m ∈ , m ≥ 2. Put a := π − m , I := {sm: s ∈ }. Then the following result is valid.

Lemma 2. For an arbitrary sequence {λk}k∈I = {λsm}s∈Z of not necessarily distinct complex numbers, satisfying the relation

λsm = sm + κs, {κs} ∈ l2, (20) the system SI , deﬁned by (19), is a Riesz basis in H.

Proof. According to the asymptotic relation (20), the system SI is l2-close to the system sin smt , cos smt s∈Z which is an orthogonal basis in H. Thus, it remains to prove, that SI is complete in H. Suppose h1, h2 ∈ L2(0, b) are such functions, that

b dj dj h (t) sin λt + h (t) cos λt dt =0, k ∈J , j = 0, m − 1. 1 dλj 2 dλj k Z0 λ=λk λ=λk ! This means that the entire function b H(λ) := (h1(t) sin λt + h2(t) cos λt) dt Z0 has zeros λk with the multiplicities mk for k ∈J . On the other hand, in view of (20) and Theorem 1, the sequence {λsm}s∈Z is the spectrum π of some boundary value problem D of the form (1)-(2) on the interval 0, m instead of (0, π). Its characteristic function admits the following representation (see (5)): b ∆(λ) = sin λb + (w1(t) sin λt + w2(t) cos λt) dt, w1,w2 ∈ L2(0, b). Z0 H(λ) Clearly, the function is entire and bounded in the whole λ-plane. By the Liouville’s ∆(λ) theorem, H(λ) ≡ C∆(λ), where C is a constant. However, lim H(λ) = 0, while lim ∆(λ) λ→+∞ λ→+∞ does not exist. Thus, C = 0 and H(λ) ≡ 0. Consequently, h1(t)= h2(t)=0 a.e. on (0, b), and the system SI is complete. The following theorem provides suﬃcient (and also necessary) conditions for the solvability π Z of Inverse Problem 2 for a = π − m , I = {sm: s ∈ }. π Theorem 3. Let p and q be arbitrary functions from L2(0, a), a = π − m , and let {λk}k∈I = {λsm}s∈Z be arbitrary complex numbers, satisfying (20). Then the functions p and q can be uniquely continued to the interval (a, π), so that p, q ∈ L2,π(0, π) and {λk}k∈I are eigenvalues of the boundary value problem D = D(p, q) with the corresponding multiplicities.

6 Proof. Since the system SI , constructed from the given {λk}k∈I, is a Riesz basis by Lemma 2, one can construct the functions p, q ∈ L2,π(a, π), using Algorithm 1. By construction, {λk}k∈I are eigenvalues of the problem D = D(p, q). The uniqueness of D follows from Theorem 2. For m = 2, we get a half inverse problem, when the functions p and q are known a priori π π Z on the interval 0, 2 and can be recovered on 2 , π by the half of the spectrum {λ2k}k∈ . Acknowledgment. This work was supported by Grant 17-11-01193 of the Russian Science Foundation.

References

[1] Marchenko, V. A. Sturm-Liouville Operators and their Applications, Naukova Dumka, Kiev (1977) (Russian); English transl., Birkhauser (1986).

[2] Levitan, B. M. Inverse Sturm-Liouville Problems, Nauka, Moscow (1984) (Russian); En- glish transl., VNU Sci. Press, Utrecht (1987).

[3] Levitan, B. M.; Sargsyan I. S. Sturm-Liouville and Dirac Operators, Nauka, Moscow (1988); English transl., Kluwer Academic Publishers, Dordrecht (1991).

[4] P¨oschel, J.; Trubowitz, E. Inverse Spectral Theory, New York, Academic Press (1987).

[5] Freiling, G.; Yurko, V. Inverse Sturm-Liouville problems and their applications. Hunting- ton, NY: Nova Science Publishers, 305 p. (2001).

[6] Malamud, M. M. On some inverse problems, Boundary Value Problems of Mathematical Physics, Kiev (1979), 116-124.

[7] Eremin, M. S. An inverse problem for a second-order integro-diﬀerential equation with a singularity, Diﬀer. Uravn. 24:2 (1988), 350–351.

[8] Yurko, V. A. Inverse problem for integrodiﬀerential operators, Mathematical Notes 50:5 (1991), 1188–1197.

[9] Buterin, S. A. The inverse problem of recovering the Volterra convolution operator from the incomplete spectrum of its rank-one perturbation, Inverse Problems 22 (2006), 2223–2236.

[10] Buterin, S. A. On an inverse spectral problem for a convolution integro-diﬀerential operator, Res. Math. 50 (2007), no.3–4, 73–181.

[11] Buterin, S. A. On the reconstruction of a convolution perturbation of the Sturm-Liouville operator from the spectrum, Diﬀ. Uravn. 46 (2010), 146–149 (Russian); English transl. in Diﬀ. Eqns. 46 (2010), 150–154.

[12] Buterin, S. A.; Choque Rivero, A. E. On inverse problem for a convolution integrodiﬀer- ential operator with Robin boundary conditions, Appl. Math. Lett. 48 (2015) 150–155.

[13] Bondarenko, N.; Buterin, S. On Recovering the Dirac Operator with an Integral Delay from the Spectrum, Results Math. 71 (2017), 1521–1529.

[14] Hochstadt, H.; Lieberman, B. An inverse Sturm-Liouville problem with mixed given data, SIAM J. Appl. Math. 34 (1978), 676–680.

7 [15] Gesztesy, F.; Simon, B. Inverse spectral analysis with partial information on the potential, II. The case of discrete spectrum, Trans. AMS 352:6 (2000), 2765–2787.

[16] Sakhnovich, L. Half-inverse problems on the ﬁnite interval, Inverse Problems 17 (2001), 527–532.

[17] Hryniv, R.O.; Mykytyuk Ya. V. Half-inverse spectral problems for SturmLiouville operators with singular potentials, Inverse Probl. 20 (2004), 1423–1444.

[18] Shieh, C.-T.; Yurko, V. A. Inverse nodal and inverse spectral problems for discontinuous boundary value problems, J. Math. Anal. Appl. 347 (2008), 266-272.

[19] Martinyuk, O.; Pivovarchik, V. On the Hochstadt-Lieberman theorem, Inverse Problems 26 (2010), 035011 (6pp).

[20] Shieh, C.-T.; Buterin, S. A.; Ignatiev, M. Yu. On Hochstadt-Lieberman theorem for Sturm- Liouville operators, Far East J. Appl. Math. 52:2 (2011), 131–146.

[21] Buterin, S. A. On half inverse problem for diﬀerential pencils with the spectral parameter in boundary conditions, Tamkang J. Math. 42:3 (2011), 355–364.

[22] Buterin, S.A.; Shieh C.-T. Incomplete inverse spectral and nodal problems for diﬀerential pencils, Results Math. 62:1 (2012), 167–179.

[23] Yang, C.-F. Inverse spectral problems for the Sturm-Liouville operator on a d-star graph, J. Math. Anal. Appl. 365 (2010), 742–749.

[24] Yang, C.-F. Hochstadt–Lieberman theorem for Dirac operator with eigenparameter de- pendent boundary conditions, Nonlinear Analysis: Theory, Methods and Applications 74:7 (2011), 2475–2484.

[25] Yang, C.-F.; Zettl, A. Half inverse problems for quadratic pencils of Sturm-Liouville operators, Taiwanese J. Math. 16:5 (2012), 1829–1846.

[26] Wei, G.; Xu, H.-K. On the missing eigenvalue problem for an inverse SturmLiouville problem, Journal de Mathematiques Pures et Appliquees, 91:5 (2009), 468–475.

[27] Wei, G.; Xu, H.-K. Inverse spectral problem with partial information given on the potential and norming constants, Trans. Amer. Math. Soc. 364 (2012), 3265–3288.

[28] Bondarenko, N.P. Partial inverse problems for the Sturm-Liouville operator on a star- shaped graph with mixed boundary conditions, J. Inverse and Ill-Posed Probl. (2018), Vol. 26, Issue 1. Pp. 1–12.

[29] Bondarenko, N.P. A partial inverse problem for the Sturm-Liouville operator on a star- shaped graph, Analysis and Mathematical Physics, published online 24 April 2017, pp. 1-14. DOI: 10.1007/s13324-017-0172-x

[30] Buterin, S.; Sat, M. On the half inverse spectral problem for an integro-diﬀerential operator, Inverse Problems in Science and Engineering (2017), Vol 25, Issue 10. Pp. 1508–1518.

[31] Bondarenko, N.P. An inverse problem for an integro-diﬀerential operator on a star-shaped graph, Mathematical Methods in the Applied Sciences (2018), Vol. 41, Issue 4. Pp. 1697- 1702.

8 Natalia Pavlovna Bondarenko 1. Department of Applied Mathematics, Samara National Research University, Moskovskoye Shosse 34, Samara 443086, Russia, 2. Department of Mechanics and Mathematics, Saratov State University, Astrakhanskaya 83, Saratov 410012, Russia, e-mail: [email protected]