Researcher profile

Steven J. Miller

Steven J. Miller contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
56works
0followers
15topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

56 published item(s)

preprint2022arXiv

$k$-Diophantine $m$-tuples in Finite Fields

In this paper, we define a $k$-Diophantine $m$-tuple to be a set of $m$ positive integers such that the product of any $k$ distinct positive integers is one less than a perfect square. We study these sets in finite fields $\mathbb{F}_p$ for odd prime $p$ and guarantee the existence of a $k$-Diophantine m-tuple provided $p$ is larger than some explicit lower bound. We also give a formula for the number of 3-Diophantine triples in $\mathbb{F}_p$ as well as an asymptotic formula for the number of $k$-Diophantine $k$-tuples.

preprint2022arXiv

Bounding Vanishing at the Central Point of Cuspidal Newforms

The Katz-Sarnak Density Conjecture states that zeros of families of $L$-functions are well-modeled by eigenvalues of random matrix ensembles. For suitably restricted test functions, this correspondence yields upper bounds for the families' order of vanishing at the central point. We generalize previous results on the $n$\textsuperscript{th} centered moment of the distribution of zeros to allow arbitrary test functions. On the computational side, we use our improved formulas to obtain significantly better bounds on the order of vanishing for cuspidal newforms, setting world records for the quality of the bounds. We also discover better test functions that further optimize our bounds. We see improvement as early as the $5$\textsuperscript{th} order, and our bounds improve rapidly as the rank grows (more than one order of magnitude better for rank 10 and more than four orders of magnitude for rank 50).

preprint2022arXiv

Class Numbers and Pell's Equation $x^2 + 105y^2 = z^2$

Two well-studied Diophantine equations are those of Pythagorean triples and elliptic curves, for the first we have a parametrization through rational points on the unit circle, and for the second we have a structure theorem for the group of rational solutions. Recently, Yekutieli discussed a connection between these two problems, and described the group structure of Pythagorean triples and the number of triples for a given hypotenuse. In arXiv:2112.03663 we generalized these methods and results to Pell's equation. We find a similar group structure and count on the number of solutions for a given $z$ to $x^2 + Dy^2 = z^2$ when $D$ is 1 or 2 modulo 4 and the class group of $\mathbb{Q}[\sqrt{-D}]$ is a free $\mathbb{Z}_2$ module, which always happens if the class number is at most 2. In this paper, we discuss the main results of arXiv:2112.03663 using some concrete examples in the case of $D=105$.

preprint2022arXiv

Modeling Random Walks to Infinity on Primes in $\mathbb{Z}[\sqrt{2}]$

An interesting question, known as the Gaussian moat problem, asks whether it is possible to walk to infinity on Gaussian primes with steps of bounded length. Our work examines a similar situation in the real quadratic integer ring $\mathbb{Z}[\sqrt{2}]$ whose primes cluster near the asymptotes $y = \pm x/\sqrt{2}$ as compared to Gaussian primes, which cluster near the origin. We construct a probabilistic model of primes in $\mathbb{Z}[\sqrt{2}]$ by applying the prime number theorem and a combinatorial theorem for counting the number of lattice points whose absolute values of their norms are at most $r^2$. We then prove that it is impossible to walk to infinity if the walk remains within some bounded distance from the asymptotes. Lastly, we perform a few moat calculations to show that the longest walk is likely to stay close to the asymptotes; hence, we conjecture that there is no walk to infinity on $\mathbb{Z}[\sqrt{2}]$ primes with steps of bounded length.

preprint2022arXiv

Walking to Infinity on the Fibonacci Sequence

An interesting open problem in number theory asks whether it is possible to walk to infinity on primes, where each term in the sequence has one more digit than the previous. In this paper, we study its variation where we walk on the Fibonacci sequence. We prove that all walks starting with a Fibonacci number and the following terms are Fibonacci numbers obtained by appending exactly one digit at a time to the right have a length of at most two. In the more general case where we append at most a bounded number of digits each time, we give a formula for the length of the longest walk.

preprint2021arXiv

A Refined Conjecture for the Variance of Gaussian Primes Across Sectors

We derive a refined conjecture for the variance of Gaussian primes across sectors, with a power saving error term, by applying the L-functions Ratios Conjecture. We observe a bifurcation point in the main term, consistent with the Random Matrix Theory (RMT) heuristic previously proposed by Rudnick and Waxman. Our model also identifies a second bifurcation point, undetected by the RMT model, that emerges upon taking into account lower order terms. For sufficiently small sectors, we moreover prove an unconditional result that is consistent with our conjecture down to lower order terms.

preprint2021arXiv

Distribution of Eigenvalues of Matrix Ensembles arising from Wigner and Palindromic Toeplitz Blocks

Random Matrix Theory (RMT) has successfully modeled diverse systems, from energy levels of heavy nuclei to zeros of $L$-functions; this correspondence has allowed RMT to successfully predict many number theoretic behaviors. However there are some operations which to date have no RMT analogue. Our motivation is to find an RMT analogue of Rankin-Selberg convolution, which constructs a new $L$-functions from an input pair. We report one such attempt; while it does not appear to model convolution, it does create new ensembles with properties hybridizing those of its constituents. For definiteness we concentrate on the ensemble of palindromic real symmetric Toeplitz (PST) matrices and the ensemble of real symmetric matrices, whose limiting spectral measures are the Gaussian and semi-circular distributions, respectively; these were chosen as they are the two extreme cases in terms of moment calculations. For a PST matrix $A$ and a real symmetric matrix $B$, we construct an ensemble of random real symmetric block matrices whose first row is $\lbrace A, B \rbrace$ and whose second row is $\lbrace B, A \rbrace$. By Markov's Method of Moments and the use of free probability, we show this ensemble converges weakly and almost surely to a new, universal distribution with a hybrid of Gaussian and semi-circular behaviors. We extend this construction by considering an iterated concatenation of matrices from an arbitrary pair of random real symmetric sub-ensembles with different limiting spectral measures. We prove that finite iterations converge to new, universal distributions with hybrid behavior, and that infinite iterations converge to the limiting spectral measure of the dominant component matrix.

preprint2021arXiv

Limiting Spectral Distributions of Families of Block Matrix Ensembles

We introduce a new matrix operation on a pair of matrices, $\text{swirl}(A,X),$ and discuss its implications on the limiting spectral distribution. In a special case, the resultant ensemble converges almost surely to the Rayleigh distribution. In proving this, we provide a novel combinatorial proof that the random matrix ensemble of circulant Hankel matrices converges almost surely to the Rayleigh distribution, using the method of moments.

preprint2021arXiv

Tinkering with Lattices: A New Take on the Erdős Distance Problem

