Stability of nonlocal difference schemes in a subspace

We consider a family of two-layer difference schemes for the heat equation with nonlocal boundary conditions containing the parameter γ. In some interval γ ∈ (1, γ ₊), the spectrum of the main difference operator contains a unique eigenvalue λ ₀ in the left complex half-plane, while the remaining eigenvalues λ ₁, λ ₂, …, λ _N?1 lie in the right half-plane. The corresponding grid space H _N is represented as the direct sum H _N = H ₀⊕H _N?1 of a one-dimensional subspace and the subspace H _N?1 that is the linear span of eigenvectors µ⁽¹⁾, µ⁽²⁾, …, µ^(N?1). We introduce the notion of stability in the subspace H _N?1 and derive a stability criterion.     

Series by Haar system with monotone coefficients

Let χ = {χ _n } _n=0 ^∞ be the Haar system normalized in L ²(0, 1) and M = {M _s } _s=1 ^∞ be an arbitrary, increasing sequence of nonnegative integers. For any subsystem of χ of the form {φ _k } = χS = {χ _n } _n∈S , where S = S(M) = {n _k } _k=1 ^∞ = {n ∈ V[p]: p ∈ M}, V[0] = {1, 2} and V[p] = {2 ^p + 1, 2 ^p + 2, …, 2 ^p+1} for p = 1, 2, … a series of the form Σ _i=1 ^∞ a _i φ _i with a _i ↘ 0 is constructed, that is universal with respect to partial series in all classes L ^r (0, 1), r ∈ (0, 1), in the sense of a.e. convergence and in the metric ofL ^r (0, 1). The constructed series is universal in the class of all measurable, finite functions on [0, 1] in the sense of a.e. convergence. It is proved that there exists a series by Haar system with decreasing coefficients, which has the following property: for any ? > 0 there exists a measurable function µ(x), x ∈ [0, 1], such that 0 ≤ µ(x) ≤ 1 and |{x ∈ [0, 1], µ(x) ≠ = 1}| < ?, and the series is universal in the weighted space L _µ[0, 1] with respect to subseries, in the sense of convergence in the norm of L _µ[0, 1].     

On homogenization for non-self-adjoint locally periodic elliptic operators

In this note we consider the homogenization problem for a matrix locally periodic elliptic operator on R ^d of the form A ^ε = ?divA(x, x/ε)?. The function A is assumed to be Hölder continuous with exponent s ∈ [0, 1] in the “slow” variable and bounded in the “fast” variable. We construct approximations for (A ^ε ? μ)^?1, including one with a corrector, and for (?Δ) ^s/2(A ^ε ? μ)^?1 in the operator norm on L ₂(R ^d ) ⁿ . For s ≠ 0, we also give estimates of the rates of approximation.     

Perturbations of strongly continuous operator semigroups,and matrix Muckenhoupt weights

Let A and A ₀ be linear continuously invertible operators on a Hilbert space ? such that A ^?1 ? A ₀ ^?1 has finite rank. Assuming that σ(A ₀) = ? and that the operator semigroup V ₊(t) = exp{iA ₀ t}, t ≥ 0, is of class C ₀, we state criteria under which the semigroups U _±(t) = exp{±iAt}, t ≥ 0, are of class C ₀ as well. The analysis in the paper is based on functional models for nonself-adjoint operators and techniques of matrix Muckenhoupt weights.     

Paranormal Elements in Normed Algebra

For a normed algebra A and natural numbers k we introduce and investigate the ∥ · ∥ closed classes P _k (A). We show that P₁(A) is a subset of P _k (A) for all k. If T in P₁(A), then Tⁿ lies in P₁(A) for all natural n. If A is unital, U, V ∈ A are such that ∥U∥ = ∥V∥ = 1, VU = I and T lies in P _k (A), then UTV lies in P _k (A) for all natural k. Let A be unital, then 1) if an element T in P₁(A) is right invertible, then any right inverse element T^?1 lies in P₁(A); 2) for ßßIßß = 1 the class P₁(A) consists of normaloid elements; 3) if the spectrum of an element T, T ∈ P₁(A) lies on the unit circle, then ∥TX∥ = ∥X∥ for all X ∈ A. If A = B(H), then the class P₁(A) coincides with the set of all paranormal operators on a Hilbert space H.     

Ranks of matrices with few distinct entries

It is consistent that P(ω ₁) is the union of less than \({2^{{\aleph _1}}}\) parts such that if A ₀,..., A _n?1, B ₀,..., B _m?1 are distinct elements of the same part, then |A ₀ ∩ · · · ∩ A _n?1 ∩ (ω ₁ ? B ₀) ∩ · · ·∩ (ω ₁ ? B _m?1)| = N₁.     

On Templeman averages and variation functions

