Advances in random matrix theory, zeta functions, and sphere packing

Over four hundred years ago, Sir Walter Raleigh asked his mathematical assistant to find formulas for the number of cannonballs in regularly stacked piles. These investigations aroused the curiosity of the astronomer Johannes Kepler and led to a problem that has gone centuries without a solution: why is the familiar cannonball stack the most efficient arrangement possible? Here we discuss the solution that Hales found in 1998. Almost every part of the 282-page proof relies on long computer verifications. Random matrix theory was developed by physicists to describe the spectra of complex nuclei. In particular, the statistical fluctuations of the eigenvalues (“the energy levels”) follow certain universal laws based on symmetry types. We describe these and then discuss the remarkable appearance of these laws for zeros of the Riemann zeta function (which is the generating function for prime numbers and is the last special function from the last century that is not understood today.) Explaining this phenomenon is a central problem. These topics are distinct, so we present them separately with their own introductory remarks.

Titration of chaos with added noise

Deterministic chaos has been implicated in numerous natural and man-made complex phenomena ranging from quantum to astronomical scales and in disciplines as diverse as meteorology, physiology, ecology, and economics. However, the lack of a definitive test of chaos vs. random noise in experimental time series has led to considerable controversy in many fields. Here we propose a numerical titration procedure as a simple “litmus test” for highly sensitive, specific, and robust detection of chaos in short noisy data without the need for intensive surrogate data testing. We show that the controlled addition of white or colored noise to a signal with a preexisting noise floor results in a titration index that: (i) faithfully tracks the onset of deterministic chaos in all standard bifurcation routes to chaos; and (ii) gives a relative measure of chaos intensity. Such reliable detection and quantification of chaos under severe conditions of relatively low signal-to-noise ratio is of great interest, as it may open potential practical ways of identifying, forecasting, and controlling complex behaviors in a wide variety of physical, biomedical, and socioeconomic systems.

Critical transport properties of random metals in large magnetic fields

The threshold behavior of the transport properties of a random metal in the critical region near a metal–insulator transition is strongly affected by the measuring electromagnetic fields. In spite of the randomness, the electrical conductivity exhibits striking phase-coherent effects due to broken symmetry, which greatly sharpen the transition compared with the predictions of effective medium theories, as previously explained for electrical conductivities. Here broken symmetry explains the sign reversal of the T → 0 magnetoconductance of the metal–insulator transition in Si(B,P), also previously not understood by effective medium theories. Finally, the symmetry-breaking features of quantum percolation theory explain the unexpectedly very small electrical conductivity temperature exponent α = 0.22(2) recently observed in Ni(S,Se)2 alloys at the antiferromagnetic metal–insulator transition below T = 0.8 K.

Ultrasonic tagging of light: Theory

A theory is provided for the detection efficiency of diffuse light whose frequency is modulated by an acoustical wave. We derive expressions for the speckle pattern of the modulated light, as well as an expression for the signal-to-noise ratio for the detector. The aim is to develop a new imaging technology for detection of tumors in humans. The acoustic wave is focused into a small geometrical volume, which provides the spatial resolution for the imaging. The wavelength of the light wave can be selected to provide information regarding the kind of tumor.

Spectral bifurcations in dispersive wave turbulence

Dispersive wave turbulence is studied numerically for a class of one-dimensional nonlinear wave equations. Both deterministic and random (white noise in time) forcings are studied. Four distinct stable spectra are observed—the direct and inverse cascades of weak turbulence (WT) theory, thermal equilibrium, and a fourth spectrum (MMT; Majda, McLaughlin, Tabak). Each spectrum can describe long-time behavior, and each can be only metastable (with quite diverse lifetimes)—depending on details of nonlinearity, forcing, and dissipation. Cases of a long-live MMT transient state dcaying to a state with WT spectra, and vice-versa, are displayed. In the case of freely decaying turbulence, without forcing, both cascades of weak turbulence are observed. These WT states constitute the clearest and most striking numerical observations of WT spectra to date—over four decades of energy, and three decades of spatial, scales. Numerical experiments that study details of the composition, coexistence, and transition between spectra are then discussed, including: (i) for deterministic forcing, sharp distinctions between focusing and defocusing nonlinearities, including the role of long wavelength instabilities, localized coherent structures, and chaotic behavior; (ii) the role of energy growth in time to monitor the selection of MMT or WT spectra; (iii) a second manifestation of the MMT spectrum as it describes a self-similar evolution of the wave, without temporal averaging; (iv) coherent structures and the evolution of the direct and inverse cascades; and (v) nonlocality (in k-space) in the transferral process.

