Nanoparticle Assembly and Oriented Attachment: Correlating Controlling Factors to the Resulting Structures.

Nanoparticle assembly and attachment are common pathways of crystal growth by which particles organize into larger scale materials with hierarchical structure and long-range order. In particular, oriented attachment (OA), which is a special type of particle assembly, has attracted great attention in recent years because of the wide range of material structures that result from this process, such as one-dimensional (1D) nanowires, two-dimensional (2D) sheets, three-dimensional (3D) branched structures, twinned crystals, defects, etc. Utilizing in situ transmission electron microscopy techniques, researchers observed orientation-specific forces that act over short distances (∼1 nm) from the particle surfaces and drive the OA process. Integrating recently developed 3D fast force mapping via atomic force microscopy with theories and simulations, researchers have resolved the near-surface solution structure, the molecular details of charge states at particle/fluid interfaces, inhomogeneity of surface charges, and dielectric/magnetic properties of particles that influence short- and long-range forces, such as electrostatic, van der Waals, hydration, and dipole-dipole forces. In this review, we discuss the fundamental principles for understanding particle assembly and attachment processes, and the controlling factors and resulting structures. We review recent progress in the field via examples of both experiments and modeling, and discuss current developments and the future outlook.

Influence of Concentration Gradients on Electroconvection at a Cation-Exchange Membrane Surface.

Membrane-based systems, such as electrodialysis, play an important role in desalination and industrial separation processes. Electrodialysis uses alternating anion- and cation-exchange membranes with a perpendicular electric field to generate concentrated and diluate streams from a feed solution. It is known that under overlimiting current conditions, reduced charge and mass transfer at the membrane interface leads to regions of high ion depletion generating instability and vortices termed electroconvection. While electroconvective mixing is known to directly impact the separation efficiency of electrodialysis, the influence of ion concentration gradients across the membrane experienced in a functional electrodialysis system is not known. Here, we report the influence of ion concentration gradients across a cation exchange membrane (Nafion) that is both aligned with and opposed to the applied electric field. Experiments were conducted by coflowing NaCl solutions of different concentrations (0.1-100 mM) on each side of the membrane, and electroconvection was visualized with a fluorescence dye (Rhodamine 6G). We obtained concentration profiles from fluorescence image data and systematically measured the thickness of the depletion boundary layer dBL under different conditions. We found smaller dBL values at a higher flow rate both with and without concentration gradients. Our results show that electroconvection is enhanced when the electric field is opposite to the direction of the concentration gradient.

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.

Towards predictive control of reversible nanoparticle assembly with solid-binding proteins.

Although a broad range of ligand-functionalized nanoparticles and physico-chemical triggers have been exploited to create stimuli-responsive colloidal systems, little attention has been paid to the reversible assembly of unmodified nanoparticles with non-covalently bound proteins. Previously, we reported that a derivative of green fluorescent protein engineered with oppositely located silica-binding peptides mediates the repeated assembly and disassembly of 10-nm silica nanoparticles when pH is toggled between 7.5 and 8.5. We captured the subtle interplay between interparticle electrostatic repulsion and their protein-mediated short-range attraction with a multiscale model energetically benchmarked to collective system behavior captured by scattering experiments. Here, we show that both solution conditions (pH and ionic strength) and protein engineering (sequence and position of engineered silica-binding peptides) provide pathways for reversible control over growth and fragmentation, leading to clusters ranging in size from 25 nm protein-coated particles to micrometer-size aggregate. We further find that the higher electrolyte environment associated with successive cycles of base addition eventually eliminates reversibility. Our model accurately predicts these multiple length scales phenomena. The underpinning concepts provide design principles for the dynamic control of other protein- and particle-based nanocomposites.

The Role of Ion Rotation in Ion Mobility: Ultrahigh-Precision Prediction of Ion Mobility Dependence on Ion Mass Distribution and Translational to Rotational Energy Transfer.

