Understanding the mechanisms of anisotropic dissolution in metal oxides by applying radiolysis simulations to liquid-phase TEM.

Iron-based redox-active minerals are ubiquitous in soils, sediments, and aquatic systems. Their dissolution is of great importance for microbial impacts on carbon cycling and the biogeochemistry of the lithosphere and hydrosphere. Despite its widespread significance and extensive prior study, the atomic-to-nanoscale mechanisms of dissolution remain poorly understood, particularly the interplay between acidic and reductive processes. Here, we use in situ liquid-phase-transmission electron microscopy (LP-TEM) and simulations of radiolysis to probe and control acidic versus reductive dissolution of akaganeite (ß-FeOOH) nanorods. Informed by crystal structure and surface chemistry, the balance between acidic dissolution at rod tips and reductive dissolution at rod sides was systematically varied using pH buffers, background chloride anions, and electron beam dose. We find that buffers, such as bis-tris, effectively inhibited dissolution by consuming radiolytic acidic and reducing species such as superoxides and aqueous electrons. In contrast, chloride anions simultaneously suppressed dissolution at rod tips by stabilizing structural elements while promoting dissolution at rod sides through surface complexation. Dissolution behaviors were systematically varied by shifting the balance between acidic and reductive attacks. The findings show LP-TEM combined with simulations of radiolysis effects can provide a unique and versatile platform for quantitatively investigating dissolution mechanisms, with implications for understanding metal cycling in natural environments and the development of tailored nanomaterials.

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.

In situ imaging of amorphous intermediates during brucite carbonation in supercritical CO2.

Progress in understanding crystallization pathways depends on the ability to unravel relationships between intermediates and final crystalline products at the nanoscale, which is a particular challenge at elevated pressure and temperature. Here we exploit a high-pressure atomic force microscope to directly visualize brucite carbonation in water-bearing supercritical carbon dioxide (scCO2) at 90 bar and 50 °C. On introduction of water-saturated scCO2, in situ visualization revealed initial dissolution followed by nanoparticle nucleation consistent with amorphous magnesium carbonate (AMC) on the surface. This is followed by growth of nesquehonite (MgCO3·3H2O) crystallites. In situ imaging provided direct evidence that the AMC intermediate acts as a seed for crystallization of nesquehonite. In situ infrared and thermogravimetric-mass spectrometry indicate that the stoichiometry of AMC is MgCO3·xH2O (x = 0.5-1.0), while its structure is indicated to be hydromagnesite-like according to density functional theory and X-ray pair distribution function analysis. Our findings thus provide insight for understanding the stability, lifetime and role of amorphous intermediates in natural and synthetic systems.

Spiers Memorial Lecture: Assembly-based pathways of crystallization.

Solution crystallization of materials ranging from simple salts to complex supramolecular assemblies has long been viewed through the lens of classical nucleation and growth theories in which monomeric building blocks assemble into ordered structures through inherent thermal fluctuations that overcome a free energy barrier and continue to grow by the addition of such units to atomic steps. However, recent observations have revealed a rich set of hierarchical pathways during both nucleation and growth involving species of a higher order than monomers. While many studies have investigated and deduced the mechanisms underlying hierarchical nucleation pathways, much less research has been directed towards the development of a mechanistic picture of growth by the assembly of more complex units. Here, we review recent investigations into crystal growth by particle attachment, with an emphasis on oriented attachment. We discuss the relationship between interfacial structure, interparticle forces, and attachment dynamics, discuss the consequences of size dependent phase stability, and examine the impact of the ligand-functionalization of primary particles.

Interplay between Short- and Long-Ranged Forces Leading to the Formation of Ag Nanoparticle Superlattice.

