Equilibrium phase behavior of polydisperse hard spheres

Physicochemical properties and their impact on ice nucleation efficiency of respiratory viral RNA and proteins

Ice nucleation processes in the earth's atmosphere are critical for cloud formation, radiation, precipitation, and climate change. We investigated the physicochemical properties and ice nucleation potential of selected viral aerosols, including their RNA and proteins, using advanced techniques such as scanning-transmission electron microscopy (S/TEM), small angle X-ray scattering (SAXS), particle analyzers, and a peltier chamber. The experiments revealed that RNA particles obtained from MS2 bacteriophage had a mean freezing point of -13.9 ± 0.3 °C, comparable to the average ice nucleation temperature of global dust particles, which is approximatively -15 °C. RNA from MS2, Influenza, SARS-CoV-1 and SARS-CoV-2 demonstrated average ice nucleation temperatures of -13.9 ± 0.3 °C, -13.7 ± 0.3 °C, -13.7 ± 0.3 °C, and -15.9 ± 0.4 °C, respectively. SAXS analysis indicated a high local crystallinity value of 0.5 of MS2 RNA particles, hinting that high crystalline nature may contribute to their effectiveness as ice nuclei. Dilution experiments show that viral RNA consistently catalyzes ice nucleation. The addition of dust-containing particles, such as Fe2O3, CuO, and TiO2, to MS2 bacteriophage droplets enhanced ice nucleation, as did UV radiation. We herein discuss the implications of this work on ice nucleation and freezing processes.

On the Role of Starchy Grains in Ice Nucleation Processes

Little is known about the role of starchy food on climate change processes like ice nucleation. Here, we investigate the ice nucleation efficiency (INE) of eight different starchy food materials, namely, corn (CO), potato (PO), barley (BA), brown rice (BR), white rice (WR), oats (OA), wheat (WH), and sweet potato (SP), in immersion freezing mode under mixed-phase cloud conditions. Notably, among all these food materials, PO and BA exhibit the highest ice nucleation efficiency with ice nucleation temperatures as high as -4.3 °C (T50 ∼ -7.0 ± 0.5 °C) and -6.5 °C (T50 ∼ -7.2 ± 0.2 °C), respectively. We also explore the effect of environmentally relevant physicochemical conditions on ice nucleation efficiency, including different pH, temperature, UV/O3/NOx exposure, and various cocontaminants. The change in shape, size, surface properties, hydrophobicity, and crystallinity of materials accounted for the altered INE. The increase in shape, size, and hydrophobicity of the sample generally reduces the INE, whereas an increase in crystallinity enhances the INE of the sample under our experimental conditions. The results suggest that environmentally relevant concentrations slightly alter INE, indicating their role as catalysts in environmental matrices. The outcome of studies on the ice nucleation properties of these food-containing aerosols might help in the physicochemical understanding of other biomolecule-induced ice nucleation, which is still an underdeveloped research area.

Topological phases in nanoparticle monolayers: can crystalline, hexatic, and isotropic-fluid phases coexist in the same monolayer?

Topological phases are stable configurations of matter in 2-dimensions (2D) formed via spontaneous symmetry breaking. These play a crucial role in determining the system properties. Though a number of fundamental studies on topological phase transitions and topological defect dynamics have been conducted with model colloidal systems (typically microns in size), the same is lacking on nanoparticle monolayers (NPMLs, typically made of ligand-coated sub-ten nanometer particles). Here, we show that in an evaporation-driven self-assembly process, the three topological phases, namely crystalline, hexatic, and isotropic-fluid phases, can coexist within the same NPML. We associate this coexistence with the local variation in particle size, which can be described by a unique frequency parameter (p25), quantifying the fraction of NPs that has size deviation greater than or equal to 25% of the mean size (where the deviation,ζ is defined as ζ = ((|Size-mean|)/mean)). The p25-values for the three phases are distinctly different: crystalline arrangement occurs when p25 < ∼0.02, while a hexatic phase exists for 0.02 ≤ p25 ≤ 0.1. For p25 ≥ 0.1, the isotropic-fluid phase occurs. Following KTHNY-theory, we further numerically extrapolate the occurrence of each phase to the accumulated excess planar strain in the NPML due to the presence of various topological defects in the structures.

Even Strong Energy Polydispersity Does Not Affect the Average Structure and Dynamics of Simple Liquids

Size-polydisperse liquids have become standard models for avoiding crystallization, thereby enabling studies of supercooled liquids and glasses formed, e.g., by colloidal systems. Purely energy-polydisperse liquids have been studied much less, but provide an interesting alternative. We here study numerically the difference in structure and dynamics obtained by introducing these two kinds of polydispersity into systems of particles interacting via the Lennard-Jones and EXP pair potentials. To a very good approximation, the average pair structure and dynamics are unchanged even for strong energy polydispersity, which is not the case for size-polydisperse systems. When the system at extreme energy polydispersity undergoes a continuous phase separation into lower and higher particle-energy regions whose structure and dynamics are different from the average, the average structure and dynamics are still virtually the same as for the monodisperse system. Our findings are consistent with the fact that the distribution of forces on the individual particles do not change when energy polydispersity is introduced, while they do change in the case of size polydispersity. A theoretical explanation remains to be found, however.

