Cwikel's bound reloaded

Abstract: There are a couple of proofs by now for the famous Cwikel--Lieb--Rozenblum (CLR) bound, which is a semiclassical bound on the number of bound states for a Schrödinger operator, proven in the 1970s. Of the rather distinct proofs by Cwikel, Lieb, and Rozenblum, the one by Lieb gives the best constant, the one by Rozenblum does not seem to yield any reasonable estimate for the constants, and Cwikel's proof is said to give a constant which is at least about 2 orders of magnitude off the truth. This situation did not change much during the last 40+ years.It turns out that this common belief, i.e, Cwikel's approach yields bad constants, is not set in stone: We give a drastic simplification of Cwikel's original approach which leads to an astonishingly good bound for the constant in the CLR inequality. Our proof is also quite flexible and leads to rather precise bounds for a large class of Schrödinger-type operators with generalized kinetic energies. Moreover, it highlights a natural but overlooked connection of the CLR bound with bounds for maximal Fourier multipliers from harmonic analysis.

\section{Introduction}\label{introduction}
We want to find natural bounds, with the right semi-classical behavior, for the number of negative eigenvalues of Schr\"odinger  operators $P^2+V$ with $P=-i\nabla$, the momentum operator, or more general operators like the polyharmonic Schr\"odinger operators $|P|^{2\alpha}+V$, including the ultra-relativistic operator $|P|+V$. We will also consider operator-valued potentials $V$.

For the one-particle Schr\"odinger operator $P^2+V$ with $P=-i\nabla$ the momentum operator and $V$ a real-valued potential, this type of bound goes back to Cwikel, Lieb, and Rozenblum \cite{Cwikel,Lieb,Ro1,Ro2}, with very different proofs. They prove
\begin{align}\label{eq:original-CLR}
	N(P^2+V) \le L_{0,d}  \int_{{\mathbb{R}}^d} V_-(x)^{d/2}\, dx
\end{align}
for the number of negative eigenvalues of a Schr\"odinger operator.
This bound is a semi-classical bound since a simple scaling argument shows that the classical phase-space volume of the region of negative energy is given by
\begin{align}
	N^\text{cl}(\eta^2+V)= \iint_{\eta^2+V(x)<0} 1\,\frac{d\eta\, dx}{(2\pi)^d} =  \frac{|B_1^d|}{(2\pi)^d}\int_{{\mathbb{R}}^d} V_-(x)^{d/2}\, dx \, .
\end{align}
where $|B_1^d|$ is the volume of the unit ball in ${\mathbb{R}}^d$.

The intuition is that the uncertainty principle forces a quantum particle to occupy roughly a classical phase-space
volume $(2\pi)^d$. Thus  $N^\text{cl}(\eta^2+V)$, counting the
volume where  the classical Hamiltonian energy
$H(\eta,x) = \eta^2+V(x)$ is negative, should control $N(P^2+V)$.
The CLR bound  \eqref{eq:original-CLR} shows that this is the case modulo the factor\footnote{We write $L_{0,d}$ etc, since
there are a class of inequalities for the $\gamma^\text{th}$ moment
of the negative eigenvalues with associated constants
$L_{\gamma,d}$, see \cite{LT1,LT2} and the reviews \cite{LaWe-review,Hu-review}.} $C_{0,d}=L_{0,d}\frac{(2\pi)^d}{|B_1^d|}$ .

The original bounds on $C_{0,d}$ in \cite{Cwikel} and \cite{Lieb}
were explicitly dimension dependent with a considerable growth in the dimension $d$. The bound due to Lieb grows like $C_{0,d}= \sqrt{\pi d}(1+O(d^{-1}))$.
See \cite{Simon} or \cite[Chapter 3.4]{Roepstorff} for an excellent discussion of Lieb's method and Remark \ref{rem:explicit-numbers} below for some explicit numbers.
However, it is expected that semi-classical arguments work better in high  dimensions. In particular, the constant $C_{0,d}$ should not grow in $d$.
The first dimension independent bound  $C_{0,d}\le 81$ was derived
by extending Cwikel's method to operator-valued potentials in 2002 in \cite{Hu}. This work extended an induction in the dimension argument by Laptev and Weidl\footnote{See also \cite{La} for some indication of the induction in dimension trick.} \cite{LaWe}, who were the first to derive Lieb--Thirring bounds with the sharp classical Lieb--Thirring constant in all dimensions in some cases.
Although the upper bound from \cite{Hu} is dimension independent, it is certainly too large for small dimensions.

For the last $40$-plus years it has been believed that any approach
based on Cwikel's method cannot, in low dimensions, yield any bounds on $C_{0,d}$ which are comparable to the ones obtained by Lieb.
This is wrong, as we will show by drastically simplifying and, at the same time, generalizing the important ideas of Cwikel.
A typical result which can be easily achieved with our method is
\begin{theorem}\label{thm:CLR}
	The number $N(P^2+V)$ of negative energy bound states of  $P^2+V$ obeys the semiclassical bound
	\begin{align}
		N(P^2+V)\le C_{0,d} \frac{|B^d|}{(2\pi)^d} \int_{{\mathbb{R}}^d} V_-(x)^{d/2}\, \d x
	\end{align}
	for all $d\ge 3$, where $B^d$ is the unit ball in ${\mathbb{R}}^d$, $|B^d|$ its volume, and the constant $C_{0,d}$ given in Table \ref{tb:better} below.

	Moreover, the same bounds with the same constants also hold in the operator-valued case, see Theorem \ref{thm:great-operator-valued}.
\end{theorem}
\begin{remarks}\label{rem:explicit-numbers}
\refstepcounter{theoremi} {\rm(\roman{theoremi})}{~} Table \ref{tb:better} below compares the upper bounds on $C_{0,d}$,  obtained with our method, with the best known ones so far for scalar and operator-valued potentials,
  \begin{table}[ht!]
  \begin{tabular}{| l | c | c | c|}
    \hline
    \multirow{2}{*}{$d$} & Our results  & \multicolumn{2}{c|}{Best known so far} \\
        & scalar and operator valued  & scalar & operator-valued         \\\hline
    3   & 7.55151 & 6.86924  &  \multirow{7}{*}{10.332}\\
    4   & 6.32791 & 6.03398  &  \\
    5   & {5.95405} & 5.96677 &   \\
    6   & {5.77058} & 6.07489 &  \\
    7   & {5.67647} & 6.24464 &  \\
    8   & {5.63198} & 6.43921 &  \\
    9   & {5.62080}  & 6.64378 &  \\ \hline

  \end{tabular}
  \vskip1ex
  \caption{}
  \label{tb:better}
  \end{table}
  All bounds on $C_{0,d}$ in the third column of the table were obtained already in the original work of Lieb more than 40 years ago\footnote{The numbers are taken from Roepstorff's book \cite[Table 3.1]{Roepstorff}}.
 Our bounds on $C_{0,d}$ also hold in the operator-valued case, see Section \ref{sec:operator-valued} below. The last column is due to Frank, Lieb and
 Seiringer \cite{FLS} and holds for all $d\ge 3$.
 Our result also gives the bound $C_{0,d}\le 5.62080 $ for $d\ge 9$, see the discussion in Appendix \ref{sec:induction-in-dimension}.  For dimensions $d=3,\ldots,9$ our upper bounds are compared with the values of the lower bound \eqref{eq:CLR-lower-bound} achievable by our method in Table \ref{tb:lower-bound} below. \\[0.2em]
\refstepcounter{theoremi} {\rm(\roman{theoremi})}{~} There have been several previous attempts to improve on Lieb's result, for example, due to Conlon \cite{Co}, Li and Yau \cite{LY}, Frank \cite{Fr}, and Weidl \cite{Weidl1,Weidl2}. All these very much different proofs shed a new light on the  Cwikel--Lieb--Rozenblum bound, but failed to give better bounds on the involved constants than already achieved by Lieb. \\[0.2em]
\refstepcounter{theoremi} {\rm(\roman{theoremi})}{~} From the point of view of physics, the other important case is $\alpha=1/2$, which corresponds to an ultra-relativistic Schr\"odinger operator $|P|+V$. In three dimensions we get the upper bound
  \begin{align}
  	N(|P|+V) \le 5.77058 \int_{{\mathbb{R}}^3} V_-(x)^{3}\, dx
  \end{align}
  which improves the result of Daubechies \cite{Daubechies}, who gets $N(|P|+V) \le 6.08 \int_{{\mathbb{R}}^3} V_-(x)^{3}\, dx$.
 \end{remarks}

For more general so-called polyharmonic Schr\"odinger-type operators our method yields the following for scalar potentials.
A similar result, with the same constants, also holds for operator-valued potentials, see Theorem \ref{thm:poly-operator-valued}.

\begin{theorem}\label{thm:poly}
	Let $P=-i\nabla$ be the momentum operator, $V=V_+-V_-$ be a real-valued potential with positive part $V_+\in L^1_{\text{\rm{loc}}}$ and negative part
	$V_-\in L^{d/\alpha}({\mathbb{R}}^d)$ with $0<\alpha<d/2$, and $P^{2\alpha}+V$ the Schr\"odinger--type operator defined via quadratic form methods on $L^2({\mathbb{R}}^d)$. Then the number $N(P^{2\alpha}+V)$ of negative energy bound states of  $P^{2\alpha}+V$ is bounded by
	\begin{align}\label{eq:poly}
		N(P^{2\alpha}+V)
			\le C_{d/\alpha}\,
			\frac{|B_1^d|}{(2\pi)^d} \int_{{\mathbb{R}}^d} V_-(x)^{\frac{d}{2\alpha}}\, dx
	\end{align}
	with constant
	\begin{align}\label{eq:C-gamma}
		C_{\gamma} = \frac{\gamma^{\gamma +1}}{4\left( \gamma-2 \right)^{\gamma-2}}
						   M_\gamma \ ,
	\end{align}
  where
  \begin{align}\label{eq:M-gamma}
  	M_\gamma= \inf\left\{
  						\left( \|m_1\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})}
  								\|m_1\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})}
  						\right)^{\gamma-2}
  							\int_0^\infty (1-t^{-1}m(t))^2 \, t^{1-\gamma}\, dt \,
  				  \right\},
  \end{align}
  the infimum taken over all $m_1,m_2\in L^2({\mathbb{R}}_+,\frac{ds}{s})$, and $m= m_1*m_2$ denotes the convolution of $m_1, m_2$ on ${\mathbb{R}}^+$ with measure $ds/s$.
\end{theorem}
The minimization problem for $M_\gamma$ in \eqref{eq:M-gamma} seems to be new. As Theorem \ref{thm:poly}  shows, it has considerable implications for the spectral theory of Schr\"odinger operators.
For the constant $M_\gamma$ above, we note the following estimate. \begin{proposition}\label{prop:simpel-bound-M-gamma}
	For all $\gamma>2$
	\begin{align}\label{eq:easy-estimates}
		\frac{2}{\gamma(\gamma-1)(\gamma-2)}
		\le M_\gamma
		\le \frac{8 }{\gamma(\gamma-2)(\gamma+ 2) }
	\end{align}
\end{proposition}
For the proof see Sections \ref{sec:lower-bound} and Remark \ref{rem:simple-choice} below.

\begin{remarks}
\refstepcounter{theoremi} {\rm(\roman{theoremi})}{~}
	Even this simple upper bound on $M_\gamma$ yields result better than available results so far in the literature: Using ideas from Rumin \cite{Ru1, Ru2}, Frank \cite{Fr} got the bound
	\begin{align*}
		N(P^{2\alpha}+V)
			\le \left(
					\frac{d(d +2\alpha)}{(d-2\alpha)^2}
				\right)^{(d-2\alpha)/(2\alpha)} \frac{d}{d -2\alpha}
			\frac{|B_1^d|}{(2\pi)^d} \int_{{\mathbb{R}}^d} V_-(x)^{\frac{d}{2\alpha}}\, dx\, .
	\end{align*}
   Computing the ratio of the constants in Frank's bound and the one from \eqref{eq:poly}, using the upper bound in \eqref{eq:easy-estimates}, one sees that our bound from Theorem
   \ref{thm:poly}  is better in the whole allowed range of $0<\alpha<d/2$. \\[0.2em]
