Approximations of the exponential of an operator with periodic coefficients

We consider the operator exponential e ^−tA , t > 0, where A is a selfadjoint positive definite operator corresponding to the diffusion equation in \mathbbRⁿ {\mathbb{R}^n} with measurable 1-periodic coefficients, and approximate it in the operator norm ||   ·   ||_{L²( \mathbbRⁿ ) ? L²( \mathbbRⁿ )} {\left\| {\; \cdot \;} \right\|_{{{L^2}\left( {{\mathbb{R}^n}} \right) \to {L^2}\left( {{\mathbb{R}^n}} \right)}}} with order O( t ^{- \fracm2} ) O\left( {{t^{{ - \frac{m}{2}}}}} \right) as t → ∞, where m is an arbitrary natural number. To construct approximations we use the homogenized parabolic equation with constant coefficients, the order of which depends on m and is greater than 2 if m > 2. We also use a collection of 1-periodic functions N _α (x), x ? \mathbbRⁿ x \in {\mathbb{R}^n} , with multi-indices α of length | a| \leqslant m \left| \alpha \right| \leqslant m , that are solutions to certain elliptic problems on the periodicity cell. These results are used to homogenize the diffusion equation with ε-periodic coefficients, where ε is a small parameter. In particular, under minimal regularity conditions, we construct approximations of order O(ε ^m ) in the L ²-norm as ε → 0. Bibliography: 14 titles.     

Anharmonic Oscillators in the Complex Plane, $\boldsymbol{\mathcal{PT}}$-symmetry,and Real Eigenvalues

For integers m ≥ 3 and 1 ≤ ℓ ≤ m − 1, we study the eigenvalue problems − u ^″(z) + [( − 1)^ℓ(iz) ^m  − P(iz)]u(z) = λu(z) with the boundary conditions that u(z) decays to zero as z tends to infinity along the rays argz=-\fracp2±\frac(l+1)pm+2\arg z=-\frac{\pi}{2}\pm \frac{(\ell+1)\pi}{m+2} in the complex plane, where P is a polynomial of degree at most m − 1. We provide asymptotic expansions of the eigenvalues λ _n . Then we show that if the eigenvalue problem is PT\mathcal{PT}-symmetric, then the eigenvalues are all real and positive with at most finitely many exceptions. Moreover, we show that when gcd(m,l)=1\gcd(m,\ell)=1, the eigenvalue problem has infinitely many real eigenvalues if and only if one of its translations or itself is PT\mathcal{PT}-symmetric. Also, we will prove some other interesting direct and inverse spectral results.     

Application of a Bernstein type inequality to rational interpolation in the Dirichlet space

We prove a Bernstein type inequality involving the Bergman and Hardy norms for rational functions in the unit disk \mathbb D {\mathbb D} that have at most n poles all of which are outside the disk \frac1r \mathbb D \frac{1}{r} {\mathbb D} , 0 < r < 1. The asymptotic sharpness of this inequality is shown as n → ∞ and r → 1—. We apply our Bernstein type inequality to an efficient Nevanlinna–Pick interpolation problem in the standard Dirichlet space constrained by the H2-nom. Bibliography: 14 titles.     

On uniqueness problems related to the Fokker–Planck–Kolmogorov equation for measures

We survey recent results related to uniqueness problems for parabolic equations for measures. We consider equations of the form ∂ _t μ = L ^* μ for bounded Borel measures on ℝ ^d  × (0, T), where L is a second order elliptic operator, for example, Lu = D_xu + ( b,?_xu ) Lu = {\Delta_x}u + \left( {b,{\nabla_x}u} \right) , and the equation is understood as the identity

Semidiscrete and asymptotic approximations for the nonstationary radiative–conductive heat transfer problem in a periodic system of grey heat shields

