BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Frank H. Stillinger

Frank H. Stillinger appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.stat-mechcond-mat.softcond-mat.mtrl-scimath.MGmath-phmath.MPphysics.atm-clus

Catalog footprint

What is connected

21works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

21 published item(s)

preprint2022arXiv

Characterization of Void Space, Large-Scale Structure, and Transport Properties of Maximally Random Jammed Packings of Superballs

Dense, disordered packings of particles are useful models of low-temperature amorphous phases of matter, biological systems, granular media, and colloidal systems. The study of dense packings of nonspherical particles enables one to ascertain how rotational degrees of freedom affect packing behavior. Here, we study superballs, a large family of deformations of the sphere, defined in three dimensions by $|x_1|^{2p}+|x_2|^{2p}+|x_3|^{2p}\leq1$, where $p>0$ is a deformation parameter indicating to what extent the shape deviates from a sphere. As $p$ increases from the sphere point ($p=1$), the superball attains cubic symmetry, and attains octahedral symmetry when $p<1$. Previous characterization of superball packings has shown that they have a maximally random jammed (MRJ) state, whose properties (e.g., packing fraction $ϕ$) vary nonanalytically as $p$ diverges from unity. Here, we use an event-driven molecular dynamics algorithm to produce MRJ superball packings. We characterize their large-scale structure by examining the small-$Q$ behavior of their structure factors and spectral densities, which indicate these packings are effectively hyperuniform. We show that the mean width $\bar{w}$ is a useful length scale to make distances dimensionless in order to compare systematically superballs of different shape. We also compute the complementary cumulative pore-size distribution $F(δ)$ and find the pore sizes tend to decrease as $|1-p|$ increases. From $F(δ)$, we estimate the fluid permeability, mean survival time, and principal diffusion relaxation time of the packings. Additionally, we compute the diffusion spreadability and find the long-time power-law scaling indicates these packings are hyperuniform. Our results can be used to help inform the design of granular materials with targeted densities and transport properties.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Static Structural Signatures of Nearly Jammed Disordered and Ordered Hard-Sphere Packings: Direct Correlation Function

Dynamical signatures are known to precede jamming in hard-particle systems, but static structural signatures have proven more elusive. The observation that compressing hard-particle packings towards jamming causes growing hyperuniformity has paved the way for the analysis of jamming as an "inverted critical point" in which the direct correlation function $c(r)$ diverges. We establish quantitative relationships between various singularities in $c(r)$ and the total correlation function $h(r)$ that provide a concrete means of identifying features that must be expressed in $c(r)$ if one hopes to reproduce details in the pair correlation function accurately. We also analyze 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) and Torquato-Jiao (TJ) packing algorithms, which both generate disordered packings, but can show perceptible structural differences. We identify a short-ranged scaling $c(r) \propto -1/r$ as $r \rightarrow 0$ and show that this, along with the developing delta function at $c(D)$, is a consequence of the growing long-rangedness in $c(r)$. Near the freezing density, we identify qualitative differences in the structure factor $S(k)$ as well as $c(r)$ between TJ- and LS-generated configurations and link them to differences in the protocols' packing dynamics. Configurations from both algorithms have structure factors that approach zero in the low-wavenumber limit as jamming is approached and are shown to exhibit a corresponding power-law decay in $c(r)$ for large $r$ as a consequence. Our work advances the notion that static signatures are exhibited by hard-particle packings as they approach jamming and underscores the utility of the direct correlation function as a means of monitoring for an incipient rigid network.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Superior Descriptor of Random Textures and Its Predictive Capacity

Two-phase random textures abound in a host of contexts, porous and composite media, ecological structures, biological media and astrophysical structures. Questions surrounding the spatial structure of such textures continue to pose many theoretical challenges. For example, can two-point correlation functions be identified that can be both manageably measured and yet reflect nontrivial higher-order structural information about the textures? We present a novel solution to this question by probing the information content of the widest class of different types of two-point functions examined to date using inverse "reconstruction" techniques. This enables us to show that a superior descriptor is the two-point cluster function $C_2({\bf r})$, which is sensitive to topological {\it connectedness} information. We demonstrate the utility of $C_2({\bf r})$ by accurately reconstructing textures drawn from materials science, cosmology and granular media, among other examples. Our work suggests an entirely new theoretical pathway to predict the bulk physical properties of random textures, and also has important ramifications for atomic and molecular systems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Densest binary sphere packings