The Erdős distance problem concerns the least number of distinct distances that can be determined by $N$ points in the plane. The integer lattice with $N$ points is known as \textit{near-optimal}, as it spans $Θ(N/\sqrt{\log(N)})$ distinct distances, the lower bound for a set of $N$ points (Erdős, 1946). The only previous non-asymptotic work related to the Erdős distance problem that has been done was for $N \leq 13$. We take a new non-asymptotic approach to this problem in a model case, studying the distance distribution, or in other words, the plot of frequencies of each distance of the $N\times N$ integer lattice. In order to fully characterize this distribution, we adapt previous number-theoretic results from Fermat and Erdős in order to relate the frequency of a given distance on the lattice to the sum-of-squares formula. We study the distance distributions of all the lattice's possible subsets; although this is a restricted case, the structure of the integer lattice allows for the existence of subsets which can be chosen so that their distance distributions have certain properties, such as emulating the distribution of randomly distributed sets of points for certain small subsets, or emulating that of the larger lattice itself. We define an error which compares the distance distribution of a subset with that of the full lattice. The structure of the integer lattice allows us to take subsets with certain geometric properties in order to maximize error; we show these geometric constructions explicitly. Further, we calculate explicit upper bounds for the error when the number of points in the subset is $4$, $5$, $9$ or $\left \lceil N^2/2\right\rceil$ and prove a lower bound in cases with a small number of points.

preprint2020arXiv

Bounds on Zeckendorf Games

Zeckendorf proved that every positive integer $n$ can be written uniquely as the sum of non-adjacent Fibonacci numbers. We use this decomposition to construct a two-player game. Given a fixed integer $n$ and an initial decomposition of $n=n F_1$, the two players alternate by using moves related to the recurrence relation $F_{n+1}=F_n+F_{n-1}$, and whoever moves last wins. The game always terminates in the Zeckendorf decomposition; depending on the choice of moves the length of the game and the winner can vary, though for $n\ge 2$ there is a non-constructive proof that Player 2 has a winning strategy. Initially the lower bound of the length of a game was order $n$ (and known to be sharp) while the upper bound was of size $n \log n$. Recent work decreased the upper bound to of size $n$, but with a larger constant than was conjectured. We improve the upper bound and obtain the sharp bound of $\frac{\sqrt{5}+3}{2}\ n - IZ(n) - \frac{1+\sqrt{5}}{2}Z(n)$, which is of order $n$ as $Z(n)$ is the number of terms in the Zeckendorf decomposition of $n$ and $IZ(n)$ is the sum of indices in the Zeckendorf decomposition of $n$ (which are at most of sizes $\log n$ and $\log^2 n$ respectively). We also introduce a greedy algorithm that realizes the upper bound, and show that the longest game on any $n$ is achieved by applying splitting moves whenever possible.

preprint2020arXiv

Central Limit Theorems for Compound Paths on the 2-Dimensional Lattice

Zeckendorf proved that every integer can be written uniquely as a sum of non-consecutive Fibonacci numbers $\{F_n\}$, and later researchers showed that the distribution of the number of summands needed for such decompositions of integers in $[F_n, F_{n+1})$ converges to a Gaussian as $n\to\infty$. Decomposition problems have been studied extensively for a variety of different sequences and notions of a legal decompositions; for the Fibonacci numbers, a legal decomposition is one for which each summand is used at most once and no two consecutive summands may be chosen. Recently, Chen et al. [CCGJMSY] generalized earlier work to $d$-dimensional lattices of positive integers; there, a legal decomposition is a path such that every point chosen had each component strictly less than the component of the previous chosen point in the path. They were able to prove Gaussianity results despite the lack of uniqueness of the decompositions; however, their results should hold in the more general case where some components are identical. The strictly decreasing assumption was needed in that work to obtain simple, closed form combinatorial expressions, which could then be well approximated and led to the limiting behavior. In this work we remove that assumption through inclusion-exclusion arguments. These lead to more involved combinatorial sums; using generating functions and recurrence relations we obtain tractable forms in $2$ dimensions and prove Gaussianity again; a more involved analysis should work in higher dimensions.

preprint2020arXiv

Constructions of Generalized MSTD Sets in Higher Dimensions

Let $A$ be a set of finite integers, define $$A+A \ = \ \{a_1+a_2: a_1,a_2 \in A\}, \ \ \ A-A \ = \ \{a_1-a_2: a_1,a_2 \in A\},$$ and for non-negative integers $s$ and $d$ define $$sA-dA\ =\ \underbrace{A+\cdots+A}_{s} -\underbrace{A-\cdots-A}_{d}.$$ A More Sums than Differences (MSTD) set is an $A$ where $|A+A| > |A-A|$. It was initially thought that the percentage of subsets of $[0,n]$ that are MSTD would go to zero as $n$ approaches infinity as addition is commutative and subtraction is not. However, in a surprising 2006 result, Martin and O'Bryant proved that a positive percentage of sets are MSTD, although this percentage is extremely small, about $10^{-4}$ percent. This result was extended by Iyer, Lazarev, Miller, ans Zhang [ILMZ] who showed that a positive percentage of sets are generalized MSTD sets, sets for $\{s_1,d_1\} \neq \{s_2, d_2\}$ and $s_1+d_1=s_2+d_2$ with $|s_1A-d_1A| > |s_2A-d_2A|$, and that in $d$-dimensions, a positive percentage of sets are MSTD. For many such results, establishing explicit MSTD sets in $1$-dimensions relies on the specific choice of the elements on the left and right fringes of the set to force certain differences to be missed while desired sums are attained. In higher dimensions, the geometry forces a more careful assessment of what elements have the same behavior as $1$-dimensional fringe elements. We study fringes in $d$-dimensions and use these to create new explicit constructions. We prove the existence of generalized MSTD sets in $d$-dimensions and the existence of $k$-generational sets, which are sets where $|cA+cA|>|cA-cA|$ for all $1\leq c \leq k$. We then prove that under certain conditions, there are no sets with $|kA+kA|>|kA-kA|$ for all $k \in \mathbb{N}.$

preprint2020arXiv

Deterministic Zeckendorf Games

Zeckendorf proved that every positive integer can be written uniquely as the sum of non-adjacent Fibonacci numbers. We further explore a two-player Zeckendorf game introduced in Baird-Smith, Epstein, Flint, and Miller: Given a fixed integer $n$ and an initial decomposition of $n = nF_1$, players alternate using moves related to the recurrence relation $F_{n+1} = F_n + F_{n_1}$, and the last player to move wins. We improve the upper bound on the number of moves possible and show that it is of the same order in $n$ as the lower bound; this is an improvement by a logarithm over previous work. The new upper bound is $3n - 3Z(n) - IZ(n) + 1$, and the existing lower bound is sharp at $n - Z(n)$ moves, where $Z(n)$ is the number of terms in the Zeckendorf decomposition of $n$ and $IZ(n)$ is the sum of indices in the same Zeckendorf decomposition of $n$. We also studied four deterministic variants of the game, where there was a fixed order on which available move one takes: Combine Largest, Split Largest, Combine Smallest and Split Smallest. We prove that Combine Largest and Split Largest realize the lower bound. Split Smallest has the largest number of moves over all possible games, and is close to the new upper bound. For Combine Split games, the number of moves grows linearly with $n$.

