Assembly by solvent evaporation: equilibrium structures and relaxation times

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.

High-Precision Calculation of Nanoparticle (Nanocrystal) Potentials of Mean Force and Internal Energies

Thermodynamic stability assessment of nanocrystal systems requires precise free energy calculations. This study highlights the importance of meticulous control over various factors, including the thermostat, time step, potential cutoff, initial configuration, sampling method, and overall simulation duration. Free energy computations in dry (solvent-free) systems are on the order of several hundred kBT but can be obtained with consistent accuracy. However, calculation of internal energies becomes challenging, as they are typically much larger in magnitude than free energies and exhibit significant noise and reduced reliability. To address this limitation, we propose a new internal energy estimate that drastically reduces the noise. We also present formulas that enable the optimization of the parameters of the harmonic bias potential for optimal convergence. Finally, we discuss the implications of these findings for the computation of free energies in nanocrystal clusters and superlattices.

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.

Ligand Effects in Assembly of Cubic and Spherical Nanocrystals: Applications to Packing of Perovskite Nanocubes

We establish the formula representing cubic nanocrystals (NCs) as hard cubes taking into account the role of the ligands and describe how these results generalize to any other NC shapes. We derive the conditions under which the hard cube representation breaks down and provide explicit expressions for the effective size. We verify the results from the detailed potential of mean force calculations for two nanocubes in different orientations as well as with spherical nanocrystals. Our results explicitly demonstrate the relevance of certain ligand conformations, i.e., "vortices", and show that edges and corners provide natural sites for their emergence. We also provide both simulations and experimental results with single component cubic perovskite nanocrystals assembled into simple cubic superlattices, which further corroborate theoretical predictions. In this way, we extend the Orbifold Topological Model (OTM) accounting for the role of ligands beyond spherical nanocrystals and discuss its extension to arbitrary nanocrystal shapes. Our results provide detailed predictions for recent superlattices of perovskite nanocubes and spherical nanocrystals. Problems with existing united atom force fields are discussed.

Body centered tetragonal nanoparticle superlattices: why and when they form?

Body centered tetragonal (BCT) phases are structural intermediates between body centered cubic (BCC) and face centered cubic (FCC) structures. However, BCC ↔ FCC transitions may or may not involve a stable BCT intermediate. Interestingly, nanoparticle superlattices usually crystallize in BCT structures, but this phase is much less frequent for colloidal crystals of micrometer-sized particles. Two origins have been proposed for the formation of BCT NPSLs: (i) the influence of the substrate on which the nanoparticle superlattice is deposited, and (ii) non-spherical nanoparticle shapes, combined with the fact that different crystal facets have different ligand organizations. Notably, none of these two mechanisms alone is able to explain the set of available experimental observations. In this work, these two hypotheses were independently tested using a recently developed molecular theory for nanoparticle superlattices that explicitly captures the degrees of freedom associated with the ligands on the nanoparticle surface and the crystallization solvent. We show that the presence of a substrate can stabilize the BCT structure for spherical nanoparticles, but only for very specific combinations of parameters. On the other hand, a truncated-octahedron nanoparticle shape strongly stabilizes BCT structures in a wide region of the phase diagram. In the latter case, we show that the stabilization of BCT results from the geometry of the system and it does not require different crystal facets to have different ligand properties, as previously proposed. These results shed light on the mechanisms of BCT stabilization in nanoparticle superlattices and provide guidelines to control its formation.

Ligand structure and adsorption free energy of nanocrystals on solid substrates

We present an investigation on the absorption of alkylthiolated nanocrystals on a solid substrate. We calculate adsorption free energies and report a number of effects induced by the substrate. Nearest neighbor distances and bonding free energies are significantly different than for a free floating case, there is a weakening of bonding free energies among nanocrystals, and the adsorption is manifestly anisotropic, i.e., stronger along certain directions of the nanocrystal core. We contend that this last result accounts for the Bain transition (fcc → bcc) observed in experimental results. We report the presence of vortices induced by the substrate, which explain the increased nearest neighbor distance among nanocrystals, which is in excellent quantitative agreement with experimental results and with the predictions of the Orbifold Topological Model. Implications for the assembly of nanostructures and future experiments are also discussed.

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.