Simple effective rule to estimate the jamming packing fraction of polydisperse hard spheres

Searching for the maximal packing fraction of hard disks confined by a circular cavity through replica exchange/event-chain Monte Carlo

We utilized a blend of replica exchange and event-chain Monte Carlo techniques to generate candidate configurations, aiming for a maximal packing fraction of hard disks within a circular enclosure. Our investigation encompassed systems comprising N particles, with N ranging from 300 to 720. Through our analysis, we identified 108 novel maximal packings, with some surpassing existing configurations by over 0.001 in packing fraction. As such, Monte Carlo methods demonstrate their efficacy in tackling optimization challenges of this nature.

Rheology of debris flow materials is controlled by the distance from jamming

Debris flows are fast-flowing and dangerous slurries of soil and water that often form when intense rainfall soaks hillsides burned by wildfire. As climate change intensifies this hazard, models capable of predicting failure and flow behaviors are needed. Here we capitalize on recent progress in the physics of dense suspensions, to determine how the physical and chemical composition of natural hillslope soils controls the viscosity and yield stress of debris flows. We show how a simple flow model—previously developed for idealized suspensions—can be extended to highly heterogeneous, natural debris flow materials. This model reconciles previously contradictory observations and could help to improve computer models that assess the hazard potential of debris flows in the field.

Debris flows are dense and fast-moving complex suspensions of soil and water that threaten lives and infrastructure. Assessing the hazard potential of debris flows requires predicting yield and flow behavior. Reported measurements of rheology for debris flow slurries are highly variable and sometimes contradictory due to heterogeneity in particle composition and volume fraction (ϕ) and also inconsistent measurement methods. Here we examine the composition and flow behavior of source materials that formed the postwildfire debris flows in Montecito, CA, in 2018, for a wide range of ϕ that encapsulates debris flow formation by overland flow. We find that shear viscosity and yield stress are controlled by the distance from jamming, Δϕ=ϕm−ϕ, where the jamming fraction ϕm is a material parameter that depends on grain size polydispersity and friction. By rescaling shear and viscous stresses to account for these effects, the data collapse onto a simple nondimensional flow curve indicative of a Bingham plastic (viscoplastic) fluid. Given the highly nonlinear dependence of rheology on Δϕ, our findings suggest that determining the jamming fraction for natural materials will significantly improve flow models for geophysical suspensions such as hyperconcentrated flows and debris flows.

Finding the grain size distribution that produces the densest arrangement in frictional sphere packings: Revisiting and rediscovering the century-old Fuller and Thompson distribution

By means of discrete-element methods, we investigate the joint effects of the grain size distribution (GSD) and contact friction on the structure of three-dimensional samples composed of spherical grains. Specifically, we compress these systems isotropically until jamming and then analyze their structure in terms of density, connectivity, coefficients of uniformity and curvature, and parameters of grading entropy. Our study focuses on power-law GSDs and particularly on the Fuller and Thompson distribution, proposed over a century ago. First, we show that, among the set of GSDs investigated, this particular distribution produces the densest and best-connected systems, falsifying a conjecture recently posed in the literature. Second, we find that the jamming packing fraction can be accurately predicted as a function of simple descriptors of the GSD, but among these descriptors the granular entropy concept proves to be the most useful. This allows for an alternative interpretation of both jamming and grading entropy concepts. Finally, we compare the Fuller and Thompson distribution with two well-known GSDs: that of the Apollonian sphere packing and that towards which granular systems evolve after intensive grain fracturing. Surprisingly, we find that these three GSDs are practically coincident in the limit of large size spans, despite having been introduced or discovered in different scientific contexts (i.e., engineering, mathematics, and earth sciences, respectively).

Optimal packing in 2D and 3D granular systems: Density, connectivity, and force distributions

In this work, we explore the influence of the grain size distribution (GSD) on density, connectivity and internal forces distributions, for both 2D and 3D granular packings built mechanically. For power law GSDs, we show that there is an exponent for which density and connectivity are optimized, and this exponent is close to those that characterize other well known GSDs such as the Fuller and Thompson distribution and the Appollonian packing. In addition, we studied the distributions of normal forces, finding that these can be well described by a power-law tail, specially for the GSDs with large size span. These results highlight the role of the GSD on internal structure and suggest important consequences on macroscopic properties.