preprint2020arXiv

Distribution of missing differences in diffsets

Lazarev, Miller and O'Bryant investigated the distribution of $|S+S|$ for $S$ chosen uniformly at random from $\{0, 1, \dots, n-1\}$, and proved the existence of a divot at missing 7 sums (the probability of missing exactly 7 sums is less than missing 6 or missing 8 sums). We study related questions for $|S-S|$, and shows some divots from one end of the probability distribution, $P(|S-S|=k)$, as well as a peak at $k=4$ from the other end, $P(2n-1-|S-S|=k)$. A corollary of our results is an asymptotic bound for the number of complete rulers of length $n$.

preprint2020arXiv

Extensions of Autocorrelation Inequalities with Applications to Additive Combinatorics

In a 2019 paper, Barnard and Steinerberger show that for $f\in L^1(\mathbf{R})$, the following autocorrelation inequality holds: \begin{equation*} \min_{0 \leq t \leq 1} \int_\mathbf{R} f(x) f(x+t)\ \mathrm{d}x \ \leq\ 0.411 ||f||_{L^1}^2, \end{equation*} where the constant $0.411$ cannot be replaced by $0.37$. In addition to being interesting and important in their own right, inequalities such as these have applications in additive combinatorics where some problems, such as those of minimal difference basis, can be encapsulated by a convolution inequality similar to the above integral. Barnard and Steinerberger suggest that future research may focus on the existence of functions extremizing the above inequality (which is itself related to Brascamp-Lieb type inequalities). We show that for $f$ to be extremal under the above, we must have \begin{equation*} \max_{x_1 \in \mathbf{R} }\min_{0 \leq t \leq 1} \left[ f(x_1-t)+f(x_1+t) \right] \ \leq\ \min_{x_2 \in \mathbf{R} } \max_{0 \leq t \leq 1} \left[ f(x_2-t)+f(x_2+t) \right] . \end{equation*} Our central technique for deriving this result is local perturbation of $f$ to increase the value of the autocorrelation, while leaving $||f||_{L^1}$ unchanged. These perturbation methods can be extended to examine a more general notion of autocorrelation. Let $d,n \in \mathbb{Z}^+$, $f \in L^1$, $A$ be a $d \times n$ matrix with real entries and columns $a_i$ for $1 \leq i \leq n$, and $C$ be a constant. For a broad class of matrices $A$, we prove necessary conditions for $f$ to extremize autocorrelation inequalities of the form \begin{equation*} \min_{ \mathbf{t} \in [0,1]^d } \int_{\mathbf{R}} \prod_{i=1}^n\ f(x+ \mathbf{t} \cdot a_i)\ \mathrm{d}x\ \leq\ C ||f||_{L^1}^n. \end{equation*}

preprint2020arXiv

On the sum of $k$-th powers in terms of earlier sums

For $k$ a positive integer let $S_k(n) = 1^k + 2^k + \cdots + n^k$, i.e., $S_k(n)$ is the sum of the first $k$-th powers. Faulhaber conjectured (later proved by Jacobi) that for $k$ odd, $S_k(n)$ could be written as a polynomial of $S_1(n)$; for example $S_3(n) = S_1(n)^2$. We extend this result and prove that for any $k$ there is a polynomial $g_k(x,y)$ such that $S_k(n) = g(S_1(n), S_2(n))$. The proof yields a recursive formula to evaluate $S_k(n)$ as a polynomial of $n$ that has roughly half the number of terms as the classical one.

preprint2012arXiv

First Order Approximations of the Pythagorean Won-Loss Formula for Predicting MLB Teams' Winning Percentages

We mathematically prove that an existing linear predictor of baseball teams' winning percentages (Jones and Tappin 2005) is simply just a first-order approximation to Bill James' Pythagorean Won-Loss formula and can thus be written in terms of the formula's well-known exponent. We estimate the linear model on twenty seasons of Major League Baseball data and are able to verify that the resulting coefficient estimate, with 95% confidence, is virtually identical to the empirically accepted value of 1.82. Our work thus helps explain why this simple and elegant model is such a strong linear predictor.

preprint2011arXiv

An elliptic curve test of the L-Functions Ratios Conjecture

We compare the L-Function Ratios Conjecture&#39;s prediction with number theory for the family of quadratic twists of a fixed elliptic curve with prime conductor, and show agreement in the 1-level density up to an error term of size X^{-(1-sigma)/2} for test functions supported in (-sigma, sigma); this gives us a power-savings for σ<1. This test of the Ratios Conjecture introduces complications not seen in previous cases (due to the level of the elliptic curve). Further, the results here are one of the key ingredients in the companion paper [DHKMS2], where they are used to determine the effective matrix size for modeling zeros near the central point for this family. The resulting model beautifully describes the behavior of these low lying zeros for finite conductors, explaining the data observed by Miller in [Mil3]. A key ingredient in our analysis is a generalization of Jutila&#39;s bound for sums of quadratic characters with the additional restriction that the fundamental discriminant be congruent to a non-zero square modulo a square-free integer M. This bound is needed for two purposes. The first is to analyze the terms in the explicit formula corresponding to characters raised to an odd power. The second is to determine the main term in the 1-level density of quadratic twists of a fixed form on GL_n. Such an analysis was performed by Rubinstein [Rub], who implicitly assumed that Jutila&#39;s bound held with the additional restriction on the fundamental discriminants; in this paper we show that assumption is justified.

preprint2011arXiv

Finding and Counting MSTD sets

We review the basic theory of More Sums Than Differences (MSTD) sets, specifically their existence, simple constructions of infinite families, the proof that a positive percentage of sets under the uniform binomial model are MSTD but not if the probability that each element is chosen tends to zero, and &#39;explicit&#39; constructions of large families of MSTD sets. We conclude with some new constructions and results of generalized MSTD sets, including among other items results on a positive percentage of sets having a given linear combination greater than another linear combination, and a proof that a positive percentage of sets are $k$-generational sum-dominant (meaning $A$, $A+A$, $...$, $kA = A + ...+A$ are each sum-dominant).

preprint2011arXiv

From Fibonacci Numbers to Central Limit Type Theorems