\refstepcounter{theoremi} {\rm(\roman{theoremi})}{~} For the constant $C_\gamma$ in \eqref{eq:C-gamma},  the lower bound from \eqref{eq:easy-estimates} yields
   \begin{align*}
 	  C_\gamma \ge \frac{\gamma^{\gamma}}{2(\gamma-1)\left( \gamma-2 \right)^{\gamma-1}} \eqqcolon C^\text{lower}_\gamma\, ,
   \end{align*}
  where $C^\text{lower}_\gamma$ is a, probably non-sharp, lower bound for the best possible constant achievable by our method\footnote{Which is not necessarily the best possible constant.}.
   Thus the upper bound on $C_\gamma$ shows
  \begin{align*}
  	\frac{C_\gamma}{C^\text{lower}_\gamma}\le  4\frac{\gamma-1}{\gamma+2}<4\, ,
  \end{align*}
  where $\gamma=d/\alpha> 2$.
  So our upper bound is less than a factor of 4 off the lower bound.  \\
\refstepcounter{theoremi} {\rm(\roman{theoremi})}{~} The above lower bound also gives the lower bound
\begin{align}\label{eq:CLR-lower-bound}
	C^\text{lower}_{0,d} = C^\text{lower}_d =  \frac{d^{d}}{2(d-1)\left( d-2 \right)^{d-1}}
\end{align}
for the constant in Theorem \ref{thm:CLR}, achievable by our method. For dimensions $3\le d \le 9$ this gives the lower bounds
  \begin{table}[h!]
  \begin{tabular}{| l | c | c |}
    \hline $d$ & Our results  & lower bound    \\\hline
    3   & 7.55151 & 6.75000  \\
    4   & 6.32791 & 5.33333  \\
    5   & {5.95405} & 4.82253  \\
    6   & {5.77058} & 4.55625   \\
    7   & {5.67647} & 4.39229   \\
    8   & {5.63198} & 4.28088   \\
    9   & {5.62080}  & 4.20028   \\ \hline
  \end{tabular}
  \vskip1ex
  \caption{}
  \label{tb:lower-bound}
  \end{table}

 In addition,
 \begin{align*}
 	C^\text{lower}_{0,d} = \frac{d^2}{2(d-1)(d-2)}\left( 1+\frac{2}{d-2} \right)^{d-2} \rightarrow \frac{e^2}{2}\ge 3.69452\, .
 \end{align*}
 This comparison shows that there is not too much room to improve
 on the upper bounds we obtained, even if one finds the sharp value in the minimization problem for $M_\gamma$ in \eqref{eq:M-gamma}.

\refstepcounter{theoremi} {\rm(\roman{theoremi})}{~} It is known that if $\alpha\ge d/2$, the operator $P^{2\alpha}-U$ always has bound states for nontrivial $U\ge 0$, so a quantitative bound of the form $N(P^{2\alpha}-U)\lesssim \int_{{\mathbb{R}}^d}U(x)^{d/\alpha}$ cannot hold if $\alpha\ge d$.
For $\alpha=1$ see \cite{Simon2} or \cite[Problem 2 in \S 45]{LL}.
For more general cases, see \cite{NW96,LSW,Pankrashkin}, and \cite{HHRW} for a simple proof of how the existence/ non-existence of a CLR type bound for operators of the form $T(P)+V$ for a large class of functions $T:{\mathbb{R}}^d\to [0,\infty)$ is related to the behavior of the symbol $T$ close to its zero-set.
\end{remarks}

As we mentioned before, our method can be generalized to  operator valued potentials, leading to the same results.
To formulate this, we need some additional notation.
An operator-valued potential $V$ is a map $V:{\mathbb{R}}^d\mapsto \mathcal{G}$ with $V(x):\mathcal{G}\to \mathcal{G}$ a bounded self-adjoint operator on an auxiliary Hilbert space\footnote{In the following, all Hilbert-spaces are considered to be separable, unless said otherwise ;-).  Physically, this auxiliary Hilbert-space corresponds to other of degrees of freedom, for example spin.} $\mathcal{G}$ for almost all $x\in{\mathbb{R}}^d$. We denote by $\mathcal{B}(\mathcal{G})$ the bounded operators on $\mathcal{G}$ and by ${\mathcal{S}}_\alpha(\mathcal{G})$ the von Neumann--Schatten ideal of compact operators on $\mathcal{G}$ with $\alpha$ summable singular values, see for example \cite{Simon-trace-ideals} for a background on  von Neumann--Schatten ideals.

\begin{theorem}[Operator-valued version of Theorem \ref{thm:poly}]\label{thm:poly-operator-valued}
 Let $\mathcal{G}$ be a Hilbert space  and $V:{\mathbb{R}}^d\to \mathcal{B}(\mathcal{G})$ an operator valued potential with positive part $V_+\in L^1_{\text{loc}}({\mathbb{R}}^d,\mathcal{B}(\mathcal{G}))$ and negative part  $V_-\in L^{d/(2\alpha)}({\mathbb{R}}^d, {\mathcal{S}}_{d/(2\alpha)}(\mathcal{G}))$. Then the number of negative energy bound states of $P^{2\alpha}\otimes\idG+V$ is bounded by
 \begin{align}
 	N(P^{2\alpha}\otimes\idG+V)
			\le C_{d/\alpha}\,
			\frac{|B_1^d|}{(2\pi)^d} \int_{{\mathbb{R}}^d} \mathop{\mathrm{tr}} \nolimits_\mathcal{G}[ V_-(x)^{\frac{d}{2\alpha}}]\, dx
 \end{align}
 with the same constant $C_\gamma$ as in Theorem \ref{thm:poly}.
\end{theorem}

For the physically most interesting case $\alpha=1$ this enables us to get  considerable  improvements on the constants in the Cwikel--Lieb--Rozenblum bound.

\begin{theorem}[Operator-valued version of Theorem \ref{thm:CLR}]\label{thm:CLR-op-valued}
	Let $\mathcal{G}$ be a Hilbert space  and $V:{\mathbb{R}}^d\to \mathcal{B}(\mathcal{G})$ an operator valued potential with positive part $V_+\in L^1_{\text{loc}}({\mathbb{R}}^d,\mathcal{B}(\mathcal{G}))$ and negative part  $V_-\in L^{d/2}({\mathbb{R}}^d, {\mathcal{S}}_{d/2}(\mathcal{G}))$. Then the number of negative energy bound states of $P^{2}\otimes\idG+V$ is bounded by
 \begin{align}
 	N(P^{2}\otimes\idG+V)
			\le C_{0,d}^{\mathrm{op}}\,
			\frac{|B_1^d|}{(2\pi)^d} \int_{{\mathbb{R}}^d} \mathop{\mathrm{tr}} \nolimits_\mathcal{G}[ V_-(x)^{\frac{d}{2}}]\, dx
 \end{align}
 with
 \begin{align}\label{eq:constants-op-valued}
 	C_{0,d}^{\mathrm{op}} = \min_{3\leq n \leq d} C_{0,n}^{\mathrm{op}} \leq \min_{3\leq n \leq d} C_n,
 \end{align}
 where $C_n$ is given by \eqref{eq:C-gamma} for $\gamma=n$.
\end{theorem}

\begin{remark}
	Table \ref{tb:better} lists	upper bounds on $C_{0,d}^{\mathrm{op}}$ for dimensions $d=3\ldots 9$, see also Appendix \ref{sec:numerics}. The constant for $d=9$ is also an upper bound on $C_{0,d}^{\mathrm{op}}$ in any dimension $d\geq 10$ by \eqref{eq:constants-op-valued}.
\end{remark}

The structure of the paper is as follows. In Section \ref{sec:CLR} we present the main ideas of our method in the case of a standard non-relativistic Schr\"odinger operator. The extension to more general kinetic energies is done in Section \ref{sec:gen-kin-energy}.

In Section \ref{sec:max-F-mult} we explain the surprising connection of semiclassical bounds and maximal Fourier
multiplier estimates, which is probably the most important new part
of our method.

Although we cannot explicitly find minimizers of the variational problem from Theorem \ref{thm:poly}, there is a natural lower bound, which is discussed in Section \ref{sec:lower-bound}. The numerical study to find reasonable  upper bounds for this variational problem is presented in Appendix \ref{sec:numerics}.

The extension to the operator-valued setting is done in Sections \ref{sec:operator-valued} and \ref{sec:trace-ideals}.
In particular, in Section \ref{sec:trace-ideals} we prove a
fully operator-valued version of Cwikel's original weak trace ideal bound.

\section{The splitting trick}\label{sec:CLR}
The main idea in the proof of Theorem \ref{thm:CLR} is quite simple.
Let $U\coloneqq V_-\ge 0$. As quadratic forms $P^2+V\ge P^2-U$.
This and the  Birman--Schwinger principle show that
\begin{align*}
	N(P^2+V)\le N(P^2-U) = n(U^{1/2}|P|^{-2}U^{1/2};1)
\end{align*}
where $N(P^2+V)$ are the number of negative eigenvalues of $P^2-U$ and $n(A;\kappa)$ the number of singular values $(s_j(A))_{j\in{\mathbb{N}} }$ greater than $\kappa>0$ of a compact operator $A$.

We denote by $\mathcal{F}$ the Fourier transform and by  $\mathcal{F}^{-1}$ its inverse, by $M_h$ the operator of  multiplication with a function $h$,
and  $A=A_{f,g} = M_f\mathcal{F}^{-1}M_g$ for $f,g$ non-negative (measurable) functions on ${\mathbb{R}}^d$. When $f(x)=U(x)^{1/2}$ and $g(\eta)= |\eta|^{-1/2}$ then $A A^* = U^{1/2}|P|^{-2}U^{1/2}$, which has the same non-zero eigenvalues as $A^*A$. Thus
\begin{align*}
	N(P^2-U) = n(A_{f,g};1)
\end{align*}
In particular, the Chebyshev--Markov inequality gives
\begin{align*}
	N(P^2-U)
	  & = n(A_{f,g};1)
	  	\le \sum_{j} \frac{(s_j(A_{f,g})-\mu)_+^2}{(1-\mu)^2}
\end{align*}
for any $0<\mu<1$.
Let's drop the dependence of $A$ on $f$ and $g$ for the moment.
We want to split $A= B+ H$, where $B$ is bounded,  and note that Ky Fan's inequality for the singular values  yields
\begin{align*}
	s_j(A) &= s_j(B+ H)
	  \le s_1(B) + s_j(H) \\
	&\le \|B\| + s_j(H)
\end{align*}
for all $j\in{\mathbb{N}} $. So if $\|B_{f,g}\|\le \mu<1$ we get
\begin{align}\label{eq:first-upper-bound}
	N(P^2-U) \le (1-\mu)^{-2} \sum_{j\in{\mathbb{N}}} s_j(H_{f,g})^2
	  = (1-\mu)^{-2} \|H_{f,g}\|_{HS}^2,
\end{align}
where $\|H\|_{HS}$ denotes the Hilbert--Schmidt norm of an operator $H$.

In order to make the above argument work, one has to be able to split $A_{f,g}= B_{f,g}+H_{f,g}$ in such a way that the Hilbert-Schmidt norm of $H_{f,g}$ is easy to calculate and one has a good bound on the operator norm of $B_{f,g}$.
It will turn out, see the calculation below, that $\|H_{f,g}\|_{HS}^2= c \int_{R^d}f(x)^d\, dx$, so the right hand side of \eqref{eq:first-upper-bound} has exactly the right (semi-classical) scaling in $f$.
But, in order to use this in \eqref{eq:first-upper-bound}, it also enforces that the upper bound $\mu$ on the operator norm of $B_{f,g}$ has to be \emph{independent} of $f$. Since for a given $\varphi\in L^2$ one can freely choose $f\ge 0$ as to make $|B_{f,g,m}\varphi|$ as big as possible, this leads naturally to the associated maximal operator $B^*_{g,m}\coloneqq\sup_{f\ge 0}|B_{f,g,m}\varphi|$.
Although this is not explicitly written in the paper by Cwikel, getting a useful bound on such a type of maximal operator is exactly what he achieved in \cite{Cwikel}, using  a dyadic decomposition in the ranges of $f$ and $g$ and collecting suitable terms.
We will do this in a much simpler and more effective way. This enables us to get a constant which is more than 10 times smaller than the original constant by Cwikel.

