🔢 mathematics

Lower bounds for self-adjoint Sturm-Liouville operators

This paper provides estimates for the lower bound of the self-adjoint operator associated with a three-coefficient Sturm-Liouville differential expression within the weighted $L^2$ -Hilbert space $L^2(\mathbb R; rdx)$ .

Original authors: Jussi Behrndt, Fritz Gesztesy, Philipp Schmitz, Carsten Trunk

Published 2026-10-21

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Jussi Behrndt, Fritz Gesztesy, Philipp Schmitz, Carsten Trunk

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

In the vast landscape of mathematical physics, there is a class of problems that act as the fundamental grammar for understanding how energy behaves in one-dimensional systems. Imagine a long, thin wire or a single path through space where a particle, like an electron, is free to move but is constantly nudged by an invisible landscape of forces. Mathematicians describe this landscape using a specific type of equation that balances the particle's tendency to spread out against the pull of these forces. For decades, scientists have been able to predict the lowest possible energy state of such a system when the path is uniform and the forces are simple. This lowest energy state is crucial because it determines whether the system is stable or if it might collapse into chaos. If the energy can drop too low, the system becomes unpredictable; if there is a firm floor beneath the energy, the system remains safe and manageable.

However, the real world is rarely uniform. The path itself might change thickness, or the medium through which the particle travels might become denser or thinner in unpredictable ways. When these variations are introduced, the simple rules that worked for the uniform case break down. The question then becomes: even with these messy, changing conditions, can we still guarantee that the energy of the system has a bottom limit? Can we find a number that says, "No matter how the path twists or how the forces fluctuate, the energy will never fall below this point"? This is the central puzzle that a team of mathematicians set out to solve. They focused on a complex version of these equations where the path's properties and the forces acting on it are allowed to vary wildly, provided they follow certain basic rules of continuity.

The researchers, working with a general mathematical expression that describes these varying conditions, developed a new method to calculate this safety floor. Instead of trying to solve the equation for every specific scenario, which would be impossible given the infinite variety of ways the path could change, they looked for a universal rule. They treated the changing path and the forces as ingredients that could be measured in specific ways. By carefully weighing the "roughness" of the path and the strength of the negative forces that try to pull the energy down, they derived a formula that acts as a safety net. This formula does not tell you the exact energy of a specific system, but it provides a guaranteed lower bound. It tells you that as long as the variations in the path and the forces stay within certain measurable limits, the energy cannot plunge into an infinite abyss.

A key part of their discovery involves how they handle the weight of the path itself. In some regions, the path might become so thin or light that it threatens to destabilize the whole system. The team realized that to find a safe lower limit, they had to account for these weak spots. They introduced a way to measure the size of the regions where the path is particularly light and combined this with the strength of the forces in those areas. If the weak spots are too large or the forces too strong in those specific zones, the safety floor drops. But their method allows them to calculate exactly how low it can go before the system loses its stability. They showed that even if the path becomes extremely thin in some places, as long as the forces do not become infinitely strong in a way that overwhelms the path's structure, a limit exists.

The team proved that their estimates hold true under a wide range of conditions. They demonstrated that if the forces trying to pull the energy down are not too intense relative to the stiffness of the path, the system remains stable. They provided specific calculations for different types of variations, showing how the lower bound changes depending on whether the variations are spread out smoothly or concentrated in sharp spikes. In one particularly clear case, they showed that if the path is uniform and the forces are bounded, the lower limit is directly related to the strength of those forces. But in the more complex, general case, the limit depends on a delicate balance between the size of the weak regions and the intensity of the forces within them.

Their work is a rigorous proof, meaning it is a logical certainty derived from the rules of mathematics, not a guess based on computer simulations or experimental data. They did not just suggest that a limit exists; they constructed a mathematical argument that guarantees it. They also identified a specific condition where the system is guaranteed to be completely stable, with an energy floor of zero or higher. This happens when the product of the total "weakness" of the path and the total strength of the negative forces is small enough. If this product stays below a certain threshold, the system is safe from collapsing.

This research fills a gap in our understanding of these mathematical models. While scientists have long known how to handle simple, uniform cases, this paper provides the first systematic way to handle the messy, variable cases that are more common in real-world applications. It offers a toolkit for physicists and engineers who work with systems where materials are not perfect and conditions are not constant. By providing these estimates, the authors have given a way to predict the stability of complex systems without needing to solve the impossible task of calculating every single detail. The result is a clearer picture of the boundaries of stability in a changing world, ensuring that even in the most variable conditions, there is always a floor beneath the energy.

Technical Summary: Lower Bounds for Self-Adjoint Sturm–Liouville Operators

Problem Statement
The paper addresses the problem of deriving explicit estimates for the lower bound of the spectrum, $\inf(\sigma(T))$ , for a general self-adjoint Sturm–Liouville operator $T$ . This operator is associated with the three-coefficient differential expression:
$\tau = \frac{1}{r} \left( -\frac{d}{dx} p \frac{d}{dx} + q \right)$
acting in the weighted $L^2$ -Hilbert space $L^2(\mathbb{R}; r dx)$ . The coefficients $p, q, r$ are real-valued, with $r$ and $p$ positive almost everywhere and locally integrable (specifically $r, 1/p, q \in L^1_{\text{loc}}(\mathbb{R})$ ).