Nanoparticle (NP) superlattices have attracted increasing attention due to their unique physicochemical properties. However, key questions persist regarding the correlation between short- and long-range driving forces for nanoparticle assembly and resultant capability to predict the transient and final superlattice structure. Here the self-assembly of Ag NPs in aqueous solutions is investigated by employing in situ liquid cell transmission electron microscopy, combined with atomic force microscopy-based force measurements, and theoretical calculations. Despite the NPs exhibiting instantaneous Brownian motion, it is found that the dynamic behavior of NPs is correlated with the van der Waals force, sometimes unexpectedly over relatively large particle separations. After the NPs assemble into clusters, a delicate balance between the hydration and van der Waals forces results in a distinct distribution of particle separation, which is ascribed to layers of hydrated ions adsorbed on the NP surface. The study demonstrates pivotal roles of the complicated correlation between interparticle forces; potentially enabling the control of particle separation, which is critical for tailoring the properties of NP superlattices.

Polyelectrolyte complex films influence the formation of polycrystalline micro-structures.

Silica-carbonate biomorphs are inorganic materials composed of thousands of crystalline nanorods that assemble complex morphologies such as helices, vessels, and sheets. We investigate the effect on biomorph crystallization of polyelectrolyte complex films that are prepared using the layer-by-layer deposition technique and post-processed to obtain three stable, chemically distinct films. Biomorph growth on poly(diallyldimethylammonium)-dominated substrates (cationic) shows polycrystalline helical and sheet structures bounded by large witherite prisms. Crystallization on poly(styrenesulfonate)-dominated (anionic) and stoichiometric substrates follows a qualitatively different pathway. We observe islands of radial mineral films that over several days extend at a remarkably constant velocity of 0.48 µm h-1 and eventually mineralize the whole substrate. Our work opens exciting avenues for the use of polyelectrolyte films as tunable substrates for biomimetic crystallization.

Systematic characterization of polycrystalline silica-carbonate helices.

Biomorphs are complex, life-like structures that emerge from the precipitation of barium carbonate and amorphous silica in alkaline media. Despite their inorganic nature, these microstructures have non-crystallographic morphologies such as helices and cardioid sheets. At the nanoscale, biomorphs arrange thousands of crystalline nanorods as hierarchical assemblies that resemble natural biominerals suggesting novel approaches towards the production of biomimetic materials. We report the synthesis of silica-carbonate biomorphs in single-phase, gradient-free solutions that differ markedly from the typical solution-gas or gel-solution setups. Our experimental approach significantly increases the duration of biomorph growth and hence assembles networks in which individual helices extend to several millimeters. These unusually long biomorphs allow the first quantitative measurements of mesoscopic parameters such as the helix wavelength, period, width, and linear as well as tangential growth velocities. We find that the latter quantities are system-specific and tightly conserved during many hours of growth. Moreover, the average double helix wavelength of (19 ± 3) µm and width of (9.6 ± 0.8) µm vary by less than 12% when the initial carbonate concentration is changed by three orders of magnitude. We also delineate the single helix growth mechanism and report the occurrence of ribbon-like structures and highly regular "superhelices". Our experiments clearly demonstrate the robustness and consistency of biomorph growth under stable chemical conditions.

Do dissolving objects converge to a universal shape?

Surprisingly, macroscopic objects such as melting ice cubes and growing stalactites approach nonintuitive geometric ideals. Here we investigate the shape of dissolving cylinders in a large volume of water. The cylinders are oriented vertically and consist of amorphous glucose or poly(ethylene glycol). The dissolution causes density differences in the surrounding fluid, which induce gravity-driven convection downward along the object. The resulting concentration gradient shapes the cylinder according to fast dissolution at the tip and slow dissolution at the base. The contour of the object approaches a power law of the form z ∝ R(2), where z is the vertical distance from the tip and R is the corresponding radius. We suggest that this paraboloidal shape is the geometric attractor for the dissolution of noncrystalline objects in the presence of gravity.

Effect of Solvent Composition on Non-DLVO Forces and Oriented Attachment of Zinc Oxide Nanoparticles.

