---
title: Fractional Kirchhoff–Choquard System
url: https://www.emergentmind.com/topics/fractional-kirchhoff-choquard-system
type: topic
---

# Fractional Kirchhoff–Choquard System

Searching arXiv for recent and foundational papers on fractional Kirchhoff–Choquard systems.
A fractional Kirchhoff–Choquard system is a class of doubly nonlocal elliptic problems in which the principal operator is modulated by a Kirchhoff coefficient depending on a global energy, while the reaction term contains a Choquard or Hartree convolution with a Riesz-type kernel. In the fractional setting, the diffusion is governed by a fractional Laplacian, a fractional \(p\)-Laplacian, a weighted fractional \(p\)-Laplacian, or variable-order fractional operators; in several works the model is studied either as a single equation or as a genuinely coupled system for \((u,v)\). Across these variants, the common structure is the coexistence of nonlocal diffusion, nonlocal global coupling through the Kirchhoff term, and nonlocal long-range interaction through the Choquard term, often at Hardy–Littlewood–Sobolev critical growth [1808.07996], [2401.08310], [2509.07597].

## 1. Defining structure and canonical forms

The defining Kirchhoff feature is that the coefficient multiplying the principal operator depends on a global norm or seminorm of the solution. In the whole-space fractional \(p\)-Kirchhoff–Choquard equation, the model takes the form
\[
M\!\left(\|u\|_{W}^{p}\right)(-\Delta)^s_p u + V(x)|u|^{p-2}u = \lambda \,(I_\mu * F(u))\, f(u) \quad \text{in } \mathbb{R}^N,
\]
with
\[
\|u\|_{W}^{p} = [u]_{s,p}^{p} + \int_{\mathbb{R}^N} V(x)|u|^p\,dx,
\qquad
(I_\mu * F(u))(x)=\int_{\mathbb{R}^N} \frac{F(u(y))}{|x-y|^\mu}\,dy
\]
[1808.07996]. This formulation exhibits the standard threefold nonlocality: the fractional operator, the energy-dependent Kirchhoff coefficient, and the Choquard convolution.

A critical Heisenberg-group variant replaces Euclidean geometry by \(\mathbb H^N\) and the Euclidean Riesz kernel by the group kernel \(|\eta^{-1}\xi|^{-\lambda}\). The equation studied there is
\[
M\!\left(\|u\|_\mu^{p}\right) \Big(\mu\,(-\Delta)^s_p u + V(\xi)|u|^{p-2}u\Big)
=
f(\xi,u)
+
\left(\int_{\mathbb H^N}\frac{|u(\eta)|^{Q_\lambda^*}}{|\eta^{-1}\xi|^\lambda}\,d\eta\right)
|u(\xi)|^{Q_\lambda^*-2}u(\xi),
\]
with
\[
\|u\|_\mu^p = \mu [u]_{s,p}^p+\int_{\mathbb H^N}V(\xi)|u|^p\,d\xi
\]
[2401.08310].

A normalized coupled system on \(\mathbb R^N\) introduces two unknowns and \(L^2\)-mass constraints:
\[
\begin{cases}
\left(M(\|u\|_s^2)\right)(-\Delta)^s u = \lambda_1 u +(I_\mu*|v|^{2^*_{\mu,s}})\,|u|^{2^*_{\mu,s}-2}u +\alpha p\,(I_\mu*|v|^q)\,|u|^{p-2}u,\\[4pt]
\left(M(\|v\|_s^2)\right)(-\Delta)^s v = \lambda_2 v +(I_\mu*|u|^{2^*_{\mu,s}})\,|v|^{2^*_{\mu,s}-2}v +\alpha q\,(I_\mu*|u|^p)\,|v|^{q-2}v,\\[4pt]
\int_{\mathbb R^N}|u|^2=d_1^2,\qquad \int_{\mathbb R^N}|v|^2=d_2^2,
\end{cases}
\]
with \(M(t)=a+bt\), \(a,b>0\), and \(\lambda_1,\lambda_2\) as Lagrange multipliers [2509.07597].

