A generalized numerical range: the range of a constrained sesquilinear form

Let M_n denote the algebra of all nxn complex matrices. For a given q?C with ∣Q∣≤1, we define and denote the q-numerical range of A?M_n by

W_q (A)={x ^? Ay:x,y?C ⁿ , _{x ^? x?y ^? y=1,x ^? y=q} }

The q-numerical radius is then given by r_q (A)=sup{∣z∣:z?W _q (A)}. When q=1,W _q (A) and r _q (A) reduce to the classical numerical range of A and the classical numerical radius of A, respectively. when q≠0, another interesting quantity associated with W _q (A) is the inner q-numerical radius defined by [rtilde] _q (A)=inf{∣z∣:z?W _q (A)}

In this paper, we describe some basic properties of W _q (A), extending known results on the classical numerical range. We also study the properties of r_q considered as a norm (seminorm if q=0) on M_n .Finally, we characterize those linear operators L on M_n that leave W_q ,r_q of [rtilde]_q invariant. Extension of some of our results to the infinite dimensional case is discussed, and open problems are mentioned.     

On the well-definedness of the order of an ordinary differential equation

It is proved that if the differential equations y ^{( n )}=f(x, y, y′, …, y ^{( n ?1 )}) and y ^{( m )}=g(x, y, y′, …, y ^{( n ?1 )}) have the same particular solutions in a suitable region where f and g are continuous real-valued functions with continuous partial derivatives (alternatively, continuous functions satisfying the classical Lipschitz condition), then n?=?m and the functions f and g are equal. This note could find classroom use in a course on differential equations as enrichment material for the unit on the existence and uniqueness theorems for solutions of initial value problems.     

Uniqueness implies existence and uniqueness conditions for nonlocal boundary value problems for nth order differential equations

For the nth order differential equation, y(n)=f(x,y,y^′,…,y(n−1)), we consider uniqueness implies existence results for solutions satisfying certain nonlocal (k+2)-point boundary conditions, 1?k?n−1. Uniqueness of solutions when k=n−1 is intimately related to uniqueness of solutions when 1?k?n−2. These relationships are investigated as well.     

Some Results for Decomposition Numbers for the Hecke Algebra of Type A with q = −1

Let ? = ?_{?, ?1}(𝔖 _n ) be the Hecke algebra of the symmetric group 𝔖 _n . For partitions λ and ν with ν 2 ? regular, define the Specht module S(λ) and the irreducible module D(ν). Define d _λν = [S(λ): D(ν)] to be the composition multiplicity of D(ν) in S(λ). In this paper we compute the decomposition numbers d _λν for all partitions of the form λ = (a, c, 1 ^b ) and ν 2 ? regular.     

Constructive existence theorems for problems of Thomas-Fermi Type

A minimal positive solution of the Thomas-Fermi problem ? = λt^?1/2 w^3/2, w(0) = 1, w(1) = w(1) is shown to exist for each λ > 0. It is proved that all positive solutions, for a given value of λ, are strictly ordered and that the minimal positive solution w_λ is a decreasing function of λ. Upper and lower analytic bounds for w _λ are given and these bounds are shown to initiate sequences of Picard and Newton iterates which converge monotonically to w _λ. A comparative analysis of the efficiency of the iteration schemes is presented. The methods used are of a general nature and can be applied to a variety of nonlinear boundary value problems of convex type [14].     

Optimal error estimates of the local discontinuous Galerkin method for nonlinear second‐order elliptic problems on Cartesian grids

In this paper, we study the local discontinuous Galerkin (LDG) methods for two‐dimensional nonlinear second‐order elliptic problems of the type u_xx + u_yy = f(x, y, u, u_x, u_y) , in a rectangular region Ω with classical boundary conditions on the boundary of Ω . Convergence properties for the solution and for the auxiliary variable that approximates its gradient are established. More specifically, we use the duality argument to prove that the errors between the LDG solutions and the exact solutions in the L² norm achieve optimal (p + 1)th order convergence, when tensor product polynomials of degree at most p are used. Moreover, we prove that the gradient of the LDG solution is superclose with order p + 1 toward the gradient of Gauss–Radau projection of the exact solution. The results are valid in two space dimensions on Cartesian meshes using tensor product polynomials of degree p ≥ 1 , and for both mixed Dirichlet–Neumann and periodic boundary conditions. Preliminary numerical experiments indicate that our theoretical findings are optimal.     

