Sustainable thermal power equipment supplier selection by Einstein prioritized linear Diophantine fuzzy aggregation operators

1.   Introduction

In this paper we deal with the following fractional Choquard equation

where $N\geq 2$, $\mu > 0$, $s \in (0, 1)$, $\alpha \in (0, N)$, $(-\Delta)^s$ and $I_{\alpha}$ denote respectively the fractional Laplacian and the Riesz potential defined by

where $C_{N, s}: = \frac{4^s \Gamma(\frac{N+2s}{2})}{\pi^{N/2}|\Gamma(-s)|}$ and $A_{N, \alpha} : = \frac{\Gamma(\frac{N-\alpha}{2})}{2^\alpha\pi^{N/2} \Gamma(\frac{\alpha}{2}) }$ are two suitable positive constants and the integral is in the principal value sense. Finally $F: \mathbb{R} \to \mathbb{R}$, $F' = f$ is a nonlinearity satisfying general assumptions specified below.

When dealing with double nonlocalities, important applications arise in the study of exotic stars: minimization properties related to (1.1) play indeed a fundamental role in the mathematical description of the gravitational collapse of boson stars ^[31,53] and white dwarf stars ^[37]. In fact, the study of the ground states to (1.1) gives information on the size of the critical initial conditions for the solutions of the corresponding pseudo-relativistic equation ^[48]. Moreover, when $s = \frac{1}{2}$, $N = 3$, $\alpha = 2$ and $F(t) = \frac{1}{r} |t|^r$, we obtain

related to the well-known massless boson stars equation ^[29,39,50], where the pseudorelativistic operator $\sqrt{-\Delta + m}$ collapses to the square root of the Laplacian. Other applications can be found in relativistic physics and in quantum chemistry ^[1,22,38] and in the study of graphene ^[56], where the nonlocal nonlinearity describes the short time interactions between particles.

In the limiting local case $s = 1$, when $N = 3$, $\alpha = 2$ and $F(t) = \frac{1}{2} |t|^2$, the equation has been introduced in 1954 by Pekar in ^[63] to describe the quantum theory of a polaron at rest. Successively, in 1976 it was arisen in the work ^[51] suggested by Choquard on the modeling of an electron trapped in its own hole, in a certain approximation to Hartree-Fock theory of one-component plasma (see also ^[32,33,69]). In 1996 the same equation was derived by Penrose in his discussion on the self-gravitational collapse of a quantum mechanical wave-function ^{[59,64,65,66]} (see also ^[70,71]) and in that context it is referred as Schrödinger-Newton system. Variational methods were also employed to derive existence and qualitative results of standing wave solutions for more generic values of $\alpha \in (0, N)$ and of power type nonlinearities $F(t) = \frac{1}{r} |t|^r$ ^[60] (see also ^{[14,20,49,58,61,62]}). The case of general functions $F$, almost optimal in the sense of Berestycki-Lions ^[5], has been treated in ^[19,61].

The fractional power of the Laplacian appearing in (1.1), when $s \in (0, 1)$, has been introduced instead by Laskin ^[47] as an extension of the classical local Laplacian in the study of nonlinear Schrödinger equations, replacing the path integral over Brownian motions with Lévy flights. This operator arises naturally in many contexts and concrete applications in various fields, such as optimization, finance, crystal dislocations, charge transport in biopolymers, flame propagation, minimal surfaces, water waves, geo-hydrology, anomalous diffusion, neural systems, phase transition and Bose-Einstein condensation (see ^{[6,25,31,45,46,55]} and references therein). Equations involving the fractional Laplacian together with local nonlinearities have been largely investigated, and some fundamental contributions can be found in ^[9,10,30]. In particular, existence and qualitative properties of the solutions for general classes of fractional NLS equations with local sources have been studied in ^{[7,11,28,41,42]}.

Mathematically, doubly nonlocal equations have been treated in ^[23,24] in the case of pure power nonlinearities (see also ^[13] for some orbital stability results and ^[12] for a Strichartz estimates approach), obtaining existence and qualitative properties of the solutions. Other results can be found in ^[4,57,67] for superlinear nonlinearities, in ^[36] for $L^2$-supercritical Cauchy problems, in ^[35] for bounded domains and in ^[72] for concentration phenomena with strictly noncritical and monotone sources.

In the present paper we address the study of (1.1) when $f$ satisfies the following set of assumptions of Berestycki-Lions type ^[5]:

(f1) $f \in C(\mathbb{R}, \mathbb{R})$;

(f2) we have

(f3) $F(t) = \int_0^t f(\tau) d\tau$ satisfies

(f4) there exists $t_0 \in \mathbb{R}$, $t_0 \neq 0$ such that $F(t_0) \neq 0$.

We observe that (f3) implies that we are in a noncritical setting: indeed the exponents $\frac{N+\alpha}{N}$ and $\frac{N+\alpha}{N-2s}$ have been addressed in ^[60] as critical for Choquard-type equations when $s = 1$, and then generalized to $s\in (0, 1)$ in ^[23]; we will assume the noncriticality in order to obtain the existence of a solution, while most of the qualitative results will be given in a possibly critical setting. This kind of general nonlinearities include some particular cases such as pure powers $f(t) \sim t^r$, cooperating powers $f(t) \sim t^r + t^h$, competing powers $f(t)\sim t^r - t^h$ and saturable functions $f(t) \sim \frac{t^3}{1+t^2}$ (which arise, for instance, in nonlinear optics ^[27]).

We deal first with existence of a ground state for (1.1), obtaining the following result.

Theorem 1.1. Assume (f1)–(f4). Then there exists a radially symmetric weak solution $u$ of (3.1), which satisfies the Pohozaev identity:

This solution is of Mountain Pass type and minimizes the energy among all the solutions satisfying (1.2).

We refer to Section 3 for the precise meaning of weak solution, of Mountain Pass type and energy, according to a variational formulation of the problem.

We point out some difficulties which arise in this framework. Indeed, the presence of the fractional power of the Laplacian does not allow to use the fact that every solution satisfies the Pohozaev identity to conclude that a Mountain Pass solution is actually a (Pohozaev) ground state, as in ^[44] (see Remark 3.5). On the other hand, the presence of the Choquard term, which scales differently from the $L^2$-norm term, does not allow to implement the classical minimization argument by ^[5,21]. Finally, the nonhomogeneity of the nonlinearity $f$ obstructs the minimization approach of ^[23,61]. Thus, we need a new approach to get existence of solutions, in the spirit of ^[16,17,18].

