Computing the arithmetic genus of Hilbert modular fourfolds

The Hilbert modular fourfold determined by the totally real quartic number field is a desingularization of a natural compactification of the quotient space , where PSL acts on by fractional linear transformations via the four embeddings of into . The arithmetic genus, equal to one plus the dimension of the space of Hilbert modular cusp forms of weight , is a birational invariant useful in the classification of these varieties. In this work, we describe an algorithm allowing for the automated computation of the arithmetic genus and give sample results.

     

Optimal logarithmic energy points on the unit sphere

We study minimum energy point charges on the unit sphere in , , that interact according to the logarithmic potential , where is the Euclidean distance between points. Such optimal -point configurations are uniformly distributed as . We quantify this result by estimating the spherical cap discrepancy of optimal energy configurations. The estimate is of order . Essential is an improvement of the lower bound of the optimal logarithmic energy which yields the second term in the asymptotical expansion of the optimal energy. Previously, this was known for the unit sphere in only. Furthermore, we present an upper bound for the error of integration for an equally-weighted numerical integration rule with the nodes forming an optimal logarithmic energy configuration. For polynomials of degree at most this bound is as . For continuous functions of satisfying a Lipschitz condition with constant the bound is as .

     

Comparison theorems of Kolmogorov type and exact values of -widths on Hardy-Sobolev classes

Let be a strip in complex plane. denotes those -periodic, real-valued functions on which are analytic in the strip and satisfy the condition , . Osipenko and Wilderotter obtained the exact values of the Kolmogorov, linear, Gel'fand, and information -widths of in , , and 2-widths of in , , .

In this paper we continue their work. Firstly, we establish a comparison theorem of Kolmogorov type on , from which we get an inequality of Landau-Kolmogorov type. Secondly, we apply these results to determine the exact values of the Gel'fand -width of in , . Finally, we calculate the exact values of Kolmogorov -width, linear -width, and information -width of in , , .

     

Even moments of generalized Rudin-Shapiro polynomials

We know from Littlewood (1968) that the moments of order of the classical Rudin-Shapiro polynomials satisfy a linear recurrence of degree . In a previous article, we developed a new approach, which enables us to compute exactly all the moments of even order for . We were also able to check a conjecture on the asymptotic behavior of , namely , where , for even and . Now for every integer there exists a sequence of generalized Rudin-Shapiro polynomials, denoted by . In this paper, we extend our earlier method to these polynomials. In particular, the moments have been completely determined for and , for and and for and . For higher values of and , we formulate a natural conjecture, which implies that , where is an explicit constant.

     

Data-sparse approximation to the operator-valued functions of elliptic operator

In previous papers the arithmetic of hierarchical matrices has been described, which allows us to compute the inverse, for instance, of finite element stiffness matrices discretising an elliptic operator The required computing time is up to logarithmic factors linear in the dimension of the matrix. In particular, this technique can be used for the computation of the discrete analogue of a resolvent

In the present paper, we consider various operator functions, the operator exponential negative fractional powers , the cosine operator function and, finally, the solution operator of the Lyapunov equation. Using the Dunford-Cauchy representation, we get integrals which can be discretised by a quadrature formula which involves the resolvents mentioned above. We give error estimates which are partly exponentially, partly polynomially decreasing.

     

Bounds on differences of adjacent zeros of Bessel functions and iterative relations between consecutive zeros

Bounds for the distance between adjacent zeros of cylinder functions are given; and are such that ; stands for the th positive zero of the cylinder (Bessel) function , , .

These bounds, together with the application of modified (global) Newton methods based on the monotonic functions and , give rise to forward ( ) and backward ( ) iterative relations between consecutive zeros of cylinder functions.

The problem of finding all the positive real zeros of Bessel functions for any real and inside an interval , 0$">, is solved in a simple way.

     

Weak coupling of solutions of first-order least-squares method

A theoretical analysis of a first-order least-squares finite element method for second-order self-adjoint elliptic problems is presented. We investigate the coupling effect of the approximate solutions for the primary function and for the flux . We prove that the accuracy of the approximate solution for the primary function is weakly affected by the flux . That is, the bound for is dependent on , but only through the best approximation for multiplied by a factor of meshsize . Similarly, we provide that the bound for is dependent on , but only through the best approximation for multiplied by a factor of the meshsize . This weak coupling is not true for the non-selfadjoint case. We provide the numerical experiment supporting the theorems in this paper.

     

