Inherent structure landscape connection between liquids, granular materials, and the jamming phase diagram

Thermodynamic properties of quasi-one-dimensional fluids

We calculate thermodynamic and structural quantities of a fluid of hard spheres of diameter σ in a quasi-one-dimensional pore with accessible pore width W smaller than σ by applying a perturbative method worked out earlier for a confined fluid in a slit pore [Franosch et al. Phys. Rev. Lett. 109, 240601 (2012)]. In a first step, we prove that the thermodynamic and a certain class of structural quantities of the hard-sphere fluid in the pore can be obtained from a purely one-dimensional fluid of rods of length σ with a central hard core of size σW=σ2-W2 and a soft part at both ends of length (σ - σW)/2. These rods interact via effective k-body potentials veff(k) (k ≥ 2). The two- and the three-body potential will be calculated explicitly. In a second step, the free energy of this effective one-dimensional fluid is calculated up to leading order in (W/σ)2. Explicit results for, e.g., the perpendicular pressure, surface tension, and the density profile as a function of density, temperature, and pore width are presented and partly compared with results from Monte-Carlo simulations and standard virial expansions. Despite the perturbative character of our approach, it encompasses the singularity of the thermodynamic quantities at the jamming transition point.

Structural Transformation between a Nematic Loose Packing and a Randomly Stacked Close Packing of Granular Disks

Packing structures of granular disks are reconstructed using magnetic resonance imaging techniques. As packing fraction increases, the packing structure transforms from a nematic loose packing to a dense packing with randomly oriented stacks. According to our model based on Edwards' volume ensemble, stack structures are statistically favored when the effective temperature decreases, which has a lower structural anisotropy than single disks, and brings down the global orientational order consequently. This mechanism identified in athermal granular materials can help us understand the nonergodic characteristics of disklike particle assemblies such as discotic mesogens and clays.

Is there a granular potential?

Granular materials, such as sand or grain, exhibit many structural and dynamic characteristics similar to those observed in molecular systems, despite temperature playing no role in their properties. This has led to an effort to develop a statistical mechanics for granular materials that has focused on establishing an equivalent to the microcanonical ensemble and a temperature-like thermodynamic variable. Here, we expand on these ideas by introducing a granular potential into the Edwards ensemble, as an analogue to the chemical potential, and explore its properties using a simple model of a granular system. A simple kinetic Monte Carlo simulation of the model shows the effect of mass transport leading to equilibrium and how this is connected to the redistribution of volume in the system. An exact analytical treatment of the model shows that the compactivity and the ratio of the granular potential to the compactivity determine the equilibrium between two open systems that are able to exchange volume and particles, and that mass moves from high to low values of this ratio. Analysis of the granular potential shows that adding a particle to the system increases the entropy at high compactivity, but decreases the entropy at low compactivity. Finally, we demonstrate the use of a small system thermodynamics method for the measurement of granular potential differences.

Elastogranular columns and beams

String and grains can be combined to create structures capable of bearing significant loads. In this work, we prepare columns and beams through a layer-by-layer deposition of granular matter and loops of fiber strings, and characterize their mechanical properties. The loops cause the grains to jam, and the inter-grain contact leads to a Hertzian-like constitutive response. Initially, one force chain that propagates vertically through the column bears most of the compressive load. As the magnitude of the load is increased, more force chains form in the column, which act in parallel to increase its stiffness, akin to a "super-Hertzian" regime. Applying a compressive prestress enables the structures to withstand shear, enabling the fabrication of cantilevered beams. This work provides a mechanical framework to use elastogranular jamming to create rapid, reusable infrastructure components, such as columns, beams, and arches from inexpensive, commonplace materials, such as rocks and string.

Jammed disks of two sizes and weights in a channel: Alternating sequences

