---
title: Generalized Rellich's Lemmas
url: https://www.emergentmind.com/topics/generalized-rellich-s-lemmas
type: topic
---

# Generalized Rellich's Lemmas

Searching arXiv for the cited works and closely related papers on generalized Rellich lemmas and Rellich-type results.
“Generalized Rellich’s lemmas” denotes a family of extensions of classical results associated with Franz Rellich rather than a single theorem. In current usage, the term covers at least four major lineages: higher-order Hardy–Rellich inequalities and their weighted refinements; uniqueness statements for Helmholtz and scattering problems on unbounded domains; compactness theorems of Rellich–Kondrachov type in generalized Sobolev frameworks; and smooth selection results for kernels of parametrized singular linear systems. Recent arXiv literature develops these themes in continuous, discrete, degenerate, Riemannian, many-body, and complex-frequency settings [1710.06955] [1401.4531] [1602.07493] [2507.04242].

## 1. Classical lineages and the scope of the term

A first classical lineage is inequality-theoretic. Birman’s sequence \(\{I_n\}_{n\in\mathbb N}\) extends the Hardy inequality (\(n=1\)) and the Rellich inequality (\(n=2\)) to all orders, with
\[
\int_0^\infty |f^{(n)}(x)|^2\,dx
\;\ge\;
\frac{[(2n-1)!!]^2}{2^{2n}}
\int_0^\infty \frac{|f(x)|^2}{x^{2n}}\,dx.
\]
For \(n=1\) this gives the classical Hardy constant \(1/4\), and for \(n=2\) it gives the Rellich constant \(9/16\) on the half-line [1710.06955].

A second lineage is uniqueness for Helmholtz-type equations. In this setting, a Rellich lemma states that an outgoing solution whose far-field or large-radius boundary behavior is sufficiently small must vanish identically. This perspective is central both in half-space and cone uniqueness theorems for \((\Delta+\lambda)u=f\) on unbounded domains and in recent extensions to complex wavenumbers associated with scattering poles [1401.4531] [2507.04242].

A third lineage is compactness. The classical Rellich–Kondrachov theorem appears as a special case of abstract compact embeddings for generalized Sobolev spaces defined by nonnegative quadratic forms \(Q(x,\xi)\), including degenerate and weighted cases [1110.6907].

A fourth lineage is algebraic and analytic rather than PDE-theoretic. Rellich’s 1969 lemma concerns analytic families of singular matrices \(D(\varepsilon)\) with \(\det D(\varepsilon)=0\), asserting the existence of a locally analytic unit vector \(x(\varepsilon)\) satisfying \(D(\varepsilon)x(\varepsilon)=0\). This has been generalized to \(C^\ell\) and multiparameter settings, to inhomogeneous systems, and to orthonormal frame fields of kernel vectors [2301.13164].

This multiplicity of meanings is a recurrent source of confusion. The modern literature does not use “generalized Rellich’s lemma” for a single canonical extension; instead, it refers to a class of results sharing a common role: converting weak decay, weak regularity, or structural singularity into rigidity, compactness, or sharp coercive control.

## 2. Higher-order Hardy–Rellich inequalities and weighted identities

In the higher-order inequality setting, Gesztesy, Littlejohn, Michael, and Wellman place Birman’s inequalities in a natural Hilbert-space framework. On the half-line,
\[
H_n([0,\infty))
=
\{f:\ f^{(j)}\in AC_{\mathrm{loc}},\ j=0,\dots,n-1,\ f^{(n)}\in L^2(0,\infty),\ f^{(j)}(0)=0\},
\]
with inner product \((f,g)_{H_n}=\int_0^\infty f^{(n)}\overline{g^{(n)}}\), one has the key inclusion \(f\in H_n([0,\infty))\Rightarrow f'\in H_{n-1}([0,\infty))\). Repeated application of the classical Hardy inequality to \(f,f',f'',\dots\) yields the entire Birman sequence, and the constants \(\frac{[(2n-1)!!]^2}{2^{2n}}\) are sharp. Equality occurs only for the trivial function. The same inequalities hold on finite intervals in the standard Sobolev space \(H_0^n((0,b))\), and the constants are again sharp [1710.06955].

