Evaporation-induced assembly of colloidal crystals

Nanocrystal Assemblies: Current Advances and Open Problems

We explore the potential of nanocrystals (a term used equivalently to nanoparticles) as building blocks for nanomaterials, and the current advances and open challenges for fundamental science developments and applications. Nanocrystal assemblies are inherently multiscale, and the generation of revolutionary material properties requires a precise understanding of the relationship between structure and function, the former being determined by classical effects and the latter often by quantum effects. With an emphasis on theory and computation, we discuss challenges that hamper current assembly strategies and to what extent nanocrystal assemblies represent thermodynamic equilibrium or kinetically trapped metastable states. We also examine dynamic effects and optimization of assembly protocols. Finally, we discuss promising material functions and examples of their realization with nanocrystal assemblies.

Solvent Isotherms and Structural Transitions in Nanoparticle Superlattice Assembly

We introduce a Molecular Theory for Compressible Fluids (MOLT-CF) that enables us to compute free energies and other thermodynamic functions for nanoparticle superlattices with any solvent content, including the dry limit. Quantitative agreement is observed between MOLT-CF and united-atom molecular dynamics simulations performed to assess the reliability and precision of the theory. Among other predictions, MOLT-CF shows that the amount of solvent within the superlattice decreases approximately linearly with its vapor pressure and that in the late stages of drying, solvent-filled voids form at lattice interstitials. Applied to single-component superlattices, MOLT-CF predicts fcc-to-bcc Bain transitions for decreasing vapor pressure and for increasing ligand length, both in agreement with experimental results. We explore the stability of other single-component phases and show that the C14 Frank-Kasper phase, which has been reported in experiments, is not a global free-energy minimum. Implications for precise assembly and prediction of multicomponent nanoparticle systems are discussed.

Defects in lead halide perovskite light-emitting diodes under electric field: from behavior to passivation strategies

Lead halide perovskites (LHPs) are emerging semiconductor materials for light-emitting diodes (LEDs) owing to their unique structure and superior optoelectronic properties. However, defects that initiate degradation of LHPs through external stimuli and prompt internal ion migration at the interfaces remain a significant challenge. The electric field (EF), which is a fundamental driving force in LED operation, complicates the role of these defects in the physical and chemical properties of LHPs. A deeper understanding of EF-induced defect behavior is crucial for optimizing the LED performance. In this review, the origins and characterization of defects are explored, indicating the influence of EF-induced defect dynamics on LED performance and stability. A comprehensive overview of recent defect passivation approaches for LHP bulk films and nanocrystals (NCs) is also provided. Given the ubiquity of EF, a summary of the EF-induced defect behavior can enhance the performance of perovskite LEDs and related optoelectronic devices.

Structure Formation in Supraparticles Composed of Spherical and Elongated Particles

We studied the evaporation-induced formation of supraparticles from dispersions of elongated colloidal particles using experiments and computer simulations. Aqueous droplets containing a dispersion of ellipsoidal and spherical polystyrene particles were dried on superamphiphobic surfaces at different humidity values that led to varying evaporation rates. Supraparticles made from only ellipsoidal particles showed short-range lateral ordering at the supraparticle surface and random orientations in the interior regardless of the evaporation rate. Particle-based simulations corroborated the experimental observations in the evaporation-limited regime and showed an increase in the local nematic ordering as the diffusion-limited regime was reached. A thin shell of ellipsoids was observed at the surface when supraparticles were made from binary mixtures of ellipsoids and spheres. Image analysis revealed that the supraparticle porosity increased with an increasing aspect ratio of the ellipsoids.

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.

Dynamic density functional theory for drying colloidal suspensions: Comparison of hard-sphere free-energy functionals

Dynamic density functional theory (DDFT) is a promising approach for predicting the structural evolution of a drying suspension containing one or more types of colloidal particles. The assumed free-energy functional is a key component of DDFT that dictates the thermodynamics of the model and, in turn, the density flux due to a concentration gradient. In this work, we compare several commonly used free-energy functionals for drying hard-sphere suspensions, including local-density approximations based on the ideal-gas, virial, and Boublík–Mansoori–Carnahan–Starling–Leland (BMCSL) equations of state as well as a weighted-density approximation based on fundamental measure theory (FMT). To determine the accuracy of each functional, we model one- and two-component hard-sphere suspensions in a drying film with varied initial heights and compositions, and we compare the DDFT-predicted volume fraction profiles to particle-based Brownian dynamics (BD) simulations. FMT accurately predicts the structure of the one-component suspensions even at high concentrations and when significant density gradients develop, but the virial and BMCSL equations of state provide reasonable approximations for smaller concentrations at a reduced computational cost. In the two-component suspensions, FMT and BMCSL are similar to each other but modestly overpredict the extent of stratification by size compared to BD simulations. This work provides helpful guidance for selecting thermodynamic models for soft materials in nonequilibrium processes, such as solvent drying, solvent freezing, and sedimentation.