Disks of two sizes and weights in alternating sequence are confined to a long and narrow channel. The axis of the channel is horizontal and its plane vertical. The channel is closed off by pistons that freeze jammed microstates out of loose disk configurations subject to moderate pressure, gravity, and random agitations. Disk sizes and channel width are such that under jamming no disk remains loose and all disks touch one wall. We present exact results for the characterization of jammed macrostates including volume and entropy. The rigorous analysis divides the disk sequences of jammed microstates into overlapping tiles from which we construct a small number of species of statistically interacting particles. Jammed macrostates depend on dimensionless control parameters inferred from ratios between measures of expansion work against the pistons, gravitational potential energy, and intensity of random agitations. These control parameters enter the configurational statistics via the activation energies prior to jamming of the particles. The range of disk weights naturally divides into regimes where qualitatively different features come into play. We sketch a path toward generalizations that include random sequences under a modified jamming protocol.

Inherent structure landscape of hard spheres confined to narrow cylindrical channels

The inherent structure landscape for a system of hard spheres confined to a hard cylindrical channel, such that spheres can only contact their first and second neighbors, is studied using an analytical model that extends previous results [Phys. Rev. Lett. 115, 025702 (2015)PRLTAO0031-900710.1103/PhysRevLett.115.025702] to provide a comprehensive picture of jammed packings over a range of packing densities. In the model, a packing is described as an arrangement of k helical sections, separated by defects, that have alternating helical twist directions and where all spheres satisfy local jamming constraints. The structure of each helical section is determined by a single helical twist angle, and a jammed packing is obtained by minimizing the length of the channel per particle with respect to the k helical section angles. An analysis of a small system of N=20 spheres shows that the basins on the inherent structure landscape associated with these helical arrangements split into a number of distinct jammed states separated by low barriers giving rise to a degree of hierarchical organization. The model accurately predicts the geometric properties of packings generated using the Lubachevsky and Stillinger compression scheme (N=10^{4}) and provides insight into the nature of the probability distribution of helical section lengths.

Marginally jammed states of hard disks in a one-dimensional channel

We have studied a class of marginally jammed states in a system of hard disks confined in a narrow channel-a quasi-one-dimensional system-whose exponents are not those predicted by theories valid in the infinite dimensional limit. The exponent γ which describes the distribution of small gaps takes the value 1 rather than the infinite dimensional value 0.41269⋯. Our work shows that there exist jammed states not found within the tiling approach of Ashwin and Bowles. The most dense of these marginal states is an unusual state of matter that is asymptotically crystalline.

Buckling of a linear chain of hard spheres in a harmonic confining potential: Numerical and analytical results for low and high compression

We extend a previous analysis of the buckling properties of a linear chain of hard spheres between hard walls under transverse harmonic confinement. Two regimes are distinguished-low compression, for which the entire chain buckles, and higher compression, for which there is localized buckling. With further increase of compression, second-neighbor contacts occur; beyond this compression the structure is no longer planar, and is not treated here. A continuous model is developed which is amenable to analytical solution in the low compression regime. This is helpful in understanding the scaling properties of both finite and infinite chains.

Jamming Below Upper Critical Dimension

Extensive numerical simulations in the past decades proved that the critical exponents of the jamming of frictionless spherical particles are the same in two and three dimensions. This implies that the upper critical dimension is d_{u}=2 or lower. In this Letter, we study the jamming transition below the upper critical dimension. We investigate a quasi-one-dimensional system: disks confined in a narrow channel. We show that the system is isostatic at the jamming transition point as in the case of standard jamming transition of the bulk systems in two and three dimensions. Nevertheless, the scaling of the excess contact number shows the linear scaling. Furthermore, the gap distribution remains finite even at the jamming transition point. These results are qualitatively different from those of the bulk systems in two and three dimensions.

Mapping diffusivity of narrow channels into one dimension

