Atomistic Mismatch Defines Energy-Structure Relationships during Oriented Attachment of Nanoparticles

Ion Correlations Decrease Particle Aggregation Rate by Increasing Hydration Forces at Interfaces

The connection between solution structure, particle forces, and emergent phenomena at solid-liquid interfaces remains ambiguous. In this case study on boehmite aggregation, we established a connection between interfacial solution structure, emerging hydration forces between two approaching particles, and the resulting structure and kinetics of particle aggregation. In contrast to expectations from continuum-based theories, we observed a nonmonotonic dependence of the aggregation rate on the concentration of sodium chloride, nitrate, or nitrite, decreasing by 15-fold in 4 molal compared to 1 molal solutions. These results are accompanied by an increase in repulsive hydration forces and interfacial oscillatory features from 0.27-0.31 nm in 0.01 molal to 0.38-0.52 nm in 2 molal. Moreover, molecular dynamics (MD) simulations indicated that these changes correspond to enhanced ion correlations near the interface and produced loosely bound aggregates that retain electrolyte between the particles. We anticipate that these results will enable the prediction of particle aggregation, attachment, and assembly, with broad relevance to interfacial phenomena.

Resolving intermediates during the growth of aluminum deuteroxide (Hydroxide) polymorphs in high chemical potential solutions

Aluminum hydroxide polymorphs are of widespread importance yet their kinetics of nucleation and growth remain beyond the reach of current models. Here we attempt to unveil the reaction processes underlying the polymorphs formation at high chemical potential. We examine their formation in-situ from supersaturated alkaline sodium aluminate solutions using deuteration and time-resolved neutron pair distribution function analyses, which indicate the formation of individual Al(OD)3 layers as an intermediate particle phase. These layers ultimately stack to form gibbsite- or bayerite-like layered heterostructures. Ex-situ characterization of the recovered precipitates using 27Al magic angle spinning nuclear magnetic resonance spectroscopy, Raman, X-ray diffraction, and scanning electron microscopy, suggests the presence of additional intermediate states, an amorphous compound bearing both tetrahededrally- and penta-coordinated Al3+. These observations reveal the complex pathways to form Al(OD)3 monolayers via either transient pentacoordinate species or amorphous-to-ordered transitions. The subsequent crystallization of admixed gibbsite/bayerite is followed by an Al(OD)3 monolayer attachment process.

Role of Solvents in Oriented Attachment of Ag Nanoparticles: Insights from Molecular Dynamics Simulations and Topological Analysis

Although the solvation force is considered one of the key forces behind the oriented attachment (OA), the precise roles of solvents in this process remain incompletely elucidated. In this study, we examined the effect of solvent polarities (water, acetone, and chloroform) on the attachment of silver nanoparticles by calculating the free energy curves for the OA process. The observed magnitudes of the binding energies and approaching and dissociation energy barriers are commensurate with the respective solvent polarities. Consequently, OA is more likely to occur in acetone with an intermediate permittivity relative to that of water and chloroform. Additionally, we identified a topological descriptor, namely, the Euler characteristic, of the solvent network, especially the water network, between two approaching surfaces, which manifests a linear correlation with the observed free energy profiles. This descriptor holds promise as a quantitative tool for predicting interactions between nanoparticles in solvent environments featuring hydrogen bond networks.

Energetics of water expulsion from intervening space between two particles during aggregation

Solvent expulsion away from an intervening region between two approaching particles plays important roles in particle aggregation yet remains poorly understood. In this work, we use metadynamics molecular simulations to study the free energy landscape of removing water molecules from gibbsite and pyrophyllite slit pores representing the confined spaces between two approaching particles. For gibbsite, removing water from the intervening region is both entropically and enthalpically unfavorable. The closer the particles approach each other, the harder it is to expel water molecules. For pyrophyllite, water expulsion is spontaneous, which is different from the gibbsite system. A smaller pore makes the water removal more favorable. When water is being drained from the intervening region, single chains of water molecules are observed in gibbsite pore, while in pyrophyllite pore water cluster is usually observed. Water-gibbsite hydrogen bonds help stabilize water chains, while water forms clusters in pyrophyllite pore to maximize the number of hydrogen bonds among themselves. This work provides the first assessment into the energetics and structure of water being drained from the intervening region between two approaching particles during oriented attachment and aggregation.

Molecular Mechanisms of Sorbed Ion Effects during Boehmite Particle Aggregation