A beautiful theorem of Zeckendorf states that every integer can be written uniquely as a sum of non-consecutive Fibonacci numbers $\{F_n\}_{n=1}^{\infty}$. Lekkerkerker proved that the average number of summands for integers in $[F_n, F_{n+1})$ is $n/(ϕ^2 + 1)$, with $ϕ$ the golden mean. This has been generalized to the following: given nonnegative integers $c_1,c_2,...,c_L$ with $c_1,c_L>0$ and recursive sequence $\{H_n\}_{n=1}^{\infty}$ with $H_1=1$, $H_{n+1} =c_1H_n+c_2H_{n-1}+...+c_nH_1+1$ $(1\le n< L)$ and $H_{n+1}=c_1H_n+c_2H_{n-1}+...+c_LH_{n+1-L}$ $(n\geq L)$, every positive integer can be written uniquely as $\sum a_iH_i$ under natural constraints on the $a_i$&#39;s, the mean and the variance of the numbers of summands for integers in $[H_{n}, H_{n+1})$ are of size $n$, and the distribution of the numbers of summands converges to a Gaussian as $n$ goes to the infinity. Previous approaches used number theory or ergodic theory. We convert the problem to a combinatorial one. In addition to re-deriving these results, our method generalizes to a multitude of other problems (in the sequel paper \cite{BM} we show how this perspective allows us to determine the distribution of gaps between summands in decompositions). For example, it is known that every integer can be written uniquely as a sum of the $\pm F_n$&#39;s, such that every two terms of the same (opposite) sign differ in index by at least 4 (3). The presence of negative summands introduces complications and features not seen in previous problems. We prove that the distribution of the numbers of positive and negative summands converges to a bivariate normal with computable, negative correlation, namely $-(21-2ϕ)/(29+2ϕ) \approx -0.551058$.

preprint2011arXiv

Gaussian Behavior in Generalized Zeckendorf Decompositions

A beautiful theorem of Zeckendorf states that every integer can be written uniquely as a sum of non-consecutive Fibonacci numbers $\{F_n\}_{n=1}^{\infty}$; Lekkerkerker proved that the average number of summands for integers in $[F_n, F_{n+1})$ is $n/(ϕ^2 + 1)$, with $ϕ$ the golden mean. Interestingly, the higher moments seem to have been ignored. We discuss the proof that the distribution of the number of summands converges to a Gaussian as $n \to \infty$, and comment on generalizations to related decompositions. For example, every integer can be written uniquely as a sum of the $\pm F_n$&#39;s, such that every two terms of the same (opposite) sign differ in index by at least 4 (3). The distribution of the numbers of positive and negative summands converges to a bivariate normal with computable, negative correlation, namely $-(21-2ϕ)/(29+2ϕ) \approx -0.551058$.

preprint2011arXiv

Generalized More Sums Than Differences Sets

A More Sums Than Differences (MSTD, or sum-dominant) set is a finite set $A\subset \mathbb{Z}$ such that $|A+A|<|A-A|$. Though it was believed that the percentage of subsets of $\{0,...,n\}$ that are sum-dominant tends to zero, in 2006 Martin and O&#39;Bryant \cite{MO} proved a positive percentage are sum-dominant. We generalize their result to the many different ways of taking sums and differences of a set. We prove that $|ε_1A+...+ε_kA|>|δ_1A+...+δ_kA|$ a positive percent of the time for all nontrivial choices of $ε_j,δ_j\in \{-1,1\}$. Previous approaches proved the existence of infinitely many such sets given the existence of one; however, no method existed to construct such a set. We develop a new, explicit construction for one such set, and then extend to a positive percentage of sets. We extend these results further, finding sets that exhibit different behavior as more sums/differences are taken. For example, notation as above we prove that for any $m$, $|ε_1A + ... + ε_kA| - |δ_1A + ... + δ_kA| = m$ a positive percentage of the time. We find the limiting behavior of $kA=A+...+A$ for an arbitrary set $A$ as $k\to\infty$ and an upper bound of $k$ for such behavior to settle down. Finally, we say $A$ is $k$-generational sum-dominant if $A$, $A+A$, ...,$kA$ are all sum-dominant. Numerical searches were unable to find even a 2-generational set (heuristics indicate the probability is at most $10^{-9}$, and almost surely significantly less). We prove the surprising result that for any $k$ a positive percentage of sets are $k$-generational, and no set can be $k$-generational for all $k$.

preprint2011arXiv

Low-lying Zeros of Cuspidal Maass Forms

The Katz-Sarnak Density Conjecture states that the behavior of zeros of a family of $L$-functions near the central point (as the conductors tend to zero) agree with the behavior of eigenvalues near 1 of a classical compact group (as the matrix size tends to infinity). Using the Petersson formula, Iwaniec, Luo and Sarnak \cite{ILS} proved that the behavior of zeros near the central point of holomorphic cusp forms agree with the behavior of eigenvalues of orthogonal matrices for suitably restricted test functions. We prove a similar result for level 1 cuspidal Maass forms, the other natural family of ${\rm GL}_2$ $L$-functions. We use the explicit formula to relate sums of our test function at scaled zeros to sums of the Fourier transform at the primes weighted by the $L$-function coefficients, and then use the Kuznetsov trace formula to average the Fourier coefficients over the family. There are numerous technical obstructions in handling the terms in the trace formula, which are surmounted through the use of smooth weight functions for the Maass eigenvalues and results on Kloosterman sums and Bessel and hyperbolic functions.

preprint2011arXiv

The Limiting Spectral Measure for Ensembles of Symmetric Block Circulant Matrices

Given an ensemble of NxN random matrices, a natural question to ask is whether or not the empirical spectral measures of typical matrices converge to a limiting spectral measure as N --> oo. While this has been proved for many thin patterned ensembles sitting inside all real symmetric matrices, frequently there is no nice closed form expression for the limiting measure. Further, current theorems provide few pictures of transitions between ensembles. We consider the ensemble of symmetric m-block circulant matrices with entries i.i.d.r.v. These matrices have toroidal diagonals periodic of period m. We view m as a &#34;dial&#34; we can &#34;turn&#34; from the thin ensemble of symmetric circulant matrices, whose limiting eigenvalue density is a Gaussian, to all real symmetric matrices, whose limiting eigenvalue density is a semi-circle. The limiting eigenvalue densities f_m show a visually stunning convergence to the semi-circle as m tends to infinity, which we prove. In contrast to most studies of patterned matrix ensembles, our paper gives explicit closed form expressions for the densities. We prove that f_m is the product of a Gaussian and a degree 2m-2 polynomial; the formula equals that of the m x m Gaussian Unitary Ensemble (GUE). The proof is by the moments. The new feature, which allows us to obtain closed form expressions, is converting the central combinatorial problem in the moment calculation into an equivalent counting problem in algebraic topology. We end with a generalization of the m-block circulant pattern, dropping the assumption that the m random variables be distinct. We prove that the limiting spectral distribution exists and is determined by the pattern of the independent elements within an m-period, depending on not only the frequency at which each element appears, but also the way the elements are arranged.

preprint2011arXiv