The diffusion of particles trapped in long narrow channels occurs predominantly in one dimension. Here, a molecular-dynamics simulation is used to study the inertial dynamics of two-dimensional hard disks confined to long, narrow, structureless channels with hard walls in the no-passing regime. We show that the diffusion coefficient obtained from the mean-squared displacement can be mapped onto the exact results for the diffusion of the strictly-one-dimensional hard rod system through an effective occupied volume fraction obtained from either the equation of state or a geometric projection of the particle interaction diameters.

Gardner Transition in Physical Dimensions

The Gardner transition is the transition that at mean-field level separates a stable glass phase from a marginally stable phase. This transition has similarities with the de Almeida-Thouless transition of spin glasses. We have studied a well-understood problem, that of disks moving in a narrow channel, which shows many features usually associated with the Gardner transition. We show that some of these features are artifacts that arise when a disk escapes its local cage during the quench to higher densities. There is evidence that the Gardner transition becomes an avoided transition, in that the correlation length becomes quite large, of order 15 particle diameters, even in our quasi-one-dimensional system.

Positional ordering of hard adsorbate particles in tubular nanopores

The phase behavior and structural properties of a monolayer of hard particles is examined in such a confinement where the adsorbed particles are constrained to the surface of a narrow hard cylindrical pore. The diameter of the pore is chosen such that only first- and second-neighbor interactions occur between the hard particles. The transfer operator method of [Percus and Zhang, Mol. Phys. 69, 347 (1990)MOPHAM0026-897610.1080/00268979000100241] is reformulated to obtain information about the structure of the monolayer. We have found that a true phase transition is not possible in the examined range of pore diameters. The monolayer of hard spheres undergoes a structural change from fluidlike order to a zigzaglike solid one with increasing surface density. The case of hard cylinders is different in the sense that a layering takes place continuously between a low-density one-row and a high-density two-row monolayer. Our results reveal a clear discrepancy with classical density functional theories, which do not distinguish smecticlike ordering in bulk from that in narrow periodic pores.

Critical behavior of hard squares in strong confinement

We examine the phase behavior of a quasi-one-dimensional system of hard squares with side-length σ, where the particles are confined between two parallel walls and only nearest-neighbor interactions occur. As in our previous work [Gurin, Varga, and Odriozola, Phys. Rev. E 94, 050603 (2016)]2470-004510.1103/PhysRevE.94.050603, the transfer operator method is used, but here we impose a restricted orientation and position approximation to yield an analytic description of the physical properties. This allows us to study the parallel fluid-like to zigzag solid-like structural transition, where the compressibility and heat capacity peaks sharpen and get higher as H→H_{c}=2sqrt[2]-1≈1.8284 and p→p_{c}=∞. Here H is the width of the channel measured in σ units and p is the pressure. We have found that this structural change becomes critical at the (p_{c},H_{c}) point. The obtained critical exponents belong to the universality class of the one-dimensional Ising model. We believe this behavior holds for the unrestricted orientational and positional case.

Anomalous structural transition of confined hard squares

Structural transitions are examined in quasi-one-dimensional systems of freely rotating hard squares, which are confined between two parallel walls. We find two competing phases: one is a fluid where the squares have two sides parallel to the walls, while the second one is a solidlike structure with a zigzag arrangement of the squares. Using transfer matrix method we show that the configuration space consists of subspaces of fluidlike and solidlike phases, which are connected with low probability microstates of mixed structures. The existence of these connecting states makes the thermodynamic quantities continuous and precludes the possibility of a true phase transition. However, thermodynamic functions indicate strong tendency for the phase transition and our replica exchange Monte Carlo simulation study detects several important markers of the first order phase transition. The distinction of a phase transition from a structural change is practically impossible with simulations and experiments in such systems like the confined hard squares.

Phase behaviour and correlations of parallel hard squares: from highly confined to bulk systems