Approaching the hard sphere limit in colloids suitable for confocal microscopy - the end of a decade lasting quest

PMMA-PHSA particles serve as the hard sphere model system since the 1980s. We investigate the fluid structure of fluorescent ones in three different solvents by laser scanning confocal microscopy: a decalin-tetrachloroethylene (TCE)-mixture and a decalin-cyclohexylbromide (CHB)-mixture with and without tetrabutylammoniumbromide (TBAB). The experimental 3D radial distribution functions are modeled by analytical theory and computer simulations taking polydispersity and the experimental position uncertainty into account. The quantitative comparison between experiment and simulation/theory establishes hard sphere like behavior for particles in decalin-TCE for a wide range of particle packing fractions. To the best of our knowledge, we present the first experimental dataset of a fluid structure that agrees convincingly with Percus-Yevick over a wide concentration range. Furthermore, charged sphere behavior is confirmed both for the decalin-CHB and the decalin-CHB-TBAB solvents, and it is demonstrated that a finite particle concentration reduces screening in the decalin-CHB-TBAB system compared to the bulk solvent.

Structure ordering and glass transition in size-asymmetric ternary mixtures of hard spheres: Variation from fragile to strong glasses

We investigate the structure and activated dynamics of a binary mixture of colloidal particles dispersed in a solvent of much smaller-sized particles. The solvent degrees of freedom are traced out from the grand partition function of the colloid-solvent mixture which reduces the system from ternary to effective binary mixture of colloidal particles. In the effective binary mixture colloidal particles interact via effective potential that consists of bare potential plus the solvent-induced interaction. Expressions for the effective potentials and pair correlation functions are derived. We used the result of pair correlation functions to determine the number of particles in a cooperatively reorganizing cluster (CRC) in which localized particles form "long-lived" nonchemical bonds with the central particle. For an event of relaxation to take place these bonds have to reorganize irreversibly, the energy involved in the processes is the effective activation energy of relaxation. Results are reported for hard sphere colloidal particles dispersed in a solvent of hard sphere particles. Our results show that the concentration of solvent can be used as a control parameter to fine-tune the microscopic structural ordering and the size of CRC that governs the glassy dynamics. We show that a small variation in the concentration of solvent creates a bigger change in the kinetic fragility which highlights a wide variation in behavior, ranging from fragile to strong glasses. We conclude that the CRC which is determined from the static pair correlation function and the fluctuations embedded in the system is probably the sole player in the physics of glass transition.

Sedimentation path theory for mass-polydisperse colloidal systems

Both polydispersity and the presence of a gravitational field are inherent to essentially any colloidal experiment. While several theoretical works have focused on the effect of polydispersity on the bulk phase behavior of a colloidal system, little is known about the effect of a gravitational field on a polydisperse colloidal suspension. We extend here the sedimentation path theory to study sedimentation-diffusion-equilibrium of a mass-polydisperse colloidal system: the particles possess different buoyant masses but they are otherwise identical. The model helps to understand the interplay between gravity and polydispersity on sedimentation experiments. Since the theory can be applied to any parent distribution of buoyant masses, it can also be used to study the sedimentation of monodisperse colloidal systems. We find that mass-polydispersity has a strong influence in colloidal systems near density matching for which the bare density of the colloidal particles equals the solvent density. To illustrate the theory, we study crystallization in sedimentation-diffusion-equilibrium of a suspension of mass-polydisperse hard spheres.

2D Colloidal Crystals with Anisotropic Impurities

We report a study of 2D colloidal crystals with anisotropic ellipsoid impurities using video microscopy. It is found that at low impurity densities, the impurity particles behave like floating disorder with which the quasi-long-range orientational order survives and the elasticity of the system is actually enhanced. There is a critical impurity density above which the 2D crystal loses the quasi-long-range orientational order. At high impurity densities, the 2D crystal breaks into polycrystalline domains separated by grain boundaries where the impurity particles aggregate. This transition is accompanied by a decrease in the elastic moduli, and it is associated with strong heterogeneous dynamics in the system. The correlation length vs impurity density in the disordered phase exhibits an essential singularity at the critical impurity density.

Vitrification is a spontaneous non-equilibrium transition driven by osmotic pressure

