Boundary estimates for solutions of non-homogeneous boundary value problems on graphs SONJA CURRIE University of the Witwatersrand School of Mathematics Private bag 3 P.O. Wits 2050 South Africa [email protected]

BRUCE A. WATSON University of the Witwatersrand School of Mathematics Private bag 3 P.O. Wits 2050 South Africa [email protected]

Abstract: We consider the Sturm-Liouville operator on a graph and give bounds for the norms of the boundary values of solutions to the non-homogeneous boundary value problem in terms of the norm of the non-homogeneity. In addition the eigenparameter dependence of these bounds is studied. Key–Words: Differential Operators on Graphs, Boundary Estimates

1 Introduction We consider the Sturm-Liouville equation ly := −

d2 y + q(x)y = λy, dx2

(1)

where q is real-valued and essentially bounded (with respect to Lebesgue measure), on a weighted graph G with formally self-adjoint boundary conditions at the nodes. For characterisations of self-adjoint boundary value problems on graphs and associated boundary conditions we refer the reader to [5] and [11]. In [10], it was shown that the geometry of a noncommensurate simple graph is uniquely dependent on the spectrum of the Laplacian on the graph. I.e. for zero potential they reconstructed the boundary conditions (of a specific type) from a single spectrum. In [6] spectral asymptotics were given for l on graphs where all edges are of equal length while in [7] and [8] eigenvalue asymptotics were given for l on general compact graphs via matrix Pr¨ufer angle techniques and Dirichlet-Neumann bracketing respectively. Variational aspects of boundary value problems on graphs were studied in [2], [8] and [20], and on trees in [19]. Sturmian oscillation theory was extended to Sturm-Liouville operators on graphs by Pokornyi and Pryadiev, and Pokornyi, Pryadiev and Al-Obeid, in [15] and [16]. Sturm-Liouville problems on finite graphs arise naturally in quantum mechanics and circuit theory, [3, 9]. In [13] and the bibliography thereof one can find an extensive collection of physical systems that give rise to Sturm-Liouville problems on graphs.

Here we consider solutions of non-homogeneous Sturm-Liouville problems on graphs and their a priori estimates. Particular attention is given to the relationship between the norm of the solutions to the non-homogeneous boundary value problem on the boundary of the graph and the norm of the nonhomogeneous term on the graph, see Theorem 3.2. In addition the eigenparameter dependence of this relationship is explored. To complete the paper an example is provided in Section 4, illustrating Theorem 3.2. The results obtained in this paper rely on an ability to make the transition between local results (on each edge) and global results (on the whole graph). Thus the method employed involves two main steps: establishing a local result on each edge; and the transition from the local results to a global result on the whole graph. It should be noted that for Sturmian systems it is only possible to find boundary estimates, of the form given in this paper, for two specific types of systems. Firstly for non-weighted systems with general, essentially bounded, Hermitian potential (not necessarily diagonal) and secondly for diagonal systems, which are equivalent to considering Sturm-Liouville equations on graphs as is done in this paper, see [7] for the equivalence. In [17, 18], Schauder considers interior estimates and estimates near the boundary for solutions of second order elliptic boundary value problems. His estimates near the boundary are for solutions of the Dirichlet problem. Estimates near the boundary for other than Dirichlet boundary conditions have been obtained by Miranda, [14], for second order elliptic

boundary value problems and by Agmon, Douglis, Nirenberg and Browder, [1, 4], for arbitrary order elliptic operators. In the above references it should be noted that the estimates are given in a region near the boundary whereas our results provide estimates on the boundary.

2 Preliminaries Let G denote a directed graph with a finite number of nodes and edges, with each edge parametrized by path-length and having finite length. Each edge, ei , of length say li can thus be considered as the interval [0, li ]. Having made this identification, it is possible to consider the differential equation (1) on the graph G to be the collection of differential equations d2 yi (2) − 2 + qi (x)yi = λyi , x ∈ [0, li ], dx for i = 1, ..., K, where qi and yi denote the restriction of q and y to ei , respectively. It is now possible, at each node, ν, to specify boundary conditions in terms of the values of y and y ′ at ν on each of the incident edges. In particular if the edges which originate at node ν are ei , i ∈ Λs (ν) and the edges which terminate at node ν are ei , i ∈ Λe (ν) then the boundary conditions at ν are of the form j∈Λs (ν)

