Quasilinear equations with natural growth in the gradients in spaces of Sobolev multipliers

We study the existence problem for a class of nonlinear elliptic equations whose prototype is of the form \(-\Delta _p u = |\nabla u|^p + \sigma \) in a bounded domain \(\Omega \subset \mathbb {R}^n\). Here \(\Delta _p\), \(p>1\), is the standard p-Laplacian operator defined by \(\Delta _p u=\mathrm{div}\, (|\nabla u|^{p-2}\nabla u)\), and the datum \(\sigma \) is a signed distribution in \(\Omega \). The class of solutions that we are interested in consists of functions \(u\in W^{1,p}_0(\Omega )\) such that \(|\nabla u|\in M(W^{1,p}(\Omega )\rightarrow L^p(\Omega ))\), a space pointwise Sobolev multipliers consisting of functions \(f\in L^{p}(\Omega )\) such that

$$\begin{aligned} \int _{\Omega } |f|^{p} |\varphi |^p dx \le C \int _{\Omega } (|\nabla \varphi |^p + |\varphi |^p) dx \quad \forall \varphi \in C^\infty (\Omega ), \end{aligned}$$

for some \(C>0\). This is a natural class of solutions at least when the distribution \(\sigma \) is nonnegative and compactly supported in \(\Omega \). We show essentially that, with only a gap in the smallness constants, the above equation has a solution in this class if and only if one can write \(\sigma =\mathrm{div}\, F\) for a vector field F such that \(|F|^{\frac{1}{p-1}}\in M(W^{1,p}(\Omega )\rightarrow L^p(\Omega ))\). As an important application, via the exponential transformation \(u\mapsto v=e^{\frac{u}{p-1}}\), we obtain an existence result for the quasilinear equation of Schrödinger type \(-\Delta _p v = \sigma \, v^{p-1}\), \(v\ge 0\) in \(\Omega \), and \(v=1\) on \(\partial \Omega \), which is interesting in its own right.

     

Estimates for large solutions of a quasilinear elliptic equation