Persistent dynamics in colloidal glasses suggest the existence of a non-equilibrium driving force for structural relaxation during glassy aging. But the implicit assumption in the literature that colloidal glasses form within the metastable state bypasses the search for a driving force for vitrification and glassy aging and its connection with a metastable state. The natural relation of osmotic pressure to number-density gradients motivates us to investigate the osmotic pressure as this driving force. We use dynamic simulation to quench a polydisperse hard-sphere colloidal liquid into the putative glass region while monitoring structural relaxation and osmotic pressure. Following quenches to various depths in volume fractionϕ(whereϕ_RCP≈ 0.678 for 7% polydispersity), the osmotic pressure overshoots its metastable value, then decreases with age toward the metastable pressure, driving redistribution of coordination number and interparticle voids that smooths structural heterogeneity with age. For quenches to 0.56 ⩽ϕ⩽ 0.58, accessible post-quench volume redistributes with age, allowing the glass to relax into a strong supercooled liquid and easily reach a metastable state. At higher volume fractions, 0.59 ⩽ϕ< 0.64, this redistribution encounters a barrier that is subsequently overcome by osmotic pressure, allowing the system to relax toward the metastable state. But forϕ⩾ 0.64, the overshoot is small compared to the high metastable pressure; redistribution of volume stops as particles acquire contacts and get stuck, freezing the system far from the metastable state. Overall, the osmotic pressure drives structural rearrangements responsible for both vitrification and glassy age-relaxation. The connection of energy, pressure, and structure identifies the glass transition, 0.63 <ϕ_g⩽ 0.64. We leverage the connection of osmotic pressure to energy density to put forth the mechanistic view that relaxation of structural heterogeneity in colloidal glasses occurs via individual particle motion driven by osmotic pressure, and is a spontaneous energy minimization process that drives the glass off and back to the metastable state.

Crystal-to-Crystal Transitions in Binary Mixtures of Soft Colloids

In this article, we demonstrate a method for inducing reversible crystal-to-crystal transitions in binary mixtures of soft colloidal particles. Through a controlled decrease of salinity and increasingly dominating electrostatic interactions, a single sample is shown to reversibly organize into entropic crystals, electrostatic attraction-dominated crystals, or aggregated gels, which we quantify using microscopy and image analysis. We furthermore analyze crystalline structures with bond order analysis to discern between two crystal phases. We observe the different phases using a sample holder geometry that allows both in situ salinity control and imaging through confocal laser scanning microscopy and apply a synthesis method producing particles with high resolvability in microscopy with control over particle size. The particle softness provides for an enhanced crystallization speed, while altering the re-entrant melting behavior as compared to hard sphere systems. This work thus provides several tools for use in the reproducible manufacture and analysis of binary colloidal crystals.

Effect of particle size distribution on polydisperse hard disks

Using Monte Carlo simulations, we systematically investigate the effect of particle size distribution on the phase behavior of polydisperse hard disks. Compared with the commonly used Gaussian-like polydisperse hard disks [P. Sampedro Ruiz, Q.-l. Lei, and R. Ni, Commun. Phys. 2, 70 (2019)], we find that the phase behavior of polydisperse hard-disk systems with lognormal and triangle distributions is significantly different. In polydisperse hard-disk systems of lognormal distributions, although the phase diagram appears similar to that of Gaussian-like polydisperse hard disks, the re-entrant melting of the hexatic or solid phase cannot be observed in sedimentation experiments. For polydisperse hard-disk systems of triangle distributions, the phase behavior is qualitatively different from the Gaussian-like and lognormal distributions, and we cannot reach any system of true polydispersity larger than 0.06, which is due to the special shape of the triangle distribution. Our results suggest that the exact particle size distribution is of primary importance in determining the phase behavior of polydisperse hard disks, and we do not have a universal phase diagram for different polydisperse hard-disk systems.

Slowing down of dynamics and orientational order preceding crystallization in hard-sphere systems

Despite intensive studies in the past decades, the local structure of disordered matter remains widely unknown. We show the results of a coherent x-ray scattering study revealing higher-order correlations in dense colloidal hard-sphere systems in the vicinity of their crystallization and glass transition. With increasing volume fraction, we observe a strong increase in correlations at both medium-range and next-neighbor distances in the supercooled state, both invisible to conventional scattering techniques. Next-neighbor correlations are indicative of ordered precursor clusters preceding crystallization. Furthermore, the increase in such correlations is accompanied by a marked slowing down of the dynamics, proving experimentally a direct relation between orientational order and sample dynamics in a soft matter system. In contrast, correlations continuously increase for nonequilibrated, glassy samples, suggesting that orientational order is reached before the sample slows down to reach (quasi-)equilibrium.

Synchronized fractionation and phase separation in binary colloids

Fractionation is necessary for self-assembly in multicomponent mixtures. Here, reversible fractionation and crystallization are realized and studied in two-dimensional binary colloids which are supersaturated by enhancing the attraction between colloidal particles. As a deep and fast supersaturation results in gels with a uniform distribution of binary particles, a gradual quasistatic supersaturation process leads to a two-step crystallization in which small particles and large particles are fractionated as coexisting crystal and liquid phases respectively. Fractionation occurs as well in the quasistatic melting of gels. We show that the synchronized fractionation and phase separation arises from the competition between the size-dependent repulsion and the tunable attraction. The results in this study demonstrate a robust mechanism of fractionation via phase separation, and have important implication in understanding the reversible formation of membraneless organelles in living cells.

Dense hard disk ordering: influence of bidispersity and quenched disorder

Using Monte Carlo simulations, the impact on structural ordering in two-dimensional systems via the interplay of size bidispersity and quenched disorder in the form of an externally applied spatially random potential, is studied for a system of hard disks. By scanning across a wide range of dense packing fractions, size ratios and roughness of the applied potential, the phase diagram is constructed, which demonstrates that both quenched and size disorders shift the onset of translational order to higher packings, while maintaining the presence of the intermediate hexatic phase. At larger disorder strengths, the signatures of structural order are absent within the range of investigated packing fractions. Further, the dynamics with increasing potential strength is analysed for the mono-component system to obtain a spatio-temporal description of the melting process. Finally, the influence of the externally rough field on the Mermin-Wagner fluctuations, characteristic to two-dimensional systems, is investigated.