We consider semidiscrete and asymptotic approximations to a solution to the nonstationary nonlinear initial-boundary-value problem governing the radiative–conductive heat transfer in a periodic system consisting of n grey parallel plate heat shields of width ε = 1/n, separated by vacuum interlayers. We study properties of special semidiscrete and homogenized problems whose solutions approximate the solution to the problem under consideration. We establish the unique solvability of the problem and deduce a priori estimates for the solutions. We obtain error estimates of order O( ?{e} ) O\left( {\sqrt {\varepsilon } } \right) and O(ε) for semidiscrete approximations and error estimates of order O( ?{e} ) O\left( {\sqrt {\varepsilon } } \right) and O(ε ^3/4) for asymptotic approximations. Bibliography: 9 titles.     

On some imbedding relations between certain sequence spaces

We introduce a space of sequences l_p^l \ell_p^\lambda of nonabsolute type, which is a p-normed space and a BK-space in the cases of 0 < p < 1 and 1 ≤ p < ∞; respectively. Further, we deduce some imbedding relations and construct a basis for the space l_p^l \ell_p^\lambda , where 1 ≤ p < ∞.     

The primes contain arbitrarily long polynomial progressions

We establish the existence of infinitely many polynomial progressions in the primes; more precisely, given any integer-valued polynomials P ₁, …, P _k  ∈ Z[m] in one unknown m with P ₁(0) = … = P _k (0) = 0, and given any ε > 0, we show that there are infinitely many integers x and m, with 1 \leqslant m \leqslant x^e1 \leqslant m \leqslant x^\varepsilon, such that x + P ₁(m), …, x + P _k (m) are simultaneously prime. The arguments are based on those in [18], which treated the linear case P _j  = (j − 1)m and ε = 1; the main new features are a localization of the shift parameters (and the attendant Gowers norm objects) to both coarse and fine scales, the use of PET induction to linearize the polynomial averaging, and some elementary estimates for the number of points over finite fields in certain algebraic varieties.     

Discrete Fourier Analysis on Fundamental Domain and Simplex of A d Lattice in d-Variables

A discrete Fourier analysis on the fundamental domain Ω _d of the d-dimensional lattice of type A _d is studied, where Ω₂ is the regular hexagon and Ω₃ is the rhombic dodecahedron, and analogous results on d-dimensional simplex are derived by considering invariant and anti-invariant elements. Our main results include Fourier analysis in trigonometric functions, interpolation and cubature formulas on these domains. In particular, a trigonometric Lagrange interpolation on the simplex is shown to satisfy an explicit compact formula and the Lebesgue constant of the interpolation is shown to be in the order of (log n) ^d . The basic trigonometric functions on the simplex can be identified with Chebyshev polynomials in several variables already appeared in literature. We study common zeros of these polynomials and show that they are nodes for a family of Gaussian cubature formulas, which provides only the second known example of such formulas.     

On the trigonometric interpolation and the entire interpolation

In this paper, we study a kind of interpolation problems on a given nodal set by trigonometric polynomials of order n and entire functions of exponential type according as the nodal set is respectively. We established some equivalent conditions and found the explicit forms of some interpolation functions on the interpolation problems. As a special case, the explicit forms of fundamential functions of (0,m)-interpolat on by trigonometric case or entire functions case (in B² _σ) respectively, if they exist, may follow from our results. Besides, we also considered the convergence of the interpolation functions at above stated. Suported by the Natural Youth Science Foundation of Beijing Normal University.     

Heat flows for a nonconvex Signorini type problem in $$ {\mathbb{R}^{N}} $$