Let \(\Omega \) be a bounded smooth domain of \(R^{n}\). We study the asymptotic behaviour of the solutions to the equation \(\triangle u-|Du|^{q}=f(u)\) in \(\Omega , 1<q<2,\) which satisfy the boundary condition \(u(x)\rightarrow \infty \) as \(x\rightarrow \partial \Omega \). These solutions are called large or blowup solutions. Near the boundary we give lower and upper bounds for the ratio \(\psi (u)/\delta \), where \(\psi (u) = \int _{u}^{\infty }1/\sqrt{2F}dt\), \(F'=f\), \(\delta =dist(x,\partial \Omega )\) or for the ratio \(u/\delta ^{(2-q)/(1-q)}\). When in particular the ratio \(f/F^{q/2}\)is regular at infinity, we find again known results (Bandle and Giarrusso, in Adv Diff Equ 1, 133–150, 1996; Giarrusso, in Comptes Rendus de l’Acad Sci 331, 777–782 2000).     

Lower order eigenvalues of a system of equations of the drifting Laplacian on the metric measure spaces

Let \(\Omega \) be a bounded domain with smooth boundary in an n-dimensional metric measure space \((\mathbb {R}^n, \langle ,\rangle , e^{-\phi }dv)\) and let \(\mathbf {u}=(u^1, \ldots , u^n)\) be a vector-valued function from \(\Omega \) to \(\mathbb {R}^n\). In this paper, we investigate the Dirichlet eigenvalue problem of a system of equations of the drifting Laplacian: \(\mathbb {L}_{\phi } \mathbf {u} + \alpha [ \nabla (\mathrm {div}\mathbf { u}) -\nabla \phi \mathrm {div} \mathbf {u}]= - \widetilde{\sigma } \mathbf {u}\), in \( \Omega \), and \(u|_{\partial \Omega }=0,\) where \(\mathbb {L}_{\phi } = \Delta - \nabla \phi \cdot \nabla \) is the drifting Laplacian and \(\alpha \) is a nonnegative constant. We establish some universal inequalities for lower order eigenvalues of this problem on the metric measure space \((\mathbb {R}^n, \langle ,\rangle , e^{-\phi }dv)\) and the Gaussian shrinking soliton \((\mathbb {R}^n, \langle ,\rangle _{\mathrm {can}}, e^{-\frac{|x|^2}{4}}dv, \frac{1}{2})\). Moreover, we give an estimate for the upper bound of the second eigenvalue of this problem in terms of its first eigenvalue on the gradient product Ricci soliton \((\Sigma \times \mathbb {R}, \langle ,\rangle , e^{-\frac{\kappa t^2}{2}}dv, \kappa )\), where \( \Sigma \) is an Einstein manifold with constant Ricci curvature \(\kappa \).     

Improved Moser–Trudinger inequality of Tintarev type in dimension <Emphasis Type="Italic">n</Emphasis> and the existence of its extremal functions

Let \(\Omega \) be a smooth bounded domain in \(\mathbb R^n\) with \(n\ge 2\), \(W^{1,n}_0(\Omega )\) be the usual Sobolev space on \(\Omega \) and define \(\lambda _1(\Omega ) = \inf \nolimits _{u\in W^{1,n}_0(\Omega )\setminus \{0\}}\frac{\int _\Omega |\nabla u|^n \mathrm{d}x}{\int _\Omega |u|^n \mathrm{d}x}\). Based on the blow-up analysis method, we shall establish the following improved Moser–Trudinger inequality of Tintarev type

$$\begin{aligned} \sup _{u\in W^{1,n}_0(\Omega ), \int _\Omega |\nabla u|^n \mathrm{{d}}x-\alpha \int _\Omega |u|^n \mathrm{{d}}x \le 1} \int _\Omega \exp (\alpha _{n} |u|^{\frac{n}{n-1}}) \mathrm{{d}}x < \infty , \end{aligned}$$

for any \(0 \le \alpha < \lambda _1(\Omega )\), where \(\alpha _{n} = n \omega _{n-1}^{\frac{1}{n-1}}\) with \(\omega _{n-1}\) being the surface area of the unit sphere in \(\mathbb R^n\). This inequality is stronger than the improved Moser–Trudinger inequality obtained by Adimurthi and Druet (Differ Equ 29:295–322, 2004) in dimension 2 and by Yang (J Funct Anal 239:100–126, 2006) in higher dimension and extends a result of Tintarev (J Funct Anal 266:55–66, 2014) in dimension 2 to higher dimension. We also prove that the supremum above is attained for any \(0< \alpha < \lambda _{1}(\Omega )\). (The case \(\alpha =0\) corresponding to the Moser–Trudinger inequality is well known.)

     

On the <Emphasis Type="Italic">p</Emphasis>-Laplacian with Robin boundary conditions and boundary trace theorems

Let \(\Omega \subset \mathbb {R}^\nu \), \(\nu \ge 2\), be a \(C^{1,1}\) domain whose boundary \(\partial \Omega \) is either compact or behaves suitably at infinity. For \(p\in (1,\infty )\) and \(\alpha >0\), define

$$\begin{aligned} \Lambda (\Omega ,p,\alpha ):=\inf _{\begin{array}{c} u\in W^{1,p}(\Omega )\\ u\not \equiv 0 \end{array}}\dfrac{\displaystyle \int _\Omega |\nabla u|^p \mathrm {d} x - \alpha \displaystyle \int _{\partial \Omega } |u|^p\mathrm {d}\sigma }{\displaystyle \int _\Omega |u|^p\mathrm {d} x}, \end{aligned}$$

where \(\mathrm {d}\sigma \) is the surface measure on \(\partial \Omega \). We show the asymptotics

$$\begin{aligned} \Lambda (\Omega ,p,\alpha )=-(p-1)\alpha ^{\frac{p}{p-1}} - (\nu -1)H_\mathrm {max}\, \alpha + o(\alpha ), \quad \alpha \rightarrow +\infty , \end{aligned}$$

where \(H_\mathrm {max}\) is the maximum mean curvature of \(\partial \Omega \). The asymptotic behavior of the associated minimizers is discussed as well. The estimate is then applied to the study of the best constant in a boundary trace theorem for expanding domains, to the norm estimate for extension operators and to related isoperimetric inequalities.

     

Multilicity results for some nonlinear ellitic roblems with asymtotically $${{\varvec{}}}$$-linear terms

Taking any \(p > 1\), we consider the asymptotically p-linear problem

$$\begin{aligned} \left\{ \begin{array}{ll} - {{\mathrm{div}}}(a(x,u,\nabla u)) + A_t(x,u,\nabla u)\ = \ \lambda ^\infty |u|^{p-2}u + g^\infty (x,u) &{}\quad \hbox {in}\;\Omega ,\\ u\ = \ 0 &{}\quad \hbox {on}\;\partial \Omega , \end{array} \right. \end{aligned}$$

where \(\Omega \) is a bounded domain in \(\mathbb R^N\), \(N\ge 2\), \(A(x,t,\xi )\) is a real function on \(\Omega \times \mathbb R\times \mathbb R^N\) which grows with power p with respect to \(\xi \) and has partial derivatives \(A_t(x,t,\xi ) = \frac{\partial A}{\partial t}(x,t,\xi )\), \(a(x,t,\xi ) = \nabla _\xi A(x,t,\xi )\). If \(A(x,t,\xi ) \rightarrow A^\infty (x,t)\) and \(\frac{g^\infty (x,t)}{|t|^{p-1}} \rightarrow 0\) as \(|t| \rightarrow +\infty \), suitable assumptions, variational methods and either the cohomological index theory or its related pseudo-index one, allow us to prove the existence of multiple nontrivial bounded solutions in the non-resonant case, i.e. if \(\lambda ^\infty \) is not an eigenvalue of the operator associated to \(\nabla _\xi A^\infty (x,\xi )\). In particular, while in [14] the model problem \(A(x,t,\xi ) = \mathcal{A}(x,t) |\xi |^p\) with \(p > N\) is studied, here our goal is twofold: extending such results not only to a more general family of functions \(A(x,t,\xi )\), but also to the more difficult case \(1 < p \le N\).

     

On Lipschitz continuity of solutions of hyperbolic Poisson’s equation

In this paper, we investigate solutions of the hyperbolic Poisson equation \(\Delta _{h}u(x)=\psi (x)\), where \(\psi \in L^{\infty }(\mathbb {B}^{n}, {\mathbb R}^n)\) and

$$\begin{aligned} \Delta _{h}u(x)= (1-|x|^2)^2\Delta u(x)+2(n-2)\left( 1-|x|^2\right) \sum _{i=1}^{n} x_{i} \frac{\partial u}{\partial x_{i}}(x) \end{aligned}$$

is the hyperbolic Laplace operator in the n-dimensional space \(\mathbb {R}^n\) for \(n\ge 2\). We show that if \(n\ge 3\) and \(u\in C^{2}(\mathbb {B}^{n},{\mathbb R}^n) \cap C(\overline{\mathbb {B}^{n}},{\mathbb R}^n )\) is a solution to the hyperbolic Poisson equation, then it has the representation \(u=P_{h}[\phi ]-G_{ h}[\psi ]\) provided that \(u\mid _{\mathbb {S}^{n-1}}=\phi \) and \(\int _{\mathbb {B}^{n}}(1-|x|^{2})^{n-1} |\psi (x)|\,d\tau (x)<\infty \). Here \(P_{h}\) and \(G_{h}\) denote Poisson and Green integrals with respect to \(\Delta _{h}\), respectively. Furthermore, we prove that functions of the form \(u=P_{h}[\phi ]-G_{h}[\psi ]\) are Lipschitz continuous.

     

Eigenvalue Counting Function for Robin Laplacians on Conical Domains

We study the discrete spectrum of the Robin Laplacian \(Q^{\Omega }_\alpha \) in \(L^2(\Omega )\), \(u\mapsto -\Delta u, \quad D_n u=\alpha u \text { on }\partial \Omega \), where \(D_n\) is the outer unit normal derivative and \(\Omega \subset {\mathbb {R}}^{3}\) is a conical domain with a regular cross-section \(\Theta \subset {\mathbb {S}}^2\), n is the outer unit normal, and \(\alpha >0\) is a fixed constant. It is known from previous papers that the bottom of the essential spectrum of \(Q^{\Omega }_\alpha \) is \(-\alpha ^2\) and that the finiteness of the discrete spectrum depends on the geometry of the cross-section. We show that the accumulation of the discrete spectrum of \(Q^\Omega _\alpha \) is determined by the discrete spectrum of an effective Hamiltonian defined on the boundary and far from the origin. By studying this model operator, we prove that the number of eigenvalues of \(Q^{\Omega }_\alpha \) in \((-\infty ,-\alpha ^2-\lambda )\), with \(\lambda >0\), behaves for \(\lambda \rightarrow 0\) as

$$\begin{aligned} \dfrac{\alpha ^2}{8\pi \lambda } \int _{\partial \Theta } \kappa _+(s)^2\mathrm {d}s +o\left( \frac{1}{\lambda }\right) , \end{aligned}$$

where \(\kappa _+\) is the positive part of the geodesic curvature of the cross-section boundary.

     

On a Fractional $$ p \& q$$ Laplacian Problem with Critical Sobolev–Hardy Exponents

We consider the following fractional \( p \& q\) Laplacian problem with critical Sobolev–Hardy exponents

$$\begin{aligned} \left\{ \begin{array}{ll} (-\Delta )^{s}_{p} u + (-\Delta )^{s}_{q} u = \frac{|u|^{p^{*}_{s}(\alpha )-2}u}{|x|^{\alpha }}+ \lambda f(x, u) &{} \text{ in } \Omega \\ u=0 &{} \text{ in } \mathbb {R}^{N}{\setminus } \Omega , \end{array} \right. \end{aligned}$$

where \(0<s<1\), \(1\le q<p<\frac{N}{s}\), \((-\Delta )^{s}_{r}\), with \(r\in \{p,q\}\), is the fractional r-Laplacian operator, \(\lambda \) is a positive parameter, \(\Omega \subset \mathbb {R}^{N}\) is an open bounded domain with smooth boundary, \(0\le \alpha <sp\), and \(p^{*}_{s}(\alpha )=\frac{p(N-\alpha )}{N-sp}\) is the so-called Hardy–Sobolev critical exponent. Using concentration-compactness principle and the mountain pass lemma due to Kajikiya [23], we show the existence of infinitely many solutions which tend to be zero provided that \(\lambda \) belongs to a suitable range.

     

Mean curvature bounds and eigenvalues of Robin Laplacians

We consider the Laplacian with attractive Robin boundary conditions,

$$\begin{aligned} Q^\Omega _\alpha u=-\Delta u, \quad \dfrac{\partial u}{\partial n}=\alpha u \text { on } \partial \Omega , \end{aligned}$$

in a class of bounded smooth domains \(\Omega \in \mathbb {R}^\nu \); here \(n\) is the outward unit normal and \(\alpha >0\) is a constant. We show that for each \(j\in \mathbb {N}\) and \(\alpha \rightarrow +\infty \), the \(j\)th eigenvalue \(E_j(Q^\Omega _\alpha )\) has the asymptotics

$$\begin{aligned} E_j(Q^\Omega _\alpha )=-\alpha ^2 -(\nu -1)H_\mathrm {max}(\Omega )\,\alpha +{\mathcal O}(\alpha ^{2/3}), \end{aligned}$$

where \(H_\mathrm {max}(\Omega )\) is the maximum mean curvature at \(\partial \Omega \). The discussion of the reverse Faber-Krahn inequality gives rise to a new geometric problem concerning the minimization of \(H_\mathrm {max}\). In particular, we show that the ball is the strict minimizer of \(H_\mathrm {max}\) among the smooth star-shaped domains of a given volume, which leads to the following result: if \(B\) is a ball and \(\Omega \) is any other star-shaped smooth domain of the same volume, then for any fixed \(j\in \mathbb {N}\) we have \(E_j(Q^B_\alpha )>E_j(Q^\Omega _\alpha )\) for large \(\alpha \). An open question concerning a larger class of domains is formulated.

     

Existence of Solution for a Singular Elliptic Equation of Kirchhoff Type

We study both the existence and uniqueness of nonnegative solution to a singular elliptic problem of Kirchhoff type, whose model is:

$$\begin{aligned} {\left\{ \begin{array}{ll} -B\left( \dfrac{1}{2}\displaystyle \int _\Omega |\nabla u|^2\mathrm {d}x\right) \Delta u=\dfrac{h(x)}{u^\gamma }, &{}\quad x\in \Omega ,\\ u>0, &{}\quad x\in \Omega ,\\ u=0, &{}\quad x\in \partial \Omega , \end{array}\right. } \end{aligned}$$

where \(\Omega \subset \mathbb {R}^n(n\ge 1)\) is a smooth bounded domain, \(\gamma >1\), \(h\in L^1(\Omega )\) is positive (i.e., \(h(x)>0\) a.e. in \(\Omega \)), \(B : \mathbb {R}^+\rightarrow \mathbb {R}^+\) is a \(C^1\)-continuous function with positive lower bound. A necessary and sufficient condition will be given for the existence of weak solution of the general nonlocal singular elliptic with strong singularity. In addition, we prove that the solution is unique under some suitable conditions.

     

Global $$C^{1,\al pha }$$ regularity and e xistence of multi ple solutions for singular p( x)-La placian equations

In this paper we study the following singular p(x)-Laplacian problem

$$\begin{aligned} \left\{ \begin{array}{l@{\quad }l} - \text{ div } \left( |\nabla u|^{p(x)-2} \nabla u\right) =\frac{ \lambda }{u^{\beta (x)}}+u^{q(x)}, &{} \text{ in }\quad \Omega , \\ u>0, &{} \text{ in }\quad \Omega , \\ u=0, &{} \text{ on }\quad \partial \Omega , \end{array}\right. \end{aligned}$$

where \(\Omega \) is a bounded domain in \(\mathbb {R}^N\), \(N\ge 2\), with smooth boundary \(\partial \Omega \), \(\beta \in C^1(\bar{\Omega })\) with \( 0< \beta (x) <1\), \(p\in C^1(\bar{\Omega })\), \(q \in C(\bar{\Omega })\) with \(p(x)>1\), \(p(x)< q(x) +1 <p^*(x)\) for \(x \in \bar{\Omega }\), where \( p^*(x)= \frac{Np(x)}{N-p(x)} \) for \(p(x) <N\) and \( p^*(x)= \infty \) for \( p(x) \ge N\). We establish \(C^{1,\alpha }\) regularity of weak solutions of the problem and strong comparison principle. Based on these two results, we prove the existence of multiple (at least two) positive solutions for a certain range of \(\lambda \).

     

Stability of Standing Waves for a Generalized Choquard Equation with Potential

We are going to study the standing waves for the generalized Choquard equation with potential:

$$ -i\partial_{t} u-\Delta u+V(x)u=\bigl(|x|^{-1}\ast |u|^{p}\bigr)|u|^{p-2}u,\quad \hbox{in } \mathbb{R}\times \mathbb{R}^{3}, $$

where \(V(x)\) is a real function, \(p>2\) is close to 2 and ? standards for the convolution in \(\mathbb{R}^{3}\). The stability of the standing waves \(u(x)=e^{i\omega t}\varphi (x)\) is investigated under suitable assumptions on the potential and the frequency \(\omega \).

     

Existence results for a singular quasilinear elliptic equation

Let \(\Omega \subset \mathbb R^N\) be a bounded domain with smooth boundary. Existence of a positive solution to the quasilinear equation

$$\begin{aligned} -\text {div}\left[ \left( a(x)+|u|^\theta \right) \nabla u\right] +\frac{\theta }{2}|u|^{\theta -2}u|\nabla u|^2=|u|^{p-2}u \quad \text {in}\ \Omega \end{aligned}$$

with zero Dirichlet boundary condition is proved. Here \(\theta >0\) and a(x) is a measurable function satisfying \(0<\alpha \le a(x)\le \beta \). The equation involves singularity when \(0<\theta \le 1\). As a main novelty with respect to corresponding results in the literature, we only assume \(\theta +2<p<\frac{2^*}{2}(\theta +2)\). The proof relies on a perturbation method and a critical point theory for E-differentiable functionals.

     

Hölder continuity for continuous solutions of the singular minimal surface equation with arbitrary zero set

In this paper we prove the following theorem: Let \(\Omega \subset \mathbb {R}^{n}\) be a bounded open set, \(\psi \in C_{c}^{2}(\mathbb {R}^{n})\), \(\psi > 0\) on \(\partial \Omega \), be given boundary values and u a nonnegative solution to the problem

$$\begin{aligned}&u \in C^{0}(\overline{\Omega }) \cap C^{2}(\{u> 0\}) \\&u = \psi \quad \text { on } \; \partial \Omega \\&{\text {div}} \left( \frac{Du}{\sqrt{1 + |Du|^{2}}}\right) = \frac{\alpha }{u \sqrt{1 + |Du|^{2}}} \quad \text { in } \; \{u > 0\} \end{aligned}$$

where \(\alpha > 0\) is a given constant. Then \(u \in C^{0, \frac{1}{2}} (\overline{\Omega })\). Furthermore we prove strict mean convexity of the free boundary \(\partial \{u = 0\}\) provided \(\partial \{u = 0\}\) is assumed to be of class \(C^{2}\) and \(\alpha \ge 1\).

     

Book reviews

We consider the following singularly perturbed nonlocal elliptic problem

$$\begin{aligned} -\left( \varepsilon ^{2}a+\varepsilon b\displaystyle \int _{\mathbb {R}^{3}}|\nabla u|^{2}dx\right) \Delta u+V(x)u=\displaystyle \varepsilon ^{\alpha -3}(W_{\alpha }(x)*|u|^{p})|u|^{p-2}u, \quad x\in \mathbb {R}^{3}, \end{aligned}$$

where \(\varepsilon >0\) is a parameter, \(a>0,b\ge 0\) are constants, \(\alpha \in (0,3)\), \(p\in [2, 6-\alpha )\), \(W_{\alpha }(x)\) is a convolution kernel and V(x) is an external potential satisfying some conditions. By using variational methods, we establish the existence and concentration of positive ground state solutions for the above equation.

     

Regularity for Inhomogeneous Dirichlet Problems of Some Schrödinger Equations on Domains

Let \(n\ge 3, \Omega \) be a bounded, simply connected and semiconvex domain in \({\mathbb {R}}^n\) and \(L_{\Omega }:=-\Delta +V\) a Schrödinger operator on \(L^2 (\Omega )\) with the Dirichlet boundary condition, where \(\Delta \) denotes the Laplace operator and the potential \(0\le V\) belongs to the reverse Hölder class \(RH_{q_0}({\mathbb {R}}^n)\) for some \(q_0\in (\max \{n/2,2\},\infty ]\). Assume that the growth function \(\varphi :\,{\mathbb {R}}^n\times [0,\infty ) \rightarrow [0,\infty )\) satisfies that \(\varphi (x,\cdot )\) is an Orlicz function and \(\varphi (\cdot ,t)\in {\mathbb {A}}_{\infty }({\mathbb {R}}^n)\) (the class of uniformly Muckenhoupt weights). Let \(H_{\varphi ,\,L_{{\mathbb {R}}^n},\,r}(\Omega )\) be the Musielak–Orlicz–Hardy space whose elements are restrictions of elements of the Musielak–Orlicz–Hardy space, associated with \(L_{{\mathbb {R}}^n}:=-\Delta +V\) on \({\mathbb {R}}^n\), to \(\Omega \). In this article, the authors show that the operators \(VL^{-1}_\Omega \) and \(\nabla ^2L^{-1}_\Omega \) are bounded from \(L^1(\Omega )\) to weak-\(L^1(\Omega )\), from \(L^p(\Omega )\) to itself, with \(p\in (1,2]\), and also from \(H_{\varphi ,\,L_{{\mathbb {R}}^n},\,r}(\Omega )\) to the Musielak–Orlicz space \(L^\varphi (\Omega )\) or to \(H_{\varphi ,\,L_{{\mathbb {R}}^n},\,r}(\Omega )\) itself. As applications, the boundedness of \(\nabla ^2{\mathbb {G}}_D\) on \(L^p(\Omega )\), with \(p\in (1,2]\), and from \(H_{\varphi ,\,L_{{\mathbb {R}}^n},\,r}(\Omega )\) to \(L^\varphi (\Omega )\) or to \(H_{\varphi ,\,L_{{\mathbb {R}}^n},\,r}(\Omega )\) itself is obtained, where \({\mathbb {G}}_D\) denotes the Dirichlet Green operator associated with \(L_\Omega \). All these results are new even for the Hardy space \(H^1_{L_{{\mathbb {R}}^n},\,r}(\Omega )\), which is just \(H_{\varphi ,\,L_{{\mathbb {R}}^n},\,r}(\Omega )\) with \(\varphi (x,t):=t\) for all \(x\in {\mathbb {R}}^n\) and \(t\in [0,\infty )\).     

Semilinear elliptic equations with Hardy potential and subcritical source term

Let \(\Omega \) be a smooth bounded domain in \({\mathbb {R}}^N\) (\(N>2\)) and \(\delta (x):=\text {dist}\,(x,\partial \Omega )\). Assume \(\mu \in {\mathbb {R}}_+, \nu \) is a nonnegative finite measure on \(\partial \Omega \) and \(g \in C(\Omega \times {\mathbb {R}}_+)\). We study positive solutions of

$$\begin{aligned} -\Delta u - \frac{\mu }{\delta ^2} u = g(x,u) \text { in } \Omega , \qquad \text {tr}^*(u)=\nu . \end{aligned}$$

(P)

Here \(\text {tr}^*(u)\) denotes the normalized boundary trace of u which was recently introduced by Marcus and Nguyen (Ann Inst H Poincaré Anal Non Linéaire, 34, 69–88, 2017). We focus on the case \(0<\mu < C_H(\Omega )\) (the Hardy constant for \(\Omega \)) and provide qualitative properties of positive solutions of (P). When \(g(x,u)=u^q\) with \(q>0\), we prove that there is a critical value \(q^*\) (depending only on \(N, \mu \)) for (P) in the sense that if \(q<q^*\) then (P) possesses a solution under a smallness assumption on \(\nu \), but if \(q \ge q^*\) this problem admits no solution with isolated boundary singularity. Existence result is then extended to a more general setting where g is subcritical [see (1.28)]. We also investigate the case where g is linear or sublinear and give an existence result for (P).