h

X

h

i

i



γij yj (lj ) + δij y j (lj ) = 0,

j∈Λe (ν)

(3)

Remark It should be noted that by setting αij = 0 = βij for i = 1, ..., N (ν) with j 6∈ Λs (ν) and γij = 0 = δij for i = 1, ..., N (ν) with j 6∈ Λe (ν), after relabelling the conditions (3) may be written as

K h X

i

γij yj (lj ) + δij y ′ j (lj ) = 0,

j=1

(f, g) =

i=1 0

li

f |ei g¯|ei dt,

Theorem 3.1 Let λ = −k2 , k > 0, then for y a solution of the boundary value problem (2)-(3), 1 1 ||y||L2 (G) = √ ||y|∂G ||L2 (∂G) 1 + O k 2k

 



(5)

as k → ∞, where ∂G denotes the boundary of G. Proof: Consider the second order Sturm-Liouville problem on the interval [0, li ] given by −yi′′ + qyi = λyi

(6)

with non-homogeneous Dirichlet boundary conditions and

yi (li ) = βi (k).

(7)

Let λi0 denote the least eigenvalue of (6) on (0, li ) with Dirichlet boundary conditions, yi (0) = 0 = yi (li ). Taking λ < Λ := min λi0 we have that (6)-(7) has a unique solution for each αi (k), βi (k) and each i = 1, . . . , K. From [12, Appendix A1], the fundamental solutions of (6) obeying the boundary conditions u1 (0) = 1 = u′2 (0), u′1 (0) = 0 = u2 (0)

u1 (t) = cosh kt + O (4)

for i = 1, ..., N, where N is the total number of linearly independent boundary conditions. Define L2 (G) to be the set of all f : G → C with fi ∈ L2 (0, li ) and inner product K Z X

3 Boundary Estimates

(8) (9)

are given asymptotically for large k > 0, by

i

αij yj (0) + βij y ′ j (0) +

j=1

with domain D(L) = {f | f, f ′ ∈ AC, l(f ) ∈ L2 (G), f obeying (3) }. In this setting, the formal self-adjointness of (2)(3) ensures that the operator L on L2 (G) is a selfadjoint operator, see [21, p. 77-78].

i=1,...,K

for i = 1, ..., N (ν) where N (ν) is the number of linearly independent boundary conditions at node ν.

K h X

Lf = −f ′′ + qf

yi (0) = αi (k)

αij yj (0) + βij y ′ j (0) +

X

then L2 (G) becomes a Hilbert space. The above boundary value problem on G can be reformulated as an operator eigenvalue problem, [5], by setting

u2 (t) =

1 sinh kt + O k

ekt k

!

ekt k2

,

(10)

!

(11)

with corresponding derivatives u′1 (t) = k sinh kt + O(ekt ), u′2 (t) = cosh kt + O

ekt k

!

(12) ,

(13)

uniformly with respect to t. Note that the Wronskian of u1 (t) and u2 (t) is equal to 1 for all t, i.e. u1 (t)u′2 (t) − u2 (t)u′1 (t) = 1,

u1 (li )αi (k) − βi (k) u2 (t) u2 (li ) [−αi (k)(−u1 (t)u2 (li ) + u1 (li )u2 (t))] u2 (li ) βi (k)u2 (t) . (14) u2 (li )

yi2 (t)

= +

w(t) := −u1 (t)u2 (li ) + u1 (li )u2 (t), then w is the solution of (6) with

+

R li

Thus from [12, Appendix A1], for large k > 0, !

,

u22 (li )



1 1+O k

 

(15)

.

(16)

=

αi (k) sinh k(li − t) + O





 

 



e2kli 8k3 u22 (li )

=

e2kli − 1 2k



1+O

 

1 k

 

×

βi2 (k) + α2i (k) + 4αi (k)βi (k)kli e−kli

1 e2kli (α2 + βi2 ) 1 + O k 8k3 u22 (li ) i

ekt k

!!

 



i

ek(li −t) k

α2i (k) + βi2 (k) ≥ |2αi (k)βi (k)|. +