Efficient CM-constructions of elliptic curves over finite fields

We present an algorithm that, on input of an integer together with its prime factorization, constructs a finite field and an elliptic curve over for which has order . Although it is unproved that this can be done for all , a heuristic analysis shows that the algorithm has an expected run time that is polynomial in , where is the number of distinct prime factors of . In the cryptographically relevant case where is prime, an expected run time can be achieved. We illustrate the efficiency of the algorithm by constructing elliptic curves with point groups of order and nextprime.

     

Strong tractability of multivariate integration using quasi-Monte Carlo algorithms

We study quasi-Monte Carlo algorithms based on low discrepancy sequences for multivariate integration. We consider the problem of how the minimal number of function evaluations needed to reduce the worst-case error from its initial error by a factor of depends on and the dimension . Strong tractability means that it does not depend on and is bounded by a polynomial in . The least possible value of the power of is called the -exponent of strong tractability. Sloan and Wozniakowski established a necessary and sufficient condition of strong tractability in weighted Sobolev spaces, and showed that the -exponent of strong tractability is between 1 and 2. However, their proof is not constructive.

In this paper we prove in a constructive way that multivariate integration in some weighted Sobolev spaces is strongly tractable with -exponent equal to 1, which is the best possible value under a stronger assumption than Sloan and Wozniakowski's assumption. We show that quasi-Monte Carlo algorithms using Niederreiter's -sequences and Sobol sequences achieve the optimal convergence order for any 0$"> independent of the dimension with a worst case deterministic guarantee (where is the number of function evaluations). This implies that strong tractability with the best -exponent can be achieved in appropriate weighted Sobolev spaces by using Niederreiter's -sequences and Sobol sequences.

     

On fundamental domains of arithmetic Fuchsian groups

Let be a totally real algebraic number field and an order in a quaternion algebra over . Assume that the group of units in with reduced norm equal to is embedded into as an arithmetic Fuchsian group. It is shown how Ford's algorithm can be effectively applied in order to determine a fundamental domain of as well as a complete system of generators of .

     

Lower bounds for the condition number of a real confluent Vandermonde matrix

Lower bounds on the condition number of a real confluent Vandermonde matrix are established in terms of the dimension , or and the largest absolute value among all nodes that define the confluent Vandermonde matrix and the interval that contains the nodes. In particular, it is proved that for any modest (the largest multiplicity of distinct nodes), behaves no smaller than , or than if all nodes are nonnegative. It is not clear whether those bounds are asymptotically sharp for modest .

     

Limiting set of second order spectra

Let be a self-adjoint operator acting on a Hilbert space . A complex number is in the second order spectrum of relative to a finite-dimensional subspace iff the truncation to of is not invertible. This definition was first introduced in Davies, 1998, and according to the results of Levin and Shargorodsky in 2004, these sets provide a method for estimating eigenvalues free from the problems of spectral pollution. In this paper we investigate various aspects related to the issue of approximation using second order spectra. Our main result shows that under fairly mild hypothesis on the uniform limit of these sets, as increases towards , contain the isolated eigenvalues of of finite multiplicity. Therefore, unlike the majority of the standard methods, second order spectra combine nonpollution and approximation at a very high level of generality.

     

An analysis of the Rayleigh-Ritz method for approximating eigenspaces

This paper concerns the Rayleigh-Ritz method for computing an approximation to an eigenspace of a general matrix from a subspace that contains an approximation to . The method produces a pair that purports to approximate a pair , where is a basis for and . In this paper we consider the convergence of as the sine of the angle between and approaches zero. It is shown that under a natural hypothesis--called the uniform separation condition--the Ritz pairs converge to the eigenpair . When one is concerned with eigenvalues and eigenvectors, one can compute certain refined Ritz vectors whose convergence is guaranteed, even when the uniform separation condition is not satisfied. An attractive feature of the analysis is that it does not assume that has distinct eigenvalues or is diagonalizable.

     