Structural and thermodynamic properties of hard-sphere fluids

This Perspective article provides an overview of some of our analytical approaches to the computation of the structural and thermodynamic properties of single-component and multicomponent hard-sphere fluids. For the structural properties, they yield a thermodynamically consistent formulation, thus improving and extending the known analytical results of the Percus-Yevick theory. Approximate expressions linking the equation of state of the single-component fluid to the one of the multicomponent mixtures are also discussed.

Yuste S. Equation of State of Four- and Five-Dimensional Hard-Hypersphere Mixtures

New proposals for the equation of state of four- and five-dimensional hard-hypersphere mixtures in terms of the equation of state of the corresponding monocomponent hard-hypersphere fluid are introduced. Such proposals (which are constructed in such a way so as to yield the exact third virial coefficient) extend, on the one hand, recent similar formulations for hard-disk and (three-dimensional) hard-sphere mixtures and, on the other hand, two of our previous proposals also linking the mixture equation of state and the one of the monocomponent fluid but unable to reproduce the exact third virial coefficient. The old and new proposals are tested by comparison with published molecular dynamics and Monte Carlo simulation results and their relative merit is evaluated.

i-Rheo: determining the linear viscoelastic moduli of colloidal dispersions from step-stress measurements

We report on the application of a Fourier transform-based method, 'i-Rheo', to evaluate the linear viscoelastic moduli of hard-sphere colloidal dispersions, both in the fluid and glass states, from a direct analysis of raw step-stress (creep) experimental data. We corroborate the efficacy of i-Rheo by comparing the outputs of creep tests performed on homogenous complex fluids to conventional dynamic frequency sweeps. A similar approach is adopted for a number of colloidal suspensions over a broad range of volume fractions. For these systems, we test the limits of the method by varying the applied stress across the materials' linear and non-linear viscoelastic regimes, and we show that the best results are achieved for stress values close to the upper limit of the materials' linear viscoelastic regime, where the signal-to-noise ratio is at its highest and the non-linear phenomena have not appeared yet. We record that, the range of accessible frequencies is controlled at the higher end by the relative weight between the inertia of the instrument and the elasticity of the complex material under investigation; whereas, the lowest accessible frequency is dictated by the extent of the materials' linear viscoelastic regime. Nonetheless, despite these constrains, we confirm the effectiveness of i-Rheo for gaining valuable information on the materials' linear viscoelastic properties even from 'creep ringing' data, confirming its potency and general validity as an accurate method for determining the material's rheological behaviour for a variety of complex systems.

Massive replica exchange Monte Carlo algorithm: a tool to access high pressure thermodynamics of hard systems

We have explored the idea of producing the equilibrium equation of state, i.e. the pressure as a function of packing fraction, βP(φ), of a confined system up to very high pressures to yield the configuration that leads to the maximum packing fraction. For this purpose we have massively implemented the replica exchange Monte Carlo algorithm in graphics processing units (GPUs), in such a way that each GPU core handles a single simulation cell. This yields a very easy scheme to implement parallelization for a very large amount of replicas (thousands), which densely sample configuration space. We have tested this idea with a very well studied system, i.e. discs confined in a circular cavity, for a number of particles N ≤ 125. In all cases, our outcomes for configurations having maximum packing fractions are in perfect agreement with those already reported and conjectured optimal in the literature, for which there is no formal mathematical proof, strongly suggesting that they are indeed optimal configurations. Furthermore, in most cases, we have obtained the same function βP(φ), by compressing loose random configurations and by decompressing copies of the configuration having the largest packing fraction. This reveals numerically that the so obtained maximum packing configurations are the correct answer.

Fabrication of optimized skin biomimics for improved interfacial retention of cosmetic emulsions