Virus Dynamics on Starlike Graphs

The field of epidemiology has presented fascinating and relevant questions for mathematicians, primarily concerning the spread of viruses in a community. The importance of this research has greatly increased over time as its applications have expanded to also include studies of electronic and social networks and the spread of information and ideas. We study virus propagation on a non-linear hub and spoke graph (which models well many airline networks). We determine the long-term behavior as a function of the cure and infection rates, as well as the number of spokes n. For each n we prove the existence of a critical threshold relating the two rates. Below this threshold, the virus always dies out; above this threshold, all non-trivial initial conditions iterate to a unique non-trivial steady state. We end with some generalizations to other networks.

preprint2010arXiv

A unitary test of the Ratios Conjecture

The Ratios Conjecture of Conrey, Farmer and Zirnbauer predicts the answers to numerous questions in number theory, ranging from n-level densities and correlations to mollifiers to moments and vanishing at the central point. The conjecture gives a recipe to generate these answers, which are believed to be correct up to square-root cancelation. These predictions have been verified, for suitably restricted test functions, for the 1-level density of orthogonal and symplectic families of L-functions. In this paper we verify the conjecture&#39;s predictions for the unitary family of all Dirichlet $L$-functions with prime conductor; we show square-root agreement between prediction and number theory if the support of the Fourier transform of the test function is in (-1,1), and for support up to (-2,2) we show agreement up to a power savings in the family&#39;s cardinality.

preprint2010arXiv

Distribution of Eigenvalues of Highly Palindromic Toeplitz Matrices

Consider the ensemble of real symmetric Toeplitz matrices whose entries are i.i.d random variables chosen from a fixed probability distribution p of mean 0, variance 1 and finite higher moments. Previous work [BDJ,HM] showed that the limiting spectral measures (the density of normalized eigenvalues) converge weakly and almost surely to a universal distribution almost that of the Gaussian, independent of p. The deficit from the Gaussian distribution is due to obstructions to solutions of Diophantine equations and can be removed (see [MMS]) by making the first row palindromic. In this paper, we study the case where there is more than one palindrome in the first row of a real symmetric Toeplitz matrix. Using the method of moments and an analysis of the resulting Diophantine equations, we show that the moments of this ensemble converge to an universal distribution with a fatter tail than any previously seen limiting spectral measure.

preprint2010arXiv

Effective equidistribution and the Sato-Tate law for families of elliptic curves

Extending recent work of others, we provide effective bounds on the family of all elliptic curves and one-parameter families of elliptic curves modulo p (for p prime tending to infinity) obeying the Sato-Tate Law. We present two methods of proof. Both use the framework of Murty-Sinha; the first involves only knowledge of the moments of the Fourier coefficients of the L-functions and combinatorics, and saves a logarithm, while the second requires a Sato-Tate law. Our purpose is to illustrate how the caliber of the result depends on the error terms of the inputs and what combinatorics must be done.

preprint2010arXiv

Low-lying Zeros of Number Field $L$-functions

One of the most important statistics in studying the zeros of L-functions is the 1-level density, which measures the concentration of zeros near the central point. Fouvry and Iwaniec [FI] proved that the 1-level density for L-functions attached to imaginary quadratic fields agrees with results predicted by random matrix theory. In this paper, we show a similar agreement with random matrix theory occurring in more general sequences of number fields. We first show that the main term agrees with random matrix theory, and similar to all other families studied to date, is independent of the arithmetic of the fields. We then derive the first lower order term of the 1-level density, and see the arithmetic enter.

preprint2010arXiv

Modeling Convolutions of $L$-Functions

A number of mathematical methods have been shown to model the zeroes of $L$-functions with remarkable success, including the Ratios Conjecture and Random Matrix Theory. In order to understand the structure of convolutions of families of $L$-functions, we investigate how well these methods model the zeros of such functions. Our primary focus is the convolution of the $L$-function associated to Ramanujan&#39;s tau function with the family of quadratic Dirichlet $L$-functions, for which J.B. Conrey and N.C. Snaith computed the Ratios Conjecture&#39;s prediction. Our main result is performing the number theory calculations and verifying these predictions for the one-level density for suitably restricted test functions up to square-root error term. Unlike Random Matrix Theory, which only predicts the main term, the Ratios Conjecture detects the arithmetic of the family and makes detailed predictions about their dependence in the lower order terms. Interestingly, while Random Matrix Theory is frequently used to model behavior of L-functions (or at least the main terms), there has been little if any work on the analogue of convolving families of L-functions by convolving random matrix ensembles. We explore one possibility by considering Kronecker products; unfortunately, it appears that this is not the correct random matrix analogue to convolving families.

preprint2010arXiv

On the number of summands in Zeckendorf decompositions

Zeckendorf proved that every positive integer has a unique representation as a sum of non-consecutive Fibonacci numbers. Once this has been shown, it&#39;s natural to ask how many summands are needed. Using a continued fraction approach, Lekkerkerker proved that the average number of such summands needed for integers in $[F_n, F_{n+1})$ is $n / (φ^2 + 1) + O(1)$, where $φ= \frac{1+\sqrt{5}}2$ is the golden mean. Surprisingly, no one appears to have investigated the distribution of the number of summands; our main result is that this converges to a Gaussian as $n\to\infty$. Moreover, such a result holds not just for the Fibonacci numbers but many other problems, such as linear recurrence relation with non-negative integer coefficients (which is a generalization of base $B$ expansions of numbers) and far-difference representations. In general the proofs involve adopting a combinatorial viewpoint and analyzing the resulting generating functions through partial fraction expansions and differentiating identities. The resulting arguments become quite technical; the purpose of this paper is to concentrate on the special and most interesting case of the Fibonacci numbers, where the obstructions vanish and the proofs follow from some combinatorics and Stirling&#39;s formula; see [MW] for proofs in the general case.

preprint2009arXiv

An Orthogonal Test of the $L$-functions Ratios Conjecture, II

Recently Conrey, Farmer, and Zirnbauer developed the L-functions Ratios conjecture, which gives a recipe that predicts a wealth of statistics, from moments to spacings between adjacent zeros and values of L-functions. The problem with this method is that several of its steps involve ignoring error terms of size comparable to the main term; amazingly, the errors seem to cancel and the resulting prediction is expected to be accurate up to square-root cancellation. We prove the accuracy of the Ratios Conjecture&#39;s prediction for the 1-level density of families of cuspidal newforms of constant sign (up to square-root agreement for support in (-1,1), and up to a power savings in (-2,2)), and discuss the arithmetic significance of the lower order terms. This is the most involved test of the Ratios Conjecture&#39;s predictions to date, as it is known that the error terms dropped in some of the steps do not cancel, but rather contribute a main term! Specifically, these are the non-diagonal terms in the Petersson formula, which lead to a Bessel-Kloosterman sum which contributes only when the support of the Fourier transform of the test function exceeds (-1, 1).

