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.

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.

Tracking nitrite's deviation from Stokes-Einstein predictions with pulsed field gradient 15N NMR spectroscopy.

Predicting the behavior of oxyanions in radioactive waste stored at the Department of Energy legacy nuclear sites requires the development of novel analytical methods. This work demonstrates 15N pulsed field gradient nuclear magnetic resonance spectroscopy to quantify the diffusivity of nitrite. Experimental results, supported by molecular dynamics simulations, indicate that the diffusivity of free hydrated nitrite exceeds that of free hydrated sodium despite the greater hydrodynamic radius of nitrite. Investigations are underway to understand how the compositional and dynamical heterogeneities of the ion networks at high concentrations affect rheological and transport properties.

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.

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.

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.

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.

SimELIT: A Novel GUI-Based Comprehensive Ion Trajectory Simulation Software for Mass Spectrometry.

Ion trajectory simulation in mass spectrometry systems from injection to detection is technically challenging but very important for better understanding the ion dynamics in instrument development. Here, we present SimELIT (Simulator of Eulerian and Lagrangian Ion Trajectories), a novel ion trajectory simulation platform. SimELIT is built upon a suite of multiphysics solvers compiled into OpenFOAM (an open-source numerical solver library particularly used for computational mechanics), with a simple web-based graphical user interface (GUI) allowing users to define the details of OpenFOAM cases and run simulations. SimELIT is a modular program and can provide extensions of physics (e.g., gas flows, electrodynamic fields) and thus enable ion trajectory simulations from the ion source to detector. The current version (SimELIT) provides two numerical solvers for ion trajectory simulations─(1) a Lagrangian particle tracker in vacuum and (2) a Eulerian ion density solver in background gas in the presence of electric fields. Here, we describe the architecture of SimELIT, including its use of Docker and the React Framework, and demonstrate the computation of ion trajectories of multiple m/z values in a static/linear voltage drop in vacuum (across a 1 m long flight tube). Further, the drift motion of ions under 1 Torr pressure conditions in a static background (N2) gas through a 20 V/cm static electric field is shown. The results produced from SimELIT were compared with SIMION and theoretical estimates. In addition, we report the computation of ion trajectories in electrodynamic fields within a planar FAIMS device operating at atmospheric pressure.

Evaluation of Atopic Dermatitis and Cutaneous Infectious Disorders Using Sequential Pattern Mining: A Nationwide Population-Based Cohort Study.

According to previous studies, the increased risk of cutaneous infectious disorders in patients with atopic dermatitis (AD) is related to impaired epidermal function, abnormal systemic immune function, and lower antimicrobial peptides. In this study, we analyzed the association between AD and cutaneous infectious disorders in the real world using sequential pattern mining (SPM). We analyzed National Health Insurance data from 2010-2013 using SPM to identify comorbid cutaneous infectious diseases and the onset durations of comorbidities. Patients with AD were at greater risk for molluscum contagiosum (adjusted odds ratio (aOR), 5.273), impetigo (aOR, 2.852), chickenpox (aOR, 2.251), otitis media (aOR, 1.748), eczema herpeticum (aOR, 1.292), and viral warts (aOR, 1.105). In SPM analysis, comorbidity of 1.06% shown in molluscum contagiosum was the highest value, and the duration of 77.42 days documented for molluscum contagiosum was the shortest onset duration among all the association rules. This study suggests that AD is associated with an increased risk of cutaneous infectious disorders. In particular, care should be taken regarding its high relevance with impetigo, molluscum contagiosum, and otitis media, which may help in preventing AD from worsening through appropriately preventing and managing the condition.

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.

Predictive Theoretical Framework for Dynamic Control of Bioinspired Hybrid Nanoparticle Self-Assembly.

At-will tailoring of the formation and reconfiguration of hierarchical structures is a key goal of modern nanomaterial design. Bioinspired systems comprising biomacromolecules and inorganic nanoparticles have potential for new functional material structures. Yet, consequential challenges remain because we lack a detailed understanding of the temporal and spatial interplay between participants when it is mediated by fundamental physicochemical interactions over a wide range of scales. Motivated by a system in which silica nanoparticles are reversibly and repeatedly assembled using a homobifunctional solid-binding protein and single-unit pH changes under near-neutral solution conditions, we develop a theoretical framework where interactions at the molecular and macroscopic scales are rigorously coupled based on colloidal theory and atomistic molecular dynamics simulations. We integrate these interactions into a predictive coarse-grained model that captures the pH-dependent reversibility and accurately matches small-angle X-ray scattering experiments at collective scales. The framework lays a foundation to connect microscopic details with the macroscopic behavior of complex bioinspired material systems and to control their behavior through an understanding of both equilibrium and nonequilibrium characteristics.

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.

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.

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.

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.

Correlation function approach for diffusion in confined geometries.

This paper describes a formalism for extracting spatially varying transport coefficients from simulations of a molecular fluid in a nanochannel. This approach is applied to self-diffusion of a Lennard-Jones fluid confined between two parallel surfaces. A numerical grid is laid over the domain confining the fluid, and fluid properties are projected onto the grid cells. The time correlation functions between properties in different grid cells are calculated and can be used as the basis for a fitting procedure for extracting spatially varying diffusion coefficients from the simulation. Results for the Lennard-Jones system show that transport behavior varies sharply near the liquid-solid boundary and that the changes depend on the details of the liquid-solid interaction. A quantitative difference between the reduced and detailed models is discussed. It is found that the difference could be associated with assumptions about the form of the transport equations at molecular scales in lieu of problems with the method itself. The study suggests that this approach to fitting molecular simulations to continuum equations may guide the development of appropriate coarse-grained equations to model transport phenomena at nanometer scales.

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.