!!#

Using (16) gives .

Z

0

li

yi2 (t) dt = ||yi ||2L2 (0,li )

Squaring this gives

=

yi2 (t) =

1 k

 

since by Schwartz’ inequality

1 [βi (k)u2 (t) − αi (k)w(t)] u2 (li ) 1 βi (k) sinh kt + O ku2 (li )

1+O

1 2li αi (k)βi (k)ekli 1+O 4 k #    2kli 2 αi (k) 1 e −1 1+O 4 k 2k

+

=

yi (t)

"



2αi (k)βi (k)ekli 1 1+O 4 k   # 2k(l −t) 2 i 1 αi (k)e 1+O . 4 k

"



Substituting (11) and (15) in (14) we obtain

=

"

1 βi2 (k)e2kt 2 2 4 k u2 (li )

βi2 (k) 1 1 1+O 2 2 4 k k u2 (li )

+ e2kli = 4k2

!!2  .

yi2 (t) dt

=

uniformly in t. Now

!!

Integrating from 0 to li gives 0

ek(li −t) k2

×

Hence for large k > 0, yi2 (t) is bounded on (0, li ) so Lebesgues dominated convergence theorem can be used, and only the pointwise limit of yi2 (t) needs to be considered for t ∈ (0, li ). For t ∈ (0, li ) and k → ∞,

Let

w(li ) = 0 w′ (li ) = −u′1 (t)u2 (li ) + u1 (li )u′2 (t) = 1.

!!

ek(li −t) k

α2i (k) sinh k(li − t) + O

+

1 w(t) = sinh k(t − li ) + O k

ek(li −t) k

sinh k(li − t) + O

yi (t) = αi u1 (t) −

+

2αi (k)βi (k) sinh kt + O

for all t.

From (7), (8) and (9)

=

+

ekt k

1 α2i + βi2 1+O 2k k 

 

.

Therefore 1 k2 u22 (li )



β 2 (k) sinh kt + O i

ekt k

!!2

||yi ||2L2 (0,li )

  ||yi |∂(0,li ) ||2L2 (∂(0,li ))  1 1+O . = 2k k

Summing over i = 1, . . . , K proves the theorem. The following corollary gives bounds for the boundary norm of solutions to the non-homogeneous boundary value problem in terms of the nonhomogeneous term. Corollary 3.2 There exists a constant C > 0 such that for k > 0 large C k

3 2

||f ||L2 (G) ≥ ||y|∂G ||L2 (∂G)

(17)

for all f ∈ L2 (G) and y the solution of −y ′′ + qy = λy + f,

(18)

where λ = −k2 , obeying the boundary conditions (3). Proof: Let Gλ denote the Green’s operator of the boundary value problem (2)-(3) and let GD λ denote the Green’s operator of the boundary value problem (2) but with Dirichlet boundary conditions at every node (i.e. y at all nodes is zero). We note that (l − λ)(Gλ − GD λ )f = f − f = 0 for f ∈ L2 (G) and where l is as given in (1). Thus (Gλ −GD λ )f is a solution of (2) and from Theorem 3.1 obeys (5) hence we obtain that, since (1 + O( k1 )) ≥ √1 for large k, 2 ||(Gλ − GD λ )f ||L2 (G) ≥ 1 √ ||[(Gλ − GD λ )f ]|∂G ||L2 (∂G) , 2 k

since GD λ vanishes on the boundary of G. Taking y = Gλ f gives (17). Remark For the system −Y ′′ + QY = λW Y + F,

(21)

with general self-adjoint boundary conditions of the form AY (0) + BY ′ (0) + CY (1) + DY ′ (1) = 0, (22) for A, B, C and D constant matrices, where either (i) Q ∈ L∞ (0, 1) is Hermitian (not necessarily diagonal), W = I and F ∈ L2 (0, 1) is Hermitian, i.e. a non-weighted system with general, Hermitian, L∞ potential or (ii) Q ∈ L∞ (0, 1) is real valued and diagonal, W is constant, real valued and diagonal and F ∈ L2 (0, 1) is diagonal, i.e. a Sturm-Liouville boundary value problem on a graph, see [7], the following result, corresponding to the above corollary, is obtained. There exists a constant C > 0 such that for k > 0 large, C 3

