Phase diagram and dynamics of Yukawa systems

Coupled dynamics in binary mixtures of model colloidal Yukawa systems

The dynamical behavior of binary mixtures consisting of highly charged colloidal particles is studied by means of Brownian dynamics simulations. We investigate differently sized, but identically charged particles with nearly identical interactions between all species in highly dilute suspensions. Different short-time self-diffusion coefficients induce, mediated by electrostatic interactions, a coupling of both self and collective dynamics of differently sized particles: the long-time self-diffusion coefficients of a larger species are increased by the presence of a more mobile, smaller species and vice versa. Similar coupling effects are observed in collective dynamics where both time constant and functional form of intermediate scattering functions' initial decay are influenced by the presence of a differently sized species. We provide a systematic analysis of coupling effects in dependence on the ratio of sizes, number densities, and the strength of electrostatic interactions.

Accessing the free expansion of a crystalline colloidal drop by optical experiments

We study poly-crystalline spherical drops of an aqueous suspension of highly charged colloidal spheres exposed to a colloid-free aqueous environment. Crystal contours were obtained from standard optical imaging. The crystal spheres first expand to nearly four times their initial volume before slowly shrinking due to dilution-induced melting. Exploiting coherent multiple-scattering by (110) Bragg reflecting crystals, time-dependent density profiles were recorded within the drop interior. These show a continuously flattening radial density gradient and a decreasing central density. Expansion curves and density profiles are qualitatively consistent with theoretical expectations based on dynamical density functional theory for the expansion of a spherical crystallite made of charged Brownian spheres. We anticipate that our study opens novel experimental access to density determination in turbid crystals.

Modeling order-disorder boundaries of colloidal dispersions in organic solvents using interaction force measurements

HYPOTHESIS

The formation of soft colloidal crystals, which are nonclose-packed ordered arrays of colloidal particles suspended in a solvent, is dictated by a single physical factor that yields a fixed threshold at order-disorder boundaries for different experimental conditions such as ion concentration, solvent type, and particle size. Identifying the determinant factor and its threshold value should enable the prediction of the critical concentrations of colloidal particles to form soft colloidal crystals.

EXPERIMENTS

Soft colloidal crystals were fabricated using a series of monohydric alcohols as dispersion media and reflectance spectra were measured to locate order-disorder boundaries. The interaction forces acting between particles were also measured by employing atomic force microscopy.

FINDINGS

The interparticle forces at the order-disorder boundaries exhibited a universal threshold that was independent of the solvent types including alcohols and water. Therefore, the determinant factor for the formation of soft colloidal crystals was determined to be the force acting between the particles. Furthermore, a priori calculation of this critical force and consequently the critical particle concentration in colloidal systems was demonstrated by referring to the pressure at the liquid-to-solid transition in a hard sphere system (Alder transition).

Transformations in crystals of DNA-functionalized nanoparticles by electrolytes

Colloidal crystals have applications in water treatments, including water purification and desalination technologies. It is, therefore, important to understand the interactions between colloids as a function of electrolyte concentration. We study the assembly of DNA-grafted gold nanoparticles immersed in concentrated electrolyte solutions. Increasing the concentration of divalent Ca2+ ions leads to the condensation of nanoparticles into face-centered-cubic (FCC) crystals at low electrolyte concentrations. As the electrolyte concentration increases, the system undergoes a phase change to body-centered-cubic (BCC) crystals. This phase change occurs as the interparticle distance decreases. Molecular dynamics analysis suggests that the interparticle interactions change from strongly repulsive to short-range attractive as the divalent-electrolyte concentration increases. A thermodynamic analysis suggests that increasing the salt concentration leads to significant dehydration of the nanoparticle environment. We conjecture that the intercolloid attractive interactions and dehydrated states favour the BCC structure. Our results gain insight into salting out of colloids such as proteins as the concentration of salt increases in the solution.

Porous crystals in charged sphere suspensions by aggregate-driven phase separation

The kinetics of phase transition processes often governs the resulting material microstructure. Using optical microscopy, we here investigate the formation and stabilization of a porous crystalline microstructure forming in low-salt suspensions of charged colloidal spheres containing aggregates comprising some 5-10 of these colloids. We observe the transformation of an initially crystalline colloidal solid with homogeneously incorporated aggregates to individual, compositionally refined crystallites of perforated morphology coexisting with an aggregate-enriched fluid phase filling the holes and separating individual crystallites. A preliminary kinetic characterization suggests that the involved processes follow power laws. We show that this route to porous materials is neither restricted to nominally single component systems nor to a particular microstructure to start from. However, it necessitates an early rapid solidification stage during which the aggregates become trapped in the bulk of the host-crystals. The thermodynamic stability of the reconstructed crystalline scaffold against melting under increased salinity was found comparable to that of pure phase crystallites grown very slowly from a melt. Future implications of this novel route to porous colloidal crystals are discussed.