Writing out the Fourier transform, one sees that  $A_{f,g}$ has kernel
\begin{align}
	A_{f,g}(x,\eta)= (2\pi)^{-d/2} e^{ix\cdot\eta} f(x)g(\eta),
\end{align}
that is,
\begin{align}
	A_{f,g}\varphi(x)= f(x)\mathcal{F}^{-1}(g\varphi)(x)
	= (2\pi)^{-d/2}  \int_{{\mathbb{R}}^d} e ^{ix\cdot\eta} f(x)g(\eta)\varphi(\eta)\, \d\eta ,
\end{align}
at least for nice enough $\varphi$.

In order to write $A_{f,g}$ as a sum of a bounded and a Hilbert-Schmidt operator, set $t=f(x) g(\eta)$, split $t= m(t) +t-m(t)$ for some function $m:[0,\infty)\to {\mathbb{R}} $, and define $B=B_{f,g,m}$ and $H_{f,g,m}$ via their kernels
\begin{align}\label{eq:kernel-B}
	B_{f,g,m}(x,\eta) &= (2\pi)^{-d/2} e^{ix\cdot\eta}m(f(x)g(\eta)) \\
	H_{f,g,m}(x,\eta) &=  (2\pi)^{-d/2} e^{ix\cdot\eta} \left(f(x)g(\eta)-m(f(x)g(\eta))\right) \label{eq:kernel-H}
\end{align}
It is then clear that $A_{f,g}= B_{f,g,m}+H_{f,g,m}$. Our starting point is that the Hilbert--Schmidt norm of $H_{f,g,m}$ is easy to calculate and it is not too hard to get an explicit bound on the operator norm of $B_{f,g,m}$ on $L^2$ under a suitable assumption on $m$.
\begin{theorem}\label{thm:great}
	The Hilbert--Schmidt norm of $H_{f,g,m}$ is given by
	\begin{align}
		\|H_{f,g,m}\|_{HS}^2= \int_{{\mathbb{R}}^d} G_{g,m}(f(x)) \, dx
	\end{align}
	where $G_{g,m}$ is given by
	\begin{align}
		G_{g,m}(u)= \int_{{\mathbb{R}}^d} |ug(\eta)-m(ug(\eta))|^2 \frac{\d \eta}{(2\pi)^d} .
	\end{align}
    If, moreover, $m$ is given by a convolution,  that is,
    \begin{align}
    	m(t) = m_1*m_2(t)=\int_0^\infty m_1(t/s)m_2(s)\frac{\d s}{s}
    \end{align}
    then for all measurable non-negative functions $f$ and $g$ the operator $B_{f,g,m}$ is bounded on $L^2({\mathbb{R}}^d)$ with
    \begin{align}\label{eq:B-bound}
    	\|B_{f,g,m}\varphi\|_{2} \le \left( \int_0^\infty |m_1(s)|^{2} \frac{d s}{s}\right)^{1/2}\left( \int_0^\infty |m_2(s)|^{2} \frac{d s}{s}\right)^{1/2}\|\varphi\|_2
    \end{align}
    for all $\varphi\in L^2({\mathbb{R}}^d)$.
\end{theorem}
\begin{remark}
	 We stress the fact, that the bound on the operator norm of $B_{f,g,m}$ is \emph{independent} of the choice of $f$, as it has to be, and also of $g$. This will turn out to be a natural consequence of
	 the convolution structure of $m$,
	 see Section \ref{sec:max-F-mult} below.
\end{remark}
\begin{proof}[Proof of the first half of Theorem \ref{thm:great}: Evaluating the Hilbert--Schmidt norm]
 Since the operator $H_{f,g,m}$ has a kernel given by the right-hand side of \eqref{eq:kernel-H}, one computes its Hilbert-Schmidt norm as
  \begin{align*}
  	\|H_{f,g,m}\|_{HS}^2
  		&= \iint_{{\mathbb{R}}^d\times{\mathbb{R}}^d} |H_{f,g,m}(x,\eta)|^2  \d x d\eta
  			= \iint_{{\mathbb{R}}^d\times{\mathbb{R}}^d} \left|f(x)g(\eta)-m(f(x)g(\eta))\right|^2  \, \frac{\d x d\eta}{(2\pi)^d} \\
  		&= \int_{{\mathbb{R}}^d} G_{g,m}(f(x))\, dx
  \end{align*}
  using the Fubini-Tonelli Theorem and the definition of $G_{g,m}$.
  The proof of the second half of Theorem \ref{thm:great} will be given later in Section \ref{sec:max-F-mult} below.
\end{proof}

In the rest of this section we will discuss how Theorem \ref{thm:great} and the bound \eqref{eq:first-upper-bound} easily lead to the Cwikel--Lieb--Rozenblum bound for a non-relativistic single-particle Schr\"odinger operator:
In this case
$g(\eta)=|\eta|^{-1}$. A simple scaling argument gives
\begin{align*}
	\|H_{f,g,m}\|_{HS}^2
		&= \iint_{{\mathbb{R}}^\times{\mathbb{R}}^d} \left(\frac{f(x)}{|\eta|}- m\left(\frac{f(x)}{|\eta|}\right)\right)^2 \frac{dxd\eta}{(2\pi)^d} \\
		&= \int_{{\mathbb{R}}^d} f(x)^d\, dx \int_{{\mathbb{R}}^d} (|\eta|^{-1}-m(|\eta|^{-1}))^2\, \frac{d\eta}{(2\pi)^d}
\end{align*}
Going to spherical coordinates shows
\begin{align*}
	\int_{{\mathbb{R}}^d} (|\eta|^{-1}-m(|\eta|^{-1}))^2\, \frac{d\eta}{(2\pi)^d}
		&= \frac{|S^{d-1}|}{(2\pi)^d} \int_0^\infty \left(r^{-1}-m(r^{-1})\right)^2 r^{d-1}\, dr \\
		&= \frac{ d|B_1^d|}{(2\pi)^d} \int_0^\infty (1-t^{-1}m(t))^2 t^{1-d}\, dt\, ,
\end{align*}
where $|S^{d-1}|$ is the surface area of the unit sphere in ${\mathbb{R}}^d$ and $|B_1^d|=|S^{d-1}|/d$ is the volume of the unit ball in ${\mathbb{R}}^d$.

In particular, scaling $f$ by  $\kappa>0$, using
$\kappa A_{f,g} =  A_{\kappa f,g}= B_{\kappa f,g,m}+ H_{\kappa f,g,m}$, the argument leading to \eqref{eq:first-upper-bound} now leads to
\begin{align}
	N(P^2-U)
	  & = n(A_{\kappa  f,g};\kappa)
	      \le (\kappa-\mu)^{-2}\sum_{j}\, \|H_{\kappa f,g,m}\|_{HS}^2 \label{eq:general-upper-bound-1} \\
	  &=  \frac{\kappa^d}{(\kappa-\mu)^{2}} \frac{ d|B_1^d|}{(2\pi)^d} \int_0^\infty (1-t^{-1}m(t))^2 t^{1-d}\, dt\, \int_{{\mathbb{R}}^d} U(x)^{d/2}\, dx \ ,  \label{eq:general-upper-bound-2}
\end{align}
as long as $\kappa> \mu\ge \|B_{\kappa f,g,m}\|$.
It is important to note here that the last factor on the right hand side of the above bound has the correct dependence on the potential $U$. That is, the factor in front of it, which depends on the upper bound $\mu$ on the operator norm of $B_{f,g,m}$, \emph{has to be independent} of $f=\sqrt{U}$.
Thanks to the second part of Theorem \ref{thm:great}, we can use
$\mu = \|m_1(s)\|_{L^2({\mathbb{R}}_+, \frac{ds}{s})} \|m_1(s)\|_{L^2({\mathbb{R}}_+, \frac{ds}{s})}$ as an upper bound for $\|B_{f,g,m}\|$, which is independent of $f$ (and $g$), so the same bound holds for $\|B_{\kappa f,g,m}\|$ for any $\kappa>0$. Using this, we can now freely optimize in $\kappa>\mu$ in \eqref{eq:general-upper-bound-2}, to get

\begin{align}\label{eq:super}
	N(P^2-U) \le C \frac{|B_1^d|}{(2\pi)^d} \int_{{\mathbb{R}}^d} U(x)^{d/2}\, dx
\end{align}
with the constant
\begin{equation}\label{eq:super-constant}
  \begin{split}
	C&=C_{d,m}= \frac{d^{d+1}}{4(d-2)^{d-2}}\mu^{d-2} \int_0^\infty (1-t^{-1}m(t))^2 t^{1-d}\, dt  \\
		&= \frac{d^{d+1}}{4(d-2)^{d-2}}\left(\|m_1\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})}\|m_2\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})}\right)^{d-2}\int_0^\infty (1-t^{-1}m(t))^2 t^{1-d}\, dt
  \end{split}
\end{equation}
where $m(t)= m_1*m_2(t)=\int_0^\infty m_1(t/s)m_2(s)\, \frac{ds}{s}$.

This gives most of the main ideas of our proof of Theorem \ref{thm:CLR}. The last new idea, which is crucially important for the proof of \eqref{eq:B-bound}, is the connection between the bound on the norm of the operator $B_{f,g,m}$ and bounds for maximal Fourier multipliers on $L^2$. This is explained in Section \ref{sec:max-F-mult}.

\begin{remark}\label{rem:simple-choice}
 In order to get good bounds on $C$, we have to find good candidates for $m_1$ and $m_2$. A simple, but not optimal, choice is $m_1(s)= s\id_{\{0<s\le 1\}}$ and $m_2(s)= s^{-1}\id_{\{s>1\}}$, in which case  $\|m_1\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})}=\|m_2\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})}=1$ and  $m(t)= m_1*m_2(t)= \min(t,t^{-1})$, so
 \begin{align*}
 	\int_0^\infty (1-t^{-1}m(t))^2 t^{1-d}\, dt
		= \int_1^\infty (1-t^{-2})^2 t^{1-d}\, dt
			= \frac{8}{(d-2)d(d+2)} \, .
 \end{align*}
 This gives
 \begin{align*}
 	C_{0,d}= \frac{2\, d^d}{(d-2)^{d-1}(d+2)}
 \end{align*}
  as a possible constant in the CLR inequality and yields $C_{0,3}\le 10.8$, already an order of a magnitude smaller than Cwikel's bound. Moreover, combining this with `stripping-off-dimensions' ideas with the help of similar bounds for operator-valued potentials, one can get this bound also uniformly in  the dimension.
To get the uniform bound claimed in Theorem \ref{thm:CLR} we have to choose better candidates for $m_1$ and $m_2$. We can achieve this in small dimensions and extend the bounds with the help of  bounds for operator-valued potentials, see Appendix \ref{sec:induction-in-dimension} and \ref{sec:numerics} and  Section \ref{sec:operator-valued}.
\end{remark}

Before we do this let us note a simple consequence of our approach for more general kinetic energies.

\section{General kinetic energies}\label{sec:gen-kin-energy}
First we consider the case where
$P^2$ is replaced by $P^{2\alpha}$ and give the
\begin{proof}[Proof of Theorem \ref{thm:poly}]
  Replacing $g(\eta)=|\eta|^{-1}$ by $g(\eta)=|\eta|^{-\alpha}$ one simply reruns the argument from the previous section. Calculating, again by scaling,
  \begin{align*}
	\|H_{f,g,m}\|_{HS}^2
		&= \iint_{{\mathbb{R}}^\times{\mathbb{R}}^d} \left(\frac{f(x)}{|\eta|^\alpha}- m\left(\frac{f(x)}{|\eta|^\alpha}\right)\right)^2 \frac{dxd\eta}{(2\pi)^d} \\
		&= \int_{{\mathbb{R}}^d} f(x)^{d/\alpha}\, dx \int_{{\mathbb{R}}^d} (|\eta|^{-\alpha}-m(|\eta|^{-\alpha}))^2\, \frac{d\eta}{(2\pi)^d}