Under (f1)–(f4) it is moreover possible to state the existence of a constant sign solution (see Proposition 3.6). This motivates the investigation of qualitative properties for general positive solutions; in this case we consider weaker or stronger assumptions in substitution to (f1)–(f3), depending on the result. In particular, we observe that (f1)–(f2) alone imply

and

where we notice that the last inequality is weaker than (f3); some of the qualitative results are still valid when $F$ has this possible critical growth. Consider finally the following stronger assumption in the origin:

$(f5)$ $\limsup_{t \to 0} \frac{|tf(t)|}{|t|^{2}} < +\infty$,

and observe that

The main qualitative results that we obtain are the following ones.

Theorem 1.2. Assume (f1)–(f2). Let $u\in H^s(\mathbb{R}^N)$ be a weak positive solution of (1.1).Then $u \in L^1(\mathbb{R}^N) \cap L^{\infty}(\mathbb{R}^N)$. The same conclusion holds for generally signed solutions by assuming also (f5).

The condition in zero of the function $f$ assumed in (f5) leads also to the following polynomial decay of the solutions.

Theorem 1.3. Assume (f1)– (f2) and (f5). Let $u\in H^s(\mathbb{R}^N)$ be a positive weak solution of (1.1). Then there exists $C', C'' > 0$ such that

The previous results generalize some of the ones in ^[23] to the case of general, not homogeneous, nonlinearities; in particular, we do not even assume $f$ to satisfy Ambrosetti-Rabinowitz type conditions nor monotonicity conditions. We observe in addition that the information $u\in L^1(\mathbb{R}^N)\cap L^2(\mathbb{R}^N)$ is new even in the power-type setting: indeed in ^[23] the authors assume the nonlinearity to be not lower critical, while here we include the possibility of criticality. Moreover, we improve the results in ^[57,67] since we do not assume $f$ to be superlinear, and we have no restriction on the parameter $\alpha$. Finally, we extend some of the results in ^[61] to the fractional framework, and some of the results in ^[7] to Choquard nonlinearities.

The paper is organized as follows. We start with some notations and recalls in Section 2. In Section 3 we obtain the existence of a ground state in a noncritical setting, and in addition the existence of a positive solution. Section 4 is dedicated to the study of the boundedness of positive solutions, while in Section 5 we investigate the asymptotic decay. Finally in the Appendix A we obtain the boundedness of general signed solutions under some more restrictive assumption.

2.   Preliminaries

Let $N\geq 2$ and $s \in (0, 1)$. Recalled the definition of the fractional Laplacian ^[25]

for every $s \in (0, 1)$, we set the fractional Sobolev space as

endowed with

In particular, we consider the subspace of radially symmetric functions $H^s_r(\mathbb{R}^N)$, and recall the continuous embedding [^[25], Theorem 3.5]

for every $p \in [2, 2^*_s]$, $2^*_s = \frac{2N}{N-2s}$ critical Sobolev exponent, and the compact embedding ^[54]

for every $p \in (2, 2^*_s)$. In addition we have the following embedding of the homogeneous space [^[25], Theorem 6.5] for some $\mathcal{{S}} > 0$

Moreover the following relation with the Gagliardo seminorm holds [^[25], Proposition 3.6], for some $C(N, s) > 0$

Thanks to this last formulation, we obtain that if $u\in H^s(\mathbb{R}^N)$ and $h: \mathbb{R} \to \mathbb{R}$ is a Lipschitz function with $h(0) = 0$, then $h(u) \in H^s(\mathbb{R}^N)$. Indeed

and

We further have the following relation with the Fourier transform [^[25], Proposition 3.3]

notice that this last expression is suitable for defining the fractional Sobolev space $W^{s, p}(\mathbb{R}^N)$ also for $s\geq 1$ and $p\geq 1$, by ^[28]

Finally, set $\alpha \in (0, N)$, we recall the following standard estimates for the Riesz potential [^[52], Theorem 4.3].

Proposition 2.1 (Hardy-Littlewood-Sobolev inequality). Let $\alpha \in (0, N)$, and let $r, \, h \in (1, +\infty)$ be such that $\frac{1}{r} - \frac{1}{h} = \frac{\alpha}{N}$. Then the map

is continuous. In particular, if $r, \, t \in (1, +\infty)$ verify $\frac 1 r + \frac 1 t = \frac{N+\alpha}N$, then there exists a constant $C = C(N, \alpha, r, t) > 0$ such that

for all $g\in L^r(\mathbb{R}^N)$ and $h\in L^t(\mathbb{R}^N)$.

3.   Existence of ground states

In this section we search for solutions to the fractional Choquard equation

by variational methods on the subspace of radially symmetric functions $H^s_r(\mathbb{R}^N)$. We recall that $F' = f$ and we assume (f1)–(f2) in order to have well defined functionals. We set $\mathcal{{D}}: H^s_r(\mathbb{R}^N) \to \mathbb{R}$ as

and define the $C^1$-functional $\mathcal{{J}}_{\mu}: H^s_r(\mathbb{R}^N) \to \mathbb{R}$ associated to (3.1) by

We notice that, by the Principle of Symmetric Criticality of Palais, the critical points of $\mathcal{{J}}_{\mu}$ are weak solutions of (3.1). Moreover, inspired by the Pohozaev identity

we define also the Pohozaev functional $\mathcal{{P}}_{\mu}: H^s_r(\mathbb{R}^N) \to \mathbb{R}$ by

Furthermore we introduce the set of paths

and the Mountain Pass (MP for short) value

Finally we set

the least energy of $\mathcal{{J}}_{\mu}$ on the Pohozaev set.

Remark 3.1. Since of key importance in the good definition of the functionals, as well as in bootstrap argument in the rest of the paper, we write here in which spaces lie the considered quantities. Let $u\in H^s(\mathbb{R}^N) \subset L^2(\mathbb{R}^N) \cap L^{2^*_s}(\mathbb{R}^N)$. By (f2) we have

Thus by the Hardy-Littlewood-Sobolev inequality we obtain

Finally, by the Hölder inequality, we have