Spherical Silver Nanocrystals Arranged in a Metastable Square Pattern

The present article demonstrates the development of two-dimensional (2D) assembly of spherical nanocrystals (NCs) in the square arrangement through the delicate balance between repulsive ligand interactions and attractive van der Waals interactions of NCs, respectively, instead of the otherwise stable hexagonal arrangement. The experimental packing efficiency values matched quite well with the theoretically calculated square arrangement patterns. The above fact indicates that the formation of the 2D square arrangement of silver NCs can be explained by introducing the concept of softness to NCs in the hard sphere model.

Fingerprinting Brownian Motions of Polymers under Flow

We present a quantitative approach to the self-dynamics of polymers under steady flow by employing a set of complementary reference frames and extending the spherical harmonic expansion technique to dynamic density correlations. Application of this method to nonequilibrium molecular dynamics simulations of polymer melts reveals a number of universal features. For both unentangled and entangled melts, the center-of-mass motions in the flow frame are described by superdiffusive, anisotropic Gaussian distributions, whereas the isotropic component of monomer self-dynamics in the center-of-mass frame is strongly suppressed. Spatial correlation analysis shows that the heterogeneity of monomer self-dynamics increases significantly under flow.

Computer simulations of self-assembly of anisotropic colloids

Computer simulations have played a significant role in understanding the physics of colloidal self-assembly, interpreting experimental observations, and predicting novel mesoscopic and crystalline structures. Recent advances in computer simulations of colloidal self-assembly driven by anisotropic or orientation-dependent inter-particle interactions are highlighted in this review. These interactions are broadly classified into two classes: entropic and enthalpic interactions. They mainly arise due to shape anisotropy, surface heterogeneity, compositional heterogeneity, external field, interfaces, and confinements. Key challenges and opportunities in the field are discussed.

Salt-concentration-dependent nucleation rates in low-metastability colloidal charged sphere melts containing small amounts of doublets

We determined bulk crystal nucleation rates in aqueous suspensions of charged spheres at low metastability. Experiments were performed in dependence on electrolyte concentration and for two different particle number densities. The time-dependent nucleation rate shows a pronounced initial peak, while postsolidification crystal size distributions are skewed towards larger crystallite sizes. At each concentration, the nucleation rate density initially drops exponentially with increasing salt concentration. The full data set, however, shows an unexpected scaling of the nucleation rate densities with metastability times the number density of particles. Parameterization of our results in terms of classical nucleation theory reveals unusually low interfacial free energies of the nucleus surfaces and nucleation barriers well below the thermal energy. We tentatively attribute our observations to the presence of doublets introduced by the employed conditioning technique. We discuss the conditions under which such small seeds may induce nucleation.

Excess entropy and Stokes-Einstein relation in simple fluids

The Stokes-Einstein (SE) relation between the self-diffusion and shear viscosity coefficients operates in sufficiently dense liquids not too far from the liquid-solid phase transition. By considering four simple model systems with very different pairwise interaction potentials (Lennard-Jones, Coulomb, Debye-Hückel or screened Coulomb, and the hard sphere limit) we identify where exactly on the respective phase diagrams the SE relation holds. It appears that the reduced excess entropy s_{ex} can be used as a suitable indicator of the validity of the SE relation. In all cases considered the onset of SE relation validity occurs at approximately s_{ex}≲-2. In addition, we demonstrate that the line separating gaslike and liquidlike fluid behaviours on the phase diagram is roughly characterized by s_{ex}≃-1.

Foam film stratification studies probe intermicellar interactions

Ultrathin foam films containing supramolecular structures like micelles in bulk and adsorbed surfactant at the liquid-air interface undergo drainage via stratification. At a fixed surfactant concentration, the stepwise decrease in the average film thickness of a stratifying micellar film yields a characteristic step size that also describes the quantized thickness difference between coexisting thick-thin flat regions. Even though many published studies claim that step size equals intermicellar distance obtained using scattering from bulk solutions, we found no reports of a direct comparison between the two length scales. It is well established that step size is inversely proportional to the cubic root of surfactant concentration but cannot be estimated by adding micelle size to Debye length, as the latter is inversely proportional to the square root of surfactant concentration. In this contribution, we contrast the step size obtained from analysis of nanoscopic thickness variations and transitions in stratifying foam films using Interferometry Digital Imaging Optical Microscopy (IDIOM) protocols, that we developed, with the intermicellar distance obtained using small-angle X-ray scattering. We find that stratification driven by the confinement-induced layering of micelles within the liquid-air interfaces of a foam film provides a sensitive probe of non-DLVO (Derjaguin-Landau-Verwey-Overbeek) supramolecular oscillatory structural forces and micellar interactions.