Oriented attachment (OA) occurs when nanoparticles in solution align their crystallographic axes prior to colliding and subsequently fuse into single crystals. Traditional colloidal theories such as DLVO provide a framework for evaluating OA but fail to capture key particle interactions due to the atomistic details of both the crystal structure and the interfacial solution structure. Using zinc oxide as a model system, we investigated the effect of the solvent on short-ranged and long-ranged particle interactions and the resulting OA mechanism. In situ TEM imaging showed that ZnO nanocrystals in toluene undergo long-range attraction comparable to 1kT at separations of 10 nm and 3kT near particle contact. These observations were rationalized by considering non-DLVO interactions, namely, dipole-dipole forces and torques between the polar ZnO nanocrystals. Langevin dynamics simulations showed stronger interactions in toluene compared to methanol solvents, consistent with the experimental results. Concurrently, we performed atomic force microscopy measurements using ZnO-coated probes for the short-ranged interaction. Our data are relevant to another type of non-DLVO interaction, namely, the repulsive solvation force. Specifically, the solvation force was stronger in water compared to ethanol and methanol, due to the stronger hydrogen bonding and denser packing of water molecules at the interface. Our results highlight the importance of non-DLVO forces in a general framework for understanding and predicting particle aggregation and attachment.

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.

Effects of particle shape and surface roughness on van der Waals interactions and coupling to dynamics in nanocrystals.

The van der Waals interaction between colloids and nanoparticles is one of the key components to understanding particle aggregation, attachment, and assembly. While the ubiquity of anisotropic particle shapes and surface roughness is well-recognized in nanocrystalline materials, the effects of both on van der Waals forces and torques have not been adequately investigated. In this study, we develop a numerical scheme to determine the van der Waals forces and torques between cubic particles with multiple configurations and relative orientations. Our results show that the van der Waals torque due to anisotropic particle shapes is appreciable at nearly all configurations and mutual angles, outcompeting Brownian torque for various materials systems and conditions. Surface roughness enhances this particle shape effect, resulting in stronger van der Waals interactions ascribed to protrusions on the surfaces. Moreover, a scaling analysis indicates that the surface roughness alters the separation dependence of the van der Waals force and, more importantly, significantly influences the dynamics of two approaching particles. Our results clearly demonstrate that surface roughness and anisotropic shape play a crucial role in the energetics and kinetics of various particle-scale and emergent phenomena, such as crystal growth by oriented attachment, nanomaterials synthesis and assembly, mud flow rheology, as well as the deposition of natural nanocrystals within the subsurface.

Discrete Orientations of Interfacial Waters Direct Crystallization of Mica-Binding Proteins.

Understanding the basis of templated molecular assembly on a solid surface requires a fundamental comprehension of both short- and long-range aqueous response to the surface under a variety of solution conditions. Herein we provide a detailed picture of how the molecular-scale response to different mica surfaces yields distinct solvent orientations that produce quasi-static directional potentials onto which macromolecules can adsorb. We connect this directionality to observed (a)symmetric epitaxial alignment of designed proteins onto these surfaces, corroborate our findings with 3D atomic force microscopy experiments, and identify slight differences in surface structure as the origin of this effect. Our work provides a detailed picture of the intrinsic electrolyte response in the vicinity of mineral interfaces, with clear predictions for experiment, and highlights the role of solvent on the predictive assembly of hierarchical materials on mineral surfaces.

Predicting Outcomes of Nanoparticle Attachment by Connecting Atomistic, Interfacial, Particle, and Aggregate Scales.