Retention of hydrophobic active agents on human skin following the use of skin-care formulations is an important indication of the performance of the deposited product. We have developed a novel system which replicates the interaction between human skin and a cosmetic emulsion to systematically establish and characterize the key parameters driving the retention process at the interface. This included a comprehensive study of the skin's biology and physical properties which influenced the process, the fabrication of advanced, improved skin biomimics, the formulation of a cosmetic model-system emulsion, comprising a hydrophobic active agent i.e. petrolatum, commonly used in cosmetic products, the development of a dedicated and highly consistent deposition rig with a corresponding cleaning set-up and the systematic characterization of retention processes on the developed mimics. This study further explores the interplay of petrolatum with skin biomimics and studies the mechanisms that give rise to improved interfacial retention. Petrolatum has been found to create an occlusive layer on the skin mimic, displaying high coverage from emulsion formulations. The large particle size emulsions yielded improved retention on the developed skin biomimics due to the microstructure of the emulsion and the counter effect of the surfactant.

Chemical potential of a test hard sphere of variable size in hard-sphere fluid mixtures

A detailed comparison between the Boublík-Mansoori-Carnahan-Starling-Leland (BMCSL) equation of state of hard-sphere mixtures is made with Molecular Dynamics (MD) simulations of the same compositions. The Labík and Smith simulation technique [S. Labík and W. R. Smith, Mol. Simul. 12, 23-31 (1994)] was used to implement the Widom particle insertion method to calculate the excess chemical potential, βμ₀^ex, of a test particle of variable diameter, σ₀, immersed in a hard-sphere fluid mixture with different compositions and values of the packing fraction, η. Use is made of the fact that the only polynomial representation of βμ₀^ex which is consistent with the limits σ₀ → 0 and σ₀ → ∞ has to be of the cubic form, i.e., c₀(η)+c¯₁(η)σ₀/M₁+c¯₂(η)(σ₀/M₁)²+c¯₃(η)(σ₀/M₁)³, where M₁ is the first moment of the distribution. The first two coefficients, c₀(η) and c¯₁(η), are known analytically, while c¯₂(η) and c¯₃(η) were obtained by fitting the MD data to this expression. This in turn provides a method to determine the excess free energy per particle, βa^ex, in terms of c¯₂, c¯₃, and the compressibility factor, Z. Very good agreement between the BMCSL formulas and the MD data is found for βμ₀^ex, Z, and βa^ex for binary mixtures and continuous particle size distributions with the top-hat analytic form. However, the BMCSL theory typically slightly underestimates the simulation values, especially for Z, differences which the Boublík-Carnahan-Starling-Kolafa formulas and an interpolation between two Percus-Yevick routes capture well in different ranges of the system parameter space.

Equation of state of polydisperse hard-disk mixtures in the high-density regime

A proposal to link the equation of state of a monocomponent hard-disk fluid to the equation of state of a polydisperse hard-disk mixture is presented. Event-driven molecular dynamics simulations are performed to obtain data for the compressibility factor of the monocomponent fluid and of 26 polydisperse mixtures with different size distributions. Those data are used to assess the proposal and to infer the values of the compressibility factor of the monocomponent hard-disk fluid in the metastable region from those of mixtures in the high-density region. The collapse of the curves for the different mixtures is excellent in the stable region. In the metastable regime, except for two mixtures in which crystallization is present, the outcome of the approach exhibits a rather good performance. The simulation results indicate that a (reduced) variance of the size distribution larger than about 0.01 is sufficient to avoid crystallization and explore the metastable fluid branch.

Subjamming transition in binary sphere mixtures

We study the influence of particle-size asymmetry on structural evolution of randomly jammed binary sphere mixtures with varying large-sphere and small-sphere composition. Simulations of jammed packings are used to assess the transition from large-sphere dominant to small-sphere dominant mixtures. For weakly asymmetric particle sizes, packing properties evolve smoothly, but not monotonically, with increasing small-sphere composition, f. Our simulations reveal that at high values of ratio α of large- to small-sphere radii (α≥α_{c}≈5.75), evolution of structural properties, such as packing density, fraction of jammed spheres, and contact statistics with f, exhibit features that suggest a sharp transition, either through discontinuities in structural measures or their derivatives. We argue that this behavior is related to the singular, composition dependence of close-packing fraction predicted in infinite aspect ratio mixtures α→∞ by the Furnas model, but occurring for finite valued range of α above a critical value, α_{c}≈5.75. The existence of a sharp transition from small- to large-f values for α≥α_{c} can be attributed to the existence of a subjamming transition of small spheres within the interstices of jammed large spheres along the line of compositions f_{sub}(α). We argue that the critical value of finite-size asymmetry α_{c}≃5.75 is consistent with the geometric criterion for the transmission of small-sphere contacts between neighboring tetrahedrally close-packed interstices of large spheres, facilitating a cooperative subjamming transition of small spheres confined within the disjoint volumes.