These models admit several extensions. Weighted singular versions combine Choquard terms with Hardy-type and Sobolev-critical singular weights [2410.05185]. Variable-exponent and variable-order formulations replace constant \(p\), \(s\), and \(\mu\) by spatially dependent exponents and kernels [2005.00617], [2005.09221]. Exponential-growth versions replace polynomial reaction terms by Trudinger–Moser type nonlinearities [1908.11285]. This suggests that the term “fractional Kirchhoff–Choquard system” is best understood as a structural category rather than a single canonical equation.

## 2. Fractional operators, Kirchhoff coefficients, and Choquard interactions

The fractional diffusion component is represented in several inequivalent but related ways. In the constant-exponent fractional \(p\)-Kirchhoff theory, the operator is the fractional \(p\)-Laplacian
\[
(-\Delta)^s_p \varphi(x) = 2\,\text{P.V.}\int_{\mathbb{R}^N} \frac{|\varphi(x)-\varphi(y)|^{p-2}\big(\varphi(x)-\varphi(y)\big)} {|x-y|^{N+ps}}\,dy,
\]
with Gagliardo seminorm
\[
[u]_{s,p}^{p} = \int_{\mathbb{R}^{2N}} \frac{|u(x)-u(y)|^p}{|x-y|^{N+ps}}\,dx\,dy
\]
[1808.07996]. On the Heisenberg group, the analogous nonlocal seminorm uses \(|\eta^{-1}\xi|\) and the homogeneous dimension \(Q=2N+2\) [2401.08310]. In weighted singular problems, the operator becomes
\[
(-\Delta)_{p,\theta}^{s}u(x) = \text{p.v.}\int_{\mathbb{R}^N} \frac{|u(x)-u(y)|^{p-2}(u(x)-u(y))} {|x|^{\theta_1}|x-y|^{N+sp}|y|^{\theta_2}} \,dy
\]
[2410.05185]. In variable-order settings, the principal term is
\[
(-\Delta)^{s(\cdot,\cdot)}_{p(\cdot,\cdot)}u(x) = \mathrm{P.V.}\int_{\mathbb{R}^N} \frac{|u(x)-u(y)|^{p(x,y)-2}(u(x)-u(y))} {|x-y|^{N+s(x,y)p(x,y)}}\,dy
\]
[2005.09221].

The Kirchhoff coefficient is the second defining ingredient. Several prototypes recur: nondegenerate lower-bounded coefficients such as \(M(t)=m_0+bt^{\theta-1}\) or \(M(t)=a+bt\) [1808.07996], [2509.07597]; power-type degenerate coefficients such as \(M(t)=t^{\theta-1}\), which satisfy \(M(0)=0\) and create additional analytical difficulties [2106.10856]; and generalized coefficients controlled by inequalities of the form
\[
M(t)\,t \le \theta\,\mathcal M(t),\qquad \mathcal M(t)=\int_0^t M(\tau)\,d\tau
\]
or
\[
m(t)\,t\le \xi M(t)
\]
[1808.07996], [2410.05185]. In all of these cases, the diffusion strength depends on the total energy rather than pointwise data, which is the classical Kirchhoff effect in a fractional nonlocal setting.

The Choquard term introduces long-range interaction. In Euclidean models it is usually written with the Riesz potential \(I_\mu(x)=|x|^{-\mu}\), \(0<\mu<N\), so that
\[
(I_\mu*f)(x)=\int_{\mathbb{R}^N}\frac{f(y)}{|x-y|^\mu}\,dy
\]
[2509.07597]. The resulting nonlinearities include generalized Hartree terms \((I_\mu*F(u))f(u)\) [1808.07996], critical terms of the form
\[
\left(\int_\Omega \frac{|u(y)|^{2^*_{\mu,s}}}{|x-y|^\mu}\,dy\right)|u|^{2^*_{\mu,s}-2}u
\]
[2203.06471], and coupled cross-interaction terms involving both \(u\) and \(v\) [2509.07597]. The Heisenberg-group analogue replaces \(|x-y|^{-\mu}\) by \(|\eta^{-1}\xi|^{-\lambda}\) and the Euclidean critical exponent by the corresponding HLS critical exponent on \(\mathbb H^N\) [2401.08310].

## 3. Critical exponents, geometry, and noncompactness

A central organizing principle is the distinction between subcritical and critical Choquard growth. In the fractional Euclidean setting, the upper Hardy–Littlewood–Sobolev critical exponent is
\[
2^*_{\mu,s}:=\frac{2N-\mu}{N-2s},
\]
while the lower exponent is
\[
2_{\mu,*}:=\frac{2N-\mu}{N}
\]
[2509.07597]. These exponents determine whether the convolution term lies in a compact regime or at the threshold where concentration phenomena occur.

