BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Christoph Koutschan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
Reconstructing a hypothetical recurrence equation from the first terms of an infinite sequence is a classical and well-known technique in experimental mathematics. We propose a variation of this technique which can succeed with fewer input terms.
A minimally rigid graph, also called Laman graph, models a planar framework which is rigid for a general choice of distances between its vertices. In other words, there are finitely many ways, up to isometries, to realize such a graph in the plane. Using ideas from algebraic and tropical geometry, we derive a recursive formula for the number of such realizations. Combining computational results with the construction of new rigid graphs via gluing techniques, we can give a new lower bound on the maximal possible number of realizations for graphs with a given number of vertices.
For an (irreducible) recurrence equation with coefficients from $\mathbb Z[n]$ and its two linearly independent rational solutions $u_n,v_n$, the limit of $u_n/v_n$ as $n\to\infty$, when exists, is called the Apéry limit. We give a construction that realises certain quotients of $L$-values of elliptic curves as Apéry limits.
In this memorial tribute to Joe Gillis, who taught us that Special Functions count, we show how the seminal Even-Gillis integral formula for the number of derangements of a multiset, in terms of Laguerre polynomials, can be used to efficiently compute not only the number of the title, but much harder ones, when it is interfaced with Wilf-Zeilberger algorithmic proof theory.
There are only aleph-zero rational numbers, while there are 2 to the power aleph-zero real numbers. Hence the probability that a randomly chosen real number would be rational is 0. Yet proving rigorously that any specific, natural, real constant, is irrational is usually very hard, witness that there are still no proofs of the irrationality of the Euler-Mascheroni constant, the Catalan constant, or Zeta(5). Inspired by Frits Beukers' elegant rendition of Apery's seminal proofs of the irrationality of Zeta(2) and Zeta(3), and heavily using algorithmic proof theory, we systematically searched for other similar integrals, that lead to irrationality proofs. We found quite a few candidates for such proofs, including the square-root of Pi times Gamma(7/3)/Gamma(-1/6) and Gamma(19/6)/Gamma(8/3) divided by the square-root of Pi.
In theory and practice of inverse problems, linear operator equations $Tx=y$ with compact linear forward operators $T$ having a non-closed range $\mathcal{R}(T)$ and mapping between infinite dimensional Hilbert spaces plays some prominent role. As a consequence of the ill-posedness of such problems, regularization approaches are required, and due to its unlimited qualification spectral cut-off is an appropriate method for the stable approximate solution of corresponding inverse problems. For this method, however, the singular system $\{σ_i(T),u_i(T),v_i(T)\}_{i=1}^\infty$ of the compact operator $T$ is needed, at least for $i=1,2,...,N$, up to some stopping index $N$. In this note we consider $n$-fold integration operators $T=J^n\;(n=1,2,...)$ in $L^2([0,1])$ occurring in numerous applications, where the solution of the associated operator equation is characterized by the $n$-th generalized derivative $x=y^{(n)}$ of the Sobolev space function $y \in H^n([0,1])$. Almost all textbooks on linear inverse problems present the whole singular system $\{σ_i(J^1),u_i(J^1),v_i(J^1)\}_{i=1}^\infty$ in an explicit manner. However, they do not discuss the singular systems for $J^n,\;n \ge 2$. We will emphasize that this seems to be a consequence of the fact that for higher $n$ the eigenvalues $σ^2_i(J^n)$ of the associated ODE boundary value problems obey transcendental equations, the complexity of which is growing with $n$. We present the transcendental equations for $n=2,3,...$ and discuss and illustrate the associated eigenfunctions and some of their properties.
We study zero-forcing detection (ZF) for multiple-input/multiple-output (MIMO) spatial multiplexing under transmit-correlated Rician fading for an N_R X N_T channel matrix with rank-1 line-of-sight (LoS) component. By using matrix transformations and multivariate statistics, our exact analysis yields the signal-to-noise ratio moment generating function (m.g.f.) as an infinite series of gamma distribution m.g.f.'s and analogous series for ZF performance measures, e.g., outage probability and ergodic capacity. However, their numerical convergence is inherently problematic with increasing Rician K-factor, N_R , and N_T. We circumvent this limitation as follows. First, we derive differential equations satisfied by the performance measures with a novel automated approach employing a computer-algebra tool which implements Groebner basis computation and creative telescoping. These differential equations are then solved with the holonomic gradient method (HGM) from initial conditions computed with the infinite series. We demonstrate that HGM yields more reliable performance evaluation than by infinite series alone and more expeditious than by simulation, for realistic values of K , and even for N_R and N_T relevant to large MIMO systems. We envision extending the proposed approaches for exact analysis and reliable evaluation to more general Rician fading and other transceiver methods.
With the exception of q-hypergeometric summation, the use of computer algebra packages implementing Zeilberger's "holonomic systems approach" in a broader mathematical sense is less common in the field of q-series and basic hypergeometric functions. A major objective of this article is to popularize the usage of such tools also in these domains. Concrete case studies showing software in action introduce to the basic techniques. An application highlight is a new computer-assisted proof of the celebrated Ismail-Zhang formula, an important q-analog of a classical expansion formula of plane waves in terms of Gegenbauer polynomials.
In the convergence analysis of numerical methods for solving partial differential equations (such as finite element methods) one arrives at certain generalized eigenvalue problems, whose maximal eigenvalues need to be estimated as accurately as possible. We apply symbolic computation methods to the situation of square elements and are able to improve the previously known upper bound, given in "p- and hp-finite element methods" (Schwab, 1998), by a factor of 8. More precisely, we try to evaluate the corresponding determinant using the holonomic ansatz, which is a powerful tool for dealing with determinants, proposed by Zeilberger in 2007. However, it turns out that this method does not succeed on the problem at hand. As a solution we present a variation of the original holonomic ansatz that is applicable to a larger class of determinants, including the one we are dealing with here. We obtain an explicit closed form for the determinant, whose special form enables us to derive new and tight upper resp. lower bounds on the maximal eigenvalue, as well as its asymptotic behaviour.
Designing mechanical devices, called linkages, that draw a given plane curve has been a topic that interested engineers and mathematicians for hundreds of years, and recently also computer scientists. Already in 1876, Kempe proposed a procedure for solving the problem in full generality, but his constructions tend to be extremely complicated. We provide a novel algorithm that produces much simpler linkages, but works only for parametric curves. Our approach is to transform the problem into a factorization task over some noncommutative algebra. We show how to compute such a factorization, and how to use it to construct a linkage tracing a given curve.
Continuing a series of articles in the past few years on creative telescoping using reductions, we develop a new algorithm to construct minimal telescopers for algebraic functions. This algorithm is based on Trager's Hermite reduction and on polynomial reduction, which was originally designed for hyperexponential functions and extended to the algebraic case in this paper.
Continuing a series of articles in the past few years on creative telescoping using reductions, we adapt Trager's Hermite reduction for algebraic functions to fuchsian D-finite functions and develop a reduction-based creative telescoping algorithm for this class of functions, thereby generalizing our recent reduction-based algorithm for algebraic functions, presented at ISSAC 2016.
Transportation processes, which play a prominent role in the life and social sciences, are typically described by discrete models on lattices. For studying their dynamics a continuous formulation of the problem via partial differential equations (PDE) is employed. In this paper we propose a symbolic computation approach to derive mean-field PDEs from a lattice-based model. We start with the microscopic equations, which state the probability to find a particle at a given lattice site. Then the PDEs are formally derived by Taylor expansions of the probability densities and by passing to an appropriate limit as the time steps and the distances between lattice sites tend to zero. We present an implementation in a computer algebra system that performs this transition for a general class of models. In order to rewrite the mean-field PDEs in a conservative formulation, we adapt and implement symbolic integration methods that can handle unspecified functions in several variables. To illustrate our approach, we consider an application in crowd motion analysis where the dynamics of bidirectional flows are studied. However, the presented approach can be applied to various transportation processes of multiple species with variable size in any dimension, for example, to confirm several proposed mean-field models for cell motility.
We study multi-parameter solutions of the inhomogeneous paraxial wave equation in a linear and quadratic approximation which include oscillating laser beams in a parabolic waveguide, spiral light beams, and other important families of propagation-invariant laser modes in weakly varying media. A "smart" lens design and a similar effect of superfocusing of particle beams in a thin monocrystal film are also discussed. In the supplementary electronic material, we provide a computer algebra verification of the results presented here, and of some related mathematical tools that were stated without proofs in the literature. We also demonstrate how computer algebra can be used to derive some of the presented formulas automatically, which is highly desirable as the corresponding hand calculations are very tedious. In numerical simulations, some of the new solutions reveal quite exotic properties which deserve further investigation including an experimental observation.
We propose a differential analog of the notion of integral closure of algebraic function fields. We present an algorithm for computing the integral closure of the algebra defined by a linear differential operator. Our algorithm is a direct analog of van Hoeij's algorithm for computing integral bases of algebraic function fields.
For multiple-input multiple-output (MIMO) spatial-multiplexing transmission, zero-forcing detection (ZF) is appealing because of its low complexity. Our recent MIMO ZF performance analysis for Rician--Rayleigh fading, which is relevant in heterogeneous networks, has yielded for the ZF outage probability and ergodic capacity infinite-series expressions. Because they arose from expanding the confluent hypergeometric function $ {_1\! F_1} (\cdot, \cdot, σ) $ around 0, they do not converge numerically at realistically-high Rician $ K $-factor values. Therefore, herein, we seek to take advantage of the fact that $ {_1\! F_1} (\cdot, \cdot, σ) $ satisfies a differential equation, i.e., it is a \textit{holonomic} function. Holonomic functions can be computed by the \textit{holonomic gradient method} (HGM), i.e., by numerically solving the satisfied differential equation. Thus, we first reveal that the moment generating function (m.g.f.) and probability density function (p.d.f.) of the ZF signal-to-noise ratio (SNR) are holonomic. Then, from the differential equation for $ {_1\! F_1} (\cdot, \cdot, σ) $, we deduce those satisfied by the SNR m.g.f. and p.d.f., and demonstrate that the HGM helps compute the p.d.f. accurately at practically-relevant values of $ K $. Finally, numerical integration of the SNR p.d.f. produced by HGM yields accurate ZF outage probability and ergodic capacity results.
The Apagodu-Zeilberger algorithm can be used for computing annihilating operators for definite sums over hypergeometric terms, or for definite integrals over hyperexponential functions. In this paper, we propose a generalization of this algorithm which is applicable to arbitrary $\partial$-finite functions. In analogy to the hypergeometric case, we introduce the notion of proper $\partial$-finite functions. We show that the algorithm always succeeds for these functions, and we give a tight a priori bound for the order of the output operator.
We derive the recurrence relations for relativistic Coulomb integrals directly from the integral representations with the help of computer algebra methods. In order to manage the computational complexity of this problem, we employ holonomic closure properties in a sophisticated way.
We prove three conjectures concerning the evaluation of determinants, which are related to the counting of plane partitions and rhombus tilings. One of them was posed by George Andrews in 1980, the other two were by Guoce Xin and Christian Krattenthaler. Our proofs employ computer algebra methods, namely, the holonomic ansatz proposed by Doron Zeilberger and variations thereof. These variations make Zeilberger's original approach even more powerful and allow for addressing a wider variety of determinants. Finally, we present, as a challenge problem, a conjecture about a closed-form evaluation of Andrews's determinant.
We employ computer algebra algorithms to prove a collection of identities involving Bessel functions with half-integer orders and other special functions. These identities appear in the famous Handbook of Mathematical Functions, as well as in its successor, the DLMF, but their proofs were lost. We use generating functions and symbolic summation techniques to produce new proofs for them.
The aim of this article is twofold: on the one hand it is intended to serve as a gentle introduction to the topic of creative telescoping, from a practical point of view; for this purpose its application to several problems is exemplified. On the other hand, this chapter has the flavour of a survey article: the developments in this area during the last two decades are sketched and a selection of references is compiled in order to highlight the impact of creative telescoping in numerous contexts.
Our goal is to compute the minimal-order recurrence of the colored Jones polynomial of the 7_4 knot, as well as for the first four double twist knots. As a corollary, we verify the AJ Conjecture for the simplest knot 7_4 with reducible non-abelian SL(2,C) character variety. To achieve our goal, we use symbolic summation techniques of Zeilberger's holonomic systems approach and an irreducibility criterion for q-difference operators. For the latter we use an improved version of the qHyper algorithm of Abramov-Paule-Petkovsek to show that a given q-difference operator has no linear right factors. En route, we introduce exterior power Adams operations on the ring of bivariate polynomials and on the corresponding affine curves.
We study the face-centered cubic lattice (fcc) in up to six dimensions. In particular, we are concerned with lattice Green's functions (LGF) and return probabilities. Computer algebra techniques, such as the method of creative telescoping, are used for deriving an ODE for a given LGF. For the four- and five-dimensional fcc lattices, we give rigorous proofs of the ODEs that were conjectured by Guttmann and Broadhurst. Additionally, we find the ODE of the LGF of the six-dimensional fcc lattice, a result that was not believed to be achievable with current computer hardware.
We consider the numerical discretization of the time-domain Maxwell's equations with an energy-conserving discontinuous Galerkin finite element formulation. This particular formulation allows for higher order approximations of the electric and magnetic field. Special emphasis is placed on an efficient implementation which is achieved by taking advantage of recurrence properties and the tensor-product structure of the chosen shape functions. These recurrences have been derived symbolically with computer algebra methods reminiscent of the holonomic systems approach.
The conjecture that the orbit-counting generating function for totally symmetric plane partitions can be written as an explicit product formula, has been stated independently by George Andrews and David Robbins around 1983. We present a proof of this long-standing conjecture.
We study q-holonomic sequences that arise as the colored Jones polynomial of knots in 3-space. The minimal-order recurrence for such a sequence is called the (non-commutative) A-polynomial of a knot. Using the "method of guessing", we obtain this polynomial explicitly for the K_p = (-2, 3, 3+2p) pretzel knots for p = -5, ..., 5. This is a particularly interesting family since the pairs (K_p, -K_{-p}) are geometrically similar (in particular, scissors congruent) with similar character varieties. Our computation of the non-commutative A-polynomial (a) complements the computation of the A-polynomial of the pretzel knots done by the first author and Mattman, (b) supports the AJ Conjecture for knots with reducible A-polynomial and (c) numerically computes the Kashaev invariant of pretzel knots in linear time. In a later publication, we will use the numerical computation of the Kashaev invariant to numerically verify the Volume Conjecture for the above mentioned pretzel knots.
We prove an integrability criterion of order 3 for a homogeneous potential of degree -1 in the plane. Still, this criterion depends on some integer and it is impossible to apply it directly except for families of potentials whose eigenvalues are bounded. To address this issue, we use holonomic and asymptotic computations with error control of this criterion and apply it to the potential of the form V(r,θ)=r^{-1} h(\exp(iθ)) with h a polynomial of degree less than 3. We find then all meromorphically integrable potentials of this form.
A sequence $f_n(q)$ is $q$-holonomic if it satisfies a nontrivial linear recurrence with coefficients polynomials in $q$ and $q^n$. Our main theorems state that $q$-holonomicity is preserved under twisting, i.e., replacing $q$ by $ωq$ where $ω$ is a complex root of unity, and under the substitution $q \to q^α$ where $α$ is a rational number. Our proofs are constructive, work in the multivariate setting of $\partial$-finite sequences and are implemented in the Mathematica package HolonomicFunctions. Our results are illustrated by twisting natural $q$-holonomic sequences which appear in quantum topology, namely the colored Jones polynomial of pretzel knots and twist knots. The recurrence of the twisted colored Jones polynomial can be used to compute the asymptotics of the Kashaev invariant of a knot at an arbitrary complex root of unity.
A variation of Zeilberger's holonomic ansatz for symbolic determinant evaluations is proposed which is tailored to deal with Pfaffians. The method is also applicable to determinants of skew-symmetric matrices, for which the original approach does not work. As Zeilberger's approach is based on the Laplace expansion (cofactor expansion) of the determinant, we derive our approach from the cofactor expansion of the Pfaffian. To demonstrate the power of our method, we prove, using computer algebra algorithms, some conjectures proposed in the paper "Pfaffian decomposition and a Pfaffian analogue of q-Catalan Hankel determinants" by Ishikawa, Tagawa, and Zeng. A minor summation formula related to partitions and Motzkin paths follows as a corollary.
In this note we reinvestigate the task of computing creative telescoping relations in differential-difference operator algebras. Our approach is based on an ansatz that explicitly includes the denominators of the delta parts. We contribute several ideas of how to make an implementation of this approach reasonably fast and provide such an implementation. A selection of examples shows that it can be superior to existing methods by a large factor.
In computer algebra there are different ways of approaching the mathematical concept of functions, one of which is by defining them as solutions of differential equations. We compare different such approaches and discuss the occurring problems. The main focus is on the question of determining possible branch cuts. We explore the extent to which the treatment of branch cuts can be rendered (more) algorithmic, by adapting Kahan's rules to the differential equation setting.
We consider algebras over a field K, generated by two variables x and y subject to the single relation yx = qxy + ax + by + c for q in K^* and a, b, c in K. We prove, that among such algebras there are precisely five isomorphism classes. The representatives of these classes, which are ubiquitous operator algebras, are called model algebras. We derive explicit multiplication formulas for y^m*x^n in terms of standard monomials x^i*y^j for many algebras of the considered type. Such formulas are used in establishing formulas of binomial type and in implementing non-commutative multiplication in a computer algebra system. By using the formulas we also study centers and ring-theoretic properties of the non-commutative model algebras.
In order to appreciate how well off we mathematicians and scientists are today, with extremely fast hardware and lots and lots of memory, as well as with powerful software, both for numeric and symbolic computation, it may be a good idea to go back to the early days of electronic computers and compare how things went then. We have chosen, as a case study, a problem that was considered a huge challenge at the time. Namely, we looked at C.L. Pekeris's seminal 1958 work on the ground state energies of two-electron atoms. We went through all the computations ab initio with today's software and hardware, with a special emphasis on the symbolic computations which in 1958 had to be made by hand, and which nowadays can be automated and generalized.
For a polynomial P, we consider the sequence of iterated integrals of ln P(x). This sequence is expressed in terms of the zeros of P(x). In the special case of ln(1 + x^2), arithmetic properties of certain coefficients arising are described. Similar observations are made for ln(1 + x^3).
The SL_3 colored Jones polynomial of the trefoil knot is a $q$-holonomic sequence of two variables with natural origin, namely quantum topology. The paper presents an explicit set of generators for the annihilator ideal of this $q$-holonomic sequence as a case study. On the one hand, our results are new and useful to quantum topology: this is the first example of a rank 2 Lie algebra computation concerning the colored Jones polynomial of a knot. On the other hand, this work illustrates the applicability and computational power of the employed computer algebra methods.
We present a computer-aided, yet fully rigorous, proof of Ira Gessel's tantalizingly simply-stated conjecture that the number of ways of walking $2n$ steps in the region $x+y \geq 0, y \geq 0$ of the square-lattice with unit steps in the east, west, north, and south directions, that start and end at the origin, equals $16^n\frac{(5/6)_n(1/2)_n}{(5/3)_n(2)_n}$ .