The role of ion rotation in determining ion mobilities is explored using the subtle gas phase ion mobility shifts based on differences in ion mass distributions between isotopomer ions that have been observed with ion mobility spectrometry (IMS) measurements. These mobility shifts become apparent for IMS resolving powers on the order of ∼1500 where relative mobilities (or alternatively momentum transfer collision cross sections; Ω) can be measured with a precision of ∼10 ppm. The isotopomer ions have identical structures and masses, differing only in their internal mass distributions, and their Ω differences cannot be predicted by widely used computational approaches, which ignore the dependence of Ω on the ion's rotational properties. Here, we investigate the rotational dependence of Ω, which includes changes to its collision frequency due to thermal rotation as well as the coupling of translational to rotational energy transfer. We show that differences in rotational energy transfer during ion-molecule collisions provide the major contribution to isotopomer ion separations, with only a minor contribution due to an increase in collision frequency due to ion rotation. Modeling including these factors allowed for differences in Ω to be calculated that precisely mirror the experimental separations. These findings also highlight the promise of pairing high-resolution IMS measurements with theory and computation for improved elucidation of subtle structural differences between ions.

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.

Radiolysis and Radiation-Driven Dynamics of Boehmite Dissolution Observed by In Situ Liquid-Phase TEM.

Over the last several decades, there have been several studies examining the radiation stability of boehmite and other aluminum oxyhydroxides, yet less is known about the impact of radiation on boehmite dissolution. Here, we investigate radiation effects on the dissolution behavior of boehmite by employing liquid-phase transmission electron microscopy (LPTEM) and varying the electron flux on the samples consisting of either single nanoplatelets or aggregated stacks. We show that boehmite nanoplatelets projected along the [010] direction exhibit uniform dissolution with a strong dependence on the electron dose rate. For nanoplatelets that have undergone oriented aggregation, we show that the dissolution occurs preferentially at the particles at the ends of the stacks that are more accessible to bulk solution than at the others inside the aggregate. In addition, at higher dose rates, electrostatic repulsion and knock-on damage from the electron beam causes delamination of the stacks and dissolution at the interfaces between particles in the aggregate, indicating that there is a threshold dose rate for electron-beam enhancement of dissolution of boehmite aggregates.

27 Al NMR diffusometry of Al13 Keggin nanoclusters.

Although nanometer-sized aluminum hydroxide clusters (i.e., ϵ-Al13 , [Al13 O4 (OH)24 (H2 O)12 ]7+ ) command a central role in aluminum ion speciation and transformations between minerals, measurement of their translational diffusion is often limited to indirect methods. Here, 27 Al pulsed field gradient stimulated echo nuclear magnetic resonance (PFGSTE NMR) spectroscopy has been applied to the AlO4 core of the ϵ-Al13 cluster with complementary theoretical simulations of the diffusion coefficient and corresponding hydrodynamic radii from a boundary element-based calculation. The tetrahedral AlO4 center of the ϵ-Al13 cluster is symmetric and exhibits only weak quadrupolar coupling, which results in favorable T1 and T2 27 Al NMR relaxation coefficients for 27 Al PFGSTE NMR studies. Stokes-Einstein relationship was used to relate the 27 Al diffusion coefficient of the ϵ-Al13 cluster to the hydrodynamic radius for comparison with theoretical simulations, dynamic light scattering from literature, and previously published 1 H PFGSTE NMR studies of chelated Keggin clusters. This first-of-its-kind observation proves that 27 Al PFGSTE NMR diffusometry can probe symmetric Al environments in polynuclear clusters of greater molecular weight than previously considered.

Shear stress dependence of force networks in 3D dense suspensions.

The geometric organization and force networks of 3D dense suspensions that exhibit both shear thinning and thickening have been examined as a function of varying strength of interparticle attractive interactions using lubrication flow discrete element simulations. Significant rearrangement of the geometric topology does not occur at either the local or global scale as these systems transition across the shear thinning and shear thickening regimes. In contrast, massive rearrangements in the balance of attractive, lubrication, and contact forces are observed with interesting behavior of network growth and competition. In agreement with prior work, in shear thinning regions the attractive force is dominant, however as the shear thickening region is approached there is growth of lubrication forces. Lubrication forces oppose the attraction forces, but as viscosity continues to increase under increasing shear stress, the lubrication forces are dominated by contact forces that also resist attraction. Contact forces are the dominant interactions during shear thickening and are an order of magnitude higher than their values in the shear-thinning regime. At high attractive interaction strength, contact networks can form even under shear thinning conditions, however high shear stress is still required before contact networks become the driving mechanism of shear thickening. Analysis of the contact force network during shear thickening generally indicates a uniformly spreading network that rapidly forms across empty domains; however the growth patterns exhibit structure that is significantly dependent upon the strength of interparticle interactions, indicating subtle variations in the mechanism of shear thickening.