A complete description of Golay pairs for lengths up to 100

In his 1961 paper, Marcel Golay showed how the search for pairs of binary sequences of length with complementary autocorrelation is at worst a problem. Andres, in his 1977 master's thesis, developed an algorithm which reduced this to a search and investigated lengths up to 58 for existence of pairs. In this paper, we describe refinements to this algorithm, enabling a search at length 82. We find no new pairs at the outstanding lengths 74 and 82. In extending the theory of composition, we are able to obtain a closed formula for the number of pairs of length generated by a primitive pair of length . Combining this with the results of searches at all allowable lengths up to 100, we identify five primitive pairs. All others pairs of lengths less than 100 may be derived using the methods outlined.

     

On the minimal polynomial of Gauss periods for prime powers

For a positive integer , set and let denote the group of reduced residues modulo . Fix a congruence group of conductor and of order . Choose integers to represent the cosets of in . The Gauss periods

corresponding to are conjugate and distinct over with minimal polynomial

To determine the coefficients of the period polynomial (or equivalently, its reciprocal polynomial is a classical problem dating back to Gauss. Previous work of the author, and Gupta and Zagier, primarily treated the case , an odd prime, with fixed. In this setting, it is known for certain integral power series and , that for any positive integer

holds in for all primes except those in an effectively determinable finite set. Here we describe an analogous result for the case , a prime power ( ). The methods extend for odd prime powers to give a similar result for certain twisted Gauss periods of the form

where denotes the usual Legendre symbol and .

     

Jacobi sums and new families of irreducible polynomials of Gaussian periods

Let 2$">, an -th primitive root of 1, mod a prime number, a primitive root modulo and . We study the Jacobi sums , , where is the least nonnegative integer such that mod . We exhibit a set of properties that characterize these sums, some congruences they satisfy, and a MAPLE program to calculate them. Then we use those results to show how one can construct families , , of irreducible polynomials of Gaussian periods, , of degree , where is a suitable set of primes mod . We exhibit examples of such families for several small values of , and give a MAPLE program to construct more of them.

     

On the embedding problem for representations

Let denote the double cover of corresponding to the element in where transpositions lift to elements of order and the product of two disjoint transpositions to elements of order . Given an elliptic curve , let denote its -torsion points. Under some conditions on elements in correspond to Galois extensions of with Galois group (isomorphic to) . In this work we give an interpretation of the addition law on such fields, and prove that the obstruction for having a Galois extension with gives a homomorphism . As a corollary we can prove (if has conductor divisible by few primes and high rank) the existence of -dimensional representations of the absolute Galois group of attached to and use them in some examples to construct modular forms mapping via the Shimura map to (the modular form of weight attached to) .

     

On strong tractability of weighted multivariate integration

We prove that for every dimension and every number of points, there exists a point-set whose -weighted unanchored discrepancy is bounded from above by independently of provided that the sequence has for some (even arbitrarily large) . Here is a positive number that could be chosen arbitrarily close to zero and depends on but not on or . This result yields strong tractability of the corresponding integration problems including approximation of weighted integrals over unbounded domains such as . It also supplements the results that provide an upper bound of the form when .

     

Hermite and Hermite-Fejér interpolation for Stieltjes polynomials

Let and be the ultraspherical polynomials with respect to . Then we denote by the Stieltjes polynomials with respect to satisfying

In this paper, we show uniform convergence of the Hermite-Fejér interpolation polynomials and based on the zeros of the Stieltjes polynomials and the product for and , respectively. To prove these results, we prove that the Lebesgue constants of Hermite-Fejér interpolation operators for the Stieltjes polynomials and the product are optimal, that is, the Lebesgue constants and have optimal order . In the case of the Hermite-Fejér interpolation polynomials for , we prove weighted uniform convergence. Moreover, we give some convergence theorems of Hermite-Fejér and Hermite interpolation polynomials for in weighted norms.

     

On the polynomial representation for the number of partitions with fixed length

In this paper, it is shown that the number of partitions of a nonnegative integer with parts can be described by a set of polynomials of degree in , where denotes the least common multiple of the integers and denotes the quotient of when divided by . In addition, the sets of the polynomials are obtained and shown explicitly for and .