preprint2009arXiv

Nuclei, Primes and the Random Matrix Connection

In this article, we discuss the remarkable connection between two very different fields, number theory and nuclear physics. We describe the essential aspects of these fields, the quantities studied, and how insights in one have been fruitfully applied in the other. The exciting branch of modern mathematics, random matrix theory, provides the connection between the two fields. We assume no detailed knowledge of number theory, nuclear physics, or random matrix theory; all that is required is some familiarity with linear algebra and probability theory, as well as some results from complex analysis. Our goal is to provide the inquisitive reader with a sound overview of the subjects, placing them in their historical context in a way that is not traditionally given in the popular and technical surveys.

preprint2008arXiv

An Orthogonal Test of the L-Functions Ratios Conjecture

We test the predictions of the L-functions Ratios Conjecture for the family of cuspidal newforms of weight k and level N, with either k fixed and N --> oo through the primes or N=1 and k --> oo. We study the main and lower order terms in the 1-level density. We provide evidence for the Ratios Conjecture by computing and confirming its predictions up to a power savings in the family&#39;s cardinality, at least for test functions whose Fourier transforms are supported in (-2, 2). We do this both for the weighted and unweighted 1-level density (where in the weighted case we use the Petersson weights), thus showing that either formulation may be used. These two 1-level densities differ by a term of size 1 / log(k^2 N). Finally, we show that there is another way of extending the sums arising in the Ratios Conjecture, leading to a different answer (although the answer is such a lower order term that it is hopeless to observe which is correct).

preprint2008arXiv

Chains of distributions, hierarchical Bayesian models and Benford&#39;s Law

Kossovsky recently conjectured that the distribution of leading digits of a chain of probability distributions converges to Benford&#39;s law as the length of the chain grows. We prove his conjecture in many cases, and provide an interpretation in terms of products of independent random variables and a central limit theorem. An interesting consequence is that in hierarchical Bayesian models priors tend to satisfy Benford&#39;s Law as the number of levels of the hierarchy increases, which allows us to develop some simple tests (based on Benford&#39;s law) to test proposed models. We give explicit formulas for the error terms as sums of Mellin transforms, which converges extremely rapidly as the number of terms in the chain grows. We may interpret our results as showing that certain Markov chain Monte Carlo processes are rapidly mixing to Benford&#39;s law.

preprint2008arXiv

Explicit constructions of infinite families of MSTD sets

We explicitly construct infinite families of MSTD (more sums than differences) sets. There are enough of these sets to prove that there exists a constant C such that at least C / r^4 of the 2^r subsets of {1,...,r} are MSTD sets; thus our family is significantly denser than previous constructions (whose densities are at most f(r)/2^{r/2} for some polynomial f(r)). We conclude by generalizing our method to compare linear forms epsilon_1 A + ... + epsilon_n A with epsilon_i in {-1,1}.

preprint2008arXiv

Order Statistics and Benford&#39;s Law

Fix a base B and let zeta have the standard exponential distribution; the distribution of digits of zeta base B is known to be very close to Benford&#39;s Law. If there exists a C such that the distribution of digits of C times the elements of some set is the same as that of zeta, we say that set exhibits shifted exponential behavior base B (with a shift of log_B C \bmod 1). Let X_1, >..., X_N be independent identically distributed random variables. If the X_i&#39;s are drawn from the uniform distribution on [0,L], then as N\to\infty the distribution of the digits of the differences between adjacent order statistics converges to shifted exponential behavior (with a shift of \log_B L/N \bmod 1). By differentiating the cumulative distribution function of the logarithms modulo 1, applying Poisson Summation and then integrating the resulting expression, we derive rapidly converging explicit formulas measuring the deviations from Benford&#39;s Law. Fix a delta in (0,1) and choose N independent random variables from any compactly supported distribution with uniformly bounded first and second derivatives and a second order Taylor series expansion at each point. The distribution of digits of any N^δconsecutive differences \emph{and} all N-1 normalized differences of the order statistics exhibit shifted exponential behavior. We derive conditions on the probability density which determine whether or not the distribution of the digits of all the un-normalized differences converges to Benford&#39;s Law, shifted exponential behavior, or oscillates between the two, and show that the Pareto distribution leads to oscillating behavior.

preprint2008arXiv

The Distribution of the Largest Non-trivial Eigenvalues in Families of Random Regular Graphs

Recently Friedman proved Alon&#39;s conjecture for many families of d-regular graphs, namely that given any epsilon > 0 `most&#39; graphs have their largest non-trivial eigenvalue at most 2 sqrt{d-1}+epsilon in absolute value; if the absolute value of the largest non-trivial eigenvalue is at most 2 sqrt{d-1} then the graph is said to be Ramanujan. These graphs have important applications in communication network theory, allowing the construction of superconcentrators and nonblocking networks, coding theory and cryptography. As many of these applications depend on the size of the largest non-trivial positive and negative eigenvalues, it is natural to investigate their distributions. We show these are well-modeled by the beta=1 Tracy-Widom distribution for several families. If the observed growth rates of the mean and standard deviation as a function of the number of vertices holds in the limit, then in the limit approximately 52% of d-regular graphs from bipartite families should be Ramanujan, and about 27% from non-bipartite families (assuming the largest positive and negative eigenvalues are independent).

preprint2008arXiv

When almost all sets are difference dominated

We investigate the relationship between the sizes of the sum and difference sets attached to a subset of {0,1,...,N}, chosen randomly according to a binomial model with parameter p(N), with N^{-1} = o(p(N)). We show that the random subset is almost surely difference dominated, as N --> oo, for any choice of p(N) tending to zero, thus confirming a conjecture of Martin and O&#39;Bryant. The proofs use recent strong concentration results. Furthermore, we exhibit a threshold phenomenon regarding the ratio of the size of the difference- to the sumset. If p(N) = o(N^{-1/2}) then almost all sums and differences in the random subset are almost surely distinct, and in particular the difference set is almost surely about twice as large as the sumset. If N^{-1/2} = o(p(N)) then both the sum and difference sets almost surely have size (2N+1) - O(p(N)^{-2}), and so the ratio in question is almost surely very close to one. If p(N) = c N^{-1/2} then as c increases from zero to infinity (i.e., as the threshold is crossed), the same ratio almost surely decreases continuously from two to one according to an explicitly given function of c. We also extend our results to the comparison of the generalized difference sets attached to an arbitrary pair of binary linear forms. For certain pairs of forms f and g, we show that there in fact exists a sharp threshold at c_{f,g} N^{-1/2}, for some computable constant c_{f,g}, such that one form almost surely dominates below the threshold, and the other almost surely above it. The heart of our approach involves using different tools to obtain strong concentration of the sizes of the sum and difference sets about their mean values, for various ranges of the parameter p.

