The Inverse Spectral Transform for the Conservative Camassa Holm Flow with Decaying Initial Data

Transcription

1 Arch. ational Mech. Anal. 4 (7 5 Digital Object Identifier (DOI.7/s z The Inverse Spectral Transform for the Conservative Camassa Holm Flow with Decaying Initial Data Jonathan Eckhardt Communicated by C. De Lellis Abstract We establish the inverse spectral transform for the conservative Camassa Holm flow with decaying initial data. In particular, it is employed to prove the eistence of weak solutions for the corresponding Cauchy problem.. Introduction The inverse scattering transform [,,39] is a powerful tool for solving initial value problems for certain nonlinear partial differential equations. Originally developed for the Korteweg de Vries equation by Gardner, Greene, Kruskal and Miura [4], this method has since been successfully etended to various other completely integrable equations. It is the aim of the present article to establish the corresponding transform for the Camassa Holm equation u t u t = u u 3uu + uu (. with decaying initial data. Due to the vast amount of literature on this equation, we only refer to a selection of articles [5,4,5,9,,6,7,43,5,59,6,63,7] containing further information. The relevance of the Camassa Holm equation stems from the fact that it constitutes a model for unidirectional wave propagation on shallow water [6,4,5,53]. Unlike the Korteweg de Vries equation, it allows for smooth solutions to blow up in finite time in a way that resembles wave-breaking. This process has been described in detail [,,6,63] and is known to only happen when the quantity ω = u u changes sign. Compared to the rather tame sign-definite case (which shares a lot of similarities with the Korteweg de Vries equation [4,58,64], indefiniteness of ω causes serious complications (for eample esearch supported by the Austrian Science Fund (FWF under Grant No. J3455.

2 Jonathan Eckhardt noticeable from the discussions in [5,54,57]. Many of these problems are already apparent for the prototypical eample of a peakon-antipeakon collision [6,,68], which is a special case of so-called multi-peakon solutions [5]. On the other side, the Camassa Holm equation is known to be formally completely integrable in the sense that there is an associated isospectral problem f + f = z ωf, (. 4 where z is a comple spectral parameter. Solving the corresponding inverse problem is essentially equivalent to solving the initial value problem for the Camassa Holm equation. For this reason, it is not surprising that the complications encountered due to wave-breaking for indefinite ω reoccur when dealing with this inverse problem. In fact, despite a large amount of articles, very little is known about the inverse problem for (. in the indefinite case and almost the entire literature on this subject is restricted to strictly positive and smooth ω (in which case the spectral problem can be transformed into a standard form that is known from the Korteweg de Vries equation [4,58,64]. Apart from the eplicitly solvable finite dimensional case [5,33], only insufficient partial uniqueness results [7,9,3,34,37] have been obtained so far for the inverse problem in the indefinite case. Making use of recent progress on the inverse spectral problem for indefinite strings in [35], we will be able to overcome these difficulties and establish the inverse spectral transform for the more general two-component Camassa Holm system u t + uu + P =, μ t + (uμ = (u 3 Pu, (.3 where the auiliary function P satisfies P P = u + μ, (.4 for a class E of decaying initial data to be defined below; this should be compared to the definition of the set D in [5, Section 6], [5, Definition 3.], [45, Definition 4.]. Note that the two-component Camassa Holm system (.3 is equivalent to its more commonly used incarnation in [8,3,4,45,5] (cf.[5, Equation (4.], which is known to contain the Camassa Holm equation as a special case. Definition.. The set D consists of all pairs (u,μ such that u is a real-valued function in H ( and μ is a non-negative finite Borel measure on with μ(b u( + u ( d (.5 for every Borel set B. B

3 The Inverse Spectral Transform 3 Associated with each pair (u,μ D is a distribution ω in H ( defined by ω(h = u(h( d + u (h ( d, h H (, (.6 so that ω = u u in a distributional sense, as well as a non-negative finite Borel measure υ on defined such that μ(b = u( + u ( d + υ(b (.7 B for every Borel set B. Let us point out that it is always possible to uniquely recover the pair (u,μfrom the distribution ω and the Borel measure υ. In fact, the function u at any point can be written as u( = ω(δ, δ (s = e s, s, (.8 which then allows us to determine the Borel measure μ as well from (.7. Now the class E is defined by imposing an additional growth restriction on pairs in D. Definition.. The set E ± consists of all pairs (u,μ D such that e ± ( u ( u( d + e ± dυ( <. (.9 Furthermore, the set E is defined as the intersection of E + and E. Justified by the La pair formulation of the two-component Camassa Holm system in [8,3,5] and the results in [33], we will consider the spectral problem f + 4 f = z ωf + z υ f. (. Due to the low regularity of the coefficients, it is not clear how this differential equation has to be interpreted, which is why we clarify this matter in Appendi A.Basic properties of the spectral problem (like realness and discreteness of the spectrum σ for eample will be discussed in Section along with some further necessary conditions. In the following section, we will solve the corresponding inverse problem for the class of coefficients corresponding to E, giving a complete characterization of all possible spectral data. More precisely, we will establish a one-to-one correspondence between E and a class of spectral data eplicitly described by Theorem 3., which will be shown to be a homeomorphism with respect to suitable topologies. Finally, in Section 4 we will introduce the conservative Camassa Holm flow as a certain isospectral flow on E and show that its integral curves define weak solutions of the two-component Camassa Holm system.

4 4 Jonathan Eckhardt As an immediate consequence of the solution of the inverse spectral problem in Theorem 3., we will see that the phase space E decomposes into a foliation of invariant isospectral sets Iso(σ, each of which can be parametrized by the set { κ σ λ e κλ λẇ (λ < }, (. where W is the entire function defined by the infinite product in (.6 and the dot denotes differentiation. In terms of these coordinates, the conservative Camassa Holm flow on each set Iso(σ becomes the simple linear flow given by κ λ =, λ σ. (. λ These facts are reminiscent of the fact that the conservative Camassa Holm flow can be viewed as a completely integrable infinite dimensional Hamiltonian system. Notation. For integrals of a function f which is locally integrable with respect to a Borel measure ν on an interval I, we will employ the convenient notation y [,y f dν, y >, f dν =, y =,, y I, (.3 [y, f dν, y <, rendering the integral left-continuous as a function of y. If the function f is locally absolutely continuous on I and g denotes a left-continuous distribution function of the Borel measure ν, then we have the integration by parts formula y f dν = gf y y g(s f (s ds,, y I, (.4 which will be used frequently throughout this article. Given a discrete set σ of nonzero reals, we denote with n σ (r the number of all λ σ with modulus not greater than r. Furthermore, we introduce the notation λ = lim r λ, (.5 λ r provided that the limit eists. Similarly, subject to eistence, we shall write ( z ( = lim z λ r λ λ r, z C, (.6 where the limit is meant to be taken in the topology of locally uniform convergence. The limit in (.6 eists if and only if the limit in (.5 eists and the sum We denote with σ the set of all real-valued sequences κ ={κ λ } indeed by σ.

5 The Inverse Spectral Transform 5 λ (.7 is finite. In this case, upon denoting the entire function in (.6 with W,wehave = Ẇ (, λ λ = Ẇ ( Ẅ (. (.8. The Direct Spectral Problem In the present section, we are going to introduce the spectral quantities that will linearize the conservative Camassa Holm flow on E and derive their basic properties. To this end, we fi an arbitrary pair (u,μ E and first recall the definition of the distribution ω in (.6 as well as the one of the Borel measure υ in (.7. We introduce the spectrum σ associated with the pair (u,μas the set of all those numbers z C for which there is a nontrivial solution f of the differential equation f + 4 f = z ωf + z υ f (. that belongs to H (. egarding the precise meaning and some basic properties of this differential equation we refer to the discussion in Appendi A. A first consequence of the growth restriction on (u,μin (.9 is the eistence of particular solutions of the differential equation (. with prescribed asymptotics. Theorem.. For every z C there is a unique solution φ ± (z, of the differential equation (. with the asymptotics φ ± (z, e, ±. (. The derivative of φ ± (z, is integrable and square integrable near ±. Proof. To begin with, let us introduce the diffeomorphism η ± : by η ± (ξ = ln( ξ, ξ, where denotes the open negative/positive semi-ais. Furthermore, let a ± be a measurable function on such that (we set α ± = u ± u as in Appendi A a ± (ξ = α ±(η ± (ξ ξ = u (η ± (ξ ± u(η ± (ξ ξ (.3 for almost all ξ (note that the right-hand side is well-defined almost everywhere and define the non-negative Borel measure β ± on via setting β ± (B = ξ dυ η ±(ξ = e ± dυ( (.4 B η ± (B

6 6 Jonathan Eckhardt for every Borel set B. The growth restriction on (u,μin (.9 implies that the function a ± is square integrable and that the measure β ± is finite (thus β ± can be etended to a Borel measure on the closure of by setting β ± ({} =. Consequently (see, for eample, [3, Section.8], [8, Theorem.], [38, Theorem A.], there is a unique matri solution Y ± (z, on of the integral equation ( ξ Y ± (z,ξ= + z ξ + z ( ( a± (s a ± (s a ± (s Y ± (z, s ds Y ± (z, s dβ ± (s, ξ, (.5 for each z C. It follows from a version of Gronwall s inequality (see, for eample, [8, Lemma.3], [38, Lemma A.] that the function Y ± (z, satisfies the estimate Y ± (z,ξ e ±(ξ z, ξ, (.6 where denotes the ma norm and the non-negative function ± is given by ξ ( ± (ξ = ma, a ± (s ξ ds + dβ ±, ξ. Plugging this estimate back into the integral equation (.5, we furthermore get ( Y ±(z,ξ z ± (ξe ±(ξ z, ξ. (.7 For the top-right entry of Y ± (z,, we are actually able to sharpen the estimate to Y±, (z,ξ ( ξ z e ± (ξ z a ± (s ds ξ, ξ. Upon plugging this back into the integral equation (.5 one more time, we obtain Y±, (z,ξ+ z ξ z ξ ± (ξe ±(ξ z, ξ. (.8 Net, we introduce the matri function U ± (z, on by ( e ± U ± (z, = e Y ± (z,η± (,. From the integral equation (.5, a substitution (use, for eample, [3, Theorem 3.6. and Corollary 3.7.] as well as the definitions of a ± and β ±, one sees that Y ± (z,η± ( y y ( = z α± (s e s α ± (s e ±s Y α ± (s ± (z,η± (s ds y ( + z e ±s Y ± (z,η± (s dυ(s for all, y. Upon employing the integration by parts formula (.4, one sees that U ± (z, is a solution of the integral equation (A.5 and thus the system (A.4.

7 The Inverse Spectral Transform 7 Now under the additional assumption that z is non-zero, let θ ± (z, and φ ± (z, be the solutions of the differential equation (. such that ( θ± (z, zφ ± (z, U ± (z, =± z θ ± [±],, (z, φ[±] ± (z, guaranteed to eist by Lemma A.. For future purposes, let us note the bounds θ± (z, e, θ [] ± (z, e, φ± (z, e ±, φ [] ± (z, e± (.9 e 4 ±( e z,, that follow from the estimates (.6 and (.8 fory ± (z,, also using that ± u( = e s α ± (s ds ( ± e,. Furthermore, the inequalities in (.7 and (.8 turn into θ± (z, e ±, z θ ± [±] (z, e±, φ± (z, e ±, φ [±] ± (z, e± ± ( z ± e e ±( e z,. (. In particular, we see that φ ± (z, has the required asymptotics (. aswellas θ ± (z, e ±, ±, which implies that φ ± (z, is indeed uniquely determined by the asymptotics in (.. Finally, the fact that the derivative of φ ± (z, is integrable and square integrable near ± follows from the last bound in (.9. It remains to set φ ± (, = e, θ± (, = e ±, for all and note that the claim is obvious in the case when z is zero. For every z C we introduce the comple number W (z as the Wronskian of the two solutions φ + (z, and φ (z,, that is, in such a way that W (z = φ + (z, φ (z, φ + (z, φ (z, (. for almost all ; see Corollary A.4. It follows readily that the set of zeros of W coincides with the spectrum σ. Thus, the net result implies that σ is a discrete set of nonzero reals with a convergence eponent of at most one. Corollary.. The functions φ ± (, and W are real entire of Cartwright class with only nonzero and real roots for each fied. Moreover, the function zφ (z, φ + (z,, z C\, (. W (z is a meromorphic Herglotz Nevanlinna function.

8 8 Jonathan Eckhardt Proof. Since the entries of Y ± (,ξare real entire of Cartwright class with only real roots (cf. [55, Section ], the same holds for φ ± (, upon noting that ( lim φ ±(z, =±e ± Y ±, z, e lim = e = φ± (,. z z z Because the function φ ± [] (, is real entire of Cartwright class as well, so is W as W (z = φ + (z, φ [] (z, φ[] + (z, φ (z,, z C. Net, we observe that the meromorphic function ± φ[] ± (z, zφ ± (z,, z C\, is a Herglotz Nevanlinna function. In fact, to this end one just needs to evaluate its imaginary part and use (A. as well as the vanishing asymptotics of the solution φ ± (z, near ±. Since we may write W (z zφ (z, φ + (z, = φ[] + (z, zφ + (z, φ[] (z, zφ (z,, z C\, this shows that the function in (. is a Herglotz Nevanlinna function as well. In particular, this guarantees that W has only nonzero and real roots indeed. emark.3. Although we will not prove this here, let us mention that it is possible to show that the eponential type of the entire function W is simply given by ρ( d, (.3 where ρ is the square root of the adon Nikodým derivative of the absolutely continuous part of the Borel measure υ (with respect to the Lebesgue measure. As the spectrum alone will not be enough, we need to introduce further spectral quantities. Since the solutions φ + (λ, and φ (λ, are linearly dependent for each fied λ σ in view of Corollary A.4, there is a unique nonzero c λ such that φ (λ, = c λ φ + (λ,,, (.4 which will be referred to as the coupling constant associated with λ. It will often be more convenient to work with the logarithmic coupling constant κ λ, defined by κ λ = ln c λ (.5 instead. We will see below that one can always recover (the sign of the coupling constant c λ from the quantity κ λ, provided the spectrum σ is known as well. Furthermore, we introduce the right/left (modified norming constant γλ,± via γλ,± = ω( φ ± (λ, + λ φ ± (λ, dυ(. (.6

9 The Inverse Spectral Transform 9 Upon employing (A. and the asymptotics from Theorem., one sees that λγλ,± = φ ± (λ, d + φ ± (λ, d + λ φ ± (λ, dυ( >. 4 (.7 The following result gives a relation between all our spectral quantities. Lemma.4. For each λ σ we have the relation Ẇ (λ = c λ ± γ λ,±. (.8 Proof. Let us fi an arbitrary z C and introduce the function W ± (z, = φ ± (z, φ [] (z, φ ± [] (z, φ (z,,, where the dot denotes differentiation with respect to the spectral parameter. From the bounds in (. and Cauchy s integral formula, we obtain the estimates φ ± (z, e ±, φ ± [±] (z, e± ( ( z + ± e e ±( e ( z +,. Upon writing φ (z, as a linear combination of φ ± (z, and θ ± (z,, wesee from the bounds in (.9 that the function W ± (z, tends to zero as ±. Net, we note that from (A.6 and (A.8 one gets φ ± (z, y = y φ [] ± (z, y = y zu (sφ ± (z, s + φ [] ± ( 4 zu(s z u (s y z φ ± (z, s dυ(s, (z, s ds, φ ± (z, s zu (sφ ± [] (z, s ds for all, y, and after differentiating with respect to z also φ ± (z, y y = zu (s φ ± (z, s + φ ± [] (z, s + u (sφ ± (z, s ds, y ( φ ± [] (z, y = 4 zu(s z u (s φ ± (z, s zu (s φ ± [] (z, s ds y ( u(s + zu (s φ ± (z, s + u (sφ ± [] (z, s ds y y z φ ± (z, s dυ(s z φ ± (z, s dυ(s. In conjunction with the integration by parts formula (.4, this gives W ± (z, y = y + z u(sφ + (z, sφ (z, s + u (s (φ + (z, φ (z, (s ds y φ + (z, sφ (z, s dυ(s.

10 3 Jonathan Eckhardt Upon letting y ± in this equation, we obtain Ẇ (z = W (z, W + (z, = u(sφ + (z, sφ (z, s + u (s (φ + (z, φ (z, (s ds (.9 + z φ + (z, sφ (z, s dυ(s, which readily yields the claimed identity (one should also note that the functions u and u are integrable due to the growth restriction in (.9 and thus the bounds in (.9 guarantee that all integrals eist indeed. In particular, the previous result shows that all zeros of the entire function W are simple. As a Cartwright class function, the Wronskian W thus admits a product representation (see, for eample, [6, Section 7.] of the form W (z = ( z, z C. (. λ Upon invoking the identities in (.8, this fact allows us to read off trace formulas for the spectrum σ from the derivatives of the function W at zero. Proposition.5. The first two trace formulas are: λ = u( d, λ = dμ. (. Proof. In view of (.8, we immediately obtain the first identity from (.9 and we are left to compute the second derivative of W at zero. For this purpose, we first introduce the entire function V ± by V ± (z = θ ± (z, φ [] [] (z, θ ± (z, φ (z,,, z C, so that we may write φ ± (z, =±W (zθ (z, + V (zφ (z,,, z C. (. In conjunction with (.9 and (., this allows us to estimate φ ± (z, e ± W (z, φ ± [±] (z, e± ± W (z 4e 5 (±e ± z ( ( zw(z ±e ± ( z e ± + + V (z e ± (.3 for all and z C. From this we see that each of the three products φ + (z, φ (z,, φ + (z, φ [] (z,, φ[] + (z, φ (z,,

11 The Inverse Spectral Transform 3 is bounded uniformly in and locally uniformly in z C. Thus, we may differentiate (.9 under the integral and in order to evaluate the derivative at zero, we first note that Ẏ ± (,ξ= ξ ( a± (s ξ a ± (s ds + a ± (s ( dβ ± (s, ξ, (.4 which follows from (.5. Moreover, for the top-right entry we even have Ÿ ±, (,ξ= ξ ξ a ± (s ds 4 ξ s After performing substitutions, this turns into the identities φ ± (, = e ± α (s ds, φ [] ± (, =± e a ± (r dr ds, ξ. ± α ± (s ds, for all. Plugging them into the differentiated integral in (.9, we arrive at ( Ẅ ( = u( + u ( d + dυ u( d, which proves the remaining claim in view of (.8. Before we are able to solve the inverse spectral problem, we need to derive one more necessary condition. More precisely, we will show that the growth restriction on the pair (u,μin (.9 implies certain asymptotic behavior of the spectral data. Proposition.6. We have the identity e ± ( u ( u( d + e ± dυ( = λ λγ. (.5 λ,± Proof. To begin with, we introduce the function m ± by m ± (z = V ±(z zw(z, z C\. From the estimates in (.3, we see that φ ± (z, e ± W (z, φ [±] ± (z, e± W (z,, where the convergence is locally uniform in z C. In much the same manner as above (upon writing θ ± as in (. and using (.9 aswellas(. to estimate, one obtains inequalities similar to (.3 forθ ± and concludes that θ ± (z, e ± ±V± (z, θ ± [±] (z, e± V± (z,, where the convergence is again locally uniform in z C. Thus, we may write m ± (z =± lim θ ± (z, zφ ± (z, = lim ξ Y ±, (z,ξ, z C\, (.6 Y ±, (z,ξ

12 3 Jonathan Eckhardt that is, the function m ± is the Weyl Titchmarsh function for (.5; cf. [35, Equation (6.]. In particular, it is a meromorphic Herglotz Nevanlinna function [35, Lemma 5.] with simple poles at all points λ σ (note that zero is not a pole since V ± vanishes at zero with residues given by res λ m ± = V ±(λ λẇ (λ = λγλ,±, λ σ, in view of Lemma.4. After a substitution, we obtain from (.4 that θ ± (, = e ± ± α ± (s ds, θ [] ± (, = e± ± α (s ds, for every, which implies that V ± ( =. Upon recalling the integral representation formula for Herglotz Nevanlinna functions (see, for eample, [66, Section 5.3], this guarantees that the function m ± admits the representation (use also [35, Lemma 7.] to conclude that there is no linear term present m ± (z = z λ(λ z λγ, z C\. (.7 λ,± Clearly, we may just as well write the function m ± as the limit m ± (z =± lim θ ± [±] (z, zφ ± [±] = lim (z, ξ Y ±, (z,ξ, z C\. (.8 Y ±, (z,ξ Now we observe that (.4 givesforeveryξ the epansion Y ±, (z,ξ ξ ξ Y ±, (z,ξ = z a ± (s ds + z dβ ± + O(z, z. Since the convergence in (.8 is uniform for all z C\ that lie in a small neighborhood of zero, we infer, after a substitution, that m ± has the epansion m ± (z = z e ± ( u ( u( d + z e ± dυ( + O(z, z. Upon differentiating (.7 and letting z, we obtain the identity in (.5. In this contet, let us also mention the following result which characterizes the subclass of E that gives rise to purely positive/negative spectrum. Proposition.7. The spectrum σ is positive/negative if and only if the Borel measure υ vanishes identically and the distribution ω is non-negative/non-positive. In this case, the distribution ω can be represented by a non-negative/non-positive finite Borel measure on (for simplicity denoted with ω as well and λ = dω. (.9

13 The Inverse Spectral Transform 33 Proof. If υ vanishes identically and ω is non-negative/non-positive, then (.6 and (.7 show that the spectrum is positive/negative. Conversely, if the spectrum is positive/negative, then [35, Lemma 7.] shows that the measures β + and β vanish identically as well as that the functions a + and a have non-decreasing/nonincreasing representatives. Clearly, the measures β + and β vanish identically if and only if so does υ. Givensomeh H ( with compact support, we set h ± (ξ = ξh(η ± (ξ, ξ, (.3 and note that ω(h = h + (ξa +(ξ dξ h (ξa (ξ dξ. + This shows that ω is a non-negative/non-positive distribution and therefore can be represented by a non-negative/non-positive Borel measure on. Now for every k N let h k H ( be the piecewise linear function such that h k is equal to one on [ k, k], equal to zero outside of [ k, k + ] and linear in between. Then from h k dω = u(h k ( d + u (h k ( d we see that ω is a finite measure as well as (.9 upon letting k. An identity similar to (.9 also holds under the sole assumption that the distribution ω can be represented by a real-valued Borel measure on (again denoted with ω as well with finite total variation. More precisely, one has the equality λ = dω (.3 in this case, which follows readily from the first trace formula in Proposition The Inverse Spectral Problem We are now going to solve the corresponding inverse spectral problem for the class E. Let us point out that the given sufficient conditions on the spectral data are also necessary in view of Corollary. (in conjunction with well-known properties of entire functions of Cartwright class; for eample, see [, Chapter 8] or [6, Theorem 7..] as well as Proposition.6 and Lemma.4. Thus, we indeed obtain a complete characterization of all possible spectral data for the class E. Theorem 3.. Let σ be a discrete set of nonzero reals such that the limit n σ (r lim r r (3. ecall that n σ (r denotes the number of all λ σ with modulus not greater than r.

14 34 Jonathan Eckhardt eists in [, and such that the entire function W is well-defined by 3 W (z = Moreover, for each λ σ let κ λ such that the sum e κλ λ λẇ (λ ( z, z C. (3. λ (3.3 is finite. Then there is a unique pair (u,μ E such that the associated spectrum coincides with σ and the logarithmic coupling constants are κ λ for each λ σ. Proof. Uniqueness. Since the given spectral data uniquely determines the Weyl Titchmarsh function m in view of (.7 and Lemma.4, the uniqueness part in [35, Theorem 6.] shows that the function a and the Borel measure β are uniquely determined as well. It follows readily from the definition of β in (.4 that this also uniquely determines the Borel measure υ. Furthermore, because of ξ u(η ± (ξ = a ± (ss ds, ξ, (3.4 ξ we conclude that the function u is uniquely determined too and thus so is μ. Eistence. Let us introduce the meromorphic Herglotz Nevanlinna function m(z = z e κ λ λ(λ z λẇ (λ, z C\. From the eistence part of [35, Theorem 6.], we obtain a real-valued and locally square integrable function a on [, and a non-negative Borel measure β on [, with β({} = such that the function m is the Weyl Titchmarsh function for ( ξ ( a(s Y (z,ξ= + z a(s Y (z, s ds a(s ξ ( (3.5 + z Y (z, s dβ(s, ξ [,, that is, if Y (z, denotes the unique solution of the integral equation (3.5 for each z C, then the function m is given by Y (z,ξ m(z = lim ξ Y (z,ξ, z C\. In order to state a required fact from [35] in a concise way, we introduce the set (, that consists of all ξ (, such that the function a is not equal 3 If the limit in (3.eistsin[,, then this is the case if and only if the sum λ eists.

15 The Inverse Spectral Transform 35 to a constant almost everywhere near ξ or such that ξ belongs to the topological support of β.let f, f L (, be orthogonal to all h L (, that satisfy ξ h(s ds =, ξ, (3.6 and functions g, g on [, be square integrable with respect to β. Upon setting Y (λ, ξ F i (λ = f i (ξ dξ + Y (λ, ξg i (ξ dβ(ξ, λ σ, i =,, λ we have the identity F (λf (λ e κ λ λẇ (λ = f (ξ f (ξ dξ + g (ξg (ξ dβ(ξ (3.7 since the spectral transform introduced in [35, Section 5] is a partial isometry. Motivated by the relation (3.4, we define the real-valued function u on via ξ u(ln(ξ = a(ss ds, ξ (,, (3.8 ξ so that u belongs to Hloc ( and satisfies a(ξ = u (ln(ξ u(ln(ξ ξ for almost all ξ (,. Furthermore, we define the Borel measure μ on by (.7, where υ is the Borel measure on such that that we have β(b = e dυ( ln(b for every Borel set B (,. Since the function a is square integrable near zero and the Borel measure β is finite near zero, it follows readily that c e ( u ( + u( c d + e dυ( < for every c. In particular, this guarantees that the function u lies in H ( near (more precisely, one sees from (3.8 that u decays eponentially near and since u + u is clearly square integrable near, we conclude that u is as well and that the measure μ is finite near. We are now left to verify that the pair (u,μactually belongs to E. In fact, in this case it is readily seen that a coincides with the function a as introduced by (.3 and that β coincides with the Borel measure β as introduced by (.4. From (.6, we then see that the function m coincides with the Weyl Titchmarsh function m for this pair as introduced in the proof of Proposition.6. Upon taking (.7 and Lemma.4 into account, one then may conclude that the pair (u,μindeed gives rise to the desired spectral data.

16 36 Jonathan Eckhardt In order to show that the pair (u,μ belongs to the set E, let us first consider the special case when a is equal to some constant a almost everywhere near and β vanishes near. From this we readily infer that ( za (ξ ξ Y (z,ξ= z(ξ ξ za (ξ ξ Y (z,ξ + za (ξ ξ as long as ξ and ξ are close enough to. Thus, the function m is given by m(z = a Y (z,ξ + Y (z,ξ a Y (z,ξ + Y (z,ξ, z C\, and evaluating the limit as z in this equation shows that a has to be zero, which immediately implies that the pair (u,μ belongs to E. In particular, this proves the claim in the case when the set σ is finite (since then the set is finite as well. Thus, back in the general case, for every k N we may find a (u k,μ k E such that the corresponding Weyl Titchmarsh function m k, is given by m k, (z = λ k e κ λ z λ(λ z λẇ (λ, z C\. With the functions a k, defined as in (.3, we infer from [35, Proposition 6.] that ξ lim k a k, (s ds = ξ a(s ds, locally uniformly for all ξ [,. Inviewof(3.4 and (3.8, this shows that the functions u k converge pointwise to u and since they are uniformly bounded in H ( by Proposition.5, a compactness argument shows that u lies in H (. To verify also the remaining growth restrictions on (u,μin the general case, we now may assume that there is an increasing sequence n such that ξ n = e n belongs to for every n N.Wefisomen N and define the function J n by J n (ζ, z = Z (z,ξ n Z (ζ, ξ n Z (z,ξ n Z (ζ, ξ n z ζ, ζ, z C\, where Z(z, is the Weyl solution given by ( zw(z Z(z,ξ= Y (z,ξ, ξ [,, V (z for every z C, and V is the entire function defined in such a way that V (z = zw(zm(z, z C\. Because the function Z(z, is a Weyl solution for every z C\, we note that J n can be rewritten as (see the proof of [35, Lemma 5.] J n (ζ, z = Z ξ n (z, s z Z (ζ, s ζ ds + ξ n Z (z, sz (ζ, s dβ(s

17 The Inverse Spectral Transform 37 for all ζ, z C\. In particular, this implies that the entire function E n defined by E n (z = zψ(z, n iψ [+] (z, n, z C, is a de Branges function, where ψ(z, is the solution of the differential equation (. so that (compare the proof of Theorem. ( ( zψ(z, e ψ [ ] = (z, e Z(z, e,, for each z C, with the quantities ω and υ defined as in (A. and (.7. The reproducing kernel K n in the corresponding de Branges space B n is given by K n (ζ, z = J n (ζ, z + ψ(z, n ψ(ζ, n, ζ, z C\. An integration by parts and the integral equation (3.5 show that for every z C\ the function Z (z, restricted to (ξ n, is orthogonal to all functions h L (, that satisfy (3.6. Thus, the identity (3.7 for the particular functions { f i (s = Z (z i,ξ n e n, s [,ξ n, z i Z (z i, s, s [ξ n,, g i (s = {, s [,ξ n, Z (z, s, s [ξ n,, with z = z and z = ζ gives (note that Z (λ, = V (λy (λ, for all λ σ K n (z,λ V (λ K n (ζ, λ V (λ e κλ λẇ (λ = K n(z,ζ = K n (z,, K n (ζ, Bn, ζ, z C\. The values of V on the set σ are readily evaluated and we obtain F(λ e κ λ λẇ (λ = F B n, F B n, (3.9 after employing simple linearity, continuity and density arguments. Now let H be a locally integrable, trace normed, real, symmetric and nonnegative definite matri function on [, such that if M(z, denotes the unique solution of the integral equation M(z, t = ( + z t ( H(sM(z, s ds, for every z C, then we have M (z, t lim t M (z, t = z e κ λ λ(λ z λẇ (λ, t [,, z C\. Such a function H is guaranteed to eist by the solution of the inverse problem for canonical systems due to de Branges; see [3, Theorem XII], [69, Theorem.4]. Finiteness of the sum in (3.3 implies(see[56, Lemma 5.5], [66, Theorem 6.7] and

18 38 Jonathan Eckhardt note that eistence of (3. and (3. together establish [56, Condition (C]; cf. [, 8.3.6] that the entire function W belongs to the Cartwright class. Moreover, since m is a Herglotz Nevanlinna function, so does the function V and thus also E n for every n N.Inviewof[3, Theorem VII] in conjunction with (3.9 and upon employing [9, Theorem I], we obtain a t n (, so that ( zψ(z, n ψ [+] (z, n = Ɣ n ( M (z, t n M (z, t n, z C, for a real matriɣ n with det Ɣ n =. Evaluating at zero shows that ( e n Ɣ n = ε n e n for some ε n. We note that the sequence t n may be chosen in such a way that it is non-increasing because the sequence of de Branges spaces B n is non-increasing with respect to inclusion. Thus, for every n N we have on the one side ( zψ(z, ψ [+] (z, = Ɣ M(z, t M(z, t n Ɣ n ( zψ(z, n ψ [+] (z, n, z C, and on the other side (due to Lemma A.also ( ( zψ(z, zψ(z, n ψ [+] = (z, + (z,, n ψ [+], z C, (z, n where ± (z,, n denotes the matri valued solution of the system (A.4 such that ± (z, n, n is the identity matri. In view of [8, Problem ], this gives M(z, t M(z, t n = Ɣ +(z,, n Ɣ n, z C. Differentiating with respect to z and evaluating at zero, we obtain t ( n H(s ds = e s ε ds e ε n e n ε e t n ε n e n n ( + u (s u(s ds ε e + ε n e n n + e s ( u (s u(s ( n ( ds + e s dυ(s. Since H is integrable near zero, the top-right entry shows that the integrals n u (s u(s ds are uniformly bounded for all n N. Moreover, from the bottom-left entry we now infer that the quantities ε n e n are uniformly bounded for all n N as well. Finally, the top-left entry implies that the sum n e s ( u (s u(s n ds + e s dυ(s

19 The Inverse Spectral Transform 39 is uniformly bounded for all n N, which concludes the proof since n. For a discrete set σ of nonzero reals such that the limit (3.eistsin[, and such that the entire function W is well-defined by (3., we define the isospectral set Iso(σ as the set of all those pairs (u,μ E whose associated spectra coincide with the set σ. It is an immediate consequence of Theorem 3. that the isospectral set Iso(σ is in one-to-one correspondence with the set 4 { σ = κ σ by means of the bijection given by λ } e κλ λẇ (λ <, (3. (u,μ {κ λ }. (3. Of course, the solution of the inverse spectral problem can easily be formulated in terms of the right/left norming constants as well. However, the picture (that is, the condition corresponding to (3.3 on the asymptotic behavior looks less symmetric. Corollary 3.. Let σ be a discrete set of nonzero reals such that the limit (3. eists in [, and such that the entire function W is well-defined by (3.. Moreover, for each λ σ let γλ,± such that λγ λ,± > and the sums λ λγλ,±, γλ,± λ λẇ (λ, (3. are finite. Then there is a unique pair (u,μ E such that the associated spectrum coincides with σ and the right/left norming constants are γλ,± for each λ σ. In order to avoid misunderstandings, let us point out that this corollary has to be read as two separate statements according to the plus-minus alternative, that is, we are only able to prescribe either the right or the left norming constants. emark 3.3. If the set σ is finite, then all conditions in Theorem 3. and Corollary 3. are trivially satisfied. The solution of the inverse spectral problem in this case has a particular simple form and can be written down eplicitly in terms of the spectral data; see [33, Section 4]. We will see net that these kinds of solutions can be used to approimate solutions in the general case in a certain way. We conclude this section by establishing a continuity property for the inverse spectral transform. In order to state it, let (u k,μ k belong to E for every k N and denote all corresponding quantities in an obvious way with an additional subscript. 4 To be precise, this has to be interpreted appropriately in the trivial case when σ is empty.

20 4 Jonathan Eckhardt Proposition 3.4. We have 5 λ( + λ λγ λ( + λ k k,λ,± λγ, (3.3 λ,± χ(λ + λ λγ χ(λ + λ k k,λ,± λγ, (3.4 λ,± as k for all functions χ C( such that the limit of χ(λ as λ eists and is finite if and only if the functions u k converge to u pointwise and ( s e s e ±r ( u k (r u k(r s dr + e ±r dυ k (r ds ± ± ± ( s e s e ±r ( u (r u(r s dr + e ±r dυ(r ds ± ± (3.5 as k for all. In this case, the functions u k converge to u even locally uniformly and the Borel measures μ k converge to μ in the sense of distributions. Proof. It follows from [35, Proposition 6.] that the first condition on convergence of (3.3 and (3.4 is equivalent to ξ ( s a k,± (r dr + ξ s a k,± (s ds dβ k,± ds ξ ξ ± a ± (s ds, (3.6 ( s s a ± (r dr + dβ ± ds, (3.7 as k for all ξ, in which case the convergence is uniform as long as ξ stays bounded. In conjunction with (3.4, one sees that this implies that the sequence u k ( converges to u( for all. Moreover, upon employing a substitution in (3.7, we readily infer that (3.5 holds as k for all as well. Both of these convergences are uniform as long as stays away from. In addition, for every smooth function h on that vanishes near and such that h has compact support we have the identity h dμ k = h(α k,± ( d ± h(u k (u k ( d + h dυ k = h ± (ξa k,± (ξ dξ h (u k ( d + h ± dβ k,± ξ ( s s = h ± (ξ a k,± (r dr + dβ k,± ds dξ h (u k ( d, 5 Let us mention that the following condition is equivalent to locally uniform convergence m k,± m ± of the corresponding Weyl Titchmarsh functions given as in (.7.

21 The Inverse Spectral Transform 4 with h ± given by (.3. Taking the limit k on the right-hand side gives h dμ k h dμ, k, (3.8 which shows that the Borel measures μ k converge to μ in the sense of distributions. For the converse direction, we now assume that the functions u k converge to u pointwise and that (3.5 holds as k for every. As observed above, the latter condition is equivalent to the fact that (3.7 holds as k for all ξ. Upon epressing the integrals in (3.6 in terms of u k and u by using (.3aswellas a substitution, we infer (use Lebesgue s dominated convergence theorem and note that convergence of (3.7 implies a uniform bound on the monotone integrands which yields a sufficient estimate for the functions u k that (3.6 holds as k for all ξ, which establishes the equivalence in the claim. We endow E with the unique first countable topology 6 such that the sequence (u k,μ k converges to (u,μif and only if the functions u k converge to u pointwise and both of the plus-minus alternatives in (3.5 hold as k for all. Corollary 3.5. If the sequence (u k,μ k converges to (u,μin E, then the functions u k converge to u uniformly as well as weakly in H (, the Borel measures μ k converge to μ in the weak topology 7 and the sequence μ k ( converges to μ(. Proof. We have seen in the first part of the proof of Proposition 3.4 that the given assumption implies that the functions u k converge to u uniformly. After choosing a suitable partition of unity (two functions are sufficient we infer from (3.8 that the sequence μ k ( converges to μ(. In particular, the sequence μ k ( is bounded which allows to deduce the remaining claims in a straightforward manner. 4. The Conservative Camassa Holm Flow Let us define the conservative Camassa Holm flow on E as a mapping : E E (4. in the following way: given a pair (u,μ E with associated spectrum σ and logarithmic coupling constants κ λ for every λ σ,aswellasat, the corresponding image t (u,μunder is defined as the unique pair in E (guaranteed to eist by Theorem 3. for which the associated spectrum coincides with σ and the logarithmic coupling constants are κ λ + t, λ σ. (4. λ 6 For eample, take the topology inherited from the injection (u,μ (m +, m when the space of Herglotz Nevanlinna functions is equipped with the topology of locally uniform convergence. 7 We regard the space of finite Borel measures on as the dual space of the Banach space C (.

22 4 Jonathan Eckhardt Of course, this definition is motivated by the well-known simple time evolution of spectral data for classical solutions of the Camassa Holm equation, which constitutes the essence of our method of solution and is due to the underlying completely integrable structure; see [4, Section 6]. It is clear that the isospectral sets Iso(σ are invariant under the conservative Camassa Holm flow. Moreover, the bijection in (3. takes the conservative Camassa Holm flow on Iso(σ to a simple linear flow on σ whose solutions are eplicitly given by (4.. In conjunction with these facts, the trace formulas in Proposition.5 give rise to conserved quantities for the flow; cf. [4,57]. Proposition 4.. The two functionals (u,μ u( d, (u,μ dμ, (4.3 on E are invariant under the conservative Camassa Holm flow. Since the Wronskian W is uniquely determined by the spectrum σ, it is invariant under the conservative Camassa Holm flow. Therefore, we see that the functional (u,μ ρ( d (4.4 on E is invariant as well (recall the notation from emark.3. In particular, this shows that the property of the function ρ vanishing almost everywhere (or equivalent, the Borel measure υ being singular with respect to the Lebesgue measure is preserved. This case corresponds to the subclass of global conservative solutions of the Camassa Holm equation in [5,5], whereas the function ρ is present for the general class of global conservative solutions of the two-component Camassa Holm system in [45]. Proposition 4.. The conservative Camassa Holm flow is continuous. Proof. Let t, t k for every k N such that t k t as k and suppose that the sequence (u k,μ k converges to (u,μin E so that (3.4 holds as k for all functions χ C( such that the limit of χ(λ as λ eists and is finite. In conjunction with Proposition.5, we infer from Corollary 3.5 that there is an ε> such that the intersections ( ε, ε σ and ( ε, ε σ k are empty for every k N. Given any continuous function τ on \{} such that the limit of τ(λ as λ eists and is finite, we choose a function χ ± C( such that χ ± (λ = τ(λe ± t λ, λ ε, as well as constants K, T such that τ(λ K for all λ with λ ε and such that t k T for all k N. We then may estimate

23 The Inverse Spectral Transform 43 τ(λ + λ τ(λ k λγk,λ,± t e k λ + λ λγλ,± e λ t K ε e T ε tk t + λ λγ k k,λ,± + χ ± (λ + λ λγk,λ,± χ ± (λ + λ λγ λ,± k for every k N. Since it follows readily from our assumptions that the right-hand side always converges to zero as k, we infer that λ( + λ k k λγ k,λ,± e t k λ χ(λ + λ λγk,λ,± t e k λ λ( + λ λγ λ,± e t λ χ(λ + λ, λγλ,± e λ t as k for all functions χ C( such that the limit of χ(λ as λ eists and is finite. In view of Proposition 3.4 and the definition of the flow,this implies that the corresponding images t k (u k,μ k converge to t (u,μ. For the global conservative solutions of the two-component Camassa Holm system in [45], continuity results similar to Proposition 4. were obtained in [45, Theorem 5.] and [46, Theorem 6.7]. However, the topologies used in [45,46] are defined in a much more intricate and less eplicit way; cf. [45, Lemma 6.4 and Lemma 6.5]. As the main result of this section, we are now going to show that the integral curve t t (u,μ for any fied initial data (u,μ E defines a weak solution of the two-component Camassa Holm system (.3. To this end, let us denote the pair t (u,μ with (u(, t, μ(, t for every t, so that u can be regarded as a function on. Note that the integral curve t t (u,μ is continuous by Proposition 4., which guarantees that the function u is at least continuous. Theorem 4.3. The pair (u,μis a weak solution of the two-component Camassa Holm system (.3 in the sense that for every test function ϕ Cc ( we have ( u(, t u(, tϕ t (, t + + P(, t ϕ (, t d dt =, (4.5 ϕ t (, t + u(, tϕ (, t dμ(, t dt = ( u(, t u(, t P(, t ϕ (, t dt,, (4.6

24 44 Jonathan Eckhardt where the function P is given by P(, t = e s u(s, t ds + e s dμ(s, t,, t. ( Proof. Let us denote all the quantities corresponding to (u,μ in an obvious way with an additional subscript. If the spectrum σ is a finite set, then it follows from [33, Section 5] that the pair (u,μis a weak solution of the two-component Camassa Holm system (.3. Otherwise, we define a pair (u k,,μ k, E for every k N in such a way that the associated spectrum coincides with σ [ k, k] and the right/left norming constants are γ,λ,± for every λ σ [ k, k]. Since the intersection σ [ k, k] is a finite set, the corresponding pairs (u k,μ k obtained from the integral curves t t (u k,,μ k, satisfy ( uk (, t u k (, tϕ t (, t + + P k (, t ϕ (, t d dt =, (4.8 ϕ t (, t + u k (, tϕ (, t dμ k (, t dt ( (4.9 uk (, t = u k (, t P k (, t ϕ (, t d dt, for every test function ϕ Cc (, where the function P k is given by P k (, t = e s u k (s, t ds + e s dμ k (s, t,, t. 4 4 Moreover, our definitions and Proposition.5 guarantee the bounds u k (, t μ k (, t,, t, k N, λ and from Proposition 3.4 we infer that the functions u k (, t converge to u(, t locally uniformly and that the Borel measures μ k (, t converge to μ(, t in the weak topology (to conclude this, we also used that the sequence μ k (, t is bounded for every fied t. In particular, this shows that the sequence of functions P k is bounded as well and converges to P at least pointwise. Finally, upon passing to the limit k in (4.8 and (4.9, we infer that the pair (u,μis a weak solution of the two-component Camassa Holm system (.3. When the Borel measure υ corresponding to the pair (u,μ vanishes identically and the distribution ω is non-negative/non-positive, then this property is preserved by the conservative Camassa Holm flow in view of Proposition.7; cf. [, Corollary 3.3]. In this case, the function u is a weak solution of the Camassa Holm equation (. in the sense that for every test function ϕ C c ( we have ( u(, t u(, tϕ t (, t + + P(, t ϕ (, t d dt =, (4.

25 The Inverse Spectral Transform 45 where the function P is given by P(, t = e s u(s, t ds + e s u (s, t ds,, t. 4 (4. Under these assumptions, global weak solutions in this sense have been obtained before in [,6]. In the more challenging general case, eistence of global weak solutions as in Theorem 4.3 has been established in [5,45,5] by means of an elaborate transformation to Lagrangian coordinates. emark 4.4. If the spectrum σ associated with the pair (u,μ is a finite set, then we recover the special class of conservative multi-peakon solutions [49], which can be written down eplicitly in terms of the spectral data [33]. We have seen in the proof of Theorem 4.3 that in the general case, our weak solution (u,μ of the two-component Camassa Holm system (.3 can be approimated by a sequence of conservative multi-peakon solutions in a certain way; cf. [47,48]. We conclude this section with a comment on the long-time behavior of the conservative Camassa Holm flow. Upon employing the method introduced in [36], it is possible to show that our weak solution (u,μof the two-component Camassa Holm system (.3 asymptotically splits into a (in general infinite train of single peakons, each corresponding to an eigenvalue λ σ of the underlying spectral problem. For classical solutions of the Camassa Holm equation (., this behavior was anticipated by McKean [63] (in accordance with numerical observations in [7] and first proved in [37] (see also [5,33] for the multi-peakon case and [6] for a particular class of low-regularity solutions. Appendi A. The Basic Differential Equation Throughout this appendi, let u be a real-valued function in Hloc ( and υ be a non-negative Borel measure on. We define the distribution ω in Hloc ( by8 ω(h = u(h( d + u (h ( d, h Hc (, (A. so that ω = u u formally, and consider the ordinary differential equation f + 4 f = z ωf + z υ f, (A. where z is a comple spectral parameter. Of course, this differential equation has to be understood in a distributional sense in general; cf. [35,44,67]. 8 The space H c ( consists of all those functions in H ( which have compact support.

26 46 Jonathan Eckhardt Definition A.. A solution of (A.isafunction f Hloc ( such that f (h ( d + f (h( d = z ω( fh + z fhdυ 4 for every function h Hc (. (A.3 In order to make this notion of solution more accessible, we will first show that the differential equation (A. is equivalent to the first order system ( ( ( ( F α ± = F ± + z ± α ± υ F α ± ± + z F ±, (A.4 where we introduced α ± = u ± u. Since υ is allowed to be a genuine Borel measure, this system has to be understood as a measure differential equation [3,8, 38,65] in general. Solutions of such an equation are not necessarily continuous but only of locally bounded variation and, in order to guarantee the unique solvability of initial value problems, one has to impose a suitable normalization. Here we will require solutions to be left-continuous, which is implicitly contained in the following definition: a solution of the system (A.4 is a function F ± : C with locally bounded variation such that F ± y = y ( y ( α± (s F ± (s ds + z α ± (s F α ± (s ± (s ds y ( + z F ± dυ (A.5 for all, y. In this case, the first component of F ± is clearly locally absolutely continuous and the second component is left-continuous; cf. (.3. With this notion of solutions, initial value problems for the system (A.4 are always uniquely solvable; see, for eample, [3, Section.8], [8, Theorem.], [38, Theorem A.]. Lemma A.. If the function f is a solution of the differential equation (A., then there is a unique left-continuous function f [] such that f [] ( = f ( zu ( f ( (A.6 for almost all and the function ( zf f [] ± ( zu (A.7 f is a solution of the system (A.4. Conversely, if the function F ± is a solution of the system (A.4, then its first component is a solution of the differential equation (A.. Proof. First one notes that for every c, Equation (A.3 takes the form ( h ( f ( zu ( f ( f (s ds + z u(s f (s + u (s f (s ds 4 c c + z f dυ d = c

27 The Inverse Spectral Transform 47 upon integrating by parts. Now from this we infer that a function f Hloc ( is a solution of (A. if and only if there is a c and a constant d C such that f ( zu ( f ( = d + 4 c f (s ds z u(s f (s + u (s f (s ds c z f dυ (A.8 for almost all.soif f is a solution of (A., then this guarantees that there is a unique left-continuous function f [] such that (A.6 holds for almost all. Let us denote the second component in (A.7 with f [±] so that f [±] ( = f ( f ( + zα ±( f ( for almost all. Then it is straightforward to show that y y ± f [±] (s ds + z α ± (s zf(s α ± (s f [±] (s ds ( =± zu y f y y f (s ds + z u(s f (s + u (s f (s ds 4 for all, y. In combination with (A.8, we end up with f [±] y =± + z y y f [±] (s ds + z zf dυ y c α ± (s zf(s α ± (s f [±] (s ds for all, y, which shows that (A.7 is a solution of the system (A.4. Now suppose that F ± is a solution of the system (A.4 and denote the respective components with subscripts. The first component of (A.5 shows that F ±, belongs to Hloc ( with zf ±, ( = F ±, ( ± F ±,( zα ± (F ±, ( (A.9 for almost all. In combination with the second component of (A.5thisgives F ±, ( =± F ±,( zα ± (F ±, ( zf ±, (c ± zf ±, (s ds c z α ± (s F ±, (s ds z α ± (s zf ±, (s ds z F ±, dυ c for some c and almost all. Plugging (A.9 twice into this equation, we see that (A.8 holds with f replaced by F ±, for some constant d. As noted before, this guarantees that F ±, is a solution of the differential equation (A.. The auiliary function f [] introduced in Lemma A. allows us to state the following basic eistence and uniqueness result for the differential equation (A.. c c

28 48 Jonathan Eckhardt Corollary A.3. For every c and d,d C there is a unique solution f of the differential equation (A. with the initial conditions f (c = d and f [] (c = d. (A. If d,d and z are real, then the solution f is real-valued as well. Proof. If z is zero, then the solutions of the differential equation (A. aregiven eplicitly and the claim is obvious. Otherwise, it follows from unique solvability of initial value problems for the system (A.4 in conjunction with Lemma A.. We are now left to introduce the Wronskian W ( f, g of two solutions f, g of the differential equation (A.. Although the following result is essentially a consequence of Lemma A. and the fact that the traces of the matrices in (A.5 vanish, we provide an independent direct proof. Corollary A.4. For any two solutions f, g of the differential equation (A. there is a unique comple number W ( f, g such that W ( f, g = f (g ( f (g( (A. for almost all. This number is zero if and only if the functions f and g are linearly dependent. Proof. First of all, note that if f is a solution of the differential equation (A. and h Hloc (, then (A.8 and the integration by parts formula (.4 show that f [] h y y = f (sh (s ds + y f (sh(s ds 4 y y (A. z u(s f (sh(s + u (s( fh (s ds z fhdυ for all, y. Ifg is another solution of the differential equation (A., then f (g ( f (g( = f (g [] ( f [] (g( (A.3 for almost all. Employing (A., one sees that the right-hand side of (A.3 is constant, which proves eistence of the number W ( f, g. The remaining claim follows in a standard way upon utilizing the uniqueness part in Corollary A.3. Let us mention that one can easily get rid of the constant potential term in the first order system (A.4 by scaling solutions with the matri function ( e e ±,, (A.4 at the cost of rescaling the remaining matri coefficients in (A.4 as well. It is even possible to further transform (A.4 into a canonical system in standard form ( G ± = zh ±G ±, (A.5 with a locally integrable, non-negative and trace normed matri function H ±. Since this is of no need for our purposes, we will not do this here but refer to [35, Proof of Theorem 6.], where such a transformation has been employed.