Crystallization of charged gold particles mediated by nonadsorbing like-charged polyelectrolyte

We report that the aqueous dispersions of negatively charged submicron-sized colloidal Au particles formed non-close-packed colloidal crystals by the addition of a like-charged linear polyelectrolyte, sodium polyacrylate (NaPAA). Au particles often form irregular aggregates in dispersions because of a strong van der Waals force acting between them. To prevent aggregation, we introduced negative electric charges on particle surfaces. By the addition of NaPAA, colloidal crystals were formed on the bottom of a sample cell because of the supply of Au particles by sedimentation and 2D diffusion even under very dilute conditions. Interparticle potential calculations demonstrated that the addition of NaPAA caused depletion attraction between the particles as well as a significant reduction in the interparticle repulsion because of the electrostatic screening effect. However, the electrostatic repulsion was strong enough to prevent the direct contact of particles in the excluded region between Au particles. Large-area crystals could be obtained by tilting the sample cell. By drying the sample, the Au particles came into contact and the non-space-filling crystals changed into closest packed crystals. These closest packed crystals exhibited a significant enhancement of Raman scattering intensity because of high hot-spot density.

Defects in crystals of soft colloidal particles

In this paper we use computer simulations to examine point defects in systems of "soft" colloidal particles including Hertzian spheres, and star polymers. We use Monte Carlo simulations to determine the deformation of the different crystals associated with vacancies and interstitials and use thermodynamic integration to predict the equilibrium concentrations of such defects. We find that the nature of the lattice distortion is mainly determined by the crystal structure and not by the specifics of the interaction potential. We can distinguish one-, two-, and three-dimensional lattice distortions and find that the range of the distortion generally depends on the dimensionality. We find that in both model systems the deformation of the body-centered cubic (BCC) crystal caused by an interstitial is one-dimensional and we show that its structure is well described as a crowdion. Similarly, we show that the one-dimensional deformation of the hexagonal (H) crystal of Hertzian spheres caused by a vacancy can be characterized as a voidion. Interestingly, with the exception of the FCC crystal in the Hertzian sphere model, in all cases we find that the interstitial concentration is higher than the vacancy concentration. Most noteworthy, the concentration of interstitials in the BCC crystals can reach up to 1%.

Electrostatic energy of Coulomb crystals with polarized electron background

Outer crusts of neutron stars and interiors of cool white dwarfs consist of bare atomic nuclei, arranged in a crystal lattice and immersed in a Fermi gas of degenerate electrons. We study electrostatic properties of such Coulomb crystals, taking into account the polarizability of the electron gas and considering different structures, which can form the ground state: body-centered-cubic (bcc), face-centered-cubic (fcc), hexagonal close-packed (hcp), and MgB_{2}-like lattices. At zero temperature the electrostatic energy provides a fundamental contribution to the total energy of the classical Coulomb crystal, which allows us to study structural transitions in the neutron-star crusts and crystallized white-dwarf interiors. To take the electron background polarization into account, we use the linear response theory with the electron dielectric function given either by the Thomas-Fermi approximation or by the random-phase approximation (RPA). We compare the widely used nonrelativistic (Lindhard) version of the RPA with the more general, relativistic (Jancovici) version. The results of the different approximations are compared to assess the importance of going beyond the Thomas-Fermi or Lindhard approximations. We also include the contribution of zero-point vibrations of ions into the ground-state energy. We show that the bcc lattice forms the ground state for any charge number Z of the atomic nuclei at the densities where the electrons are relativistic (ρ≳10^{6}g cm^{-3}), while at the nonrelativistic densities (ρ≲10^{6}g cm^{-3}) the fcc and hcp lattices can form the ground state. The MgB_{2}-like lattice never forms the ground state at realistic densities in the crystallized regions of degenerate stars. The RPA corrections strongly affect the boundaries between the phases. As a result, transitions between different ground-state structures depend on Z in a nontrivial way. The relativistic and quantum corrections produce less dramatic effects, moderately shifting the phase boundaries.

Dynamical properties across different coarse-grained models for ionic liquids