For bounded-domain problems with scalar unknown, the critical Choquard exponent is again
\[
2^*_{\mu,s}=\frac{2N-\mu}{N-2s},
\]
and the associated term is critical in the sense of the Hardy–Littlewood–Sobolev inequality [2203.06471], [2106.10856]. In the Heisenberg-group framework, the critical exponent is written
\[
Q_\lambda^*=\frac{2Q-\lambda}{Q-2s},
\qquad Q=2N+2,
\]
and marks the critical threshold for the fractional Sobolev embedding on \(\mathbb H^N\) [2401.08310].

Weighted singular models introduce a weighted critical Sobolev exponent
\[
p_s^*(\beta,\theta)=\frac{p(N-\beta)}{N-sp-\theta}
\]
and also weighted Choquard windows
\[
p_s^{\flat}(\delta,\mu)=\frac{(N-\delta-\mu/2)p}{N},
\qquad
p_s^{\sharp}(\delta,\theta,\mu)=\frac{(N-\delta-\mu/2)p}{N-sp-\theta}
\]
[2410.05185]. Variable-exponent settings replace these fixed exponents by pointwise critical bounds \(p^*(x)\) and variable HLS balance conditions [2005.00617], [2005.09221].

The analytic consequence of criticality is loss of compactness. Several papers state this explicitly: standard embeddings fail to be compact at the critical Sobolev or HLS level [2401.08310], [2410.05185], [2203.06471]. On unbounded domains, translation invariance and escape of mass add a second source of noncompactness [1808.07996], [2410.05185]. A plausible implication is that the Kirchhoff coefficient intensifies the difficulty, because weak convergence of \(u_n\) alone does not directly control the nonlinear factor \(M(\|u_n\|)\); this point is stated directly in the bounded-domain critical Choquard analysis [2203.06471].

## 4. Variational formulation and principal analytical tools

The dominant framework is variational. The energy functional usually has a Kirchhoff part minus a Choquard part, with additional local terms when present. For the whole-space generalized model,
\[
J_\lambda(u) = \frac{1}{p}\mathcal M\!\left(\|u\|_W^p\right) - \frac{\lambda}{2} \int_{\mathbb{R}^{2N}} \frac{F(u(x))F(u(y))}{|x-y|^\mu}\,dx\,dy
\]
[1808.07996]. For the Heisenberg-group critical problem,
\[
I_\mu(u) = \frac1p\,M(\|u\|_\mu^p) - \frac{1}{2Q_\lambda^*} \int_{\mathbb H^N}\int_{\mathbb H^N} \frac{|u(\xi)|^{Q_\lambda^*}|u(\eta)|^{Q_\lambda^*}}{|\eta^{-1}\xi|^\lambda} \,d\eta\,d\xi - \int_{\mathbb H^N}F(\xi,u)\,d\xi
\]
[2401.08310]. For the normalized two-component system,
\[
J(u,v) = \frac{a}{2}\big([u]_s^2+[v]_s^2\big) +\frac{b}{4}\big([u]_s^4+[v]_s^4\big) -\frac12\int_{\mathbb{R}^N}(I_\mu*|u|^{2^*_{\mu,s}})|v|^{2^*_{\mu,s}} -\alpha\int_{\mathbb{R}^N}(I_\mu*|u|^p)|v|^q
\]
on the constraint manifold
\[
S(d_1,d_2)= \left\{(u,v)\in H^s(\mathbb{R}^N)\times H^s(\mathbb{R}^N): \|u\|_2^2=d_1^2,\ \|v\|_2^2=d_2^2 \right\}
\]
[2509.07597].

The principal tools are recurrent across the literature. The Hardy–Littlewood–Sobolev inequality controls the convolution term in essentially every formulation [1808.07996], [2203.06471], [2106.10856]. When singular weights are present, a doubly weighted Stein–Weiss inequality plays the same role [2410.05185]. Fractional Sobolev embeddings provide the local compactness input in subcritical regimes [1808.07996], while concentration–compactness handles critical loss of mass and concentration at points or at infinity [2401.08310], [2203.06471], [2003.05194]. Mountain-pass geometry is the standard existence mechanism for nontrivial critical points [1808.07996], [2410.05185]. Symmetry-based multiplicity uses Krasnosel’skii genus, Benci pseudo-index, Fountain theorem, or Dual Fountain theorem, depending on the parity structure and the functional setting [2401.08310], [2005.09221], [2003.05194].