We study a fluid of two-dimensional parallel hard squares in bulk and under confinement in channels, with the aim of evaluating the performance of fundamental-measure theory (FMT). To this purpose, we first analyse the phase behaviour of the bulk system using FMT and Percus-Yevick (PY) theory, and compare the results with molecular dynamics and Monte Carlo simulations. In a second step, we study the confined system and check the results against those obtained from the transfer matrix method and from our own Monte Carlo simulations. Squares are confined to channels with parallel walls at angles of 0° or 45° relative to the diagonals of the parallel hard squares, respectively, which allows for an assessment of the effect of the external-potential symmetry on the fluid structural properties. In general FMT overestimates bulk correlations, predicting the existence of a columnar phase (absent in simulations) prior to crystallization. The equation of state predicted by FMT compares well with simulations, although the PY approach with the virial route is better in some range of packing fractions. The FMT is highly accurate for the structure and correlations of the confined fluid due to the dimensional crossover property fulfilled by the theory. Both density profiles and equations of state of the confined system are accurately predicted by the theory. The highly non-uniform pair correlations inside the channel are also very well described by FMT.

Glasslike behavior of a hard-disk fluid confined to a narrow channel

Disks moving in a narrow channel have many features in common with the glassy behavior of hard spheres in three dimensions. In this paper we study the caging behavior of the disks that sets in at characteristic packing fraction ϕ(d). Four-point overlap functions similar to those studied when investigating dynamical heterogeneities have been determined from event-driven molecular dynamics simulations and the time-dependent dynamical length scale has been extracted from them. The dynamical length scale increases with time and, on the equilibration time scale, it is proportional to the static length scale associated with the zigzag ordering in the system, which grows rapidly above ϕ(d). The structural features responsible for the onset of caging and the glassy behavior are easy to identify as they show up in the structure factor, which we have determined exactly from the transfer-matrix approach.

Understanding the ideal glass transition: lessons from an equilibrium study of hard disks in a channel

We use an exact transfer-matrix approach to compute the equilibrium properties of a system of hard disks of diameter σ confined to a two-dimensional channel of width 1.95σ at constant longitudinal applied force. At this channel width, which is sufficient for next-nearest-neighbor disks to interact, the system is known to have a great many jammed states. Our calculations show that the longitudinal force (pressure) extrapolates to infinity at a well-defined packing fraction ϕ(K) that is less than the maximum possible ϕ(max), the latter corresponding to a buckled crystal. In this quasi-one-dimensional problem there is no question of there being any real divergence of the pressure at ϕ(K). We give arguments that this avoided phase transition is a structural feature, the remnant in our narrow channel system of the hexatic to crystal transition, but that it has the phenomenology of the (avoided) ideal glass transition. We identify a length scale ξ̃(3) as our equivalent of the penetration length for amorphous order: In the channel system, it reaches a maximum value of around 15σ at ϕ(K), which is larger than the penetration lengths that have been reported for three-dimensional systems. It is argued that the α-relaxation time would appear on extrapolation to diverge in a Vogel-Fulcher manner as the packing fraction approaches ϕ(K).

Inherent structures, fragility, and jamming: insights from quasi-one-dimensional hard disks

We study a quasi-one-dimensional system of hard disks confined between hard lines to explore the relationship between the inherent structure landscape, the thermodynamics, and the dynamics of the fluid. The transfer matrix method is used to obtain an exact description of the landscape, equation of state, and provide a mapping of configurations of the equilibrium fluid to their local jammed structures. This allows us to follow how the system samples the landscape as a function of occupied volume fraction ϕ. Configurations of the ideal gas map to the maximum in the distribution of inherent structures, with a jamming volume fraction ϕ(J)(*), and sample more dense basins with increasing ϕ. This suggests jammed states with a density below ϕ(J)(*) are inaccessible from the equilibrium fluid. The configurational entropy of the fluid decreases rapidly at intermediate ϕ before plateauing at a low value and going to zero as the most dense packing is approached. This leads to the appearance of a maximum in both the isobaric heat capacity and the inherent structure pressure. We also show that the system exhibits a crossover from fragile to strong fluid behavior, located at the heat capacity maximum. Structural relaxation in the fragile fluid is shown to be controlled by the presence of high order saddle points caused by neighboring defects that are unstable with respect to jamming and spontaneously rearrange to form a stable local environment. In the strong fluid, the defect concentration is low so that defects do not interact and relaxation occurs through the hopping of isolated defects between stable local packing environments.