Visualization of Aluminum Ions at the Mica Water Interface Links Hydrolysis State-to-Surface Potential and Particle Adhesion.

When hydrolyzable cations such as aluminum interact with solid-water interfaces, macroscopic interfacial properties (e.g., surface charge and potential) and interfacial phenomena (e.g., particle adhesion) become tightly linked with the microscopic details of ion adsorption and speciation. We use in situ atomic force microscopy to directly image individual aluminum ions at a mica-water interface and show how adsorbate populations change with pH and aluminum activity. Complementary streaming potential measurements then allow us to build a triple layer model (TLM) that links surface potentials to adsorbate populations, via equilibrium binding constants. Our model predicts that hydrolyzed species dominate the mica-water interface, even when unhydrolyzed species dominate the solution. Ab initio molecular dynamics (AIMD) simulations confirm that aluminum hydrolysis is strongly promoted at the interface. The TLM indicates that hydrolyzed adsorbates are responsible for surface-potential inversions, and we find strong correlations between hydrolyzed adsorbates and particle-adhesion forces, suggesting that these species mediate adhesion by chemical bridging.

Transport of Colloidal Particles in Microscopic Porous Medium Analogues with Surface Charge Heterogeneity: Experiments and the Fundamental Role of Single-Bead Deposition.

Understanding colloid transport in subsurface environments is challenging because of complex interactions among colloids, groundwater, and porous media over several length scales. Here, we report a versatile method to assemble bead-based microfluidic porous media analogues with chemical heterogeneities of different configurations. We further study the transport of colloidal particles through a family of porous media analogues that are randomly packed with oppositely charged beads with different mixing ratios. We recorded the dynamics of colloidal particle deposition at the level of single grains. From these, the maximum surface coverage (θmax = 0.051) was measured directly. The surface-blocking function and the deposition coefficient (kpore = 3.56 s-1) were obtained. Using these pore-scale parameters, the transport of colloidal particles was modeled using a one-dimensional advection-dispersion-deposition equation under the assumption of irreversible adsorption between oppositely charged beads and colloids, showing very good agreement with experimental breakthrough curves and retention profiles at the scale of the entire porous medium analogue. This work presents a new approach to fabricate chemically heterogeneous porous media in a microfluidic device that enables the direct measurement of pore-scale colloidal deposition. Compared with the conventional curve-fitting method for deposition constant, our approach allows quantitative prediction of colloidal breakthrough and retention via coupling of direct pore-scale measurements and an advection-dispersion-deposition model.

Trends in mica-mica adhesion reflect the influence of molecular details on long-range dispersion forces underlying aggregation and coalignment.

Oriented attachment of nanocrystalline subunits is recognized as a common crystallization pathway that is closely related to formation of nanoparticle superlattices, mesocrystals, and other kinetically stabilized structures. Approaching particles have been observed to rotate to achieve coalignment while separated by nanometer-scale solvent layers. Little is known about the forces that drive coalignment, particularly in this "solvent-separated" regime. To obtain a mechanistic understanding of this process, we used atomic-force-microscopy-based dynamic force spectroscopy with tips fabricated from oriented mica to measure the adhesion forces between mica (001) surfaces in electrolyte solutions as a function of orientation, temperature, electrolyte type, and electrolyte concentration. The results reveal an ∼60° periodicity as well as a complex dependence on electrolyte concentration and temperature. A continuum model that considers the competition between electrostatic repulsion and van der Waals attraction, augmented by microscopic details that include surface separation, water structure, ion hydration, and charge regulation at the interface, qualitatively reproduces the observed trends and implies that dispersion forces are responsible for establishing coalignment in the solvent-separated state.

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.