k2

||F ||L2 (0,1) ≥ ||Y|∂(0,1) ||L2 (∂(0,1))

(23)

for Y the solution of (21), with λ = −k2 , obeying the boundary conditions (22).

(19)

for all f ∈ L2 (G). But [(Gλ − GD λ )f ]|∂G = [Gλ f ]|∂G giving

4 Example In this section we provide an example to illustrate Theorem 3.2. We also show that (17) is the best possible estimate that can be obtained. Consider the second order differential equation

1 ||(Gλ − GD λ )f ||L2 (G) ≥ √ ||[Gλ f ]|∂G ||L2 (∂G) . 2 k

−y ′′ − λy = eαt

Now as Gλ and GD λ are both resolvent operators we have

on [0, π] where α is a constant, i.e. a graph with a single edge of length π, with the boundary conditions

||f ||L2 (G) ||f ||L2 (G) =C |λ| k2 (20) for λ → −∞, where C > 0 is a constant. Hence, combining (19) and (20), we obtain that

y(0) = 0, y ′ (π) = 0.

||(Gλ − GD λ )f ||L2 (G) ≤ C

2C

||f ||L2 (G) k2

≥ =

1 √ ||[(Gλ − GD λ )f ]|∂G ||L2 (∂G) k 1 √ ||[Gλ f ]|∂G ||L2 (∂G) k

(24)

(25)

Then the solutions of (24) on the interval [0, π] are of the form y=

eαt + aekt + be−kt , − k2

α2

where a and b are constants, and λ = −k2 .

(26)

From (26) and the boundary conditions (25), the constants a and b are given as follows αeπα + ke−πk , k(α2 − k2 )(ekπ + e−kπ ) αeπα − ke−πk . k(α2 − k2 )(ekπ + e−kπ )

a = − b =

Substituting the constants back into (26) and evaluating at 0 and π gives y(0) = 0, (k − α)eπ(k+α) + (k + α)eπ(α−k) − 2k . k(α2 − k2 )(ekπ + e−kπ )

y(π) =

We now look at the case of α = 2k. Then y(π) =

−e2πk −ke3πk + 3keπk − 2k ≃ 3k3 (ekπ + e−kπ ) 3k2

giving

||f ||L2 (G) =

s

[6] R. C ARLSON , Spectral theory and spectral gaps for periodic Schr¨odinger operators on graphs, Waves Random Media, 14 (2004), S29-S45. [7] S. C URRIE , B.A. WATSON , Eigenvalue asymptotics for differential operators on graphs, J. Com. Appl. Math., 182 (2005), 13-31. [8] S. C URRIE , B.A. WATSON , Dirichlet Neumann bracketing for boundary value problems on graphs, Electronic J. Differential Equations, 2005 (2005), 1–11 [9] N. G ERASIMENKO , B. PAVLOV, Scattering problems on non-compact graphs, Theoretical and Mathematical Physics, 74 (1988), 230-240.

e2πk ≃ . 3k2

[10] B. G UTKIN , U. S MILANSKY, Can one hear the shape of a graph, J. Phys. A: Math. Gen., 34 (2001), 6061-6068.

e4πk − 1 e2πk ≃ √ , 4k 2 k

[11] M. H ARMER , Hermitian symplectic geometry and extension theory, J. Phys. A: Math. Gen., 33 (2000), 9193-9203.

||y|∂G ||L2 (∂G) Also

[5] R. C ARLSON , Adjoint and self-adjoint differential operators on graphs, Electronic J. Differential Equations, 1998 (1998), No. 06, 1–10.

and therefore ||f ||L2 (G) 3 3 ≃ k2. ||y|∂G ||L2 (∂G) 2 Thus

3

3 ||f ||L2 (G) ≥ ||y|∂G ||L2 (∂G) .

2k 2 Showing the power term in k to be optimal. References: [1] S. AGMON , A. D OUGLIS , L. N IRENBERG Estimates near the Boundary for Solutions of Elliptic Partial Differential Equations Satisfying General Boundary Conditions. I, Comm. of Pure and Applied Math., XII (1959), 623–727.