Room-temperature ionic liquids (RTILs) stand out among molecular liquids for their rich physicochemical characteristics, including structural and dynamic heterogeneity. The significance of electrostatic interactions in RTILs results in long characteristic length- and timescales, and has motivated the development of a number of coarse-grained (CG) simulation models. In this study, we aim to better understand the connection between certain CG parameterization strategies and the dynamical properties and transferability of the resulting models. We systematically compare five CG models: a model largely parameterized from experimental thermodynamic observables; a refinement of this model to increase its structural accuracy; and three models that reproduce a given set of structural distribution functions by construction, with varying intramolecular parameterizations and reference temperatures. All five CG models display limited structural transferability over temperature, and also result in various effective dynamical speedup factors, relative to a reference atomistic model. On the other hand, the structure-based CG models tend to result in more consistent cation-anion relative diffusion than the thermodynamic-based models, for a single thermodynamic state point. By linking short- and long-timescale dynamical behaviors, we demonstrate that the varying dynamical properties of the different CG models can be largely collapsed onto a single curve, which provides evidence for a route to constructing dynamically-consistent CG models of RTILs.

Point defects in crystals of charged colloids

Charged colloidal particles-on both the nano and micron scales-have been instrumental in enhancing our understanding of both atomic and colloidal crystals. These systems can be straightforwardly realized in the lab and tuned to self-assemble into body-centered-cubic (BCC) and face-centered-cubic (FCC) crystals. While these crystals will always exhibit a finite number of point defects, including vacancies and interstitials-which can dramatically impact their material properties-their existence is usually ignored in scientific studies. Here, we use computer simulations and free-energy calculations to characterize vacancies and interstitials in FCC and BCC crystals of point-Yukawa particles. We show that, in the BCC phase, defects are surprisingly more common than in the FCC phase, and the interstitials manifest as so-called crowdions: an exotic one-dimensional defect proposed to exist in atomic BCC crystals. Our results open the door to directly observe these elusive defects in the lab.

Spontaneous organization of supracolloids into three-dimensional structured materials

Periodic nano- or microscale structures are used to control light, energy and mass transportation. Colloidal organization is the most versatile method used to control nano- and microscale order, and employs either the enthalpy-driven self-assembly of particles at a low concentration or the entropy-driven packing of particles at a high concentration. Nonetheless, it cannot yet provide the spontaneous three-dimensional organization of multicomponent particles at a high concentration. Here we combined these two concepts into a single strategy to achieve hierarchical multicomponent materials. We tuned the electrostatic attraction between polymer and silica nanoparticles to create dynamic supracolloids whose components, on drying, reorganize by entropy into three-dimensional structured materials. Cryogenic electron tomography reveals the kinetic pathways, whereas Monte Carlo simulations combined with a kinetic model provide design rules to form the supracolloids and control the kinetic pathways. This approach may be useful to fabricate hierarchical hybrid materials for distinct technological applications.

Stiffness of the interface between a colloidal body-centered cubic crystal and its liquid

Equilibrium interfaces were established between body-centered cubic (BCC) crystals and their liquid using charged colloidal particles in an electric bottle. By measuring a time series of interfacial positions and computing the average power spectrum, their interfacial stiffness was determined according to the capillary fluctuation method. For the (100) and the (114) interfaces, the stiffnesses were 0.15 and 0.18 [Formula: see text] (σ: particle diameter), respectively, and were isotropic in the plane of the interface. For comparison, similar charged colloids were used to create an interface between a face-centered cubic (FCC) crystal and its liquid. Its stiffness was significantly larger: 0.26 [Formula: see text] This result gives experimental support to the explanations offered for the preferential nucleation of BCC over FCC in metallic alloys.

Dispersion relations of Yukawa fluids at weak and moderate coupling

In this paper we compare different theoretical approaches to describe the dispersion of collective modes in Yukawa fluids when the interparticle coupling is relatively weak, so that the kinetic and potential contributions to the dispersion relation compete with each other. A thorough comparison with the results from molecular dynamics simulation allows us to conclude that, in the investigated regime, the best description is provided by the sum of the generalized excess bulk modulus and the Bohm-Gross kinetic term.

Colloidal crystallites under external oscillation

We study the two-dimensional assemblies of interacting colloidal particles in a loosely focussed optical trap. As the optical confinement is increased, the system becomes ordered and we investigate how these crystallites maintain their order under externally imposed oscillation. For small amplitudes, the crystalline order remains intact and the system behaves like a rigid body as predicted by numerical simulations. However, the rigidity breaks at large amplitudes, which we infer to be caused by the anharmonic component of the confinement potential. These studies are general enough to be applied to other physical systems comprising ordered finite-sized assemblies under external dynamic perturbation.