In particular we observe that $(I_{\alpha}*F(u))f(u)$ does not lie in $L^2(\mathbb{R}^N)$, generally. On the other hand, if $\varphi \in H^s(\mathbb{R}^N) \subset L^2(\mathbb{R}^N) \cap L^{2^*_s}(\mathbb{R}^N)$, we notice that the found summability of $(I_{\alpha}*F(u))f(u)$ is enough to have

well defined.

We present now an existence result for (3.1).

Theorem 3.2. Assume ${\rm{(f1)–(f4)}}$. Let $\mu > 0$ be fixed. Then there exists a Mountain Pass solution $u$ of (3.1), that is

Moreover, the found solution satisfies the Pohozaev identity

Proof. We split the proof in some steps.

Step 1. We first show that $\mathcal{{J}}_{\mu}$ satisfies the Palais-Smale-Pohozaev condition at every level $b \in \mathbb{R}$, that is each sequence $u_n$ in $H^s_r(\mathbb{R}^N)$ satisfying

converges up to a subsequence. Indeed (3.4) and (3.6) imply

Thus we obtain that $b\geq 0$ and $u_n$ is bounded in $H^s_r(\mathbb{R}^N)$.

Step 2. After extracting a subsequence, denoted in the same way, we may assume that $u_n \rightharpoonup u_0$ weakly in $H^s_r(\mathbb{R}^N)$. Taking into account the assumptions (f1)–(f3), we obtain

and

Thus we derive that $\langle \mathcal{{J}}_{\mu}'(u_n), u_n \rangle \to 0$ and $\langle \mathcal{{J}}_{\mu}'(u_n), u_0 \rangle \to 0$, and hence

which implies $u_n \to u_0$ strongly in $H^s_r(\mathbb{R}^N)$.

Step 3. Denote by

the sublevel of $\mathcal{{J}}_{\mu}$ and by

the set of critical points of $\mathcal{{J}}_{\mu}$ satisfying the Pohozaev identity. Then, by Steps 1–2, $K_b$ is compact. Arguing as in [^[40], Proposition 4.5] (see also [^[43], Proposition 3.1 and Corollary 4.3]), we obtain for any $b\in \mathbb{R}$, $\bar{{\varepsilon}} > 0$ and any $U$ open neighborhood of $K_b$, that there exist an ${\varepsilon} \in (0, \bar {\varepsilon})$ and a continuous map $\eta:[0, 1] \times H^s_r(\mathbb{R}^N) \to H^s_r(\mathbb{R}^N)$ such that

$(1^o)$ $\eta(0, u) = u \quad \forall u \in H^s_r(\mathbb{R}^N)$;

$(2^o)$ $\eta(t, u) = u \quad \forall (t, u) \in [0, 1] \times [\mathcal{{J}}_{\mu}\leq b- \bar{\varepsilon}]$;

$(3^o)$ $\mathcal{{J}}_{\mu}(\eta(t, u)) \leq \mathcal{{J}}_{\mu}(u) \quad \ \forall (t, u) \in [0, 1] \times H^s_r(\mathbb{R}^N)$;

$(4^o)$ $\eta (1, [\mathcal{{J}}_{\mu} \leq b+ {\varepsilon}] \setminus U) \subset [\mathcal{{J}}_{\mu} \leq b- {\varepsilon}]$;

$(5^o)$ $\eta (1, [\mathcal{{J}}_{\mu} \leq b+ {\varepsilon}]) \subset [\mathcal{{J}}_{\mu} \leq b - {\varepsilon}] \cup U$;

$(6^o)$ if $K_b = \emptyset$, then $\eta (1, [\mathcal{{J}}_{\mu} \leq b+ {\varepsilon}]) \subset [\mathcal{{J}}_{\mu} \leq b- {\varepsilon}].$

Step 4. By exploiting (f4) and arguing as in [^[61], Proposition 2.1], we obtain the existence of a function $v \in H^s_r(\mathbb{R}^N)$ such that $\mathcal{{D}}(v) > 0$. Thus defined $\gamma(t) : = v(\cdot/t)$ for $t > 0$ and $\gamma(0): = 0$ we have $\mathcal{{J}}(\gamma(t)) < 0$ for $t$ large and $\mathcal{{J}}(\gamma(t)) > 0$ for $t$ small; this means, after a suitable rescaling, that $l(\mu)$ is finite and strictly positive. In particular we observe that $0\notin K_{l(\mu)}$.

Step 5. By applying the deformation result at level $b = l(\mu) > 0$, the existence of a Mountain Pass solution $u$ is then obtained classically. Moreover, $u \in K_{l(\mu)}$ by construction, thus $u \not \equiv 0$ and $\mathcal{{P}}_{\mu}(u) = 0$.

We prove now that the found solution is actually a ground state over the Pohozaev set.

Proposition 3.3. The Mountain Pass level and the Pohozaev minimum level coincide, that is

In particular, the solution found in Theorem 3.2 is a Pohozaev minimum.

Proof. Let $u \in H^s_r(\mathbb{R}^N)\setminus\{0\}$ such that $\mathcal{{P}}_{\mu}(u) = 0$; observe that $\mathcal{{D}}(u) > 0$. We define $\gamma(t): = u(\cdot/t)$ for $t \neq 0$ and $\gamma(0): = 0$ so that $t\in (0, +\infty) \mapsto \mathcal{{J}}_{\mu}(\gamma(t))$ is negative for large values of $t$, and it attains the maximum in $t = 1$. After a suitable rescaling we have $\gamma \in \Gamma_{\mu}$ and thus

Passing to the infimum in Eq (3.7) we have $p(\mu) \geq l(\mu)$. Let now $\gamma \in \Gamma_{\mu}$. By definition we have $\mathcal{{J}}_{\mu}(\gamma(1)) < 0$, thus by

we obtain $\mathcal{{P}}_{\mu}(\gamma(1)) < 0$. In addition, since $\mathcal{{D}}(u) = o(\|u_\|{H^s}^2)$ as $u \to 0$ and $\gamma(t)\to 0$ as $t \to 0$ in $H^s_r(\mathbb{R}^N)$, we have

Thus there exists a $t^*$ such that $\mathcal{{P}}_{\mu}(\gamma(t^*)) = 0$, and hence

passing to the infimum we come up with $p(\mu) \leq l(\mu)$, and hence the claim.

Proof of Theorem 1.1. We obtain the result by matching Theorem 3.2 and Proposition 3.3.