The same paper gives a spectral-operator interpretation through the continuous \(n\)-fold Cesàro operator
\[
(T_n f)(x)
=
\frac1{x^n}
\underbrace{\int_0^x\int_0^{t_1}\cdots\int_0^{t_{n-1}}}_{n\ \mathrm{integrals}}
f(u)\,du\,dt_{n-1}\cdots dt_1,
\]
whose norm satisfies \(\|T_n\|_{L^2\to L^2}=2^n(2n-1)!!\). Moreover, \(T_n^{-1}=\frac{d^n}{dx^n}[x^n\cdot]\) on a suitable dense domain, and Mellin-transform methods show that \(T_n\) is normal with purely absolutely-continuous spectrum [1710.06955].

A different but complementary generalization is the weighted \(L^p\)-Rellich identity for quasilinear second-order degenerate elliptic operators. For vector fields \(X_i=\sum_j \sigma_{ij}(x)\partial_{x_j}\), \(\nabla_{\mathcal L}=\sigma\nabla\), \(\mathcal L=-\nabla_{\mathcal L}^*\nabla_{\mathcal L}\), and \(\mathcal L_pu=-\nabla_{\mathcal L}^*(|\nabla_{\mathcal L}u|^{p-2}\nabla_{\mathcal L}u)\), the paper “Sharp remainder formulae for general weighted Hardy and Rellich type inequalities for \(1<p<\infty\)” proves an exact identity of the form
\[
\int_\Omega V\,|\mathcal L u|^p
=
\lambda\int_\Omega W\,|u|^p
+\text{nonnegative remainder terms},
\]
under the structural condition
\[
\mathcal L\!\bigl(V\,|\mathcal L\phi|^{p-2}\mathcal L\phi\bigr)
=
\lambda\,W\,|\phi|^{p-2}\phi
\quad\text{in }\mathcal D'(\Omega),
\]
with \(-\frac{\mathcal L\phi}{\phi}\ge0\). The remainder is organized through the nonnegative functional
\[
C_p(\xi,\eta)=|\xi|^p-|\xi-\eta|^p-p|\xi-\eta|^{p-2}\Re[(\xi-\eta)\cdot\overline\eta]\ge0,
\]
which allows the argument to work for all \(1<p<\infty\), not only \(p\ge2\). Discarding the nonnegative terms yields a sharp weighted Rellich inequality, and equality in the vanishing-remainder sense occurs only in the virtual case \(u=c\phi\) [2603.02381].

In the Euclidean Laplacian case, taking \(\sigma=I_N\), \(V\equiv1\), and \(\phi(x)=|x|^{-\alpha}\) with \(\alpha=\frac{N-2p}{p}\), one obtains explicit sharp identities. For \(p=2\) and \(N\ge4\), the paper derives a new \(L^2\)-Rellich identity with coefficient \(C_N=\frac{N(N-4)}4\), and dropping the nonnegative remainders recovers the classical sharp Rellich inequality [2603.02381].

A Riemannian analogue appears in work on complete noncompact manifolds. For a positive weight \(p:M\to(0,\infty)\) satisfying \(|\nabla p|\ge c>0\) almost everywhere and suitable lower bounds on \(\Delta p\) or on \(\mathrm{div}(a(x)|\nabla p|^{p-2}\nabla p)\), the generalized Rellich inequality on \((M,g)\) takes the form
\[
\int_M \rho^{\,a-4+2p}|\Delta\varphi|^p\,dV
\ge
K_p(a)\int_M \rho^{\,a-4}|\varphi|^p\,dV
+
L_p(a)\int_M \rho^{\,a+p-2}|\nabla\varphi|^p\,dV
\]
for \(\varphi\in C_0^\infty(M\setminus\{x_0\})\), where \(\rho=d(x,x_0)\), \(|\nabla\rho|\equiv1\), and
\[
K_p(a)=p^p\Bigl|\frac{a+4-n}{p}\Bigr|^p,
\qquad
L_p(a)=p^p\Bigl|\frac{n-a-2p}{p}\Bigr|^{p-1}(a+2-n).
\]
When \(M=\mathbb R^n\) and \(\rho(x)=|x|\), this recovers the Euclidean weighted Rellich inequality, and the constants agree with the Euclidean literature in the flat case [2411.07260].

