Controlling the Short-Range Order and Packing Densities of Many-Particle Systems

Questions surrounding the spatial disposition of particles in various condensed-matter systems continue to pose many theoretical challenges. This paper explores the geometric availability of amorphous many-particle configurations that conform to a given pair correlation function g(r). Such a study is required to observe the basic constraints of non-negativity for g(r) as well as for its structure factor S(k). The hard sphere case receives special attention, to help identify what qualitative features play significant roles in determining upper limits to maximum amorphous packing densities. For that purpose, a five-parameter test family of g's has been considered, which incorporates the known features of core exclusion, contact pairs, and damped oscillatory short-range order beyond contact. Numerical optimization over this five-parameter set produces a maximum-packing value for the fraction of covered volume, and about 5.8 for the mean contact number, both of which are within the range of previous experimental and simulational packing results. However, the corresponding maximum-density g(r) and S(k) display some unexpected characteristics. A byproduct of our investigation is a lower bound on the maximum density for random sphere packings in $d$ dimensions, which is sharper than a well-known lower bound for regular lattice packings for d >= 3.Comment: Appeared in Journal of Physical Chemistry B, vol. 106, 8354 (2002). Note Errata for the journal article concerning typographical errors in Eq. (11) can be found at http://cherrypit.princeton.edu/papers.html However, the current draft on Cond-Mat (posted on August 8, 2002) is correct

Ground states of stealthy hyperuniform potentials. II. Stacked-slider phases

Stealthy potentials, a family of long-range isotropic pair potentials, produce infinitely degenerate disordered ground states at high densities and crystalline ground states at low densities in d-dimensional Euclidean space R^d. In the previous paper in this series, we numerically studied the entropically favored ground states in the canonical ensemble in the zero-temperature limit across the first three Euclidean space dimensions. In this paper, we investigate using both numerical and theoretical techniques metastable stacked-slider phases, which are part of the ground-state manifold of stealthy potentials at densities in which crystal ground states are favored entropically. Our numerical results enable us to devise analytical models of this phase in two, three, and higher dimensions. Utilizing this model, we estimated the size of the feasible region in configuration space of the stacked-slider phase, finding it to be smaller than that of crystal structures in the infinite-system-size limit, which is consistent with our recent previous work. In two dimensions, we also determine exact expressions for the pair correlation function and structure factor of the analytical model of stacked-slider phases and analyze the connectedness of the ground-state manifold of stealthy potentials in this density regime. We demonstrate that stacked-slider phases are distinguishable states of matter; they are nonperiodic, statistically anisotropic structures that possess long-range orientational order but have zero shear modulus. We outline some possible future avenues of research to elucidate our understanding of this unusual phase of matter

Inverse Statistical Mechanics: Probing the Limitations of Isotropic Pair Potentials to Produce Ground-State Structural Extremes

Inverse statistical-mechanical methods have recently been employed to design optimized short-ranged radial (isotropic) pair potentials that robustly produce novel targeted classical ground-state many-particle configurations. The target structures considered in those studies were low-coordinated crystals with a high degree of symmetry. In this paper, we further test the fundamental limitations of radial pair potentials by targeting crystal structures with appreciably less symmetry, including those in which the particles have different local structural environments. These challenging target configurations demanded that we modify previous inverse optimization techniques. Using this modified optimization technique, we have designed short-ranged radial pair potentials that stabilize the two-dimensional kagome crystal, the rectangular kagome crystal, and rectangular lattices, as well as the three-dimensional structure of CaF$_2$ crystal inhabited by a single particle species. We verify our results by cooling liquid configurations to absolute zero temperature via simulated annealing and ensuring that such states have stable phonon spectra. Except for the rectangular kagome structure, all of the target structures can be stabilized with monotonic repulsive potentials. Our work demonstrates that single-component systems with short-ranged radial pair potentials can counterintuitively self-assemble into crystal ground states with low symmetry and different local structural environments. Finally, we present general principles that offer guidance in determining whether certain target structures can be achieved as ground states by radial pair potentials