Classical theories of particle aggregation, such as Derjaguin-Landau-Verwey-Overbeek (DLVO), do not explain recent observations of ion-specific effects or the complex concentration dependence for aggregation. Thus, here, we probe the molecular mechanisms by which selected alkali nitrate ions (Na+, K+, and NO3-) influence aggregation of the mineral boehmite (γ-AlOOH) nanoparticles. Nanoparticle aggregation was analyzed using classical molecular dynamics (CMD) simulations coupled with the metadynamics rare event approach for stoichiometric surface terminations of two boehmite crystal faces. Calculated free energy landscapes reveal how electrolyte ions alter aggregation on different crystal faces relative to pure water. Consistent with experimental observations, we find that adding an electrolyte significantly reduces the energy barrier for particle aggregation (∼3-4×). However, in this work, we show this is due to the ions disrupting interstitial water networks, and that aggregation between stoichiometric (010) basal-basal surfaces is more favorable than between (001) edge-edge surfaces (∼5-6×) due to the higher interfacial water densities on edge surfaces. The interfacial distances in the interlayer between aggregated particles with electrolytes (∼5-10 Å) are larger than those in pure water (a few Ångströms). Together, aggregation/disaggregation in salt solutions is predicted to be more reversible due to these lower energy barriers, but there is uncertainty on the magnitudes of the energies that lead to aggregation at the molecular scale. By analyzing the peak water densities of the first monolayer of interstitial water as a proxy for solvent ordering, we find that the extent of solvent ordering likely determines the structures of aggregated states as well as the energy barriers to move between them. The results suggest a path for developing a molecular-level basis to predict the synergies between ions and crystal faces that facilitate aggregation under given solution conditions. Such fundamental understanding could be applied extensively to the aggregation and precipitation utilization in the biological, pharmaceutical, materials design, environmental remediation, and geological regimes.

Kinetics of Oriented Attachment of Mica Crystals

Oriented attachment (OA), that is, the coalescence of crystals through attachment on coaligned crystal faces, is a nonclassical crystal growth process. Before attachment, a mesocrystal consisting of coaligned parallel crystals but with liquid separating them was observed. Fundamental questions such as why OA is kinetically favored and whether a mesocrystal stage is a prerequisite for OA are raised. Through combining brute-force molecular dynamics simulations and path samplings based on extensive umbrella simulations, we address these questions with a case study on the OA of a mica nanocrystal onto a mica crystal substrate in water. Brute-force simulations show that if two mica crystals are attached but largely misaligned, coalignment hardly appears. Thus, if OA is possible, then coalignment must appear before the attachment between crystals. Electrophoresis of the nanocrystal toward the substrate surface is spontaneous, but mesocrystal formation is occasional, also shown by brute-force simulations. Free energies along different pathways show that OA is spontaneous and kinetically favored over non-OA, and a mesocrystal formation is just a bifurcation in the pathway. OA is through a pathway in which the nanocrystal is tilted with respect to the substrate. Part of the nanocrystal is attached to the substrate first, and then, OA is gradually completed. Once a mesocrystal is occasionally formed, then a jump event is needed for the nanocrystal to get back to the OA pathway. The sampling technique here can hopefully guide the design of nanostructured materials facilitated by OA.

Effects of Size and Shape on the Tolerances for Misalignment and Probabilities for Successful Oriented Attachment of Nanoparticles

Oriented attachment (OA) of nanoparticles is an important pathway of crystal growth, but there is a lack of tools to model OA. Here, we present several simple models that relate the probability of achieving OA to basic geometric parameters, such as particle size, shape, and lattice periodicity. A Moiré-domain model is applied to understand twist misorientations between parallel surfaces, and it predicts that the range of twist angles yielding perfect OA is inversely related to the width of the contact area. This idea is explored further through a surface functional model, which investigates how patterns of crystallographic registration can drive the emergence of complex orientational energy landscapes. The energy landscapes are predicted to possess multiple local minima that can trap particles in imperfect alignments, and these local minima become deeper and more numerous as the contact area increases, which makes OA more challenging for large particles. A second set of models is presented to understand the sequence of events by which two crystallographic faces become coplanar after the collision. We use a central force approximation to predict the odds that two particle faces will attain coalignment when the particles collide with random misalignments, and we show that in the absence of special biasing forces, the probability of attaining alignment on a given face is roughly proportional to its solid angle as viewed from the center of the particle. The model thus predicts that OA is most favorable between well-faceted particles and becomes exceedingly unlikely for large spherical particles that express many microfacets.