Influence of Softness on the Stability of Binary Colloidal Crystals

Mixtures of two types of nanoparticles can self-assemble into a wide variety of binary colloidal crystals (also called binary nanoparticle superlattices), which are interesting for their structural diversity and potential applications. Although so-called packing models-which usually treat the particles as "hard" with only excluded volume interactions-seem to explain many reported dense crystalline phases, these models often fail to predict the right structure. Here, we examine the role of soft repulsive interparticle interactions on binary colloidal crystals comprising two sizes of spherical particles; such "softness" can arise due to ligand shells or screened electrostatics. We determine the ground state phase diagram of binary systems of particles interacting with an additive inverse power law potential using a basin hopping algorithm to calculate the enthalpy of an extremely large pool of candidate structures. We find that a surprisingly small amount of softness can destabilize dense packings in favor of less densely packed structures, which provides further evidence that considerations beyond packing are necessary for describing many of the observed phases of binary colloidal crystals. Importantly, we find that several of the phases stabilized by softness, which are characterized by relatively few interparticle contacts and a tendency for local icosahedral order, are more likely to be observed experimentally than those predicted by packing models. We also report a previously unknown dense AB₄ phase and conduct free energy calculations to examine how the stability of several crystals will vary with temperature. Our results further our understanding of why particular binary colloidal crystals form and will be useful as a reference for experimentalists working with softly repulsive colloids.

The EXP pair-potential system. III. Thermodynamic phase diagram

This paper determines the thermodynamic phase diagram of the EXP system of particles interacting by the purely repulsive exponential pair potential. The solid phase is face-centered cubic (fcc) at low densities and pressures. At higher densities and pressures, the solid phase is body-centered cubic (bcc) with a re-entrant liquid phase at the highest pressures simulated. The investigation first identifies the phase diagram at zero temperature at which the following four crystal structures are considered: fcc, bcc, hexagonal close packed, and cubic diamond. There is a T = 0 phase transition at pressure 2.651 × 10^-3 with the thermodynamically stable structure being fcc below and bcc above this pressure. The densities of the two crystal structures at the phase transition are 1.7469 × 10^-2 (fcc) and 1.7471 × 10^-2 (bcc). At finite temperatures, the fcc-bcc, fcc-liquid, and bcc-liquid coexistence lines are determined by numerical integration of the Clausius-Clapeyron equation and validated by interface-pinning simulations at selected state points. The bcc-fcc phase transition is a weak first-order transition. The liquid-fcc-bcc triple point, which is determined by the interface-pinning method, has temperature 5.9 × 10^-5 and pressure 2.5 × 10^-6; the triple-point densities are 1.556 × 10^-3 (liquid), 1.583 × 10^-3 (bcc), and 1.587 × 10^-3 (fcc).

Hydrodynamic simulations of charge-regulation effects in colloidal suspensions

Self-organization of charged soft matter is of crucial importance in biology. However, it is an extremely complex phenomenon due to dynamical couplings between the hydrodynamic flow, ions, and charges of soft matter. For colloidal suspensions, the coupling between the former two has already been studied by numerical simulations, with the colloid surface charge being fixed. However, the self-organization of colloids and/or the application of an external electric field make the electrostatic environment of each colloid inhomogeneous in both space and time. Thus, this leads to inhomogenisation of the surface charge of each colloid under ionisation equilibrium conditions. This effect is known as "charge regulation" and is of great importance in various electrostatic and electrokinetic phenomena not only in colloid suspensions but also in solutions of biomolecules. However, there has so far been no success in taking the charge regulation effect into account in numerical simulations of colloidal electrokinetics. Here, we extend the fluid particle dynamics (FPD) method to incorporate the charge regulation effect. We present a theoretical formulation of the method and its application to two types of problems, where charge regulation plays an important role: (i) cluster formation of colloid particles and (ii) a single colloid particle under an external field. By these examples, we show not only the importance of considering charge-regulation effects in the self-organization of charged systems but also the applicability of our simulation method to more complex problems of charged soft matter systems.

Thermodynamics of two-dimensional Yukawa systems across coupling regimes