While lower bounds for the special case of the Schrödinger operator ( $p=r=1$ ) are well-established in the literature (e.g., via Lieb–Thirring inequalities or specific estimates for $q \in L^s$ ), a systematic study for the general three-coefficient case has been lacking. The authors note that the weight function $r$ and the coefficient $p$ fundamentally influence the lower bound. Specifically, if $r$ is non-constant or unbounded, the $L^\infty$ -norm of $1/r$ (or a related function neutralizing small values of $r$ ) and the norm of $1/p$ become critical components in the estimates.

Methodology
The authors employ a variational approach combined with functional analytic inequalities to bound the quadratic form associated with $T$ . The methodology proceeds as follows:

Reduction to Negative Spectrum: The analysis is restricted to the set $D_-(T) = \{f \in D(T) \mid (Tf, f)_r \leq 0\}$ . If this set is trivial, the operator is non-negative.
Decomposition of Potential: The potential $q$ is decomposed into its positive and negative parts, $q = q_+ - q_-$ . The lower bound is determined by the behavior of the negative part $q_-$ .
Integration by Parts and Lemma 4.2: Using integration by parts and boundary conditions derived from the limit-point case at $\pm \infty$ , the authors establish that for $f \in D_-(T)$ , the kinetic energy term $\|\sqrt{p}f'\|_2^2$ is controlled by the potential term $\|q_- f^2\|_1$ .
Interpolation and Norm Estimates: The core of the technical contribution involves deriving $L^\infty$ $L^{\infty}$ bounds for eigenfunctions (or test functions) in terms of their $L^2$ $L^{2}$ norms and the coefficients $p$ $p$ and $q_-$ $q_{-}$ . The authors utilize:
- Hölder's Inequality: To handle the interaction between $q_-$ and $f^2$ in various $L^s$ spaces.
- Sobolev-type Inequalities: Specifically, estimates for $\|f\|_\infty$ derived in [2] (Lemma 4.4), which relate the $L^\infty$ norm to $\|\sqrt{p}f'\|_2$ and $\|f\|_2$ under assumptions on $1/p$ (either $L^\infty$ or $L^\eta$ ).
- Young's Inequality: Used in the warm-up section for the Schrödinger case to derive explicit constants.
Handling the Weight Function: A key technical device is the introduction of a non-negative function $g \in L^\infty(\mathbb{R})$ and the set $\Omega_g = \{x \in \mathbb{R} \mid r(x)g(x) < 1\}$ . The measure $\mu(\Omega_g)$ appears in the denominator of the final estimates, allowing the authors to manage cases where $1/r$ is not in $L^\infty(\mathbb{R})$ .

Key Contributions and Results
The paper provides a series of theorems (3.2–3.5) offering lower bound estimates under different integrability assumptions on the coefficients.

Theorem 3.2: Assumes $1/p \in L^\infty(\mathbb{R})$ and $q_- \in L^1_u(\mathbb{R})$ (uniformly locally integrable). It provides a lower bound involving $\|q_-\|_u$ , $\|1/p\|_\infty$ , and a function $g$ such that $\mu(\Omega_g)\beta < 1$ .
Theorem 3.3: Assumes $1/p \in L^\infty(\mathbb{R})$ and $q_- \in L^s(\mathbb{R})$ for $s \in [1, \infty]$ . This yields explicit constants depending on $\|q_-\|_s$ and $\|1/p\|_\infty$ . The authors note that while these bounds are not optimal (comparing them to the known Schrödinger case $p=r=1$ ), they provide a systematic extension to the general setting.
Theorem 3.5: Generalizes the result to the case where $1/p \in L^\eta(\mathbb{R})$ for $\eta \in [1, \infty)$ and $q_- \in L^s(\mathbb{R})$ , provided $\eta + s > 2$ . This allows for coefficients $p$ that are not bounded away from zero in the $L^\infty$ sense, provided they satisfy specific integrability conditions.
Proposition 3.6: Provides a sufficient condition for the non-negativity of $T$ : if $1/p, q_- \in L^1(\mathbb{R})$ and $\|1/p\|_1 \|q_-\|_1 < 1$ , then $\inf(\sigma(T)) \geq 0$ .

In the special case where $1/r \in L^\infty(\mathbb{R})$ , the authors show that one can choose $g = 1/r$ , making $\mu(\Omega_g) = 0$ , which simplifies the estimates to explicit formulas involving $\|1/r\|_\infty$ .

Significance and Claims
The authors state that the main objective is to "initiate and contribute to this circle of problems" regarding lower bounds for general Sturm–Liouville operators, as a systematic study was previously missing.

Motivation: The work is motivated by the importance of spectral estimates in quantum mechanics and the specific influence of the weight $r$ and coefficient $p$ on the spectrum, which cannot be simply scaled as in the constant coefficient case.
Methodological Inspiration: The proofs are strongly inspired by the techniques in [2], which dealt with bounds on nonreal eigenvalues of indefinite Sturm–Liouville operators.
Scope: The paper focuses on deriving these estimates and does not claim to provide optimal constants (as evidenced by Remark 3.4, which compares their constants to the known optimal ones for the Schrödinger case).
Applicability: The results are established for the whole line $\mathbb{R}$ , but Remark 3.7 notes that the estimates remain valid for the half-line $(0, \infty)$ with Dirichlet boundary conditions at the regular endpoint.

The paper concludes by acknowledging that while the bounds are not necessarily sharp, they provide the first systematic framework for estimating the lower bound of the spectrum for the general three-coefficient Sturm–Liouville operator in weighted spaces.

More like this