The densest binary sphere packings in the alpha-x plane of small to large sphere radius ratio alpha 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 omega), 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 [S. Torquato and Y. 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 [A. B. Hopkins, Y. Jiao, F. H. Stillinger, and S. Torquato, 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 novel 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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Designed Diamond Ground State via Optimized Isotropic Monotonic Pair Potentials

We apply inverse statistical-mechanical methods to find a simple family of optimized isotropic, monotonic pair potentials, under certain constraints, whose ground states for a wide range of pressures is the diamond crystal. These constraints include desirable phonon spectra and the widest possible pressure range for stability. We also ascertain the ground-state phase diagram for a specific optimized potential to show that other crystal structures arise for other pressures. Cooling disordered configurations interacting with our optimized potential to absolute zero frequently leads to the desired diamond crystal ground state, revealing that the capture basin for the global energy minimum is large and broad relative to the local energy minima basins.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Exotic Ground States of Directional Pair Potentials via Collective-Density Variables

Collective-density variables have proved to be a useful tool in the prediction and manipulation of how spatial patterns form in the classical many-body problem. Previous work has employed properties of collective-density variables along with a robust numerical optimization technique to find the classical ground states of many-particle systems subject to radial pair potentials in one, two and three dimensions. That work led to the identification of ordered and disordered classical ground states. In this paper, we extend these collective-coordinate studies by investigating the ground states of directional pair potentials in two dimensions. Our study focuses on directional potentials whose Fourier representations are non-zero on compact sets that are symmetric with respect to the origin and zero everywhere else. We choose to focus on one representative set which has exotic ground-state properties: two circles whose centers are separated by some fixed distance. We obtain ground states for this "two-circle" potential that display large void regions in the disordered regime. As more degrees of freedom are constrained the ground states exhibit a collapse of dimensionality characterized by the emergence of filamentary structures and linear chains. This collapse of dimensionality has not been observed before in related studies.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Nonequilibrium static diverging length scales on approaching a prototypical model glassy state

Maximally random jammed states of hard spheres are prototypical glasses. We study the small wavenumber $k$ behavior of the structure factor $S(k)$ of overcompressed million-sphere packings as a function of density up to the jammed state. We find both a precursor to the glassy jammed state evident long before the jamming density is reached and an associated growing length scale, extracted from the volume integral of the direct correlation function $c(r)$, which diverges at the "critical" jammed state. We also define a nonequilibrium index $X$ and use it to demonstrate that the packings studied are intrinsically nonequilibrium in nature well before the critical point is reached.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Nonequilibrium static growing length scales in supercooled liquids on approaching the glass transition

The small wavenumber $k$ behavior of the structure factor $S(k)$ of overcompressed amorphous hard-sphere configurations was previously studied for a wide range of densities up to the maximally random jammed state, which can be viewed as a prototypical glassy state [A. Hopkins, F. H. Stillinger and S. Torquato, Phys. Rev. E, 86, 021505 (2012)]. It was found that a precursor to the glassy jammed state was evident long before the jamming density was reached as measured by a growing nonequilibrium length scale extracted from the volume integral of the direct correlation function $c(r)$, which becomes long-ranged as the critical jammed state is reached. The present study extends that work by investigating via computer simulations two different atomic models: the single-component Z2 Dzugutov potential in three dimensions and the binary-mixture Kob-Andersen potential in two dimensions. Consistent with the aforementioned hard-sphere study, we demonstrate that for both models a signature of the glass transition is apparent well before the transition temperature is reached as measured by the length scale determined from from the volume integral of the direct correlation function in the single-component case and a generalized direct correlation function in the binary-mixture case. The latter quantity is obtained from a generalized Orstein-Zernike integral equation for a certain decoration of the atomic point configuration. We also show that these growing length scales, which are a consequence of the long-range nature of the direct correlation functions, are intrinsically nonequilibrium in nature as determined by an index $X$ that is a measure of deviation from thermal equilibrium. It is also demonstrated that this nonequilibrium index, which increases upon supercooling, is correlated with a characteristic relaxation time scale.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Inherent Structures for Soft Long-Range Interactions in Two-Dimensional Many-Particle Systems

We generate inherent structures, local potential-energy minima, of the "$k$-space overlap potential" in two-dimensional many-particle systems using a cooling and quenching simulation technique. The ground states associated with the $k$-space overlap potential are stealthy ({\it i.e.,} completely suppress single scattering of radiation for a range of wavelengths) and hyperuniform ({\it i.e.,} infinite wavelength density fluctuations vanish). However, we show via quantitative metrics that the inherent structures exhibit a range of stealthiness and hyperuniformity depending on the fraction of degrees of freedom that are constrained. Inherent structures in two dimensions typically contain five-particle rings, wavy grain boundaries, and vacancy-interstitial defects. The structural and thermodynamic properties of inherent structures are relatively insensitive to the temperature from which they are sampled, signifying that the energy landscape is relatively flat and devoid of deep wells. Using the nudged-elastic-band algorithm, we construct paths from ground-state configurations to inherent structures and identify the transition points between them. In addition, we use point patterns generated from a random sequential addition (RSA) of hard disks, which are nearly stealthy, and examine the particle rearrangements necessary to make the configurations absolutely stealthy. We introduce a configurational proximity metric to show that only small local, but collective, particle rearrangements are needed to drive initial RSA configurations to stealthy disordered ground states. These results lead to a more complete understanding of the unusual behaviors exhibited by the family of "collective-coordinate" potentials to which the $k$-space overlap potential belongs.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Non-Universality of Density and Disorder in Jammed Sphere Packings

We show for the first time that collectively jammed disordered packings of three-dimensional monodisperse frictionless hard spheres can be produced and tuned using a novel numerical protocol with packing density $ϕ$ as low as 0.6. This is well below the value of 0.64 associated with the maximally random jammed state and entirely unrelated to the ill-defined ``random loose packing'' state density. Specifically, collectively jammed packings are generated with a very narrow distribution centered at any density $ϕ$ over a wide density range $ϕ\in [0.6,~0.74048\ldots]$ with variable disorder. Our results support the view that there is no universal jamming point that is distinguishable based on the packing density and frequency of occurence. Our jammed packings are mapped onto a density-order-metric plane, which provides a broader characterization of packings than density alone. Other packing characteristics, such as the pair correlation function, average contact number and fraction of rattlers are quantified and discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Phase diagram and structural diversity of the densest binary sphere packings

The densest binary sphere packings have historically been very difficult to determine. The only rigorously known packings in the alpha-x plane of sphere radius ratio alpha and relative concentration x are at the Kepler limit alpha = 1, where packings are monodisperse. Utilizing an implementation of the Torquato-Jiao sphere-packing algorithm [S. Torquato and Y. Jiao, Phys. Rev. E 82, 061302 (2010)], we present the most comprehensive determination to date of the phase diagram in (alpha,x) for the densest binary sphere packings. Unexpectedly, we find many distinct new densest packings.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Densest local packing diversity. II. Application to three dimensions

The densest local packings of N three-dimensional identical nonoverlapping spheres within a radius Rmin(N) of a fixed central sphere of the same size are obtained for selected values of N up to N = 1054. In the predecessor to this paper [A.B. Hopkins, F.H. Stillinger and S. Torquato, Phys. Rev. E 81 041305 (2010)], we described our method for finding the putative densest packings of N spheres in d-dimensional Euclidean space Rd and presented those packings in R2 for values of N up to N = 348. We analyze the properties and characteristics of the densest local packings in R3 and employ knowledge of the Rmin(N), using methods applicable in any d, to construct both a realizability condition for pair correlation functions of sphere packings and an upper bound on the maximal density of infinite sphere packings. In R3, we find wide variability in the densest local packings, including a multitude of packing symmetries such as perfect tetrahedral and imperfect icosahedral symmetry. We compare the densest local packings of N spheres near a central sphere to minimal-energy configurations of N+1 points interacting with short-range repulsive and long-range attractive pair potentials, e.g., 12-6 Lennard-Jones, and find that they are in general completely different, a result that has possible implications for nucleation theory. We also compare the densest local packings to finite subsets of stacking variants of the densest infinite packings in R3 (the Barlow packings) and find that the densest local packings are almost always most similar, as measured by a similarity metric, to the subsets of Barlow packings with the smallest number of coordination shells measured about a single central sphere, e.g., a subset of the FCC Barlow packing. We additionally observe that the densest local packings are dominated by the spheres arranged with centers at precisely distance Rmin(N) from the fixed sphere's center.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Duality Relations for the Classical Ground States of Soft-Matter Systems

Bounded interactions are particularly important in soft-matter systems, such as colloids, microemulsions, and polymers. We derive new duality relations for a class of soft potentials, including three-body and higher-order functions, that can be applied to ordered and disordered classical ground states. These duality relations link the energy of configurations associated with a real-space potential to the corresponding energy of the dual (Fourier-transformed) potential. We apply the duality relations by demonstrating how information about the classical ground states of short-ranged potentials can be used to draw new conclusions about the ground states of long-ranged potentials and vice versa. The duality relations also lead to bounds on the T=0 system energies in density intervals of phase coexistence. Additionally, we identify classes of "self-similar" potentials, for which one can relate low- and high-density ground-state energies. We analyze the ground state configurations and thermodynamic properties of a one-dimensional system previously thought to exhibit an infinite number of structural phase transitions and comment on the known ground states of purely repulsive monotonic potentials in the context of our duality relations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Geometrical Ambiguity of Pair Statistics. II. Heterogeneous Media

In a previous paper [Jiao, Stillinger and Torquato, PRE 81, 011105 (2010)], we considered the geometrical ambiguity of pair statistics associated with point configurations. Here we focus on the analogous problem for heterogeneous media (materials). The complex structures of heterogeneous media are usually characterized via statistical descriptors, such as the $n$-point correlation function $S_n$. An intricate inverse problem of practical importance is to what extent a medium can be reconstructed from the two-point correlation function $S_2$ of a target medium. Recently, general claims of the uniqueness of reconstructions using $S_2$ have been made based on numerical studies. Here, we provide a systematic approach to characterize the geometrical ambiguity of $S_2$ for both continuous two-phase heterogeneous media and their digitized representations in a mathematically precise way. In particular, we derive the exact conditions for the case where two distinct media possess identical $S_2$. The degeneracy conditions are given in terms of integral and algebraic equations for continuous media and their digitized representations, respectively. By examining these equations and constructing their rigorous solutions for specific examples, we conclusively show that in general $S_2$ is indeed not sufficient information to uniquely determine the structure of the medium, which is consistent with the results of our recent study on heterogeneous media reconstruction [Jiao, Stillinger and Torquato, PNAS 106, 17634 (2009)]. The uniqueness issue of multiphase media reconstructions and additional structural information required to characterize heterogeneous media are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Jammed Hard-Particle Packings: From Kepler to Bernal and Beyond

This review describes the diversity of jammed configurations attainable by frictionless convex nonoverlapping (hard) particles in Euclidean spaces and for that purpose it stresses individual-packing geometric analysis. A fundamental feature of that diversity is the necessity to classify individual jammed configurations according to whether they are locally, collectively, or strictly jammed. Each of these categories contains a multitude of jammed configurations spanning a wide and (in the large system limit) continuous range of intensive properties, including packing fraction $ϕ$, mean contact number $Z$, and several scalar order metrics. Application of these analytical tools to spheres in three dimensions (an analog to the venerable Ising model) covers a myriad of jammed states, including maximally dense packings (as Kepler conjectured), low-density strictly-jammed tunneled crystals, and a substantial family of amorphous packings. With respect to the last of these, the current approach displaces the historically prominent but ambiguous idea of ``random close packing" (RCP) with the precise concept of ``maximally random jamming" (MRJ). This review also covers recent advances in understanding jammed packings of polydisperse sphere mixtures, as well as convex nonspherical particles, e.g., ellipsoids, ``superballs", and polyhedra. Because of their relevance to error-correcting codes and information theory, sphere packings in high-dimensional Euclidean spaces have been included as well. We also make some remarks about packings in (curved) non-Euclidean spaces. In closing this review, several basic open questions for future research to consider have been identified.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Novel Ground-State Crystals with Controlled Vacancy Concentrations: From Kagomé to Honeycomb to Stripes

We introduce a one-parameter family, $0 \leq H \leq 1$, of pair potential functions with a single relative energy minimum that stabilize a range of vacancy-riddled crystals as ground states. The "quintic potential" is a short-ranged, nonnegative pair potential with a single local minimum of height $H$ at unit distance and vanishes cubically at a distance of $\rt$. We have developed this potential to produce ground states with the symmetry of the triangular lattice while favoring the presence of vacancies. After an exhaustive search using various optimization and simulation methods, we believe that we have determined the ground states for all pressures, densities, and $0 \leq H \leq 1$. For specific areas below $3\rt/2$, the ground states of the "quintic potential" include high-density and low-density triangular lattices, kagomé and honeycomb crystals, and stripes. We find that these ground states are mechanically stable but are difficult to self-assemble in computer simulations without defects. For specific areas above $3\rt/2$, these systems have a ground-state phase diagram that corresponds to hard disks with radius $\rt$. For the special case of H=0, a broad range of ground states is available. Analysis of this case suggests that among many ground states, a high-density triangular lattice, low-density triangular lattice, and striped phases have the highest entropy for certain densities. The simplicity of this potential makes it an attractive candidate for experimental realization with application to the development of novel colloidal crystals or photonic materials.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Optimized Monotonic Convex Pair Potentials Stabilize Low-Coordinated Crystals

We have previously used inverse statistical-mechanical methods to optimize isotropic pair interactions with multiple extrema to yield low-coordinated crystal classical ground states (e.g., honeycomb and diamond structures) in d-dimensional Euclidean space R^d. Here we demonstrate the counterintuitive result that no extrema are required to produce such low-coordinated classical ground states. Specifically, we show that monotonic convex pair potentials can be optimized to yield classical ground states that are the square and honeycomb crystals in R^2 over a non-zero number density range. Such interactions may be feasible to achieve experimentally using colloids and polymers.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Phase Behavior of Colloidal Superballs: Shape Interpolation from Spheres to Cubes

The phase behavior of hard superballs is examined using molecular dynamics within a deformable periodic simulation box. A superball's interior is defined by the inequality $|x|^{2q} + |y|^{2q} + |z|^{2q} \leq 1$, which provides a versatile family of convex particles ($q \geq 0.5$) with cube-like and octahedron-like shapes as well as concave particles ($q < 0.5$) with octahedron-like shapes. Here, we consider the convex case with a deformation parameter q between the sphere point (q = 1) and the cube (q = 1). We find that the asphericity plays a significant role in the extent of cubatic ordering of both the liquid and crystal phases. Calculation of the first few virial coefficients shows that superballs that are visually similar to cubes can have low-density equations of state closer to spheres than to cubes. Dense liquids of superballs display cubatic orientational order that extends over several particle lengths only for large q. Along the ordered, high-density equation of state, superballs with 1 < q < 3 exhibit clear evidence of a phase transition from a crystal state to a state with reduced long-ranged orientational order upon the reduction of density. For $q \geq 3$, long-ranged orientational order persists until the melting transition. The width of coexistence region between the liquid and ordered, high-density phase decreases with q up to q = 4.0. The structures of the high-density phases are examined using certain order parameters, distribution functions, and orientational correlation functions. We also find that a fixed simulation cell induces artificial phase transitions that are out of equilibrium. Current fabrication techniques allow for the synthesis of colloidal superballs, and thus the phase behavior of such systems can be investigated experimentally.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Superdense Crystal Packings of Ellipsoids

Particle packing problems have fascinated people since the dawn of civilization, and continue to intrigue mathematicians and scientists. Resurgent interest has been spurred by the recent proof of Kepler's conjecture: the face-centered cubic lattice provides the densest packing of equal spheres with a packing fraction $ϕ\approx0.7405$ \cite{Kepler_Hales}. Here we report on the densest known packings of congruent ellipsoids. The family of new packings are crystal (periodic) arrangements of nearly spherically-shaped ellipsoids, and always surpass the densest lattice packing. A remarkable maximum density of $ϕ\approx0.7707$ is achieved for both prolate and oblate ellipsoids with aspect ratios of $\sqrt{3}$ and $1/\sqrt{3}$, respectively, and each ellipsoid has 14 touching neighbors. Present results do not exclude the possibility that even denser crystal packings of ellipsoids could be found, and that a corresponding Kepler-like conjecture could be formulated for ellipsoids.

0 citations0 reviewsNo code linkedOpen access