Thermodynamics of two-dimensional Yukawa (screened Coulomb or Debye-Hückel) systems is studied systematically using molecular dynamics (MD) simulations. Simulations cover very broad parameter range spanning from weakly coupled gaseous states to strongly coupled fluid and crystalline states. Important thermodynamic quantities, such as internal energy and pressure, are obtained and accurate physically motivated fits are proposed. This allows us to put forward simple practical expressions to describe thermodynamic properties of two-dimensional Yukawa systems. For crystals, in addition to numerical simulations, the recently developed shortest-graph interpolation method is applied to describe pair correlations and hence thermodynamic properties. It is shown that the finite-temperature effects can be accounted for by using simple correction of peaks in the pair correlation function. The corresponding correction coefficients are evaluated using MD simulation. The relevance of the obtained results in the context of colloidal systems, complex (dusty) plasmas, and ions absorbed to interfaces in electrolytes is pointed out.

Molecular dynamics simulation of charged colloids confined between hard walls: pre-melting and pre-freezing across the BCC-fluid coexistence

Molecular dynamics (MD) computer simulations are used to study the structure of hard-core Yukawa systems confined between two parallel hard walls. States around the coexistence between a fluid and a body-centered cubic (BCC) crystal are considered. In all cases a pronounced layering in the vicinity of the walls is observed. Using a thermodynamic integration scheme, we determine the wall-fluid interfacial free energy γ which is negative and monotonically decreasing with increasing bulk density of the fluid. In the case of the fluid, the layers next to the walls undergo a transition from a fluid to a hexagonal structure. This pre-freezing transition occurs well below the coexistence bulk density of the fluid. The confined BCC crystal in (111) orientation shows melted regions between crystalline face-centered cubic (FCC) layers close to the wall and the BCC bulk region.

Structure formation in soft nanocolloids: liquid-drop model

Using a model where soft nanocolloids such as spherical polymer brushes and star polymers are viewed as compressible liquid drops, we theoretically explore contact interactions between such particles. By numerically minimizing the phenomenological free energy consisting of bulk and surface terms, we find that at small deformations the drop-drop interaction is pairwise additive and described by a power law. We also propose a theory to describe the small-deformation regime, and the agreement is very good at all drop compressibilities. The large-deformation regime, which is dominated by many-body interactions, is marked by a rich phase diagram which includes the face- and body-centered-cubic, σ, A15, and simple hexagonal lattice as well as isostructural and re-entrant transitions. Most of these features are directly related to the non-convex deformation free energy emerging from many-body effects in the partial-faceting regime. The phase diagram, which depends on just two model parameters, contains many of the condensed phases observed in experiments. We also provide statistical-mechanical arguments that relate the two model parameters to the molecular architecture of the polymeric nanocolloids, chain rigidity, and solvent quality. The model represents a generic framework for the overarching features of the phase behavior of polymeric nanocolloids at high compressions.

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.

Depinning and nonequilibrium dynamic phases of particle assemblies driven over random and ordered substrates: a review

We review the depinning and nonequilibrium phases of collectively interacting particle systems driven over random or periodic substrates. This type of system is relevant to vortices in type-II superconductors, sliding charge density waves, electron crystals, colloids, stripe and pattern forming systems, and skyrmions, and could also have connections to jamming, glassy behaviors, and active matter. These systems are also ideal for exploring the broader issues of characterizing transient and steady state nonequilibrium flow phases as well as nonequilibrium phase transitions between distinct dynamical phases, analogous to phase transitions between different equilibrium states. We discuss the differences between elastic and plastic depinning on random substrates and the different types of nonequilibrium phases which are associated with specific features in the velocity-force curves, fluctuation spectra, scaling relations, and local or global particle ordering. We describe how these quantities can change depending on the dimension, anisotropy, disorder strength, and the presence of hysteresis. Within the moving phase we discuss how there can be a transition from a liquid-like state to dynamically ordered moving crystal, smectic, or nematic states. Systems with periodic or quasiperiodic substrates can have multiple nonequilibrium second or first order transitions in the moving state between chaotic and coherent phases, and can exhibit hysteresis. We also discuss systems with competing repulsive and attractive interactions, which undergo dynamical transitions into stripes and other complex morphologies when driven over random substrates. Throughout this work we highlight open issues and future directions such as absorbing phase transitions, nonequilibrium work relations, inertia, the role of non-dissipative dynamics such as Magnus effects, and how these results could be extended to the broader issues of plasticity in crystals, amorphous solids, and jamming phenomena.

Role of Entropy in the Expulsion of Dopants from Optically Trapped Colloidal Assemblies