\end{align*}
and
  \begin{align*}
	\int_{{\mathbb{R}}^d} (|\eta|^{-\alpha}-m(|\eta|^{-\alpha}))^2\, \frac{d\eta}{(2\pi)^d}
		&= \frac{|S^{d-1}|}{(2\pi)^d} \int_0^\infty \left(r^{-\alpha}-m(r^{-\alpha})\right)^2 r^{d-1}\, dr \\
		&= \frac{ d|B_1^d|}{\alpha (2\pi)^d} \int_0^\infty (1-t^{-1}m(t))^2 t^{1-\frac{d}{\alpha}}\, dt\, ,
  \end{align*}
  one sees that the argument leading to \eqref{eq:general-upper-bound-2} remains virtually unchanged, only  $d$ gets replaced by by $d/\alpha$. Thus
  \begin{align*}
  	N(P^{2\alpha}+V)\le C \frac{ d|B_1^d|}{\alpha (2\pi)^d} \int_{{\mathbb{R}}^d} V_-(x)^{\frac{d}{2\alpha}}\, dx
  \end{align*}
  with constant
  \begin{align*}
	C&=\frac{(\frac{d}{\alpha})^{\frac{d}{\alpha}+1}}{4(\frac{d}{\alpha}-2)^{\frac{d}{\alpha}-2}}
		\left(\|m_1\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})}\|m_2\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})}\right)^{\frac{d}{\alpha}-2}\int_0^\infty (1-t^{-1}m(t))^2 t^{1-\frac{d}{\alpha}}\, dt
  \end{align*}
  For $m_1$ and $m_2$ we make the simple choice from Remark \ref{rem:simple-choice}. Then  $m(t)=m_1*m_2(t)= \min(t,t^{-1})$ and
 $\mu=\|m_1\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})}\|m_2\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})}=1$. Hence,
 \begin{align*}
 	\int_0^\infty (1-t^{-1}m(t))^2 t^{1-\frac{d}{\alpha}}\, dt
 		= \int_1^\infty (1-t^{-2})^2 t^{1-\frac{d}{\alpha}}\, dt
 		= \frac{8}{(\frac{d}{\alpha}-2)\frac{d}{\alpha}(\frac{d}{\alpha}+2)}
 \end{align*}
 and collecting terms finishes the proof of Theorem \ref{thm:poly}.
\end{proof}
\begin{remark}
	For the number of negative energy bound states of $P^{2\alpha}+U$ the so-far best bounds are due to Frank \cite{Fr,FrReview}. Using ideas from Rumin \cite{Ru1,Ru2}, he got the bound
	\begin{align*}
		N(P^{2\alpha}+V)
			\le \left(
					\frac{(\frac{d}{\alpha}(\frac{d}{\alpha+2})}{(\frac{d}{\alpha}-2)^2}
				\right)^{\frac{d}{2\alpha}-1} \frac{\frac{d}{\alpha}}{\frac{d}{\alpha}-2}
			\frac{|B_1^d|}{(2\pi)^d} \int_{{\mathbb{R}}^d} V_-(x)^{\frac{d}{2\alpha}}\, dx\, .
	\end{align*}
Even with the non-optimal choice of $m_1$ and $m_2$ above, a simple calculation shows that the bound from Theorem \ref{thm:poly}  is better as long as $2< \left( 1+2\alpha/d\right)^{d/(2\alpha)}$. Since $0<\delta\mapsto \left( 1+1/\delta\right)^{\delta} $ is strictly increasing, this is the case as soon as  $d>2\alpha$, that is, the whole range of allowed values of $\alpha$.
\end{remark}

For more general kinetic energies of the form $T(P)$ with $T$ a suitable non-negative measurable function obeying some mild growth condition at infinity, we have:

\begin{theorem}\label{thm:gen-kin-energy}
	The number of negative energy bound states of a Schr\"odinger--type operator $T(P)+V$, defined suitably with the help of quadratic form methods on $L^2$, obeys the  bound
	\begin{align}
		N(T(P)+V) \le \lambda^{-2} \int_{{\mathbb{R}}^d} G_{T}\big((\lambda+1)^2 V_-(x)\big)\, \d x
	\end{align}
	for any $\lambda>0$, with
	\begin{align*}
		G_{T}(u) = \int_{\frac{u}{T}> 1} \left[\Big(\frac{u}{T(\eta)}\Big)^{1/2}- \Big(\frac{u}{T(\eta)}\Big)^{-1/2}\right]^2\, \frac{\d \eta}{(2\pi)^d}
	\end{align*}
	and $ V_-=\max(-V,0)$, the negative part of $V$.
\end{theorem}
\begin{proof}
In this case we use $g(\eta)= T(\eta)^{-1/2}$, $f(x)= V_-(x)$, and make again the  the choice
 $m_1(s)=2s\id_{\{0<s\le 1\}}$ and $m_2(s)= s^{-1}\id_{s\ge 1}$. So
 $\mu=\|m_1\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})}\|m_2\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})}=1$ and with $\lambda=\kappa-\mu=\kappa-1$, the same argument leading to  \eqref{eq:general-upper-bound-1} now gives
 \begin{align*}
 	N(T(P)+V)\le N(T(P)-V_-) \le \lambda^{-2} \ \|H_{(\lambda+1)f,g,m}\|_{HS}^2\, .
 \end{align*}
 for any $\lambda>0$.
 Using Theorem \ref{thm:great} to calculate the Hilbert--Schmidt norm shows
 \begin{align*}
 	\|H_{(\lambda+1)f,g,m}\|_{HS}^2 = \int_{{\mathbb{R}}^d} G_{T}\big((\lambda+1)^2 V_-(x)\big)\, \d x \ ,
 \end{align*}
 since $m(t)=m_1*m_2(t)=\min(t,t^{-1})$.
\end{proof}
\begin{remarks}
\refstepcounter{theoremi} {\rm(\roman{theoremi})}{~} The bound given in Theorem \ref{thm:gen-kin-energy} improves the bound from \cite{HHRW}, which was based on Cwikel's original method.
  The result proven in \cite{HHRW} shows that under some rather mild general conditions on the kinetic energy symbol $T$ the operator
  $T(P)+V$ has weakly coupled bound states for any non-trivial potential $V$, no matter how small, if $\eta\mapsto \frac{1}{T(\eta)}$ is
  not integrable over the set $\{T<u\}$ for all small $u>0$, which is
  equivalent to $G_T(u)=\infty$ for all small $u>0$. This shows that
  the bound given by Theorem \ref{thm:CLR} is quite natural.  \\[0.2em]
\refstepcounter{theoremi} {\rm(\roman{theoremi})}{~} In all cases where one can explicitly calculate or find explicit upper bounds for $G_T$, Theorem \ref{thm:gen-kin-energy} gives an upper bound of the form
  \begin{align*}
  	N(T(P)+V)\lesssim N^\text{cl}(T +V)
  \end{align*}
  with $ N^\text{cl}(T+V) $ the classical phase-space volume of the set where $H(\eta,x)=T(\eta)+V(x)<0$. See \cite{Solomyak,BL,BLS} and the discussion in Section 6 of \cite{HHRW}, where it is also shown that logarithmic corrections to the classical phase-space guess appear in critical cases.
  \end{remarks}
\section{The connection with maximal Fourier multipliers}
\label{sec:max-F-mult}
In this section we finish the proof of Theorem \ref{thm:great}. Recall that given functions $f,g:{\mathbb{R}}^d\to[0,\infty)$ and a bounded function $m:{\mathbb{R}}_+\to{\mathbb{R}}_+ $, the operator $B_{f,g,m}$ is given by
\begin{align}
	B_{f,g,m}\varphi(x) = (2\pi)^{-d/2} \int_{{\mathbb{R}}^d} e^{ix\eta} m(f(x)g(\eta)) \varphi(\eta)\, d\eta \ ,
\end{align}
at least for nice enough $\varphi$, say from the set of Schwartz functions. We would like to conclude that $B_{f,g,m}$ is a bounded operator on $L^2({\mathbb{R}}^d)$, which might suggest to look for results when a pseudo-differential operator with symbol $a(x,\eta)= m(f(x)g(\eta))$ is bounded. But such results need enough differentiability of the symbol $a$, which we do not have. More importantly, we need an estimate independent of $f$, which one cannot get without looking more closely into the structure of the problem.
To see how the product structure $f(x) g(\eta)$ helps in the operator bound, we rewrite $B_{f,g,m}$ as
\begin{equation}
  \begin{split}
	B_{f,g,m}\varphi(x)
		&= (2\pi)^{-d/2} \int_{{\mathbb{R}}^d} e^{ix\eta} m(tg(\eta)) \varphi(\eta)\, d\eta \, \Big\vert_{t=f(x)} \\
		&= \mathcal{F}^{-1}\left[ m(tg(\cdot))\varphi(\cdot) \right](x)\Big\vert_{t=f(x)}
  \end{split}
\end{equation}
This suggest to look at the Fourier multiplier $B_{t,g,m}$ defined by
\begin{align}\label{eq:F-mult}
	B_{t,g,m}\varphi \coloneqq \mathcal{F}^{-1}\left[ m(tg(\cdot))\varphi(\cdot) \right]
\end{align}
and the associated maximal operator\footnote{To be more careful, one should take the supremum over a dense subset of ${\mathbb{R}}_+$, to ensure measurablity, but for $\varphi$ in the Schwarz class this makes no difference}
\begin{align}\label{eq:max-F-mult}
	B^*_{g,m}\varphi(x)\coloneqq \sup_{t>0}|B_{t,g,m}\varphi (x)|\,.