We pass to investigate more in details Pohozaev minima, showing that it is a general fact that they are solutions of the Eq (3.1).

Proposition 3.4. Every Pohozaev minimum is a solution of (3.1), i.e.,

imply

As a consequence

Proof. Let $u$ be such that $\mathcal{{J}}_{\mu}(u) = p(\mu)$ and $\mathcal{{P}}_{\mu}(u) = 0$. In particular, considered $\gamma(t) = u(\cdot/t)$, we have that $\mathcal{{J}}_{\mu}(\gamma(t))$ is negative for large values of $t$ and its maximum value is $p(\mu)$ attained only in $t = 1$.

Assume by contradiction that $u$ is not critical. Let $I: = [1-\delta, 1+\delta]$ be such that $\gamma(I) \cap K_{p(\mu)} = \emptyset$, and set $\bar{{\varepsilon}} : = p(\mu) - \max_{t \notin I} \mathcal{{J}}_{\mu}(\gamma(t)) > 0$. Let now $U$ be a neighborhood of $K_{p(\mu)}$ verifying $\gamma(I) \cap U = \emptyset$: by the deformation lemma presented in the proof of Theorem 3.2 there exists an $\eta:[0, 1]\times H^s_r(\mathbb{R}^N) \to H^s_r(\mathbb{R}^N)$ at level $p(\mu)\in \mathbb{R}$ with properties $(1^o)$-$(6^o)$. Define then $\tilde{\gamma}(t): = \eta(1, \gamma(t))$ a deformed path.

For $t \notin I$ we have $\mathcal{{J}}_{\mu}(\gamma(t)) < p(\mu) - \bar{{\varepsilon}}$, and thus by $(2^o)$ we gain

Let now $t \in I$: we have $\gamma(t) \notin U$ and $\mathcal{{J}}_{\mu}(\gamma(t))\leq p(\mu) \leq p(\mu)+ {\varepsilon}$, thus by $(4^o)$ we obtain

Joining (3.8) and (3.9) we have

which is an absurd, since after a suitable rescaling it results that $\tilde{\gamma} \in \Gamma_{\mu}$, thanks to $(3^o)$.

Remark 3.5. We point out that it is not known, even in the case of local nonlinearities ^[7], if

On the other hand, by assuming that every solution of (3.1) satisfies the Pohozaev identity (see e.g., [^[67], Proposition 2] and [^[23], Eq (6.1)]), the claim holds true.

We show now that, under the same assumptions of Theorem 3.2, we can find a solution with constant sign.

Proposition 3.6. Assume ${\rm{(f1)–(f4)}}$ and that $F\not \equiv 0$ on $(0, +\infty)$ (i.e., $t_0$ in assumption $\text{(f4)}$ can be chosen positive). Then there exists a positive radially symmetric solution of (3.1), which is minimum over all the positive functions on the Pohozaev set.

Proof. Let us define

We have that $g$ still satisfies (f1)–(f4). Thus, by Theorem 3.2 there exists a solution $u$ of

where $G(t): = \int_0^t g(\tau) d \tau$. We show now that $u$ is positive. We start observing the following: by (2.2) we have

thus

In particular, written $u = u_+ - u_-$, by the previous argument we have $u_- = \frac{|u|-u}{2} \in H^s_r(\mathbb{R}^N)$. Thus, chosen $u_-$ as test function, we obtain

By definition of $g$ and (2.2) we have

Splitting the domain, we gain

Thus we obtain that the left-hand side of (3.10) is sum of non positive pieces, thus $u_- \equiv 0$, that is $u\geq 0$. Hence $g(u) = f(u)$ and $G(u) = F(u)$, which imply that $u$ is a (positive) solution of (3.1).

4.   Regularity

In this section we prove some regularity results for (3.1). We split the proof of Theorem 1.2 in different steps.

We start from the following lemma, that can be found in [^[61], Lemma 3.3].

Lemma 4.1 (^[61]). Let $N \geq 2$ and $\alpha \in (0, N)$. Let $\lambda \in [0, 2]$ and $q, r, h, k \in [1, +\infty)$ be such that

Let $\theta \in (0, 2)$ satisfying

Let $H \in L^h(\mathbb{R}^N)$, $K \in L^k(\mathbb{R}^N)$ and $u \in L^q(\mathbb{R}^N) \cap L^r(\mathbb{R}^N)$. Then

for some $C > 0$ (depending on $\theta$).

By a proper use of Lemma 4.1 we obtain now an estimate on the Choquard term depending on $H^s$-norm of the function.

Lemma 4.2. Let $N\geq 2$, $s \in (0, 1)$ and $\alpha \in (0, N)$. Let moreover $\theta \in (\frac{\alpha}{N}, 2- \frac{\alpha}{N})$ and $H, K \in L^{\frac{2N}{\alpha}}(\mathbb{R}^N) + L^{\frac{2N}{\alpha+2s}}(\mathbb{R}^N)$. Then for every ${\varepsilon} > 0$ there exists $C_{{\varepsilon}, \theta} > 0$ such that

for every $u \in H^s(\mathbb{R}^N)$.

Proof. Observe that $2-\theta \in (\frac{\alpha}{N}, 2- \frac{\alpha}{N})$ as well. We write

We split $\int_{ \mathbb{R}^N} \left(I_{\alpha}*\big(H|u|^{\theta}\big)\right)K|u|^{2-\theta} \, dx$ in four pieces and choose

in Lemma 4.1, to obtain

Recalled that $\frac{2N}{N-2s} = 2^*_s$ and the Sobolev embedding (2.1), we obtain

where $\lesssim$ denotes an inequality up to a constant. We want to show now that, since $\frac{2N}{\alpha} > \frac{2N}{\alpha+2s}$, we can choose the decomposition of $H$ and $K$ such that the $L^{\frac{2N}{\alpha+2s}}$-pieces are arbitrary small (see [^[8], Lemma 2.1]). Indeed, let

be a first decomposition. Let $M > 0$ to be fixed, and write

Since $H_2 \chi_{\{|H_2| \leq M\}} \in L^{\frac{2N}{\alpha+2s}}(\mathbb{R}^N) \cap L^{\infty}(\mathbb{R}^N)$ and $\frac{2N}{\alpha} \in (\frac{2N}{\alpha+2s}, \infty)$, we have $H_2 \chi_{\{|H_2| \leq M\}} \in L^{\frac{2N}{\alpha}}(\mathbb{R}^N)$, and thus

