Glass-transition properties of Yukawa potentials: from charged point particles to hard spheres

Shoving model and the glass transition in one-component plasma

A modified shoving model is applied to estimate the location of the glass transition in a one-component plasma. The estimated value of the coupling parameter Γ ≃ 570 at the glass transition is compared with other predictions available in the literature.

Bridge functions of classical one-component plasmas

In a recent paper, Lucco Castello et al. [arXiv:2107.03537] performed systematic extractions of classical one-component plasma bridge functions from molecular dynamics simulations and provided an accurate parametrization that was incorporated in their isomorph-based empirically modified hypernetted chain approach for Yukawa one-component plasmas. Here the extraction technique and parametrization strategy are described in detail, while the deficiencies of earlier efforts are discussed. The structural and thermodynamic predictions of the updated version of the integral equation theory approach are compared with extensive available simulation results revealing a truly unprecedented level of accuracy in the entire dense liquid region of the Yukawa phase diagram.

Structure and thermodynamics of two-dimensional Yukawa liquids

The thermodynamic and structural properties of two-dimensional dense Yukawa liquids are studied with molecular dynamics simulations. The "exact" thermodynamic properties are simultaneously employed in an advanced scheme for the determination of an equation of state that shows an unprecedented level of accuracy for the internal energy, pressure, and isothermal compressibility. The "exact" structural properties are utilized to formulate a novel empirical correction to the hypernetted-chain approach that leads to a very high accuracy level in terms of static correlations and thermodynamics.

Theoretical Estimate of the Glass Transition Line of Yukawa One-Component Plasmas

The mode coupling theory of supercooled liquids is combined with advanced closures to the integral equation theory of liquids in order to estimate the glass transition line of Yukawa one-component plasmas from the unscreened Coulomb limit up to the strong screening regime. The present predictions constitute a major improvement over the current literature predictions. The calculations confirm the validity of an existing analytical parameterization of the glass transition line. It is verified that the glass transition line is an approximate isomorphic curve and the value of the corresponding reduced excess entropy is estimated. Capitalizing on the isomorphic nature of the glass transition line, two structural vitrification indicators are identified that allow a rough estimate of the glass transition point only through simple curve metrics of the static properties of supercooled liquids. The vitrification indicators are demonstrated to be quasi-universal by an investigation of hard sphere and inverse power law supercooled liquids. The straightforward extension of the present results to bi-Yukawa systems is also discussed.

Formation of a transient amorphous solid in low density aqueous charged sphere suspensions

Colloidal glasses formed from hard spheres, nearly hard spheres, ellipsoids and platelets or their attractive variants, have been studied in great detail. Complementing and constraining theoretical approaches and simulations, the many different types of model systems have significantly advanced our understanding of the glass transition in general. Despite their early prediction, however, no experimental charged sphere glasses have been found at low density, where the competing process of crystallization prevails. We here report the formation of a transient amorphous solid formed from charged polymer spheres suspended in thoroughly deionized water at volume fractions of 0.0002-0.01. From optical experiments, we observe the presence of short-range order and an enhanced shear rigidity as compared to the stable polycrystalline solid of body centred cubic structure. On a density dependent time scale of hours to days, the amorphous solid transforms into this stable structure. We further present preliminary dynamic light scattering data showing the evolution of a second slow relaxation process possibly pointing to a dynamic heterogeneity known from other colloidal glasses and gels. We compare our findings to the predicted phase behaviour of charged sphere suspensions and discuss possible mechanisms for the formation of this peculiar type of colloidal glass.

Equilibrium fluid-crystal interfacial free energy of bcc-crystallizing aqueous suspensions of polydisperse charged spheres

The interfacial free energy is a central quantity in crystallization from the metastable melt. In suspensions of charged colloidal spheres, nucleation and growth kinetics can be accurately measured from optical experiments. In previous work, from these data effective nonequilibrium values for the interfacial free energy between the emerging bcc nuclei and the adjacent melt in dependence on the chemical potential difference between melt phase and crystal phase were derived using classical nucleation theory (CNT). A strictly linear increase of the interfacial free energy was observed as a function of increased metastability. Here, we further analyze these data for five aqueous suspensions of charged spheres and one binary mixture. We utilize a simple extrapolation scheme and interpret our findings in view of Turnbull's empirical rule. This enables us to present the first systematic experimental estimates for a reduced interfacial free energy, σ(0,bcc), between the bcc-crystal phase and the coexisting equilibrium fluid. Values obtained for σ(0,bcc) are on the order of a few k(B)T. Their values are not correlated to any of the electrostatic interaction parameters but rather show a systematic decrease with increasing size polydispersity and a lower value for the mixture as compared to the pure components. At the same time, σ(0) also shows an approximately linear correlation to the entropy of freezing. The equilibrium interfacial free energy of strictly monodisperse charged spheres may therefore be still greater.

Glass transition of charged particles in two-dimensional confinement

The glass transition of mesoscopic charged particles in two-dimensional confinement is studied by mode-coupling theory. We consider two types of effective interactions between the particles, corresponding to two different models for the distribution of surrounding ions that are integrated out in coarse-grained descriptions. In the first model, a planar monolayer of charged particles is immersed in an unbounded isotropic bath of ions, giving rise to an isotropically screened Debye-Hückel (Yukawa)-type effective interaction. The second, experimentally more relevant system is a monolayer of negatively charged particles that levitate atop a flat horizontal electrode, as frequently encountered in laboratory experiments with complex (dusty) plasmas. A steady plasma current toward the electrode gives rise to an anisotropic effective interaction potential between the particles, with an algebraically long-ranged in-plane decay. In a comprehensive parameter scan that covers the typical range of experimentally accessible plasma conditions, we calculate and compare the mode-coupling predictions for the glass transition in both kinds of systems.

Soft mean spherical approximation for dusty plasma liquids: One-component Yukawa systems with plasma shielding

The structure and thermodynamics of strongly coupled dusty plasmas are investigated with the soft mean spherical approximation. This integral theory approach is analytically solvable for Yukawa pair interactions yielding a closed-form solution for the direct correlation function. The pair correlation function, the structure factor, and basic thermodynamic quantities are calculated for a wide range of parameters. Exact consistency between the "energy"-"virial" thermodynamic routes and approximate consistency between the "energy"-"compressibility" paths is demonstrated. Comparison with extensive molecular dynamics results is carried out and a remarkable agreement from the Coulomb limit to the strongly screened limit is revealed. The soft mean spherical approximation is concluded to be particularly well suited for the study of dusty plasma liquids, uniquely combining simplicity and accuracy.