Stieltjes Spectral Parameter Space

• Stieltjes Spectral Parameter Space is a framework that uses positive measures and Stieltjes transforms to encode spectral information in fields such as causal nonlocal gravity and operator theory.
• It encapsulates key parameters like the characteristic spectral scale (μ_char) and total spectral weight (M₀), which govern wave propagation, dispersive corrections, and short-range gravity modifications.
• The framework unifies methodologies from Bayesian ringdown analyses, modified dispersion relations, advection–diffusion homogenization, and random-matrix spectral detection.

Searching arXiv for papers directly relevant to “Stieltjes spectral parameter space” and closely related Stieltjes spectral frameworks. “Stieltjes spectral parameter space” most specifically denotes the $(\mu_{\rm char},M_0)$ plane introduced for causal nonlocal gravity models whose nonlocality is encoded by a positive spectral density $\rho(\mu)$ in a retarded Stieltjes-type kernel (Balfagon, 17 Mar 2026). In a broader mathematical and physical literature, closely related Stieltjes spectral frameworks appear whenever a positive measure, a self-adjoint operator, or a resolvent-type analytic function organizes the dependence of an observable on a spectral variable. In that wider sense, the phrase refers not to a single universal construction but to a family of parameterizations built from Stieltjes transforms, spectral measures, and positivity-preserving analytic structures (Murphy et al., 2024).

1. Definition in causal nonlocal gravity

In the gravity setting, the nonlocal operator is defined by the retarded Stieltjes-type kernel

$$\mathcal{K}^{-1} = \int_0^\infty \rho(\mu)\, (-\Box_g + \mu)^{-1}_R\,d\mu, \qquad \rho(\mu)\geq 0.$$

Here $(-\Box_g+\mu)^{-1}_R$ is the retarded Green operator for a massive Klein–Gordon operator on a curved spacetime, and the full kernel is a positive superposition of such retarded propagators. The positivity condition $\rho(\mu)\ge 0$ is stated to preserve causality, unitarity, and absence of ghost-like excitations at linear order; the appendix makes the causal statement explicit through

$$\mathrm{supp}\,\mathcal{K}^{-1}f \subset J^+(\mathrm{supp}\,f).$$

The associated transfer function is a completely monotone Stieltjes function satisfying

$$(-1)^n \frac{d^n m}{d(\omega^2)^n} \geq 0 \quad \forall\, n\in\mathbb{N}_0,\;\omega^2>0$$

(Balfagon, 17 Mar 2026).

The parameter space itself is built from moments of $\rho(\mu)$. The total spectral weight is

$M_0 = \int_0^\infty \rho(\mu)\,d\mu,$

and inverse moments are defined by

$$\mu_{\rm eff}^{-n} = M_0^{-1}\int\rho(\mu)\,\mu^{-n}\,d\mu.$$

For observational discussions, the paper uses a characteristic spectral scale $\rho(\mu)$0, appropriate when $\rho(\mu)$1 is concentrated around one mass scale, so that

$\rho(\mu)$2

This gives the geometric interpretation of the plane: the horizontal axis records where the spectral density sits, and the vertical axis records how much total spectral weight the kernel carries (Balfagon, 17 Mar 2026).

The same framework introduces the causal scale $\rho(\mu)$3 through

$\rho(\mu)$4

and for a sharply peaked spectrum,

$\rho(\mu)$5

The fiducial range is

$\rho(\mu)$6

This places the intended causal nonlocal regime at short distances rather than in the infrared (Balfagon, 17 Mar 2026).

2. Kernel, transfer function, and observable meaning

The basic frequency-domain object is the Stieltjes transfer function

$\rho(\mu)$7

This is the spectral representation that carries the information in $\rho(\mu)$8 into wave propagation. When the gravitational-wave frequency lies well below the spectral support, $\rho(\mu)$9, the low-frequency behavior is suppressed by the characteristic scale of the measure, so large $$\mathcal{K}^{-1} = \int_0^\infty \rho(\mu)\, (-\Box_g + \mu)^{-1}_R\,d\mu, \qquad \rho(\mu)\geq 0.$$0 implies small low-energy effects (Balfagon, 17 Mar 2026).

The leading dispersive correction is written as