Let S be a countable semigroup acting in a measure-preserving fashion (g ? T _g ) on a measure space (Ω, A, µ). For a finite subset A of S, let |A| denote its cardinality. Let (A _k ) _k=1 ^∞ be a sequence of subsets of S satisfying conditions related to those in the ergodic theorem for semi-group actions of A. A. Tempelman. For A-measureable functions f on the measure space (Ω, A, μ) we form for k ≥ 1 the Templeman averages \(\pi _k (f)(x) = \left| {A_k } \right|^{ - 1} \sum\nolimits_{g \in A_k } {T_g f(x)}\) and set V _q f(x) = (Σ _k≥1|π _k+1(f)(x) ? π _k (f)(x)|q)^1/q when q ∈ (1, 2]. We show that there exists C > 0 such that for all f in L ¹(Ω, A, µ) we have µ({x ∈ Ω: V _q f(x) > λ}) ≤ C(∫_Ω | f | dµ/λ). Finally, some concrete examples are constructed.     

Averaging of systems of differential inclusions with slow and fast variables

We consider the system of differential inclusions

$$\dot x \in \mu F(t, x, y, \mu ), x(0) = x_0 , \dot y \in G(t, x, y, \mu ), y(0) = y_0 $$

, where F,G: D → Kυ (\(R^{m_1 } \)), Kυ (\(R^{m_2 } \)) are mappings into the sets of nonempty convex compact sets in the Euclidean spaces \(R^{m_1 } \) and \(R^{m_2 } \), respectively, D = R ₊ × \(R^{m_1 } \) × \(R^{m_2 } \) × [0, a], a > 0, and µ is a small parameter. The functions F and G and the right-hand side of the averaged problem \(\dot u\) ∈ µF ₀(u), u(0) = x ₀, F ₀(u) ∈ Kυ (\(R^{m_1 } \)), satisfy the one-sided Lipschitz condition with respect to the corresponding phase variables. Under these and some other conditions, we prove that, for each ? > 0, there exists a µ > 0 such that, for an arbitrary µ ∈ (0, µ₀] and any solution x _µ(·), y _µ(·) of the original problem, there exists a solution u _µ(·) of the averaged problem such that ∥x _µ(t) ? y _µ(t) ∥ ≤ ? for t ∈ [0, 1/µ]. Furthermore, for each solution u _µ(·)of the averaged problem, there exists a solution x _µ(·), y _µ(·) of the original problem with the same estimate.

     

Inversion and subspaces of a finite field

Consider two F _q -subspaces A and B of a finite field, of the same size, and let A ^?1 denote the set of inverses of the nonzero elements of A. The author proved that A ^?1 can only be contained in A if either A is a subfield, or A is the set of trace zero elements in a quadratic extension of a field. Csajbók refined this to the following quantitative statement: if A ^?1 ? B, then the bound |A ^?1∩B| ≤ 2|B|/q ? 2 holds. He also gave examples showing that his bound is sharp for |B| ≤ q ³. Our main result is a proof of the stronger bound |A ^?1 ∩ B| ≤ |B|/q · (1 + O _d (q ^?1/2)), for |B| = q ^d with d > 3. We also classify all examples with |B| ≤ q ³ which attain equality or near-equality in Csajbók’s bound.     

Geometry and inequalities of geometric mean

We study some geometric properties associated with the t-geometric means A ?_tB:= A^1/2(A^?1/2BA^?1/2)^tA^1/2 of two n × n positive definite matrices A and B. Some geodesical convexity results with respect to the Riemannian structure of the n × n positive definite matrices are obtained. Several norm inequalities with geometric mean are obtained. In particular, we generalize a recent result of Audenaert (2015). Numerical counterexamples are given for some inequality questions. A conjecture on the geometric mean inequality regarding m pairs of positive definite matrices is posted.     

On idempotent <Emphasis Type="Italic">τ</Emphasis>-measurable operators affiliated to a von Neumann algebra

Let τ be a faithful normal semifinite trace on a von Neumann algebra M, let p, 0 < p < ∞, be a number, and let L_p(M, τ) be the space of operators whose pth power is integrable (with respect to τ). Let P and Q be τ-measurable idempotents, and let A ≡ P ? Q. In this case, 1) if A ≥ 0, then A is a projection and QA = AQ = 0; 2) if P is quasinormal, then P is a projection; 3) if Q ∈ M and A ∈ L_p(M, τ), then A² ∈ L_p(M, τ). Let n be a positive integer, n > 2, and A = Aⁿ ∈ M. In this case, 1) if A ≠ 0, then the values of the nonincreasing rearrangement μ_t(A) belong to the set {0} ∪ [‖A^n?2‖^?1, ‖A‖] for all t > 0; 2) either μ_t(A) ≥ 1 for all t > 0 or there is a t₀ > 0 such that μ_t(A) = 0 for all t > t₀. For every τ-measurable idempotent Q, there is aunique rank projection P ∈ M with QP = P, PQ = Q, and PM = QM. There is a unique decomposition Q = P + Z, where Z² = 0, ZP = 0, and PZ = Z. Here, if Q ∈ L_p(M, τ), then P is integrable, and τ(Q) = τ(P) for p = 1. If A ∈ L₁(M, τ) and if A = A³ and A ? A² ∈ M, then τ(A) ∈ R.     