On the other hand

which can be made arbitrary small for $M\gg 0$. In particular we choose the decomposition so that

and thus

In the last term of (4.1) we use the generalized Young's inequality $ab \leq \frac{\delta}{2} a^2 + \frac{1}{2\delta} b^2$, with

so that

Merging the pieces, we have the claim.

The following technical result can be found in [^[35], Lemma 3.5].

Lemma 4.3 (^[35]). Let $a, b \in \mathbb{R}$, $r \geq 2$ and $k\geq 0$. Set $T_k: \mathbb{R} \to [-k, k]$ the truncation in $k$, that is

and write $a_k: = T_k(a)$, $b_k: = T_k(b)$. Then

Notice that the (optimal) Sobolev embedding tells us that $H^s(\mathbb{R}^N) \hookrightarrow L^{2^*_s}(\mathbb{R}^N)$. In the following we show that $u$ belongs to some $L^r(\mathbb{R}^N)$ with $r > 2^*_s = \frac{2N}{N-2s}$; we highlight that we make no use of the Caffarelli-Silvestre $s$-harmonic extension method, and work directly in the fractional framework.

Proposition 4.4. Let $H, K \in L^{\frac{2N}{\alpha}}(\mathbb{R}^N) + L^{\frac{2N}{\alpha+2s}}(\mathbb{R}^N)$. Assume that $u \in H^s(\mathbb{R}^N)$ solves

in the weak sense. Then

Moreover, for each of these $r$, we have

with $C_r > 0$ not depending on $u$.

Proof. By Lemma 4.2 there exists $\lambda > 0$ (that we can assume large) such that

Let us set

and observe that

and

We thus define the bilinear form

for every $\varphi, \psi \in H^s(\mathbb{R}^N)$. Since, by (4.3) and (4.2), we have

for each $\varphi \in H^s(\mathbb{R}^N)$, we obtain that $a_n$ is coercive. Set

we obtain by Lax-Milgram theorem that, for each $n \in \mathbb{N}$, there exists a unique $u_n \in H^s(\mathbb{R}^N)$ solution of

that is

in the weak sense; moreover the theorem tells us that

(since $1/2$ appears as coercivity coefficient in (4.4)), and thus $u_n$ is bounded. Hence $u_n \rightharpoonup \bar{u}$ in $H^s(\mathbb{R}^N)$ up to a subsequence for some $\bar{u}$. This means in particular that $u_n \to \bar{u}$ almost everywhere pointwise.

Thus we can pass to the limit in

we need to check only the Choquard term. We first see by the continuous embedding that $u_n \rightharpoonup \bar{u}$ in $L^q(\mathbb{R}^N)$, for $q \in [2, 2^*_s]$. Split again $H = H^*+ H_*$, $K = K^* + K_*$ and work separately in the four combinations; we assume to work generally with $\tilde{H} \in \{H^*, H_*\}$, $\tilde{H}\in L^{\beta}(\mathbb{R}^N)$ and $\tilde{K} \in \{K^*, K_*\}$, $\tilde{K}\in L^{\gamma}(\mathbb{R}^N)$, where $\beta, \gamma \in \{ \frac{2N}{\alpha}, \frac{2N}{\alpha+2s}\}$. Then one can easily prove that $\tilde{H}_n u_n \rightharpoonup \tilde{H} \bar{u}$ in $L^r(\mathbb{R}^N)$ with $\frac{1}{r} = \frac{1}{\beta} + \frac{1}{q}$. By the continuity and linearity of the Riesz potential we have $I_{\alpha} * (H_n u_n) \rightharpoonup I_{\alpha} * (H \bar{u})$ in $L^h(\mathbb{R}^N)$, where $\frac{1}{h} = \frac{1}{r} - \frac{\alpha}{n}$. As before, we obtain $\left(I_{\alpha}*\big(H_n u_n\big)\right) K_n \rightharpoonup \left(I_{\alpha}*\big(H \bar{u}\big)\right) K$ in $L^k(\mathbb{R}^N)$, where $\frac{1}{k} = \frac{1}{\gamma} + \frac{1}{h}$. Simple computations show that if $\beta = \gamma = \frac{2N}{\alpha}$ and $q = 2$, then $k' = 2$; if $\beta = \frac{2N}{\alpha}$, $\gamma = \frac{2N}{\alpha+2s}$ (or viceversa) and $q = 2$, then $k' = 2^*_s$; if $\beta = \gamma = \frac{2N}{\alpha+2s}$ and $q = 2^*_s$, then $k' = 2^*_s$. Therefore $H^s(\mathbb{R}^N) \subset L^{k'}(\mathbb{R}^N)$ and we can pass to the limit in all the four pieces, obtaining

Therefore, $\bar{u}$ satisfies

as well as $u$. But we can see this problem, similarly as before, with a Lax-Milgram formulation and obtain the uniqueness of the solution. Thus $\bar{u} = u$ and hence

and almost everywhere pointwise. Let now $k\geq 0$ and write

where $T_k$ is the truncation introduced in Lemma 4.3. Let $r\geq 2$. We have $|u_{n, k}|^{r/2} \in H^s(\mathbb{R}^N)$, by exploiting (2.2) and the fact that $h(t): = (T_k(t))^{r/2}$ is a Lipschitz function with $h(0) = 0$. By (2.2) and by Lemma 4.3 we have

Set

it results that $\varphi \in H^s(\mathbb{R}^N)$, since again $h(t): = T_k(t) |T_k(t)|^{r-2}$ is a Lipschitz function with $h(0) = 0$. Thus we can choose it as a test function in (4.5) and obtain, by polarizing the identity (2.2),

and since $u_n \varphi \geq |u_{n, k}|^r$ we gain

Focus on the Choquard term on the right-hand side. We have

Focus on $(I)$. Consider $r \in [2, \frac{2N}{\alpha})$, so that $\theta: = \frac{2}{r} \in (\frac{\alpha}{N}, 2-\frac{\alpha}{N})$. Choose moreover $v: = |u_{n, k}|^{r/2} \in H^s(\mathbb{R}^N)$ and ${\varepsilon}^2: = \frac{2(r-1)}{r^2} > 0$. Thus, observed that if a function belongs to a sum of Lebesgue spaces then its absolute value does the same ([^[3], Proposition 2.3]), by Lemma 4.2 we obtain