Strong size selectivity in the self-assembly of rounded nanocubes into 3D mesocrystals

The self-assembly of nanoparticles into highly ordered crystals is largely influenced by variations in the size and shape of the constituent particles, with crystallization generally not observed if their polydispersity is too large. Here, we report on size selectivity in the self-assembly of rounded cubic maghemite nanoparticles into three-dimensional mesocrystals. Different X-ray scattering techniques are used to study and compare a nanoparticle dispersion that is used later for self-assembly, an ensemble of mesocrystals grown on a substrate, as well as an individual mesocrystal. The individual μm-sized mesocrystal is isolated using a focused-ion-beam-based technique and investigated by the diffraction of a micro-focused X-ray beam. Structural analysis reveals that individual mesocrystals have a drastically smaller size dispersity of nanoparticles than that in the initial dispersion, implying very strong size selectivity during self-assembly. The small size dispersity of the nanoparticles within individual mesocrystals is accompanied by a very narrow lattice parameter distribution. In contrast, the lattice parameter distribution within all mesocrystals of an ensemble is about four times wider than that of individual mesocrystals, indicating significant size fractionalization between mesocrystals during self-assembly. The small size dispersity within each mesocrystal has important implications for their physical properties.

Systems with Size and Energy Polydispersity: From Glasses to Mosaic Crystals

We use Langevin dynamics simulations to study dense 2d systems of particles with both size and energy polydispersity. We compare two types of bidisperse systems which differ in the correlation between particle size and interaction parameters: in one system big particles have high interaction parameters and small particles have low interaction parameters, while in the other system the situation is reversed. We study the different phases of the two systems and compare them to those of a system with size but not energy bidispersity. We show that, depending on the strength of interaction between big and small particles, cooling to low temperatures yields either homogeneous glasses or mosaic crystals. We find that systems with low mixing interaction, undergo partial freezing of one of the components at intermediate temperatures, and that while this phenomenon is energy-driven in both size and energy bidisperse systems, it is controlled by entropic effects in systems with size bidispersity only.

Packing Polydisperse Colloids into Crystals: When Charge-Dispersity Matters

Monte Carlo simulations, fully constrained by experimental parameters, are found to agree well with a measured phase diagram of aqueous dispersions of nanoparticles with a moderate size polydispersity over a broad range of salt concentrations, c_{s}, and volume fractions, ϕ. Upon increasing ϕ, the colloids freeze first into coexisting compact solids then into a body centered cubic phase (bcc) before they melt into a glass forming liquid. The surprising stability of the bcc solid at high ϕ and c_{s} is explained by the interaction (charge) polydispersity and vibrational entropy.

Molecular heterogeneity drives reconfigurable nematic liquid crystal drops

With few exceptions1-3, polydispersity or molecular heterogeneity in matter tends to impede self-assembly and state transformation. For example, shape transformations of liquid droplets with monodisperse ingredients have been reported in equilibrium4-7 and non-equilibrium studies8,9, and these transition phenomena were understood on the basis of homogeneous material responses. Here, by contrast, we study equilibrium suspensions of drops composed of polydisperse nematic liquid crystal oligomers (NLCOs). Surprisingly, molecular heterogeneity in the polydisperse drops promotes reversible shape transitions to a rich variety of non-spherical morphologies with unique internal structure. We find that variation of oligomer chain length distribution, temperature, and surfactant concentration alters the balance between NLCO elastic energy and interfacial energy, and drives formation of nematic structures that range from roughened spheres to 'flower' shapes to branched filamentous networks with controllable diameters. The branched structures with confined liquid crystal director fields can be produced reversibly over areas of at least one square centimetre and can be converted into liquid crystal elastomers by ultraviolet curing. Observations and modelling reveal that chain length polydispersity plays a crucial role in driving these morphogenic phenomena, via spatial segregation. This insight suggests new routes for encoding network structure and function in soft materials.

Refractive index matched, nearly hard polymer colloids

Refractive index matched particles serve as essential model systems for colloid scientists, providing nearly hard spheres to explore structure and dynamics. The poly(methyl methacrylate) latexes typically used are often refractive index matched by dispersing them in binary solvent mixtures, but this can lead to undesirable changes, such as particle charging or swelling. To avoid these shortcomings, we have synthesized refractive index matched colloids using polymerization-induced self-assembly (PISA) rather than as polymer latexes. The crucial difference is that these diblock copolymer nanoparticles consist of a single core-forming polymer in a single non-ionizable solvent. The diblock copolymer chosen was poly(stearyl methacrylate)-poly(2,2,2-trifluoroethyl methacrylate) (PSMA-PTFEMA), which self-assembles to form PTFEMA core spheres in n-alkanes. By monitoring scattered light intensity, n-tetradecane was found to be the optimal solvent for matching the refractive index of such nanoparticles. As expected for PISA syntheses, the diameter of the colloids can be controlled by varying the PTFEMA degree of polymerization. Concentrated dispersions were prepared, and the diffusion of the PSMA-PTFEMA nanoparticles as a function of volume fraction was measured. These diblock copolymer nanoparticles are a promising new system of transparent spheres for future colloidal studies.