\end{align}
It is clear that, as operators on $L^2$, one has $\|B_{f,g,m}\|\le \|B^*_{g,m}\|$ for the corresponding operator norms.
On the other hand, choosing $f(x)$ in such a way as to make $|B_{f,g,m}\varphi(x)|$ arbitrarily close to $B^*_{g,m}\varphi(x)$, shows the `reverse bound'
$\sup_{f\ge 0}\|B_{f,g,m}\varphi\|_2 \ge \|B^*_{g,m}\varphi\|_2$ for a given fixed Schwartz function $\varphi$.
 Thus $\sup_{f\ge 0}\|B_{f,g,m}\|= \|B^*_{g,m}\|$, so having a bound for $B_{f,g,m}$ which is \emph{uniform} in the choice of the function  $f$ is \emph{ equivalent} to having a bound for the maximal Fourier multiplier $B^*_{g,m}$.

 This is our starting point for the proof of the second half of Theorem \ref{thm:great}.
\begin{theorem}\label{thm:max-op-bound}
	Let $g$ be a measurable non-negative function on ${\mathbb{R}}^d$ and assume that $m:{\mathbb{R}}_+\to{\mathbb{R}} $ is given by a convolution,
	\begin{align*}
		m(t)= m_1*m_2(t)= \int_0^\infty m_1(t/s)m_2(s)\frac{ds}{s}
	\end{align*}
	with $m_1,m_2\in L^2({\mathbb{R}}_+,\frac{ds}{s})$.
	Then the maximal Fourier multiplier $B^*_{g,m}$, defined in \eqref{eq:max-F-mult}, is bounded  on $L^2({\mathbb{R}}^d)$ with bound
	\begin{align*}
		\|B^*_{g,m}\| \le \|m_1\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})} \|m_2\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})}
	\end{align*}
	for its operator norm.
\end{theorem}

\begin{remark}
	There are several different but related proofs of boundedness of maximal Fourier multipliers available in the literature, see, e.g., \cite{Carbery,DappaTrebels,Rubio}. These works concentrate on getting $L^p$ bounds and do not care much about the involved constants. For us the $L^2$ boundedness is important, with good bounds on the operator norm.
\end{remark}

\begin{proof}
  The proof is easy and uses nothing more than the Cauchy-Schwarz inequality: since $m$ is given by a convolution type integral, we have
  \begin{align*}
  	B_{t,g,m}\varphi(x)
  		& =	\int_0^\infty \mathcal{F}^{-1}\left[ m_1(tg/s)\varphi \right](x) \, m_2(s)\, \frac{ds}{s}.
  \end{align*}
Interchanging the integrals is certainly fine for nice enough $\varphi$ from a dense subset of $L^2$, say the Schwartz functions. Applying the triangle and then the Cauchy-Schwarz inequality for the $ds/s$ integration yields
 \begin{align}
   |B_{t,g,m}\varphi(x)|
  		& \le 	\int_0^\infty \left|\mathcal{F}^{-1}\left[ m_1(tg/s)\varphi \right](x)\right| \, |m_2(s)|\, \frac{ds}{s}\nonumber\\
  		&\le \left( \int_0^\infty \left|\mathcal{F}^{-1}\left[ m_1(tg/s)\varphi \right](x)\right|^2 \frac{ds}{s} \right)^{1/2}
  			\|m_2\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})}\label{eq:punch}
 \end{align}
 The punchline is that the measure $ds/s$ is invariant under scaling, so we can scale $s$ by a fixed factor $t$ to see that
 \begin{align*}
 	\int_0^\infty \left|\mathcal{F}^{-1}\left[ m_1(tg/s)\varphi \right](x)\right|^2 \frac{ds}{s}
 	= \int_0^\infty \left|\mathcal{F}^{-1}\left[ m_1(g/s)\varphi \right](x)\right|^2 \frac{ds}{s}\, ,
 \end{align*}
 that is, the right hand side of \eqref{eq:punch} is independent of $t>0$.
 So
 \begin{align*}
 	B^*_{g,m}\varphi(x) = \sup_{t>0}|B_{t,g,m}\varphi(x)|
 		\le \left( \int_0^\infty \left|\mathcal{F}^{-1}\left[ m_1(g/s)\varphi \right](x)\right|^2 \frac{ds}{s} \right)^{1/2}
  			\|m_2\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})}.
 \end{align*}
 In particular,
 \begin{align*}
 	\|B^*_{g,m}\varphi\|_2^2
 		&\le  \|m_2\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})}^2
 				\int_{{\mathbb{R}}^d} \int_0^\infty \left|\mathcal{F}^{-1}\left[ m_1(g/s)\varphi \right](x)\right|^2 \, \frac{ds}{s}\, dx\, .
 \end{align*}
Using Fubini--Tonelli to interchange the integrals and  Plancherel's theorem for the $L^2$ norm of the Fourier transform, one sees that
\begin{align*}
	\int_{{\mathbb{R}}^d} \int_0^\infty & \left|\mathcal{F}^{-1}\left[ m_1(g/s)\varphi \right](x)\right|^2 \, \frac{ds}{s}\, dx\,
		= \int_0^\infty \int_{{\mathbb{R}}^d}  |m_1(g(\eta)/s)|^2|\varphi(\eta)|^2 \, d\eta\,  \frac{ds}{s}\, \\
	&=  \int_0^\infty \int_{{\mathbb{R}}^d}  |m_1(s^{-1})|^2|\varphi(\eta)|^2 \, d\eta\,  \frac{ds}{s}\,
		= \| m_1\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})}^2
			\|\varphi(\eta)\|_2^2
\end{align*}
where we also used the same scaling argument\footnote{Strictly speaking one gets an inequality if $g(\eta)=0$ since  $m_1(0)=0$.} as before to scale out $g(\eta)$ and the invariance of $ds/s$ under inversion $s\mapsto s^{-1}$. Thus
\begin{align*}
	\|B^*_{g,m}\varphi\|_2\le \| m_1\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})} \| m_2\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})} \|\varphi\|_2
\end{align*}
and we are done.
\end{proof}
The next result, which also proves the second half of Theorem \ref{thm:great}, is an easy consequence of Theorem \ref{thm:max-op-bound}.
\begin{corollary}
	Let $f,g$ be measurable non-negative functions on ${\mathbb{R}}^d$ and assume that $m:{\mathbb{R}}_+\to{\mathbb{R}} $ is given by a  convolution,
	\begin{align*}
		m(t)= m_1*m_2(t)= \int_0^\infty m_1(t/s)m_2(s)\frac{ds}{s}
	\end{align*}
	with $m_1,m_2\in L^2({\mathbb{R}}_+,\frac{ds}{s})$.
	Then the operator $B_{f,g,m}$, defined by \eqref{eq:kernel-B}, i.e.,  given by the kernel
	\begin{align*}
		B_{f,g,m}(x,\eta) &= (2\pi)^{-d/2} e^{ix\cdot\eta}m(f(x)g(\eta))
	\end{align*}
	is bounded  on $L^2({\mathbb{R}}^d)$ with bound
	\begin{align*}
		\sup_{g\ge0}\big\|\sup_{f\ge0}|B_{f,g,m}\varphi|\big\|_2 \le \|m_1\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})} \|m_2\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})}\|\varphi\|_2 \, .
	\end{align*}
\end{corollary}
\begin{proof} By definition of the maximal Fourier multiplier we have  $|B_{f,g,m}\varphi(x)|\le B^*_{g,m}\varphi(x)$ and thus also
$\sup_{f\ge 0}|B_{f,g,m}\varphi(x)|\le B^*_{g,m}\varphi(x)$ for almost every $x\in{\mathbb{R}}^d$.

Since the $L^2$--bound from Theorem \ref{thm:max-op-bound} is independent of $g\ge 0$, we can also take the supremum in $g\ge 0$ , after taking the $L^2$--norm.
\end{proof}

\section{A lower bound for the variational problem}
\label{sec:lower-bound}
Recall that the variational problem, which comes up in a natural way in our bound on the number of bound states is
\begin{align}
	P_\gamma= \inf \left\{ (\|m_1\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})}\|m_2\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})})^{\gamma-2} \int_0^\infty (1-t^{-1} m_1*m_2(t))^2 t^{1-\gamma}\,dt \right\},
\end{align}
where the convolution $m_1*m_2$ is on ${\mathbb{R}}_+$ with its scaling invariant measure $\frac{ds}{s}$, and the infimum is taken over all functions $m_1,m_2:{\mathbb{R}}_+\to{\mathbb{R}} $ .

\begin{theorem}
We have the lower bound
  \begin{align*}
  	P_\gamma\ge \frac{2}{(\gamma-2)(\gamma-1)\gamma}
  \end{align*}
  for all $\gamma>2$.
\end{theorem}
\begin{proof}
  The proof is straightforward, once one notices that, for example by Cauchy--Schwarz, one has $\|m\|_\infty \le \|m_1\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})}\|m_2\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})}$  for $m =m_1*m_2$.
  Thus
  \begin{align*}
  	P_\gamma &\ge \inf_{m} \left\{ \|m\|_\infty^{\gamma-2} \int_0^\infty (1-t^{-1} m(t))^2 t^{1-\gamma}\,dt \right\} \\
  	&= \inf_{l>0}\inf_{\|m\|_\infty=l} \left\{ l^{\gamma-2} \int_0^\infty (1-t^{-1} m(t))^2 t^{1-\gamma}\,dt \right\}
  \end{align*}
 To make $\int_0^\infty (1-t^{-1} m(t))^2 t^{1-\gamma}\,dt$ as small as possible under the constraint $\|m\|_\infty=l$, one simply chooses
 $m(t)=\min(t,l)$. Then
 \begin{align*}
 	\int_0^\infty (1-t^{-1} m(t))^2 t^{1-\gamma}\,dt
 	= \int_l^\infty (1-t^{-1}l)^2 t^{1-\gamma}\,dt
 	= l^{2-\gamma} \frac{2}{(\gamma-2)(\gamma-1)\gamma}
 \end{align*}
 and one sees that
 \begin{align*}
 	\inf_{\|m\|_\infty=l} \left\{ l^{\gamma-2} \int_0^\infty (1-t^{-1} m(t))^2 t^{1-\gamma}\,dt \right\}
 	= \frac{2}{(\gamma-2)(\gamma-1)\gamma}.
 \end{align*}
\end{proof}

\section{Extension to operator--valued potentials}\label{sec:operator-valued}
In this section we extend our method to operator--valued potentials and give the proof of Theorem \ref{thm:poly-operator-valued}, i.e. we prove that the number of negative bound states of $P^{2\alpha}\otimes \idG + V$ is bounded by
 \begin{align*}
 	N(P^{2\alpha}\otimes\idG+V)
			\le C_{d/\alpha}\,
			\frac{|B_1^d|}{(2\pi)^d} \int_{{\mathbb{R}}^d} \mathop{\mathrm{tr}} \nolimits_\mathcal{G}[ V_-(x)^{\frac{d}{2\alpha}}]\, dx \ ,
 \end{align*}
 where $V:{\mathbb{R}}^d\to \mathcal{B}(\mathcal{G})$ is an operator valued potential with positive part $V_+\in L^1_{\text{loc}}({\mathbb{R}}^d,\mathcal{B}(\mathcal{G}))$ and negative part  $V_-\in L^{d/(2\alpha)}({\mathbb{R}}^d, {\mathcal{S}}_{d/(2\alpha)}(\mathcal{G}))$.

 Let $U(x)= V(x)_-$, the negative part of $V(x)$ defined by spectral calculus.
The Birman--Schwinger operator corresponding to $|P|^{2\alpha}\otimes\idG -U$ is given by
\begin{align*}
	K= \sqrt{U} (|P|^{-2\alpha}\otimes\idG) \sqrt{U}
\end{align*}
and we again have
\begin{align*}
	N(|P|^{2\alpha}\otimes\idG +V) \le N(|P|^{2\alpha}\otimes\idG -U)
		= n(K; 1).
\end{align*}
Now we factor $K$ as $K= \widetilde  {A}_{f,g}^* \widetilde  {A}_{f,g}$ where $\widetilde  {A}_{f,g}$ has kernel
\begin{align*}
	\widetilde  {A}_{f,g}\varphi (\eta)
		= (2\pi)^{-d/2} \int_{{\mathbb{R}}^d} e^{-i\eta \cdot x} g(\eta) f(x)\varphi(x)\, dx \ ,
\end{align*}
$g(\eta)= |\eta|^{-\alpha}$ is real--valued (even positive), and $f(x)=\sqrt{U(x)}$ takes values in the self-adjoint positive operators on $\mathcal{G}$. We split this as
\begin{align*}
	\widetilde  {A}_{f,g}= \widetilde  {B}_{f,g,m} + \widetilde  {H}_{f,g,m}
\end{align*}
with a function $m:[0,\infty)\to {\mathbb{R}} $, so that
\begin{align}\label{eq:kernel-B-operator-valued}
	\widetilde  {B}_{f,g,m}\varphi(\eta) &= (2\pi)^{-d/2} \int_{{\mathbb{R}}^d} e^{-i\eta \cdot x} m(g(\eta) f(x))\varphi(x)\, dx
		= \mathcal{F}\left[ m(tf )\varphi \right](\eta)\Big|_{t=g(\eta)}
\end{align}
and
\begin{align}\label{eq:kernel-H-operator-valued}
	\widetilde  {H}_{f,g,m}\varphi(\eta) &= (2\pi)^{-d/2} \int_{{\mathbb{R}}^d} e^{-i\eta \cdot x} \left[g(\eta)f(x) -m(g(\eta) f(x))\right]\varphi(x)\, dx \ ,
\end{align}
where $\varphi$ is a function from a nice dense subset of $L^2({\mathbb{R}}^d,\mathcal{G})$, so that the integrals converge and $m(t f(x))$ is an operator on $\mathcal{G}$ defined via functional calculus.
\begin{remark}
  With a slight abuse of notation, we write $\mathcal{F}$ in the definition of $\widetilde  {B}_{f,g,m}$, which strictly speaking denotes the Fourier transform on $L^2({\mathbb{R}}^d)$, instead of $\mathcal{F}\otimes\idG$, the Fourier transform on $L^2({\mathbb{R}}^d,\mathcal{G}) = L^2({\mathbb{R}}^d)\otimes\mathcal{G}$.  In addition, in the definition of $\widetilde  {B}_{f,g,m}$ and
	$\widetilde  {H} _{f,g,m}$ above we swapped the role of $f$ and $g$ compared to the discussion in Section \ref{sec:max-F-mult}. This is convenient, since by assumption $g(\eta)$ is a  multiplication operator  on $\mathcal{G}$, and this makes a maximal Fourier multiplier estimate, now with $g$ instead of $f$, easier.
	The general case can be reduced to this setting, see Section \ref{sec:trace-ideals} below.
\end{remark}

The following theorem is the replacement of Theorem \ref{thm:great} in the operator-valued setting.

\begin{theorem}\label{thm:great-operator-valued}
	 $\widetilde  {H}_{f,g,m}$ is  Hilbert--Schmidt operator on $\mathcal{H}=L^2({\mathbb{R}}^d,\mathcal{G})$ with Hilbert--Schmidt norm given by
	\begin{align}\label{eq:operator-valued-HS-norm}
		\|\widetilde  {H}_{f,g,m}\|_{{\mathcal{S}}_2(\mathcal{H})}^2
		= \int_{{\mathbb{R}}^d} \mathop{\mathrm{tr}} \nolimits_\mathcal{G}\left[G_{g,m}(f(x))\right]\, dx \ ,
	\end{align}
	where $G_{g,m}$ is again given by
	\begin{align}
		G_{g,m}(u)= \int_{{\mathbb{R}}^d} |ug(\eta)-m(ug(\eta))|^2 \frac{\d \eta}{(2\pi)^d} .
	\end{align}
    If, moreover, $m=m_1*m_2$
    then for all measurable non-negative functions $g$ and non-negative operator-valued functions $f$ the operator $\widetilde  {B}_{f,g,m}$ is bounded on $\mathcal{H}$ with
    \begin{align}
    	\|\widetilde  {B}_{f,g,m}\varphi\|_{\mathcal{H}}
    		\le  \|m_1\|_{L^2({\mathbb{R}}_+, \frac{ds}{s})}
    			\|m_2\|_{L^2({\mathbb{R}}_+, \frac{ds}{s})}
    				\|\varphi\|_{\mathcal{H}}
    \end{align}
    for all $\varphi\in \mathcal{H}$.
\end{theorem}
\begin{proof}
  To prove \eqref{eq:operator-valued-HS-norm}, we note that the Hilbert--Schmidt operators on $\mathcal{H}= L^2({\mathbb{R}}^d,\mathcal{G})$ are isomorphic to operators with kernels in $L^2({\mathbb{R}}^d\times{\mathbb{R}}^d, {\mathcal{S}}_2(\mathcal{G}))$ and
  \begin{align*}
  	\|\widetilde  {H}\|_{{\mathcal{S}}_2(\mathcal{H})}^2 &= \mathop{\mathrm{tr}} \nolimits_{\mathcal{H}}\big[ \widetilde  {H}^* \widetilde  {H} \big]
  		= \iint_{{\mathbb{R}}^d\times{\mathbb{R}}^d} \|\widetilde  {H}(\eta,x)\|_{{\mathcal{S}}_2(\mathcal{G})}^2 \, d x d\eta \ ,
  \end{align*}
see Lemma \ref{lem:hilbert-schmidt}.

	Using the explicit form of the `kernel' of $\widetilde  {H}_{f,g,m}$ given in \eqref{eq:kernel-H-operator-valued} this shows
 \begin{align*}
 	\|\widetilde  {H}\|_{{\mathcal{S}}_2(\mathcal{H})}^2
 		&= (2\pi)^{-d}\int_{{\mathbb{R}}^d}\int_{{\mathbb{R}}^d} \mathop{\mathrm{tr}} \nolimits_\mathcal{G}\big[ |g(\eta)f(x)- m(g(\eta f(x)))|^2 \big]\, d\eta d x \\
 		&= \int_{{\mathbb{R}}^d} \mathop{\mathrm{tr}} \nolimits_\mathcal{G}\big[ G_{g,m}(f(x)) \big]\, d x
 \end{align*}
 by the definition of $G_{g,m}$ and the spectral theorem.

 Concerning the boundedness of $\widetilde  {B}_{f,g,m}$ we recall \eqref{eq:kernel-B-operator-valued} and, if $m=m_1*m_2$,
 \begin{align*}
 	\widetilde  {B}_{f,t,m}\varphi(\eta)
 		\coloneqq \mathcal{F}\left[ m(tf)\varphi \right](\eta)
 		= \int_0^\infty \mathcal{F}\left[ m_1(f /s)\varphi \right](\eta) \, m_2(ts)\, \frac{ds}{s}.
 \end{align*}
 Thus,
 \begin{align*}
 	\big\| \widetilde  {B}_{f,t,m}\varphi(\eta) \big\|_\mathcal{G}
 		&\le  \int_0^\infty \big\| \mathcal{F}\left[ m_1(f /s)\varphi \right](\eta) \big\|_\mathcal{G}\, |m_2(ts)|\, \frac{ds}{s} \\
 		&\le \left( \int_0^\infty \big\| \mathcal{F}\left[ m_1(f /s)\varphi \right](\eta) \big\|_\mathcal{G}^2 \,\frac{ds}{s}\right)^{1/2}
 			\left( \int_0^\infty |m_2(ts)|^2\, \frac{ds}{s} \right)^{1/2} \\
 		&= \left( \int_0^\infty \big\| \mathcal{F}\left[ m_1(f /s)\varphi \right](\eta) \big\|_\mathcal{G}^2 \,\frac{ds}{s}\right)^{1/2} \|m_2\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})}
 \end{align*}
 due to the scaling invariance of $ds/s$. We therefore have a maximal operator bound
 \begin{align*}
 	\widetilde  {B}_{f,m}^* \varphi(\eta)
 		\coloneqq \sup_{t>0}
 					\big\| \widetilde  {B}_{f,t,m}\varphi(\eta) \big\|_\mathcal{G}
 					\le \left( \int_0^\infty \big\| \mathcal{F}\left[ m_1(f /s)\varphi \right](\eta) \big\|_\mathcal{G}^2 \,\frac{ds}{s}\right)^{1/2} \|m_2\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})} \ .
 \end{align*}
