---
title: Double Operator Integrals
url: https://www.emergentmind.com/topics/double-operator-integrals
type: topic
---

# Double Operator Integrals

Double operator integrals are operator-valued integrals of the form
\[
\iint \Phi(x,y)\,dE_A(x)\,T\,dE_B(y),
\]
where \(E_A\) and \(E_B\) are spectral measures, \(T\) is an operator, and \(\Phi\) is a scalar kernel. In the Birman–Solomyak framework they convert spectral data of two operators into a linear transformation on operators, and thereby provide a basic tool for perturbation theory, operator Lipschitz estimates, trace formulas, and functional calculus. In the classical self-adjoint setting the theory starts on Hilbert–Schmidt operators and extends, for suitable kernels, to bounded operators and Schatten ideals; later work develops multiple operator integrals, applies them to noncommuting pairs, and proposes generalized forms for non-Hermitian and continuous-spectrum settings [1509.02803][2510.25789][2503.16239].

## 1. Foundational construction

Let \(A\) and \(B\) be self-adjoint operators on a Hilbert space \(\mathscr H\), with spectral measures \(E_A\) and \(E_B\). The formal double operator integral
\[
\iint \Phi(x,y)\,dE_A(x)\,T\,dE_B(y)
\]
is first defined for Hilbert–Schmidt operators \(T\in \boldsymbol S_2\). The standard construction introduces a spectral measure \(\mathcal E\) on the product space by
\[
\mathcal E(\Lambda_1\times\Lambda_2)(T)=E_A(\Lambda_1)\,T\,E_B(\Lambda_2),
\]
and then sets
\[
\iint \Phi(x,y)\,dE_A(x)\,T\,dE_B(y)
=
\left(\int \Phi\,d\mathcal E\right)T.
\]
This gives a bounded map on \(\boldsymbol S_2\) with
\[
\left\|\iint \Phi\,dE_A\,T\,dE_B\right\|_{\boldsymbol S_2}
\le
\|\Phi\|_{L^\infty}\,\|T\|_{\boldsymbol S_2}
\]
[1509.02803][1505.07173].

A recent foundations paper reformulates the same starting point in terms of a product projection-valued measure \(G\) acting on the Hilbert space of Hilbert–Schmidt operators. For projection-valued measures \(E\) and \(F\) on Hilbert spaces \(H\) and \(K\), it constructs \(G\) on \(X\times Y\) so that
\[
G(\Gamma\times\Delta)=E(\Gamma)F(\Delta),
\]
and defines, for \(T\in (K,H)\),
\[
\iint \phi(x,y)\,dE(x)\,T\,dF(y)
:=
\int_{X\times Y}\phi\,dG\,T.
\]
For separated kernels \(\phi(x,y)=\alpha(x)\beta(y)\), this reduces to
\[
\iint \alpha(x)\beta(y)\,dE(x)\,T\,dF(y)
=
\left(\int \alpha\,dE\right)\,T\,\left(\int \beta\,dF\right).
\]
That paper also gives a new proof of the existence of the product of two projection-valued measures and a variant approach to arbitrary, not necessarily separable, Hilbert spaces [2510.25789].

In finite dimensions the same object becomes a finite double sum. If
\[
A=\sum_i \lambda_i P_i,\qquad B=\sum_j \mu_j Q_j,
\]
then
\[
T_\Phi^{A,B}(T)=\sum_{i,j}\Phi(\lambda_i,\mu_j)\,P_i\,T\,Q_j.
\]
This finite-dimensional formula is the discrete model behind much of the abstract theory [2503.16239].

## 2. Schur multipliers and tensor-product characterizations

The central boundedness question is: for which kernels \(\Phi\) does the transformation
\[
T\mapsto \iint \Phi(x,y)\,dE_A(x)\,T\,dE_B(y)
\]
extend from \(\boldsymbol S_2\) to \(\mathcal B(\mathscr H)\), or to a Schatten ideal \(\boldsymbol S_p\)? Such kernels are Schur multipliers. In the classical formulation, \(\Phi\) is a Schur multiplier with respect to \((E_A,E_B)\) precisely when it belongs to the corresponding multiplier class \(\mathfrak M(E_A,E_B)\) [1509.02803].