Modeling Solution Drying by Moving a Liquid-Vapor Interface: Method and Applications

A method of simulating the drying process of a soft matter solution with an implicit solvent model by moving the liquid-vapor interface is applied to various solution films and droplets. For a solution of a polymer and nanoparticles, we observe “polymer-on-top” stratification, similar to that found previously with an explicit solvent model. Furthermore, “polymer-on-top” is found even when the nanoparticle size is smaller than the radius of gyration of the polymer chains. For a suspension droplet of a bidisperse mixture of nanoparticles, we show that core-shell clusters of nanoparticles can be obtained via the “small-on-outside” stratification mechanism at fast evaporation rates. “Large-on-outside” stratification and uniform particle distribution are also observed when the evaporation rate is reduced. Polymeric particles with various morphologies, including Janus spheres, core-shell particles, and patchy particles, are produced from drying droplets of polymer solutions by combining fast evaporation with a controlled interaction between the polymers and the liquid-vapor interface. Our results validate the applicability of the moving interface method to a wide range of drying systems. The limitations of the method are pointed out and cautions are provided to potential practitioners on cases where the method might fail.

A Deep Learning Framework Discovers Compositional Order and Self-Assembly Pathways in Binary Colloidal Mixtures

Binary colloidal superlattices (BSLs) have demonstrated enormous potential for the design of advanced multifunctional materials that can be synthesized via colloidal self-assembly. However, mechanistic understanding of the three-dimensional self-assembly of BSLs is largely limited due to a lack of tractable strategies for characterizing the many two-component structures that can appear during the self-assembly process. To address this gap, we present a framework for colloidal crystal structure characterization that uses branched graphlet decomposition with deep learning to systematically and quantitatively describe the self-assembly of BSLs at the single-particle level. Branched graphlet decomposition is used to evaluate local structure via high-dimensional neighborhood graphs that quantify both structural order (e.g., body-centered-cubic vs face-centered-cubic) and compositional order (e.g., substitutional defects) of each individual particle. Deep autoencoders are then used to efficiently translate these neighborhood graphs into low-dimensional manifolds from which relationships among neighborhood graphs can be more easily inferred. We demonstrate the framework on in silico systems of DNA-functionalized particles, in which two well-recognized design parameters, particle size ratio and interparticle potential well depth can be adjusted independently. The framework reveals that binary colloidal mixtures with small interparticle size disparities (i.e., A- and B-type particle radius ratios of r _A/r _B = 0.8 to r _A/r _B = 0.95) can promote the self-assembly of defect-free BSLs much more effectively than systems of identically sized particles, as nearly defect-free BCC-CsCl, FCC-CuAu, and IrV crystals are observed in the former case. The framework additionally reveals that size-disparate colloidal mixtures can undergo nonclassical nucleation pathways where BSLs evolve from dense amorphous precursors, instead of directly nucleating from dilute solution. These findings illustrate that the presented characterization framework can assist in enhancing mechanistic understanding of the self-assembly of binary colloidal mixtures, which in turn can pave the way for engineering the growth of defect-free BSLs.

Shaping in the Third Direction; Fabrication of Hemispherical Micro-Concavity Array by Using Large Size Polystyrene Spheres as Template for Direct Self-Assembly of Small Size Silica Spheres

Silica and polystyrene spheres with a small size ratio (r = 0.005) form by sequential hanging drop self-assembly, a binary colloidal crystal through which calcination transforms in a silica-ordered concavity array. These arrays are capable of light Bragg diffraction and shape dependent optical phenomena, and they can be transformed into inverse-opal structures. Hierarchical 2D and 3D super-structures with ordered concavities as structural units were fabricated in this study.

Mild stratification in drying films of colloidal mixtures

Size stratification of bidisperse colloidal mixtures during vertical drying was investigated using the implicit solvent Langevin dynamics (LD) simulation and the explicit solvent lattice Boltzmann (LB) method. Simulations were performed for the Péclet number (Pe) over a wide range of 1-1000. In the case of a low size ratio of 2, mild stratification was observed in both simulation methods, in contrast to distinct stratification with thick "small-on-top" or "large-on-top" layers. The LD simulations exhibited a "small-on-top" stratification or mixed state. In contrast, the LB simulations exhibited a "large-on-top" or mixed state, according to the variation in Pe. The results demonstrated that the explicit solvent reduced the collective diffusion under moderate Pe conditions. This suppressed the steep concentration gradient of small particles in the packed region of particles near the air-solvent interface. Thus, distinguishable stratification patterns were obtained for the implicit and explicit solvent models.

Out-of-equilibrium processes in crystallization of organic-inorganic perovskites during spin coating