Several papers replace the Palais–Smale condition by the Cerami condition, especially when growth is slower or compactness is weak. This occurs in the generalized whole-space equation without Ambrosetti–Rabinowitz [1808.07996] and in the weighted singular critical problem on \(\mathbb R^N\) [2410.05185]. In normalized problems, the Pohozaev identity and Pohozaev manifold are central:
\[
\mathcal{P}(d_1,d_2):= \{(u,v)\in S(d_1,d_2): P_\alpha(u,v)=0\}
\]
[2509.07597]. In singular problems, truncation is required because the functional is not \(C^1\); the regularized problem
\[
M([u]_{X_0}^2)(-\Delta)^s u = \lambda \big((u^+) + \tfrac1n\big)^{-\gamma} + \left(\int_\Omega \frac{|u(y)|^{2^{*}_{\mu,s}}}{|x-y|^\mu}\,dy\right)|u|^{2^{*}_{\mu,s}-2}u
\]
is introduced precisely for that purpose [2106.10856].

## 5. Existence, multiplicity, and normalized solutions

The existence theory spans subcritical, critical, singular, weighted, and constrained regimes. In the generalized whole-space fractional \(p\)-Kirchhoff equation with Choquard nonlinearity, there is a nontrivial weak solution for every \(\lambda>0\) under assumptions \((V)\), \((M_1)\), \((M_2)\), and \((F_1)\)–\((F_4)\), without assuming the Ambrosetti–Rabinowitz condition [1808.07996]. A similar all-\(\lambda\) existence theorem holds for the weighted critical singular equation in \(\mathbb R^N\): if \(0<\mu<ps+\theta<N\) and \((V_1)\), \((m_1)\)-\((m_2)\), \((F_1)\)-\((F_4)\) hold, then the problem admits at least one nontrivial weak solution for every \(\lambda>0\) [2410.05185].

Critical bounded-domain Choquard problems yield richer multiplicity patterns. For
\[
M(t)=a+bt^{\theta-1},\qquad a,b>0,\ \theta\ge1,
\]
the bounded-domain fractional Kirchhoff–Choquard problem has at least one positive solution for small \(\lambda\) when \(2<q<2\theta\), at least one positive solution for large \(\lambda\) when \(2\theta\le q<2_s^*\), at least two positive solutions for sufficiently large \(\lambda\) in the superlinear case \(\theta\ge 2^*_{\mu,s}\), and at least two positive solutions in the concave regime \(1<q<2\) under the stated restrictions on \(\mu\) [2203.06471]. For the degenerate singular critical problem with prototype \(M(t)=t^{\theta-1}\), there exists \(\lambda_0>0\) such that for every \(\lambda\in(0,\lambda_0)\) the problem has at least two distinct positive weak solutions, and every weak solution is bounded and belongs to \(L^\infty(\Omega)\cap C^{0,s}(\mathbb R^N)\) when \(\mu<\min\{N,4s\}\) [2106.10856].

The Heisenberg-group critical equation exhibits a parameter-dependent dichotomy. In the critical case \(\tau=Q_\lambda^*/p\), the problem has infinitely many solutions for \(\mu\) sufficiently large by means of the Krasnosel’skii genus theorem and a Palais–Smale condition above a threshold
\[
\mu > \frac{2p}{m_0}H_{Q_\lambda^*}^{-Q_\lambda^*/p}
\]
[2401.08310]. In the subcritical case \(\tau\in\left(1,\frac{Q_\lambda^*}{p}\right)\), for every \(m\in\mathbb N\) there exists \(\mu_m>0\) such that for all \(0<\mu<\mu_m\), the problem has at least \(m\) pairs of solutions [2401.08310].