Size-Dependent Localization in Polydisperse Colloidal Glasses

We have investigated concentrated suspensions of polydisperse hard spheres and have determined the dynamics and sizes of individual particles using confocal microscopy. With increasing concentration, the dynamics of the small and large particles start to differ. The large particles exhibit slower dynamics and stronger localization. Moreover, as the particle size increases, the local volume fraction ϕ_{loc} also increases. In the glass state, the localization length significantly decreases beyond ϕ_{loc}≈0.67. This suggests a link between local crowding and dynamical heterogeneities. However dynamical arrest of subpopulations seems not directly linked to a large value of ϕ_{loc}, indicating the importance of collective effects.

A precise algorithm to detect voids in polydisperse circle packings

Computer simulations are the primary tool for studying polydisperse particle packings quanti- tatively. For the problem of packing N unequal circles in a larger container circle, nothing is known a priori about the optimal packing (i.e. the packing with the highest packing fraction). Simulations usually start from a random initial configuration with the aim to finish with a dense final packing. Unfortunately, smaller circles often get stuck in trapped positions and prevent the rest of the packing from growing larger. Hence, the knowledge of the structure of unoccupied areas or holes inside a packing is important to be able to move trapped circles into free circular places or voids . A novel algorithm is proposed for detecting such voids in two-dimensional arbitrary circle packings by a decomposition of the contact graph. Combined with a clever object jumping strategy and together with other heuristic methods like swaps and shifts, this approach increases the packing fraction ϕ significantly. Its effectiveness for jumping across the maximally random jammed barrier ( ϕ MRJ ≈0.8575 in the large- N limit) for small benchmark instances as well as for large problem sizes (up to N ≈10 3 ) is demonstrated.

Extension of the BMCSL equation of state for hard spheres to the metastable disordered region: Application to the SAFT approach

A simple modification of the Boublík-Mansoori-Carnahan-Starling-Leland equation of state is proposed for an application to the metastable disordered region. The new model has a positive pole at the jamming limit and can accurately describe the molecular simulation data of pure hard in the stable fluid region and along the metastable branch. The new model has also been applied to binary mixtures hard spheres, and an excellent description of the fluid and metastable branches can be obtained by adjusting the jamming packing fraction. The new model for hard sphere mixtures can be used as the repulsive term of equations of state for real fluids. In this case, the modified equations of state give very similar predictions of thermodynamic properties as the original models, and one can remove the multiple liquid density roots observed for some versions of the Statistical Associating Fluid Theory (SAFT) at low temperature without any modification of the dispersion term.

Influence of particle size distribution on random close packing of spheres

The densest amorphous packing of rigid particles is known as random close packing. It has long been appreciated that higher densities are achieved by using collections of particles with a variety of sizes. For spheres, the variety of sizes is often quantified by the polydispersity of the particle size distribution: the standard deviation of the radius divided by the mean radius. Several prior studies quantified the increase of the packing density as a function of polydispersity. A particle size distribution is also characterized by its skewness, kurtosis, and higher moments, but the influence of these parameters has not been carefully quantified before. In this work, we numerically generate many sphere packings with different particle radii distributions, varying polydispersity and skewness independently of one another. We find that the packing density can increase significantly with increasing skewness and in some cases skewness can have a larger effect than polydispersity. However, the packing fraction is relatively insensitive to the higher moment value of the kurtosis. We present a simple empirical formula for the value of the random close packing density as a function of polydispersity and skewness.