Solution of large‐scale weighted least‐squares problems

A sequence of least‐squares problems of the form min_y∥G^1/2(A^T y?h)∥₂, where G is an n×n positive‐definite diagonal weight matrix, and A an m×n (m?n) sparse matrix with some dense columns; has many applications in linear programming, electrical networks, elliptic boundary value problems, and structural analysis. We suggest low‐rank correction preconditioners for such problems, and a mixed solver (a combination of a direct solver and an iterative solver). The numerical results show that our technique for selecting the low‐rank correction matrix is very effective. Copyright © 2002 John Wiley & Sons, Ltd.     

One application of Floquet's theory to Lp–Lq estimates for hyperbolic equations with very fast oscillations

In this paper we consider the strictly hyperbolic equation u_tt−λ²(t)b²(t)Δu=0. The coefficient consists of an increasing function λ=λ(t) and a non‐constant periodic function b=b(t). We study the question for the influence of these parts on L_p–L_q decay estimates for the solution of the Cauchy problem. A fairly wide class of equations will be described for which the influence of the oscillating part dominates. This implies, on the one hand, that there exist no L_p–L_q decay estimates and, on the other hand, that the energy estimate from Gronwall's inequality is near to an optimal one. Copyright © 1999 John Wiley & Sons, Ltd.     

A result concerning nilpotent values with generalized skew derivations on Lie ideals

ABSTRACT

Let n≥1 be a fixed integer, R a prime ring with its right Martindale quotient ring Q, C the extended centroid, and L a non-central Lie ideal of R. If F is a generalized skew derivation of R such that (F(x)F(y)?yx)ⁿ = 0 for all x,y∈L, then char(R) = 2 and R?M₂(C), the ring of 2×2 matrices over C.     

Existence of global solutions with prescribed asymptotic behavior for nonlinear ordinary differential equations