In particular,
 \begin{align*}
 	\|\widetilde  {B}_{f,m}^* \varphi\|_{L^2({\mathbb{R}}^d)}^2
 	\le \|m_2\|_{L^2({\mathbb{R}}_+.\frac{ds}{s})}^2 \int_{{\mathbb{R}}^d} \int_0^\infty \big\| \mathcal{F}\left[ m_1(f /s)\varphi \right](\eta) \big\|_\mathcal{G}^2 \,\frac{ds}{s}\, d\eta \ ,
 \end{align*}
and
 \begin{align*}
 	\int_{{\mathbb{R}}^d} \int_0^\infty &\big\| \mathcal{F}\left[ m_1(f /s)\varphi \right](\eta) \big\|_\mathcal{G}^2 \,\frac{ds}{s}\, d\eta
 	 = \int_0^\infty \int_{{\mathbb{R}}^d} \left\langle \mathcal{F}[m_1(f/s)\varphi](\eta), \mathcal{F}[m_1(f/s)\varphi](\eta) \right\rangle_\mathcal{G} \, d\eta\, \frac{d s}{s} \\
 	&= \int_0^\infty \int_{{\mathbb{R}}^d} \left\langle m_1(f(x)/s)\varphi(x), m_1(f(x)/s)\varphi(x) \right\rangle_\mathcal{G} \, d x\, \frac{d s}{s} \\
 	&= \int_{{\mathbb{R}}^d} \left\langle \varphi(x), \int_0^\infty m_1(f(x)/s)^2\, \frac{ds}{s} \varphi(x) \right\rangle_\mathcal{G}\, dx \\
 	 &= \|m_1\|_{L^2({\mathbb{R}}_+, \frac{ds}{s})}^2 \int_{{\mathbb{R}}^d} \|\varphi(x)\|_\mathcal{G}^2\, dx
 	 	= \|m_1\|_{L^2({\mathbb{R}}_+, \frac{ds}{s})}^2
 	 		\|\varphi\|_\mathcal{H}^2
 \end{align*}
 where we again used that, by scaling $\int_0^\infty m_1(r/s)^2\, \frac{ds}{s}= \|m_1\|_{L^2({\mathbb{R}}_+, \frac{ds}{s})}^2$ for all $r>0$, so by functional calculus
 \begin{align*}
	\int_0^\infty m_1(f(x)/s)^2\, \frac{ds}{s}
	= \|m_1\|_{L^2({\mathbb{R}}_+, \frac{ds}{s})}^2\id_\mathcal{G}\, .
 \end{align*}
 Altogether, we get the operator-valued version of our previous maximal Fourier multiplier bound in the form
 \begin{align*}
 	\|\widetilde  {B}_{f,m}^* \varphi\|_{L^2({\mathbb{R}}^d)}^2
 	\le \|m_1\|_{L^2({\mathbb{R}}_+, \frac{ds}{s})}
 		\|m_2\|_{L^2({\mathbb{R}}_+, \frac{ds}{s})}
 	 		\|\varphi\|_\mathcal{H} \ ,
 \end{align*}
 and it is easy to see that
 \begin{align*}
 	\|\widetilde  {B}_{f,g,m}\varphi\|_{\mathcal{H}} \le \|\widetilde  {B}_{f,m}^* \varphi\|_{L^2({\mathbb{R}}^d)} \ ,
 \end{align*}
 which completes the proof of Theorem \ref{thm:great-operator-valued}.
\end{proof}

The proof of Theorem \ref{thm:poly-operator-valued} is now easy: one simply does the same steps as in the scalar case with \eqref{eq:general-upper-bound-1} replaced by
\begin{align*}
		N(P^{2\alpha}\otimes\idG-U)
	  & = n(\widetilde  {A}_{\kappa f, g};\kappa)
	      \le (\kappa-\mu)^{-2}\sum_{j}\, \|\widetilde  {H}_{\kappa f,g,m}\|_{{\mathcal{S}}_2(\mathcal{H})}^2 \ ,
\end{align*}
where now $\mu \ge \|\widetilde  {B}_{\kappa f,g,m}\varphi\|_{\mathcal{H}}$.
As before, Theorem \ref{thm:great-operator-valued} gives a $\kappa$-independent bound for $\|\widetilde  {B}_{\kappa f,g,m}\varphi\|_{\mathcal{H}}$,  in particular, we can take any  $\mu\ge \|m_1\|_{L^2({\mathbb{R}}_+, \frac{ds}{s})} \|m_2\|_{L^2({\mathbb{R}}_+, \frac{ds}{s})} $. It also allows us to calculate the Hilbert-Schmidt norm.
For $g(\eta)=|\eta|^{-\alpha}$ we get
\begin{align*}
	G_{g,m}(u) = u^{d/\alpha} \int_{{\mathbb{R}}^d} (|\eta|^{-\alpha}-m(|\eta|^{-\alpha}))^2\, \frac{d\eta}{(2\pi)^d} \ ,
\end{align*}
so
\begin{align*}
	\|\widetilde  {H}_{\kappa f,g,m}\|_{{\mathcal{S}}_2(\mathcal{H})}^2
	= \kappa^{d/\alpha} \int_{{\mathbb{R}}^d} (|\eta|^{-\alpha}-m(|\eta|^{-\alpha}))^2\, \frac{d\eta}{(2\pi)^d}
		\int_{{\mathbb{R}}^d} \mathop{\mathrm{tr}} \nolimits_\mathcal{G}\left[ f(x)^{d/\alpha} \right]\, dx \ .
\end{align*}
Using this in the above bound for $N(P^{2\alpha}\otimes\idG-U)$ and
minimizing over $\kappa$, as in the scalar case, finishes the proof of Theorem \ref{thm:poly-operator-valued}.

\section{Trace ideal estimates}\label{sec:trace-ideals}

In this section we show how the ideas developed so far give a simple proof of a fully operator-valued version of Cwikel's theorem. Such an inequality was first proved in \cite{Fr}.