Ground states of stealthy hyperuniform potentials: I. Entropically favored configurations

Systems of particles interacting with "stealthy" pair potentials have been shown to possess infinitely degenerate disordered hyperuniform classical ground states with novel physical properties. Previous attempts to sample the infinitely degenerate ground states used energy minimization techniques, introducing algorithmic dependence that is artificial in nature. Recently, an ensemble theory of stealthy hyperuniform ground states was formulated to predict the structure and thermodynamics that was shown to be in excellent agreement with corresponding computer simulation results in the canonical ensemble (in the zero-temperature limit). In this paper, we provide details and justifications of the simulation procedure, which involves performing molecular dynamics simulations at sufficiently low temperatures and minimizing the energy of the snapshots for both the high-density disordered regime, where the theory applies, as well as lower densities. We also use numerical simulations to extend our study to the lower-density regime. We report results for the pair correlation functions, structure factors, and Voronoi cell statistics. In the high-density regime, we verify the theoretical ansatz that stealthy disordered ground states behave like "pseudo" disordered equilibrium hard-sphere systems in Fourier space. These results show that as the density decreases from the high-density limit, the disordered ground states in the canonical ensemble are characterized by an increasing degree of short-range order and eventually the system undergoes a phase transition to crystalline ground states. We also provide numerical evidence suggesting that different forms of stealthy pair potentials produce the same ground-state ensemble in the zero-temperature limit. Our techniques may be applied to sample this limit of the canonical ensemble of other potentials with highly degenerate ground states

Negative thermal expansion in single-component systems with isotropic interactions

We have devised an isotropic interaction potential that gives rise to negative thermal expansion (NTE) behavior in equilibrium many-particle systems in both two and three dimensions over a wide temperature and pressure range (including zero pressure). An optimization procedure is used in order to find a potential that yields a strong NTE effect. A key feature of the potential that gives rise to this behavior is the softened interior of its basin of attraction. Although such anomalous behavior is well known in material systems with directional interactions (e.g., zirconium tungstate), to our knowledge this is the first time that NTE behavior has been established to occur in single-component many-particle systems for isotropic interactions. Using constant-pressure Monte Carlo simulations, we show that as the temperature is increased, the system exhibits negative, zero and then positive thermal expansion before melting (for both two- and three-dimensional systems). The behavior is explicitly compared to that of a Lennard-Jones system, which exhibits typical expansion upon heating for all temperatures and pressures.Comment: 21 pages, 13 figure

Hyperuniformity Order Metric of Barlow Packings

The concept of hyperuniformity has been a useful tool in the study of large-scale density fluctuations in systems ranging across the natural and mathematical sciences. One can rank a large class of hyperuniform systems by their ability to suppress long-range density fluctuations through the use of a hyperuniformity order metric $\bar{\Lambda}$. We apply this order metric to the Barlow packings, which are the infinitely degenerate densest packings of identical rigid spheres that are distinguished by their stacking geometries and include the commonly known fcc lattice and hcp crystal. The "stealthy stacking" theorem implies that these packings are all stealthy hyperuniform, a strong type of hyperuniformity which involves the suppression of scattering up to a wavevector $K$. We describe the geometry of three classes of Barlow packings, two disordered classes and small-period packings. In addition, we compute a lower bound on $K$ for all Barlow packings. We compute $\bar{\Lambda}$ for the aforementioned three classes of Barlow packings and find that to a very good approximation, it is linear in the fraction of fcc-like clusters, taking values between those of least-ordered hcp and most-ordered fcc. This implies that the $\bar{\Lambda}$ of all Barlow packings is primarily controlled by the local cluster geometry. These results indicate the special nature of anisotropic stacking disorder, which provides impetus for future research on the development of anisotropic order metrics and hyperuniformity properties.Comment: 13 pages, 7 figure