Focus on $(II)$. Assuming $r < \min\{\frac{2N}{\alpha}, \frac{2N}{N-2s}\}$, we have $u_n \in L^r(\mathbb{R}^N)$ and $H_n \in L^{\frac{2N}{\alpha}}(\mathbb{R}^N)$, thus

for the Hölder inequality. Similarly

Thus, since $\frac{1}{a} + \frac{1}{b} = \frac{N+\alpha}{N}$, we have by the Hardy-Littlewood-Sobolev inequality (see Proposition 2.1) that

With respect to $k$, the second factor on the right-hand side is bounded, while the first factor goes to zero thanks to the dominated convergence theorem, thus

Joining (4.6), (4.8), (4.9), (4.10) we obtain

That is, by Sobolev inequality (2.1)

Letting $k\to +\infty$ by the monotone convergence theorem (since $u_{n, k}$ are monotone with respect to $k$ and $u_{n, k} \to u_n$ pointwise) we have

and thus $u_n \in L^{\frac{r}{2}2^*_s}(\mathbb{R}^N)$. Notice that $\frac{r}{2} \in \big[1, \min\{\frac{N}{\alpha}, \frac{N}{N-2s}\}\big)$. If $N-2s < \alpha$ we are done. Otherwise, set $r_1: = r$, we can now repeat the argument with

Again, if $\frac{2N}{\alpha} < 2\left(\frac{N}{N-2s}\right)^2$ we are done, otherwise we repeat the argument. Inductively, we have

thus $\frac{2N}{\alpha} < 2\left(\frac{N}{N-2s}\right)^m$ after a finite number of steps. For such $r = r_m$, consider again (4.11): by the almost everywhere convergence of $u_n$ to $u$ and Fatou's lemma

Being $u_n$ equibounded in $H^s(\mathbb{R}^N)$ and thus in $L^{2^*_s}(\mathbb{R}^N)$, by the iteration argument we have that it is equibounded also in $L^r(\mathbb{R}^N)$; in particular, the bound is given by $\|{u}\|_2$ times a constant $C(r)$. Thus the right-hand side is a finite quantity, and we gain $u \in L^{\frac{r}{2}2^*_s}(\mathbb{R}^N)$, which is the claim.

The following Lemma states that $I_{\alpha}*g \in L^{\infty}(\mathbb{R}^N)$ whenever $g$ lies in $L^q(\mathbb{R}^N)$ with $q$ in a neighborhood of $\frac{N}{\alpha}$ (in particular, it generalizes Proposition 2.1 to the case $h = \infty$ and $r \approx \frac{N}{\alpha}$).

In addition, it shows the decay at infinity of the Riesz potential, which will be useful in Section 5.

Proposition 4.5. Assume that $\text{(f1)–(f2)}$ hold. Let $u\in H^s(\mathbb{R}^N)$ be a solution of (3.1). Then $u\in L^q(\mathbb{R}^N)$ for $q \in \big[2, \frac{N}{\alpha} \frac{2N}{N-2s}\big)$, and

that is, continuous and zero at infinity. In particular,

and

Proof. We first check to be in the assumptions of Proposition 4.4. Indeed, by (f1)–(f2) and the fact that $u\in H^s(\mathbb{R}^N)\subset L^2(\mathbb{R}^N) \cap L^{2^*_s}(\mathbb{R}^N)$ we obtain that

lie in $L^{\frac{2N}{\alpha}}(\mathbb{R}^N) + L^{\frac{2N}{\alpha+2s}}(\mathbb{R}^N)$, since bounded by functions in this sum space (see e.g., [^[3], Proposition 2.3]). Now by Proposition 4.4 we have $u \in L^q(\mathbb{R}^N)$ for $q \in [2, \frac{N}{\alpha} \frac{2N}{N-2s})$.

To gain the information on the convolution, we want to use Young's Theorem, which states that if $g, h$ belong to two Lebesgue spaces with conjugate (finite) indexes, then $g*h \in C_0(\mathbb{R}^N)$. We first split

where

We need to show that $F(u) \in L^{q_1}(\mathbb{R}^N)\cap L^{q_2}(\mathbb{R}^N)$ for some $q_i$ satisfying

that is

or equivalently $q_2 < \frac{N}{\alpha} < q_1$. Recall that

Note that $u \in L^q(\mathbb{R}^N)$ for $q \in [2, \frac{N}{\alpha} \frac{2N}{N-2s})$ implies

for some $q_2 < \frac{N}{\alpha} < q_1$. Thus we have the claim.

Once obtained the boundedness of the Choquard term, we can finally gain the boundedness of the solution.

Proposition 4.6. Assume that $\text{(f1)–(f2)}$ hold. Let $u\in H^s(\mathbb{R}^N)$ be a positive solution of (3.1). Then $u\in L^{\infty}(\mathbb{R}^N)$.

Proof. By Lemma 4.5 we obtain

Thus $u$ satisfies the following nonautonomous problem, with a local nonlinearity

with $a$ bounded. In particular

where

Set $\gamma: = \max\{1, \frac{\alpha+2s}{N-2s}\} \in [1, 2^*_s)$, we thus have

Hence we are in the assumptions of [^[26], Proposition 5.1.1] and we can conclude.

We observe that a direct proof of the boundedness for generally signed solutions, but assuming also (f5), can be found in Appendix A.

Gained the boundedness of the solutions, we obtain also some additional regularity, which will be implemented in some bootstrap argument for the $L^1$-summability.

Proposition 4.7. Assume that $\text{(f1)–(f2)}$ hold. Let $u\in H^s(\mathbb{R}^N)\cap L^{\infty}(\mathbb{R}^N)$ be a weak solution of (3.1). Then $u \in H^{2s}(\mathbb{R}^N) \cap C^{0, \gamma}(\mathbb{R}^N)$ for any $\gamma \in (0, \min\{1, 2s\})$. Moreover $u$ satisfies (3.1) almost everywhere.

Proof. By Proposition 4.6, Proposition 4.5 and (f2) we have that $u\in L^{\infty}(\mathbb{R}^N)$ satisfies