## 3. Discrete Rellich inequalities and factorization methods

The discrete theory replaces derivatives by lattice differences and the Laplacian by the discrete Dirichlet Laplacian. For finitely supported sequences \(A=\{A_n\}_{n\in\mathbb N_0}\) with \(A_0=A_1=0\), the one-dimensional operator is
\[
(-\Delta A)_n=2A_n-A_{n+1}-A_{n-1},
\]
and the bi-Laplacian is \((-\Delta)^2\). The classical one-dimensional discrete Rellich inequality states that, for sequences supported in \(n\ge2\),
\[
\sum_{n=2}^\infty |(-\Delta)^2A|_n^2
\ge
\frac1{16}\sum_{n=2}^\infty n^{-4}|A_n|^2,
\]
with sharp constant \(1/16\) [2309.04923].

The paper “On improvements of the Hardy, Copson and Rellich inequalities” strengthens this by a factorization identity. For \(A\in C_c(\mathbb N_0)\) with \(A_0=A_1=0\), there exists a nonnegative remainder operator \(R\) and a weight sequence \(p^{(2)}=\{p_n^{(2)}\}_{n\ge2}\) such that
\[
\sum_{n=2}^\infty |(-\Delta)^2A|_n^2
=
\sum_{n=2}^\infty p_n^{(2)}|A_n|^2
+
\sum_{n=0}^\infty |(RA)_n|^2.
\]
Hence
\[
\sum_{n=2}^\infty |(-\Delta)^2A|_n^2
\ge
\sum_{n=2}^\infty p_n^{(2)}|A_n|^2,
\qquad
p_n^{(2)}=\frac{((-\Delta)^2u)_n}{u_n},
\quad
u_n=n^{3/2},
\]
and an explicit recursion shows \(p_n^{(2)}>\frac1{16n^4}\) for all \(n\ge2\) [2309.04923].

The proof proceeds by seeking a decomposition
\[
(-\Delta)^2-\mathrm{diag}(p_n^{(2)})=RR^T
\]
and writing
\[
(RA)_n=\lambda_nA_n-\mu_nA_{n+1}+\nu_nA_{n+2}.
\]
The coefficients are chosen so that \(R\) annihilates the trial sequence \(u_n=n^{3/2}\), which yields a three-term recursion for the \(\lambda_n\). The paper then proves positivity and quantitative bounds,
\[
n^{-1}p_n^{(2)}<\lambda_n<p_n^{(2)}p_{n+1}^{(2)},
\qquad n\ge2,
\]
establishing the existence of the factorization [2309.04923].

The optimality issue is subtle. Gerhat–Krejčiřík–Štampach had shown that \(p_n^{(2)}\) cannot be enlarged pointwise for all but finitely many \(n\) without violating the inequality, and the present paper emphasizes that the same bounds reappear in infinitely many factorization presentations. In this sense, \(p^{(2)}\) is presented as the best possible pointwise Rellich weight in the discrete half-line setting. The same discrete Rellich improvement feeds directly into sharpened Knopp inequalities of orders \(a=2\) and \(a\ge1\), as well as improvements of generalized Hardy-type inequalities [2309.04923].

## 4. Uniqueness on unbounded domains and scattering at complex frequencies

In PDE and scattering theory, generalized Rellich lemmas are uniqueness statements for solutions of Helmholtz-type equations under decay and support conditions. Vesalainen studies
\[
(\Delta+\lambda)u=f
\quad\text{in }\mathbb R^n,
\qquad
u\in B_2,
\]
with \(B^*\) and \(B_2\) the Agmon–Hörmander spaces, and proves several generalizations of the classical uniqueness theorem on unbounded domains. In the half-space case, if \(f\) is super-exponentially decaying and supported in \(\{x_n\ge0\}\), then the scattered wave \(u\) vanishes in the lower half-space \(H=\mathbb R^{n-1}\times(-\infty,0)\). The proof uses a Carleman estimate extracted from Sylvester–Uhlmann,
\[
\|e^{-Tx_n}u\|_{L^2}
\le
C\,\frac1T\|e^{-Tx_n}(\Delta+\lambda)u\|_{L^2},
\]
together with cutoff arguments and unique continuation [1401.4531].