In this setting let $(X,dx)$ and $(Y,dy)$ be sigma-finite measure spaces and $\mathcal{H}, \mathcal{G}$ (separable) Hilbert spaces. We denote by $L^p(X,{\mathcal{S}}_p(\mathcal{H}))$ the set of measurable functions $f: X\to {\mathcal{S}}_p(\mathcal{H})$, where ${\mathcal{S}}_p(\mathcal{H})$ is the space of $p$-summable compact operators, i.e. the von Neumann--Schatten class, on $\mathcal{H}$, such that
\begin{align*}
	\|f\|_{L^p(X,{\mathcal{S}}_p(\mathcal{H}))}^p \coloneqq \int_{X} \|f(x)\|_{{\mathcal{S}}_p(\mathcal{H})}^p\, dx <\infty \ .
\end{align*}
Similarly, we denote by $L^p_\text{w}(Y, \mathcal{B}(\mathcal{G}))$ the set of of all measurable functions $g:Y\to \mathcal{B}(\mathcal{G})$, with values in the bounded operators on $\mathcal{G}$, such that
\begin{align*}
	\|g\|_{L^p_\text{w}(Y, \mathcal{B}(\mathcal{G}))}^p
		\coloneqq \sup_{t>0} t^p\left|\left\{ y\in Y:\, \|g(y)\|_{\mathcal{B}(\mathcal{G})>t}\right\} \right| <\infty
\end{align*}
A map $A: L^2(X,\mathcal{H})\to L^2(Y,\mathcal{G})$ is in the weak trace--ideal ${\mathcal{S}}_{p,\text{w}}= {\mathcal{S}}_{p,\text{w}}(L^2(X,\mathcal{H}), L^2(Y,\mathcal{G}))$ if
\begin{align}
	\left\| f\Phi^* g \right\|_{p,\text{w}}
		\coloneqq \Big(\sup_{n\in{\mathbb{N}}} n^p s_n(A)\Big)^{1/p}
\end{align}
where $s_n(A)$ are teh singular--values of $A$, i.e. the eigenvalues of $A^* A:L^2(X,\mathcal{H})\to L^2(X,\mathcal{H})$.

\begin{theorem}[Fully operator valued version of Cwikel's theorem]\label{thm:Cwikel}
  Let $\Phi: L^2(X,\mathcal{H})\to L^2(Y,\mathcal{G})$ be a unitary operator,
  which is also bounded from  $L^1(X,\mathcal{H})$ into $L^\infty(Y,\mathcal{G})$.
  If $p>2$ and $f\in L^p(X,{\mathcal{S}}_p(\mathcal{H}))$ and $g\in L^p_\text{w}(Y, \mathcal{B}(\mathcal{G}))$, then $f\Phi^* g$ is in the weak trace ideal
  ${\mathcal{S}}_p(L^2(X,\mathcal{H}),L^2(Y,\mathcal{G}))$ and
  \begin{align*}
  \begin{split}
  	\left\| f\Phi^* g \right\|_{p,\text{w}}^p
  		&\coloneqq \sup_{n\in{\mathbb{N}}} n^p s_n(f\Phi^* g) \\
  		&\le \frac{2(p-2)}{p+2} \, \left(\frac{p}{p-2} \right)^{p} \|\Phi\|_{L^1\to L^\infty}^2 \|f\|_{L^p(X,{\mathcal{S}}_p(\mathcal{H}))}^p \|g\|_{L^p_\text{w}(Y, \mathcal{B}(\mathcal{G}))}^p \,.
  \end{split}
  \end{align*}
\end{theorem}

\begin{remark}
	Theorem \ref{thm:Cwikel} improves the result of Frank in \cite{Fr},
	\begin{align*}
		\left\| f\Phi^* g \right\|_{p,\text{w}}^p
			\le \frac{p}{2} \left(\frac{p}{p-2} \right)^{p-1}
				\|\Phi\|_{L^1\to L^\infty}^2 \|f\|_{L^p(X,{\mathcal{S}}_p(\mathcal{H}))}^p \|g\|_{L^p_\text{w}(Y, \mathcal{B}(\mathcal{G}))}^p \, ,
	\end{align*}
 by a factor of $(p+2)/2$. In addition, his bound in the scalar case, when $\Phi$ is the usual Fourier-transform,  is worse than Theorem \ref{thm:Cwikel} by a factor of $\frac{1}{2}(1+p/2)^{p/2}>1$ in the allowed range $p>2$.
\end{remark}

\begin{proof}
 First we note that one can reduce the result to the case when $g$ is pointwise a positive multiple of the identity operator on $\mathcal{G}$. As operators on $\mathcal{G}$ one has  $g(y) g(y)^* \le \|g(y)\|_{\mathcal{B}(\mathcal{G})}^2\id_{\mathcal{G}}$. Thus
  with $A_1= f\Phi^* g$ we have
  \begin{align*}
  	A_1 A_1^* = f\Phi^* g g^* \Phi f^* \le  f\Phi^* (\|g\|_{\mathcal{B}(\mathcal{G})}\id_{\mathcal{G}})^2 \Phi f^* = A_2 A_2^*
  \end{align*}
 with $A_2= f\Phi^* \|g\|_{\mathcal{B}(\mathcal{G})}\id_{\mathcal{G}}= f\Phi^* \|g\|_{\mathcal{B}(\mathcal{G})} $ where, for simplicity, we   wrote
 $\|g\|_{\mathcal{B}(\mathcal{G})}$ for $\|g\|_{\mathcal{B}(\mathcal{G})}\id_\mathcal{G}$.
 Since   the singular values of $A_1$ are the square roots of the eigenvalues of $A_1^* A_1$, which has the same non-zero-eigenvalues as $ A_1 A_1^*$ we see that the nonzero singular values of $A_1$ obey the bound
 $s_n(A_1) \le s_n(A_2)$.

 Similarly, $|f(x)|\coloneqq \sqrt{f(x)^*f(x)}$  is a non negative operator on
 $\mathcal{H}$ and
 \begin{align*}
 	A_2^* A_2 = \|g\|_{\mathcal{B}(\mathcal{G})} \Phi^* f^*f\Phi^* \|g\|_{\mathcal{B}(\mathcal{G})}
 	= \|g\|_{\mathcal{B}(\mathcal{G})} \Phi^* |f|^2\Phi^* \|g\|_{\mathcal{B}(\mathcal{G})}
 	= A_3^* A_3
 \end{align*}
 with $ A_3= |f|\Phi^* \|g\|_{\mathcal{B}(\mathcal{G})}$. So the singular
 values of $A_2$ are the same as the singular values of $A_3$ and
 without loss of generality, we can assume that $g$ is a non-negative
 function and $f$ takes values in the non-negative operators on $\mathcal{H}$.
 By scaling, we can also assume that $\|f\|_{L^p(X,{\mathcal{S}}_p(\mathcal{H}))}= \|g\|_{L^p_\text{w}(Y)}^p=1$.

 Since $\Phi: L^1(X,\mathcal{H})\to L^\infty(Y,\mathcal{G})$ is bounded, Lemma \ref{lem:kernel-op-valued} shows that it has a kernel $\Phi(\cdot,\cdot)$ such that for all $f\in L^2(X,\mathcal{H})$,
 \begin{align*}
 \Phi f(y) = \int_X \Phi(y,x) f(x) \, dx
 \end{align*}
 for almost all $y\in Y$. Moreover, $\sup_{(y,x)\in Y\times X}\|\Phi(y,x)\|_{\mathcal{B}(\mathcal{H},\mathcal{G})}= \|\Phi\|_{L^1\to L^\infty}$.
 Having reduced the estimate to scalar non-negative functions $g$ and non-negative operator-valued functions $f$ we  can rewrite  $A_{f,g}= g\Phi f$ as
 \begin{align}
 	A_{f,g}\varphi(y) = \int_{X} g(y) \Phi(y,x) f(x) \varphi(x)\, dx
 		= \int_{X}  \Phi(y,x) g(y) f(x) \varphi(x)\, dx
 \end{align}
 using that $g(y)$ is now a non-negative scalar. Thus, we can take again an arbitrary function $m:{\mathbb{R}}_+\to {\mathbb{R}} $ with $m(0)=0$ and split
  \begin{align}
  	 B_{f,g,m}\varphi(y)
  	 	& \coloneqq
  	 		\int_{X}  \Phi(y,x) m\big( g(y) f(x) \big) \varphi(x)\, dx
  	 			\label{eq:def-B-fully-operator-valued} \\
  	 H_{f,g,m}\varphi(y)
  	 	& \coloneqq
  	 		\int_{X}  \Phi(y,x) \big[g(y)f(x)-m\big( g(y) f(x) \big)\big] \varphi(x)\, dx \label{eq:def-H-fully-operator-valued}
  \end{align}
  The above expression are well-defined by the spectral theorem, since $g$ is a non-negative function and $f$ takes values in the non-negative operators on $\mathcal{H}$, so $m(g(y)f(x))$ is a bounded operator on $\mathcal{H}$ for almost all $y$ and $x$, when $m$ is bounded. Thus the integrals in \eqref{eq:def-B-fully-operator-valued} and \eqref{eq:def-B-fully-operator-valued}  converge for all $\varphi$ from a dense subset of $L^2(X,\mathcal{H})$, for example the piecewise constant functions.

  Scaling in $f$ by $\kappa>0$, we get from Ky Fan's inequality
  \begin{equation}
    \begin{split}
    	s_n(g\Phi f)
    		&= \kappa^{-1}s_n(A_{\kappa,f,g})\le \kappa^{-1}\left[ \|B_{\kappa f,g,m}\| + s_n(H_{\kappa f,g,m}) \right] \\
    		& \le \kappa^{-1} \left[ \mu + n^{-1/2} \|H_{\kappa, f,g,m}\|_{HS} \right]
    \end{split}
  \end{equation}
      where we take $\mu = \|m_1\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})} \|m_2\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})}$, the upper bound on the norm of $B_{\kappa f,g,m}$ from Lemma \ref{lem:great-fully-operator-valued} below and we used
      $s_n(H)\le n^{-1}\sum_{j=1}^n s_j(H)^2\le n^{-1}\|H\|_{HS}^2$, for any Hilbert-Schmidt operator, due to the monotonicity of its singular values.
  Thus using the bound \eqref{eq:bound-H-fully-operator-valued} one gets
  \begin{align*}
  	    	s_n(g\Phi f)
    			 \le \kappa^{-1} \left[ \mu + n^{-1/2} D^{1/2} \kappa^{p/2} \right]
  \end{align*}
 with $D= p \, \|\Phi\|_{L^1\to L^\infty}^2 \int_0^\infty (1-t^{-1}m(t))^2 t^{1-p}\, dt$, and minimizing this over $\kappa>0$ we have
 \begin{align*}
 	s_n(g\Phi f)
 		\le \frac{p}{p-2} \left( \frac{p-2}{2} \right)^{2/p}
 			(\mu^{p-2}D)^{1/p}\, n^{-1/p}
 \end{align*}
 for the singular values for all $n\in{\mathbb{N}}$.

 Making again the simplest choice $m_1(s)= 2s\id_{\{0\le s\le 1\}}$ and $m_2(s)= s^{-1}\id_{\{s\ge 1\}}$ one checks that $m=m_1*m_2= \min(t,t^{-1})$ is allowed since $m'\le 1$. Calculating the numbers finishes the proof of Theorem \ref{thm:Cwikel}.
\end{proof}