Entropy-Induced Self-Assembly of Colloidal Crystals with High Reflectivity and Narrow Reflection Bandwidth

Cracks and defects, which could result in lower reflectivity and larger full width at half maximum (FWHM), are the major obstacles for obtaining highly ordered structures of colloidal crystals (CCs). The high-quality CCs with high reflectivity (more than 90%) and 9.2 nm narrow FWHM have been successfully fabricated using a fixed proportion of a soft matter system composed of silica particles (SPs), polyethylene glycol diacrylate (PEGDA), and ethanol. The influences of refractivity difference, volume fractions, and particle dimension on FWHM were illuminated. Firstly, we clarified the influences of the planar interface and the bending interface on the self-assembly. The CCs had been successfully fabricated on the planar interface and presented unfavorable results on the bending interface. Secondly, a hard sphere system consisting of SPs, PEGDA, and ethanol was established, and the entropy-driven phase transition mechanism of a polydisperse system was expounded. The FWHM and reflectivity of CCs showed an increasing trend with increasing temperature. Consequently, high-quality CCs were obtained by adjusting temperatures (ordered structure formed at 90 °C and solidified at 0 °C) based on the surface phase rule of the system. We acquired a profound understanding of the principle and process of self-assembly, which is significant for preparation and application of CCs such as optical filters.

Nonisotropic Self-Assembly of Nanoparticles: From Compact Packing to Functional Aggregates

Quantum strongly correlated systems that exhibit interesting features in condensed matter physics often need an unachievable temperature or pressure range in classical materials. One solution is to introduce a scaling factor, namely, the lattice parameter. Synthetic heterostructures named superlattices or supracrystals are synthesized by the assembling of colloidal atoms. These include semiconductors, metals, and insulators for the exploitation of their unique properties. Most of them are currently limited to dense packing. However, some of desired properties need to adjust the colloidal atoms neighboring number. Here, the current state of research in nondense packing is summarized, discussing the benefits, outlining possible scenarios and methodologies, describing examples reported in the literature, briefly discussing the challenges, and offering preliminary conclusions. Penetrating such new and intriguing research fields demands a multidisciplinary approach accounting for the coupling of statistic physics, solid state and quantum physics, chemistry, computational science, and mathematics. Standard interactions between colloidal atoms and emerging fields, such as the use of Casimir forces, are reported. In particular, the focus is on the novelty of patchy colloidal atoms to meet this challenge.

Formation of Laves phases in buoyancy matched hard sphere suspensions

Colloidal Laves phases (LPs) are promising precursors for photonic materials. Laves phases have not yet been observed to form in experiments on colloidal suspensions of hard spheres (HS), even though they have been reported in computer simulations. LP formation so far has been achieved only for binary mixtures of colloidal charged spheres or ligand-stabilized nano-particles after drying. Using static light scattering, we monitored LP formation and annealing in a binary mixture of buoyant hard sphere approximants (size ratio Γ = 0.77, number or molar fraction of small spheres xS = 0.76) for volume fractions in the fluid-crystal coexistence regions. All samples spontaneously formed MgZn2 type LPs on the time scale of weeks to months via bulk nucleation and growth. Irrespective of the initial suspension volume fractions, the LP volume fraction at coexistence is ΦCOEX = 0.59 which is significantly below the close packing limit ΦMAX = 0.615 and remarkably close to the expectation from simulation. At low volume fractions, crystals anneal to high quality during coarsening which is in line with recent theoretical expectations for the thermodynamic stability of different LP types. At large volume fractions, however, the diffractograms evolve towards a more MgCu2-like appearance which we attribute to the formation of randomly stacked LPs. Such structures are not known from atomic systems.

Domain size polydispersity effects on the structural and dynamical properties in lipid monolayers with phase coexistence

In lipid monolayers with phase coexistence, domains of the liquid-condensed phase always present size polydispersity. However, very few theoretical works consider size distribution effects on the monolayer properties. Because of the difference in surface densities, domains have excess dipolar density with respect to the surrounding liquid expanded phase, originating a dipolar inter-domain interaction. This interaction depends on the domain area, and hence the presence of a domain size distribution is associated with interaction polydispersity. Inter-domain interactions are fundamental to understanding the structure and dynamics of the monolayer. For this reason, it is expected that polydispersity significantly alters monolayer properties. By means of Brownian dynamics simulations, we study the radial distribution function (RDF), the average mean square displacement and the average time-dependent self-diffusion coefficient, D(t), of lipid monolayers with normally distributed size domains. For this purpose, we vary the relevant system parameters, polydispersity and interaction strength, within a range of experimental interest. We also analyze the consequences of using a monodisperse model to determine the interaction strength from an experimental RDF. We find that polydispersity strongly affects the value of the interaction strength, which is greatly underestimated if polydispersity is not considered. However, within a certain range of parameters, the RDF obtained from a polydisperse model can be well approximated by that of a monodisperse model, by suitably fitting the interaction strength, even for 40% polydispersities. For small interaction strengths or small polydispersities, the polydisperse systems obtained from fitting the experimental RDF have an average mean square displacement and D(t) in good agreement with that of the monodisperse system.

