Is Random Close Packing of Spheres Well Defined?

Despite its long history, there are many fundamental issues concerning random packings of spheres that remain elusive, including a precise definition of random close packing (RCP). We argue that the current picture of RCP cannot be made mathematically precise and support this conclusion via a molecular dynamics study of hard spheres using the Lubachevsky-Stillinger compression algorithm. We suggest that this impasse can be broken by introducing the new concept of a maximally random jammed state, which can be made precise.Comment: 6 pages total, 2 figure

Nuclear Shadowing in DIS: Numerical Solution of the Evolution Equation for the Green Function

Within a light-cone QCD formalism based on the Green function technique incorporating color transparency and coherence length effects we study nuclear shadowing in deep-inelastic scattering at moderately small Bjorken x_{Bj}. Calculations performed so far were based only on approximations leading to an analytical harmonic oscillatory form of the Green function. We present for the first time an exact numerical solution of the evolution equation for the Green function using realistic form of the dipole cross section and nuclear density function. We compare numerical results for nuclear shadowing with previous predictions and discuss differences.Comment: 21 pages including 3 figures; a small revision of the tex

Anisotropy in the helicity modulus of a 3D XY-model: application to YBCO

We present a Monte Carlo study of the helicity moduli of an anisotropic classical three-dimensional (3D) XY-model of YBCO in superconducting state. It is found that both the ab-plane and the c-axis helicity moduli, which are proportional to the inverse square of the corresponding magnetic field penetration depth, vary linearly with temperature at low temperatures. The result for the c-axis helicity modulus is in disagreement with the experiments on high quality samples of YBCO. Thus we conclude that purely classical phase fluctuations of the superconducting order parameter cannot account for the observed c-axis electrodynamics of YBCO.Comment: 7 pages, 1 figur

Kinetic Inductance of Josephson Junction Arrays: Dynamic and Equilibrium Calculations

We show analytically that the inverse kinetic inductance $L^{-1}$ of an overdamped junction array at low frequencies is proportional to the admittance of an inhomogeneous equivalent impedance network. The $ij^{th}$ bond in this equivalent network has an inverse inductance $J_{ij}\cos(\theta_i^0-\theta_j^0-A_{ij})$, where $J_{ij}$ is the Josephson coupling energy of the $ij^{th}$ bond, $\theta_i^0$ is the ground-state phase of the grain $i$, and $A_{ij}$ is the usual magnetic phase factor. We use this theorem to calculate $L^{-1}$ for square arrays as large as $180\times 180$. The calculated $L^{-1}$ is in very good agreement with the low-temperature limit of the helicity modulus $\gamma$ calculated by conventional equilibrium Monte Carlo techniques. However, the finite temperature structure of $\gamma$, as a function of magnetic field, is \underline{sharper} than the zero-temperature $L^{-1}$, which shows surprisingly weak structure. In triangular arrays, the equilibrium calculation of $\gamma$ yields a series of peaks at frustrations $f = \frac{1}{2}(1-1/N)$, where $N$ is an integer $\geq 2$, consistent with experiment.Comment: 14 pages + 6 postscript figures, 3.0 REVTe

Universal Short-Time Dynamics in the Kosterlitz-Thouless Phase

We study the short-time dynamics of systems that develop ``quasi long-range order'' after a quench to the Kosterlitz-Thouless phase. With the working hypothesis that the ``universal short-time behavior'', previously found in Ising-like systems, also occurs in the Kosterlitz-Thouless phase, we explore the scaling behavior of thermodynamic variables during the relaxational process following the quench. As a concrete example, we investigate the two-dimensional $6$-state clock model by Monte Carlo simulation. The exponents governing the magnetization, the second moment, and the autocorrelation function are calculated. From them, by means of scaling relations, estimates for the equilibrium exponents $z$ and $\eta$ are derived. In particular, our estimates for the temperature-dependent anomalous dimension $\eta$ that governs the static correlation function are consistent with existing analytical and numerical results and, thus, confirm our working hypothesis.Comment: 16 pages, 9 postscript figures, REVTEX 3.0, submitted to Phys. Rev.

Spin glass behavior of frustrated 2-D Penrose lattice in the classical planar model

Via extensive Monte Carlo studies we show that the frustrated XY Hamiltonian on a 2-D Penrose lattice admits of a spin glass phase at low temperature. Studies of the Edwards-Anderson order parameter, spin glass susceptibility, and local (linear) susceptibility point unequivocally to a paramagnetic to spin glass transition as the temperature is lowered. Specific heat shows a rounded peak at a temperature above the spin glass transition temperature, as is commonly observed in spin glasses. Our results strongly suggest that the critical point exponents are the same as obtained by Bhatt and Young in the ${\pm}J$ Ising model on a square lattice. However, unlike in the latter case, the critical temperature is clearly finite (nonzero). The results imply that a quasiperiodic 2-D array of superconducting grains in a suitably chosen transverse magnetic field should behave as a superconducting glass at low temperature.Comment: RevTex, 4 pages Including 4 figures. To appear in the June 1 1996 issue of Phys. Rev. B (Rapid Communications). Revised/replaced edition contains an erratum at the end of the paper, also to appear in Phys. Rev.