A basic representation uses the integral projective tensor product. If
\[
\Phi(x,y)=\int_\Omega \varphi(x,\omega)\psi(y,\omega)\,d\lambda(\omega)
\]
with
\[
\int_\Omega
\|\varphi(\cdot,\omega)\|_{L^\infty(E_A)}
\|\psi(\cdot,\omega)\|_{L^\infty(E_B)}
\,d\lambda(\omega)<\infty,
\]
then the DOI is bounded, and the kernel lies in
\[
L^\infty(E_A)\,\hat\otimes_{\mathrm i}\,L^\infty(E_B).
\]
Equivalently, Schur multipliers are described by the Haagerup tensor product
\[
L^\infty(E_A)\otimes_{\rm h}L^\infty(E_B),
\]
that is, by expansions
\[
\Phi(x,y)=\sum_{n\ge0}\varphi_n(x)\psi_n(y)
\]
with \(\{\varphi_n\}\) and \(\{\psi_n\}\) bounded in \(L^\infty(\ell^2)\). In the separable setting the Birman–Solomyak–Peller representation theorem identifies these descriptions of \(\mathfrak M(E_A,E_B)\) [1509.02803].

Wavelet methods sharpen this picture on \(\mathbb R^2\). A 2021 paper proves that if a bounded symbol \(k\) lies in
\[
B^{\frac{1}{p}-\frac{1}{2}}_{\frac{2p}{2-p},\,p}(\mathbb R,L_\infty(\mathbb R)),
\qquad 0<p<2,
\]
then \(k\) is an \(S_p\)-Schur multiplier, with
\[
\|k\|_{M_p}
\lesssim
\|k\|_{B^{\frac{1}{p}-\frac{1}{2}}_{\frac{2p}{2-p},\,p}(\mathbb R,L_\infty(\mathbb R))}
+
\|k\|_{L_\infty(\mathbb R^2)}.
\]
That argument extends Birman–Solomyak’s result to symbols without compact domain and recasts the multiplier theorem in vector-valued Besov and wavelet terms [2104.06567].

A later extension enlarges the admissible integrand class from measurable tensor-decomposable functions to a closure space obtained as limits of projective or integral projective tensor products. In that framework, any continuous function on a compact subset of \(\mathbb R^m\) can be uniformly approximated by finite sums of products of linear functions, and random DOIs and MOIs are defined by the same spectral-decomposition paradigm [2210.09392].

## 3. Perturbation theory for one operator

The most classical DOI formula expresses a nonlinear operator difference through a linear integral transform. If \(A\) and \(B\) are self-adjoint and
\[
f^{[1]}(s,t)=\frac{f(s)-f(t)}{s-t},
\]
then under standard hypotheses,
\[
f(A)-f(B)
=
\iint f^{[1]}(s,t)\,dE_A(s)\,(A-B)\,dE_B(t).
\]
Birman–Solomyak established this formula first when \(A-B\in \boldsymbol S_2\), and Peller’s survey presents the corresponding bounded-difference version when the divided difference is a Schur multiplier [1509.02803].

This DOI representation is the mechanism behind operator Lipschitz and operator differentiability results. A sufficient condition is
\[
f\in B_{\infty,1}^1(\mathbb R).
\]
In that case \(f\) is operator Lipschitz, and
\[
\|f(A)-f(B)\|
\le
C\,\|f\|_{B_{\infty,1}^1}\,\|A-B\|.
\]
The same Besov condition yields Fréchet differentiability:
\[
\frac{d}{dt}f(A+tK)\Big|_{t=0}
=
\iint \frac{f(s)-f(t)}{s-t}\,dE_A(s)\,K\,dE_A(t).
\]
For Hölder classes \(\Lambda_\alpha\), \(0<\alpha<1\), one obtains
\[
\|f(A)-f(B)\|
\le
C\,\|f\|_{\Lambda_\alpha}\,\|A-B\|^\alpha,
\]
and there are corresponding Schatten estimates such as
\[
A-B\in \boldsymbol S_p
\quad\Longrightarrow\quad
f(A)-f(B)\in \boldsymbol S_{p/\alpha}
\]
for \(f\in \Lambda_\alpha\) and \(1<p<\infty\) [1509.02803].

The same apparatus yields trace formulas. If \(A-B\in \boldsymbol S_1\), the DOI trace identity implies the Lifshits–Krein formula
\[
\operatorname{trace}(f(B)-f(A))
=
\int_{\mathbb R} f'(t)\,\xi(t)\,dt
\]
for appropriate \(f\), where \(\xi\) is the spectral shift function. The DOI proof proceeds by differentiating \(f(A+t(B-A))\), taking the trace of the derivative DOI, and integrating in \(t\) [1509.02803].