Predicting nanoparticle aggregation and attachment phenomena requires a rigorous understanding of the interplay among crystal structure, particle morphology, surface chemistry, solution conditions, and interparticle forces, yet no comprehensive picture exists. We used an integrated suite of experimental, theoretical, and simulation methods to resolve the effect of solution pH on the aggregation of boehmite nanoplatelets, a case study with important implications for the environmental management of legacy nuclear waste. Real-time observations showed that the particles attach preferentially along the (010) planes at pH 8.5 and the (101) planes at pH 11. To rationalize these results, we established the connection between key physicochemical phenomena across the relevant length scales. Starting from molecular-scale simulations of surface hydroxyl reactivity, we developed an interfacial-scale model of the corresponding electrostatic potentials, with subsequent particle-scale calculations of the resulting driving forces allowing successful prediction of the attachment modes. Finally, we scaled these phenomena to understand the collective structure at the aggregate-scale. Our results indicate that facet-specific differences in surface chemistry produce heterogeneous surface charge distributions that are coupled to particle anisotropy and shape-dependent hydrodynamic forces, to play a key role in controlling aggregation behavior.

Fractal structures in two-metal electrodeposition systems II: Cu and Zn.

In this second part of our study on fractal co-electrochemical deposition, we investigate the Cu-Zn system. Macroscopic and microscopic inspection shows a sensitive dependence of the morphology of the final pattern on initial concentrations. The pattern is seen to undergo a transition from classical dendrites to randomly ramified deposits, with each slight increase in [Cu(2+)](0), while [Zn(2+)](0) is maintained constant. The variational trends in chemical composition, growth velocity, and fractal dimension with increasing [Cu(2+)](0) are analyzed. The latter is seen to generally increase with copper (II) ion concentration. In contrast, the growth rate of the deposits is seen to decrease with increasing concentration of Cu(2+) ions. A new probe of dense ramified morphology, the pattern density, is introduced and seen to increase with [Cu(2+)](0). XRD measurements reveal that the observed properties correlate with the birth of copper-rich nuclei, which disrupt the crystalline anisotropy of the two-metal alloy.

Fractal structures in two-metal electrodeposition systems I: Pb and Zn.

Pattern formation in two-metal electrochemical deposition has been scarcely explored in the chemical literature. In this paper, we report new experiments on zinc-lead fractal co-deposition. Electrodeposits are grown in special cells at a fixed large value of the zinc ion concentration, while that of the lead ion is increased gradually. A very wide diversity of morphologies are obtained and classified. Most of the deposited domains are almost exclusively Pb or Zn. But certain regions originating at the base cathode, ranging from a short grass alley to dense, grown-up bushes or shrubs, manifest a combined Pb-Zn composition. Composition is determined using scanning electron microscopy/energy dispersive x ray measurements as well atomic absorption spectroscopy. Pb domains are characterized by shiny leaf-like and dense deposits as well as flowers with round, balloon-like corollas. The Zn zones display a greater variety of morphologies such as thick trunks and thin and fine branching, in addition to minute "cigar flower" structures. The various morphologies are analyzed and classified from the viewpoint of fractal nature, characterized by the box-count fractal dimension. Finally, macroscopic spatial alternation between two different characteristic morphologies is observed under certain conditions.

Visualizing Solution Structure at Solid-Liquid Interfaces using Three-Dimensional Fast Force Mapping.

Amongst the challenges for a variety of research fields are the visualization of solid-liquid interfaces and understanding how they are affected by the solution conditions such as ion concentrations, pH, ligands, and trace additives, as well as the underlying crystallography and chemistry. In this context, three-dimensional fast force mapping (3D FFM) has emerged as a promising tool for investigating solution structure at interfaces. This capability is based on atomic force microscopy (AFM) and allows the direct visualization of interfacial regions in three spatial dimensions with sub-nanometer resolution. Here we provide a detailed description of the experimental protocol for acquiring 3D FFM data. The main considerations for optimizing the operating parameters depending on the sample and application are discussed. Moreover, the basic methods for data processing and analysis are discussed, including the transformation of the measured instrument observables into tip-sample force maps that can be linked to the local solution structure. Finally, we shed light on some of the outstanding questions related to 3D FFM data interpretation and how this technique can become a central tool in the repertoire of surface science.

Oriented attachment induces fivefold twins by forming and decomposing high-energy grain boundaries.