From Yielding to Shear Jamming in a Cohesive Frictional Suspension.

Simulations are used to study the steady shear rheology of dense suspensions of frictional particles exhibiting discontinuous shear thickening and shear jamming, in which finite-range cohesive interactions result in a yield stress. We develop a constitutive model that combines yielding behavior and shear thinning at low stress with the frictional shear thickening at high stresses, in good agreement with the simulation results. This work shows that there is a distinct difference between solids below the yield stress and in the shear-jammed state, as the two occur at widely separated stress levels, with an intermediate region of stress in which the material is flowable.

Silver Nanocube and Nanobar Growth via Anisotropic Monomer Addition and Particle Attachment Processes.

Understanding the growth mechanism of noble metal nanocrystals during solution synthesis is of significant importance for shape and property control. However, much remains unknown about the growth pathways of metal nanoparticles due to the lack of direct observation. Using an in situ transmission electron microscopy technique, we directly observed Ag nanocube and nanobar growth in an aqueous solution through both classical monomer-by-monomer addition and nonclassical particle attachment processes. During the particle attachment process, Ag nanocubes and nanobars were formed via both oriented and nonoriented attachment. Our calculations, along with the dynamics of the observed attachment, showed that the van der Waals force overcomes hydrodynamic and friction forces and drives the particles toward each other at separations of 10-100 nm in our experiments. During classical growth, anisotropic growth was also revealed, and the resulting unsymmetrical shape constituted an intermediate state for Ag nanocube growth. We hypothesized that the temporary symmetry breaking resulted from different growth rates on (001) surfaces due to a local surface concentration variation caused by the imbalance between the consumption of Ag+ near the surface and the diffusion of Ag+ from the bulk to the surface.

Review of the Scientific Understanding of Radioactive Waste at the U.S. DOE Hanford Site.

This Critical Review reviews the origin and chemical and rheological complexity of radioactive waste at the U.S. Department of Energy Hanford Site. The waste, stored in underground tanks, was generated via three distinct processes over decades of plutonium extraction operations. Although close records were kept of original waste disposition, tank-to-tank transfers and conditions that impede equilibrium complicate our understanding of the chemistry, phase composition, and rheology of the waste. Tank waste slurries comprise particles and aggregates from nano to micro scales, with varying densities, morphologies, heterogeneous compositions, and complicated responses to flow regimes and process conditions. Further, remnant or changing radiation fields may affect the stability and rheology of the waste. These conditions pose challenges for transport through conduits or pipes to treatment plants for vitrification. Additionally, recalcitrant boehmite degrades glass quality and the high aluminum content must be reduced prior to vitrification for the manufacture of waste glass of acceptable durability. However, caustic leaching indicates that boehmite dissolves much more slowly than predicted given surface normalized rates. Existing empirical models based on ex situ experiments and observations generally only describe material balances and have not effectively predicted process performance. Recent advances in the areas of in situ microscopy, aberration-corrected transmission electron microscopy, theoretical modeling across scales, and experimental methods for probing the physics and chemistry at mineral-fluid and mineral-mineral interfaces are being implemented to build robustly predictive physics-based models.

Simulation of shear thickening in attractive colloidal suspensions.

