Formal orthogonal polynomials and Hankel/Toeplitz duality

For classical polynomials orthogonal with respect to a positive measure supported on the real line, the moment matrix is Hankel and positive definite. The polynomials satisfy a three term recurrence relation. When the measure is supported on the complex unit circle, the moment matrix is positive definite and Toeplitz. Then they satisfy a coupled Szeg recurrence relation but also a three term recurrence relation. In this paper we study the generalization for formal polynomials orthogonal with respect to an arbitrary moment matrix and consider arbitrary Hankel and Toeplitz matrices as special cases. The relation with Padé approximation and with Krylov subspace iterative methods is also outlined.This research was supported by the National Fund for Scientific Research (NFWO), project Lanczos, grant #2.0042.93.     

On the convergence of multipoint Padé-type approximants and quadrature formulas associated with the unit circle

We study the convergence of rational interpolants with prescribed poles on the unit circle to the Herglotz-Riesz transform of a complex measure supported on [–, ]. As a consequence, quadrature formulas arise which integrate exactly certain rational functions. Estimates of the rate of convergence of these quadrature formulas are also included.This research was performed as part of the European project ROLLS under contract CHRX-CT93-0416.     

A moment problem associated to rational Szegő functions

Leta ₁,...,a _p be distinct points in the finite complex plane ?, such that |a _j|>1,j=1,..., p and let \(b_j = 1/\bar \alpha _j ,\) j=1,..., p. Let μ₀, μ _π ^(j) , ν _π ^(j) j=1,..., p;n=1, 2,... be given complex numbers. We consider the following moment problem. Find a distribution ψ on [?π, π], with infinitely many points of increase, such that $$\begin{array}{l} \int_{ - \pi }^\pi {d\psi (\theta ) = \mu _0 ,} \\ \int_{ - \pi }^\pi {\frac{{d\psi (\theta )}}{{(e^{i\theta } - a_j )^n }} = \mu _n^{(j)} ,} \int_{ - \pi }^\pi {\frac{{d\psi (\theta )}}{{(e^{i\theta } - b_j )^n }} = v_n^{(j)} ,} j = 1,...,p;n = 1,2,.... \\ \end{array}$$ It will be shown that this problem has a unique solution if the moments generate a positive-definite Hermitian inner product on the linear space of rational functions with no poles in the extended complex plane ?^* outside {a ₁,...,a _p,b ₁,...,b _p}.     

Rational wavelets on the real line

In this note a representation of the discrete Green's function of a compact discretization of a two point boundary value problem of order n 2 is given which among other things allows a direct proff of the convergence (and divergence) properties.     

Power law statistics of rippled graphene nanoflakes

Graphene nanoflakes (GNFs) are predicted to possess novel magnetic, optical, and spintronic properties. They have recently been synthesized and a number of applications are being studied. Here we investigate the statistical properties of rippled GNFs (50–5,000 atoms) at $\text{ T}=300$  K. An adjacency matrix is calculated from the coordinates and we find that the free energy, enthalpy, entropy, and atomic displacement all show power law behavior. The vibrational energy versus the Wiener index also shows power law character. We distinguish between using Euclidean topographical indices and compare them to topological ones. These properties are determined from atomic coordinates using MATLAB routines.     

Recursive relations for block Hankel and Toeplitz systems part I: Direct recursions

A set of formulas is given for the relations that exist between the first and last block now or column of the inverse of a block Hankel or Toeplitz matrix. This is related to making arbitrary steps in a matrix Padé table.     

On computing rational Gauss-Chebyshev quadrature formulas

We provide an algorithm to compute the nodes and weights for Gauss-Chebyshev quadrature formulas integrating exactly in spaces of rational functions with arbitrary real poles outside . Contrary to existing rational quadrature formulas, the computational effort is very low, even for extremely high degrees, and under certain conditions on the poles it can be shown that the complexity is of order . This method is based on the derivation of explicit expressions for Chebyshev orthogonal rational functions, which are (thus far) the only examples of explicitly known orthogonal rational functions on with arbitrary real poles outside this interval.