A recent foundations treatment also proves the Daletskii–Krein formula for strongly differentiable perturbations of a densely-defined self-adjoint operator. For \(f\in W_1(\mathbb R)\), with divided difference
\[
\phi_f(x,y)=
\begin{cases}
\dfrac{f(x)-f(y)}{x-y},&x\ne y,\\[4pt]
f'(x),&x=y,
\end{cases}
\]
it shows that
\[
\frac{d}{ds}f(A(s))
=
\iint \phi_f(x,y)\,dE_s(x)\,\Phi'(s)\,dE_s(y)
\]
for \(A(s)=A_0+\Phi(s)\) [2510.25789].

## 4. Noncommuting pairs and triple operator integrals

If \(A\) and \(B\) commute, the joint spectral theorem defines \(f(A,B)\) directly from a joint spectral measure. If they do not commute, there is no joint spectral measure, and the natural replacement is the DOI-based definition
\[
f(A,B)=\iint f(x,y)\,dE_A(x)\,dE_B(y),
\]
provided \(f\) is a Schur multiplier with respect to \((E_A,E_B)\). For bounded self-adjoint \(A,B\), the condition
\[
f\in B_{\infty,1}^1(\mathbb R^2)
\]
is sufficient to make this definition meaningful [1505.07173][1504.01189].

Perturbation theory for such \(f(A,B)\) no longer closes at the level of double integrals. Instead,
\[
f(A_1,B_1)-f(A_2,B_2)
\]
is represented as a sum of triple operator integrals involving the partial divided differences
\[
\mathfrak D_1 f(x_1,x_2,y)=\frac{f(x_1,y)-f(x_2,y)}{x_1-x_2},
\qquad
\mathfrak D_2 f(x,y_1,y_2)=\frac{f(x,y_1)-f(x,y_2)}{y_1-y_2}.
\]
The natural tensor spaces for these kernels are not the ordinary Haagerup tensor products but the asymmetric Haagerup-like tensor products
\[
L^\infty\otimes_{\rm h}L^\infty\otimes^{\rm h}L^\infty,
\qquad
L^\infty\otimes^{\rm h}L^\infty\otimes_{\rm h}L^\infty.
\]
For \(f\in B_{\infty,1}^1(\mathbb R^2)\), \(\mathfrak D_1 f\) and \(\mathfrak D_2 f\) belong to these spaces with norms controlled by \(\|f\|_{B_{\infty,1}^1}\) [1505.07173].

This leads to the sharp Schatten estimate
\[
\|f(A_1,B_1)-f(A_2,B_2)\|_{\boldsymbol S_p}
\le
C\,\|f\|_{B_{\infty,1}^1}
\big(\|A_1-A_2\|_{\boldsymbol S_p}+\|B_1-B_2\|_{\boldsymbol S_p}\big),
\qquad 1\le p\le 2.
\]
The same line of work proves that no analogous Lipschitz estimate holds in \(\boldsymbol S_p\) for \(p>2\), and no such global Lipschitz estimate holds in operator norm, even for bounded functions with compactly supported Fourier transform. The restriction \(p\le 2\) is therefore sharp [1504.01189][1505.07173].

The same framework extends to noncommuting unitary pairs, and triple operator integrals also yield trace-class commutator estimates for almost commuting self-adjoint operators, which in turn lead to an extension of the Helton–Howe trace formula to arbitrary functions in \(B_{\infty,1}^1(\mathbb R^2)\) [1509.02803].

## 5. Multiple operator integrals and higher-order analysis

Double operator integrals are the \(m=2\) case of multiple operator integrals. For spectral measures \(E_1,\dots,E_m\) and an integrand \(\Psi\) in the integral projective tensor product
\[
L^\infty(E_1)\hat\otimes_{\mathrm i}\cdots\hat\otimes_{\mathrm i}L^\infty(E_m),
\]
the \(m\)-fold integral
\[
\underbrace{\int\cdots\int}_{m}
\Psi(x_1,\dots,x_m)\,dE_1(x_1)T_1\,dE_2(x_2)\cdots T_{m-1}\,dE_m(x_m)
\]
is well defined and bounded. Higher divided differences \(\mathfrak D^m f\) belong to these tensor products when
\[
f\in B_{\infty,1}^m(\mathbb R),
\]
and this yields formulas for higher derivatives:
\[
\frac{d^m}{dt^m}f(A+tK)\Big|_{t=0}
=
m!\underbrace{\int\cdots\int}_{m+1}
\mathfrak D^m f(s_1,\dots,s_{m+1})\,dE_A(s_1)K\cdots K\,dE_A(s_{m+1}).
\]
The same machinery controls higher-order differences and higher-order trace formulas under perturbations of class \(\boldsymbol S_m\) [1509.02803].