\begin{lemma}\label{lem:great-fully-operator-valued}
  Let $p>2$, $\mathcal{H}$ and $\mathcal{G}$ auxiliary Hilbert--spaces, $(X,dx)$ and
  $(Y,dy)$ $\sigma$--finite measure spaces, $0\le g\in L^p_{\text{w}}(Y)$,
  $0\le f\in L^p(X{\mathcal{S}}_p(\mathcal{H}))$, $\Phi: L^2(X,\mathcal{H})\to L^2(Y,\mathcal{G})$ unitary
  and also bounded from  $L^1(X,\mathcal{H})\to L^\infty(Y,\mathcal{G})$. Then
  for all continuous and piecewise differentiable bounded functions $m:{\mathbb{R}}_+\to {\mathbb{R}} $ with  $m(0)=0$ and $m'\le 1$ the
  operator $\widetilde  {H}_{f,g,m}$ defined in \eqref{eq:def-H-fully-operator-valued} is
  a Hilbert--Schmidt operator and
  \begin{equation}
  \begin{split}\label{eq:bound-H-fully-operator-valued}
  	\|\widetilde  {H}_{f,g,m}&\|_{{\mathcal{S}}_2(L^2(X,\mathcal{H})\to L^2(Y,\mathcal{G}))}^2
  		= \mathop{\mathrm{tr}} \nolimits_{L^2(X,\mathcal{H})}\left[ \widetilde  {H}_{f,g,m}^*\widetilde  {H}_{f,g,m} \right] \\
  	&\le p \, \|\Phi\|_{L^1\to L^\infty}^2 \int_0^\infty (1-t^{-1}m(t))^2 t^{1-p}\, dt \,
  				\|g\|_{L^p_\text{w}(Y)}^p \| f\|_{L^p(X,S_p(\mathcal{H}))}^p\,.
  \end{split}
  \end{equation}

  Moreover, if $m=m_1*m_2$, then the operator
  $\widetilde  {B}_{f,g,m}$ defined in \eqref{eq:def-B-fully-operator-valued} is
  bounded from $L^2(X,\mathcal{H})$ to $ L^2(Y,\mathcal{G})$ and
  \begin{align}\label{eq:bound-B-fully-operator-valued}
  	\|\widetilde  {B}_{f,g,m}\|_{L^2\to L^2} \le \|m_1\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})} \|m_2\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})}\, .
  \end{align}
\end{lemma}
\begin{remark}
  As the proof of Lemma \ref{lem:great-fully-operator-valued} shows one even has a bound on  $\widetilde  {B}_{f,g,m}$ in the form
  \begin{align*}
  	\sup_{f\ge 0} \big\| \sup_{g\ge 0}\|\widetilde  {B}_{f,g,m}\varphi\|_\mathcal{G} \big\|_{L^2(Y)}
  		\le \|m_1\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})} \|m_2\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})}^2
 		 		\|\varphi\|_{L^2(X,\mathcal{H})}
  \end{align*}
 	where the first supremum is taken over all functions $g:Y\to [0,\infty)$ and the second supremum is taken over all non-negative operator-valued functions $f:X\to \mathcal{B}(\mathcal{H})$.
\end{remark}

\begin{proof}
  For notational simplicity we set
  \begin{align*}
  	C= \|\Phi\|_{L^1(X,\mathcal{H})\to L^\infty(Y,\mathcal{G})} = \mathop{\mathrm{esssup}}\limits_{(y,x)\in Y\times X} \|\Phi(x,y)\|_{\mathcal{B}(\mathcal{H},\mathcal{G})}\, .
  \end{align*}
 and note
 \begin{align*}
 	\|\widetilde  {H}_{f,g,m}\|_{{\mathcal{S}}_2(L^2(X,\mathcal{H})\to L^2(Y,\mathcal{G}))}^2
 		&= \iint_{Y\times X} \mathop{\mathrm{tr}} \nolimits_{\mathcal{H}}\left[ \widetilde  {H}_{f,g,m}(y,x)^*\widetilde  {H}_{f,g,m}(y,x) \right]\, dy dx \, .
 \end{align*}
 Because $g$ is real-valued, even positive, and $f$ takes values in the non-negative, hence self-adjoint, operators
 \begin{align*}
 	\widetilde  {H}_{f,g,m}&(y,x)^*\widetilde  {H}_{f,g,m}(y,x) =\\
 		&= \big[g(y)f(x)-m\big( g(y) f(x) \big)\big]\Phi(y,x)^*\Phi(y,x)\big[g(y)f(x)-m\big( g(y) f(x) \big)\big] \\
 		&\le C^2 \big[g(y)f(x)-m\big( g(y) f(x) \big)\big]^2,
 \end{align*}
 so, setting $
 	G(u)\coloneqq \int_{Y}\left[u g(y) -m(ug(y)) \right]^2\, dy$,
 we have
 \begin{align*}
 	\|\widetilde  {H}_{f,g,m}\|_{{\mathcal{S}}_2(L^2(X,\mathcal{H})\to L^2(Y,\mathcal{G}))}^2
 		&\le  C^2 \int_{X} \mathop{\mathrm{tr}} \nolimits_\mathcal{H} G(f(x)) \, dx \, .
 \end{align*}
 With  $k(t)= (t-m(t))^2 $, the layer-cake principle shows
 \begin{align*}
 	G(u) = \int_0^\infty k'(ut)u |\{y\in Y: g(y)>t\}|\, dt\, .
 \end{align*}
 By definition $|\{y\in Y: g(y)>t\}|\le t^{-p}\|g\|_{L^p_\text{w}(Y)}^p$
 for all $t>0$. If $m(0)=0$ and $m'(t)\le 1$, then $k'\ge 0$. Thus
  \begin{align*}
 	G(u) &\le  u^p\,\|g\|_{L^p_\text{w}(Y)}^p
 				\int_0^\infty k'(t) t^{-p}\, dt\, .
 \end{align*}
 For any $0<\varepsilon<L$ we have
 \begin{align*}
 	\int_\varepsilon^L k'(ut)u t^{-p}\, dt
 		= p[k(t)t^{-1-p}]_\varepsilon^L +p \int_\varepsilon^L k(t) t^{-1-p}\, dt
 \end{align*}
 and if $0<t \mapsto k(t) t^{-1-p}\in L^1({\mathbb{R}}_+,dt)$,
 then there exist  sequences $\varepsilon_n\to 0$ and $L_n\to\infty$
 such that  $\lim_{n\to\infty}k(\varepsilon_n) \varepsilon_n^{-1-p}=0 = \lim_{n\to\infty}k(L_n) L_n^{-1-p}$,
 that is, the boundary term $ [k(t)t^{-1-p}]_{\varepsilon_n}^{L_n}$ vanishes in the limit $n\to\infty$. Hence integration by parts is
 justified as soon as the right hand side of \eqref{eq:bound-H-fully-operator-valued} is finite and
 \begin{align*}
 	\mathop{\mathrm{tr}} \nolimits_\mathcal{H} G(f(x)) \le
 		p \int_0^\infty k(t) t^{-1-p}\, dt\, \|g\|_{L^p_\text{w}(Y)}^p \mathop{\mathrm{tr}} \nolimits_\mathcal{H} (f(x)^p)
 \end{align*}
Integrating this over $X$ finishes the proof of \eqref{eq:bound-H-fully-operator-valued}.
 \smallskip

 To prove \eqref{eq:bound-B-fully-operator-valued} we introduce
 \begin{align}\label{eq:def-maximal-operator-fully-operator-valued}
 	\widetilde  {B}_{f,t,m}\varphi(y)
 		\coloneqq \int_{X}  \Phi(y,x) m\big( t f(x) \big) \varphi(x)\, dx
 		= \Phi[m(tf)\varphi](y)
 \end{align}
 for $t\ge 0$ (using $m(0)=0$). If $m=m_1*m_2$, convolution on ${\mathbb{R}}_+$, then a by now familiar calculation yields
 \begin{align*}
  	\widetilde  {B}_{f,t,m}\varphi(y)
 		=
 			\int_0^\infty \Phi[m_1(sf)\varphi](y) \, m_2(ts)\,\frac{ds}{s}
 \end{align*}
 and therefore the Cauchy--Schwarz inequality gives
 \begin{align*}
  	\|\widetilde  {B}_{f,t,m}\varphi(y)\|_\mathcal{G}
 		&\le
 			\int_0^\infty \|\Phi[m_1(sf)\varphi](y)\|_\mathcal{G} \, |m_2(ts)|\,\frac{ds}{s} \\
 		&\le \left(  \int_0^\infty \|\Phi[m_1(sf)\varphi](y)\|_\mathcal{G} ^2 \,\frac{ds}{s} \right)^{1/2}
 			\left(  \int_0^\infty  |m_2(ts)|^2\,\frac{ds}{s} \right)^{1/2} \, .
 \end{align*}
 By scaling, the right hand side above does not depend on $t>0$ anymore and, since $\widetilde  {B}_{f,0,m}\varphi(y)=0$, we get the bound
 \begin{align*}
  	\widetilde  {B}_{f,m}^*\varphi(y)
  		= \sup_{t> 0}\|\widetilde  {B}_{f,t,m}\varphi(y)\|_\mathcal{G}
 		\le \|m_2\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})} \left(  \int_0^\infty \|\Phi[m_1(sf)\varphi](y)\|_\mathcal{G} ^2 \,\frac{ds}{s} \right)^{1/2} .
 \end{align*}
 for the associated maximal operator
 $\widetilde  {B}_{f,m}^*\varphi(y)\coloneqq \sup_{t\ge 0}\|\widetilde  {B}_{f,t,m}\varphi(y)\|_\mathcal{G} $.
 In particular,
 \begin{align}\label{eq:max-op-fully-0p-valued-1}
 	\|\widetilde  {B}_{f,m}^*\varphi\|_{L^2(Y,dy)}^2
 		& \le \|m_2\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})}^2 \,\
 				\int_{Y}   \int_0^\infty \|\Phi[m_1(sf)\varphi](y)\|_\mathcal{G} ^2 \,\frac{ds}{s} \, dy\, .
 \end{align}
 Interchanging the integrals, the last factor on the right hand side of
 \eqref{eq:max-op-fully-0p-valued-1} is given by
 \begin{align*}
 	\int_0^\infty &\int_{Y}  \|\Phi[m_1(sf)\varphi](y)\|_\mathcal{G} ^2 \, \, dy\, \frac{ds}{s}
 		=
 			\int_0^\infty \|\Phi[m_1(sf)\varphi]\|_{L^2(Y,\mathcal{G})}^2 \, \frac{ds}{s} \\
 		&=  \int_0^\infty \|m_1(sf)\varphi\|_{L^2(X,\mathcal{H})}^2 \, \frac{ds}{s}
 			= \int_{X} \int_0^\infty \big\langle
 					m_1(sf(x))\varphi(x), m_1(sf(x))\varphi(x)
 				\big\rangle_\mathcal{H}  \, \frac{ds}{s} \, dx \\
 		&= \int_{X} \big\langle
 				\varphi(x), \int_0^\infty  m_1(sf(x))^2\, \frac{ds}{s}
 					\varphi(x)
 			  \big\rangle_\mathcal{H}   \, dx \, .
 \end{align*}
 As functions of the real variable $r\ge 0$ the scaling invariance of the measure $ds/s$ on ${\mathbb{R}}_+ $ and $m_1(0)=0$ gives
  $ \int_0^\infty  m_1(sr)^2\, \frac{ds}{s} = \|m_1\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})}^2\id_{\{r>0\}}$,  so the spectral theorem implies
 \begin{align*}
 	\big\langle
 		\varphi(x), \int_0^\infty  m_1(sf(x))^2\, \frac{ds}{s}
 					\varphi(x)
 	\big\rangle_\mathcal{H}
 		&=
 		\|m_1\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})}^2
 		 \big\langle
 			\varphi(x), \id_{\{f(x)>0\}} \varphi(x)
 		 \big\rangle_\mathcal{H}  \\
 		&\le \|m_1\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})}^2
 		 		\|\varphi(x)\|_\mathcal{H}^2 \, .
 \end{align*}
 Using this in \eqref{eq:max-op-fully-0p-valued-1} shows
 \begin{align}\label{eq:max-op-fully-0p-valued-2}
	\|\widetilde  {B}_{f,m}^*\varphi\|_{L^2(Y,dy)}
		\le \|m_1\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})} \|m_2\|_{L^2({\mathbb{R}}_+,\frac{ds}{s})}^2
 		 		\|\varphi\|_{L^2(X,\mathcal{H})} \,.
 \end{align}
 which proves \eqref{eq:bound-B-fully-operator-valued}, since $\|\widetilde  {B}_{f,g,m}\varphi(y)\|_\mathcal{G}\le \widetilde  {B}_{f,m}^*\varphi(y)$ for all $y\in Y$.
\end{proof}

\bigskip

Mathematical Physics (math-ph); Analysis of PDEs (math.AP); Spectral Theory (math.SP)

Cwikel's bound reloaded