Natural and synthetic nanoparticles composed of fivefold twinned crystal domains have distinct properties. The formation mechanism of these fivefold twinned nanoparticles is poorly understood. We used in situ high-resolution transmission electron microscopy combined with molecular dynamics simulations to demonstrate that fivefold twinning occurs through repeated oriented attachment of ~3-nanometer gold, platinum, and palladium nanoparticles. We discovered two different mechanisms for forming fivefold twinned nanoparticles that are driven by the accumulation and elimination of strain. This was accompanied by decomposition of grain boundaries and the formation of a special class of twins with a net strain of zero. These observations allowed us to develop a quantitative picture of the twinning process. The mechanisms provide guidance for controlling twin structures and morphologies across a wide range of materials.

Connecting energetics to dynamics in particle growth by oriented attachment using real-time observations.

The interplay between crystal and solvent structure, interparticle forces and ensemble particle response dynamics governs the process of crystallization by oriented attachment (OA), yet a quantitative understanding is lacking. Using ZnO as a model system, we combine in situ TEM observations of single particle and ensemble assembly dynamics with simulations of interparticle forces and responses to relate experimentally derived interparticle potentials to the underlying interactions. We show that OA is driven by forces and torques due to a combination of electrostatic ion-solvent correlations and dipolar interactions that act at separations well beyond 5 nm. Importantly, coalignment is achieved before particles reach separations at which strong attractions drive the final jump to contact. The observed barrier to attachment is negligible, while dissipative factors in the quasi-2D confinement of the TEM fluid cell lead to abnormal diffusivities with timescales for rotation much less than for translation, thus enabling OA to dominate.

Using Atom Dynamics to Map the Defect Structure Around an Impurity in Nano-Hematite.

The bulk behavior of materials is often controlled by minor impurities that create nonperiodic localized defect structures due to ionic size, symmetry, and charge balance mismatches. Here, we used transmission electron microscopy (TEM) of atom-resolved dynamics to directly map the topology of Fe vacancy clusters surrounding structurally incorporated U6+ in nanohematite (α-Fe2O3). Ab initio molecular dynamic simulations provided additional independent constraints on coupled U, Fe, and vacancy mobility in the solid. A clearer understanding of how such an apparently incompatible element can be accommodated by hematite emerged. The results were readily interpretable without the need for sophisticated data reconstruction methods, model structures, or ultrathin samples, and with the proliferation of aberration-corrected TEM facilities, the approach is accessible. Given sufficient z-contrast, the ability to observe impurity-vacancy structures by means of atom hopping can be used to directly probe the association of impurities and such defects in other materials, with promising applications across a broad range of disciplines.

Correlating inter-particle forces and particle shape to shear-induced aggregation/fragmentation and rheology for dilute anisotropic particle suspensions: A complementary study via capillary rheometry and in-situ small and ultra-small angle X-ray scattering.

HYPOTHESIS: Understanding the stability and rheological behavior of suspensions composed of anisotropic particles is challenging due to the complex interplay of hydrodynamic and colloidal forces. We propose that orientationally-dependent interactions resulting from the anisotropic nature of non-spherical sub-units strongly influences shear-induced particle aggregation/fragmentation and suspension rheological behavior. EXPERIMENTS: Wide-, small-, and ultra-small-angle X-ray scattering experiments were used to simultaneously monitor changes in size and fractal dimensions of boehmite aggregates from 6 to 10,000 Å as the sample was recirculated through an in-situ capillary rheometer. The latter also provided simultaneous suspension viscosity data. Computational fluid dynamics modeling of the apparatus provided a more rigorous analysis of the fluid flow. FINDINGS: Shear-induced aggregation/fragmentation was correlated with a complicated balance between hydrodynamic and colloidal forces. Multi-scale fractal aggregates formed in solution but the largest could be fragmented by shear. Orientationally-dependent interactions lead to a relatively large experimental suspension viscosity when the hydrodynamic force was small compared to colloidal forces. This manifests even at low boehmite mass fractions.