Determination of the local density of polydisperse nanoparticle assemblies

Quantitative characterization of the average structure of dense nanoparticle assemblies and aggregates is a common problem in nanoscience. Small-angle scattering is a suitable technique, but it is usually limited to not too big assemblies due to the limited experimental range, low concentrations to avoid interactions, and monodispersity to keep calculations tractable. In the present paper, a straightforward analysis of the generally available scattered intensity - even for large assemblies, at high concentrations - is detailed, providing information on the local volume fraction of polydisperse particles with hard sphere interactions. It is based on the identical local structure of infinite homogeneous nanoparticle assemblies and their subsets forming finite-sized clusters. This approach is extended to polydispersity, using Monte-Carlo simulations of hard and moderately sticky hard spheres. As a result, a simple relationship between the observed structure factor minimum - termed the correlation hole - and the average local volume fraction κ on the scale of neighboring particles is proposed and validated through independent aggregate simulations. This relationship shall be useful as an efficient tool for the structural analysis of arbitrarily aggregated colloidal systems.

Phase separation dynamics of polydisperse colloids: a mean-field lattice-gas theory

Strong theoretical evidence shows that dense colloidal mixtures phase-separate in two stages and the denser phase contains long-lived composition heterogeneities.

Dynamics of highly polydisperse colloidal suspensions as a model system for bacterial cytoplasm

There are various kinds of macromolecules in bacterial cell cytoplasm. The size polydispersity of the macromolecules is so significant that the crystallization and the phase separation could be suppressed, thus stabilizing the liquid state of bacterial cytoplasm. On the other hand, recent experiments suggested that the macromolecules in bacterial cytoplasm should exhibit glassy dynamics, which should be also affected significantly by the size polydispersity of the macromolecules. In this work, we investigate the anomalous and slow dynamics of highly polydisperse colloidal suspensions, of which size distribution is chosen to mimic Escherichia coli cytoplasm. We find from our Langevin dynamics simulations that the diffusion coefficient (D_{tot}) and the displacement distribution functions (P(r,t)) averaged over all colloids of different sizes do not show anomalous and glassy dynamic behaviors until the system volume fraction ϕ is increased up to 0.82. This indicates that the intrinsic polydispersity of bacterial cytoplasm should suppress the glass transition and help maintain the liquid state of the cytoplasm. On the other hand, colloids of each kind show totally different dynamic behaviors depending on their size. The dynamics of colloids of different size becomes non-Gaussian at a different range of ϕ, which suggests that a multistep glass transition should occur. The largest colloids undergo the glass transition at ϕ=0.65, while the glass transition does not occur for smaller colloids in our simulations even at the highest value of ϕ. We also investigate the distribution (P(θ,t)) of the relative angles of displacement for macromolecules and find that macromolecules undergo directionally correlated motions in a sufficiently dense system.

Self-Assembly of Colloidal Nanocrystals: From Intricate Structures to Functional Materials

Chemical methods developed over the past two decades enable preparation of colloidal nanocrystals with uniform size and shape. These Brownian objects readily order into superlattices. Recently, the range of accessible inorganic cores and tunable surface chemistries dramatically increased, expanding the set of nanocrystal arrangements experimentally attainable. In this review, we discuss efforts to create next-generation materials via bottom-up organization of nanocrystals with preprogrammed functionality and self-assembly instructions. This process is often driven by both interparticle interactions and the influence of the assembly environment. The introduction provides the reader with a practical overview of nanocrystal synthesis, self-assembly, and superlattice characterization. We then summarize the theory of nanocrystal interactions and examine fundamental principles governing nanocrystal self-assembly from hard and soft particle perspectives borrowed from the comparatively established fields of micrometer colloids and block copolymer assembly. We outline the extensive catalog of superlattices prepared to date using hydrocarbon-capped nanocrystals with spherical, polyhedral, rod, plate, and branched inorganic core shapes, as well as those obtained by mixing combinations thereof. We also provide an overview of structural defects in nanocrystal superlattices. We then explore the unique possibilities offered by leveraging nontraditional surface chemistries and assembly environments to control superlattice structure and produce nonbulk assemblies. We end with a discussion of the unique optical, magnetic, electronic, and catalytic properties of ordered nanocrystal superlattices, and the coming advances required to make use of this new class of solids.

Effect of Energy Polydispersity on the Nature of Lennard-Jones Liquids

