1. Introduction
The well-known limit circle/limit point classification introduced by Weyl tells us how many boundary conditions are necessary to define a self-adjoint realization of a three-coefficient Sturm–Liouville differential expression given by
|
|
|
(see Hypothesis 2.1 for details). Limit point implies no boundary conditions are needed at the given endpoint, while limit circle implies one boundary condition is needed. Since the introduction by Weyl, the limit circle/limit point dichotomy of Sturm–Liouville expressions has been thoroughly studied- for recent treatments of Sturm–Liouville theory with encyclopedic references, we refer the interested reader to [13] and [39].
The purpose of the present paper is to study an extension of the classic binary classification of limit circle/limit point in the case of nonoscillatory Sturm–Liouville differential expressions. This is achieved by means of the regularization index which gives a natural finer classification of certain limit point endpoints: those with finite regularization index lead to eigenvalues satisfying Weyl asymptotics (i.e., grow like ); those with infinite regularization index have eigenvalues satisfying growth quicker than (i.e., trace class resolvent); those for which the regularization index is undefined have nonempty essential spectrum and/or eigenvalues that grow slower than the above.
In order to define the regularization index for the endpoint resp., we need to assume resp., has self-adjoint realizations with trace class resolvents. This spectral condition, essential in our work, turns out to be equivalent with the following simple integrability condition on the product of the principal and nonprincipal solutions near the respective endpoint (see Theorem 4.2): Assume that the equation on , is nonoscillatory at the endpoint resp., for some , and that for some
|
|
|
where , resp., , are any principal, resp. nonprincipal solutions of near the endpoint resp., . This integrability assumption then allows one to iteratively construct a spectral parameter power series (i.e., a Taylor series in the spectral parameter ) for solutions of the Sturm–Liouville problem . This series can be equivalently viewed as a type of Born expansion, and we show that in a certain precise sense this series is well-behaved if and only if the aforementioned trace class condition is satisfied (see Remark 4.3). The regularization index at the singular endpoint (resp., ) is then defined by comparing the growth in of the coefficients of the power series of the principal and nonprincipal solutions as (resp., . See Definition 5.1.
The regularization index turns out to be well-behaved under Darboux transforms as seen in Theorem 9.1. Hence, one important implication is that a finite regularization index allows one to quantify how far certain limit point endpoints are away from being Darboux transformed to a limit circle endpoint, which in turn can be regularized in the sense of Niessen and Zettl (see [35] and [39, Thm. 8.3.1]), thus the appropriateness of the name regularization index.
In particular, our work extends the notion of regularization to include limit point nonoscillatory endpoints of Sturm–Liouville expressions with finite regularization index that can be Darboux transformed (equivalently, transformed into Schrödinger form). As a corollary, we obtain Weyl eigenvalue asymptotics for this class of problems. For more information on Darboux transforms directly related to the current study, we refer to [2], [16] (see also [26] and the extensive list of references therein) and [1], [18], [19] (in the context of exceptional orthogonal polynomials).
We also study the structure of the Weyl -function under our regularization process. In particular, we obtain an alternative proof that the Weyl -functions in the case of finite regularization index is in the subclass consisting of generalized Nevanlinna–Herglotz functions with negative squares (where is the floor function and ), no nonreal poles, and the only generalized pole of nonpositive type at infinity. This extends some of the results for specific examples studied in the series of papers [24]–[27]. In addition to these papers, for more information on singular Weyl -functions we refer to [8], [17], and [25].
The notion of the regularization index has appeared (sometimes implicitly) in various previous works. The prime example of a Sturm–Liouville operator for which the index is particularly useful is the perturbed spherical Schrödinger operator (or Bessel operator). Some earlier works include the papers of Fulton [7], Fulton-Langer [8], and Kurasov-Luger [30], where analytic perturbations where studied and the authors relied on the Frobenius method. Here the regularization index can be explicitly computed in terms of the roots of the corresponding indicial equation. Kostenko, Sakhnovich, and Teschl in a series of papers [23]–[27] included nonanalytic perturbations and used methods more inline with the present paper. In fact, our perturbation result Theorem 5.4 can be viewed as a natural generalization of the perturbative approach used in [23, Lem. 2.2] for spherical Schrödinger operators. See Remark 5.5 for more details.
The first explicit definition of an index, seemingly equivalent to the one defined in the present paper, seems to originate from the work of Kaltenbäck and Woracek on canonical systems and Pontryagin spaces of entire functions in [21, 22]. This index, denoted by the greek letter , can be associated to certain canonical Hamiltonian systems which encompass the Sturm–Liouville operators treated here as a special case. Of particular significance is the paper of Winkler and Woracek [37], where an accessible method for computing the index is provided and the paper of Langer and Woracek [32], where the special case of Sturm–Liouville operators is treated. In particular, Langer and Woracek use essentially the same recursion to define as we do to define , however subtle differences remain; see, for instance, Open Problem 6.5 and Remark 6.5. We nonetheless believe both notions to be equivalent and fully agree with the assessment of the authors of [37] that ‘limit point but finite index’ is in many respects similar to the limit circle case, which becomes even more apparent through our regularization process. Interestingly, when the index is infinite, certain examples can still share properties of the limit circle case such as Weyl asymptotics- see Section 11.3 which includes inverse quartic potentials.
We do not use the theory of canonical systems or Pontryagin spaces in our work, though we certainly believe that there are interesting connections to these areas which deserve further attention. We also avoid the notion of rigged (distributional) Hilbert spaces and supersingular perturbations (see [3], [29], [30, App. A], [33]). Instead our proofs mainly rely on classical ODE-methods for absolutely continuous functions. This is motivated by the fact that our Definition 5.1 of the regularization index relies exclusively on growth properties of classical solutions to near the endpoints, rather than on their -integrability or membership in a rigged Hilbert space.
The present paper is organized as follows:
-
•
In Section 2 we give some background and introduce the main integrability, equivalently trace class, assumptions that will be used throughout the paper.
-
•
In Section 3 we construct a spectral parameter power series representation for the principal solution to the Sturm–Liouville problem and show its convergence in Proposition 3.3.
-
•
In Section 4 we start with the crucial Theorem 4.2 giving us equivalent characterizations of the trace class resolvent condition in terms of properties of the principal and nonprincipal solutions. We also list an array of consequences of Theorem 4.2 for the properties of the entire nonprincipal solution, in particular, Corollaries 4.5 and 4.7.
-
•
In Section 5 we introduce the regularization index in Definition 5.1 and discuss some simple examples in Remark 5.1. We then proceed to prove a stability result in Theorem 5.4, generalizing earlier work on perturbed spherical Schrödinger operators.
-
•
In Section 6 we relate the regularization index to the classic limit circle/limit point classification. This is done in Theorem 6.1. We then study the relationship between the regularization index and the index used in [32], [37]. This relationship is encapsulated in Open Problem 6.5.
-
•
In Section 7 we introduce in Definition 7.1 a useful choice of normalized system of entire solutions, which we label ‘naturally normalized system’. This notion plays an important role in the subsequent sections. We also give a more intuitive characterization of this normalization in Theorem 7.6 (see in particular (7.4)).
-
•
In Section 8 we introduce additional hypotheses to study the relationship between Darboux transforms and naturally normalized systems. The crucial result of this section is Corollary 8.5 showing that the natural normalization is preserved under Darboux transforms.
-
•
In Section 9 we determine how the regularization index changes under a Darboux transform depending on the (non)principality property of the seed function (see Theorem 9.1). We also illustrate how applying a series of Darboux transforms can be viewed as a regularization process, showing that Weyl asymptotics hold for problems with finite regularization indices.
-
•
In Section 10 we use the results on Darboux transforms of the previous section to explicitly compute Weyl -functions for problems with finite regularization indices.
-
•
In Section 11 we demonstrate our results by considering a variety of examples including generalized Bessel operators, Jacobi differential operators, and Schrödinger operators on the half-line with power potentials. We also provide an example with an infinite regularization index for which Weyl asymptotics still holds (Mie-type potentials) and consider the Laguerre operator at for which our main hypothesis is not satisfied.
-
•
Appendix A contains certain technical proofs not included in the main text.
We include a few open problems throughout the paper.
2. A trace class integrability Condition
We begin by recalling the typical integrability hypotheses that we will make throughout.
Hypothesis 2.1.
Let and suppose that are Lebesgue measurable functions on
such that the following items – hold:
a.e. on , .
a.e. on , .
is real-valued a.e. on , .
Given Hypothesis 2.1, we study Sturm–Liouville operators associated with the general,
three-coefficient differential expression
(2.1) |
|
|
|
As such, the Wronskian of and , for , is defined by
|
|
|
with
|
|
|
denoting the first quasi-derivative of a function .
In the following we will drop the a.e. from equalities between functions in .
Let us now introduce maximal and minimal operators in associated with in the usual manner as follows.
Definition 2.2.
Assume Hypothesis 2.1. Given as in (2.1), the maximal operator in associated with is defined by
|
|
|
|
|
|
|
|
|
The preminimal operator in associated with is defined by
|
|
|
|
|
|
|
|
|
One can prove that is closable, and one then defines the minimal operator as the closure of . We have that .
It is known (see, e.g., [35]) that if the equation
(2.2) |
|
|
|
is nonoscillatory near (resp., ), meaning that its solutions have finitely many zeros in a vicinity of the respective endpoint, then there exists an up to constant multiples unique solution (resp., ) of (2.2) satisfying
|
|
|
for any linearly independent solution (resp., ) of (2.2). In this case (resp., ) is called the principal solution of (2.2) at (resp., ), and (resp., ) is called a nonprincipal solution of (2.2). Note that the nonoscillatory condition (2.2) near (resp., ) is equivalent to the semiboundedness of one (hence any) self-adjoint realization of (resp., ).
We now come to the main spectral condition of the present paper. We say that satisfies the trace class property at (resp., at ) if and only if every self-adjoint realization of (resp., ) for some (hence any) has trace class resolvent for some (hence any) in the resolvent set . The main goal of the present paper is to study the implications of the trace class property in the semibounded case:
Hypothesis 2.3.
Assume that satisfies the trace class property at resp., and that self–adjoint realizations of resp., with are semibounded.
As we will demonstrate, it is more practical to work instead with Hypothesis 2.4 stated below, as it is computationally more tractable. We will eventually prove in Theorem 4.2 that Hypothesis 2.4 is in fact equivalent to Hypothesis 2.3 at the respective endpoint, giving an easy criterion for the trace class property of .
Hypothesis 2.4.
Assume that the equation is nonoscillatory at the endpoint resp., for some , and that for some
(2.3) |
|
|
|
where , resp., , are any principal, resp. nonprincipal solutions of (2.2) near the endpoint resp., .
It is important to clarify that as our analysis is local, we will mostly require condition (2.3) to hold only at one of the endpoints, which we conventionally take as . In case of doubt, we will explicitly state that we require Hypothesis 2.4 at (resp., at ). We now add some context for this hypothesis.
Remark
Notice that Hypothesis 2.4 holds at limit circle nonoscillatory endpoints for every by definition.
We will show in Corollary 3.6 that if Hypothesis 2.4 holds for some it holds for all . Therefore, if in (2.1), we can choose without loss of generality in Hypothesis 2.4, and one confirms that linearly independent solutions are given by and with . So we now consider two case distinctions:
Suppose . Then so is principal and is nonprincipal.
Hence in this case Hypothesis 2.4 at is equivalent to assuming that (cf. [32, Def. 7.1])
(2.4) |
|
|
|
which in particular implies that .
Suppose . Then exists, solves , and satisfies so that is principal and is nonprincipal.
Hence in this case Hypothesis 2.4 is equivalent to (cf. [32, Def. 8.1])
(2.5) |
|
|
|
where might or might not be in . Therefore, if , then Hypothesis 2.4 holds at if and only if one of (2.4) or (2.5) holds.
If in (2.1), then Hypothesis 2.4 at is equivalent to assuming
|
|
|
(cf. [32, Def. 9.3]).
3. Properties of the principal solution
The goal of the present section is to construct a fundamental solution satisfying
(3.1) |
|
|
|
which is principal at the endpoint and entire in (a similar construction can be performed at the endpoint ). No additional requirements are needed at . The key to this construction is the following technical lemma which we will use repeatedly in the present paper. While we believe it to be known, we were unable to find this result in the literature, so we include the proof in the appendix.
Lemma 3.1.
Let , where and is open. Denote by the maximal operator associated with .
-
Let be given such that for all . Moreover, assume that for , with being holomorphic in . Then the mapping is an -valued holomorphic mapping and has locally around a series expansion
(3.2) |
|
|
|
where each is in and
(3.3) |
|
|
|
-
Assume that has locally the series representation (3.2) in the space with satisfying (3.3) in particular is an -valued analytic mapping. Then for is in and satisfies for .
We now construct the solution via the infinite power series given by
(3.4) |
|
|
|
In fact, rewriting (3.1) as
|
|
|
we see that (3.4) is the usual Born series, where we view as the coupling constant for the potential (see [15] and Remark 4.3).
We note that clearly depends on the choice of , so a more precise notation would be (see (3.2)). However, to keep the notation short we will suppress this -dependence and simply write as is customary with and . It will turn out that the choice of does not play any significant role (see Cor. 3.6).
We remark that other spectral parameter power series have been discussed in [28], specifically, the numerical aspects regarding eigenvalue problems (see also the review [31]). An equivalent construction also appeared in [32] in relation to the index mentioned in the Introduction.
Assuming Hypothesis 2.4 holds at , we define , where is a principal solution of (2.2). That is, we begin by constructing the series (3.4) about a point such that is nonoscillatory (though this can be extended to all by Corollaries 3.4–3.6). We then define iteratively
(3.5) |
|
|
|
where is a nonprincipal solution of (2.2) satisfying .
In the following, we fix a such that , have no zeros on . We remark that in this case we also have (see [35, Thm. 2.2(iii)])
(3.6) |
|
|
|
In Lemma 3.2 below we prove that the integral (3.5) indeed exists.
As we assume Hypothesis 2.4, we can define the function
(3.7) |
|
|
|
Note that for . This leads to the following:
Lemma 3.2.
Assume Hypothesis 2.1 and that Hypothesis 2.4 holds at . Let be chosen such that have no zeros on . Then the following estimates hold for
(3.8) |
|
|
|
Proof.
We proceed by induction. The estimate is trivial for and let us assume it holds for up to . Observe that as by assumption , do not have zeros for , it follow from (3.6) that
(3.9) |
|
|
|
Using now the monotonicity of together with the induction hypothesis we obtain
|
|
|
|
|
|
|
|
|
|
|
|
which finishes the proof.
∎
We now want to show that the series (3.4) is indeed entire in . From the previous lemma it follows that (3.4) is convergent on . Moreover, by definition and a direct calculation shows that
|
|
|
In particular for it follows by Lemma 3.1 that
|
|
|
To show that (3.4) converges not only in but in fact defines an entire function in for all , let us choose an , with small enough such that . Consider the entire system of solutions and of satisfying
|
|
|
Let us now define
|
|
|
Note that is holomorphic in its first argument and satisfies . By a standard uniqueness results for differential equations we must have for . Hence for any fixed it follows that can be analytically continued to a holomorphic function in the disc of radius around . Letting and thus , we see that is indeed entire in for all . To extend this result to , observe that we can write using Lemma 3.1 ,
|
|
|
with and , for . Again from the uniqueness of solutions to differential equations, we can iteratively conclude that for . We have thus shown the following:
Proposition 3.3.
The infinite series (3.4) converges for all , , and defines a function which is entire in and satisfies
|
|
|
It should be noted that while the series (3.4) converges for all , the estimate in (3.8) will not hold in general for .
We now note a few immediate corollaries from the construction of .
Corollary 3.4.
Assume Hypothesis 2.1 and that Hypothesis 2.4 holds at . Then defined by (3.4), (3.5) satisfies
(3.10) |
|
|
|
In particular, is nonoscillatory for all .
Note that this immediately implies that (3.4) can be interpreted as a perturbative Born series, in the sense that higher-order corrections , become negligible in the limit . This should be contrasted with the case of discrete spectrum but non-trace class resolvents (see Theorem 4.2 and the subsequent Remark 4.3).
The converse of the above corollary will be stated in Corollary 4.3. From the previous corollary, we also conclude the following.
Corollary 3.5.
Assume Hypothesis 2.1 and that Hypothesis 2.4 holds at . Then is principal at for all .
Proof.
We know that in the nonoscillatory case at , a solution of is principal at if and only if the function is not integrable near the endpoint (see e.g. [35, Thm. 2.2(ii)]). Due to (3.10), is not integrable near if and only if is not integrable for any .
∎
Note that as is principal at and is nonoscillatory at for all , we can obtain a nonprincipal solution of via the formula , where is chosen such that does not vanish on . In particular, the asymptotic behavior of nonprincipal solutions for is already dictated by the corresponding behavior of the principal solution (cf. the proof of Lem. 4.6). Thus, (3.10) also implies the independence of Hypothesis 2.4 from the generalized eigenvalue .
Corollary 3.6.
The Hypothesis 2.4 is independent of , that is, if it holds for one it will hold for all .
As previously pointed out, Corollaries 3.4–3.6 now imply that the choice does not play any significant role in the iterative construction of . We also obtain the following from [17, Lem. 3.2].
Corollary 3.7.
Assume Hypothesis 2.1 and that Hypothesis 2.4 holds at . Then all self-adjoint realizations of the restriction to an interval with have a purely discrete spectrum.
Remark
The inverse of Corollary 3.7 does not hold. A simple counterexample is given by the Laguerre differential expression, for which all self–adjoint realizations have purely discrete spectrum, but the principal solution at has asymptotically different behavior for different generalized eigenvalues (see (11.1)), hence Hypothesis 2.4 does not hold. The details can be found in Section 11.5.
We next turn to the properties of a second linearly independent fundamental solution .
4. Properties of the nonprincipal solution
By Corollary 3.7 and Remark 3.7 we know that the following hypothesis is weaker than Hypothesis 2.4.
Hypothesis 4.1.
Assume that all self–adjoint realizations of the restriction to an interval with and the Dirichlet boundary condition at have a purely discrete spectrum.
As shown in [17], Hypothesis 4.1 is equivalent to the existence of an entire fundamental system of solutions , of , real on the real axis, such that is principal for all and . The tilde indicates that our standard Hypothesis 2.4 is not assumed, and no additional normalization conditions on and are imposed. We now state the following:
Theorem 4.2.
Assume Hypotheses 2.1, 4.1 and let , be chosen as above. Then the following are equivalent:
-
Hypothesis 2.4 at for some
-
for all
-
for all
-
for all
-
for all and
-
Hypothesis 2.3 at .
Moreover, in conditions – ‘for all ’ can be replaced by ‘for some distinct ’.
Proof.
The equivalence between – is rather simple, however point requires more technical arguments. We provide the complete proof in Appendix A.
∎
We can now state the converse of Corollary 3.4.
Corollary 4.3.
Assume Hypothesis 2.1 and let be an entire fundamental solution of which is principal at for all . If satisfies
|
|
|
then Hypothesis 2.4 holds at and is equal to constructed via (3.4), (3.5) up to a multiplicative constant.
Remark
We already observed in Corollary 3.4 that the Born series given through (3.4) and (3.5) has the convenient property of the leading term being dominant as , that is, higher order terms can be viewed as small corrections. Theorem 4.2 further tells us that this happens if and only if has self-adjoint realizations with trace class resolvents, meaning that Hypothesis 2.4 is the most general condition under which one can expect a well-behaved Born series with the spectral parameter as the coupling constant. We find it interesting that the condition for the mere existence of an entire fundamental solution which is principal at is significantly weaker, and only requires self-adjoint realizations of to have a purely discrete spectrum (see Hypothesis 4.1). Being entire, will again have an everywhere convergent power series expansion of the form (3.4), however with not necessarily given through (3.5). As in the absence of the trace class resolvent condition the behavior of must necessarily depend on the spectral parameter due to Theorem 4.2 , it follows that higher order terms cannot be viewed as small corrections for , despite the series being convergent. See Section 11.5 for an explicit example of this phenomenon.
Returning to the normalization (3.10), we obtain the following corollary.
Corollary 4.5.
Assume Hypothesis 2.1 holds and let satisfy the normalization (3.10). Then any entire fundamental solution satisfying will also satisfy
(4.1) |
|
|
|
|
Proof.
That follows immediately from (A) and (A.3) in the proof of Theorem 4.2 in Appendix A.
∎
Note implies that is linearly independent of and hence nonprincipal.
For technical reasons, we will also need that converges locally uniformly for . This is shown next.
Lemma 4.6.
Denote by . Then as , the entire function converges locally uniformly in to the constant function .
Proof.
Note that locally in we can write
|
|
|
where is sufficiently close to such that does not vanish, and is holomorphic. Now as for we have with the error being locally uniform in , we conclude that the same is true for .
∎
As a corollary of Lemma 4.6 we can now prove the following (cf. Lemma 3.2).
Corollary 4.7.
Consider the power series expansion of with respect to ,
|
|
|
Then
(4.2) |
|
|
|
Proof.
Note that as converges locally uniformly for , we can conclude that locally uniformly for and . In particular,
|
|
|
∎
Due to Lemma 3.1 we know that and for . Hence we can write a general expression for in terms of
|
|
|
|
(4.3) |
|
|
|
|
Note that the two integrals above sum up to the usual formula involving the Green’s function, however we prefer to keep these integrals separate.
In the proof of Lemma 5.3 we will show that , while the constant will be determined later and depends additionally on the choice of . Note that it is crucial that as for all , to make sure that the first integral in (4.3) exists due to Hypothesis 2.4.
We finish this section with a few technical results on the behavior of as , which will be used in the following section.
Lemma 4.8.
The functions and are nonoscillatory as for all . Additionally, the function is nonoscillatory as for .
Proof.
We proceed by induction. The claim is clearly true for by Theorem 4.2 . Let us assume it is true up to . Then we can use the general formula
|
|
|
to obtain
|
|
|
By the induction hypothesis, the right-hand side is nonoscillatory as . Hence, is monotonic as , implying that it is nonoscillatory. To see that is nonoscillatory, observe that
|
|
|
where we already know that the right-hand side is nonoscillatory. As before this implies that is monotonic, in particular nonoscillatory near . As is nonoscillatory, the same must be true of .
The proof for is similar but no longer requires induction. For we have
|
|
|
hence is monotonic, implying that is nonoscillatory, as .
∎
Note that we have also shown in the previous proof that is monotonic as , and we will use this fact later.
Corollary 4.9.
The function is monotonic as . In particular the limit
exists in .
Analogously, using that is nonoscillatory it follows that is monotonic as for . We need this fact in the proof of Lemma 5.3.
Corollary 4.10.
The function is monotonic as for .
5. The regularization index
We now come to the main definition of the present paper. Observe that as by the principality of and the nonprincipality of , the regularization index is well-defined.
Definition 5.1 (Regularization index).
Assume Hypothesis 2.1 and that Hypothesis 2.4 holds at . Let be given via (3.4), (3.5) and take any entire nonprincipal solution satisfying which is real for . Then we define the regularization index of at the endpoint to be the smallest non-negative integer such that
|
|
|
|
|
|
|
|
in case such an integer exists. Otherwise
|
|
|
and we set .
The appropriateness of the terminology ‘regularization index’ will become more apparent in Theorem 9.2 and Section 10. For now, we note that, loosely speaking, this index (when ) allows one to quantify how far a limit point endpoint is away from being transformed to a limit circle endpoint via Darboux transforms. We will see in Theorem 6.1 that exactly corresponds to the limit circle case, which, if singular, can then be regularized in the sense of Niessen and Zettl (see [35] and [39, Thm. 8.3.1]). This has implications for the spectral theory of self-adjoint realizations, which will be discussed in Sections 9 and 10.
It turns out that the regularization index at (resp., ) depends only on as explained
in the following remark.
Remark
Note that any other system
, satisfying the assumptions of Definition 5.1 is related to , via
|
|
|
where and is an entire function that is real on the real axis. It is easy to see that the choice , lead to the same regularization index, so . The question whether the regularization index depends on the choice of is the content of Corollary 5.5. It turns out that is independent of , meaning that the regularization index at resp., depends only on the Sturm–Liouville differential expression .
Returning to Remark 2.4, we provide a few general examples demonstrating how the regularization index depends on the behavior of near singular endpoints. Here and in the following the notation (for ) will be shorthand for .
Let . Assume , , and with . By Remark 2.4 , we must restrict the powers to satisfy the integrability conditions given there. Independently of the integrability of , this condition reads . If , that is , and , then is constant and (times a possible logarithm) for . Thus one concludes in this case. If , that is , and , then and (times a possible logarithm) so that in this case. Combining these two cases gives the regularization index as in general. The logarithms show up if and either for , or for . Notice that the index is always finite and for .
Let . If one assumes and with , then Hypothesis 2.4 holds by Remark 2.4 (note is necessary as otherwise would be integrable near contradicting its principality). Furthermore, is limit circle for and (times a possible logarithm) for so that . Once again, the index is finite and for .
The following lemma is one of the main structural results of the present paper and will enable us to interpret the regularization index as an extension of the binary limit circle/limit point classification.
Lemma 5.3.
Assume Hypothesis 2.1 and that Hypothesis 2.4 holds at . Then
(5.1) |
|
|
|
Proof.
We will use the representation (4.3). Let us assume that . First, observe that
|
|
|
|
|
|
|
|
Where we defined . We know from Corollary 4.7 and Corollary 4.10 that monotonically as . In particular, for and sufficiently close to . Thus it follows that
|
|
|
Next, let us consider
|
|
|
|
|
|
|
|
where . Now observe that as , we have by Definition 5.1 that , and this convergence is monotonic due to Corollary 4.9. In particular, we can assume that is chosen close enough to such that for . An application of dominated convergence gives us
|
|
|
It now follows from (4.2) that necessarily . Moreover, as we have . Hence, (5.1) follows, finishing the proof.
∎
Note that we still have some freedom in choosing the , which is expected as is only unique up to additions of entire multiples of . We will see in Section 7 that further conditions need to be imposed on through the choice of to guarantee that (5.1) holds for all (see Theorem 7.6).
We continue this section with the following result on the stability of the regularization index under perturbations of the potential . The proof is modeled on Lemma 2.2 in [23].
Theorem 5.4.
Assume Hypothesis 2.1 and that Hypothesis 2.4 holds at . Choose , to be principal resp. nonprincipal solutions of , satisfying and define the perturbed Sturm–Liouville differential expression
|
|
|
where the potential satisfies and
(5.2) |
|
|
|
Then will satisfy Hypothesis 2.4 and the regularization indices of and at coincide, that is, .
Proof.
We will use a similar Born series construction as for . Let and define iteratively using the Green’s function
|
|
|
|
|
|
Note that we then have for and .
Let be chosen such that do not have any zeros for . Then just as in the proof of Lemma 3.2 we can show inductively that for and . Note that we have as . It follows that converges for close enough to , and . One can check that and that is not integrable, meaning that is a principal solution for . As the regularization index only depends on the behavior of the principal (or nonprincipal) solution as if the coefficients are fixed, it follows that .
∎
Provided Hypothesis 2.4 is satisfied, we can always chose and (5.2) will hold. Note that this just corresponds to a spectral shift of . Thus we have shown
Corollary 5.5.
The regularization index is independent of the choice of in (3.4).
Remark
The perturbation condition (5.2) generalizes the condition in [23, Hypo. 2.1] for perturbed spherical Schrödinger operators. In fact, the unperturbed case treated in [23] corresponds to
|
|
|
Entire principal and nonprincipal solutions are given in terms of Bessel functions and satisfy
|
|
|
|
|
|
|
|
Thus for to be integrable we need to require and
|
|
|
This is equivalent to Hypothesis 2.1 stated in [23] and allows for the inclusion of the classical Coulomb case .
6. Relation to limit circle/limit point classification
We will now describe the relationship between the regularization index and the limit circle/limit point classification of at .
Theorem 6.1.
Assume Hypothesis 2.1 and that Hypothesis 2.4 holds at . Then is in the limit circle case at if and only if .
Proof.
First, let us assume that is in the limit circle case at (recall that Hypothesis 2.4 always holds in the limit circle case). We can then write (see (4.3))
(6.1) |
|
|
|
|
|
|
|
|
with , which exists due to the limit circle assumption. Using (3.9) with and , we conclude that for close enough to ,
|
|
|
holds. In particular for we have
|
|
|
implying , showing that indeed .
Now assume that . We will again use the above integral representation of , that is, the first line of (6.1), with chosen such that , have no zeros on . As , we have that the limit
(6.2) |
|
|
|
exists. Note that due to our assumption on , both terms in the numerator have the same sign (namely the sign of close to ). In particular, this implies
|
|
|
(due to monotonicity one can substitute for ) which shows that is in the limit circle case at .
∎
As a corollary, we have the following.
Corollary 6.2.
Assume Hypothesis 2.1, that Hypothesis 2.4 holds at and theta satisfies (4.1). Then is in the limit point case at if and only if
|
|
|
This brings up the following general open problem which can now be understood as an extension of studying stability of the limit circle/limit point classification:
Open Problem 6.3.
Under which types of perturbation of the coefficient functions is a finite regularization index of at stable?
Note that Theorem 5.4 gives one class of perturbations under which the regularization index remains stable.
For general we also have the following characterization.
Proposition 6.4.
Assume Hypothesis 2.1, that Hypothesis 2.4 holds at , and . Then
|
|
|
|
(6.3) |
|
|
|
|
In case , we have for all .
Proof.
The proof is a simple adaptation of the second part of the proof of Theorem 6.1. Let us first show that . By definition we have that (cf. (6.2))
(6.4) |
|
|
|
As is nonoscillatory for by Lemma 4.8, we can assume that is chosen such that , , and have no zeros on . Thus, both terms in the numerator of (6.4) have the same sign, so we can conclude that
|
|
|
Again, as is nonoscillatory at , it follows that
|
|
|
showing that .
Now assume that . Then we can define similarly to (3.5)
|
|
|
as the integral converges. Note that differs from by at most a multiple of . From (3.6) we thus have for close enough to that
|
|
|
implying that . Hence we must have . This argument also shows the statement regarding , thus finishing the proof.
∎
Note that the above proposition implies that due to (4.2). However, it can happen that for . For example, the generalized Bessel operator in Section 11.1 and the Jacobi equation in Section 11.2 both show that this is satisfied for roughly half of . This leads to the following open problem:
Open Problem 6.5.
If , can one characterize for what one has ?
Open Problem 6.5 is closely related to the question of equivalence between and the index studied [32], as explained in the following remark.
Remark
In [32] an index related to canonical systems is considered, which is defined through an -condition. In the Schrödinger case the is index is defined through the following procedure (see [32, Eq. (9.3)]). Choose a principal solution of and define recursively the following sequence:
|
|
|
|
|
|
|
|
It turns out that for all odd . In fact, the above recursion is equivalent the recursion (4.3) with written in two steps, and is a nonprincipal solution of . This implies that one can identify .
The authors then define the index to be the smallest such that . Thus, this index is closely related to local -integrability of the . A similar condition is used in the case and either or not in , leading to the indices , respectively, . Comparing [32, Prop. 7.18, Ex. 9.5] to our Remark 5.1 it would appear that in these cases one would have , and .
Based on these considerations, we expect that the correct answer to Open Problem 6.5 would be if and only if (cf. Cor. 10.3). This can be directly checked in the setting of Remark 5.1 where and have power like behavior at .
7. natural normalization
Note that so far there is plenty of freedom in choosing the nonprincipal solution . In fact, , with being entire and real-valued on , will still satisfy all the results from the previous section. We now introduce another stronger normalization requirement by additionally including condition below, which narrows down the class of admissible ’s. It turns out that this condition is very convenient in the study of Darboux transforms in Section 8 and guarantees that the corresponding Weyl -function lies in a suitable generalized Nevanlinna class in Section 10.
Definition 7.1 (Naturally normalized system).
Assume Hypothesis 2.1 and that Hypothesis 2.4 holds at . We call the system of entire solutions , of naturally normalized at if and only if
-
is principal at for all and for
-
for all
-
for all .
As is fixed up to multiplicative constants by , we will sometimes also say that is naturally normalized if it satisfies and , where is implicitly assumed.
Note that by Corollary 4.5, if , are naturally normalized then we must also have for . Also, if condition is vacuous. A more complete description of naturally normalized systems , is given in Theorem 7.6. Moreover, in the special case of the perturbed spherical Schrödinger operator, condition defines a ‘Frobenius solution’ in the sense of [24, Def. 3.10].
Remark We remark that it would be of interest to directly compare our notion of a naturally normalized system of solutions to those used in [32]. In particular, it is shown there that in the presence of a finite index every solution of the Sturm–Liouville problem (more generally, canonical systems) attains regularized boundary values in the sense that finitely many divergent terms are discarded in defining the values (see [32, Thms. 4.2, 7.4, and 9.6]). The regularized boundary values are then used to fix a fundamental system of solutions in order to construct a singular Weyl -function, much as we do in Section 10 via a naturally normalized system.
It will follow from the proof of Lemma 7.3 that given defined via (3.4), (3.5) (which is unique up to multiplicative constants), the natural normalization condition is equivalent to the recursion
(7.1) |
|
|
|
In fact, (7.1) can also be used to simply define for all from some initial choice of with coming from the recursion (4.3), implying that at least one naturally normalized always exists. For uniqueness, see Remark 7.3.
Lemma 7.3.
Assume Hypothesis 2.1, that Hypothesis 2.4 holds at , and let for be defined inductively by (7.1). Then there exists such that , have no zeros for and
(7.2) |
|
|
|
holds, where is defined in (3.7). In particular, all integrals (7.1) exist.
Proof.
The proof is similar to the proof of Lemma 3.2. We will show the statement via induction. For and choosing close enough to such that(3.6) holds we obtain
|
|
|
|
(7.3) |
|
|
|
|
If necessary, redefining such that the last integral above is it follows that for , showing (7.2) for . Now using induction and monotonicity of , we obtain for
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
finishing the proof.
∎
One can now proceed as in Proposition 3.3 to show that will define an entire function for each and .
Remark
Note that if and are both naturally normalized, one has
|
|
|
where is a real polynomial of degree . If were of a higher degree or a non-polynomial entire function, (7.2) could not be satisfied by both and . The appearance of the polynomial comes from the freedom to choose arbitrary real constants in (4.3), together with the fact that is also fixed only up to the addition of a constant multiple of . In particular, in the limit circle case, , our normalization condition fixes up to the addition of a constant multiple of . This is not surprising, as we just recover the usual normalization for nonoscillatory singular Sturm–Liouville operators in the limit circle case, that is, and will be normalized in the sense of [9].
Next, we want to refine Lemma 5.3 in the case when the system , is naturally normalized. For this, we will make use of the following lemma.
Lemma 7.5.
Assume , are nonoscillatory near , such that
|
|
|
exist, that is, the integrals converge. Then are nonoscillatory and if , then .
Proof.
Recall that with will have no sign changes for close enough to (see (3.6)). The same is true for by assumption. From this, it follows that are monotonic, hence nonoscillatory. Thus we can estimate
|
|
|
|
|
|
|
|
finishing the proof.
∎
We state now the refinement of Lemma 5.3 for the case that is normalized according to Definition 7.1.
Theorem 7.6.
Assume Hypothesis 2.1, that Hypothesis 2.4 holds at , and let , be naturally normalized. Then for all we have
|
|
|
and
|
|
|
Proof.
Note that by (3.8), the definition of , and (7.3) we have
|
|
|
Together with (5.1) the theorem follows by Lemma 7.5.
∎
Note that Theorem 7.6 can be summarized as stating that
(7.4) |
|
|
|
where is shorthand for .
Note that in the more general setting of Corollary 4.5, that is only assuming and the standard assumption , the three limits
(7.5) |
|
|
|
always exist in the extended real numbers (use the previous theorem and Remark 5.1 ). This will be useful in the proof of Corollary 8.5 which relies on L’Hôpital’s rule.
Before continuing with an application to Darboux transforms, we show that the property of being a naturally normalized system is independent of the choice of in Definition 7.1.
Proposition 7.7.
The property of being a naturally normalized system is independent of the choice .
Proof.
Let , be naturally normalized for some fixed choice of . Choose any and write
|
|
|
We need to show that satisfies
(7.6) |
|
|
|
To this end, observe that
|
|
|
Thus as by (3.10), we have
|
|
|
Now for we have . In fact, for and for some by Lemma 7.3. Hence, we can estimate
|
|
|
Note that in the above sum and converges as long as , which is true for close enough to . Moreover, is a monotonically increasing function (for increasing), while . Thus and the claim (7.6) follows finishing the proof.
∎
8. Connection with Darboux transforms
Let us now turn to an application of the regularization index. For this we make additional regularity assumptions on our Sturm–Liouville differential expression.
Hypothesis 8.1.
In addition to Hypothesis 2.1, assume further that and .
With these assumptions the Sturm–Liouville differential expression (2.1) can be transformed into an equivalent Schrödinger differential expression given by
(8.1) |
|
|
|
via the Liouville transform (see [4] and [13, Sect. 3.5]).
As mentioned in [13, Sect. 3.5], we point out that the conditions given in Hypothesis 8.1 allow for different examples that are not included under the typical conditions assumed, namely, the conditions with on . This can be seen by considering the elementary example on for instance.
Returning to applying the transform, choose and define
|
|
|
|
(8.2) |
|
|
|
|
|
|
|
|
|
|
|
|
Then solves , if and only if
|
|
|
where one readily verifies that under the assumptions in Hypothesis 8.1, one has where .
Note in particular that the regularization index remains invariant under the above transformation.
The motivation for introducing Hypothesis 8.1 is twofold. First, defining Darboux transforms for general Sturm–Liouville differential expressions requires additional regularity assumptions on and (see [16]). Secondly, Darboux transforms applied to differential expressions in Schrödinger form (also Liouville form) have a much simpler form (cf. [18]).
Similar to the Liouville transform, Hypothesis 8.1 is weaker than the typical assumptions for Darboux transformations such as those found in [16], allowing for more general examples to be considered.
To avoid unnecessary notation, we will assume that and , so that is already in Schrödinger form for the rest of this section.
Let us now assume that is a positive solution of , meaning that
|
|
|
Such , often called the seed function, exists if and only if is nonoscillatory at both endpoints, which we will assume from now on (see, e.g., [14, Cor. 2.4]). Then as a formal differential expression, can be factorized as follows:
|
|
|
|
Note that we avoid the common notation to emphasize that and are just formal differential expressions rather than operators. We define the associated Darboux transformed differential expression by
|
|
|
|
where
|
|
|
as can be verified by a direct computation. We say that is obtained from via a Darboux transform with seed function .
Take two functions such that , and define , . Then a quick computation shows that
|
|
|
We want to study what happens if we apply to the naturally normalized system , defined in the previous sections (cf. [26, Sect. 3]). Here we assume that Hypothesis 2.4 holds without any additional assumptions on the regularization index . We have to distinguish between two cases:
Case 1: The seed function is principal at : In this case, possibly after scaling by a real non-zero constant, we have
|
|
|
that is, , as principal solutions are unique up to scalar multiples. We define
|
|
|
Here we have used the analogue of Lemma 3.1 for . Note that . In a similar manner we define
|
|
|
Case 2: The seed function is nonprincipal at : In this case, possibly after scaling by a real non-zero constant and adding a real multiple of to it, we have
|
|
|
that is, . We then define
(8.3) |
|
|
|
In a similar manner, we define
(8.4) |
|
|
|
8.1. Properties of ,
Independent of the (non)principality of the seed function , we have the equality
(8.5) |
|
|
|
This formula follows from the identity , which can be verified by direct computation (note that the Wronskian is the same for and ).
Throughout this section we will repeatedly use L’Hôpital’s rule as summarized in the following lemma.
Lemma 8.2.
Let be given such that either or . If the limit exists in the extended real numbers, then
|
|
|
Proof.
This is a simple application of L’Hôpital’s rule:
|
|
|
∎
We now proceed with the normalization properties of and . We begin with the following proposition.
Proposition 8.3.
Assume Hypothesis 2.1, that Hypothesis 2.4 holds at , and let , be naturally normalized with either or . Then
-
, for all ,
-
, for all .
In particular, satisfies Hypothesis 2.4 at .
Proof.
First observe that . Moreover, if is a solution of , , we have
(8.6) |
|
|
|
which implies that is monotonic as , as is nonoscillatory for all . It then follows that is nonoscillatory, implying that is nonoscillatory at for all . It is also nonoscillatory for , as being nonoscillatory for implies being nonoscillatory for any .
Note the limit exists in the extended real numbers as is nonoscillatory near by (8.6). In case is nonprincipal, we obtain using L’Hôpital’s rule,
|
|
|
If is principal we instead write
|
|
|
Now the asymptotics and are equivalent by Corollary 4.5 together with (8.5), finishing the proof.
∎
The next results concern the (non)principality of , .
Proposition 8.4.
Assume Hypothesis 2.1, that Hypothesis 2.4 holds at , and let , be naturally normalized such that additionally either or , where is the seed function.
-
If the seed function is nonprincipal and is in the limit circle case at , then is nonprincipal and is principal at for all .
-
In all other cases is principal and is nonprincipal at for all
Proof.
By assumption, is in the limit circle case at so that by Theorem 6.1. We know from the proof of Proposition 8.3 that , are nonoscillatory at , meaning that the limit must exist in the extended real numbers (as ). Hence, using L’Hôpital’s rule and definitions (8.3), (8.4), as is assumed to be nonprincipal, we obtain
|
|
|
Note that , allowing us to use L’Hôpital. Hence, is principal and is nonprincipal for all by Proposition 8.3.
This case is essentially identical to the previous one with the roles of and interchanged. First let us assume that is nonprincipal, that is, . Then necessarily is in the limit point case as otherwise, we are in the previous case . Thus we have and we obtain
|
|
|
In case is principal, that is, we obtain
|
|
|
Similarly to before, it follows from Proposition 8.3 that is nonprincipal and is principal for all .
∎
As a corollary we now obtain the following.
Corollary 8.5.
Assume Hypothesis 2.1, that Hypothesis 2.4 holds at , and let , be naturally normalized with either or . Then
|
|
|
In particular, the system , is naturally normalized.
Proof.
First let us assume that is principal, that is, . Then we compute using (7.5) and Lemma 8.2,
|
|
|
and
|
|
|
Note that for we indeed have , while for we have , allowing us to use Lemma 8.2.
The case of nonprincipal, that is, can be shown analogously.
∎
9. Regularization via Darboux transforms
The goal of the present section is to demonstrate how applying a series of Darboux transforms can, in certain cases, be viewed as a type of regularization procedure. More precisely, we will say that a Schrödinger differential expression is ‘regularizable via Darboux transforms at ’ if and only if there exists a finite sequence of Darboux transforms
|
|
|
such that is in the limit circle case at .
It turns out that being regularizable via Darboux transforms at is equivalent to being finite (see Theorem 9.2 below). The key to this observation is the following theorem which shows that the regularization index is well-behaved under Darboux transforms. For this, recall that will satisfy Hypothesis 2.4 in case satisfies it (see Proposition 8.3), and thus will have a well-defined regularization index .
Theorem 9.1.
Assume Hypothesis 2.1 and that Hypothesis 2.4 holds for at the endpoint . If is the regularization index of at then the regularization index of at satisfies where we interpret as
-
if the seed function is principal at ,
-
and
|
|
|
if the seed function is nonprincipal at .
Proof.
Let us choose a naturally normalized system , and let , if is principal, and , if is nonprincipal. Then using L’Hôpital’s rule we can compute
|
|
|
where , if is principal, and , if is nonprincipal. Note that the requirements for guarantee that we are in the setting of Lemma 8.2. Thus, in case is nonprincipal, the regularization index of is given by by Definition 5.1. In case is principal (so is nonprincipal), which can only happen if is nonprincipal and , we have and . Thus , finishing the proof.
∎
We can now prove the following result stating that Hypothesis 2.4 with is necessary to transform a Schrödinger differential expression into one which is in the limit circle case at through a finite number of Darboux transforms. This then gives us a complete characterization of those Schrödinger differential expressions which can be regularized via Darboux transforms.
Theorem 9.2.
Assume Hypothesis 2.1 and let be a Schrödinger differential expression which is nonoscillatory at both endpoints. Then can be transformed via a finite series of Darboux transforms to a Schrödinger differential expression which is in the limit circle case at if and only if Hypothesis 2.4 holds at and . In this case, the minimal number of Darboux transforms is and is achieved if the seed functions are always chosen to be nonprincipal at .
Proof.
First, let us remark that as is assumed to be nonoscillatory at both endpoints, nonvanishing seed functions can always be found. Indeed, as we will see in (10.2), one can explicitly write down such seed functions which are nonprincipal at both endpoints.
Hence, it follows from Theorem 9.1 that if there exists a sequence of many Darboux transforms leading to a which is in the limit circle case at . It is also immediate from Theorem 9.1 that no smaller number of Darboux transforms will achieve this.
It remains to show that if Hypothesis 2.4 does not hold for , then it also does not hold for its Darboux transform .
Let us now assume that the seed function is principal (the other case is proven in a similar manner). Then using L’Hôpital’s rule we obtain
|
|
|
Now it follows from Theorem 4.2 that if Hypothesis 2.4 does not hold for , it will also not hold for . By induction, the same is true for any obtained from through a finite series of Darboux transform, implying that must remain in the limit point case at . This finishes the proof.
∎
As a simple corollary we also obtain.
Corollary 9.3.
A Schrödinger differential expression which is nonoscillatory at both endpoints can be transformed via Darboux transforms to one which is in the limit circle case at both endpoints if and only if Hypothesis 2.4 holds and . In this case the minimal number of Darboux transforms is and is achieved by choosing seed functions which are nonprincipal at both endpoints.
In other words, is regularizable via Darboux transforms at both endpoints if and only if .
We should remark that while these results are proven for Schrödinger differential expressions, analogous statements can be made for general Sturm–Liouville differential expressions satisfying Hypothesis 8.1 through the use of the Liouville transform which leaves the regularization index invariant.
To justify the terminology ‘regularization via Darboux transforms’ we recall that Darboux transforms as above can change the spectrum only at the value of the spectral parameter of the seed function, provided correct boundary conditions are specified (see [2]). In particular, self-adjoint realizations of and self-adjoint realizations of its Darboux transform will have in general similar spectral properties. Thus we expect problems having finite regularization indices to behave similarly to regular problems. As an example, we say that satisfies Weyl asymptotics if and only if every self-adjoint realization of has a discrete spectrum such that the eigenvalues additionally satisfy the asymptotics
(9.1) |
|
|
|
Here it is implicitly assumed that the integral above is finite.
It is known that satisfies Weyl asymptotics if it is regular, or more generally limit circle nonoscillatory at both endpoints (see [35, Remark 3.1, Lem. 3.5(3)]). Hence, we obtain from Corollary 9.3 an elementary proof of the following fact.
Corollary 9.4.
Assume Hypotheses 2.1 and 8.1. If Hypothesis 2.4 is satisfied at with , then the Weyl asymptotics (9.1) hold.
Proof.
Note that given satisfying Hypothesis 8.1 we can transform it to an equivalent differential expression
in Schrödinger form via the Liouville transform. A similar isospectral transformation with holds for self-adjoint realizations of (see, e.g., [6, Sect. 3.2] for the regular case). From Prop. 10.5, and as the interval remains invariant under Darboux transforms, it follows that the eigenvalues of (hence of ), satisfy , in particular must be finite. But according to (8.2). This shows (9.1) finishing the proof.
∎
We remark that the above result is not new. In fact [37, p. 33] implies that, provided the regularization indices are finite, Weyl asymptotics hold even without Hypothesis 8.1. This result is based on [38, Thm. 4.8], which appears to use very different techniques compared to the present paper.
Proposition 9.5.
([38, Thm. 4.8], [37])
Assume Hypothesis 2.1. If Hypothesis 2.4 is satisfied at with , then the Weyl asymptotics (9.1) hold.
Proof.
This result follows from [37, p. 33] since implies used in [32] is also finite (see Remark 6.5 for instance).
∎
As our regularization procedure relies on the Liouville and Darboux transforms, Hypothesis 8.1 was necessary. This raises the following question:
Open Problem 9.6.
Assume Hypothesis 2.1. If Hypothesis 8.1 does not hold, but Hypothesis 2.4 is satisfied at with and at least one index positive, can the problem still be regularized in the sense of transforming into an associated regular problem using a different method?
Regarding Weyl asymptotics, our work relies heavily on the finiteness of the regularizations indices. However, having an infinite regularization index is still compatible with Weyl asymptotics see Section 11.3 for an example. This leads us to the following problem:
Open Problem 9.7.
Can one characterize when Weyl asymptotics will hold based on the behavior of the system for problems with infinite index?
We point out that Sections 11.3 and 11.4 show that for an infinite regularization index, Weyl-asymptotics may or may not hold. In such cases we cannot rely on Darboux transforms, so we instead use the following lemma to prove Weyl asymptotics for the example in Section 11.3 which has an infinite regularization index.
Lemma 9.8.
Let be finite and assume that satisfies Weyl asymptotics. If the potential is bounded from below, then satisfies Weyl asymptotics.
Proof.
Fix and , and assume w.l.o.g that .
Let be the Dirichlet eigenvalues of for small enough. Then these eigenvalues will satisfy Weyl asymptotics, as , as the problem is regular at both endpoints. Similarly, the Dirichlet eigenvalues of also satisfy the same asymptotics.
By the Sturm–Picone comparison theorem we must have that for all . Moreover [39, Thm. 4.4.4] (and the remark after it) shows that the Dirichlet eigenvalues and increase as a function of , while [39, Thm. 10.8.2] shows that resp. , where and are the eigenvalues of the Friedrichs realization of and , respectively. Note that as has Weyl asymptotics by assumption, we have as . Thus it follows that for all , proving Weyl asymtptotics for .
∎
10. Weyl -function in case of finite regularization index
In this section we compute Weyl -functions in the case of a finite regularization index. As it turns out, these -functions will be in the generalized Nevanlinna–Herglotz class with . More explicitly, the underlying -functions can be written in terms of -functions of limit circle problems having the familiar Nevanlinna–Herglotz property (see Propositions 10.1 and 10.5).
Again to simplify our analysis, we assume that is in the Schrödinger form (8.1) with . As we are interested in at least one endpoint being in the limit point case, we will exclude the case of both endpoints being limit circle and without loss of generality assume with .
Assuming momentarily that , we denote by the unique self-adjoint realization of (as both endpoints are in the limit point case). Consider now an arbitrary naturally normalized system , at the endpoint (and analogously , at the endpoint ). Then we can define the singular Weyl -function (see [17], [25]) satisfying the equation
(10.1) |
|
|
|
where is some holomorphic function on .
It is known that the spectrum of can be recovered from the limits , . In particular, in the presence of a purely discrete spectrum, can be extended to a meromorphic function with poles at the eigenvalues of .
As described in the previous section, we will ‘regularize’ the expression by applying a sequence of Darboux transforms which lower the regularization index to zero so that is a limit circle endpoint while simultaneously lowering . While this process will not be unique, we will choose the optimal method in the sense of the least number of transforms while also adding the fewest eigenvalues possible. To this end we are looking for a seed function solving satisfying the following properties:
-
for all
-
is nonprincipal at both endpoints.
We choose an appropriate seed function satisfying conditions , above as follows. Consider any (note that is bounded from below as is nonoscillatory at both endpoints). Then is nonprincipal at , as otherwise would be an eigenvalue of (recall that Hypothesis 2.4 implies that principal solutions are -integrable near the endpoint in question). Moreover, it follows from [14, Cor. 2.4] that has a fixed sign on . We now define the seed function as
(10.2) |
|
|
|
It is easy to see that conditions , are satisfied. Moreover, we add a constant multiple of to in order to also have
|
|
|
Introducing the notation , we arrive at the equation
(10.3) |
|
|
|
This corresponds to changing in (10.1) by an additive constant.
Applying to both sides of (10.3) yields
|
|
|
Note that is the -function of the self-adjoint realization of the Darboux transformed , with Dirichlet boundary conditions (if any) since and are principal at and , respectively. Moreover, as is nonprincipal at , it follows that , implying that , that is, one eigenvalue was added to the spectrum. Also, by Theorem 9.1, we have that both indices have been lowered by 1.
Now, since we assumed , we can repeat this procedure times as follows: we introduce the notation and shift by a constant after choosing at each step through
|
|
|
|
|
|
Hence, after steps one arrives at
|
|
|
(10.4) |
|
|
|
with
(10.5) |
|
|
|
where denotes the eigenvalue added (below the previous step’s spectrum) at the th step with and the product is empty and equal to for . If , then the products in (10.5) are considered empty and equal to 1 (i.e., no Darboux transformation was needed to make the endpoint limit circle so (10) and (10.1) agree). As and remain principal at and , respectively (see Proposition 8.4), is the -function of the self-adjoint realization of with necessarily Dirichlet boundary conditions at (as ) and either Dirichlet at if or no boundary conditions at if (as ). That is, is the Friedrichs realization of (see [9, Thm. 4.7]).
As it is known that when is a limit circle endpoint, the Weyl -function defined via (10.5) with a naturally normalized system of solutions is in fact a Nevanlinna–Herglotz function (also called a Pick function; see [9, Eq. (5.12)]), we arrive at the following result (see Remark 10.3 for weaker conditions for the theorem to hold):
Proposition 10.1.
Assume Hypotheses 2.1 and 8.1. Furthermore, assume Hypothesis 2.4 is satisfied at with , let , and denote by either the unique self-adjoint realization of if or the Friedrichs realization otherwise. Then for each choice of real numbers satisfying there is a Nevanlinna–Herglotz function and constants , such that the Weyl -function in (10.1) can be written as
|
|
|
If , the sum and products are understood as empty and equal to 0 and 1, respectively.
Moreover, can be understood as the -function for the Friedrichs realization, , of an -times Darboux transformed so that .
Proof.
Follows from the previous discussion (after a possible Liouville transform) by using the naturally normalized system of solutions at the endpoint corresponding to on the left-hand side of (10.1). Note that the solutions will still be naturally normalized after the application of Darboux transforms by Corollary 8.5.
∎
Remark
Due to Remark 7.3, a naturally normalized system at is essentially unique up to the addition of to , where is an arbitrary real polynomial of degree . Under such a transformation the -function in (10.1) would become . In particular, we can always normalize such that the new -function has the simpler form
(10.6) |
|
|
|
Notice that one implication of Proposition 10.1 is that the -function satisfying (10.1) for any naturally normalized system is not a Nevanlinna–Herglotz function unless . In fact, as a corollary we see that under the assumptions of the previous theorem, the -function of the original problem is a type of generalized Nevanlinna–Herglotz function.
Corollary 10.3.
In addition to the assumptions and notation of Proposition 10.1, let
Then any singular -function coming from (10.1) with , are naturally normalized is in the subclass consisting of generalized Nevanlinna–Herglotz functions with negative squares, no nonreal poles, and the only generalized pole of nonpositive type at infinity.
Proof.
We first point out that is the class of Nevanlinna–Herglotz functions. Also, by Remark 10.1 the -function can be written as , with of the form (10.6), and a real polynomial of degree . Furthermore, the -function in (10.5) is a Nevanlinna–Herglotz function satisfying
(10.7) |
|
|
|
as this holds for any Weyl -function for Sturm–Liouville operators with one limit circle endpoint. In particular, it has a Nevanlinna–Herglotz representation (see [13, Eq. (6.2.39)]).
Moreover, by performing one more Darboux transformation with nonprincipal seed function, one would multiply by another simple pole and still have a Nevanlinna–Herglotz function (see 10). Thus, cannot decay and must grow sublinearly by (10.7). In fact, it must grow since if it had a finite limit at infinity, one could redefine in (10) by adding this limit times to reach a contradiction since the new -function would decay.
We suspect this is well-known behavior for the Friedrichs -function of this form, but could not find an explicit statement in the literature.
The generalized Nevanlinna–Herglotz property for of the form (10.6) now follows from repeated applications of [26, Cor. 3.6] taking in [26, Eq. (3.14)] at each step and understanding the infinite limit in [26, Cor. 3.6] as unbounded growth of the ratio. In particular, the first and then every other multiplication by raises the index in by one. As the class is invariant under addition of a real polynomial of degree (see [8, p. 190]), the claim follows for the general .
∎
A few remarks are now in order.
Remark
Proposition 10.1 and Corollary 10.3 remain true by assuming Hypothesis 2.4 holds with a finite regularization index at only one endpoint, and the other endpoint is nonoscillatory, with no further restrictions.
We remark that since the seed function was chosen to be nonprincipal at both endpoints at each step, this is the optimal method to regularize the endpoint with smallest regularization index in the sense of the least number of transforms and adding the fewest eigenvalues by Theorem 9.2.
The definition of the regularization index at is closely related to the angular momentum of a perturbed Bessel operator studied in the series of papers [23]–[27]. More precisely, , that is, in our language the perturbed Bessel operator studied in [23]–[27] has a regularization index at .
So far we applied Darboux transforms until one of the endpoints (in our case ) is in the limit circle case. If the regularization indices are equal, then the Darboux transformed expression is in the limit circle nonoscillatory case at both endpoints. Otherwise, assuming , we now want to continue the procedure -times choosing seed functions with nonprincipal behavior at both endpoints as before. Care will now be needed as the principal/nonprincipal behavior of the naturally normalized system at the endpoint will swap with every further transform by Proposition 8.4.
Assume now that , and that we have completed Darboux transforms to arrive at (10.5). As and are still principal at and , respectively, the next Darboux transformation is exactly the same as previously. Hence, preceding as before, choosing , yields
|
|
|
(10.8) |
|
|
|
where
|
|
|
However, will now be nonprincipal at , while will be principal at , by Proposition 8.4. Note that will remain principal at (which will continue to be the case). This time is the -function of the self-adjoint realization of with the Neumann-type boundary condition at defined via and either Dirichlet at if or no boundary conditions at if . While this step once again adds an eigenvalue, we will have to modify the next seed function chosen since the principal/nonprincipal behavior near interchanged.
Assuming now that (so that at least one more step is needed in the regularization process), we choose the seed function with ,
|
|
|
which is nonprincipal at both endpoints (as the nonprincipality at remained unchanged through each transformation). The main difference now is that no constant shift is needed for the choice of seed function above. Therefore we need only apply to (10) on this step to arrive at
|
|
|
(10.9) |
|
|
|
where
|
|
|
(10.10) |
|
|
|
In particular, (10) shows that we did not add an eigenvalue in this step and instead added a zero to the -function (since was not an eigenvalue of the previous step). Moreover, the principal/nonprincipal behavior at once again swaps, so is the -function for the corresponding Friedrichs realization once again (with Dirichlet boundary conditions at if and only if ).
In case we have to apply another Darboux transform (from then on the pattern will repeat). Similarly to before, as the endpoint is in the limit circle case, defines a Neumann-type boundary condition at this endpoint. In fact, is the lowest eigenvalue of the self-adjoint realization of with the -boundary condition at . This mean that we can apply [14, Cor. 2.4] to conclude that for any we have that
|
|
|
is nonvanishing and nonprincipal at the endpoint (as otherwise would be a smaller eigenvalue). In particular, no constant shift in the -function is required. We can apply to both sides of (10)
which adds an eigenvalue and swaps the principal/nonprincipal behavior at .
Finally, this process can be iterated times to arrive at the case of both regularization indices becoming zero and the -function
|
|
|
|
(10.11) |
|
|
|
|
We remark that the additional poles and zeros added after the th step in (10) will necessarily interlace by construction, and a total of eigenvalues were added during the regularization process.
We now summarize this regularization process by extending Proposition 10.1:
Proposition 10.5.
Assume Hypotheses 2.1 and 8.1. Furthermore, assume Hypothesis 2.4 is satisfied at with , let , , and denote by either the unique self-adjoint realization of if or the Friedrichs realization otherwise. Then for each choice of real numbers satisfying there is a Nevanlinna–Herglotz function and constants , such that the Weyl -function in (10.1) can be written as
|
|
|
(10.12) |
|
|
|
Moreover, if is even, then can be understood as the -function for the Friedrichs realization, , of a -times Darboux transformed quasi-regular . If is odd, then can be understood as the -function for a self-adjoint realization with Dirichlet boundary condition at the endpoint with larger index, and a Neumann-type boundary condition at the other endpoint.
In both cases, where so that the number of eigenvalues added during the regularization process is .
Proof.
Follows from the previous discussion combined with Prop. 10.1.
∎
As both endpoints of are now in the limit circle nonoscillatory case, a corresponding regular expression can be found by utilizing [39, Thm. 8.3.1], effectively regularizing . Once again, since the seed function was chosen to be nonprincipal at both endpoints at each step, this is the optimal method to regularize both endpoints in the sense of the least number of transforms and adding the fewest eigenvalues by Corollary 9.3. We also point out that the analog of Remark 10.1 is true in this more general case as well, that is, the sum in (10.5) can be removed.
We end by posing the following question:
Open Problem 10.6.
If both are infinite or one is infinite and the other is not defined but the problem is still nonoscillatory for all , is the -function corresponding to naturally normalized systems in some larger special class of functions than generalized Nevanlinna–Herglotz functions?
Given the relation , we would expect the -functions when to display superpolynomial growth. Moreover, the corresponding special class of functions would have to be invariant under the addition of arbitrary entire functions which are real on the real line by Remark 7.3.