preprint2007arXiv

A Symplectic Test of the L-Functions Ratios Conjecture

Recently Conrey, Farmer and Zirnbauer conjectured formulas for the averages over a family of ratios of products of shifted L-functions. Their L-functions Ratios Conjecture predicts both the main and lower order terms for many problems, ranging from n-level correlations and densities to mollifiers and moments to vanishing at the central point. There are now many results showing agreement between the main terms of number theory and random matrix theory; however, there are very few families where the lower order terms are known. These terms often depend on subtle arithmetic properties of the family, and provide a way to break the universality of behavior. The L-functions Ratios Conjecture provides a powerful and tractable way to predict these terms. We test a specific case here, that of the 1-level density for the symplectic family of quadratic Dirichlet characters arising from even fundamental discriminants d \le X. For test functions supported in (-1/3, 1/3) we calculate all the lower order terms up to size O(X^{-1/2+epsilon}) and observe perfect agreement with the conjecture (for test functions supported in (-1, 1) we show agreement up to errors of size O(X^{-epsilon}) for any epsilon). Thus for this family and suitably restricted test functions, we completely verify the Ratios Conjecture&#39;s prediction for the 1-level density.

preprint2007arXiv

The Modulo 1 Central Limit Theorem and Benford&#39;s Law for Products

We derive a necessary and sufficient condition for the sum of M independent continuous random variables modulo 1 to converge to the uniform distribution in L^1([0,1]), and discuss generalizations to discrete random variables. A consequence is that if X_1, ..., X_M are independent continuous random variables with densities f_1, ..., f_M, for any base B as M \to \infty for many choices of the densities the distribution of the digits of X_1 * ... * X_M converges to Benford&#39;s law base B. The rate of convergence can be quantified in terms of the Fourier coefficients of the densities, and provides an explanation for the prevalence of Benford behavior in many diverse systems.

preprint2006arXiv

A Derivation of the Pythagorean Won-Loss Formula in Baseball

It has been noted that in many professional sports leagues a good predictor of a team&#39;s won-loss percentage is Bill James&#39; Pythagorean Formula RSobs^c / (RSobs^c + RAobs^c), where RSobs (resp. RAobs) is the observed average number of runs scored (allowed) per game and c is a constant for the league; for baseball the best agreement is when c is about 1.82. We provide a theoretical justification for this formula and value of c by modelling the number of runs scored and allowed in baseball games as independent random variables drawn from Weibull distributions with the same b and c but different a; the probability density f(x;a,b,c) is 0 for x < b and is (c/a) ((x-b)/a)^{c-1} exp(-((x-b)/a)^c) otherwise. This model leads to a predicted won-loss percentage of (RS-b)^c / ((RS-b)^c + (RA-b)^c); here RS (resp. RA) is the mean of the random variable corresponding to runs scored (allowed), and RS - b (resp. RA - b) is an estimator of RSobs (resp. RAobs). An analysis of the 14 American League teams from the 2004 baseball season shows that (1) given that the runs scored and allowed in a game cannot be equal, the runs scored and allowed are statistically independent; (2) the best fit Weibull parameters attained from a least squares or a maximum likelihood analysis give good fits; least squares gives a mean value of c of 1.79 with a standard deviation of .09, and maximum likelihood gives a mean value of c of 1.74 with a standard deviation of .06, which agree beautifully with the observed best value of 1.82 attained by fitting RSobs^c / (RSobs^c + RAobs^c) to the observed winning percentages.

preprint2006arXiv

Closed-Form Bayesian Inferences for the Logit Model via Polynomial Expansions

Articles in Marketing and choice literatures have demonstrated the need for incorporating person-level heterogeneity into behavioral models (e.g., logit models for multiple binary outcomes as studied here). However, the logit likelihood extended with a population distribution of heterogeneity doesn&#39;t yield closed-form inferences, and therefore numerical integration techniques are relied upon (e.g., MCMC methods). We present here an alternative, closed-form Bayesian inferences for the logit model, which we obtain by approximating the logit likelihood via a polynomial expansion, and then positing a distribution of heterogeneity from a flexible family that is now conjugate and integrable. For problems where the response coefficients are independent, choosing the Gamma distribution leads to rapidly convergent closed-form expansions; if there are correlations among the coefficients one can still obtain rapidly convergent closed-form expansions by positing a distribution of heterogeneity from a Multivariate Gamma distribution. The solution then comes from the moment generating function of the Multivariate Gamma distribution or in general from the multivariate heterogeneity distribution assumed. Closed-form Bayesian inferences, derivatives (useful for elasticity calculations), population distribution parameter estimates (useful for summarization) and starting values (useful for complicated algorithms) are hence directly available. Two simulation studies demonstrate the efficacy of our approach.

preprint2006arXiv

Distribution of Eigenvalues of Real Symmetric Palindromic Toeplitz Matrices and Circulant Matrices

Consider the ensemble of real symmetric Toeplitz matrices, each independent entry an i.i.d. random variable chosen from a fixed probability distribution p of mean 0, variance 1, and finite higher moments. Previous investigations showed that the limiting spectral measure (the density of normalized eigenvalues) converges weakly and almost surely, independent of p, to a distribution which is almost the standard Gaussian. The deviations from Gaussian behavior can be interpreted as arising from obstructions to solutions of Diophantine equations. We show that these obstructions vanish if instead one considers real symmetric palindromic Toeplitz matrices, matrices where the first row is a palindrome. A similar result was previously proved for a related circulant ensemble through an analysis of the explicit formulas for eigenvalues. By Cauchy&#39;s interlacing property and the rank inequality, this ensemble has the same limiting spectral distribution as the palindromic Toeplitz matrices; a consequence of combining the two approaches is a version of the almost sure Central Limit Theorem. Thus our analysis of these Diophantine equations provides an alternate technique for proving limiting spectral measures for certain ensembles of circulant matrices.

preprint2006arXiv

Low lying zeros of L-functions with orthogonal symmetry

We investigate the moments of a smooth counting function of the zeros near the central point of L-functions of weight k cuspidal newforms of prime level N. We split by the sign of the functional equations and show that for test functions whose Fourier transform is supported in (-1/n, 1/n), as N --> oo the first n centered moments are Gaussian. By extending the support to (-1/n-1, 1/n-1), we see non-Gaussian behavior; in particular the odd centered moments are non-zero for such test functions. If we do not split by sign, we obtain Gaussian behavior for support in (-2/n, 2/n) if 2k >= n. The nth centered moments agree with Random Matrix Theory in this extended range, providing additional support for the Katz-Sarnak conjectures. The proof requires calculating multidimensional integrals of the non-diagonal terms in the Bessel-Kloosterman expansion of the Petersson formula. We convert these multidimensional integrals to one-dimensional integrals already considered in the work of Iwaniec-Luo-Sarnak, and derive a new and more tractable expression for the nth centered moments for such test functions. This new formula facilitates comparisons between number theory and random matrix theory for test functions supported in (-1/n-1, 1/n-1) by simplifying the combinatorial arguments. As an application we obtain bounds for the percentage of such cusp forms with a given order of vanishing at the central point.