A probabilistic extension constructs random MOIs and random DOIs from random operators defined through spectral decompositions. In that setting, one obtains tail bounds for norms of higher random operator derivatives, higher random operator differences, and Taylor remainders of random operator-valued functions. The DOI formula
\[
\frac{d}{dt}\psi(X(t))
=
T_{X(t),X(t)}^{\psi^{[1]}}\!\left(\frac{dX(t)}{dt}\right)
\]
appears as the first-order instance of that random MOI calculus [2210.09392].

These developments suggest that DOI theory is not merely a two-variable perturbation device. A plausible implication is that it is the basic rank-two component of a hierarchy of operator integral formulas governing derivatives, Taylor expansions, and trace identities.

## 6. Generalized and nonclassical variants

Several recent papers propose generalized double operator integrals beyond the classical self-adjoint or normal setting. In finite dimensions, Generalized Double Operator Integrals are defined for arbitrary complex matrices by replacing the spectral-measure framework with a projector–nilpotent decomposition. If
\[
X=\sum_{k,i}\lambda_k P_{k,i}+\sum_{k,i}N_{k,i},
\]
then the generalized DOI
\[
T_\Phi^{X_1,X_2}(Y)
\]
is a finite sum involving the values and partial derivatives of \(\Phi\) at eigenvalue pairs and the operators
\[
P_{k_1,i_1}YP_{k_2,i_2},\quad
P_{k_1,i_1}YN_{k_2,i_2}^{q_2},\quad
N_{k_1,i_1}^{q_1}YP_{k_2,i_2},\quad
N_{k_1,i_1}^{q_1}YN_{k_2,i_2}^{q_2}.
\]
When \(X_1\) and \(X_2\) are Hermitian and all nilpotent parts vanish, this reduces to the classical finite-dimensional DOI. In that setting the map \(\Phi\mapsto T_\Phi^{X_1,X_2}\) is a linear homomorphism, one has perturbation formulas such as
\[
f(X_1)-f(X_2)=T_{f^{[1]}}^{X_1,X_2}(X_1-X_2),
\]
and one derives Lipschitz and Hölder estimates in Frobenius norm, with constants depending on derivatives of the divided difference and on nilpotent orders [2503.16239][2506.19971].

A related 2025 program extends generalized multiple operator integrals to operators with continuous spectra by replacing classical projection-valued spectral measures with a projector-plus-nilpotent spectral decomposition adapted to non-normal operators. In that framework generalized DOIs appear as the \(c=1\) case of generalized MOIs, inherit decomposition, continuity, norm, and Lipschitz estimates, and are used to obtain a Kreĭn-type spectral shift formula and its higher-order analogues in the continuous-spectrum setting [2507.23049].

A conceptually different construction in quantum stochastic calculus studies causal double product integrals over the triangular domain \(a<x<y<b\),
\[
\prod_{a<x<y<b}\Big(1 + i\frac{\lambda}{2}(dP_x dQ_y - dQ_x dP_y) + i \frac{\mu}{2}(dP_x dP_y + dQ_x dQ_y)\Big).
\]
That paper explicitly presents this object as playing a role similar to a DOI, while emphasizing the differences: the kernel is operator-valued, the domain is causal, and the construction proceeds by product integrals and second quantization rather than spectral measures. This suggests a broader family of two-variable operator constructions in which DOI ideas survive outside classical spectral theory [1506.04294].

Double operator integrals therefore occupy a junction between spectral integration, perturbation theory, tensor-product operator spaces, and noncommutative analysis. In the classical theory they linearize \(f(A)-f(B)\), organize trace and differentiability formulas, and control functions of noncommuting pairs through higher operator integrals. In more recent generalizations they also serve as a template for extending perturbation theory to non-Hermitian matrices, generalized continuous-spectrum operators, and stochastic noncommutative settings.

Source: https://www.emergentmind.com/topics/double-operator-integrals