$$\mathcal{K}^{-1} = \int_0^\infty \rho(\mu)\, (-\Box_g + \mu)^{-1}_R\,d\mu, \qquad \rho(\mu)\geq 0.$$1

and, in the low-frequency regime,

$$\mathcal{K}^{-1} = \int_0^\infty \rho(\mu)\, (-\Box_g + \mu)^{-1}_R\,d\mu, \qquad \rho(\mu)\geq 0.$$2

This formula is the direct map from a spectral density to a propagation constraint. It shows that $$\mathcal{K}^{-1} = \int_0^\infty \rho(\mu)\, (-\Box_g + \mu)^{-1}_R\,d\mu, \qquad \rho(\mu)\geq 0.$$3 controls the overall amplitude of the modification, while $$\mathcal{K}^{-1} = \int_0^\infty \rho(\mu)\, (-\Box_g + \mu)^{-1}_R\,d\mu, \qquad \rho(\mu)\geq 0.$$4 or $$\mathcal{K}^{-1} = \int_0^\infty \rho(\mu)\, (-\Box_g + \mu)^{-1}_R\,d\mu, \qquad \rho(\mu)\geq 0.$$5 sets the decoupling scale (Balfagon, 17 Mar 2026).

The static limit is described by a Yukawa-like potential,

$$\mathcal{K}^{-1} = \int_0^\infty \rho(\mu)\, (-\Box_g + \mu)^{-1}_R\,d\mu, \qquad \rho(\mu)\geq 0.$$6

and in the general spectral form,

$$\mathcal{K}^{-1} = \int_0^\infty \rho(\mu)\, (-\Box_g + \mu)^{-1}_R\,d\mu, \qquad \rho(\mu)\geq 0.$$7

For a delta-like spectrum $$\mathcal{K}^{-1} = \int_0^\infty \rho(\mu)\, (-\Box_g + \mu)^{-1}_R\,d\mu, \qquad \rho(\mu)\geq 0.$$8, this reduces to the standard Yukawa form. The Stieltjes parameter space therefore unifies dispersive wave propagation and short-range static modifications within a single measure-theoretic description (Balfagon, 17 Mar 2026).

3. Observational mapping and excluded regions

The gravity paper uses three observational inputs. First, a Bayesian ringdown analysis of $$\mathcal{K}^{-1} = \int_0^\infty \rho(\mu)\, (-\Box_g + \mu)^{-1}_R\,d\mu, \qquad \rho(\mu)\geq 0.$$9 BBH events in GWTC-3 constrains a universal QNM deformation parameter to

$(-\Box_g+\mu)^{-1}_R$0

so effectively

$(-\Box_g+\mu)^{-1}_R$1

The cumulative log Bayes factor is reported as $(-\Box_g+\mu)^{-1}_R$2. Within the Stieltjes framework, the deformation scales roughly as

$(-\Box_g+\mu)^{-1}_R$3

and for the fiducial scale the estimate is

$(-\Box_g+\mu)^{-1}_R$4

Accordingly, ringdown excludes models producing percent-level deviations but does not constrain the fiducial Stieltjes regime (Balfagon, 17 Mar 2026).

Second, the paper maps published GWTC-3 modified dispersion relation bounds into the spectral plane. The LVK phenomenology uses

$(-\Box_g+\mu)^{-1}_R$5

The graviton-mass bound

$(-\Box_g+\mu)^{-1}_R$6

implies

$(-\Box_g+\mu)^{-1}_R$7

which the paper describes as far below anything relevant to the fiducial model (Balfagon, 17 Mar 2026).

Third, the strongest propagation constraint comes from GW170817: $(-\Box_g+\mu)^{-1}_R$8 Using $(-\Box_g+\mu)^{-1}_R$9 Hz and the group-velocity relation, the paper states that this excludes spectral densities producing measurable dispersion whenever

$\rho(\mu)\ge 0$0

The resulting hierarchy is

$\rho(\mu)\ge 0$1

The excluded region is therefore the infrared part of the spectrum: spectral densities with support below $\rho(\mu)\ge 0$2, power-law tails $\rho(\mu)\ge 0$3 with $\rho(\mu)\ge 0$4 extending into that infrared regime, and any spectrum producing dispersion larger than the GW170817 bound in the LIGO band. By contrast,

$\rho(\mu)\ge 0$5

is consistent with current GW constraints (Balfagon, 17 Mar 2026).