preprint2006arXiv

The Low Lying Zeros of a GL(4) and a GL(6) family of L-functions

We investigate the large weight (k --> oo) limiting statistics for the low lying zeros of a GL(4) and a GL(6) family of L-functions, {L(s,phi x f): f in H_k(1)} and {L(s,phi times sym^2 f): f in H_k(1)}; here phi is a fixed even Hecke-Maass cusp form and H_k(1) is a Hecke eigenbasis for the space H_k(1) of holomorphic cusp forms of weight k for the full modular group. Katz and Sarnak conjecture that the behavior of zeros near the central point should be well modeled by the behavior of eigenvalues near 1 of a classical compact group. By studying the 1- and 2-level densities, we find evidence of underlying symplectic and SO(even) symmetry, respectively. This should be contrasted with previous results of Iwaniec-Luo-Sarnak for the families {L(s,f): f in H_k(1)} and {L(s,sym^2f): f in H_k(1)}, where they find evidence of orthogonal and symplectic symmetry, respectively. The present examples suggest a relation between the symmetry type of a family and that of its twistings, which will be further studied in a subsequent paper. Both the GL(4) and the GL(6) families above have all even functional equations, and neither is naturally split from an orthogonal family. A folklore conjecture states that such families must be symplectic, which is true for the first family but false for the second. Thus the theory of low lying zeros is more than just a theory of signs of functional equations. An analysis of these families suggest that it is the second moment of the Satake parameters that determines the symmetry group.

preprint2005arXiv

Incomplete Quadratic Exponential Sums in Several Variables

We consider incomplete exponential sums in several variables of the form S(f,n,m) = \frac{1}{2^n} \sum_{x_1 \in \{-1,1\}} ... \sum_{x_n \in \{-1,1\}} x_1 ... x_n e^{2πi f(x)/p}, where m>1 is odd and f is a polynomial of degree d with coefficients in Z/mZ. We investigate the conjecture, originating in a problem in computational complexity, that for each fixed d and m the maximum norm of S(f,n,m) converges exponentially fast to 0 as n grows to infinity. The conjecture is known to hold in the case when m=3 and d=2, but existing methods for studying incomplete exponential sums appear to be insufficient to resolve the question for an arbitrary odd modulus m, even when d=2. In the present paper we develop three separate techniques for studying the problem in the case of quadratic f, each of which establishes a different special case of the conjecture. We show that a bound of the required sort holds for almost all quadratic polynomials, a stronger form of the conjecture holds for all quadratic polynomials with no more than 10 variables, and for arbitrarily many variables the conjecture is true for a class of quadratic polynomials having a special form.

preprint2005arXiv

Investigations of Zeros Near the Central Point of Elliptic Curve L-Functions

We explore the effect of zeros at the central point on nearby zeros of elliptic curve L-functions, especially for one-parameter families of rank r over Q. By the Birch and Swinnerton Dyer Conjecture and Silverman&#39;s Specialization Theorem, for t sufficiently large the L-function of each curve E_t in the family has r zeros (called the family zeros) at the central point. We observe experimentally a repulsion of the zeros near the central point, and the repulsion increases with r. There is greater repulsion in the subset of curves of rank r+2 than in the subset of curves of rank r in a rank r family. For curves with comparable conductors, the behavior of rank 2 curves in a rank 0 one-parameter family over Q is statistically different from that of rank 2 curves from a rank 2 family. Unlike excess rank calculations, the repulsion decreases markedly as the conductors increase, and we conjecture that the r family zeros do not repel in the limit. Finally, the differences between adjacent normalized zeros near the central point are statistically independent of the repulsion, family rank and rank of the curves in the subset. Specifically, the differences between adjacent normalized zeros are statistically equal for all curves investigated with rank 0, 2 or 4 and comparable conductors from one-parameter families of rank 0 or 2 over Q.

preprint2005arXiv

Variation in the number of points on elliptic curves and applications to excess rank

Michel proved that for a one-parameter family of elliptic curves over Q(T) with non-constant j(T) that the second moment of the number of solutions modulo p is p^2 + O(p^{3/2}). We show this bound is sharp by studying y^2 = x^3 + Tx^2 + 1. Lower order terms for such moments in a family are related to lower order terms in the n-level densities of Katz and Sarnak, which describe the behavior of the zeros near the central point of the associated L-functions. We conclude by investigating similar families and show how the lower order terms in the second moment may affect the expected bounds for the average rank of families in numerical investigations.

preprint2004arXiv

Constructing Elliptic Curves over $\mathbb{Q}(T)$ with Moderate Rank

We give several new constructions for moderate rank elliptic curves over $\mathbb{Q}(T)$. In particular we construct infinitely many rational elliptic surfaces (not in Weierstrass form) of rank 6 over $\mathbb{Q}$ using polynomials of degree two in $T$. While our method generates linearly independent points, we are able to show the rank is exactly 6 \emph{without} having to verify the points are independent. The method generalizes; however, the higher rank surfaces are not rational, and we need to check that the constructed points are linearly independent.

preprint2003arXiv

1- and 2-Level Densities for Rational Families of Elliptic Curves: Evidence for the Underlying Group Symmetries

Following Katz-Sarnak, Iwaniec-Luo-Sarnak, and Rubinstein, we use the 1- and 2-level densities to study the distribution of low lying zeros for one-parameter rational families of elliptic curves over Q(t). Modulo standard conjectures, for small support the densities agree with Katz and Sarnak&#39;s predictions. Further, the densities confirm that the curves&#39; L-functions behave in a manner consistent with having r zeros at the critical point, as predicted by the Birch and Swinnerton-Dyer conjecture. By studying the 2-level densities of some constant sign families, we find the first examples of families of elliptic curves where we can distinguish SO(even) from SO(odd) symmetry.

preprint2003arXiv

Eigenvalue Spacing Distribution for the Ensemble of Real Symmetric Toeplitz Matrices

Consider the ensemble of Real Symmetric Toeplitz Matrices, each entry iidrv from a fixed probability distribution p of mean 0, variance 1, and finite higher moments. The limiting spectral measure (the density of normalized eigenvalues) converges weakly to a new universal distribution with unbounded support, independent of p. This distribution&#39;s moments are almost those of the Gaussian&#39;s; the deficit may be interpreted in terms of Diophantine obstructions. With a little more work, we obtain almost sure convergence. An investigation of spacings between adjacent normalized eigenvalues looks Poissonian, and not GOE.