Counting sets with small sumset and applications

We study the number of k-element sets A? {1,...,N} with |A+A| ≤ K|A| for some (fixed) K > 0. Improving results of the first author and of Alon, Balogh, Samotij and the second author, we determine this number up to a factor of 2 ^{o ( k )} N ^{o (1)} for most N and k. As a consequence of this and a further new result concerning the number of sets A??/N? with |A+A| ≤ c|A|², we deduce that the random Cayley graph on ?/N? with edge density ½ has no clique or independent set of size greater than (2+o(1)) log₂ N, asymptotically the same as for the Erd?s-Rényi random graph. This improves a result of the first author from 2003 in which a bound of 160log₂ N was obtained. As a second application, we show that if the elements of A ? ? are chosen at random, each with probability 1/2, then the probability that A+A misses exactly k elements of ? is equal to (2+O(1))^?k/2 as k → ∞.     

The Number of <Emphasis Type="Italic">k</Emphasis>-Sumsets in an Abelian Group

Let G be an abelian group of order n. The sum of subsets A₁,...,A_k of G is defined as the collection of all sums of k elements from A₁,...,A_k; i.e., A₁ + A₂ + · · · + A_k = {a₁ + · · · + a_k | a₁ ∈ A₁,..., a_k ∈ A_k}. A subset representable as the sum of k subsets of G is a k-sumset. We consider the problem of the number of k-sumsets in an abelian group G. It is obvious that each subset A in G is a k-sumset since A is representable as A = A₁ + · · · + A_k, where A₁ = A and A₂ = · · · = A_k = {0}. Thus, the number of k-sumsets is equal to the number of all subsets of G. But, if we introduce a constraint on the size of the summands A₁,...,A_k then the number of k-sumsets becomes substantially smaller. A lower and upper asymptotic bounds of the number of k-sumsets in abelian groups are obtained provided that there exists a summand A_i such that |A_i| = n log^qn and |A₁ +· · ·+ A_i-1 + A_i+1 + · · ·+A_k| = n log^qn, where q = -1/8 and i ∈ {1,..., k}.     

Blow-up of critical norms for the 3-D Navier-Stokes equations

Let u =(uh, u3) be a smooth solution of the 3-D Navier-Stokes equations in R3× [0, T). It was proved that if u3 ∈ L∞(0, T;˙B-1+3/p p,q(R3)) for 3 p, q ∞ and uh∈ L∞(0, T; BMO-1(R3)) with uh(T) ∈ VMO-1(R3), then u can be extended beyond T. This result generalizes the recent result proved by Gallagher et al.(2016), which requires u ∈ L∞(0, T;˙B-1+3/pp,q(R3)). Our proof is based on a new interior regularity criterion in terms of one velocity component, which is independent of interest.     

Decay estimates of discretized Green's functions for Schrdinger type operators

For a sparse non-singular matrix A, generally A~(-1)is a dense matrix. However, for a class of matrices,A~(-1)can be a matrix with off-diagonal decay properties, i.e., |A_(ij)~(-1)| decays fast to 0 with respect to the increase of a properly defined distance between i and j. Here we consider the off-diagonal decay properties of discretized Green's functions for Schr¨odinger type operators. We provide decay estimates for discretized Green's functions obtained from the finite difference discretization, and from a variant of the pseudo-spectral discretization. The asymptotic decay rate in our estimate is independent of the domain size and of the discretization parameter.We verify the decay estimate with numerical results for one-dimensional Schr¨odinger type operators.     

Convolutions of cantor measures without resonance

Denote by µ _a the distribution of the random sum \((1 - a)\sum\nolimits_{j = 0}^\infty {{w_j}{a^j}} \), where P(ω _j = 0) = P(ω _j = 1) = 1/2 and all the choices are independent. For 0 < a < 1/2, the measure µ _a is supported on C _a , the central Cantor set obtained by starting with the closed united interval, removing an open central interval of length (1 ? 2a), and iterating this process inductively on each of the remaining intervals. We investigate the convolutions µ _a * (µ _b ° S _λ ^?1 ), where S _λ (x) = λx is a rescaling map. We prove that if the ratio log b/ log a is irrational and λ ≠ 0, then