Static and dynamical properties of a hard-disk fluid confined to a narrow channel

The thermodynamic properties of disks moving in a channel sufficiently narrow that they can collide only with their nearest neighbors can be solved exactly by determining the eigenvalues and eigenfunctions of an integral equation. Using it, we have determined the correlation length ξ of this system. We have developed an approximate solution which becomes exact in the high-density limit. It describes the system in terms of defects in the regular zigzag arrangement of disks found in the high-density limit. The correlation length is then effectively the spacing between the defects. The time scales for defect creation and annihilation are determined with the help of transition-state theory, as is the diffusion coefficient of the defects, and these results are found to be in good agreement with molecular dynamics simulations. On compressing the system with the Lubachevsky-Stillinger procedure, jammed states are obtained whose packing fractions ϕJ are a function of the compression rate γ. We find a quantitative explanation of this dependence by making use of the Kibble-Zurek hypothesis. We have also determined the point-to-set length scale ξPS for this system. At a packing fraction ϕ close to its largest value ϕmax, ξPS has a simple power law divergence, ξPS∼1/(1-ϕ/ϕmax), while ξ diverges much faster, ln(ξ)∼1/(1-ϕ/ϕmax).

Detailed characterization of rattlers in exactly isostatic, strictly jammed sphere packings

We generate jammed disordered packings of 100≤N≤2000 monodisperse hard spheres in three dimensions whose strictly jammed backbones are demonstrated to be exactly isostatic with unprecedented numerical accuracy. This is accomplished by using the Torquato-Jiao (TJ) packing algorithm as a means of studying the maximally random jammed (MRJ) state. The rattler fraction of these packings converges towards 0.015 in the infinite-system limit, which is markedly lower than previous estimates for the MRJ state using the Lubachevsky-Stillinger protocol. This is because the packings that the TJ algorithm creates are closer to the true MRJ state, as shown using bond-orientational and translational order metrics. The rattler pair correlation statistics exhibit strongly correlated behavior contrary to the conventional understanding that they be randomly (Poisson) distributed. Dynamically interacting "polyrattlers" may be found imprisoned in shared cages as well as interacting through "bottlenecks" in the backbone and these clusters are mainly responsible for the sharp increase in the rattler pair correlation function near contact. We discover the surprising existence of polyrattlers with cluster sizes of up to five rattlers (which is expected to increase with system size) and present a distribution of polyrattler occurrence as a function of cluster size and system size. We also enumerate all of the rattler interaction topologies we observe and present images of several examples, showing that MRJ packings of monodisperse spheres can contain large rattler cages while still obeying the strict jamming criterion. The backbone spheres that encage the rattlers are significantly hypostatic, implying that correspondingly hyperstatic regions must exist elsewhere in these isostatic packings. We also observe that rattlers in hard-sphere packings share an apparent connection with the low-temperature two-level system anomalies that appear in real amorphous insulators and semiconductors.

Transition state theory and the dynamics of hard disks

The dynamics of two- and five-disk systems confined in a square has been studied using molecular dynamics simulations and compared with the predictions of transition state theory. We determine the partition functions Z and Z(‡) of transition state theory using a procedure first used by Salsburg and Wood for the pressure. Our simulations show this procedure and transition state theory are in excellent agreement with the simulations. A generalization of the transition state theory to the case of a large number of disks N is made and shown to be in full agreement with simulations of disks moving in a narrow channel. The same procedure for hard spheres in three dimensions leads to the Vogel-Fulcher-Tammann formula for their alpha relaxation time.