Normalized solutions introduce a different existence paradigm. For the coupled constrained system with \(M(t)=a+bt\), the \(L^2\)-subcritical case
\[
2<p+q<4+\frac{4s-2\mu}{N}
\]
yields a normalized ground state with negative energy
\[
m_\alpha(d_1,d_2)<0,
\]
negative multipliers \(\tilde\lambda_1,\tilde\lambda_2<0\), and positivity plus radial monotonicity of the solution [2509.07597]. In the \(L^2\)-supercritical case
\[
4+\frac{8s-2\mu}{N}<p+q<2^*_{\mu,s}\text{ (as used in the paper)},
\]
the solution is obtained by a mountain-pass argument at a positive energy level
\[
\sigma_\alpha(d_1,d_2)\in \left( 0,\, \frac{abHL^3}{2}+\frac{b^3HL^6}{12} +\frac{2}{3}\left(\frac{b^2HL^4}{4}+aHL\right)^{3/2} \right)
\]
[2509.07597].

Odd nonlinearities and radial symmetry frequently strengthen the results. In the variable-order problem without Ambrosetti–Rabinowitz, oddness of \(f\) yields a sequence of nontrivial weak solutions with energies tending to \(+\infty\) by the Fountain theorem and another sequence with negative critical values converging to \(0\) by the Dual Fountain theorem [2005.09221]. In the variable-exponent whole-space problem, radial symmetry is used to obtain a nontrivial radial weak solution [2005.00617].

## 6. Extensions, special regimes, and conceptual scope

The literature shows that fractional Kirchhoff–Choquard theory is not limited to the standard scalar polynomial model. One direction is geometric generalization: the Heisenberg-group equation replaces Euclidean scaling by sub-Riemannian homogeneous structure, changing the critical exponent to \(Q_\lambda^*=(2Q-\lambda)/(Q-2s)\) and requiring a fractional concentration–compactness principle on \(\mathbb H^N\) [2401.08310]. Another direction is weighted singularity: the whole-space weighted equation combines a Hardy-type potential, a Sobolev-critical weighted local term, and a weighted Choquard convolution term with critical singular weights [2410.05185].

A further extension concerns the choice of operator and growth law. Variable-exponent and variable-order models employ spaces \(W^{s,p(\cdot,\cdot)}\) and \(W^{s(\cdot,\cdot),p(\cdot,\cdot)}\), together with variable HLS inequalities and compact radial embeddings of Strauss or Lions type [2005.00617], [2005.09221]. Exponential-growth problems replace algebraic superlinearity by nonlinearities behaving like \(\exp(|u|^{n/(n-s)})\), and the key compactness threshold is then governed by the fractional Trudinger–Moser inequality rather than a polynomial Sobolev exponent [1908.11285]. Magnetic and semiclassical variants add a fractional magnetic operator \((-\Delta)^s_A\), a magnetic Gagliardo seminorm \([u]_{s,A}\), and small-\(\varepsilon\) concentration phenomena, while retaining the Kirchhoff–Choquard structure [2003.05194].

Several recurring misconceptions are corrected by the published results. First, the subject is not confined to nondegenerate Kirchhoff coefficients; degenerate cases with \(M(0)=0\) are treated explicitly in both scalar singular problems and magnetic semiclassical problems [2106.10856], [2003.05194]. Second, Ambrosetti–Rabinowitz is not structurally necessary for existence theory: both constant-exponent and variable-order models establish existence and multiplicity without that assumption, replacing it with superlinearity and monotonicity conditions tailored to Cerami compactness or Nehari geometry [1808.07996], [2005.09221]. Third, “system” need not mean only two-component coupling. Some papers use a genuinely coupled pair \((u,v)\) [2509.07597], while others study a single equation with multiple interacting nonlocal mechanisms; a plausible implication is that the term has acquired a broader structural meaning in the recent literature.

Taken together, these works define fractional Kirchhoff–Choquard systems as a research area centered on three analytic themes: nonlocal diffusion, global Kirchhoff coupling, and nonlocal HLS-type interaction. The principal technical questions are compactness at critical growth, the effect of degeneracy in \(M\), and the construction of critical points under symmetry, mass constraints, singularity, or weighted geometry. The available results cover existence, ground states, multiplicity, regularity, normalized solutions, and semiclassical states across Euclidean, weighted, magnetic, variable-exponent, and Heisenberg-group settings [2203.06471], [2106.10856], [2509.07597].

Source: https://www.emergentmind.com/topics/fractional-kirchhoff-choquard-system