Controlling an assembly of colloidal particles under external forces can be helpful in developing soft nanomaterials with novel functionalities. How external impurities organize within such confined systems is of fundamental and technological interest, especially when the system sizes are so small that even a single dopant can interact with an appreciable fraction of the system. To address this question, we use a defocused laser beam to form two-dimensional colloidal crystallites containing foreign dopants. Our studies reveal a surprising position dependence in the fate of dopants getting either spontaneously expelled or permanently internalized within the crystallite. This phenomenon arises due to the subtle interplay between the effects of external confinement and the role of entropy in the thermodynamics of small assemblies of interacting particles.

Interpolation method for pair correlations in classical crystals

Effects of anharmonicity on the pair correlation function of classical crystals are studied. The recently proposed shortest-graph approach using the Gaussian representation of the individual correlation peaks (the peak width is determined by the length of the shortest graph connecting a given pair of particles) is further improved, to account for anharmonic corrections due to finite temperatures and hard-sphere-like interactions. Two major effects are identified, leading to a modification of the correlation peaks at large or short distances: (i) the peaks at large distances, well described by Gaussians, should be calculated from the finite-temperature phonon spectra; (ii) at short distances, the correlation peaks deviate significantly from the Gaussian form due to the lattice discreteness. We propose the analytical interpolation method, based on the shortest-graph approach, which includes both effects. By employing the molecular dynamics simulations, the accuracy of the method is verified for three- and two-dimensional crystals with the Yukawa, inverse-power-law, and pseudo-hard-sphere pair interactions. The capabilities of the method are demonstrated by calculating the phase diagram of a three-dimensional Yukawa system.

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.

How do polydisperse repulsive colloids crystallize?

A modified version of the Gibbs-ensemble Monte-Carlo method reveals how polydisperse charged colloidal particles can build complex colloidal crystals. It provides general rules that are applicable to this fractionated crystallization that stems from size segregation. It explains the spontaneous formation of complex crystals with very large unit-cells in suspensions of nanoparticles with a broad size distribution.

Melting in three-dimensional and two-dimensional Yukawa systems

Solid-liquid phase transitions in three-dimensional (3D) and two-dimensional (2D) Yukawa systems were studied numerically and analytically, including the melting of the fcc and bcc 3D lattices, and of a hexagonal primitive (hp) 2D lattice. An approach is proposed for the determination of the melting lines in these systems. The suggested approach takes into account the nonlinearity (anharmonicity) of pair interaction forces and allows one to correctly predict the conditions of melting for 3D and 2D crystal systems. The obtained results are compared with the existing theoretical and numerical data.

Interplay of anisotropy in shape and interactions in charged platelet suspensions

Motivated by the intriguing phase behavior of charged colloidal platelets, we investigate the structure and dynamics of charged repulsive disks by means of Monte Carlo simulations. The electrostatic interactions are taken into account through an effective two-body potential, obtained within the nonlinear Poisson-Boltzmann formalism, which has the form of anisotropic screened Coulomb potential. Recently, we showed that the original intrinsic anisotropy of the electrostatic potential in competition with excluded volume effects leads to a rich phase behavior that not only includes various liquid-crystalline phases but also predicts the existence of novel structures composed of alternating nematic-antinematic sheets. Here, we examine the structural and dynamical signatures of each of the observed structures for both translational and rotational degrees of freedom. Finally, we discuss the influence of effective charge value and our results in relation to experimental findings on charged platelet suspensions.

Crystallization and reentrant melting of charged colloids in nonpolar solvents

We explore the crystallization of charged colloidal particles in a nonpolar solvent mixture. We simultaneously charge the particles and add counterions to the solution with aerosol-OT (AOT) reverse micelles. At low AOT concentrations, the charged particles crystallize into body-centered-cubic (bcc) or face-centered-cubic (fcc) Wigner crystals; at high AOT concentrations, the increased screening drives a thus far unobserved reentrant melting transition. We observe an unexpected scaling of the data with particle size, and account for all behavior with a model that quantitatively predicts both the reentrant melting and the data collapse.

Fluid approach to evaluate sound velocity in Yukawa systems and complex plasmas

The conventional fluid description of multicomponent plasma, supplemented by an appropriate equation of state for the macroparticle component, is used to evaluate the longitudinal sound velocity of Yukawa fluids. The obtained results are in very good agreement with those obtained earlier employing the quasilocalized charge approximation and molecular dynamics simulations in a rather broad parameter regime. Thus, a simple yet accurate tool to estimate the sound velocity across coupling regimes is proposed, which can be particularly helpful in estimating the dust-acoustic velocity in strongly coupled dusty (complex) plasmas. It is shown that, within the present approach, the sound velocity is completely determined by particle-particle correlations and the neutralizing medium (plasma), apart from providing screening of the Coulomb interaction, has no other effect on the sound propagation. The ratio of the actual sound velocity to its "ideal gas" (weak coupling) scale only weakly depends on the coupling strength in the fluid regime but exhibits a pronounced decrease with the increase of the screening strength. The limitations of the present approach in applications to real complex plasmas are briefly discussed.

