Multiple positive solutions of singular and critical elliptic problem in {\mathbb{R}^2} with discontinuous nonlinearities

Let Ω be a bounded domain in ${\mathbb{R}^2}$ with smooth boundary. We consider the following singular and critical elliptic problem with discontinuous nonlinearity: $$(P_\lambda)\left \{\begin{array}{ll} - \Delta u = \lambda \left(\frac{m(x, u) e^{\alpha{u}^2}}{|x|^{\beta}} + u^{q}g(u - a)\right),\quad{u} > 0 \quad {\rm in} \quad \Omega\\u \quad \quad = 0\quad {\rm on} \quad \partial \Omega \end{array}\right.$$ where ${0\leq q < 1 ,0< \alpha\leq4\pi}$ and ${\beta \in [0, 2)}$ such that ${\frac{\beta}{2} + \frac{\alpha}{4\pi} \leq 1}$ and ${{g(t - a) = \left\{\begin{array}{ll}1, t \leq a\\ 0, t > a.\end{array}\right.}}$ Under the suitable assumptions on m(x, t) we show the existence and multiplicity of solutions for maximal interval for λ.     

Stationary solutions of the nonlinear Schrödinger equation with fast-decay potentials concentrating around local maxima

We study positive bound states for the equation ${- \varepsilon^2 \Delta u + Vu = u^p, \quad {\rm in} \quad \mathbb{R}^N}$ , where ${\varepsilon > 0}$ is a real parameter, ${\frac{N}{N-2} < p < \frac{N+2}{N-2}}$ and V is a nonnegative potential. Using purely variational techniques, we find solutions which concentrate at local maxima of the potential V without any restriction on the potential.     

Regularity lifting of weak solutions for nonlinear sub-Laplace equations on homogeneous groups

Let G be a homogeneous group, and let X ₁, X ₂, … , X _m be left invariant real vector fields being homogeneous of degree one on G. We consider the following Dirichlet boundary value problem of the sub-Laplace equation involving the critical exponent and singular term: $$\left\{\begin{array}{ll}-\sum_{j=1}^{m}X_j^2u(x)-\frac{a}{\|x\|^\nu}u(x)=u^{\frac{Q+2}{Q-2}}(x), x\in\Omega,\\ u(x)=0, \quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\,\,\,\, x\in \partial\Omega,\end{array}\right.$$ where ${\Omega\subset G}$ is a bounded domain with smooth boundary and ${\mathbf{0}\in\Omega}$ , Q is the homogeneous dimension of G, ${a\in \mathbb{R},\ \nu <2 }$ . We boost u to ${L^p(\Omega)}$ for any ${1\leq p < \infty}$ if ${u\in S^{1,2}_0(\Omega)}$ is a weak solution of the problem above.     

On averages of Fourier coefficients of Maass cusp forms

Let φ be a primitive Maass cusp form and t _φ (n) be its nth Fourier coefficient at the cusp infinity. In this short note, we are interested in the estimation of the sums ${\sum_{n \leq x}t_{\varphi}(n)}$ and ${\sum_{n \leq x}t_{\varphi}(n^2)}$ . We are able to improve the previous results by showing that for any ${\varepsilon > 0}$ $$\sum_{n \leq x}t_{\varphi}(n) \ll\, _{\varphi, \varepsilon} x^{\frac{1027}{2827} + \varepsilon} \quad {and}\quad\sum_{n \leq x}t_{\varphi}(n^2) \ll\,_{\varphi, \varepsilon} x^{\frac{489}{861} + \varepsilon}.$$      

A sharp equivalence between H ∞ functional calculus and square function estimates

Let (T _t ) _{t?≥ 0} be a bounded analytic semigroup on L ^p (Ω), with 1?<?p?<?∞. Let ?A denote its infinitesimal generator. It is known that if A and A ^* both satisfy square function estimates ${\bigl\|\bigl(\int_{0}^{\infty} \vert A^{\frac{1}{2}} T_t(x)\vert^2 {\rm d}t \bigr)^{\frac{1}{2}}\bigr\|_{L^p} \lesssim \|x\|_{L^p}}$ and ${\bigl\|\bigl(\int_{0}^{\infty} \vert A^{*\frac{1}{2}} T_t^*(y) \vert^2 {\rm d}t \bigr)^{\frac{1}{2}}\bigr\|_{L^{p^\prime}} \lesssim \|y\|_{L^{p^\prime}}}$ for ${x\in L^p(\Omega)}$ and ${y\in L^{p^\prime}(\Omega)}$ , then A admits a bounded ${H^{\infty}(\Sigma_\theta)}$ functional calculus for any ${\theta>\frac{\pi}{2}}$ . We show that this actually holds true for some ${\theta<\frac{\pi}{2}}$ .     

Infinitely many solutions for an elliptic problem involving critical Sobolev growth and Hardy potential

In this paper, we will prove the existence of infinitely many solutions for the following elliptic problem with critical Sobolev growth and a Hardy potential: $$-\Delta u-\frac{\mu}{|x|^2}u = |u|^{2^{\ast}-2}u+a u\quad {\rm in}\;\Omega,\quad u=0 \quad {\rm on}\; \partial\Omega,\qquad (*)$$ under the assumptions that N ≥ 7, ${\mu\in \left[0,\frac{(N-2)^2}4-4\right)}$ and a > 0, where ${2^{\ast}=\frac{2N}{N-2}}$ , and Ω is an open bounded domain in ${\mathbb{R}^N}$ which contains the origin. To achieve this goal, we consider the following perturbed problem of (*), which is of subcritical growth, $$-\Delta u-\frac{\mu}{|x|^2}u = |u|^{2^{\ast}-2-\varepsilon_n}u+au \quad {\rm in}\,\Omega, \quad u=0 \quad {\rm on}\;\partial\Omega,\qquad(\ast\ast)_n$$ where ${\varepsilon_{n} > 0}$ is small and ${\varepsilon_n \to 0}$ as n → + ∞. By the critical point theory for the even functionals, for each fixed ${\varepsilon_{n} > 0}$ small, (**) _n has a sequence of solutions ${u_{k,\varepsilon_{n}} \in H^{1}_{0}(\Omega)}$ . We obtain the existence of infinitely many solutions for (*) by showing that as n → ∞, ${u_{k,\varepsilon_{n}}}$ converges strongly in ${H^{1}_{0}(\Omega)}$ to u _k , which must be a solution of (*). Such a convergence is obtained by applying a local Pohozaev identity to exclude the possibility of the concentration of ${\{u_{k,\varepsilon_n}\}}$ .     

Infinitely many solutions for the prescribed curvature problem of polyharmonic operator

We consider the following prescribed curvature problem for polyharmonic operator: $$\left\{\begin{array}{llll} D_{m} u = \tilde{K}(y)|u|^{m^*-2}u\; {\rm in}\; \mathbb{S}^N\\ u \quad\; >0\qquad\quad\quad\quad\quad{\rm on}\; \mathbb{S}^N\\ u \quad\; \in H^{m}(\mathbb{S}^N), \end{array} \right.$$ where ${m^*=\frac{2N}{N-2m}, N\geq 2m+1,m \in \mathbb{N}_{+}, \tilde{K}}$ is positive and rationally symmetric, ${\mathbb{S}^N}$ is the unit sphere with the induced Riemannian metric ${g=g_{\mathbb{S}^N},}$ and D _m is the elliptic differential operator of 2m order given by $$\begin{array}{lll}D_m={\prod\limits_{k=1}^m}{\left(-\Delta_g+\frac{1}{4}(N-2k)(N+2k-2)\right)}\end{array}$$ where Δ _g is the Laplace-Beltrami operator on ${\mathbb{S}^N}$ . We will show that problem (P) has infinitely many non-radial positive solutions, whose energy can be arbitrary large.     

The geometric Neumann problem for the Liouville equation

Let Ω denote the upper half-plane ${\mathbb{R}_+^2}$ or the upper half-disk ${D_{\varepsilon}^+\subset \mathbb{R}_+^2}$ of center 0 and radius ${\varepsilon}$ . In this paper we classify the solutions ${v\in\;C^2(\overline{\Omega}\setminus\{0\})}$ to the Neumann problem $$\left\{\begin{array}{lll}{\Delta v+2 Ke^v=0\quad {\rm in}\,\Omega\subseteq \mathbb{R}^2_+=\{(s, t)\in \mathbb{R}^2: t >0 \},}\\ {\frac{\partial v}{\partial t}=c_1e^{v/2}\quad\quad\quad{\rm on}\,\partial\Omega\cap\{s >0 \},}\\ {\frac{\partial v}{\partial t}=c_2e^{v/2}\quad\quad\quad{\rm on}\,\partial\Omega\cap\{s <0 \},}\end{array}\right.$$ where ${K, c_1, c_2 \in \mathbb{R}}$ , with the finite energy condition ${\int_{\Omega} e^v < \infty}$ As a result, we classify the conformal Riemannian metrics of constant curvature and finite area on a half-plane that have a finite number of boundary singularities, not assumed a priori to be conical, and constant geodesic curvature along each boundary arc.     

Isolated boundary singularities of semilinear elliptic equations

Given a smooth domain ${\Omega\subset\mathbb{R}^N}$ such that ${0 \in \partial\Omega}$ and given a nonnegative smooth function ?? on ???, we study the behavior near 0 of positive solutions of ???u?=?u ^q in ?? such that u =? ?? on ???\{0}. We prove that if ${\frac{N+1}{N-1} < q < \frac{N+2}{N-2}}$ , then ${u(x)\leq C |x|^{-\frac{2}{q-1}}}$ and we compute the limit of ${|x|^{\frac{2}{q-1}} u(x)}$ as x ?? 0. We also investigate the case ${q= \frac{N+1}{N-1}}$ . The proofs rely on the existence and uniqueness of solutions of related equations on spherical domains.     

Existence and Multiplicity Results for Some Elliptic Systems in Unbounded Cylinders

We study the following nonlinear elliptic system of Lane–Emden type $$\left\{\begin{array}{ll} -\Delta u = {\rm sgn}(v) |v| ^{p-1} \qquad \qquad \qquad \; {\rm in} \; \Omega , \\ -\Delta v = - \lambda {\rm sgn} (u)|u| \frac{1}{p-1} + f(x, u)\; \; {\rm in}\; \Omega , \\ u = v = 0 \qquad \qquad \qquad \quad \quad \;\;\;\;\; {\rm on}\; \partial \Omega , \end{array}\right.$$ where ${\lambda \in \mathbb{R}}$ . If ${\lambda \geq 0}$ and ${\Omega}$ is an unbounded cylinder, i.e., ${\Omega = \tilde \Omega \times \mathbb{R}^{N-m} \subset \mathbb{R}^{N}}$ , ${N - m \geq 2, m \geq 1}$ , existence and multiplicity results are proved by means of the Principle of Symmetric Criticality and some compact imbeddings in partially spherically symmetric spaces. We are able to state existence and multiplicity results also if ${\lambda \in \mathbb{R}}$ and ${\Omega}$ is a bounded domain in ${\mathbb{R}^{N}, N \geq 3}$ . In particular, a good finite dimensional decomposition of the Banach space in which we work is given.     

Multibump bound states for quasilinear Schrödinger systems with critical frequency

In this paper, we are concerned with the multibump solutions for the following quasilinear Schrödinger system in ${\mathbb{R}^N}$ : $$\left\{\begin{array}{ll}-\Delta{u} + \lambda{a(x)u} - \frac{1}{2}(\Delta|u|^2)u = \frac{2\alpha}{\alpha + \beta}|u|^{\alpha-2}|\upsilon|^\beta u, \\-\Delta{\upsilon} + \lambda{b(x)\upsilon} - \frac{1}{2}(\Delta|\upsilon|^2)\upsilon = \frac{2\beta}{\alpha + \beta}|u|^\alpha|\upsilon|^{\beta-2} \upsilon, \\u(x) \rightarrow 0, \upsilon(x) \rightarrow 0 \quad as|x| \rightarrow \infty,\end{array}\right.$$ where λ > 0 is a parameter, α, β > 2 satisfying α + β < 2 · 2*, here ${2^{*} = \frac{2N}{N-2}}$ is the critical Sobolev exponent for ${N \geq 3}$ and a(x), b(x) are nonnegative potentials. Using variational methods, we prove that if the zero sets of a(x) and b(x) have several common isolated connected components ${\Omega_{1}, . . . ,\Omega_{k}}$ such that the interior of ${\Omega_{i} (i = 1, 2, . . . , k)}$ is not empty and ${\partial\Omega_{i} (i = 1, 2, . . . , k)}$ is smooth, then for λ sufficiently large, the system admits, for any nonempty subset ${J \subset \{1, 2, . . . , k\}}$ , a solution which is trapped in a neighborhood of ${\cup_{j\epsilon{J}} \Omega_{j}}$ .     

Quenched to annealed transition in the parabolic Anderson problem

We study limit behavior for sums of the form $\frac{1}{|\Lambda_{L|}}\sum_{x\in \Lambda_{L}}u(t,x),$ where the field $\Lambda_L=\left\{x\in {\bf{Z^d}}:|x|\le L\right\}$ is composed of solutions of the parabolic Anderson equation $$u(t,x) = 1 + \kappa \mathop{\int}_{0}^{t} \Delta u(s,x){\rm d}s + \mathop{\int}_{0}^{t}u(s,x)\partial B_{x}(s). $$ The index set is a box in Z ^d , namely $\Lambda_{L} = \left\{x\in {\bf Z}^{\bf d} : |x| \leq L\right\}$ and L = L(t) is a nondecreasing function $L : [0,\infty)\rightarrow {\bf R}^{+}. $ We identify two critical parameters $\eta(1) < \eta(2)$ such that for $\gamma > \eta(1)$ and L(t) = e^{γ t} , the sums $\frac{1}{|\Lambda_L|}\sum_{x\in \Lambda_L}u(t,x)$ satisfy a law of large numbers, or put another way, they exhibit annealed behavior. For $\gamma > \eta(2)$ and L(t) = e^{γ t} , one has $\sum_{x\in \Lambda_L}u(t,x)$ when properly normalized and centered satisfies a central limit theorem. For subexponential scales, that is when $\lim_{t \rightarrow \infty} \frac{1}{t}\ln L(t) = 0,$ quenched asymptotics occur. That means $\lim_{t\rightarrow \infty}\frac{1}{t}\ln\left (\frac{1}{|\Lambda_L|}\sum_{x\in \Lambda_L}u(t,x)\right) = \gamma(\kappa),$ where $\gamma(\kappa)$ is the almost sure Lyapunov exponent, i.e. $\lim_{t\rightarrow \infty}\frac{1}{t}\ln u(t,x)= \gamma(\kappa).$ We also examine the behavior of $\frac{1}{|\Lambda_L|}\sum_{x\in \Lambda_L}u(t,x)$ for L = e ^{γ t} with γ in the transition range $(0,\eta(1))$      

1D Schrödinger Operators with Short Range Interactions: Two-Scale Regularization of Distributional Potentials

For real ${L_\infty(\mathbb{R})}$ -functions ${\Phi}$ and ${\Psi}$ of compact support, we prove the norm resolvent convergence, as ${\varepsilon}$ and ${\nu}$ tend to 0, of a family ${S_{\varepsilon \nu}}$ of one-dimensional Schrödinger operators on the line of the form $$S_{\varepsilon \nu} = -\frac{d^2}{dx^2} + \frac{\alpha}{\varepsilon^2} \Phi \left( \frac{x}{\varepsilon} \right) + \frac{\beta}{\nu} \Psi \left(\frac{x}{\nu} \right),$$ provided the ratio ${\nu/\varepsilon}$ has a finite or infinite limit. The limit operator S ₀ depends on the shape of ${\Phi}$ and ${\Psi}$ as well as on the limit of ratio ${\nu/\varepsilon}$ . If the potential ${\alpha\Phi}$ possesses a zero-energy resonance, then S ₀ describes a non trivial point interaction at the origin. Otherwise S ₀ is the direct sum of the Dirichlet half-line Schrödinger operators.     

A decomposition of the universal embedding space for the near polygon {{\mathbb H}_n}

Let ${{\mathbb H}_n, n \geq 1}$ , be the near 2n-gon defined on the 1-factors of the complete graph on 2n?+?2 vertices, and let e denote the absolutely universal embedding of ${{\mathbb H}_n}$ into PG(W), where W is a ${\frac{1}{n+2} \left(\begin{array}{c}2n+2 \\ n+1\end{array}\right)}$ -dimensional vector space over the field ${{\mathbb F}_2}$ with two elements. For every point z of ${{\mathbb H}_n}$ and every ${i \in {\mathbb N}}$ , let Δ _i (z) denote the set of points of ${{\mathbb H}_n}$ at distance i from z. We show that for every pair {x, y} of mutually opposite points of ${{\mathbb H}_n, W}$ can be written as a direct sum ${W_0 \oplus W_1 \oplus \cdots \oplus W_n}$ such that the following four properties hold for every ${i \in \{0,\ldots,n \}}$ : (1) ${\langle e(\Delta_i(x) \cap \Delta_{n-i}(y)) \rangle = {\rm PG}(W_i)}$ ; (2) ${\left\langle e \left( \bigcup_{j \leq i} \Delta_j(x) \right) \right\rangle = {\rm PG}(W_0 \oplus W_1 \oplus \cdots \oplus W_i)}$ ; (3) ${\left\langle e \left( \bigcup_{j \leq i} \Delta_j(y) \right) \right\rangle = {\rm PG}(W_{n-i}\oplus W_{n-i+1} \oplus \cdots \oplus W_n)}$ ; (4) ${\dim(W_i) = |\Delta_i(x) \cap \Delta_{n-i}(y)| = \left(\begin{array}{c}n \\ i\end{array}\right)^2 - \left(\begin{array}{c}n \\ i-1\end{array}\right) \cdot \left(\begin{array}{c}n \\ i+1\end{array}\right)}$ .     

Boundedness of Generalized Riesz Transforms on Orlicz–Hardy Spaces Associated to Operators

Let ${\Phi}$ be a continuous, strictly increasing and concave function on (0, ∞) of critical lower type index ${p_\Phi^- \in(0,\,1]}$ . Let L be an injective operator of type ω having a bounded H _∞ functional calculus and satisfying the k-Davies–Gaffney estimates with ${k \in {\mathbb Z}_+}$ . In this paper, the authors first introduce an Orlicz–Hardy space ${H^{\Phi}_{L}(\mathbb{R}^n)}$ in terms of the non-tangential L-adapted square function and then establish its molecular characterization. As applications, the authors prove that the generalized Riesz transform ${D_{\gamma}L^{-\delta/(2k)}}$ is bounded from the Orlicz–Hardy space ${H^{\Phi}_{L}(\mathbb{R}^n)}$ to the Orlicz space ${L^{\widetilde{\Phi}}(\mathbb{R}^n)}$ when ${p_\Phi^- \in (0, \frac{n}{n+ \delta - \gamma}]}$ , ${0 < \gamma \le \delta < \infty}$ and ${\delta- \gamma < n (\frac{1}{p_-(L)}-\frac{1}{p_+(L)})}$ , or from ${H^{\Phi}_{L}(\mathbb{R}^n)}$ to the Orlicz–Hardy space ${H^{\widetilde \Phi}(\mathbb{R}^n)}$ when ${p_\Phi^-\in (\frac{n}{n + \delta+ \lfloor \gamma \rfloor- \gamma},\,\frac{n}{n+ \delta- \gamma}]}$ , ${1\le \gamma \le \delta < \infty}$ and ${\delta- \gamma < n (\frac{1}{p_-(L)}-\frac{1}{p_+(L)})}$ , or from ${H^{\Phi}_{L}(\mathbb{R}^n)}$ to the weak Orlicz–Hardy space ${WH^\Phi(\mathbb{R}^n)}$ when ${\gamma = \delta}$ and ${p_\Phi=n/(n + \lfloor \gamma \rfloor)}$ or ${p_\Phi^-=n/(n + \lfloor \gamma \rfloor)}$ with ${p_\Phi^-}$ attainable, where ${\widetilde{\Phi}}$ is an Orlicz function whose inverse function ${\widetilde{\Phi}^{-1}}$ is defined by ${\widetilde{\Phi}^{-1}(t):=\Phi^{-1}(t)t^{\frac{1}{n}(\gamma- \delta)}}$ for all ${t \in (0,\,\infty)}$ , ${p_\Phi}$ denotes the strictly critical lower type index of ${\Phi}$ , ${\lfloor \gamma \rfloor}$ the maximal integer not more than ${\gamma}$ and ${(p_-(L),\,p_+(L))}$ the range of exponents ${p \in[1,\, \infty]}$ for which the semigroup ${\{e^{-tL}\}_{t >0 }}$ is bounded on ${L^p(\mathbb{R}^n)}$ .     

Invariant submanifolds of Sasakian space forms

In the present study, we consider isometric immersions ${f : M \rightarrow \tilde{M}(c)}$ of (2n + 1)-dimensional invariant submanifold M ²ⁿ⁺¹ of (2m + 1) dimensional Sasakian space form ${\tilde{M}^{2m+1}}$ of constant ${ \varphi}$ -sectional curvature c. We have shown that if f satisfies the curvature condition ${\overset{\_}{R}(X, Y) \cdot \sigma =Q(g, \sigma)}$ then either M ²ⁿ⁺¹ is totally geodesic, or ${||\sigma||^{2}=\frac{1}{3}(2c+n(c+1)),}$ or ${||\sigma||^{2}(x) > \frac{1}{3}(2c+n(c+1)}$ at some point x of M ²ⁿ⁺¹. We also prove that ${\overset{\_ }{R}(X, Y)\cdot \sigma = \frac{1}{2n}Q(S, \sigma)}$ then either M ²ⁿ⁺¹ is totally geodesic, or ${||\sigma||^{2}=-\frac{2}{3}(\frac{1}{2n}\tau -\frac{1}{2}(n+2)(c+3)+3)}$ , or ${||\sigma||^{2}(x) > -\frac{2}{3}(\frac{1}{2n} \tau (x)-\frac{1}{2} (n+2)(c+3)+3)}$ at some point x of M ²ⁿ⁺¹.     

A weighted Hardy inequality and nonexistence of positive solutions

In this article, we prove that the following weighted Hardy inequality $$\begin{array}{ll}\big(\frac{|{d-p}|}{p}\big)^{p}\, \int\limits_{\Omega}\, \frac{|{u}|^{p}}{|{x}|^{p}}\;d\mu \\ \quad \quad \le \int\limits_{\Omega}\,|{\nabla u}|^{p}\;d\mu+ \big(\frac{|{d-p}|}{p}\big)^{p-1}\,\textrm{sgn}(d-p)\,\int\limits_{\Omega}|{u}|^{p}\,\frac{(x^{t}Ax)^{p/2}}{|{x}|^{p}}\; d\mu \quad \quad \quad (1) \end{array}$$ holds with optimal Hardy constant ${\big(\frac{|d-p|}{p}\big)^{p}}$ for all ${u \in W^{1,p}_{\mu,0}(\Omega)}$ if the dimension d ≥ 2, 1 < p < d, and for all ${u \in W^{1,p}_{\mu,0}(\Omega{\setminus}\{0\})}$ if p > d ≥ 1. Here we assume that Ω is an open subset of ${\mathbb{R}^{d}}$ with ${0 \in \Omega}$ , A is a real d × d-symmetric positive definite matrix, c > 0, and $$ d \mu: = \rho(x) \,dx \qquad \textrm{with} \quad \rho(x) = c \cdot \exp(-\frac{1}{p}(x^{t}Ax)^{p/2}), \quad x \in\Omega.\quad \quad (2) $$ If p > d ≥ 1, then we can deduce from (1) a weighted Poincaré inequality on ${W^{1,p}_{\mu,0}(\Omega \setminus\{0\})}$ . Due to the optimality of the Hardy constant in (1), we can establish the nonexistence (locally in time) of positive weak solutions of a p-Kolmogorov parabolic equation perturbed by a singular potential in dimension d = 1, for 1 < p <  + ∞, and when Ω =  ]0, + ∞[.     

Regularity of semi-stable solutions to fourth order nonlinear eigenvalue problems on general domains

We examine the fourth order problem $\Delta ^2 u = \lambda f(u) $ in $ \Omega $ with $ \Delta u = u =0 $ on $ {\partial \Omega }$ , where $ \lambda > 0$ is a parameter, $ \Omega $ is a bounded domain in $\mathbb{R }^N$ and where $f$ is one of the following nonlinearities: $ f(u)=e^u$ , $ f(u)=(1+u)^p $ or $ f(u)= \frac{1}{(1-u)^p}$ where $ p>1$ . We show the extremal solution is smooth, provided $$\begin{aligned} N < 2 + 4 \sqrt{2} + 4 \sqrt{ 2 - \sqrt{2}} \approx 10.718 \text{ when} f(u)=e^u, \end{aligned}$$ and $$\begin{aligned} N < \frac{4p}{p-1} + \frac{4(p+1)}{p-1} \left( \sqrt{ \frac{2p}{p+1}} + \sqrt{ \frac{2p}{p+1} - \sqrt{ \frac{2p}{p+1}}} - \frac{1}{2} \right) \end{aligned}$$ when $ f(u)=(u+1)^p$ . New results are also obtained in the case where $ f(u)=(1-u)^{-p}$ . These are substantial improvements to various results on critical dimensions obtained recently by various authors. To do that, we derive a new stability inequality satisfied by minimal solutions of the above equation, which is more amenable to estimates as it allows a method of proof reminiscent of the second order case.