Complex phenomena are prevalent during the formation of materials, which affect their processing-structure-function relationships. Thin films of methylammonium lead iodide (CH₃NH₃PbI₃, MAPI) are processed by spin coating, antisolvent drop, and annealing of colloidal precursors. The structure and properties of transient and stable phases formed during the process are reported, and the mechanistic insights of the underlying transitions are revealed by combining in situ data from grazing-incidence wide-angle X-ray scattering and photoluminescence spectroscopy. Here, we report the detailed insights on the embryonic stages of organic-inorganic perovskite formation. The physicochemical evolution during the conversion proceeds in four steps: i) An instant nucleation of polydisperse MAPI nanocrystals on antisolvent drop, ii) the instantaneous partial conversion of metastable nanocrystals into orthorhombic solvent-complex by cluster coalescence, iii) the thermal decomposition (dissolution) of the stable solvent-complex into plumboiodide fragments upon evaporation of solvent from the complex and iv) the formation (recrystallization) of cubic MAPI crystals in thin film.

Complex phenomena are prevalent during the formation of materials, and they affect the processing structure-function relationship. Here the authors elucidate the stochastic transformation processes happening during the spin coating of perovskite colloidal precursors by multimodal characterization.

Deep learning for characterizing the self-assembly of three-dimensional colloidal systems

Creating a systematic framework to characterize the structural states of colloidal self-assembly systems is crucial for unraveling the fundamental understanding of these systems' stochastic and non-linear behavior. The most accurate characterization methods create high-dimensional neighborhood graphs that may not provide useful information about structures unless these are well-defined reference crystalline structures. Dimensionality reduction methods are thus required to translate the neighborhood graphs into a low-dimensional space that can be easily interpreted and used to characterize non-reference structures. We investigate a framework for colloidal system state characterization that employs deep learning methods to reduce the dimensionality of neighborhood graphs. The framework next uses agglomerative hierarchical clustering techniques to partition the low-dimensional space and assign physically meaningful classifications to the resulting partitions. We first demonstrate the proposed colloidal self-assembly state characterization framework on a three-dimensional in silico system of 500 multi-flavored colloids that self-assemble under isothermal conditions. We next investigate the generalizability of the characterization framework by applying the framework to several independent self-assembly trajectories, including a three-dimensional in silico system of 2052 colloidal particles that undergo evaporation-induced self-assembly.

Assembly of nanocrystal clusters by solvent evaporation: icosahedral order and the breakdown of the Maxwell regime

We carry out molecular dynamics simulations of N gold alkylthiolated nanocrystals (0 ≤ N ≤ 29) contained in liquid droplets of octane, nonane and decane coexisting with its vapor. The equilibrium structures that result when all the solvent dries up consist of highly symmetric nanocrystal clusters with different degrees of icosahedral order that are thoroughly characterized. We show that the relaxation times follow two regimes, a first for small nanocrystal packing fraction, dominated by the diffusion of vapor molecules (Maxwell regime, relaxation times independent of N) and another, for larger packing fractions, where the solvent diffuses through the cluster (with relaxation times growing like N2/3). We discuss the connection to the assembly of superlattices, prediction of lattice constants and evaporation models.

Stratification of polymer mixtures in drying droplets: Hydrodynamics and diffusion

We study the evaporation-induced stratification of a mixture of short and long polymer chains in a drying droplet using molecular simulations. We systematically investigate the effects of hydrodynamic interactions (HI) on this process by comparing hybrid simulations accounting for HI between polymers through the multiparticle collision dynamics technique with free-draining Langevin dynamics simulations neglecting the same. We find that the dried supraparticle morphologies are homogeneous when HI are included but are stratified in core-shell structures (with the short polymers forming the shell) when HI are neglected. The simulation methodology unambiguously attributes this difference to the treatment of the solvent in the two models. We rationalize the presence (or absence) of stratification by measuring phenomenological multicomponent diffusion coefficients for the polymer mixtures. The diffusion coefficients show the importance of not only solvent backflow but also HI between polymers in controlling the dried supraparticle morphology.

Mechanochromic and Thermochromic Sensors Based on Graphene Infused Polymer Opals

High quality opal-like photonic crystals containing graphene are fabricated using evaporation-driven self-assembly of soft polymer colloids. A miniscule amount of pristine graphene within a colloidal crystal lattice results in the formation of colloidal crystals with a strong angle-dependent structural color and a stop band that can be reversibly shifted across the visible spectrum. The crystals can be mechanically deformed or can reversibly change color as a function of their temperature, hence their sensitive mechanochromic and thermochromic response make them attractive candidates for a wide range of visual sensing applications. In particular, it is shown that the crystals are excellent candidates for visual strain sensors or integrated time-temperature indicators which act over large temperature windows. Given the versatility of these crystals, this method represents a simple, inexpensive, and scalable approach to produce multifunctional graphene infused synthetic opals and opens up exciting applications for novel solution-processable nanomaterial based photonics.