Noise in neurons is message dependent

Neuronal responses are conspicuously variable. We focus on one particular aspect of that variability: the precision of action potential timing. We show that for common models of noisy spike generation, elementary considerations imply that such variability is a function of the input, and can be made arbitrarily large or small by a suitable choice of inputs. Our considerations are expected to extend to virtually any mechanism of spike generation, and we illustrate them with data from the visual pathway. Thus, a simplification usually made in the application of information theory to neural processing is violated: noise is not independent of the message. However, we also show the existence of error-correcting topologies, which can achieve better timing reliability than their components.

The vibrational energy flow transition in organic molecules: Theory meets experiment

Most large dynamical systems are thought to have ergodic dynamics, whereas small systems may not have free interchange of energy between degrees of freedom. This assumption is made in many areas of chemistry and physics, ranging from nuclei to reacting molecules and on to quantum dots. We examine the transition to facile vibrational energy flow in a large set of organic molecules as molecular size is increased. Both analytical and computational results based on local random matrix models describe the transition to unrestricted vibrational energy flow in these molecules. In particular, the models connect the number of states participating in intramolecular energy flow to simple molecular properties such as the molecular size and the distribution of vibrational frequencies. The transition itself is governed by a local anharmonic coupling strength and a local state density. The theoretical results for the transition characteristics compare well with those implied by experimental measurements using IR fluorescence spectroscopy of dilution factors reported by Stewart and McDonald [Stewart, G. M. & McDonald, J. D. (1983) J. Chem. Phys. 78, 3907–3915].

Seven newly discovered intron positions in the triose-phosphate isomerase gene: evidence for the introns-late theory.

The gene encoding the glycolytic enzyme triose-phosphate isomerase (TPI; EC 5.3.1.1) has been central to the long-standing controversy on the origin and evolutionary significance of spliceosomal introns by virtue of its pivotal support for the introns-early view, or exon theory of genes. Putative correlations between intron positions and TPI protein structure have led to the conjecture that the gene was assembled by exon shuffling, and five TPI intron positions are old by the criterion of being conserved between animals and plants. We have sequenced TPI genes from three diverse eukaryotes--the basidiomycete Coprinus cinereus, the nematode Caenorhabditis elegans, and the insect Heliothis virescens--and have found introns at seven novel positions that disrupt previously recognized gene/protein structure correlations. The set of 21 TPI introns now known is consistent with a random model of intron insertion. Twelve of the 21 TPI introns appear to be of recent origin since each is present in but a single examined species. These results, together with their implication that as more TPI genes are sequenced more intron positions will be found, render TPI untenable as a paradigm for the introns-early theory and, instead, support the introns-late view that spliceosomal introns have been inserted into preexisting genes during eukaryotic evolution.

Founder-effect speciation theory: failure of experimental corroboration.

The theory of founder-effect speciation proposes that colonization by very few individuals of an empty habitat favors rapid genetic changes and the evolution of a new species. We report here the results obtained in a 10-year-long and large-scale experiment with Drosophila pseudoobscura designed to test the theory. In our experimental protocol, populations are established with variable numbers of very few individuals and allowed to expand greatly for several generations until conditions of severe competition for resources are reached and the population crashes. A few random survivors are then taken to start a new population expansion and thus initiate a new cycle of founding events, population flushes, and crashes. Our results provide no support for the theories proposing that new species are very likely to appear as by-products of founder events.