On stochastic and quasistochastic methods for solving equations

A modification of the quasi-Monte Carlo method previously suggested by the authors is discussed. This modification is successfully applied to solve integral equations of the second kind thanks to a number of advantages over the quasi-Monte Carlo method. Thus, solving these equations by the quasi-Monte Carlo method requires estimating the sum of a series whose terms are integrals with infinitely increasing constructive dimension. As is known, this is one of the main factors hindering the application of the quasi-Monte Carlo method. Another difficulty is the dominated convergence condition, which ensures the absolute convergence of the series whose sum is to be estimated.The modification under consideration makes it possible to overcome the first difficulty and relax the dominated convergence condition.

Two families of estimates are suggested, which can be applied in the modified algorithm of the quasi-Monte Carlo method. The integral equation

$\phi (x) = \int {k(x,y)\phi (y)\mu (dy) + f(x)(\bmod \mu ),} $

is considered, where x ∈ D ? ? ^s and f and k are given functions defined on the supports of the measures μ and μ ? μ. The values of ? _n (x) = ∫k(x, y)?_{n - 1}(y)μ(dy) + f(x)are estimated step by step. The first group of estimates serves to evaluate ? _n at a fixed point x’:

$\xi _1 (x') = \frac{1}{N}\sum\limits_{j = 1}^N {\xi _1^n (y_j )} ,where\xi _1^n (y) = \frac{{k(x',y)\hat \phi _{n - 1} (y)}}{{p_{n - 1} (y)}} + f(x'),$

here, the y _j are distributed with density p _n?1 and \(\hat \phi _{n - 1} (y)\) is the estimate obtained at the preceding step. The other group of estimates makes it possible to evaluate ? _n at random points.

Optimal parameters of the estimates are found and the corresponding theorems are proved.The theory has been verified on the example of a difference analogue of the Navier-Stokes equation; experimental results are presented.     

On the Darboux Integrability of the Hindmarsh–Rose Burster

We study the Hindmarsh–Rose burster which can be described by the differential system = y-x~3+ bx~2+ I-z,  = 1-5 x2~-y, z = μ(s(x-x_0)-z),where b, I, μ, s, x_0 are parameters. We characterize all its invariant algebraic surfaces and all its exponential factors for all values of the parameters. We also characterize its Darboux integrability in function of the parameters. These characterizations allow to study the global dynamics of the system when such invariant algebraic surfaces exist.     

On the Surface Integral Approximation of the Second Derivatives of the Potential of a Bulk Charge Located in a Layer of Small Thickness

We consider a bulk charge potential of the form

$$u(x) = \int\limits_\Omega {g(y)F(x - y)dy,x = ({x_1},{x_2},{x_3}) \in {\mathbb{R}^3},} $$

where Ω is a layer of small thickness h > 0 located around the midsurface Σ, which can be either closed or open, and F(x ? y) is a function with a singularity of the form 1/|x ? y|. We prove that, under certain assumptions on the shape of the surface Σ, the kernel F, and the function g at each point x lying on the midsurface Σ (but not on its boundary), the second derivatives of the function u can be represented as

$$\frac{{{\partial ^2}u(x)}}{{\partial {x_i}\partial {x_j}}} = h\int\limits_\Sigma {g(y)\frac{{{\partial ^2}F(x - y)}}{{\partial {x_i}\partial {x_j}}}} dy - {n_i}(x){n_j}(x)g(x) + {\gamma _{ij}}(x),i,j = 1,2,3,$$

where the function γ_ij(x) does not exceed in absolute value a certain quantity of the order of h², the surface integral is understood in the sense of Hadamard finite value, and the n_i(x), i = 1, 2, 3, are the coordinates of the normal vector on the surface Σ at a point x.

     

Weighted polyharmonic equation with Navier boundary conditions in a half space

We study positive solutions of the following polyharmonic equation with Hardy weights associated to Navier boundary conditions on a half space:?????(-?)~mu(x)=u~p(x)/|x|~s,in R_+~n,u(x)=-?u(x)=…=(-?)~(m-1)u(x)=0,on ?R_+~n,(0.1)where m is any positive integer satisfying 02mn.We first prove that the positive solutions of(0.1)are super polyharmonic,i.e.,(-?)~iu0,i=0,1,...,m-1.(0.2) For α=2m,applying this important property,we establish the equivalence between (0.1) and the integral equation u(x)=c_n∫R_+~n(1/|x-y|~(n-α)-1/|x~*-y|~(n-α))u~p(y)/|y|~sdy,(0.3) where x~*=(x1,...,x_(n-1),-x_n) is the reflection of the point x about the plane R~(n-1).Then,we use the method of moving planes in integral forms to derive rotational symmetry and monotonicity for the positive solution of(0.3),in whichαcan be any real number between 0 and n.By some Pohozaev type identities in integral forms,we prove a Liouville type theorem—the non-existence of positive solutions for(0.1).     

Blue sky catastrophe in relaxation systems with one fast and two slow variables

We establish conditions under which three-dimensional relaxational systems of the form

$$\dot x = f(x,y,\mu ),\varepsilon \dot y = g(x,y),x = (x_1 ,x_2 ) \in \mathbb{R}^2 ,y \in \mathbb{R},$$

where 0 ≤ ε ? 1, |µ| ? 1, and f, g ∈ C ^∞, exhibit the so-called blue sky catastrophe [the appearance of a stable relaxational cycle whose period and length tend to infinity as µ tends to some critical value µ_*(ε), µ_*(0) = 0].

     

Infinitely many solutions for Kirchhoff equations with sign-changing potential and Hartree nonlinearity

This paper is concerned with the following Kirchhoff-type equations:

$$\begin{aligned} \left\{ \begin{array}{ll} \displaystyle -\big (a+b\int _{\mathbb {R}^{3}}|\nabla u|^{2}\mathrm {d}x\big )\Delta u+ V(x)u+\mu \phi |u|^{p-2}u=f(x, u)+g(x,u), &{} \text{ in } \mathbb {R}^{3},\\ (-\Delta )^{\frac{\alpha }{2}} \phi = \mu |u|^{p}, &{} \text{ in } \mathbb {R}^{3},\\ \end{array} \right. \end{aligned}$$

where \(a>0,~b,~\mu \ge 0\) are constants, \(\alpha \in (0,3)\), \(p\in [2,3+2\alpha )\), the potential V(x) may be unbounded from below and \(\phi |u|^{p-2}u\) is a Hartree-type nonlinearity. Under some mild conditions on V(x), f(x, u) and g(x, u), we prove that the above system has infinitely many nontrivial solutions. Specially, our results cover the general Schrödinger equations, the Kirchhoff equations and the Schrödinger–Poisson system.

     

Positive Solutions of a Third-Order Boundary Value Problem with Full Nonlinearity

This paper is concerned with the existence of positive solutions of the third-order boundary value problem with full nonlinearity

$$\begin{aligned} \left\{ \begin{array}{lll} u'''(t)&{}=f(t,u(t),u'(t),u''(t)),\quad t\in [0,1],\\ u(0)&{}=u'(1)=u''(1)=0, \end{array}\right. \end{aligned}$$

where \(f:[0,1]\times \mathbb {R}^+\times \mathbb {R}^+\times \mathbb {R}^-\rightarrow \mathbb {R}^+\) is continuous. Under some inequality conditions on f as |(x, y, z)| small or large enough, the existence results of positive solution are obtained. These inequality conditions allow that f(t, x, y, z) may be superlinear, sublinear or asymptotically linear on x, y and z as \(|(x,y,z)|\rightarrow 0\) and \(|(x,y,z)|\rightarrow \infty \). For the superlinear case as \(|(x,y,z)|\rightarrow \infty \), a Nagumo-type growth condition is presented to restrict the growth of f on y and z. Our discussion is based on the fixed point index theory in cones.

     

Linear FDEs in the frame of generalized ODEs: variation-of-constants formula

We present a variation-of-constants formula for functional differential equations of the form

$$\dot y = \mathcal{L}\left( t \right)y_t + f\left( {y_t,t} \right),\;y_{t_0}= \varphi $$

, where \(\mathcal{L}\) is a bounded linear operator and φ is a regulated function. Unlike the result by G. Shanholt (1972), where the functions involved are continuous, the novelty here is that the application t \(t \mapsto f\left( {y_t,t} \right)\) is Kurzweil integrable with t in an interval of ?, for each regulated function y. This means that t \(t \mapsto f\left( {y_t,t} \right)\) may admit not only many discontinuities, but it can also be highly oscillating and yet, we are able to obtain a variation-of-constants formula. Our main goal is achieved via theory of generalized ordinary differential equations introduced by J.Kurzweil (1957). As a matter of fact, we establish a variation-of-constants formula for general linear generalized ordinary differential equations in Banach spaces where the functions involved are Kurzweil integrable. We start by establishing a relation between the solutions of the Cauchy problem for a linear generalized ODE of type

$$\frac{{dx}}{{d\tau }} = D\left[ {A\left( t \right)x} \right],x\left( {{t_0}} \right) = \tilde x$$

and the solutions of the perturbed Cauchy problem

$$\frac{{dx}}{{d\tau }} = D\left[ {A\left( t \right)x + F\left( {x,t} \right)} \right],x\left( {{t_0}} \right) = \tilde x$$

Then we prove that there exists a one-to-one correspondence between a certain class of linear generalized ODE and the Cauchy problem for a linear functional differential equations of the form

$$\dot y = \mathcal{L}\left( t \right)y_t,\;y_{t_0} = \varphi$$

, where \(\mathcal{L}\) is a bounded linear operator and φ is a regulated function. The main result comes as a consequence of such results. Finally, because of the extent of generalized ODEs, we are also able to describe the variation-of-constants formula for both impulsive FDEs and measure neutral FDEs.

     

Analytic Hypoellipticity for a New Class of Sums of Squares of Vector Fields in $${\mathcal {R}}^3$$

This paper is devoted to a substantial generalization of previous work on the analytic hypoellipticity of sums of squares \(P=\sum _1^4X^2_j\) of real vector fields with real analytic coefficient in three variables. For p(x, y) quasi-homogeneous in (x, y), consider the vector fields

$$\begin{aligned} X_1 = \frac{\partial }{\partial x}, \quad X_2=-\frac{\partial }{\partial y} + p(x,y)\frac{\partial }{\partial t}, \quad X_3=x^{n_1}\frac{\partial }{\partial t}, \quad X_4=y^{n_2}\frac{\partial }{\partial t}, \end{aligned}$$

\( n_1, n_2 \ne 0\). We show that the operator

$$\begin{aligned} P=\sum _1^4 X_j^2, \end{aligned}$$

well known to be \(C^\infty \)-hypoelliptic, is actually analytic hypoelliptic near the origin in \({\mathcal {R}}^3\).

     

A Spectral Gap Estimate and Applications

We consider the Schrödinger operator

$$ \text{-} \frac{d^{2}}{d x^{2}} + V {\text{on an interval}}~~[a,b]~{\text{with Dirichlet boundary conditions}},$$

where V is bounded from below and prove a lower bound on the first eigenvalue λ ₁ in terms of sublevel estimates: if w _V (y) = |{x ∈ [a, b] : V (x) ≤ y}|, then

$$\lambda_{1} \geq \frac{1}{250} \min\limits_{y > \min V}{\left( \frac{1}{w_{V}(y)^{2}} + y\right)}.$$

The result is sharp up to a universal constant if {x ∈ [a, b] : V(x) ≤ y} is an interval for the value of y solving the minimization problem. An immediate application is as follows: let \({\Omega } \subset \mathbb {R}^{2}\) be a convex domain and let \(u:{\Omega } \rightarrow \mathbb {R}\) be the first eigenfunction of the Laplacian ? Δ on Ω with Dirichlet boundary conditions on ?Ω. We prove

$$\| u \|_{L^{\infty}({\Omega})} \lesssim \frac{1}{\text{inrad}({\Omega})} \left( \frac{\text{inrad}({\Omega})}{\text{diam}({\Omega})} \right)^{1/6} \|u\|_{L^{2}({\Omega})},$$

which answers a question of van den Berg in the special case of two dimensions.

     

Infinitely many periodic solutions for a class of second-order Hamiltonian systems

In this paper we study the existence of infinitely many periodic solutions for second-order Hamiltonian systems

$$\left\{ {\begin{array}{*{20}c} {\ddot u(t) + A(t)u(t) + \nabla F(t,u(t)) = 0,} \\ {u(0) - u(T) = \dot u(0) - \dot u(T) = 0,} \\ \end{array} } \right.$$

, where F(t, u) is even in u, and ?F(t, u) is of sublinear growth at infinity and satisfies the Ahmad-Lazer-Paul condition.

     

Mean Square Estimates for Coefficients of Symmetric Power <Emphasis Type="Italic">L</Emphasis>-Functions

Let L(sym ^j f,s) be the jth symmetric power L-function attached to a holomorphic Hecke eigencuspform f(z) for the full modular group, and \(\lambda_{\mathrm{sym}^{j}f}(n)\) denote its nth coefficient. In this paper we are able to prove that

$\int_{1}^{x}\bigg|\sum_{n\leq y}\lambda_{\mathrm{sym}^{3}f}(n)\bigg|^{2}dy=O\bigl(x^{2}\bigr),$

and

$\int_{1}^{x}\bigg|\sum_{n\leq y}\lambda_{\mathrm{sym}^{4}f}(n)\bigg|^{2}dy=O\bigl(x^{\frac{11}{5}}\log x\bigr).$

     

Conditional exponential stability of a differential system with linear dichotomous Coppel-Conti approximation

We prove the conditional exponential stability of the zero solution of the nonlinear differential system

$$\dot y = A(t)y + f(t,y),{\mathbf{ }}y \in R^n ,{\mathbf{ }}t \geqslant 0,$$

with L ^p -dichotomous linear Coppel-Conti approximation .x = A(t)x whose principal solution matrix X _A (t), X _A (0) = E, satisfies the condition

$$\mathop \smallint \limits_0^t \left\| {X_A (t)P_1 X_A^{ - 1} (\tau )} \right\|^p d\tau + \mathop \smallint \limits_t^{ + \infty } \left\| {X_A (t)P_2 X_A^{ - 1} (\tau )} \right\|^p d\tau \leqslant C_p (A) < + \infty ,{\mathbf{ }}p \geqslant 1,{\mathbf{ }}t \geqslant 0,$$

where P ₁ and P ₂ are complementary projections of rank k ∈ {1, …, n ? 1} and rank n ? k, respectively, and with a higher-order infinitesimal perturbation f:[0, ∞) × U → R ⁿ that is piecewise continuous in t ≥ 0 and continuous in y in some neighborhood U of the origin.

     

On the Existence of Ground State Solutions for Fractional Schrödinger–Poisson Systems with General Potentials and Super-quadratic Nonlinearity

In this article, we are concerned with the following fractional Schrödinger–Poisson system:

$$\begin{aligned} \left\{ \begin{array}{ll} (-\Delta )^{s}u+V(x)u+\phi u=f(u)&{} \quad \hbox {in}~\mathbb {R}^{3},\\ (-\Delta )^{t}\phi =u^2&{} \quad \hbox {in}~\mathbb {R}^{3},\\ \end{array} \right. \end{aligned}$$

where \(0<s\le t<1\), \(2s+2t>3\), and \(f\in C(\mathbb {R},\mathbb {R})\). Under more relaxed assumptions on potential V(x) and f(x), we obtain the existence of ground state solutions for the above problem by adopting some new tricks. Our results here extend the existing study.

     

Aubry–Mather Sets for Relativistic Oscillators with Anharmonic Potentials

In this paper, by the Aubry–Mather theory, it is proved that there are many periodic solutions and usual or generalized quasiperiodic solutions for relativistic oscillator with anharmonic potentials models d/dt(x/(1-|x|~2~(1/2))+ |x|~(α-1)x=p(t),where p(t) ∈ C~0(R~1) is 1-periodic and α 0.     

Strong isochronism of a center and a focus for systems with homogeneous nonlinearities

We suggest a new approach to studying the isochronism of the system

${{dx} \mathord{\left/ {\vphantom {{dx} {dt}}} \right. \kern-\nulldelimiterspace} {dt}} = - y + p_n (x,y),{{dy} \mathord{\left/ {\vphantom {{dy} {dt}}} \right. \kern-\nulldelimiterspace} {dt}} = x + q_n (x,y),$

where p _n and q _n are homogeneous polynomials of degree n. This approach is based on the normal form

${{dX} \mathord{\left/ {\vphantom {{dX} {dt}}} \right. \kern-\nulldelimiterspace} {dt}} = - Y + XS(X,Y),{{dY} \mathord{\left/ {\vphantom {{dY} {dt}}} \right. \kern-\nulldelimiterspace} {dt}} = X + YS(X,Y)$

and its analog in polar coordinates. We prove a theorem on sufficient conditions for the strong isochronism of a center and a focus for the reduced system and obtain examples of centers with strong isochronism of degrees n = 4, 5. The present paper is the first to give examples of foci with strong isochronism for the system in question.

     

On the exceptional sets in Sylvester expansions*

For any x ?? (0, 1], let the series \( {\sum}_{n=1}^{\infty }1/{d}_n(x) \) be the Sylvester expansion of x, where {d _j (x),?j?≥?1} is a sequence of positive integers satisfying d₁(x)?≥?2 and d_j?+?1(x)?≥?d _j (x)(d _j (x)???1)?+?1 for j?≥?1. Suppose ? : ? → ?⁺ is a function satisfying ?(n+1) – ? (n) → ∞ as n → ∞. In this paper, we consider the set

$$ E\left(\phi \right)=\left\{x\kern0.5em \in \left(0,1\right]:\kern0.5em \underset{n\to \infty }{\lim}\frac{\log {d}_n(x)-{\sum}_{j=1}^{n-1}\log {d}_j(x)}{\phi (n)}=1\right\} $$

and quantify the size of the set in the sense of Hausdorff dimension. As applications, for any β > 1 and γ > 0, we get the Hausdorff dimension of the set \( \left\{x\in \kern1em \left(0,1\right]:\kern0.5em {\lim}_{n\to \infty}\left(\log {d}_n(x)-{\sum}_{j=1}^{n-1}\log {d}_j(x)\right)/{n}^{\beta }=\upgamma \right\}, \) and for any τ > 1 and η > 0, we get a lower bound of the Hausdorff dimension of the set \( \left\{x\kern0.5em \in \kern0.5em \left(0,1\right]:\kern1em {\lim}_{n\to \infty}\left(\log {d}_n(x)-{\sum}_{j=1}^{n-1}\log {d}_j(x)\right)/{\tau}^n=\eta \right\}. \)     

Dynamic control for non-linear systems with delay

We consider some class of non-linear systems of the form

$\dot x = A( \cdot )x + \sum\limits_{i = 1}^l {A_i ( \cdot )x(t - \tau _i (t)) + b( \cdot )u} ,$

where A(·) ∈ ? ^{n × n} , A _i (·) ∈ ? ^{n × n} , b(·) ∈ ? ⁿ , whose coefficients are arbitrary uniformly bounded functionals.

A special type of the Lyapunov-Krasovskii functional is used to synthesize dynamic control described by the equation

$\dot u = \rho ( \cdot )u + (m( \cdot ),x),$

where ρ(·) ∈ ?¹, m(·) ∈ ? ⁿ , which makes the system globally asymptotically stable. Also, the situation is considered where the control u enters into the system not directly but through a pulse element performing an amplitude-frequency modulation.

     

Infinitely many solutions for a nonlinear Schrödinger equation with general nonlinearity

We prove the existence of infinitely many solutions for

$$\begin{aligned} - \Delta u + V(x) u = f(u) \quad \text { in } \mathbb {R}^N, \quad u \in H^1(\mathbb {R}^N), \end{aligned}$$

where V(x) satisfies \(\lim _{|x| \rightarrow \infty } V(x) = V_\infty >0\) and some conditions. We require conditions on f(u) only around 0 and at \(\infty \).

     

Hölder seminorm preserving linear bijections and isometries

Let (X, d) be a compact metric and 0 < α < 1. The space Lip ^α (X) of Hölder functions of order α is the Banach space of all functions ? from X into \(\mathbb{K}\) such that ∥?∥ = max{∥?∥_∞, L(?)} < ∞, where

$L(f) = sup\{ \left| {f(x) - f(y)} \right|/d^\alpha (x,y):x,y \in X, x \ne y\} $

is the Hölder seminorm of ?. The closed subspace of functions ? such that

$\mathop {\lim }\limits_{d(x,y) \to 0} \left| {f(x) - f(y)} \right|/d^\alpha (x,y) = 0$

is denoted by lip ^α (X). We determine the form of all bijective linear maps from lip ^α (X) onto lip ^α (Y) that preserve the Hölder seminorm.