where $g(x): = - \mu u(x) + (I_{\alpha}*F(u))(x) f(u(x))$. We prove first that $u\in H^{2s}(\mathbb{R}^N)$. Indeed, we already know that $f(u)$, $F(u)$ and $I_{\alpha}*F(u)$ belong to $L^{\infty}(\mathbb{R}^N)$. By Remark 3.1, we obtain

In particular,

Since $u$ is a weak solution, we have, fixed $\varphi \in H^s(\mathbb{R}^N)$,

Since $g \in L^2(\mathbb{R}^N)$, we can apply Plancharel theorem and obtain

Since $H^s(\mathbb{R}^N) = \mathcal{{F}}(H^s(\mathbb{R}^N))$ and $\varphi$ is arbitrary, we gain

By definition, we obtain $u \in H^{2s}(\mathbb{R}^N)$, which concludes the proof. Observe moreover that $\mathcal{{F}}^{-1}\big((1+|\xi|^{2s})\widehat{u}\big) = u +g \in L^2(\mathbb{R}^N) \cap L^{\infty}(\mathbb{R}^N)$, thus by definition $u \in H^{2s}(\mathbb{R}^N) \cap W^{2s, \infty}(\mathbb{R}^N)$. By the embedding [^[28], Theorem 3.2] we obtain $u \in C^{0, \gamma}(\mathbb{R}^N)$ if $2s < 1$ and $\gamma \in (0, 2s)$, while $u \in C^{1, \gamma}(\mathbb{R}^N)$ if $2s > 1$ and $\gamma \in (0, 2s-1)$ (see also [^[68], Proposition 2.9]).

It remains to show that $u$ is an almost everywhere pointwise solution. Thanks to the fact that $u\in H^{2s}(\mathbb{R}^N)$, we use again (4.13), where we can apply Plancharel theorem (that is, we are integrating by parts (4.12)) and thus

Since $\varphi \in H^s(\mathbb{R}^N)$ is arbitrary, we obtain

This concludes the proof.

We observe, by the proof, that if $s\in (\tfrac{1}{2}, 1)$, then $u \in C^{1, \gamma}(\mathbb{R}^N)$ for any $\gamma \in (0, 2s-1)$, and $u$ is a classical solution, with $(-\Delta)^s u \in C(\mathbb{R}^N)$ and equation (3.1) satisfied pointwise.

We end this section by dealing with the summability of $u$ in Lebesgue spaces $L^r(\mathbb{R}^N)$ for $r < 2$.

Remark 4.8. We start noticing that, if a solution $u$ belongs to some $L^q(\mathbb{R}^N)$ with $q < 2$, then $u\in L^1(\mathbb{R}^N)$. Assume thus $q \in (1, 2)$ and let $u \in L^q(\mathbb{R}^N) \cap L^{\infty}(\mathbb{R}^N)$, then we have

Thanks to Proposition 4.7, $u$ satisfies (3.1) almost everywhere, thus we have

which equivalently means that the Bessel operator verifies

Thus by [^[2], Theorem 1.2.4] we obtain that $u$ itself lies in the same Lebesgue space, that is

If $\frac{qN}{N + \alpha(2-q)} < 1$, we mean that $(I_{\alpha}*F(u)) f(u) \in L^1(\mathbb{R}^N) \cap L^{\infty}(\mathbb{R}^N)$, and thus $u \in L^1(\mathbb{R}^N) \cap L^{\infty}(\mathbb{R}^N)$. We convey this when we deal with exponents less than $1$.

If $q < 2$, then

and we can implement a bootstrap argument to gain $u \in L^1(\mathbb{R}^N)$. More precisely

where $q_n \to 0$ (but we stop at $1$). Thus, in order to implement the argument, we need to show that $u\in L^q(\mathbb{R}^N)$ for some $q < 2$.

We show now that $u\in L^1(\mathbb{R}^N)$. It is easy to see that, if the problem is (strictly) not lower-critical, i.e., (f2) holds together with

for some $\beta\in ({N+\alpha\over N}, {N+\alpha\over N-2s})$, then $u\in L^1(\mathbb{R}^N)$. Indeed $u\in H^s(\mathbb{R}^N)\cap L^\infty(\mathbb{R}^N)\subset L^2(\mathbb{R}^N)\cap L^\infty(\mathbb{R}^N)$ and

where ${1\over q} = {\beta\over 2}-{\alpha\over 2N}$; noticed that $q < 2$, we can implement the bootstrap argument of Remark 4.8.

We will show that the same conclusion can be reached by assuming only (f2).

Proposition 4.9. Assume that $\text{(f1)}– \text{(f2)}$ hold. Let $u\in H^s(\mathbb{R}^N)\cap L^{\infty}(\mathbb{R}^N)$ be a weak solution of (3.1). Then $u \in L^1(\mathbb{R}^N)$.

Proof of Proposition 4.9. For a given solution $u\in H^s(\mathbb{R}^N)\cap L^\infty(\mathbb{R}^N)$ we set again

Since $u\in L^2(\mathbb{R}^N)\cap L^\infty(\mathbb{R}^N)$, by (f2) we have $H$, $K\in L^{2N\over \alpha}(\mathbb{R}^N)$. For $n\in \mathbb{N}$, we set

Then we have

Since ${\rm{supp}}\; (H-H_n)\subset \big\{ |x|\leq n\big\}$ is a bounded set, we have for any $\beta \in [1, {2N\over \alpha}]$

We write our equation (3.1) as

where we introduced the function $R_n$ by

Now we consider the following linear equation:

We have the following facts:

(i) The given solution $u$ solves (4.16).

(ii) By the property (4.15) with $\beta \in (\frac{2N}{N+\alpha}, {2N\over \alpha})$, there exists $q_1 \in (1, 2)$, namely ${1\over q_1} = {1\over\beta}+{{1\over 2}}-{\alpha\over 2N}$, such that $R_n\in L^{q_1}(\mathbb{R}^N)\cap L^2(\mathbb{R}^N)$.

(iii) By the property (4.14), for any $r\in (\frac{2N}{2N-\alpha}, 2] \subset (1, 2]$

is well-defined and verifies

Here $C_{r, n}$ satisfies $C_{r, n}\to 0$ as $n\to\infty$.

We show only (iii). Since $v\in L^r(\mathbb{R}^N)$, by Hardy-Littlewood-Sobolev inequality and Hölder inequality we obtain

where $C_r > 0$ is independent of $n$, $v$. Thus by (4.14) we have $C_{r, n}: = C_r\|{H_n}\|_{2N\over \alpha}\|K_\|{2N\over \alpha}\to 0$ as $n\to\infty$.

