作者:
Adam B. HopkinsFrank H. StillingerSalvatore TorquatoDepartment of Chemistry
Department of Physics Princeton Center for Theoretical Science Princeton Institute for the Science and Technology of Materials Program in Applied and Computational Mathematics Princeton University Princeton New Jersey 08544 USA
The densest binary sphere packings in the α-x plane of small to large sphere radius ratio α and small sphere relative concentration x have historically been very difficult to determine. Previous research had led to ...
详细信息
The densest binary sphere packings in the α-x plane of small to large sphere radius ratio α and small sphere relative concentration x have historically been very difficult to determine. Previous research had led to the prediction that these packings were composed of a few known “alloy” phases including, for example, the AlB2 (hexagonal ω), HgBr2, and AuTe2 structures, and to XYn structures composed of close-packed large spheres with small spheres (in a number ratio of n to 1) in the interstices, e.g., the NaCl packing for n=1. However, utilizing an implementation of the Torquato-Jiao sphere-packing algorithm [Torquato and Jiao, Phys. Rev. E 82, 061302 (2010)], we have discovered that many more structures appear in the densest packings. For example, while all previously known densest structures were composed of spheres in small to large number ratios of one to one, two to one, and very recently three to one, we have identified densest structures with number ratios of seven to three and five to two. In a recent work [Hopkins et al., Phys. Rev. Lett. 107, 125501 (2011)], we summarized these findings. In this work, we present the structures of the densest-known packings and provide details about their characteristics. Our findings demonstrate that a broad array of different densest mechanically stable structures consisting of only two types of components can form without any consideration of attractive or anisotropic interactions. In addition, the structures that we have identified may correspond to currently unidentified stable phases of certain binary atomic and molecular systems, particularly at high temperatures and pressures.
作者:
Salvatore TorquatoYang Jiao[]Department of Chemistry
Department of Physics Princeton Center for Theoretical Science Program of Applied and Computational Mathematics Princeton Institute of the Science and Technology of Materials Princeton University Princeton New Jersey 08544 USA
We have recently devised organizing principles to obtain maximally dense packings of the Platonic and Archimedean solids and certain smoothly shaped convex nonspherical particles [Torquato and Jiao, Phys. Rev. E 81, 0...
详细信息
We have recently devised organizing principles to obtain maximally dense packings of the Platonic and Archimedean solids and certain smoothly shaped convex nonspherical particles [Torquato and Jiao, Phys. Rev. E 81, 041310 (2010)]. Here we generalize them in order to guide one to ascertain the densest packings of other convex nonspherical particles as well as concave shapes. Our generalized organizing principles are explicitly stated as four distinct propositions. All of our organizing principles are applied to and tested against the most comprehensive set of both convex and concave particle shapes examined to date, including Catalan solids, prisms, antiprisms, cylinders, dimers of spheres, and various concave polyhedra. We demonstrate that all of the densest known packings associated with this wide spectrum of nonspherical particles are consistent with our propositions. Among other applications, our general organizing principles enable us to construct analytically the densest known packings of certain convex nonspherical particles, including spherocylinders, “lens-shaped” particles, square pyramids, and rhombic pyramids. Moreover, we show how to apply these principles to infer the high-density equilibrium crystalline phases of hard convex and concave particles. We also discuss the unique packing attributes of maximally random jammed packings of nonspherical particles.
Hyperuniformity characterizes a state of matter that is poised at a critical point at which density or volume-fraction fluctuations are anomalously suppressed at infinite wavelengths. Recently, much attention has been...
详细信息
Hyperuniformity characterizes a state of matter that is poised at a critical point at which density or volume-fraction fluctuations are anomalously suppressed at infinite wavelengths. Recently, much attention has been given to the link between strict jamming (mechanical rigidity) and (effective or exact) hyperuniformity in frictionless hard-particle packings. However, in doing so, one must necessarily study very large packings in order to access the long-ranged behavior and to ensure that the packings are truly jammed. We modify the rigorous linear programming method of Donev et al. [J. Comput. Phys. 197, 139 (2004)] in order to test for jamming in putatively collectively and strictly jammed packings of hard disks in two dimensions. We show that this rigorous jamming test is superior to standard ways to ascertain jamming, including the so-called “pressure-leak” test. We find that various standard packing protocols struggle to reliably create packings that are jammed for even modest system sizes of N≈103 bidisperse disks in two dimensions; importantly, these packings have a high reduced pressure that persists over extended amounts of time, meaning that they appear to be jammed by conventional tests, though rigorous jamming tests reveal that they are not. We present evidence that suggests that deviations from hyperuniformity in putative maximally random jammed (MRJ) packings can in part be explained by a shortcoming of the numerical protocols to generate exactly jammed configurations as a result of a type of “critical slowing down” as the packing's collective rearrangements in configuration space become locally confined by high-dimensional “bottlenecks” from which escape is a rare event. Additionally, various protocols are able to produce packings exhibiting hyperuniformity to different extents, but this is because certain protocols are better able to approach exactly jammed configurations. Nonetheless, while one should not generally expect exact hyperuniformity for d
作者:
Steven AtkinsonYang JiaoSalvatore TorquatoDepartment of Chemistry
Department of Physics Princeton Center for Theoretical Science Program of Applied and Computational Mathematics Princeton Institute for the Science and Technology of Materials Princeton University Princeton New Jersey 08544 USA
Dense packings of hard particles have important applications in many fields, including condensed matter physics, discrete geometry, and cell biology. In this paper, we employ a stochastic search implementation of the ...
详细信息
Dense packings of hard particles have important applications in many fields, including condensed matter physics, discrete geometry, and cell biology. In this paper, we employ a stochastic search implementation of the Torquato-Jiao adaptive-shrinking-cell (ASC) optimization scheme [Nature (London) 460, 876 (2009)] to find maximally dense particle packings in d-dimensional Euclidean space Rd. While the original implementation was designed to study spheres and convex polyhedra in d≥3, our implementation focuses on d=2 and extends the algorithm to include both concave polygons and certain complex convex or concave nonpolygonal particle shapes. We verify the robustness of this packing protocol by successfully reproducing the known putative optimal packings of congruent copies of regular pentagons and octagons, then employ it to suggest dense packing arrangements of congruent copies of certain families of concave crosses, convex and concave curved triangles (incorporating shapes resembling the Mercedes-Benz logo), and “moonlike” shapes. Analytical constructions are determined subsequently to obtain the densest known packings of these particle shapes. For the examples considered, we find that the densest packings of both convex and concave particles with central symmetry are achieved by their corresponding optimal Bravais lattice packings; for particles lacking central symmetry, the densest packings obtained are nonlattice periodic packings, which are consistent with recently-proposed general organizing principles for hard particles. Moreover, we find that the densest known packings of certain curved triangles are periodic with a four-particle basis, and we find that the densest known periodic packings of certain moonlike shapes possess no inherent symmetries. Our work adds to the growing evidence that particle shape can be used as a tuning parameter to achieve a diversity of packing structures.
作者:
Steven AtkinsonFrank H. StillingerSalvatore TorquatoDepartment of Chemistry
Department of Physics Princeton Center for Theoretical Science Program of Applied and Computational Mathematics Princeton Institute for the Science and Technology of Materials Princeton University Princeton New Jersey 08544 USA
The nonequilibrium process by which hard-particle systems may be compressed into disordered, jammed states has received much attention because of its wide utility in describing a broad class of amorphous materials. Wh...
详细信息
The nonequilibrium process by which hard-particle systems may be compressed into disordered, jammed states has received much attention because of its wide utility in describing a broad class of amorphous materials. While dynamical signatures are known to precede jamming, the task of identifying static structural signatures indicating the onset of jamming have proven more elusive. The observation that compressing hard-particle packings towards jamming is accompanied by an anomalous suppression of density fluctuations (termed “hyperuniformity”) has paved the way for the analysis of jamming as an “inverted critical point” in which the direct correlation function c(r), rather than the total correlation function h(r), diverges. We expand on the notion that c(r) provides both universal and protocol-specific information as packings approach jamming. By considering the degree and position of singularities (discontinuities in the nth derivative) as well as how they are changed by the convolutions found in the Ornstein-Zernike equation, we establish quantitative statements about the structure of c(r) with regards to singularities it inherits from h(r). These relations provide a concrete means of identifying features that must be expressed in c(r) if one hopes to reproduce various details in the pair correlation function accurately and provide stringent tests on the associated numerics. We also analyze the evolution of systems of three-dimensional monodisperse hard spheres of diameter D as they approach ordered and disordered jammed configurations. For the latter, we use the Lubachevsky-Stillinger (LS) molecular dynamics and Torquato-Jiao (TJ) sequential linear programming algorithms, which both generate disordered packings, but can show perceptible structural differences. We identify a short-ranged scaling c(r)∝−1/r as r→0 that accompanies the formation of the delta function at c(D) that indicates the formation of contacts in all cases, and show that this scaling behavior is,
ABACUS (Atomic-orbital Based Ab-initio Computation at USTC) is an open-source software for first-principles electronic structure calculations and molecular dynamics simulations. It mainly features density functional t...
详细信息
Using the Deep Potential methodology, we construct a model that reproduces accurately the potential energy surface of the SCAN approximation of density functional theory for water, from low temperature and pressure to...
详细信息
This paper presents a new Lagrangian type scheme for solving the Euler equations of compressible gas *** this new scheme the system of equations is discretized by Runge-Kutta Discontinuous Galerkin(RKDG)method,and the...
详细信息
This paper presents a new Lagrangian type scheme for solving the Euler equations of compressible gas *** this new scheme the system of equations is discretized by Runge-Kutta Discontinuous Galerkin(RKDG)method,and the mesh moves with the fluid *** scheme is conservative for the mass,momentum and total energy and maintains second-order *** scheme avoids solving the geometrical part and has free *** of some numerical tests are presented to demonstrate the accuracy and the non-oscillatory property of the scheme.
We find a fraction of atoms remain unionized after the laser pulse when the tunneled electrons are released in a certain window of initial field phase and transverse velocity. The survival window shifts with laser pol...
详细信息
We find a fraction of atoms remain unionized after the laser pulse when the tunneled electrons are released in a certain window of initial field phase and transverse velocity. The survival window shifts with laser polarization ellipticity and its width varies with respect to laser intensity and atomic ionization potential. Neutral atom yield can be calculated by summing up tunneling probabilities in the window. Our theory can quantitatively reproduce the distribution of the survival yields vs laser ellipticity observed for helium in experiment. For other atom species with smaller ionization potential such as magnesium, our theory predicts a wider distribution than the strong-field approximation model while closer to the three-dimensional semiclassical electron ensemble simulations, indicating the important role of the Coulomb effects.
To overcome the problem of inefficient computing time and unreliable results in MCNP5 calculation, a two-step method is adopted to calculate the energy deposition of prompt γ-rays in detectors for depleted uranium sp...
详细信息
To overcome the problem of inefficient computing time and unreliable results in MCNP5 calculation, a two-step method is adopted to calculate the energy deposition of prompt γ-rays in detectors for depleted uranium spherical shells under D-T neutron irradiation. In the first step, the γ-ray spectrum for energy below 7 MeV is calculated by MCNP5 code; secondly, the electron recoil spectrum in a BC501 A liquid scintillator detector is simulated based on EGSnrc Monte Carlo Code with the γ-ray spectrum from the first step as input. The comparison of calculated results with experimental ones shows that the simulations agree well with experiment in the energy region 0.4–3 MeV for the prompt γ-ray spectrum and below 4 MeVee for the electron recoil spectrum. The reliability of the two-step method in this work is validated.
暂无评论