In the companion paper [ Ingebrigtsen , T. S. ; Tanaka , H. J. Phys. Chem. B 2015 , 119 , 11052 ] the effect of size polydispersity on the nature of Lennard-Jones (LJ) liquids, which represent most molecular liquids without hydrogen bonds, was studied. More specifically, it was shown that even highly size polydisperse LJ liquids are Roskilde-simple (RS) liquids. RS liquids are liquids with strong correlation between constant volume equilibrium fluctuations of virial and potential energy and are simpler than other types of liquids. Moreover, it was shown that size polydisperse LJ liquids have isomorphs to a good approximation. Isomorphs are curves in the phase diagram of RS liquids along which structure, dynamics, and some thermodynamic quantities are invariant in dimensionless (reduced) units. In this paper, we study the effect of energy polydispersity on the nature of LJ liquids. We show that energy polydisperse LJ liquids are RS liquids. However, a tendency of particle segregation, which increases with the degree of polydispersity, leads to a loss of strong virial-potential energy correlation but is mitigated by increasing temperature and/or density. Isomorphs are a good approximation also for energy polydisperse LJ liquids, although particle-resolved quantities display a somewhat poorer scaling compared to the mean quantities along the isomorph.

Hiding in Plain View: Colloidal Self-Assembly from Polydisperse Populations

We report small-angle x-ray scattering experiments on aqueous dispersions of colloidal silica with a broad monomodal size distribution (polydispersity, 14%; size, 8 nm). Over a range of volume fractions, the silica particles segregate to build first one, then two distinct sets of colloidal crystals. These dispersions thus demonstrate fractional crystallization and multiple-phase (bcc, Laves AB_{2}, liquid) coexistence. Their remarkable ability to build complex crystal structures from a polydisperse population originates from the intermediate-range nature of interparticle forces, and it suggests routes for designing self-assembling colloidal crystals from the bottom up.

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.

Gelation in mixtures of polymers and bidisperse colloids

We investigated the effects of varying the volume fraction of large particles (r) on the linear rheology and microstructure of mixtures of polymers and bidisperse colloids, in which the ratio of the small and large particle diameters was α=0.31 or α=0.45. Suspensions formulated at a total volume fraction of ϕ_{T}=0.15 and a constant concentration of polymer in the free volume c/c^{*}≈0.7 contained solid-like gels for small r and fluids or fluids of clusters at large r. The solid-like rheology and microstructure of these suspensions changed little with r when r was small, and fluidized only when r>0.8. By contrast, dense suspensions with ϕ_{T}=0.40 and α=0.31 contained solid-like gels at all concentrations of large particles and exhibited only modest rheological and microstructural changes upon varying the volume fraction of large particles. These results suggest that the effect of particle-size dispersity on the properties of colloid-polymer mixtures are asymmetric in particle size and are most pronounced near a gelation boundary.

Effect of Size Polydispersity on the Nature of Lennard-Jones Liquids

Polydisperse fluids are encountered everywhere in biological and industrial processes. These fluids naturally show a rich phenomenology exhibiting fractionation and shifts in critical point and freezing temperatures. We study here the effect of size polydispersity on the basic nature of Lennard-Jones (LJ) liquids, which represent most molecular liquids without hydrogen bonds, via two- and three-dimensional molecular dynamics computer simulations. A single-component liquid constituting spherical particles and interacting via the LJ potential is known to exhibit strong correlations between virial and potential energy equilibrium fluctuations at constant volume. This correlation significantly simplifies the physical description of the liquid, and these liquids are now known as Roskilde-simple (RS) liquids. We show that this simple nature of the single-component LJ liquid is preserved even for very high polydispersities (above 40% polydispersity for the studied uniform distribution). We also investigate isomorphs of moderately polydisperse LJ liquids. Isomorphs are curves in the phase diagram of RS liquids along which structure, dynamics, and some thermodynamic quantities are invariant in dimensionless units. We find that isomorphs are a good approximation even for polydisperse LJ liquids. The theory of isomorphs thus extends readily to size polydisperse fluids and can be used to improve even further the understanding of these intriguing systems.

On polydispersity and the hard sphere glass transition

We investigate the dynamics of polydisperse hard spheres at high packing fractions ϕ. We use extensive numerical simulations based on an experimentally-realistic particle size distribution (PSD) and compare to commonly-used PSDs such as Gaussian or top hat distribution. We find that the mode of kinetic arrest depends on the PSD's shape and not only on its variance. For the experimentally-realistic PSD we find ageing dynamics even though the density correlators decay fully to zero for ϕ ≥ 0.59. We observe substantial decoupling of the dynamics of the smallest and largest particles. While the smallest particles remain diffusive in all our simulations, a power-law describes the largest-particle diffusion, suggesting an ideal arrest at ϕc ∼ 0.588. The latter is however averted just before ϕc, due to the presence of the mobile smallest particles. In addition, we identify that a partial aging mechanism is at work, whose effects are most pronounced for the largest particles. By comparing our results with recent experimental observations of ergodic behavior up to ϕ ∼ 0.6 in a hard-sphere system, we argue that this is an effect of polydispersity, which smears out the glass transition.

Crystallization kinetics of colloidal binary mixtures with depletion attraction

In this work the crystallization kinetics of colloidal binary mixtures with attractive interaction potential (Asakura-Oosawa) has been addressed. Parameters such as fraction of crystals, linear crystal dimension and crystal packing have been quantified in order to understand how the crystal formation is driven in terms of the depth of the attractive potential and the composition of the binary mixture (described by the number ratio). It was found that inside the eutectic triangle, crystallization is mainly governed by nucleation and the crystal packing is close to the close-packing of hard spheres. Moving out from the eutectic triangle towards small component results in the crystallization of small spheres. Enrichment of the eutectic mixture with large component results in the crystallization of both large and small spheres, however, the kinetics are completely different from those of the eutectic composition. Crosslinked polystyrene microgels with nearly hard sphere interactions were used as model systems. Attraction was introduced by addition of linear polystyrene. The time evolution of crystallization has been followed by static light scattering.