[12] H. H OCHSTADT, On Inverse Problems Associated with Sturm-Liouville Operators, J. Differential Equations, 17 (1975), 220-235. [13] P. K UCHMENT, Graph models for waves in thin structures, Waves Random Media, 12 (2002), R1-R24. [14] C. M IRANDA , Sul Problema misto per le equazioni lineari ellittiche, Ann. Mat. Pura Appl., Ser. 4, 39 (1955), 279–303. [15] Y U .V. P OKORNYI , V.L. P RYADIEV, Some problems of the qualitative Sturm-Liouville theory on a spatial network, Russian Math. Surveys, 59:3 (2004), 515-552.

[2] F. A LI M EHMETI , A characterization of a generalized C ∞ -notion on nets, Integral Equations and Operator Theory, 9 (1986), 753-766.

[16] Y U .V. P OKORNYI , V.L. P RYADIEV, A. A L O BEID, On the oscillation of the spectrum of a boundary value problem on a graph, Math. Notes, 60 (1997), 351-353.

[3] J. AVRON , Adiabatic quantum transport in multiply connected systems, Reviews of Modern Physics., 60 (1988), 873-915.

¨ [17] J. S CHAUDER , Uber lineare elliptische Differentialgleichungen zweiter Ordnung, Math. Z., 38 (1934), 257–282.

[4] F.E. B ROWDER , Estimates and existence theorems for elliptic boundary value problems, Proc. Nat. Acad. Sci. U.S.A., 45 (1959), 365–372.

[18] J. S CHAUDER , Numerische Absch¨atzungen in elliptischen linearen Differentialgleichungen, Studia Math., 5 (1934), 34–42.

[19] M. S OLOMYAK, On the spectrum of the Laplacian on regular metric trees, Waves Random Media, 14 (2004), S155-S171. [20] J. VON B ELOW, A characteristic equation associated to an eigenvalue problem on c2 -networks, Linear Algebra and Applications, 71 (1985), 309-325. [21] K. YOSHIDA , Functional Analysis, SpringerVerlag, Berlin Heidelberg New York, 1980.

Boundary estimates for solutions of non-homogeneous boundary ...

values of solutions to the non-homogeneous boundary value problem in terms of the norm of the non-homogeneity. In addition the eigenparameter dependence ...

65KB Sizes 0 Downloads 332 Views

Recommend Documents

Boundary Final.pdf
BOYD LAKE. WELD COUNTY. ROAD 62 ... POWELL. SANCTUARY. SEA GULL. GOLDCO GRUMMAN. DRAFT. HORSE ... Page 1 of 1. Boundary Final.pdf.

boundary layer
RH and OH not only depends on the strength and spatial distribution of the. RH emissions ..... Any present-day numerical model is only capable of sim- ulating 2 ...

Boundary Final.pdf
Effective. 29 May, 2016. Page 1 of 1. Boundary Final.pdf. Boundary Final.pdf. Open. Extract. Open with. Sign In. Main menu. Displaying Boundary Final.pdf.

Growing Garlic - Boundary Garlic
on time and curing properly are all important for producing bulbs with good ... Grain straw is not recommended because it can host wheat curl mite which will ...

School District Boundary Task Force
A request was made for the department to gather information about how other ... Disclaimer: Staffing of the task force by the Department of Education and ...

Boundary Lines new.pdf
17th Parallel---------North Vietnam and South Vietnam. Hindenburg Line-------- Germany and Poland. Maginot Line-----France and Germany. Oder Neisse ...

School Boundary Map.pdf
San. Benito Rd. N. ogales. Ave. Old Adobe. Way. Horizon Ln. Dolo res Ave ... Ave. River Rd. Baggins Hill Rd. San. Diego Rd. Santa Barbara. Estrada Ave.

NO-BOUNDARY EMERGENCE AND BOOK OF CHANGE1
3. WHAT NO-BOUNDARY EMERGENCE IS. 3.1. Emergence in Stephen Hawking's ..... According to Hubble law, the recession velocity of a galaxy is in direct.

Boundary Element Formulation of Harmonic ... - Semantic Scholar
that this will allow tailoring new interpolates for particular needs. Another direction for research is ... ment methods with the software library BEMLIB. Chapman &.

