Laura Filion

Entropy-Driven Formation of Binary Semiconductor-Nanocrystal Superlattices

One of the main reasons for the current interest in colloidal nanocrystals is their propensity to form superlattices, systems in which (different) nanocrystals are in close contact in a well-ordered three-dimensional (3D) geometry resulting in novel material properties. However, the principles underlying the formation of binary nanocrystal superlattices are not well understood. Here, we present a study of the driving forces for the formation of binary nanocrystal superlattices by comparing the formed structures with full free energy calculations. The nature (metallic or semiconducting) and the size-ratio of the two nanocrystals are varied systematically. With semiconductor nanocrystals, self-organization at high temperature leads to superlattices (AlB(2), NaZn(13), MgZn(2)) in accordance with the phase diagrams for binary hard-sphere mixtures; hence entropy increase is the dominant driving force. A slight change of the conditions results in structures that are energetically stabilized. This study provides rules for the rational design of 3D nanostructured binary semiconductors, materials with promises in thermoelectrics and photovoltaics and which cannot be reached by any other technology.

Crystal nucleation of hard spheres using molecular dynamics, umbrella sampling, and forward flux sampling: a comparison of simulation techniques.

Over the last number of years several simulation methods have been introduced to study rare events such as nucleation. In this paper we examine the crystal nucleation rate of hard spheres using three such numerical techniques: molecular dynamics, forward flux sampling, and a Bennett-Chandler-type theory where the nucleation barrier is determined using umbrella sampling simulations. The resulting nucleation rates are compared with the experimental rates of Harland and van Megen [Phys. Rev. E 55, 3054 (1997)], Sinn et al. [Prog. Colloid Polym. Sci. 118, 266 (2001)], Schätzel and Ackerson [Phys. Rev. E 48, 3766 (1993)], and the predicted rates for monodisperse and 5% polydisperse hard spheres of Auer and Frenkel [Nature 409, 1020 (2001)]. When the rates are examined in units of the long-time diffusion coefficient, we find agreement between all the theoretically predicted nucleation rates, however, the experimental results display a markedly different behavior for low supersaturation. Additionally, we examined the precritical nuclei arising in the molecular dynamics, forward flux sampling, and umbrella sampling simulations. The structure of the nuclei appears independent of the simulation method, and in all cases, the nuclei contains on average significantly more face-centered-cubic ordered particles than hexagonal-close-packed ordered particles.

Observation of a Ternary Nanocrystal Superlattice and Its Structural Characterization by Electron Tomography

The first genuine ternary colloidal crystal (see picture) is composed of PbSe nanocrystals of two different diameters (blue and green) and of CdSe nanocrystals (red). Electron tomography shows that the superlattice is isostructural with the atomic lattice AlMgB4.

Entropy-driven formation of large icosahedral colloidal clusters by spherical confinement

Icosahedral symmetry, which is not compatible with truly long-range order, can be found in many systems, such as liquids, glasses, atomic clusters, quasicrystals and virus-capsids. To obtain arrangements with a high degree of icosahedral order from tens of particles or more, interparticle attractive interactions are considered to be essential. Here, we report that entropy and spherical confinement suffice for the formation of icosahedral clusters consisting of up to 100,000 particles. Specifically, by using real-space measurements on nanometre- and micrometre-sized colloids, as well as computer simulations, we show that tens of thousands of hard spheres compressed under spherical confinement spontaneously crystallize into icosahedral clusters that are entropically favoured over the bulk face-centred cubic crystal structure. Our findings provide insights into the interplay between confinement and crystallization and into how these are connected to the formation of icosahedral structures.

Erasing no-man’s land by thermodynamically stabilizing the liquid–liquid transition in tetrahedral particles

One of the most controversial hypotheses for explaining the origin of the thermodynamic anomalies characterizing liquid water postulates the presence of a metastable second-order liquid-liquid critical point [1] located in the “no-man’s land” [2]. In this scenario, two liquids with distinct local structure emerge near the critical temperature. Unfortunately, since spontaneous crystallization is rapid in this region, experimental support for this hypothesis relies on significant extrapolations, either from the metastable liquid or from amorphous solid water [3, 4]. Although the liquid-liquid transition is expected to feature in many tetrahedrally coordinated liquids, including silicon [5], carbon [6] and silica, even numerical studies of atomic and molecular models have been unable to conclusively prove the existence of this transition. Here we provide such evidence for a model in which it is possible to continuously tune the softness of the interparticle interaction and the flexibility of the bonds, the key ingredients controlling the existence of the critical point. We show that conditions exist where the full coexistence is thermodynamically stable with respect to crystallization. Our work offers a basis for designing colloidal analogues of water exhibiting liquid-liquid transitions in equilibrium, opening the way for experimental confirmation of the original hypothesis.