Exposing a dynamical signature of the freezing transition through the sound propagation gap

The conventional view of freezing holds that nuclei of the crystal phase form in the metastable fluid through purely stochastic thermal density fluctuations. The possibility of a change in the character of the fluctuations as the freezing point is traversed is beyond the scope of this perspective. Here we show that this perspective may be incomplete by examination of the time autocorrelation function of the longitudinal current, computed by molecular dynamics for the hard-sphere fluid around its freezing point. In the spatial window where sound is overdamped, we identify a change in the long-time decay of the correlation function at the known freezing points of monodisperse and moderately polydisperse systems. The fact that these findings agree with previous experimental studies of colloidal systems in which particle are subject to diffusive dynamics, suggests that the dynamical signature we identify with the freezing transition is a consequence of packing effects alone.

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.

Prediction of polydisperse hard-sphere mixture behavior using tridisperse systems

How many state-variables are needed to predict the equation of state and the jamming density of polydisperse mixtures in glassy, non-equilibrium compressed states? We propose to define equivalent and maximally equivalent systems as those that match three and five moments of a given polydisperse size distribution, respectively. Fluids can be represented well by an equivalent system with only s = 2 components (bidisperse). As little as s = 3 components (tridisperse) are enough to achieve a maximally equivalent system. Those match macroscopic properties in glassy states, but also the volume fraction of rattlers, suggesting strong microstructural equivalency too. For many soft and granular systems, tridisperse, maximally equivalent systems allow for a closed analytical treatment and well-controlled industrial applications, while our proposal waits for experimental validation.

Solid-liquid transition in polydisperse Lennard-Jones systems

We study melting of a face-centered crystalline solid consisting of polydisperse Lennard-Jones spheres with Gaussian polydispersity in size. The phase diagram reproduces the existence of a nearly temperature invariant terminal polydispersity (δ(t) =/~ 0.11), with no signature of reentrant melting. The absence of reentrant melting can be attributed to the influence of the attractive part of the potential upon melting. We find that at terminal polydispersity the fractional density change approaches zero, which seems to arise from vanishingly small compressibility of the disordered phase. At constant temperature and volume fraction the system undergoes a sharp transition from crystalline solid to the disordered amorphous or fluid state with increasing polydispersity. This has been quantified by second- and third-order rotational invariant bond orientational order, as well as by the average inherent structure energy. The translational order parameter also indicates a similar sharp structural change at δ =/~ 0.09 in case of T(*) = 1.0, φ = 0.58. The free energy calculation further supports the sharp nature of the transition. The third-order rotationally invariant bond order shows that with increasing polydispersity, the local cluster favors a more icosahedral arrangement and the system loses its local crystalline symmetry. Interestingly, the value of structure factor S(k) of the amorphous phase at δ =/~ 0.10 (just beyond the solid-liquid transition density at T(*) = 1) becomes 2.75, which is below the value of 2.85 required for freezing given by the empirical Hansen-Verlet rule of crystallization, well known in the theory of freezing.

Reentrant disordered phase in a system of repulsive rods on a Bethe-like lattice

We solve exactly a model of monodispersed rigid rods of length k with repulsive interactions on the random locally tree-like layered lattice. For k≥4 we show that with increasing density, the system undergoes two phase transitions: first, from a low-density disordered phase to an intermediate density nematic phase and, second, from the nematic phase to a high-density reentrant disordered phase. When the coordination number is four, both phase transitions are continuous and in the mean field Ising universality class. For an even coordination number larger than four, the first transition is discontinuous, while the nature of the second transition depends on the rod length k and the interaction parameters.

Impurity partitioning during colloidal crystallization

We have found that an impurity partitioning takes place during growth of colloidal crystals, which was recognized by the fact that the impurity concentration in the solid (CS) was different from that in the initial solution (C0). The effective partition coefficient k(eff) (=CS/C0) was investigated for pure polystyrene and polystyrene dyed with fluorescent particles by changing the ratio of particle diameters d(imp)/d(cryst) and growth rate V. At each size ratio for the polystyrene impurity, k(eff) was less than unity and increased to unity with increasing V, whereas at a given growth rate, k(eff) increased to unity as d(imp)/d(cryst) approached unity. These results were consistent with the solute behavior analyzed using the Burton, Prim, and Slichter (BPS) model. The obtained k0, equilibrium partition coefficient, from a BPS plot increased as d(imp)/d(cryst) approached unity. In contrast, while the fluorescent particles also followed the BPS model, they showed higher k0 values than those of the same size of polystyrene particles. A k0 value greater than unity was obtained for impurities that were similar in size to the host particle. This behavior is attributed to the positive free energy of fusion associated with the incorporation of the fluorescent particles into the host matrix. Such positive free energy of fusion implies the presence of the enthalpy associated with interaction between particles.