$D({\mu _a} * ({\mu _b} \circ S_\lambda ^{ - 1})) = \min ({\dim _H}({C_a}) + {\dim _H}({C_b}),1)$

, where D denotes any of correlation, Hausdorff or packing dimension of a measure.

We also show that, perhaps surprisingly, for uncountably many values of λ the convolution µ_1/4* (µ_1/3 ° S _λ ^?1 ) is a singular measure, although dim _H (C _1/4) + dim _H (C _1/3) > 1 and log(1/3)/ log(1/4) is irrational.     

G-matrices,J-orthogonal matrices,and their sign patterns

A real matrix A is a G-matrix if A is nonsingular and there exist nonsingular diagonal matrices D₁ and D₂ such that A^?T = D₁AD₂, where A^?T denotes the transpose of the inverse of A. Denote by J = diag(±1) a diagonal (signature) matrix, each of whose diagonal entries is +1 or ?1. A nonsingular real matrix Q is called J-orthogonal if Q^TJQ = J. Many connections are established between these matrices. In particular, a matrix A is a G-matrix if and only if A is diagonally (with positive diagonals) equivalent to a column permutation of a J-orthogonal matrix. An investigation into the sign patterns of the J-orthogonal matrices is initiated. It is observed that the sign patterns of the G-matrices are exactly the column permutations of the sign patterns of the J-orthogonal matrices. Some interesting constructions of certain J-orthogonal matrices are exhibited. It is shown that every symmetric staircase sign pattern matrix allows a J-orthogonal matrix. Sign potentially J-orthogonal conditions are also considered. Some examples and open questions are provided.     

Elliptic function of level 4

The article is devoted to the theory of elliptic functions of level n. An elliptic function of level n determines a Hirzebruch genus called an elliptic genus of level n. Elliptic functions of level n are also of interest because they are solutions of the Hirzebruch functional equations. The elliptic function of level 2 is the Jacobi elliptic sine function, which determines the famous Ochanine–Witten genus. It is the exponential of the universal formal group of the form F(u, v) = (u² ? v²)/(uB(v) ? vB(u)), B(0) = 1. The elliptic function of level 3 is the exponential of the universal formal group of the form F(u, v) = (u²A(v) ? v²A(u))/(uA(v)² ? vA(u)²), A(0) = 1, A″(0) = 0. In the present study we show that the elliptic function of level 4 is the exponential of the universal formal group of the form F(u, v) = (u²A(v) ? v²A(u))/(uB(v) ? vB(u)), where A(0) = B(0) = 1 and for B′(0) = A″(0) = 0, A′(0) = A₁, and B″(0) = 2B₂ the following relation holds: (2B(u) + 3A₁u)² = 4A(u)³ ? (3A₁² ? 8B₂)u²A(u)². To prove this result, we express the elliptic function of level 4 in terms of the Weierstrass elliptic functions.     

Commutators of Riesz transforms of magnetic Schrödinger operators

Let \({A=-(\nabla-i{\vec a})\cdot (\nabla-i{\vec a}) +V}\) be a magnetic Schrödinger operator acting on \({L^2({\mathbb R}^n)}\), n ≥  1, where \({{\vec a}=(a_1, \ldots, a_n)\in L^2_{\rm loc}({\mathbb R}^n, {\mathbb R}^n)}\) and \({0\leq V\in L^1_{\rm loc}({\mathbb R}^n)}\). In this paper, we show that when a function \({b\in {\rm BMO}({\mathbb R}^n)}\), the commutators [b, T _k ]f = T _k (b f) ? b T _k f, k = 1, . . . , n, are bounded on \({L^p({\mathbb R}^n)}\) for all 1 < p < 2, where the operators T _k are Riesz transforms (?/?x _k  ? i a _k )A ^?1/2 associated with A.     

Extended phase properties and stability analysis of RKN-type integrators for solving general oscillatory second-order initial value problems

In this paper, we study in detail the phase properties and stability of numerical methods for general oscillatory second-order initial value problems whose right-hand side functions depend on both the position y and velocity y ^'. In order to analyze comprehensively the numerical stability of integrators for oscillatory systems, we introduce a novel linear test model y ^?(t) + ? ² y(t) + µ y ^'(t)=0 with µ<2?. Based on the new model, further discussions and analysis on the phase properties and stability of numerical methods are presented for general oscillatory problems. We give the new definitions of dispersion and dissipation which can be viewed as an essential extension of the traditional ones based on the linear test model y ^?(t) + ? ² y(t)=0. The numerical experiments are carried out, and the numerical results showthatthe analysisofphase properties and stability presentedinthispaper ismoresuitableforthenumericalmethodswhentheyareappliedtothe generaloscillatory second-order initial value problem involving both the position and velocity.