The same paper proves a discrete-lattice analogue on \(\mathbb Z^n\). For
\[
(-\Delta_{\mathrm{disc}}-\lambda)u=f,
\qquad \lambda\in(0,n),
\]
if \(f\) is super-exponentially decaying and vanishes in the discrete cone
\[
C=\{m\in\mathbb Z^n:\ |m_1|+\cdots+|m_{n-1}|<m_n\},
\]
then \(u(m)=0\) for all \(m\in C\). The key inputs are a discrete Paley–Wiener theorem and repeated use of the lattice equation to propagate the value \(u(0)\) arbitrarily deep into the cone, where super-exponential decay forces it to vanish [1401.4531].

Vesalainen also gives complex-variable proofs under weaker hypotheses: exponential decay of \(f\) in a half-space and a polynomial-decay case with exponentially thin support. In Fourier variables, one uses
\[
(|\xi|^2-\lambda)\widehat u(\xi)=\widehat f(\xi),
\]
analytic continuation of \(\widehat f\), vanishing on the real sphere \(\{|\xi|^2=\lambda\}\) via the Rellich–Vekua lemma, and division by \(|\xi|^2-\lambda\) to continue \(\widehat u\). These results feed into the discreteness of non-scattering energies for non-compactly supported potentials with suitable decay and support assumptions [1401.4531].

A further extension concerns complex wavenumbers \(k\in\mathbb C_-\), where outgoing fields may grow exponentially at infinity. Liu, Sun, and Zhang study exterior scattering in \(\Omega=\mathbb R^2\setminus\overline D\) for a bounded Lipschitz obstacle \(D\), defining outgoing solutions by the Green representation with \(\Phi_k(x,y)=\frac{i}{4}H_0^{(1)}(k|x-y|)\). They prove two generalized Rellich lemmas. The first assumes
\[
\lim_{r\to\infty}
e^{2|\Im(kr)|}
\int_{|x|=r}|u(x)|^2\,ds(x)=0,
\]
for any \(C^2\) solution of \(\Delta u+k^2u=0\), and concludes \(u\equiv0\). The second assumes that \(u\) is outgoing and only requires
\[
\lim_{r\to\infty}
e^{2\Im(kr)}
\int_{|x|=r}|u(x)|^2\,ds(x)=0.
\]
When \(\Im k=0\) and \(k>0\), these reduce to the classical Rellich lemma [2507.04242].

The proofs use Fourier–Bessel expansions,
\[
u(r,\phi)=\sum_{n=-\infty}^{\infty} a_n(r)Y_n(\phi),
\qquad
a_n(r)=\alpha_nH_n^{(1)}(kr)+\beta_nH_n^{(2)}(kr),
\]
together with large-argument asymptotics of the Hankel functions. The new feature is the exponential reweighting, which compensates exactly for the growth of \(H_n^{(1)}(kr)\) and \(H_n^{(2)}(kr)\) when \(\Im k<0\). This addresses a common misconception: the classical unweighted Rellich condition is not adequate for complex poles with negative imaginary part, because outgoing solutions no longer decay in the usual sense [2507.04242].

These lemmas are then used to prove uniqueness of obstacles from far-field data at non-real \(k\), and they underpin an inside-out duality for scattering poles. In the latter, the near-field operator \(\mathcal N_k\) factorizes as
\[
\mathcal N_k=G\,\mathcal S_k,
\]
where \(\mathcal S_k\) is the single-layer operator and \(G\) is an interior solution operator. Off the set of scattering poles, \(\mathcal S_k\) is injective with dense range in the appropriate trace space, so \(\mathcal N_k\) is injective; at scattering poles, injectivity fails. The linear sampling method then identifies poles through the blow-up of regularized solution norms for the equation \(\mathcal N_k g=\Phi_k(\cdot,z)|_\Gamma\) [2507.04242].

