arxiv: v1 [math.sp] 4 Sep 2015

Partial inverse problems for Sturm-Liouville operators on trees Natalia Bondarenko, Chung-Tsun Shieh arxiv:1509.01534v1 [math.sp] 4 Sep 2015 Abstract. In this paper, inverse spectral problems for Sturm-Liouville operators on a tree (a graphwithout cycles) arestudied. Weshow thatif thepotential onanedgeisknown apriori, then b 1 spectral sets uniquely determine the potential functions on a tree with b external edges. Constructive solutions, based on the method of spectral mappings, are provided for the considered inverse problems. Keywords: quantum graphs; Sturm-Liouville operators; inverse spectral problems; method of spectral mappings. AMS Mathematics Subject Classification (2010): 34A55 47E05 34B45 34L40 1. Introduction. This paper concerns the theory of inverse spectral problems for Sturm-Liouville operators on geometrical graphs. Inverse problems consist in recovering differential operators from their spectral characteristics. Differential operators on graphs (quantum graphs) have applications in various fields of science and engineering(mechanics, chemistry, electronics, nanoscale technology and others) and attract a considerable attention of mathematicians in recent years. There is an extensive literature devoted to differential operators on graphs and their applications, we mention only some research papers and surveys [1, 2, 3, 4, 5, 6, 7]. There are different kinds of inverse problems studied for quantum graphs, one of them is to recover the coefficients of the operator while some information is known a priori. This paper is focused on the reconstruction of the potential of the Sturm-Liouville operator on a tree(a graph without cycles) with a prescribed structure and standard matching conditions in the vertices. V.A. Yurko [8, 9] studied such inverse problems on trees by the Weyl vector, the system of spectra and the spectral data. These problems are generalizations of the well studied inverse problems for Sturm-Liouville operators on a finite interval (see monographs [10, 11, 12, 13] and references therein). By the method of spectral mappings [13, 14], V.A. Yurko proved uniqueness theorems and developed a constructive algorithm for solution of inverse problems on trees. Later on, this solution was generalized to the case of arbitrary graphs [7]. In this paper, we formulate and solve partial inverse problems for the Sturm-Liouville operator on the tree. We suppose that the Sturm-Liouville potential is known on the part of the graph and show that we need less data to recover the potential on the remaining part. We know the only work [15] in this direction, where the potential is known on a half of one edge and completely on the other edges of the star-shaped graph, and the author solves the Hochstadt-Lieberman-type problem [16] by a part of the spectrum. In this paper, we assume that the potential is known on one edge of a tree, then reconstruct the potential on the remaining part by the system of spectra or the Weyl functions. By developing the ideas of V.A. Yurko [8, 9], we show that one needs one less spectral set or one less Weyl function for the solution of the partial inverse problem. We consider separately the cases of boundary and internal edges, present constructive solutions and corresponding uniqueness theorems for both of them. The paper is organized as follows. In Section 2, we introduce the notation and briefly describe the solution of inverse problems on trees by V.A. Yurko [8, 9]. In Section 3, we formulate our main results and outline their constructive solutions. Section 4 contains proofs of the technical lemmas from Section 3. In Section 5 we illustrate our method by an example. 1