The influence of attractive forces between particles under conditions of large particle volume fraction, ϕ, is addressed using numerical simulations which account for hydrodynamic, Brownian, conservative and frictional contact forces. The focus is on conditions for which a significant increase in the apparent viscosity at small shear rates, and possibly the development of a yield stress, is observed. The high shear rate behavior for Brownian suspensions has been shown in recent work [R. Mari, R. Seto, J. F. Morris and M. M. Denn PNAS, 2015, 112, 15326-15330] to be captured by the inclusion of pairwise forces of two forms, one a contact frictional interaction and the second a repulsive force often found in stabilized colloidal dispersions. Under such conditions, shear thickening is observed when shear stress is comparable to the sum of the Brownian stress, kT/a3, and a characteristic stress based on the combination of interparticle force, i.e. σ ∼ F0/a2 with kT the thermal energy, F0 the repulsive force scale and a the particle radius. At sufficiently large ϕ, this shear thickening can be very abrupt. Here it is shown that when attractive interactions are present with the noted forces, the shear thickening is obscured, as the viscosity shear thins with increasing shear rate, eventually descending from an infinite value (yield stress conditions) to a plateau at large stress; this plateau is at the same level as the large-shear rate viscosity found in the shear thickened state without attractive forces. It is shown that this behavior is consistent with prior observations in shear thickening suspensions modified to be attractive through depletion flocculation [V. Gopalakrishnan and C. F. Zukoski J. Rheol., 2004, 48, 1321-1344]. The contributions of the contact, attractive, and hydrodynamics forces to the bulk stress are presented, as are the contact networks found at different attractive strengths.

Bead-Based Microfluidic Sediment Analogues: Fabrication and Colloid Transport.

Mobile colloids can act as carriers for low-solubility contaminants in the environment. However, the dominant mechanism for this colloid-facilitated transport of chemicals is unclear. Therefore, we developed a bead-based microfluidic platform of sediment analogues and measured both single and population transport of model colloids. The porous medium is assembled through a bead-by-bead injection method. This approach has the versatility to build both electrostatically homogeneous and heterogeneous media at the pore scale. A T-junction at the exit also allowed for encapsulation and enumeration of colloids effluent at single particle resolution to give population dynamics. Tortuosity calculated from pore-scale trajectory analysis and its comparison with lattice Boltzmann simulations revealed that transport of colloids was influenced by the size exclusion effect. The porous media packed by positively and negatively charged beads into two layers showed distinctive colloidal particle retention and significant remobilization and re-adsorption of particles during water flushing. We demonstrated the potential of our method to fabricate porous media with surface heterogeneities at the pore scale. With both single and population dynamics measurement, our platform has the potential to connect pore-scale and macroscale colloid transport on a lab scale and to quantify the impact of grain surface heterogeneities that are natural in the subsurface environment.

Biomolecular electrostatics and solvation: a computational perspective.

An understanding of molecular interactions is essential for insight into biological systems at the molecular scale. Among the various components of molecular interactions, electrostatics are of special importance because of their long-range nature and their influence on polar or charged molecules, including water, aqueous ions, proteins, nucleic acids, carbohydrates, and membrane lipids. In particular, robust models of electrostatic interactions are essential for understanding the solvation properties of biomolecules and the effects of solvation upon biomolecular folding, binding, enzyme catalysis, and dynamics. Electrostatics, therefore, are of central importance to understanding biomolecular structure and modeling interactions within and among biological molecules. This review discusses the solvation of biomolecules with a computational biophysics view toward describing the phenomenon. While our main focus lies on the computational aspect of the models, we provide an overview of the basic elements of biomolecular solvation (e.g. solvent structure, polarization, ion binding, and non-polar behavior) in order to provide a background to understand the different types of solvation models.

Effect of bubbles and silica dissolution on melter feed rheology during conversion to glass.

Nuclear-waste melter feeds are slurry mixtures of wastes with glass-forming and glass-modifying additives (unless prefabricated frits are used), which are converted to molten glass in a continuous electrical glass-melting furnace. The feeds gradually become continuous glass-forming melts. Initially, the melts contain dissolving refractory feed constituents that are suspended together with numerous gas bubbles. Eventually, the bubbles escape, and the melts homogenize and equilibrate. Knowledge of various physicochemical properties of the reacting melter feed is crucial for understanding the feed-to-glass conversion that occurs during melting. We studied the melter feed viscosity during heating and correlated it with the volume fractions of dissolving quartz (SiO2) particles and the gas phase. The measurements were performed with a rotating spindle rheometer on the melter feed heated at 5 K/min, starting at several different temperatures. The effects of undissolved quartz particles, gas bubbles, and compositional inhomogeneity on the melter feed viscosity were determined by fitting a linear relationship between the logarithm of viscosity and the volume fractions of suspended phases.