Sub-millimetre gravity supplies the strongest direct probe at the fiducial scale. Because

$\rho(\mu)\ge 0$6

the natural window is short-range gravity. Using Eöt-Wash bounds at $\rho(\mu)\ge 0$7,

$\rho(\mu)\ge 0$8

the paper translates the constraint to

$\rho(\mu)\ge 0$9

It describes this as the strongest existing constraint on the Stieltjes kernel at the fiducial scale. The projected gravitational-wave sensitivities remain much weaker than the fiducial prediction: LIGO O5 at $$\mathrm{supp}\,\mathcal{K}^{-1}f \subset J^+(\mathrm{supp}\,f).$$0, Einstein Telescope / Cosmic Explorer at $$\mathrm{supp}\,\mathcal{K}^{-1}f \subset J^+(\mathrm{supp}\,f).$$1, and LISA at $$\mathrm{supp}\,\mathcal{K}^{-1}f \subset J^+(\mathrm{supp}\,f).$$2, versus $$\mathrm{supp}\,\mathcal{K}^{-1}f \subset J^+(\mathrm{supp}\,f).$$3 for the fiducial nonlocal effect (Balfagon, 17 Mar 2026).

4. Spectral-measure frameworks beyond gravity

A closely related Stieltjes spectral framework appears in homogenization of advection–diffusion. For the diagonal effective diffusivity $$\mathrm{supp}\,\mathcal{K}^{-1}f \subset J^+(\mathrm{supp}\,f).$$4, the quantity is written as

$$\mathrm{supp}\,\mathcal{K}^{-1}f \subset J^+(\mathrm{supp}\,f).$$5

where $$\mathrm{supp}\,\mathcal{K}^{-1}f \subset J^+(\mathrm{supp}\,f).$$6 is a Stieltjes function,

$$\mathrm{supp}\,\mathcal{K}^{-1}f \subset J^+(\mathrm{supp}\,f).$$7

with $$\mathrm{supp}\,\mathcal{K}^{-1}f \subset J^+(\mathrm{supp}\,f).$$8 a positive spectral measure. The underlying operator is

$$\mathrm{supp}\,\mathcal{K}^{-1}f \subset J^+(\mathrm{supp}\,f).$$9

so $$(-1)^n \frac{d^n m}{d(\omega^2)^n} \geq 0 \quad \forall\, n\in\mathbb{N}_0,\;\omega^2>0$$0 is self-adjoint. For steady flows $$(-1)^n \frac{d^n m}{d(\omega^2)^n} \geq 0 \quad \forall\, n\in\mathbb{N}_0,\;\omega^2>0$$1 is a compact self-adjoint operator; for space-time periodic flows it is generally unbounded self-adjoint. Because $$(-1)^n \frac{d^n m}{d(\omega^2)^n} \geq 0 \quad \forall\, n\in\mathbb{N}_0,\;\omega^2>0$$2 is Stieltjes, Padé approximants provide nested rigorous bounds,

$$(-1)^n \frac{d^n m}{d(\omega^2)^n} \geq 0 \quad \forall\, n\in\mathbb{N}_0,\;\omega^2>0$$3

which induce upper and lower bounds on $$(-1)^n \frac{d^n m}{d(\omega^2)^n} \geq 0 \quad \forall\, n\in\mathbb{N}_0,\;\omega^2>0$$4. The paper’s main computational contribution is an iterative closed-form method for generating arbitrarily high spectral moments, thereby making the Stieltjes representation operational in both spatially and space-time periodic flows (Murphy et al., 2024).

The relaxation literature uses an analogous construction on the positive semiaxis. Spectral functions for the Havriliak–Negami, Jurlewicz–Weron–Stanislavski, and excess-wing models are treated as Stieltjes functions on $$(-1)^n \frac{d^n m}{d(\omega^2)^n} \geq 0 \quad \forall\, n\in\mathbb{N}_0,\;\omega^2>0$$5. Their Stieltjes character implies that the response function $$(-1)^n \frac{d^n m}{d(\omega^2)^n} \geq 0 \quad \forall\, n\in\mathbb{N}_0,\;\omega^2>0$$6 and the relaxation function $$(-1)^n \frac{d^n m}{d(\omega^2)^n} \geq 0 \quad \forall\, n\in\mathbb{N}_0,\;\omega^2>0$$7 are nonnegative. The memory kernel is determined by