Boundary Realignment Plan.pdf
Boundary Realignment Plan.pdf. Boundary Realignment Plan.pdf. Open. Extract. Open with. Sign In. Main menu. Displaying Boundary Realignment Plan.pdf.

Jones Boundary Map.pdf
Page 1 of 1. n. OSO. STAPLES. LIPES. YORKTOWN. C. R 41. EVERHART. CEDAR PASS. BOSTON. LOIRE. TIMBERGATE. VAIL. STRASBOURG. GRAFORD.

Boundary Element Formulation of Harmonic ... - Semantic Scholar
On a deeper level, BEM makes possible the comparison of trans- finite harmonic ... Solving a Dirichlet problem could seem a high price to pay, but the quality of the .... Euclidean space, and not just to some large yet bounded domain. This task ...

The First Boundary of Bangladesh.pdf
Abdul Karim,“Suba Bangla: Government and Politics,” History of Bangladesh, 1704-. 1971, Edited by Sirajul Islam, Dhaka: Asiatic Society of Bangladesh, ...

Stability Boundary for Haptic Rendering: Influence of ...
Senior Scientist. Carsten Preusche. Senior Scientist ... From the control point of view, a haptic system is a sampled-data controlled mechatronic device. Fig.

Query-Boundary Disputes.pdf
... your information and proper guidance. ATTY. SHERWIN MARC ARCEGA LINDO. Provincial Chief Legal Counsel. l. Page 3 of 4. Query-Boundary Disputes.pdf.

School District Boundary Task Force
School District Boundary Task Force. Study Assignment. In certain areas of the state, the battle over what are known as “minor boundary changes” for school.

Internal boundary layer model for the evolution of ...
Feb 5, 2012 - and ranging (lidar) data we computed z02 =10−2 m (Supplementary ... portantly, the analytical expression was not fit to sand flux data; all ...

numerical approximation of the boundary control for the ...
Oct 11, 2006 - a uniform time T > 0 for the control of all solutions is required. This time T depends on the size of the domain and the velocity of propagation of waves. In general, any semi-discrete dynamics generates spurious high-frequency oscilla

Stability Boundary for Haptic Rendering: Influence of ...
are supported by simulations and experimental data using the DLR Light-Weight ... applications, such as virtual prototyping [1] and maintainability analysis [2, .... Classical control tools have been applied to the linear system in order to obtain ..

Boundary scanning and complex dynamics - Mark McClure
with this notebook. .... [2] J. Glynn and T. Gray, The Beginner's Guide to Mathematica Version 4. Cambridge University Press, NY, 2000. BoundaryScanPP.nb. 9.

ocsb-georges-vanier-boundary-map.pdf
LANGFORD. HINES. KEIGHLEY. EAGLESON. MANNING. KANATA. WITHERSPOON. GOULBOURN FORCED. MARCH. LEGGET. HAZELDEAN. LEACOCK. RICHARDSON SIDE. BRADLEY SIDE. IRONSIDE. OAK CREEK. KANATA. SOLANDT. SANDHILL. CASTLEFRANK. PALLADIUM. HUNTMAR. RICHARDSON SIDE.

ocsb-holy-redeemer-boundary-map.pdf
THEN proceeding east along the middle of Hazeldean Road to Eagleson Road,. THEN proceeding north along the middle of Eagleson Road to Highway #417,. THEN proceeding west along the middle of Highway #417 to the point of commencement at the Palladium D

ocsb-good-shepherd-boundary-map.pdf
Page 1 of 264. SHEFFIELD. AVIATION. ROCKCLIFFE. REGIONAL ROAD 174. PAGE. TENTH LINE. WALL. WALKLEY EXTENSION. INDUSTRIAL. BATHGATE. MONTREAL. ROLLING BROOK. DOLMAN RIDGE. HIGHWAY 417. RAMSAYVILLE. CONROY. BLAIR. RUSSELL. HAWTHORNE. PORTOBELLO. NANTES

ocsb-st-isidore-boundary-map.pdf
Page 1 of 3. RIDGESIDE FARM. PINERIDGE. CONSTANCE BAY. BEACH. ALLBIRCH. MALEY. KLONDIKE. MARCH. GREENLAND. ARMITAGE. SHERRUBY.