Temporal evolution of concentration and microstructure of colloidal films during vertical drying: a lattice Boltzmann simulation study

We study the temporal and structural development of colloid films during vertical drying using the lattice Boltzmann (LB) simulation. The dispersed particles moving in Brownian motion have excluded volume and hydrodynamic interactions in the film. The concentrated colloidal film formed by solvent evaporation is modeled as an uniaxial compression of colloids with a planar moving interface. The simulation studies are carried out over a wide range of Péclet number (Pe), the relative ratio between the evaporation rate and the diffusion rate of colloids. The results clearly demonstrate a temporal variation of colloid concentration as the evaporation rate increases. In the case of high Pe, the increase of colloid concentration in the top layer creates structural features that can be distinguished along the height of the film, and eventually can induce a large tensile stress in the layer. However, surprisingly, the colloids are maximally crystallized in the case of moderate Pe. The LB simulation results are further compared with those from previous studies of the Brownian Dynamics (BD) simulation and the continuum model for the evaporation film. The LB and BD results match well both at low and high Pe limits. The qualitatively significant differences between LB and BD simulations at a moderate Pe indicate that hydrodynamic interactions (HIs) play an important role in this Pe. The presence of HIs induces a greater reduction of diffusion than under geometrical restriction alone, and the effect is conspicuous when particles are driven both by diffusion and by advection.

Assembly by solvent evaporation: equilibrium structures and relaxation times

We present a study describing the dynamics and equilibrium of the assembly of nanostructures by solvent evaporation. We first consider N nanocrystals stabilized by capping ligands in a spherical droplet of liquid solvent coexisting with its gas and show that, as the liquid solvent evaporates slowly, NCs crystallize into clusters of high symmetry based on tetrahedral and octahedral units: tetrahedron (N = 4), octahedron (N = 6), icosahedron (N = 13), Archimedean truncated tetrahedron (N = 16) and Z₂₀ (N = 21). We derive explicit formulas for the process and rigorously compute relaxation times, which drastically increase when the packing parameter reaches the hard-sphere liquid-solid transition η = 0.49. This result shows that contrary to what occurs in an evaporation of a single component system, the relaxation times are not determined by the diffusion constant of the vapor, but rather, are dominated by the residence time of solvent molecules trapped within the capping ligands. Our theory provides a number of predictions that enable the design of new structures while improving the control and quality of their assembly.

Control of Stratification in Drying Particle Suspensions via Temperature Gradients

A potential strategy for controlling stratification in a drying suspension of bidisperse particles is studied using molecular dynamics simulations. When the suspension is maintained at a constant temperature during fast drying, it can exhibit "small-on-top" stratification with an accumulation (depletion) of smaller (larger) particles in the top region of the drying film, consistent with the prediction of current theories based on diffusiophoresis. However, when only the region near the substrate is thermalized at a constant temperature, a negative temperature gradient develops in the suspension because of evaporative cooling at the liquid-vapor interface. Since the associated thermophoresis is stronger for larger nanoparticles, a higher fraction of larger nanoparticles migrate to the top of the drying film at fast evaporation rates. As a result, stratification is converted to "large-on-top". Very strong small-on-top stratification can be produced with a positive thermal gradient in the drying suspension. Here, we explore a way to produce a positive thermal gradient by thermalizing the vapor at a temperature higher than that of the solvent. Possible experimental approaches to realize various thermal gradients in a suspension undergoing solvent evaporation and thus to produce different stratification states in the drying film are suggested.

Ordering of colloidal hard spheres under gravity: from monolayer to multilayer

The phase behaviour of hard spheres confined by a gravitational potential to a thin layer (up to several monolayers) near a hard, flat surface is investigated using grand canonical Monte Carlo simulation. Depending on the strength of the gravitational field, the bottom monolayer of spheres may adopt uniform hexagonal order before, during, or after the growth of the second layer of particles. The crossover from ordering with a sparsely populated overlayer to ordering with almost one-third of the system's particles forming a second layer is observed upon decreasing the dimensionless Péclet number Pe = mgσ/kBT from 18 to 16. The particular sensitivity of the nature of the transition to particle size in this range is interpreted in terms of competing influences on the base layer structure by particles in the overlayer: promotion of order through increased pressure, versus stabilization of defects through occupation of low-lying sites on top of them. Simulations of grain boundaries between 2-D ordered domains of different orientation are used to correlate the degree of overlayer coverage to its effects on grain boundary stiffness as an indicator of defect free energy. Finally, we examine the structure of the ordered phases at coexistence over a range of gravitational strengths and find that orientational ordering of the second monolayer occurs along with first-order transition of the base layer at Pe = 8 but not at Pe = 10.