Now we show $u\in L^{q_1}(\mathbb{R}^N)$, where $q_1\in (1, 2)$ is given in (ii). Since $((-\Delta)^s+\mu)^{-1}:\, L^r(\mathbb{R}^N)\to L^r(\mathbb{R}^N)$ is a bounded linear operator for $r\in(1, 2]$ (see [^[2], Theorem 1.2.4]), (4.16) can be rewritten as

where

By choosing $\beta \in (2, \frac{2N}{\alpha})$ we have $q_1 \in (\frac{2N}{2N-\alpha}, 2)\subset(1, 2)$, thus we observe that for $n$ large, $T_n$ is a contraction in $L^2(\mathbb{R}^N)$ and in $L^{q_1}(\mathbb{R}^N)$. We fix such an $n$.

Since $T_n$ is a contraction in $L^2(\mathbb{R}^N)$, we can see that $u\in H^s(\mathbb{R}^N)$ is a unique fixed point of $T_n$. In particular, we have

On the the other hand, since $T_n$ is a contraction in $L^{q_1}(\mathbb{R}^N)$, $(T_n^k (0))_{k = 1}^\infty$ also converges in $L^{q_1}(\mathbb{R}^N)$. Thus the limit $u$ belongs to $L^{q_1}(\mathbb{R}^N)$.

Since $q_1 < 2$ we can use the bootstrap argument of Remark 4.8 to get $u\in L^1(\mathbb{R}^N)$, and reach the claim.

5.   Asymptotic decay

We prove now the polynomial decay of the solutions. We start from two standard lemmas, whose proofs can be found for instance in [^[15], Lemma A.1 and Lemma A.3].

Lemma 5.1 (Maximum Principle). Let $\Sigma \subset \mathbb{R}^N$, possibly unbounded, and let $u\in H^s(\mathbb{R}^N)$ be a weak subsolution of

with $a > 0$, in the sense that

for every positive $\varphi \in H^s(\mathbb{R}^N)$ with ${\mathop{{\rm{supp}}}}(\varphi) \subset \mathbb{R}^N \setminus \Sigma$. Assume moreover that

Then

Lemma 5.2 (Comparison function). Let $b > 0$. Then there exists a strictly positive continuous function $W\in H^s(\mathbb{R}^N)$ such that, for some positive constants $C', C''$ (depending on $b$), it verifies

pointwise, with $r: = b^{-1/2s}$, and

We show first some conditions which imply the decay at infinity of the solutions.

Lemma 5.3. Assume that $\text{(f1)–(f2)}$ hold. Let $u$ be a weak solution of (3.1). Assume

and

Then we have

Proof. Being $u$ solution of

where $\chi \in L^{\frac{N}{2s}}(\mathbb{R}^N) \cap L^{\infty}(\mathbb{R}^N)$, we have the representation formula

where $\mathcal{{K}}$ is the Bessel kernel; we recall that $\mathcal{{K}}$ is positive, it satisfies $\mathcal{{K}}(x) \leq \frac{C}{|x|^{N+2s}}$ for $|x| \geq 1$ and $\mathcal{{K}} \in L^q(\mathbb{R}^N)$ for $q \in [1, 1 + \tfrac{2s}{N-2s})$ (see [^[28], page 1241 and Theorem 3.3]). Let us fix $\eta > 0$; we have, for $x \in \mathbb{R}^N$,

As regards the first piece

while for the second piece, fixed a whatever $q \in (1, 1 + \tfrac{2s}{N-2s})$ and its conjugate exponent $q' > \frac{N}{2s}$, we have by Hölder inequality

where the second factor can be made small for $|x| \gg 0$. Joining the pieces, we have (5.3).

We observe that the assumptions of the Lemma are fulfilled by assuming that $u$ is bounded thanks to Proposition 4.9. We are now ready to prove the polynomial decay of the solutions.

Conclusion of the proof of Theorem 1.3. Observe that, by (f5) and Lemma 5.3, we have

Thus we obtain, by applying Proposition 4.5, that

Thus, by (5.5) and the positivity of $u$, we have for some $R'\gg 0$

Similarly

Notice that we always intend differential inequalities in the weak sense, that is tested with functions in $H^s(\mathbb{R}^N)$ with supports contained in the reference domain (e.g., $\mathbb{R}^N \setminus B_{R'}$).

In addition, by Lemma 5.2 we have that there exist two positive functions $\underline{W}'$, $\overline{W}'$ and three positive constants $R''$, $C'$ and $C''$ depending only on $\mu$, such that

and

Set $R: = \max\{ R', 2R''\}$. Let $\underline{C}_1$ and $\overline{C}_1$ be some lower and upper bounds for $u$ on $B_R$, $\underline{C}_2: = \min_{B_R} \overline{W}'$ and $\overline{C}_2: = \max_{B_R} \underline{W}'$, all strictly positive. Define

so that

Thanks to the comparison principle in Lemma 5.1, and redefining $C'$ and $C''$, we obtain

By the boundedness of $u$, we obtain the claim.

We see that, for non sublinear $f$ (that is, (f5)), the decay is essentially given by the fractional operator. It is important to remark that, contrary to the limiting local case $s = 1$ (see ^[60]), the Choquard term in case of linear $f$ does not affect the decay of the solution.

Remark 5.4. We observe that the conclusion of the proof of Theorem 1.3 can be substituted by exploiting a result in ^[30]. Indeed write $V: = -(I_{\alpha}*F(u)) \frac{f(u)}{u}$, which is bounded and zero at infinity as observed in (5.4)–(5.5), and gain

Up to dividing for $\|{u}\|_2$, we may assume $\|{u}\|_2 = 1$. Thus we are in the assumptions of [^[30], Lemma C.2] and obtain, even for changing-sign solutions of (3.1),

together with

for some $C_1, C_2 > 0$.

Acknowledgments

The first and second authors are supported by PRIN 2017JPCAPN "Qualitative and quantitative aspects of nonlinear PDEs'' and by INdAM-GNAMPA. The third author is supported in part by Grant-in-Aid for Scientific Research (19H00644, 18KK0073, 17H02855, 16K13771) of Japan Society for the Promotion of Science.

Conflict of interest

All authors declare no conflicts of interest in this paper.