$$(-1)^n \frac{d^n m}{d(\omega^2)^n} \geq 0 \quad \forall\, n\in\mathbb{N}_0,\;\omega^2>0$$8

and the Laplace exponent

$$(-1)^n \frac{d^n m}{d(\omega^2)^n} \geq 0 \quad \forall\, n\in\mathbb{N}_0,\;\omega^2>0$$9

is a complete Bernstein function. A partner memory is defined by

$\rho(\mu)$0

and the two memories satisfy the Sonine relation

$\rho(\mu)$1

Here the “parameter space” is the positive spectral axis itself, together with the Stieltjes class that constrains admissible frequency-domain response functions (Górska et al., 2021).

At a more abstract level, generalized Stieltjes transforms extend the same measure-to-function correspondence: $\rho(\mu)$2 A major structural result is that products remain generalized Stieltjes transforms: $\rho(\mu)$3 with an explicit representing measure $\rho(\mu)$4. The paper describes this as a Stieltjes convolution of measures and proves that the product measure has no singular continuous part. This gives a rigorous measure-theoretic calculus behind Stieltjes spectral parameterizations and their closure properties (Gomilko et al., 2022).

5. Random-matrix spectral variables and threshold domains

In random-matrix theory, Stieltjes spectral parameter spaces usually appear as complex or real domains on which resolvents are analytic. A generalized Wishart matrix,

$\rho(\mu)$5

is built from two $\rho(\mu)$6 matrices with

$\rho(\mu)$7

Its Stieltjes transform $\rho(\mu)$8 satisfies a cubic equation, which is the master equation for the limiting spectral density. The classical Marčenko–Pastur setting is recovered at $\rho(\mu)$9, while $M_0 = \int_0^\infty \rho(\mu)\,d\mu,$0 gives the “null” cross-correlation case; in the regime $M_0 = \int_0^\infty \rho(\mu)\,d\mu,$1, the density converges to the Wigner semicircle law. For $M_0 = \int_0^\infty \rho(\mu)\,d\mu,$2, the support undergoes a spectral phase transition at $M_0 = \int_0^\infty \rho(\mu)\,d\mu,$3: a gap around zero opens only when $M_0 = \int_0^\infty \rho(\mu)\,d\mu,$4, and a Dirac mass at $M_0 = \int_0^\infty \rho(\mu)\,d\mu,$5 appears for $M_0 = \int_0^\infty \rho(\mu)\,d\mu,$6 because the $M_0 = \int_0^\infty \rho(\mu)\,d\mu,$7 rank-two projectors cannot span the full $M_0 = \int_0^\infty \rho(\mu)\,d\mu,$8-dimensional space. In this setting, the spectral parameter space is the complex $M_0 = \int_0^\infty \rho(\mu)\,d\mu,$9-plane, with the physical branch selected by the condition $$\mu_{\rm eff}^{-n} = M_0^{-1}\int\rho(\mu)\,\mu^{-n}\,d\mu.$$0 as $$\mu_{\rm eff}^{-n} = M_0^{-1}\int\rho(\mu)\,\mu^{-n}\,d\mu.$$1 and by analyticity in the upper half-plane (Bouchaud et al., 2020).

A different random-matrix usage occurs in the spectral signal detection threshold problem. For noise covariance models with separable variance profile, the threshold $$\mu_{\rm eff}^{-n} = M_0^{-1}\int\rho(\mu)\,\mu^{-n}\,d\mu.$$2 is the right endpoint of the support of the limiting spectral distribution $$\mu_{\rm eff}^{-n} = M_0^{-1}\int\rho(\mu)\,\mu^{-n}\,d\mu.$$3,

$$\mu_{\rm eff}^{-n} = M_0^{-1}\int\rho(\mu)\,\mu^{-n}\,d\mu.$$4

The Stieltjes transform is evaluated on the real half-line

$$\mu_{\rm eff}^{-n} = M_0^{-1}\int\rho(\mu)\,\mu^{-n}\,d\mu.$$5