## 5. Compactness and spectral Rellich theorems

The generalized compactness theory begins with abstract Sobolev spaces built from possibly degenerate quadratic forms. Let \(Q(x,\xi)=\xi^TQ(x)\xi\) be a measurable family of symmetric nonnegative quadratic forms on a finite measure space \((X,\mu)\), and define
\[
W^{Q,p}(X,\mu)
=
\overline{\Lip_{Q,p}(X)}^{\|\cdot\|_{W^{Q,p}}},
\qquad
\|f\|_{W^{Q,p}}
=
\Bigl(\|f\|_{L^p}^p+\int_X Q(x,\nabla f(x))^{p/2}\,d\mu\Bigr)^{1/p}.
\]
Chua, Rodney, and Wheeden prove an abstract compact embedding theorem: if a bounded set \(S\subset L^N(X,\mu)\times\mathcal X\) satisfies localization, finite-overlap control, and a Poincaré-type estimate on finitely many small sets covering all but an \(\varepsilon\)-portion of \(X\), then the projection \(\pi(f,g)=f\) is compact from \(S\) into \(L^q(X,\mu)\) for every \(1\le q<N\). This gives the classical Rellich–Kondrachov theorem when \(Q(x,\xi)=|\xi|^2\), but it also yields compactness for weighted boundary-degenerate spaces, \(s\)-John domains with distance weights, and even quasimetric spaces without any gradient structure, provided local oscillation can be controlled in the required way [1110.6907].

A spectral version appears in the many-body Schrödinger setting. Ito and Skibsted consider a generalized \(N\)-body Hamiltonian
\[
H=-\Delta+\sum_{b\in B}V_b(x)
\quad\text{on }L^2(\Omega),
\]
where \(\Omega\) is a hard-core configuration space and the pair potentials may be singular, with a decomposition \(V_b=V_b^{(1)}+V_b^{(2)}\) satisfying the soft-potential hypotheses of Condition 1.2. They define the optimal Besov-type space
\[
B^0
=
\{\psi\in L^2_{\mathrm{loc}}(\Omega):\ \chi_{|x|<R}\psi\in L^2\ \forall R,\ \lim_{R\to\infty}\|\chi_{R<|x|<2R}\psi\|=0\},
\]
and prove the generalized Rellich theorem: if \(\psi\in B^0\) and \(H\psi=E\psi\) with \(E\) a real non-threshold spectral point, then \(\psi\in L^2(\Omega)\) [1602.07493].

The proof uses a Mourre estimate with a rescaled Graf vector field,
\[
A_R=\tfrac12(w_R\cdot p+p\cdot w_R),
\]
yielding positivity of \(i[H,A_R]\) up to compact errors on a non-threshold spectral window, together with functional-calculus localization for a propagation observable \(B\). The argument proceeds in two stages: first proving \(\psi\in L^2_{1/4}\), then bootstrapping to \(\psi\in L^2\). According to the paper, the spaces \(B^0\) and \(B^*\) are optimal, and the theorem covers Coulomb atomic and molecular models with hard-core nuclei [1602.07493].

These compactness and spectral results illustrate a broader structural point. In one direction, generalized Rellich–Kondrachov theorems convert local oscillation estimates into precompactness. In another, generalized Rellich theorems in scattering and many-body analysis convert weak generalized-eigenfunction control into genuine square integrability. The terminology is shared because both types of result express a rigidity principle beyond the classical Euclidean Sobolev setting.

## 6. Parametrized singular systems, frame fields, and sensitivity

Rellich’s name also appears in the smooth dependence of solutions to parametrized singular linear systems. Let
\[
D(\varepsilon)=[Y_{ij}(\varepsilon)]_{1\le i,j\le n}
\]
be an \(n\times n\) real matrix with analytic entries near \(\varepsilon_0\in\mathbb R\), and assume \(\det D(\varepsilon)=0\) for all \(\varepsilon\) in a neighborhood. Rellich’s original lemma asserts the existence, after possibly shrinking the neighborhood, of analytic functions \(x_i(\varepsilon)\) such that the vector \(x(\varepsilon)\) satisfies
\[
D(\varepsilon)x(\varepsilon)=0,
\qquad
|x(\varepsilon)|=1.
\]
The proof uses nonvanishing minors, cofactor vectors, a vanishing-order argument, and normalization [2301.13164].