Anisotropic elasticity of experimental colloidal Wigner crystals

Colloidal particles interacting via a long-range repulsion can, in contrast to hard-sphere systems, exhibit crystalline ordering at low volume fraction. Here we experimentally investigate the structure and properties of such "colloidal Wigner crystals." We find a body-centered-cubic crystalline phase at volume fractions of ϕ≳15%, which exhibits large fluctuations of individual particles from their average positions. We determine the three independent crystalline elastic constants and find that these crystals are very compliant and highly anisotropic.

Practical expressions for the internal energy and pressure of Yukawa fluids

Simple practical expressions that allow estimation of thermodynamic properties of Yukawa fluids in a wide range of coupling, up to the fluid-solid phase transition, are presented. These expressions demonstrate excellent agreement with the available results from numerical simulations. The approach provides simple and accurate tools to estimate thermodynamic properties of Yukawa fluids and related systems in a broad range of parameters.

Exponential approximation for one-component Yukawa plasma

A theory based on the exponential approximation of the liquid-state theory is applied to study properties of several models of one-component Yukawa plasma characterized by different values of the screening parameter z. The results of the new theory are compared to the results of a conventional theory, which is based on the first-order mean spherical approximation, and to the results of a Monte Carlo simulation. The new theory shows improvements in the predictions for the thermodynamic and structural properties of Yukawa plasmas with high and intermediate values of the screening parameter, z, and coupling parameter, Γ. For low values of z and Γ, the new theory is comparable in accuracy to the conventional theory, which in turn agrees well with the results of the Monte Carlo simulation.

Effective potential and interdiffusion in binary ionic mixtures

We calculate interdiffusion coefficients in a two-component, weakly or strongly coupled ion plasma (gas or liquid, composed of two ion species immersed into a neutralizing electron background). We use an effective potential method proposed recently by Baalrud and Daligaut [Phys. Rev. Lett. 110, 235001 (2013)]. It allows us to extend the standard Chapman-Enskog procedure of calculating the interdiffusion coefficients to the case of strong Coulomb coupling. We compute binary diffusion coefficients for several ionic mixtures and fit them by convenient expressions in terms of the generalized Coulomb logarithm. These fits cover a wide range of plasma parameters spanning from weak to strong Coulomb couplings. They can be used to simulate diffusion of ions in ordinary stars as well as in white dwarfs and neutron stars.

Design principles for Bernal spirals and helices with tunable pitch

Using the framework of potential energy landscape theory, we describe two in silico designs for self-assembling helical colloidal superstructures based upon dipolar dumbbells and Janus-type building blocks, respectively. Helical superstructures with controllable pitch length are obtained using external magnetic field driven assembly of asymmetric dumbbells involving screened electrostatic as well as magnetic dipolar interactions. The pitch of the helix is tuned by modulating the Debye screening length over an experimentally accessible range. The second design is based on building blocks composed of rigidly linked spheres with short-range anisotropic interactions, which are predicted to self-assemble into Bernal spirals. These spirals are quite flexible, and longer helices undergo rearrangements via cooperative, hinge-like moves, in agreement with experiment.

Crystallization kinetics of colloidal model suspensions: recent achievements and new perspectives

Colloidal model systems allow studying crystallization kinetics under fairly ideal conditions, with rather well-characterized pair interactions and minimized external influences. In complementary approaches experiment, analytic theory and simulation have been employed to study colloidal solidification in great detail. These studies were based on advanced optical methods, careful system characterization and sophisticated numerical methods. Over the last decade, both the effects of the type, strength and range of the pair-interaction between the colloidal particles and those of the colloid-specific polydispersity have been addressed in a quantitative way. Key parameters of crystallization have been derived and compared to those of metal systems. These systematic investigations significantly contributed to an enhanced understanding of the crystallization processes in general. Further, new fundamental questions have arisen and (partially) been solved over the last decade: including, for example, a two-step nucleation mechanism in homogeneous nucleation, choice of the crystallization pathway, or the subtle interplay of boundary conditions in heterogeneous nucleation. On the other hand, via the application of both gradients and external fields the competition between different nucleation and growth modes can be controlled and the resulting microstructure be influenced. The present review attempts to cover the interesting developments that have occurred since the turn of the millennium and to identify important novel trends, with particular focus on experimental aspects.