Summary Conditions are given for the nonlinear differential equation (1)L _n y+f(t, y, ..., ...,y ^(n–1)=0to have solutions which exist on a given interval [t₀, )and behave in some sense like specified solutions of the linear equation (2)L _n z=0as t.The global nature of these results is unusual as compared to most theorems of this kind, which guarantee the existence of solutions of (1)only for sufficiently large t. The main theorem requires no assumptions regarding oscillation or nonoscillation of solutions of (2).A second theorem is specifically applicable to the situation where (2)is disconjugate on [t ₀, ),and a corollary of the latter applies to the case where Lz=z ⁿ.     

A local normal form theorem for infinitary logic with unary quantifiers

We prove a local normal form theorem of the Gaifman type for the infinitary logic L_∞ω( Q _u)^ω whose formulas involve arbitrary unary quantifiers but finite quantifier rank. We use a local Ehrenfeucht‐Fraïssé type game similar to the one in [9]. A consequence is that every sentence of L_∞ω( Q _u)^ω of quantifier rank n is equivalent to an infinite Boolean combination of sentences of the form (?^≥iy)ψ(y), where ψ(y) has counting quantifiers restricted to the (2^n–1 – 1)‐neighborhood of y. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Regularized traces of higher-order singular differential operators

We consider singular differential operators of order 2m, m ∈ ?, with discrete spectrum in L ₂[0, + ∞). For self-adjoint extensions given by the boundary conditions y(0) = y″(0) = ? = y ^(2m?2)(0) = 0 or y′(0) = y?(0) = ? = y ^(2m?1)(0) = 0, we obtain regularized traces. We present the explicit form of the spectral function, which can be used for calculating regularized traces.     

Best Approximation of Functions like |x| exp(−A|x|)

We determine the exact order of best approximation by polynomials and entire functions of exponential type of functions like?_λ, α(x)=|x|^λ exp(−A|x|^−α). In particular, it is shown thatE(?_λ, α, _n, L_p(−1, 1))∼n^{−(2λp+αp+2)/2p(1+α)}×exp(−(1+α⁻¹)(Aα)^1/(1+α) cos απ/2(1+α) n^α/(1+α)), whereE(?_λ, α, _n, L_p(−1, 1)) denotes best polynomial approximation of?_λ, αinL_p(−1, 1),λ∈,α∈(0, 2],A>0, 1?p?∞. The problem, concerning the exact order of decrease ofE(?_0, 2, _n, L_∞(−1, 1)), has been posed by S. N. Bernstein.     

A new quadratic nonconforming finite element on rectangles

A new quadratic nonconforming finite element on rectangles (or parallelograms) is introduced. The nonconforming element consists of P₂ ⊕ Span{x²y,xy²} on a rectangle and eight degrees of freedom. Our element is essentially of seven degrees of freedom since the degree of freedom associated with the integration on rectangle is essentially of bubble‐function nature. Global basis functions are constructed for both Dirichlet and Neumann type of problems; accordingly the corresponding dimensions are counted. The local and global interpolation operators are defined. Error estimates of optimal order are derived in both broken energy and L²(Ω) norms for second‐order of elliptic problems. Brief numerical results are also shown to confirm the optimality of the presented quadratic nonconforming element. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

The scaling window of the 2‐SAT transition

We consider the random 2‐satisfiability (2‐SAT) problem, in which each instance is a formula that is the conjunction of m clauses of the form x∨y, chosen uniformly at random from among all 2‐clauses on n Boolean variables and their negations. As m and n tend to infinity in the ratio m/n→α, the problem is known to have a phase transition at α_c=1, below which the probability that the formula is satisfiable tends to one and above which it tends to zero. We determine the finite‐size scaling about this transition, namely the scaling of the maximal window W(n, δ)=(α_?(n,δ), α₊(n,δ)) such that the probability of satisfiability is greater than 1?δ for α<α_? and is less than δ for α>α₊. We show that W(n,δ)=(1?Θ(n^?1/3), 1+Θ(n^?1/3)), where the constants implicit in Θ depend on δ. We also determine the rates at which the probability of satisfiability approaches one and zero at the boundaries of the window. Namely, for m=(1+ε)n, where ε may depend on n as long as |ε| is sufficiently small and |ε|n^1/3 is sufficiently large, we show that the probability of satisfiability decays like exp(?Θ(nε³)) above the window, and goes to one like 1?Θ(n^?1|ε|^?3 below the window. We prove these results by defining an order parameter for the transition and establishing its scaling behavior in n both inside and outside the window. Using this order parameter, we prove that the 2‐SAT phase transition is continuous with an order parameter critical exponent of 1. We also determine the values of two other critical exponents, showing that the exponents of 2‐SAT are identical to those of the random graph. © 2001 John Wiley & Sons, Inc. Random Struct. Alg., 18: 201–256 2001     

Elliptic boundary value problems with λ-dependent boundary conditions

In this paper second order elliptic boundary value problems on bounded domains Ω⊂Rn with boundary conditions on ∂Ω depending nonlinearly on the spectral parameter are investigated in an operator theoretic framework. For a general class of locally meromorphic functions in the boundary condition a solution operator of the boundary value problem is constructed with the help of a linearization procedure. In the special case of rational Nevanlinna or Riesz-Herglotz functions on the boundary the solution operator is obtained in an explicit form in the product Hilbert space L²(Ω)⊕(L²m(∂Ω)), which is a natural generalization of known results on λ-linear elliptic boundary value problems and λ-rational boundary value problems for ordinary second order differential equations.     

Construction of positive definite cubature formulae and approximation of functions via Voronoi tessellations

Let Ω ⊂ ℝ ^d be a compact convex set of positive measure. In a recent paper, we established a definiteness theory for cubature formulae of order two on Ω. Here we study extremal properties of those positive definite formulae that can be generated by a centroidal Voronoi tessellation of Ω. In this connection we come across a class of operators of the form L_n[f](x): = ?_i=1ⁿ f_i(x)(f(y_i) + á?f(y_i), x-y_i?)L_n[f](\boldsymbol{x}):= \sum_{i=1}^n \phi_i(\boldsymbol{x})(f(\boldsymbol{y}_i) + \langle\nabla f(\boldsymbol{y}_i), \boldsymbol{x}-\boldsymbol{y}_i\rangle), where y₁,..., y_n\boldsymbol{y}_1,\dots, \boldsymbol{y}_n are distinct points in Ω and {ϕ ₁, ..., ϕ _n } is a partition of unity on Ω. We present best possible pointwise error estimates and describe operators L _n with a smallest constant in an L ^p error estimate for 1 ≤ p < ∞ . For a generalization, we introduce a new type of Voronoi tessellation in terms of a twice continuously differentiable and strictly convex function f. It allows us to describe a best operator L _n for approximating f by L _n [f] with respect to the L ^p norm.     

Estimates for the Hill operator, II

Consider the Hill operator Ty=-y^″+q(t)y in L²(R), where the real potential q is 1-periodic and q,q^′∈L²(0,1). The spectrum of T consists of spectral bands separated by gaps γ_n,n?1 with length |γ_n|?0. We obtain two-sided estimates of the gap lengths ∑n²|γ_n|² in terms of . Moreover, we obtain the similar two-sided estimates for spectral data (the height of the corresponding slit on the quasimomentum domain, action variables for the KdV equation and so on). In order prove this result we use the analysis of a conformal mapping corresponding to quasimomentum of the Hill operator. That makes it possible to reformulate the problems for the differential operator as the problems of the conformal mapping theory. Then the proof is based on the analysis of the conformal mapping, the embedding theorems and the identities. Furthermore, we obtain the similar two-sided estimates for potentials which have p?2 derivatives.     

Varadhan–Léandre estimates for a family of random vectors

We deal with the density p ^? (y) of a random vector's family with the following Wiener–Chaos decomposition F ^? =∑ _n=0 ^∞ ? ⁿ I _n (f _n ). We first give several results of regularity for the Wiener functional F ^? . Secondly, we analyze the asymptotic behavior of the density p ^? (y). Exactly, we establish the Varadhan–Léandre estimates using Malliavin Calculus. Finally this general result is applied to the study of the density in one example of hyperbolic stochastic partial differential equation.

     

Orthogonal Polynomial Solutions of Spectral Type Differential Equations: Magnus' Conjecture

Let τ=σ+ν be a point mass perturbation of a classical moment functional σ by a distribution ν with finite support. We find necessary conditions for the polynomials {Q_n(x)}^∞_n=0, orthogonal relative to τ, to be a Bochner–Krall orthogonal polynomial system (BKOPS); that is, {Q_n(x)}^∞_n=0 are eigenfunctions of a finite order linear differential operator of spectral type with polynomial coefficients: L_N[y](x)=∑^N_i=1 ℓ_i(x) y⁽ⁱ⁾(x)=λ_ny(x). In particular, when ν is of order 0 as a distribution, we find necessary and sufficient conditions for {Q_n(x)}^∞_n=0 to be a BKOPS, which strongly support and clarify Magnus' conjecture which states that any BKOPS must be orthogonal relative to a classical moment functional plus one or two point masses at the end point(s) of the interval of orthogonality. This result explains not only why the Bessel-type orthogonal polynomials (found by Hendriksen) cannot be a BKOPS but also explains the phenomena for infinite-order differential equations (found by J. Koekoek and R. Koekoek), which have the generalized Jacobi polynomials and the generalized Laguerre polynomials as eigenfunctions.