Bellido and Prieto-Martínez extend this in three directions. First, if the entries \(Y_{ij}\) are only \(C^\ell\), one still obtains a \(C^\ell\) null-vector field \(f(\varepsilon)\), and if \(\operatorname{rank}D(\varepsilon_0)\) equals the maximal rank in the neighborhood, then \(f\) can be normalized to a nonvanishing \(C^\ell\) unit kernel vector locally. Second, the same conclusion holds for several real parameters \(\varepsilon\in\mathbb R^N\), both in the analytic and \(C^\ell\) categories, again provided one localizes at maximal-rank points. Third, an inhomogeneous version is proved: if
\[
\rank D(\varepsilon_0)=\rank[D(\varepsilon_0)\mid b(\varepsilon_0)]
\]
and \(\rank[D(\varepsilon)\mid b(\varepsilon)]\le r\) in a neighborhood, then \(D(\varepsilon)x(\varepsilon)=b(\varepsilon)\) admits a local analytic or \(C^\ell\) solution obtained from Cramer’s rule on a suitable invertible minor [2301.13164].

The same paper develops frame fields of kernel vectors. If \(k=n-r\) is the nullity of \(D(\varepsilon)\), then under the hypotheses of the analytic theorem there exists an analytic \(k\)-frame field
\[
(x_1(\varepsilon),\dots,x_k(\varepsilon))
\]
such that \(D(\varepsilon)x_i(\varepsilon)=0\) and \(x_i(\varepsilon)\cdot x_j(\varepsilon)=\delta_{ij}\). The inductive step adjoins a rank-one term,
\[
B(\varepsilon)=D(\varepsilon)+x_k(\varepsilon)x_k(\varepsilon)^T,
\]
which increases the rank by one and reduces the nullity by one, allowing iteration [2301.13164].

A further development concerns \((n-1)\)-deficient systems, where \(\dim\ker D(\varepsilon)=1\). Differentiating
\[
D(\varepsilon)x(\varepsilon)=0,
\qquad
\|x(\varepsilon)\|=1,
\]
at \(\varepsilon_0\) gives
\[
D(\varepsilon_0)\frac{\partial x}{\partial\varepsilon_j}(\varepsilon_0)
=
-\frac{\partial D}{\partial\varepsilon_j}(\varepsilon_0)u,
\qquad
u^T\frac{\partial x}{\partial\varepsilon_j}(\varepsilon_0)=0,
\]
where \(u=x(\varepsilon_0)\). This yields a direct sensitivity method. The paper also introduces an adjoint method based on the zero-Lagrangian
\[
\mathcal L(x,p,\lambda)=F(x)+p^TD(\varepsilon)x+\lambda(x^Tx-1),
\]
leading to the adjoint system
\[
D(\varepsilon_0)^Tp
=
-\nabla_xF(x(\varepsilon_0))-\lambda u,
\]
and the sensitivity formula
\[
\frac{d}{d\varepsilon}F(x(\varepsilon))\Big|_{\varepsilon=\varepsilon_0}
=
p^T\frac{\partial D}{\partial\varepsilon}(\varepsilon_0)u.
\]
The paper emphasizes that this is asymptotically more efficient when many different scalar functionals are considered for the same system matrix [2301.13164].

The examples in that work also clarify limitations. Analyticity cannot in general be replaced by mere \(C^k\) without rank constancy, and in the multiparameter case a unit null-vector may fail to extend continuously through points where maximal-rank localization is unavailable. These obstructions are part of the modern meaning of generalized Rellich lemmas in the algebraic setting: they are not only existence theorems, but also sharp descriptions of when smooth or analytic kernel selection is possible [2301.13164].

Source: https://www.emergentmind.com/topics/generalized-rellich-s-lemmas