where it is real-valued and smooth, and this region is explicitly treated as the relevant spectral parameter domain for spike detection, denoising, and singular-value shrinkage. The paper further introduces an auxiliary implicit variable $$\mu_{\rm eff}^{-n} = M_0^{-1}\int\rho(\mu)\,\mu^{-n}\,d\mu.$$6 living in a constrained interval and develops nested Newton schemes that compute $$\mu_{\rm eff}^{-n} = M_0^{-1}\int\rho(\mu)\,\mu^{-n}\,d\mu.$$7, $$\mu_{\rm eff}^{-n} = M_0^{-1}\int\rho(\mu)\,\mu^{-n}\,d\mu.$$8, and $$\mu_{\rm eff}^{-n} = M_0^{-1}\int\rho(\mu)\,\mu^{-n}\,d\mu.$$9 to machine precision with asymptotic cost $\rho(\mu)$00 per evaluation scale (Leeb, 2019).

The local semicircle law provides the complementary complex-domain picture. For a normalized Wigner matrix, the empirical Stieltjes transform

$\rho(\mu)$01

is studied uniformly in the spectral parameter $\rho(\mu)$02. The imaginary part $\rho(\mu)$03 is the resolution scale: $\rho(\mu)$04 is a smoothed density estimate with kernel bandwidth $\rho(\mu)$05, and the local law controls

$\rho(\mu)$06

down to nearly optimal scale $\rho(\mu)$07. The paper states that, with high probability and uniformly in $\rho(\mu)$08,

$\rho(\mu)$09

for $\rho(\mu)$10. Here the spectral parameter space is the domain

$\rho(\mu)$11

and the dependence on $\rho(\mu)$12 determines both smoothing and stability (Götze et al., 2015).

6. Operator curves, holomorphic families, and scope of the concept

In operator theory, the term can refer to an explicit spectral curve. For the generalized Stieltjes operator $\rho(\mu)$13 on Sobolev-Lebesgue spaces, the spectrum is

$\rho(\mu)$14

The paper denotes the curve by

$\rho(\mu)$15

so that the spectral set is $\rho(\mu)$16. This curve is symmetric with respect to the real axis, passes through $\rho(\mu)$17 at $\rho(\mu)$18, lies inside the disk of radius $\rho(\mu)$19, and converges to $\rho(\mu)$20 as $\rho(\mu)$21. In special symmetric cases it collapses to a real interval, including $\rho(\mu)$22 for the classical Stieltjes operator on $\rho(\mu)$23 (Miana et al., 2019).

A second operator-theoretic usage concerns Stieltjes and inverse Stieltjes holomorphic families of linear relations on a separable Hilbert space. These are special Nevanlinna families defined on

$\rho(\mu)$24

with positivity or negativity on the negative real axis: $\rho(\mu)$25 for Stieltjes families, and

$\rho(\mu)$26

for inverse Stieltjes families. Their basic realization is by compressed resolvents of selfadjoint relations,

$\rho(\mu)$27

with uniqueness up to unitary equivalence under a minimality condition. The paper also shows that inversion and the map $\rho(\mu)$28 move between the Stieltjes and inverse Stieltjes classes. In this framework the spectral parameter space is the cut plane together with the order structure carried by the negative semi-axis (Arlinskiĭ et al., 2018).

The harmonic-analysis literature gives a further, more flexible perspective. In Fourier–Stieltjes algebras,

$\rho(\mu)$29

and for measures on any locally compact non-discrete Abelian group there are no non-trivial constraints between the norm, the spectral radius, and the supremum of the Fourier–Stieltjes transform, even when all convolution powers are singular. This shows that Stieltjes-related spectral data can organize a parameter space without forcing rigid algebraic relations among the most natural spectral quantities (Ohrysko et al., 2018).

The literature summarized here therefore suggests that “Stieltjes spectral parameter space” is not a single standardized object. In causal nonlocal gravity it is literally a two-parameter plane $\rho(\mu)$30; in transport theory it is the moment space of a positive spectral measure; in random-matrix theory it is typically a complex or real resolvent domain; and in operator theory it may be a Beta-function spectral curve or a holomorphic family on $\rho(\mu)$31. The common structure is the use of a positive measure or positivity-preserving analytic class to encode spectral information, together with a distinguished parameter—$\rho(\mu)$32, $\rho(\mu)$33, $\rho(\mu)$34, $\rho(\mu)$35, or $\rho(\mu)$36—that probes support, moments, thresholds, or resolvent branches.