We prove the existence of a global heat flow u : Ω ×  \mathbbR⁺ ? \mathbbR^N {\mathbb{R}^{+}} \to {\mathbb{R}^{N}}, N > 1, satisfying a Signorini type boundary condition u(∂Ω ×  \mathbbR⁺ {\mathbb{R}^{+}}) ⊂  \mathbbRⁿ {\mathbb{R}^{n}}), n \geqslant 2 n \geqslant 2 , and \mathbbR^N {\mathbb{R}^{N}}) with boundary ∂ [`(W)] \bar{\Omega } such that φ(∂Ω) ⊂ \mathbbR^N {\mathbb{R}^{N}} is given by a smooth noncompact hypersurface S. Bibliography: 30 titles.     

Piecewise-linear,discontinuous Galerkin method for a fractional diffusion equation

We use a piecewise-linear, discontinuous Galerkin method for the time discretization of a fractional diffusion equation involving a parameter in the range − 1 < α < 0. Our analysis shows that, for a time interval (0,T) and a spatial domain Ω, the error in L_￥((0,T);L₂(W))L_\infty\bigr((0,T);L_2(\Omega)\bigr) is of order k ^2 + α , where k denotes the maximum time step. Since derivatives of the solution may be singular at t = 0, our result requires the use of non-uniform time steps. In the limiting case α = 0 we recover the known O(k ²) convergence for the classical diffusion (heat) equation. We also consider a fully-discrete scheme that employs standard (continuous) piecewise-linear finite elements in space, and show that the additional error is of order h ²log(1/k). Numerical experiments indicate that our O(k ^2 + α ) error bound is pessimistic. In practice, we observe O(k ²) convergence even for α close to − 1.     

The generalized bisymmetric solutions of the matrix equation A 1 X 1 B 1 + A 2 X 2 B 2 + ? + A l X l B l = C and its optimal approximation

For fixed generalized reflection matrix P, i.e. P ^T  = P, P ² = I, then matrix X is said to be generalized bisymmetric, if X = X ^T  = PXP. In this paper, an iterative method is constructed to find the generalized bisymmetric solutions of the matrix equation A ₁ X ₁ B ₁ + A ₂ X ₂ B ₂ + ⋯ + A _l X _l B _l  = C where [X ₁,X ₂, ⋯ ,X _l ] is real matrices group. By this iterative method, the solvability of the matrix equation can be judged automatically. When the matrix equation is consistent, for any initial generalized bisymmetric matrix group , a generalized bisymmetric solution group can be obtained within finite iteration steps in the absence of roundoff errors, and the least norm generalized bisymmetric solution group can be obtained by choosing a special kind of initial generalized bisymmetric matrix group. In addition, the optimal approximation generalized bisymmetric solution group to a given generalized bisymmetric matrix group in Frobenius norm can be obtained by finding the least norm generalized bisymmetric solution group of the new matrix equation , where . Given numerical examples show that the algorithm is efficient. Research supported by: (1) the National Natural Science Foundation of China (10571047) and (10771058), (2) Natural Science Foundation of Hunan Province (06JJ2053), (3) Scientific Research Fund of Hunan Provincial Education Department(06A017).     

Regularity of the Extremal Solution of Some Nonlinear Elliptic Problems Involving the <Emphasis Type="Italic">p</Emphasis>-Laplacian

We consider the equation on a smooth bounded domain of with zero Dirichlet boundary conditions where p ≥ 2, λ > 0 and f satisfies typical assumptions in the subject of extremal solutions. We prove that, for such general nonlinearities f, the extremal solution u ^* belongs to L ^∞(Ω) if N < p + p/(p − 1) and if N < p(1 + p/(p − 1)). This work was partially supported by MCyT BMF 2002-04613-CO3-02.     

On uniqueness problems related to elliptic equations for measures

We consider equations of the form L*μ = 0 for bounded measures on _boxclose^d {\mathbb{R}^{d}} , where L is a second order elliptic operator, for example, Lu = Δu + (b,∇u), and the equation is understood as the identity

Homogenization with corrector for a multidimensional periodic elliptic operator near an edge of an inner gap

For a second order differential operator A \mathcal{A} _ε  = −div g(x/ε)∇ + ε ⁻²p(x/ε) in L ₂(ℝ ^d ) with periodic coefficients and small parameter ε > 0 we study an approximation of the resolvent of A \mathcal{A} _ε at a point close to an edge of an inner gap in the spectrum of A \mathcal{A} _ε . Under certain regularity conditions, we construct an approximation (with a first order corrector taken into account) for the resolvent with error estimate of order O(ε ²). We show that a proper effective operator and a proper corrector are associated to each (regular) edge of the gap. Bibliography: 14 titles.     

On the boundedness of one recurrent sequence in a banach space

We establish necessary and sufficient conditions under which a sequence x ₀ = y ₀ , x _n+1 = Ax _n  + y _n+1 , n ≥ 0, is bounded for each bounded sequence { y_n :n \geqslant 0 } ì { x ? è_{n = 1}^￥ D( Aⁿ ) |sup_{n \geqslant 0} || Aⁿ x || < ￥ }\left\{ {y_n :n \geqslant 0} \right\} \subset \left\{ {\left. {x \in \bigcup\nolimits_{n = 1}^\infty {D\left( {A^n } \right)} } \right|\sup _{n \geqslant 0} \left\| {A^n x} \right\| < \infty } \right\}, where A is a closed operator in a complex Banach space with domain of definition D(A) .     

Existence of solution of the pullback equation involving volume forms

Let Ω ⊂ ℝ ⁿ be a smooth, bounded domain. We study the existence and regularity of diffeomorphisms of Ω satisfying the volume form equation

Numerical approximations and solution techniques for the space-time Riesz–Caputo fractional advection-diffusion equation

In this paper, we consider a space-time Riesz–Caputo fractional advection-diffusion equation. The equation is obtained from the standard advection-diffusion equation by replacing the first-order time derivative by the Caputo fractional derivative of order α ∈ (0,1], the first-order and second-order space derivatives by the Riesz fractional derivatives of order β ₁ ∈ (0,1) and β ₂ ∈ (1,2], respectively. We present an explicit difference approximation and an implicit difference approximation for the equation with initial and boundary conditions in a finite domain. Using mathematical induction, we prove that the implicit difference approximation is unconditionally stable and convergent, but the explicit difference approximation is conditionally stable and convergent. We also present two solution techniques: a Richardson extrapolation method is used to obtain higher order accuracy and the short-memory principle is used to investigate the effect of the amount of computations. A numerical example is given; the numerical results are in good agreement with theoretical analysis.     

On the Scope of Averaging for Frankl’s Conjecture

Let be a union-closed family of subsets of an m-element set A. Let . For b ∈ A let w(b) denote the number of sets in containing b minus the number of sets in not containing b. Frankl’s conjecture from 1979, also known as the union-closed sets conjecture, states that there exists an element b ∈ A with w(b) ≥ 0. The present paper deals with the average of the w(b), computed over all b ∈ A. is said to satisfy the averaged Frankl’s property if this average is non-negative. Although this much stronger property does not hold for all union-closed families, the first author (Czédli, J Comb Theory, Ser A, 2008) verified the averaged Frankl’s property whenever n ≥ 2 ^m  − 2 ^m/2 and m ≥ 3. The main result of this paper shows that (1) we cannot replace 2 ^m/2 with the upper integer part of 2 ^m /3, and (2) if Frankl’s conjecture is true (at least for m-element base sets) and then the averaged Frankl’s property holds (i.e., 2 ^m/2 can be replaced with the lower integer part of 2 ^m /3). The proof combines elementary facts from combinatorics and lattice theory. The paper is self-contained, and the reader is assumed to be familiar neither with lattices nor with combinatorics. This research was partially supported by the NFSR of Hungary (OTKA), grant no. T 049433, T 48809 and K 60148.     

Diagnostic tests for local coalescences of variables in Boolean functions

The article studies diagnostic tests for local k -fold coalescences of variables in Boolean functions f( [(x)\tilde]ⁿ )( 1 ￡ k ￡ n,  1 ￡ t ￡ 2^2k ) f\left( {{{\tilde{x}}^n}} \right)\left( {1 \leq k \leq n,\;1 \leq t \leq {2^{{2^k}}}} \right) . Upper and lower bounds are proved for the Shannon function of the length of the diagnostic test for local k -fold coalescences generated by the system of functions F_t^k \Phi_t^k . The Shannon function of the length of a complete diagnostic test for local k -fold coalescences behaves asymptotically as 2 ^k (n − k + 1) for n → ∞, k → ∞.