Fabrication of large binary colloidal crystals with a NaCl structure

Binary colloidal crystals offer great potential for tuning material properties for applications in, for example, photonics, semiconductors and spintronics, because they allow the positioning of particles with quite different characteristics on one lattice. For micrometer-sized colloids, it is believed that gravity and slow crystallization rates hinder the formation of high-quality binary crystals. Here, we present methods for growing binary colloidal crystals with a NaCl structure from relatively heavy, hard-sphere-like, micrometer-sized silica particles by exploring the following external fields: electric, gravitational, and dielectrophoretic fields and a structured surface (colloidal epitaxy). Our simulations show that the free-energy difference between the NaCl and NiAs structures, which differ in their stacking of the hexagonal planes of the larger spheres, is very small (≈0.002 kBT). However, we demonstrate that the fcc stacking of the large spheres, which is crucial for obtaining the pure NaCl structure, can be favored by using a combination of the above-mentioned external fields. In this way, we have successfully fabricated large, 3D, oriented single crystals having a NaCl structure without stacking disorder.

Vacancy-stabilized crystalline order in hard cubes

We examine the effect of vacancies on the phase behavior and structure of systems consisting of hard cubes using event-driven molecular dynamics and Monte Carlo simulations. We find a first-order phase transition between a fluid and a simple cubic crystal phase that is stabilized by a surprisingly large number of vacancies, reaching a net vacancy concentration of approximately 6.4% near bulk coexistence. Remarkably, we find that vacancies increase the positional order in the system. Finally, we show that the vacancies are delocalized and therefore hard to detect.

Simulation of nucleation in almost hard-sphere colloids: The discrepancy between experiment and simulation persists

In this paper we examine the phase behavior of the Weeks-Chandler-Andersen (WCA) potential with βε = 40. Crystal nucleation in this model system was recently studied by Kawasaki and Tanaka [Proc. Natl. Acad. Sci. U.S.A. 107, 14036 (2010)], who argued that the computed nucleation rates agree well with experiment, a finding that contradicted earlier simulation results. Here we report an extensive numerical study of crystallization in the WCA model, using three totally different techniques (Brownian dynamics, umbrella sampling, and forward flux sampling). We find that all simulations yield essentially the same nucleation rates. However, these rates differ significantly from the values reported by Kawasaki and Tanaka and hence we argue that the huge discrepancy in nucleation rates between simulation and experiment persists. When we map the WCA model onto a hard-sphere system, we find good agreement between the present simulation results and those that had been obtained for hard spheres [L. Filion, M. Hermes, R. Ni, and M. Dijkstra, J. Chem. Phys. 133, 244115 (2010); S. Auer and D. Frenkel, Nature 409, 1020 (2001)].

Predicting patchy particle crystals: variable box shape simulations and evolutionary algorithms

We consider several patchy particle models that have been proposed in literature and we investigate their candidate crystal structures in a systematic way. We compare two different algorithms for predicting crystal structures: (i) an approach based on Monte Carlo simulations in the isobaric-isothermal ensemble and (ii) an optimization technique based on ideas of evolutionary algorithms. We show that the two methods are equally successful and provide consistent results on crystalline phases of patchy particle systems.

Stability of LS and LS2 crystal structures in binary mixtures of hard and charged spheres

We study by computer simulations the stability of various crystal structures in a binary mixture of large and small spheres interacting either with a hard sphere or a screened-Coulomb potential. In the case of hard-core systems, we consider structures that have atomic prototypes CrB, gammaCuTi, alphaIrV, HgBr2, AuTe2, Ag2Se and the Laves phases (MgCu2, MgNi2, and MgZn2) as well as a structure with space group symmetry 74. By utilizing Monte Carlo simulations to calculate Gibbs free energies, we determine composition versus pressure and constant volume phase diagrams for diameter ratios of q=0.74, 0.76, 0.8, 0.82, 0.84, and 0.85 for the small and large spheres. For diameter ratios 0.76 < or = q < or = 0.84, we find the Laves phases to be stable with respect to the other crystal structures that we considered and the fluid mixture. By extrapolating to the thermodynamic limit, we show that the MgZn2 structure is the most stable one of the Laves structures. We also calculate phase diagrams for equally and oppositely charged spheres for size ratio of 0.73 taking into consideration the Laves phases and CsCl. In the case of equally charged spheres, we find a pocket of stable Laves phases, while in the case of oppositely charged spheres, Laves phases are found to be metastable with respect to the CsCl and fluid phases.