2. Inverse problems on a tree In this section, we introduce the notation and provide the main results of V.A. Yurko on the inverse problems on trees (see works [8, 9] for more details). Consider a compact tree G with the vertices V = {v i } m+1 i=1 and edges E = {e j } m j=1. For each vertex v V, we denote the set of edges associated with v by E v and call the size of E v the degree of v. Assume that the tree G does not contain vertices of degree 2. The vertices of degree 1 are called boundary vertices. Denote the set of boundary vertices of the graph G by G. For the sake of convenience, let each boundary vertex v i be an end of the edge e i, such edges are called boundary edges. All other vertices and edges are called internal. Let the vertex v r G be the root of the tree. Each edge e j E is viewed as a segment [0,T j ] and is parametrized by the parameter x j [0,T j ]. The value x j = 0 correspond to one of the end vertices of the edge e j, and x j = T j corresponds to another one. For a boundary edge, the end x j = 0 corresponds to the boundary vertex v j. A function on the tree G can be represented as a vector function y = [y j ] m j=1, where y j = y j (x j ), x j [0,T j ], j = 1,m. Let e j = [v i,v k ], i.e. the vertex v i corresponds to the end x j = 0 and the vertex v k corresponds to x j = T j. Introduce the following notation y j (v i ) = y j (0), y j (v k ) = y j (T j ), y j(v i ) = y j(0), y j(v k ) = y j(t j ). If v i G, we omit the index of the edge and write y(v i ) and y (v i ). Consider the Sturm-Liouville equation on G: y j +q j(x j )y j = λy j, x j [0,T j ], j = 1,m 1. (1) where λ is the spectral parameter, q j L[0,T j ]. We call the function q = [q j ] m j=1 the potential on the graph G. The functions y j, y j are absolutely continuous on the segments [0,T j ] and satisfy the standard matching conditions in the internal vertices v V\ G: { y j (v) = y k (v), e j,e k E v (continuity condition), e j E v y j (v) = 0, (Kirchhoff s condition). (2) Let L 0 and L k, v k G, be the boundary value problem for system (1) with the matching conditions (2) and the following conditions in the boundary vertices: L 0 : y(v i ) = 0, v i G, (3) L k : y (v k ) = 0, y(v i ) = 0, v i G\{v k }. (4) It is well-known, that the problems L k have discrete spectra, which are the countable sets of eigenvalues Λ k = {λ ks } s=1, k = 0 or v k G. Fix a boundary vertex v k G. Let Ψ k = [ψ kj ] m j=1, ψ kj = ψ kj (x j,λ), be the solution of the system (1), satisfying the matching conditions (2) and the boundary conditions ψ kk (0,λ) = 1, ψ kj (0,λ) = 0, v j G\{v k }. Denote M k (λ) = ψ kk (0,λ). The functions Ψ k and M k are called the Weyl solution and the Weyl function of (1) with respect to the boundary vertex v k, respectively. The notion of the Weyl function for the tree generalizes the notion of the Weyl function (m-function) for the classical Sturm-Liouville operator on a finite interval [10, 13]. If the tree G consists of only one edge, then M k (λ) coincide with the classical Weyl function. Consider the following inverse problems. 2

Inverse Problem 1. Given the spectra Λ 0, Λ k, v k G\{v r }, construct the potential q on the tree G. Inverse Problem 2. Given the Weyl functions M k (λ), v k G\{v r }, construct the potential q on the tree G. Note that if the number of boundary vertices is b, then one needs b spectra or b 1 Weyl functions to recover the potential. We do not require the data associated with the root v r. There is a close relation between Inverse Problems 1 and 2. The Weyl functions can be represented in the form M k (λ) = k(λ) 0 (λ), v k G, (5) where k (λ) are characteristic functions of the boundary value problems L k. If the eigenvalues Λ k are known, one can construct characteristic functions as infinite products by Hadamard s theorem. Thus, with the system of spectra, one can obtain the Weyl functions and reduce Inverse Problem 1 to Inverse Problem 2. V.A. Yurko has proved, that Inverse Problems 1 and 2 are uniquely solvable, and provided a constructive algorithm for the solution by the method of spectral mappings [13]. In the remaining of this section, we shall briefly describe his algorithm. Let the Weyl functions M k (λ), v k G\{v r } be given. Consider the following auxiliary problem. Problem IP(k). Given M k (λ), construct the potential q k (x k ) on the edge e k. Note that this problem is not equivalent to the inverse problem on the finite interval, since the Weyl function M k (λ) contains information from the whole graph. However, it can be solved uniquely by the method of spectral mappings, and the potential on the boundary edges can be recovered. Then V.A. Yurko used so-called µ-procedure to recover the potential on the internal edges. We reformulate these ideas in the form, which is more convenient for us in the future. Theorem 1. Let v be an internal vertex, connected with the set of boundary vertices V G\{v r } and only one other vertex. Suppose the potentials q k on the edges e k are known for all v k V, as well as a Weyl function M k (λ) for at least one vertex from the set V. Denote G the graph by removing the vertices v k V together with the corresponding edges e k from the graph G. Then the Weyl function for the graph G with respect the the vertex v can be determined from the given information. Applying Theorem 1, one can cut the boundary edges off, until the potential will be recovered on the whole graph. 3. Partial inverse problems In this section, the main results of the paper are formulated. We assume that the potential is known on one edge of the tree and formulate partial inverse problems. We consider separately the cases of boundary and internal edge. The first one appears to be trivial, for the second one we describe the procedure of the constructive solution. For the convenience of the reader, the proofs of the technical lemmas are provided in Section 4. Inverse Problem 3. Let e f be a boundary edge (f r). Given the potential q f on the edge e f and the spectra Λ 0, Λ k, v k G\{v f,v r }. Construct the potential q on the tree G. The solution of Inverse Problem 3 is a slight modification of the method described in Section 2. From Λ 0, Λ k, v k G\{v f,v r }, we easily construct the potentials q k for v k G\{v f,v r }. The potential q f is known, so we can apply Theorem 1 iteratively and recover the potential on G. 3

Now let e f be an internal edge. If this edge is removed, the graph splits into two parts, call them P 1 and P 2. Let P 1 and P 2 be the sets of boundary vertices of P 1 and P 2, respectively. Fix two arbitrary vertices v r1 P 1 and v r2 P 2. Inverse Problem 4. Given the potential q f on the internal edge e f, the spectra Λ 0, Λ k v k G\{v r1,v r2 }. Construct the potential q on the tree G. Solution of Inverse Problem 4. For simplicity, we assume that the ends of the edge e f have degree 3. The general case requires minor modifications. If one splits each of the ends of e f into three vertices, the tree splits into five subtrees G i, i = 1,5, such that v r1 G 2, v r2 G 5, and G 3 contains the only edge e f (see fig. 1). Let v 1 and v 4 are arbitrary boundary vertices of the trees G 1 and G 4 (different from the ends of e f ), v r1 = v 2, v r2 = v 5, e f = [v 3,v 6 ]. G v 1 G 2 1 v2 v 3 G 3 v 6 v 4 G 4 G 5 v 5 Figure 1 Step 1. Construct the characteristic functions k (λ) by the given spectra Λ k, k = 0 and v k G\{v 2,v 5 }. Find M k (λ) by formula (5). Step 2. Consider trees G 1 and G 4. Recover the potential q on the edges of G 1 and G 4, using the solutions of the problems IP(k) for v k G 1 \{v 3 } and v k G 4 \{v 6 }, and them applying Theorem 1 iteratively. Step 3. Introduce the characteristic functions of the boundary value problems for the Sturm-Liouville equations (1) on the graphs G 1 -G 5 with the standard matching conditions (2) in internal vertices and the following conditions in the boundary vertices: graphg 1 DD 1 (λ): y(v k ) = 0, v k G 1, ND DN NN 1 (λ): y (v 1 ) = 0, y(v k ) = 0, v k G 1 \{v 1 }, 1 (λ): y (v 3 ) = 0, y(v k ) = 0, v k G 1 \{v 3 }, 1 (λ): y (v 1 ) = 0, y (v 3 ) = 0, y(v k ) = 0, v k G 1 \{v 1,v 3 }. { D graphg 2 (λ): y(v k ) = 0, v k G 2, 2 N 2 (λ): y (v 3 ) = 0, y(v k ) = 0, v k G 2 \{v 3 }. DD 3 (λ): y(v 3 ) = 0, y(v 6 ) = 0, graphg ND 3 (λ): y (v 3 ) = 0, y(v 6 ) = 0, 3 DN 3 (λ): y(v 3 ) = 0, y (v 6 ) = 0, NN 3 (λ): y (v 3 ) = 0, y (v 6 ) = 0. 4

DD 4 (λ): y(v k ) = 0, v k G 4, graphg ND 4 (λ): y (v 4 ) = 0, y(v k ) = 0, v k G 4 \{v 4 }, 4 DN 4 (λ): y (v 6 ) = 0, y(v k ) = 0, v k G 4 \{v 6 }, NN 4 (λ): y (v 4 ) = 0, y (v 3 ) = 0, y(v k ) = 0, v k G 4 \{v 4,v 6 }. { D graphg 5 (λ): y(v k ) = 0, v k G 5, 5 N 5 (λ): y (v 6 ) = 0, y(v k ) = 0, v k G 5 \{v 6 }. Lemma 1. The following relation holds DD 1 (λ) D 2 (λ) 0 0 0 0 0 D 2 (λ) 1 0 0 0 0 (λ) = DN 1 (λ) N 2 (λ) 0 1 0 0 0 0 ND 3 (λ) DD 3 (λ) DD 4 (λ) 0. (6) 0 0 0 0 DD 4 (λ) D 5 (λ) 0 0 NN 3 (λ) DN 3 (λ) DN 4 (λ) N 5 (λ) If one changes DD 1 (λ) to ND 1 (λ) and DN 1 (λ) to NN 1 (λ), he obtains the determinant equal to 1 (λ). Similarly, if one changes DD 4 (λ) to ND 4 (λ) and DN 4 (λ) to NN 4 (λ), he gets 4 (λ). Step 4. Note that the functions 0 (λ), 1 (λ), 4 (λ) are known from Step 1. Since we know the potential on the graphs G 1, G 4 (from Step 2) and G 3 (given a priori), we can easily construct the characteristic functions for these graphs. Consider the relation (6) and similar relations for 1 (λ) and 4 (λ) as a system of equations with respect to D 2 (λ), N 2 (λ), D 5 (λ) and N 5 (λ) in the following form a 11 D 2 D 5 +a 12 N 2 D 5 +a 13 D 2 N 5 +a 14 N 2 N 5 = 0, a 21 D 2 D 5 +a 22 N 2 D 5 +a 23 D 2 N 5 +a 24 N 2 N 5 = 1, (7) a 31 D 2 D 5 +a 32 N 2 D 5 +a 33 D 2 N 5 +a 34 N 2 N 5 = 4, where a ij = a ij (λ), i = 1,3, j = 1,4, are known coefficients. Step 5. Multiply the first equation of (7) by 1 and subtract the second equations, multiplyed by 0. Apply the similar trasform to the first and the third equations. Then we obtain the system { b 11 D 2 D 5 +b 12 N 2 D 5 +b 13 D 2 N 5 +b 14 N 2 N 5 = 0, where Divide both equations by D 2 D 5. where b 21 D 2 D 5 +b 22 N 2 D 5 +b 23 D 2 N 5 +b 24 N 2 N 5 = 0, b 1i = a 1i 1 a 2i 0, b 2i = a 1i 4 a 3i 0, i = 1,4. (8) b i1 +b i2 M2 +b i3 M5 +b i4 M2 M5 = 0, i = 1,2, (9) M 2 (λ) = N 2 (λ) D 2 (λ), M5 (λ) = N 5 (λ) D 5 (λ) are (up to the sign) the Weyl functions for the subtrees G 2 and G 5 associated with the vertices v 3 and v 6, respectively. Step 6. From the system (9) we easily derive M 5 = b i1 +b i2 M2 b i3 +b i4 M2, i = 1,2. 5

Hence (b 11 +b 12 M2 )(b 23 +b 24 M2 ) = (b 21 +b 22 M2 )(b 13 +b 14 M2 ). Finally, we obtain the quadratic equation with respect to M 2 (λ): with analytic coefficients A(λ), B(λ), C(λ): A(λ) M 2 2 (λ)+b(λ) M 2 (λ)+c(λ) = 0, (10) A = b 12 b 24 b 22 b 14, B = b 11 b 24 +b 12 b 23 b 21 b 14 b 22 b 13, C = b 11 b 23 b 21 b 13. (11) Step 7. Consider the Sturm-Liouville equation (1) on the tree G with the potential q = 0. Implement Steps 1 6 for this case and obtain the quadratic equation A 0 (λ) M 2 20(λ)+B 0 (λ) M 20 (λ)+c 0 (λ) = 0, (12) analogous to (10). Denote ρ = λ, Reρ 0, S δ := {ρ: Reρ 0, Imρ δ}, δ > 0, [1] = 1+O(ρ 1 ). Let f(ρ 2 ) be an analytic function and ε > 0. Denote Z ε (f) := {ρ: f(ρ 2 ) εn 0}. Lemma 2. The following asymptotic relations hold A(λ) = A 0 (λ)[1], B(λ) = B 0 (λ)[1], C(λ) = C 0 (λ)[1], ρ S δ Z ε (A 0 B 0 C 0 ), ρ. Consequently, D(λ) = D 0 (λ)[1] for ρ S δ Z ε (D 0 ), ρ, where D(λ) and D 0 (λ) are discriminants of equations (10) and (12), respectively. Lemma 3. A 0 (λ) 0, D 0 (λ) 0. It follows from Lemmas 2 and 3, that the quadratic equation (10) does not degenerate for ρ S δ Z ε (A 0 D 0 ), and two roots of (10) are different by asymptotics as ρ. One can easily find an asymptotic representation of M2 (λ) for any particular graph and choose the correct root of (10) on some region of S δ for sufficiently large ρ. Then the function M 2 (λ) can be constructed for all λ C except its singularities by analytic continuation. Similarly one can find M 5 (λ). Step 8. Consider the tree G 2 with the root v 2. Solve problems IP(k) by M k (λ), v k G 2 \{v 2,v 6 }, and by M 2 (λ) for v 3, obtain the potential on the boundary edges except e 2. Then apply the cutting of boundary edges by Theorem 1 and recover the potential q on G 2. The subtree G 5 can be treated similarly. Thus, we recovered the potential q on the whole graph G. In parallel, we have proved the following uniqueness theorem. Theorem 2. Let the potential q f on the edge e f (f r) be known. (i) If e f in a boundary edge, the spectra Λ 0, Λ k, v k G\{v f,v r }, uniquely determine the potential q on the whole graph G. (ii) If e f is an internal edge, the spectra Λ 0, Λ k v k G\{v r1,v r2 } uniquely determine the potential q on the whole graph G. Using the described methods with some technical modifications, one can solve partial inverse problems by Weyl functions. Inverse Problem 5. Let e f be a boundary edge (f r). Given the potential q f on the edge e f and the Weyl functions M k (λ), v k G\{v f,v r }. Construct the potential q on the tree G. 6

Inverse Problem 6. Given the potential q f on the internal edge e f, the Weyl functions M k (λ) v k G\{v r1,v r2 }. Construct the potential q on the tree G. Thus, ifthenumber ofboundaryedgesisbandthepotentialisknownononeedge(boundary or internal), b 2 Weyl functions are required to construct q on the whole graph. 4. Proofs 4.1. Proof of Lemma 1. Consider the Sturm-Liouville equation (1) on the tree G. Let C j (x j,λ) and S j (x j,λ) be solutions of (1) on the edge e j under initial conditions C j (0,λ) = S j(0,λ) = 1, C j(0,λ) = S j (0,λ) = 0. Any solution y = [y j ] m j=1 of the equation (1) on G admits the following representation y j (x j,λ) = M 0 j (λ)c j(x j,λ)+m 1 j (λ)s j(x j,λ), j = 1,m, x j [0,T j ]. (13) Let BC be some fixed boundary conditions in the vertices v G of the form y(v) = 0 or y (v) = 0 (for instance, we consider conditions (3) for the problem L and (4) for the problem L k ). Denote by L the boundary value problem for the Sturm-Liouville equation (1) with the standard matching conditions (2) and the boundary conditions BC. If y is a solution of a boundary value problem L, substitute (13) into (2) and BC, and obtain a linear algebraic system with respect to Mj 0(λ), M1 j (λ). It is easy to check that the determinant of this system is a characteristic function (λ) of the boundary value problem L, i.e. zeros of (λ) coincide with the eigenvalues of L. Example 1. Consider the problem L 0 for the star-type graph for m = 3. Then boundary conditions (3) yield M1 0(λ) = M0 2 (λ) = M0 3 (λ) = 0. Consequently, from (2) we obtain the system with respect to Mj 0 (λ), j = 1,2,3, with the determinant S 1 (T 1,λ) S 2 (T 2,λ) 0 0 (λ) = 0 S 2 (T 2,λ) S 3 (T 3,λ) S 1(T 1,λ) S 2(T 2,λ) S 3(T 3,λ). In the general case, the following assertion is valid. Lemma 4. Let w V and the degree of w be equal n. Splitting the vertex w, we split G into n subtrees G i, i = 1,n. For each i = 1,n, let D i (λ) and N i (λ) be characteristic functions for boundary value problems for equation (1) on tree G i with matching conditions (2), boundary conditions BC for v G G i and the Dirichlet condition y(u) = 0 for D i (λ) and the Neumann condition y (u) = 0 for N i (λ). Then the characteristic function (λ) for G with the conditions (2) and BC admits the following representation: D 1 (λ) D 2 (λ) 0... 0 0 D 2 (λ) D 3 (λ)... 0 (λ) =.... (14) 0 0 0 D n(λ) N 1 (λ) N 2 (λ) N 3 (λ)... N n (λ) Indeed, if we write the determinant for (λ) and analyze the participation of the edges of G i in this determinant, we can easily see that (λ) = D i (λ)d i(λ) + N i (λ)e i(λ), where the functions D i (λ) and E i (λ) do not depend on the subtree G i. Thus we can consider the simplest case of the star-type graph, when each G i contains only one edge, and then change 7

the multipliers, corresponding to subgraphs G i, to D i (λ) and N i (λ). Thus we directly obtain (14) from the formula for the star-type graph. Lemma 1 follows from Lemma 14 for the graph in the fig. 1. Alternatively, one can derive (6) from (21), changing characteristic functions for one-edge subtrees by general characteristic function. 4.2. Proof of Lemma 2. Together with L consider the boundary value problem L 0 for equation (1) with q 0, the matching conditions (2) and the boundary conditions BC. If some symbol γ denotes the object related to L, we denote by the symbol γ 0 the similar object related to L 0. In particular, 0 (λ) is the characteristic function of L 0. Let the symbol P(ρ) stand for different polynomials of sinρt j and cosρt j, j = 1,m. Lemma 5. The characteristic function (λ) has the following asymptotic behavior: (λ) = 0 (λ)+o(ρ d ) = P(ρ) ρ d 1 +O(ρ d ), ρ S δ, ρ, where P(ρ) 0 and d = m i n, where m is the number of the edges, i is the number of internal vertices and n is the number of boundary vertices with the Neumann boundary condition y (v) = 0. Proof. The claim of the lemma immediately follows from the standard asymptotic formulas C j (x j,λ) = cosρx j +O(ρ 1 ), C j (x j,λ) = ρsinρx j +O(1), S j (x j,λ) = sinρx j ρ +O(ρ 2 ), S j (x,λ) = cosρx j, ρ S δ, ρ, and the construction of (λ). The relation P(ρ) 0 follows from the regularity of the standard matching conditions. Applying Lemma 5 to the characteristic functions, defined on Step 3 of the algorithm, we derive asymptotic representations for the coefficients c = a ij,b ij,a,b,c in the following form: c(λ) = c 0 (λ)+o(ρ d ) = P(ρ) ρ d 1 +O(ρ d ), ρ S δ, ρ, where d stands for different integers. This relation yields Lemma 2. 4.3. Proof of Lemma 3. In this subsection, we consider only the problem L 0 with q 0, so we omit the index 0 for brevity. For simplicity, let T f = 1. Taking into account, that DD 3 = sinρ ρ, ND 3 = DN 3 = cosρ, NN 3 = ρsinρ and doing some algebra with the expressions (6), (8), (11), we derive A(λ) = F 1 (λ)f 4 (λ) 0 (λ) sin2 ρ DD ρ 2 4 (λ) D 5 (λ)χ(λ), (15) B(λ) = F 1 (λ)f 4 (λ) sinρ { } ρ 0(λ) D 5 (λ)π(λ)+ D 5 (λ)sinρ ρ ξ(λ) DD 4 (λ) N 5 (λ)sinρ ρ χ(λ), (16) C(λ) = F 1 (λ)f 4 (λ) 0 (λ) sinρ { ρ N 5 (λ) Π(λ)+ sinρ } ρ ξ(λ), (17) 8

where F i (λ) = DD i (λ) NN i (λ) DN i (λ) ND i (λ), i = 1,4, Π(λ) = 2 DD 1 (λ) DD 2 (λ) DD 4 (λ), χ(λ) and ξ(λ) are characteristic functions of the graphs G 1 G 2 G 3 and G 1 G 2 G 3 G 4, respectively. Here we mean that the copies of the vertex v 3 (and v 6 in the second graph) are joined into one vertex with the standard matching conditions (2). Lemma 6. Let v 1 and v 2 be two fixed vertices from G. Denote by DD (λ), DN (λ), ND (λ) and NN (λ) the characteristic functions for equation (1) on the tree G with the matching conditions (2), with the following boundary conditions: DD (λ): y(v 1 ) = y(v 2 ) = 0, DN (λ): y(v 1 ) = y (v 2 ) = 0, ND (λ): y (v 1 ) = y(v 2 ) = 0, NN (λ): y (v 1 ) = y (v 2 ) = 0, and with the conditions BC in the vertices v G\{v 1,v 2 }. Then DD (λ) NN (λ) DN (λ) ND (λ) 0. (18) Proof. We shall divide the proof into the following steps. 1. Let the tree G consists of the only edge [v 1,v 2 ]. Then one can check the relation (18) by direct calculation. 2. Let the vertices v 1 and v 2 be connected by edges with the same vertex v, and let there also be subtrees G i, i = 1,n, from the vertex v (see fig. 2). Denote by D i (λ) and N i (λ) the characteristic functions for G i with the matching conditions (2), the boundary conditions BC and y(v) = 0 for D i (λ) and y (v) = 0 for N i (λ). According to Lemma 4, the following relation holds where DD (λ) = sinρt 1sinρT 2 K (λ)+ 1 ρ 2 ρ (sinρt 1cosρT 2 +cosρt 1 sinρt 2 ) Π (λ), Π (λ) = n D i (λ), K (λ) = Π (λ) i=1 n i=1 N i (λ) D i (λ). Using similar representations for NN (λ), DN (λ) and ND (λ), we derive DD (λ) NN (λ) DN (λ) ND (λ) = ( Π (λ) ) 2 0. 3. Now let the vertices v 1 and v 2 be connected by the edges with v 3 and v 4, respectively. Let the tree G splits by the vertices v 3 and v 4 into the subtrees G i, i = 1,n 1, connected with v 3, the subtrees G j, j = 1,n 2, connected with v 4, the subtree G 0, including the both vertices v 3 and v 4, and the edges e 1, e 2 (see fig. 3). Denote by D i (λ), N i (λ), i = 1,n 1, and by D j (λ), N j (λ), j = 1,n 2, the characteristic functions for the subtrees G i with the Dirichlet or Neumann boundary condition in v 3 and for the subtrees G i with the Dirichlet or Neumann boundary condition in v 4, respectively. Let DD 0 (λ), DN 0 (λ), ND 0 (λ) and NN 0 (λ) be characteristic functions for the subtree G 0 with the following boundary conditions DD DN ND NN 0 (λ): y(v 3 ) = y(v 4 ) = 0, 0 (λ): y(v 3 ) = y (v 4 ) = 0, 0 (λ): y (v 3 ) = y(v 4 ) = 0, 0 (λ): y (v 3 ) = y (v 4 ) = 0, 9

G 1 v 2 v 1 e e 2 1 v G 2 G 3 Figure 2 v 1 G 1 G 2 e 1 v 3 G 0 v 4 G 1 e 2 v 2 Figure 3 G2 and the conditions BC in other boundary vertices. Denote the functions n 1 Π 1 (λ) = i=1 n 2 D i (λ), Π 2 (λ) = j=1 D j (λ), n K 1 (λ) = Π 1 (λ) N i (λ) n N D i=1 i (λ), K 2 (λ) = Π 2 (λ) j (λ) D j=1 j (λ). KK (λ) = DD 0 (λ) K 1 (λ) K 2 (λ)+ ND 0 (λ) Π 1 (λ) K 2 (λ) + DN 0 (λ) K 1 (λ) Π 2(λ)+ NN 0 (λ) Π 1(λ) Π 2(λ), ΠK (λ) = DD 0 (λ) Π 1 (λ) K 2 (λ)+ DN 0 (λ) Π 1 (λ) Π 2 (λ), KΠ (λ) = DD 0 (λ) K 1 (λ) Π 2 (λ)+ ND 0 (λ) Π 1 (λ) Π 2 (λ), ΠΠ (λ) = DD 0 (λ) Π 1(λ) Π 2(λ). In view of Lemma 4, the following relation holds (19) DD (λ) = sinρt 1sinρT 2 KK (λ)+ cosρt 1sinρT 2 ΠK (λ) ρ 2 ρ 10

+ sinρt 1cosρT 2 KΠ (λ)+cosρt 1 cosρt 2 ΠΠ (λ). ρ Together with the similar relations for DN (λ), ND (λ) and NN (λ), it yields DD (λ) NN (λ) DN (λ) ND (λ) = ΠΠ (λ) KK (λ) ΠK (λ) KΠ (λ) Taking (19) into account, we obtain ΠΠ (λ) KK (λ) ΠK (λ) KΠ (λ) = ( DD 0 (λ) NN 0 (λ) DN 0 (λ) ND 0 (λ)) ( Π 1 (λ) Π 2 (λ)) 2. By virtue of Lemma 5, Π i (λ) 0, i = 1,2. Therefore the relation (18) holds for the tree G if and only if it holds for the subtree G 0. By induction, the claim of the lemma is valid for any tree G. By virtue of Lemmas 5, 6 and (15), A(λ) 0. It follows from (15), (16), (17), that D(λ) = B 2 (λ) 4A(λ)C(λ) = F 2 1 (λ)f2 4 (λ)sin2 ρ ρ 2 0 (λ) Note that the expression in the bracket above equals to In view of Lemma 5, D 5 (λ)π(λ) = C 1r p exp(r(t 1))[1], { } 2 D 5 (λ)π(λ)+ D 5 (λ)sinρ ρ ξ(λ)+ DD 4 (λ) N 5 (λ)sinρ ρ χ(λ). D 5 (λ)π(λ)+ sinρ ρ 0(λ). sinρ ρ 0(λ) = C 2 r q exp(r(t +1))[1], where ρ = ir, r +, T = m T j, C 1, C 2, p and q are some constants. Clearly, the second j=1 term grows faster than the first one. Therefore 0 (λ) 0 implies D(λ) 0. The proof of Lemma 3 is finished. Using Lemma 5, one can also check, that B(λ) and D(λ) have the same power of ρ in the denominator, so the roots of (10) have different asymptotic behavior. 5. Example In this section, we provide the solution of Inverse Problem 4 for the example of the graph in the fig. 4. For simplicity, let T j = 1, j = 1,5. Let x 3 = 0 corresponds to the vertex v 3 and x 3 = 1 corresponds to v 6. For the boundary edges, x j = 0 correspond to the boundary vertices. The matching conditions (2) take the form v 3 : y 1 (1) = y 2 (1) = y 3 (0), y 1 (1)+y 2 (1) y 3 (0) = 0, v 6 : y 3 (1) = y 4 (1) = y 5 (1), y 3 (1)+y 4 (1)+y 5 (1) = 0. (20) For this example, each subtree G i consists of only one edge e i, i = 1,5. Let us know the spectra Λ, Λ 1, Λ 4 and the potential q 3. Using the given spectra, one can easily find the characteristic functions (λ), 1 (λ), 4 (λ) and the Weyl functions M 1 (λ), M 4 (λ). Solving problems IP(1) and IP(4), recover q 1 and q 3. 11

v 1 v 2 e 1 e 2 v 3 e 3 v 6 e 4 e 5 v 4 v 5 Figure 4. Example Consider the boundary value problem L. Represent the solution y in the form (13) and substitute it into (2) and (3). From (3), one gets M1 0(λ) = M0 2 (λ) = M0 4 (λ) = M0 5 (λ) = 0. Then matching conditions (20) yield the system S 1 S 2 0 0 0 0 M1 1 0 S 2 1 0 0 0 M 1 S 1 S 2 2 0 1 0 0 M 0 3 0 0 C 3 S 3 S 4 0 M 1 = 0. (21) 3 0 0 0 0 S 4 S 5 M 1 0 0 C 3 S 3 S 4 S 5 4 M5 1 Here we omit arguments (1, λ) and (λ) for brevity. The characteristic function (λ) equals the determinant of (21). Since we know q 1, q 3 and q 4, we can solve (1) and obtain the functions S j (x j,λ) and C j (x j,λ) for j = 1,3,4. Therefore the determinant admits the following representation = a 11 S 2 S 5 +a 12 S 2 S 5 +a 13 S 2 S 5 +a 14S 2 S 5, where a 11 = S 1 S 3 S 4 S 3 S 4 +S 1 C 3 S 4 C 3 S 4, a S 12 = S 3 S 4 1 S 3 S 4, a 13 = (S 1S 3 +S 1 C 3 )S 4, a 14 = S 1 S 3 S 4. If one change S 1 to C 1 or S 4 to C 4, he obtains analogous relations for 1 (λ) and 4 (λ), respectively. Thus we arrive at the system (7). Let q 0 on G. Then a 0 11 = sin3ρ ρ Cj 0 (x j,λ) = cosρx j, Sj 0 (x j,λ) = sinρx j, ρ, a 0 12 = a0 13 = sin2ρsinρ, a 0 ρ 2 14 = sin3 ρ, ρ 3 a 0 21 = a0 31 = cos3ρ, a0 22 = a0 33 = sin2ρcosρ, ρ a 0 23 = a0 32 = cos2ρsinρ, a 0 24 ρ = a0 34 = cosρsin2 ρ. ρ 2 12

0 = 9sin5ρ+13sin3ρ+6sinρ, 0 16ρ 3 1 = 0 4 = 9cos5ρ+7cos3ρ+2cosρ. 16ρ 2 Using (8), we obtain b 0 11 = b 0 21 = 3sin6ρ 2sin4ρ+13sin2ρ, b 0 16ρ 3 12 = b 0 23 = 3cos6ρ+6cos4ρ+3cos2ρ 6, 16ρ 4 b 0 13 = b0 22 = 3cos6ρ 10cos4ρ+13cos2ρ 6, b 0 32ρ 4 14 = b0 24 = 3sin6ρ+12sin4ρ 15sin2ρ. 32ρ 5 Substitute these formulas into (11): A 0 = 27sin12ρ+174sin10ρ 420sin8ρ+378sin6ρ+153sin4ρ 468sin2ρ 2048ρ 9, B 0 = 27cos12ρ+84cos10ρ+106cos8ρ 764cos6ρ+1099cos4ρ 344cos2ρ 154 2048ρ 8, C 0 = 27sin12ρ+48sin10ρ+140sin8ρ 336sin6ρ 71sin4ρ+512sin2ρ 1024ρ 7. Calculate the discriminant of equation (12): D 0 = B 2 0 4A 0 C 0 = (6561cos24ρ 52488cos22ρ+128628cos20ρ+83592cos18ρ 987134cos16ρ+1543976cos14ρ+702372cos12ρ 4646312cos10ρ+3755087cos8ρ +3053616cos6ρ 4805144cos4ρ 4176688cos2ρ+5393934)/(8388608ρ 16 ). We used wxmaxima 12.04.0 for calculations. Obviously, A 0 (λ) 0, D 0 (λ) 0, so according to Lemma 2, the roots of equation (10) in the general case have different asymptotics: M 2 1 ρcosρ (λ) = sinρ [1], M2 2 (λ) = 1+6cos2 ρ 3sinρcosρ [1]. Since M 2 (λ) = S 2 (1,λ), only the root S 2 (1,λ) M1 2 (λ) is the required one. Finally, one can easily find M 5 (λ) and solve classical Sturm-Liouville inverse problems by Weyl functions on the edges e 2 and e 5. Acknowledgments. This work was supported by Grant 1.1436.2014K of the Russian Ministry of Education and Science, by Grants 13-01-00134, 14-01-31042 and 15-01-04864 of Russian Foundation for Basic Research and by the Mathematics Research Promotion Center of Taiwan. References [1] Langese, J.; Leugering, G.; Schmidt, J. Modelling, analysis and control of dynamic elastic multi-link structures. Birkhäuser, Boston (1994). [2] Kottos, T; Smilansky, U. Quantum chaos on graphs. Phys. Rev. Lett. 79(1997), 4794 4797. [3] Kuchment, P. Quantum graphs. Some basic structures. Waves Random Media 14 (2004), S107 S128. [4] Belishev, M.I. Boundary spectral inverse problem on a class of graphs (trees) by the BC method. Inverse Problems 20 (2004), 647 672. 13

[5] Pokornyi, Yu.; Borovskikh, A. Differential equations on networks (geometric graphs). J. Math. Sci. (N.Y.) 119:6 (2004), 691 718. [6] Analysis on Graphs and Its Applications, edited by P. Exner, J.P. Keating, P. Kuchment, T. Sunada and Teplyaev, A. Proceedings of Symposia in Pure Mathematics, AMS, 77 (2008). [7] Yurko, V.A. Inverse spectral problems for differential operators on arbitrary compact graphs. J. Inverse and Ill-Posed Probl. 18:3 (2010), 245 261. [8] Yurko, V. Inverse spectral problems for Sturm-Liouville operators on graphs. Inverse Problems 21 (2005), 1075 1086. [9] Yurko, V.A. On recovering Sturm-Liouville operators on graphs, Mathematical Notes, 79:4 (2006), 572-582. [10] Marchenko, V. A. Sturm-Liouville Operators and their Applications, Naukova Dumka, Kiev (1977) (Russian); English transl., Birkhauser (1986). [11] Levitan, B.M. Inverse Sturm-Liouville Problems, Nauka, Moscow(1984)(Russian); English transl., VNU Sci. Press, Utrecht (1987). [12] Pöschel, J; Trubowitz, E. Inverse Spectral Theory, New York, Academic Press (1987). [13] Freiling, G.; Yurko, V. Inverse Sturm-Liouville problems and their applications. Huntington, NY: Nova Science Publishers, 305 p. (2001). [14] Yurko, V.A. Method of Spectral Mappings in the Inverse Problem Theory. Inverse and Ill-Posed Problems Series, Utrecht: VSP (2002). [15] Yang, C.-F. Inverse spectral problems for the Sturm-Liouville operator on a d-star graph, J. Math. Anal. Appl. 365 (2010), 742 749. [16] Hochstadt, H.; Lieberman, B. An inverse SturmLiouville problem with mixed given data, SIAM J. Appl. Math. 34 (1978), 676-680. Natalia Bondarenko Department of Mechanics and Mathematics Saratov State University Astrakhanskaya 83, Saratov 410012, Russia bondarenkonp@info.sgu.ru Chung-Tsun Shieh Department of Mathematics Tamkang University 151 Ying-chuan Road Tamsui, Taipei County, Taiwan 25137, R.O.C. ctshieh@mail.tku.edu.tw 14