Search NASA⌕ Search

SEARCH · Search NASA

Results for “counterions”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Combined Machine Learning and Molecular Dynamics Reveal Two States of Hydration of a Single Functional Group of Cationic Polymeric Brushes

The state of hydration of a macromolecular system regulates a plethora of different properties of such a system. In this article, we develop a novel machine learning (ML) approach, based on the unsupervised clustering algorithm, for probing the hydration behavior of the {N(CH 3 ) 3 } + functional group of the PMETAC [Poly(2-(methacryloyloxy)ethyl trimethylammonium chloride] polyelectrolyte (PE) brush system. The PE brushes and the brush-supported water molecules and counterions (chloride ions) are first described using all-atom molecular dynamics (MD) simulations. The simulation data is subsequently used in our ML framework to identify that (1) the {N(CH 3 ) 3 } + functional groups of the PMETAC brushes have two distinct hydration states with one state (state 1) being characterized by less structured water molecules and the other state (state 2) being characterized by more structured water molecules and (2) an enhancement in the brush grafting density leads to the progressive dissapparenace of state 2. An increase in the grafting density increases the number of chloride counterions in a given volume around the {N(CH 3 ) 3 } + functional group and increases the number of shared water molecules between the {N(CH 3 ) 3 } + and Cl - . The chloride counterions are associated with a hydration layer with much less structured water molecules. Therefore, with an increase in the grafting density, an increase in the percentage of shared water molecules leads to the prevalence of the hydration state [of the {N(CH 3 ) 3 } + moiety] with less structured water molecules. Finally, we explain how the present findings are commensurate with two key previous related results, namely a significantly large chloride ion mobility inside the PMETAC brush layer and the {N(CH 3 ) 3 } + -Cl - average distance remaining independent of the PMETAC brush grafting density. Furthermore, we anticipate that the combined ML-MD-simulation approach proposed in this study can be adapted to probe other soft matter systems to reveal new insights of the underlying mechanisms of emergent phenomenon.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Modeling salt adsorption in electrical double layer for capacitive deionization

Abstract This work presents an electrical double layer (EDL) model for capacitive deionization (CDI) by intercorrelating the salt adsorption, the electrode voltage, and the surface charge density. The counterion condensation is considered a crucial contribution to salt adsorption because of the high surface charge density in the charged micropores. The counterions condense when the surface charge density exceeds a critical value, which is predicted by the counterion condensation theory. The EDL model accurately correlates the NaCl salt adsorption and electrode voltage with various surface charge densities and well predicts the salt adsorption in a wide range of external salt concentrations. The EDL model reproduces the non‐monotonic relationship between the salt adsorption and the specific surface area observed in experiments and provides an explanation from the modeling perspective. Exhibiting superior accuracy and predictability for salt adsorption, the EDL model could serve as an enabling tool contributing to the development of CDI processes.

Lin, Yu‐Jeng↗
DOE OSTI2024

CO 2 hydrogenation over rhodium cluster catalyst nucleated within a manganese oxide framework

Rhodium-based manganese oxide frameworks were explored as a prototype for carbon dioxide reactive capture and conversion. Three-dimensional frameworks of MnOx were utilized as support structures to isolate Rh metal centers. V, Na, and Zn were introduced as counterions to stabilize the structure and for their beneficial effect as promoters. Here, with this multicomponent catalyst, Rh active centers with MnOxs and varied counterions, we were able to selectively tune the catalytic performance of the material via the choice of counterion and structure of the host material. With cryptomelane-type tunnel manganese oxides octahedral molecular sieve (OMS2), we found that Rh-V-OMS2 was highly stable even after 48 hours on stream with a reaction rate of around 1.5x10 -4 mol CO 2 /g Rh /s, surpassing the net reactivity of other initially more active combinations. Furthermore, during CO 2 hydrogenation, in situ XAFS showed that single Rh atoms nucleated into nanoparticles/ sub-nanometer clusters with a coordination number of 5.5 or less. Our finding of the correlation between the reaction rate and particle size offers the potential for enhanced control over the reaction rate by tuning particle size. Our activity study with control experiments demonstrates that the activities of the catalysts are proved due to the unique metal support interaction offered by the Rh-X-MnO.

10 SYNTHETIC FUELS↗
DOE OSTI2025

Interaction of Planar Brushes with End-Functionalized Carboxylic or Amine Groups: Water Structure, Hydrogen Bonds, and Electrostatic Correlations

In this work we present a detailed analysis of two parallel plates grafted with polyethylene glycol (PEG) terminated with carboxylic or amine groups, and discuss the cases of opposite and same charge in the presence of electrolytes. We analyze the screening of charged groups by electrolytes, the role of hydrogen bonds in enhancing attraction or repulsion, hydration forces, competitive binding with counterions, and the hydrogen bond network, crucial to stabilizing the system. Our simulations reveal that, for large plate separation, the system is well described as a “Stern layer”, characterized by strong counterion binding and a subsequent weak diffuse layer. For short plate separation, the system is characterized by strong correlations among end groups and counterions. We also discuss the relevance of charged carboxylic/amine groups in different problems with a particular emphasis on nanoparticle assembly.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Screening of Coulombic Interactions To Achieve a Higher Power Factor in Conjugated Polymers

Thermoelectric properties of conducting polymers typically suffer from molecular chain disordering, as charge transport is predominantly controlled by morphology. This is especially more problematic when counterions are introduced to tune the carrier concentration for optimal thermoelectric performance, which disturbs the morphology further. Here, in this work, we introduce a new avenue for enhancing thermoelectric properties without needing to regulate the morphology, namely, by controlling the coulombic interaction between polarons and counterions. We perform in situ de-doping thermoelectric experiments over 3 orders of magnitude change in electrical conductivity of three distinct thermoelectric polymers, namely, poly(3-hexylthiophene-2,5-diyl) (P3HT), poly[2,5-bis(3-dodecylthiophen-2-yl)thieno[3,2-b]thiophene] (PBTTT-C 12 ), and poly[2,5-(2-octyldodecyl)-3,6-diketopyrrolopyrrole-alt-5,5-(2,5-di(thien-2-yl)thieno[3,2-b]thiophene)] (OD-PDPP2T-TT) conjugated polymers, followed by grazing-incidence wide-angle X-ray scattering (GIWAXS) to study their respective morphologies. We demonstrate a 9-fold enhancement in the thermoelectric power factor in OD-PDPP2T-TT compared to PBTTT-C 12 and link it to the coulombic screening of charge carriers, including in the optimally doped regime. We support this hypothesis using Boltzmann transport equations and show that, in both P3HT and PBTTT-C 12 , as the polymer is doped, impurity scattering remains the dominant scattering mechanism, while in OD-PDPP2T-TT, the scattering mechanism changes from impurity to acoustic phonon limited, resulting in more effective screening of ionized counterions. Our results provide an additional knob to enhance the fundamental understanding of thermoelectric physics of conducting polymers and provide a pathway to achieve higher performance in the field of organic thermoelectrics.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Lewis acid–base pair doping of p-type organic semiconductors

Doping is required to increase the electrical conductivity of organic semiconductors for uses in electronic and energy conversion devices. The limited number of commonly used p-type dopants suggests that new dopants or doping mechanisms could improve the efficiency of doping and provide new means for processing doped polymers. Drawing on Lewis acid–base pair chemistry, we combined Lewis acid dopant B(C 6 F 5 ) 3 (BCF) with the weak Lewis base benzoyl peroxide (BPO). The detailed behavior of p-type doping of the model polymer poly(3-hexylthiophene) (P3HT) with this Lewis acid–base pair in solution was examined. Solution 19F-NMR spectra confirmed the formation of the expected counterion, as well as side products from reactions with solvent. BCF : BPO was also found to efficiently dope a range of semiconducting polymers with varying chemical structures demonstrating that the BCF : BPO combination has an effective electron affinity of at least 5.3 eV. In thin films of regioregular P3HT cast from the doped solutions, delocalized polarons formed due to the large counterions leading to a large polaron-counterion distance. At and above 0.2 eq. BCF : BPO doping, amorphous areas of the film became doped, disrupting the structural order of the films. Despite the change in structural order, thin films of regioregular P3HT doped with 0.2 eq. BCF : BPO had a conductivity of 25 S cm -1 . This study demonstrates the effectiveness of a two-component Lewis acid–base doping mechanism and suggests additional two-component Lewis acid–base chemistries should be explored.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Thiourea-Derived Chelating Ligands and Their Organometallic Compounds: Investigations into Their Anticancer Activity

Thiones have been investigated as ligands in metal complexes with catalytic and biological activity. We report the synthesis, characterization, and biological evaluation of a series of MII/III complexes of the general formulae [M II (cym)(L)Cl]X (cym = η 6 -p-cymene) or [M III (Cp*)(L)Cl]X (Cp* = η 5 -pentamethylcyclopentadienyl), where X = Cl - or PF 6 - , and L represents heterocyclic derivatives of thiourea. The thiones feature a benzyl-triazolyl pendant and they act as bidentate ligands via N,S-coordination to the metal centers. Several derivatives have been investigated by single-crystal X-ray diffraction analysis. NMR investigations showed a counterion-dependent shift of several protons due to the interaction with the counterions. These NMR investigations were complemented with X-ray diffraction analysis data and the effects of different counterions on the secondary coordination sphere were also investigated by DFT calculations. In biological studies, the Ir benzimidazole derivative was found to accumulate in the cytoplasm and it was the most cytotoxic derivative investigated.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

From +I to +IV, Alkalis to Actinides: Capturing Cations across the Periodic Table with Keggin Polyoxometalate Ligands

Coordination chemistry trends across the periodic table are often difficult to probe experimentally due to limitations in finding a versatile but consistent chelating platform that can accommodate various elements without changing its coordination mode. Herein, we present new metal/ligand systems covering a wide range of ionic radii, charges, and elements. Five different ligands derived from the Keggin structure (HBW 11 O 39 8– , PW 11 O 39 7– , SiW 11 O 39 8– , GeW 11 O 39 8– , and GaW 11 O 39 9– ) were successfully crystallized with six different cations (Na + , Sr 2+ , Ba 2+ , La 3+ , Ce 4+ , and Th 4+ ) and characterized by single-crystal X-ray diffraction. Twenty-five new compounds were obtained by using Cs + as the counterion, yielding a consistent base formula of Cs x [M(XW 11 O 39 ) 2 ]·nH 2 O. Despite having a similar first-coordination sphere geometry (i.e., 8-coordinated), the nature of the central cation was found to impact the long-range geometry of the complexes. This unique crystallographic data set shows that, despite the traditional consensus, the local geometry of the cation (i.e., metal–oxygen bond distance) is not enough to depict the full impact of the complexed metal ion. The bending and twisting of the complexes, as well as ligand–ligand distances, were all impacted by the nature of the central cation. We also observed that counterions play a critical role by stabilizing the geometry of the M(XW 11 ) 2 complex and directing complex–complex interactions in the lattice. We also define certain structural limits for this type of complex, with the large Ba 2+ ion seemingly approaching those limits. Finally, this study thus lays the foundation for capturing the coordination chemistry of other rarer elements across the periodic table such as Ra 2+ , Ac 3+ , Bk 4+ , Cf 3+ , etc.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Ionic Tunability of Conjugated Polyelectrolyte Solutions

Conjugated polyelectrolytes (CPEs), which combine a π-conjugated polymer backbone with pendant ionic functionalities, offer an opportunity for electrostatic control of materials properties. In this work, the mesoscale morphology and physical properties of a high-mobility conjugated polyelectrolyte are tuned by the addition of salt, variation of the charge-compensating counterion, and complexation with an oppositely charged polyelectrolyte containing the same π-conjugated backbone. In systems with a single polyelectrolyte species, added ions screen the electrostatic repulsions stabilizing the gel-phase, resulting in the dissolution of ionic cross-links and hydrophobic collapse. Further, exchanging the charge compensating counterion is found to enable finer structural control of the solution behavior and modulation of the self-doping behavior. Finally, novel CPE–CPE complexes resulting in dense solutions and gels of semiconductive material are produced by combination of oppositely charged polyelectrolytes. Such concentrated CPE formulations should be useful materials for mixed electronic–ionic conduction and pseudocapacitative energy storage.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Charge-Density-Specific Response of Grafted Polyelectrolytes to Electric Fields: Bending or Tilting?

The response of polyelectrolytes (PEs) to applied electric fields drive applications in energetics, diagnostics, materials development, and many more. Here we employ all-atom molecular dynamics (MD) simulations to probe the response of grafted PE brushes to axial electric fields. For PEs with large charge densities, the electric field triggers a left-right asymmetric distribution of counterions around the PE backbone: consequently, depending on the location (left or right), there is an unequal screening and an unequal force on the PE functional groups causing a bending-driven brush height reduction. However, for the weakly charged PEs, the electric field causes a uniform distribution of the counterions across the brush enforcing a uniform (without left-right asymmetry) partially unscreened PE brush: therefore, all the brush segments experience similar force causing a brush tilting-driven brush height reduction. As a result, such bending versus tilting responses is commensurate with the electric-field-driven increase (decrease) in the flexibility of strongly (weakly) charged PEs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Highly Altered State of Proton Transport in Acid Pools in Charged Reverse Micelles

Transport mechanisms of solvated protons of 1 M HCl acid pools, confined within reverse micelles (RMs) containing the negatively charged surfactant sodium bis(2-ethylhexyl) sulfosuccinate (NaAOT) or the positively charged cetyltrimethylammonium bromide (CTABr), are analyzed with reactive force field simulations to interpret dynamical signatures from TeraHertz absorption and dielectric relaxation spectroscopy. We find that the forward proton hopping events for NaAOT are further suppressed compared to a nonionic RM, while the Grotthuss mechanism ceases altogether for CTABr. We attribute the sluggish proton dynamics for both charged RMs as due to headgroup and counterion charges that expel hydronium and chloride ions from the interface and into the bulk interior, thereby increasing the pH of the acid pools relative to the nonionic RM. For charged NaAOT and CTABr RMs, the localization of hydronium near a counterion or conjugate base reduces the Eigen and Zundel configurations that enable forward hopping. Thus, localized oscillatory hopping dominates, an effect that is most extreme for CTABr in which the proton residence time increases dramatically such that even oscillatory hopping is slow.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Similar but Different: Structural and Spectroscopic Characterization of Series of Neodymium, Europium, Americium, and Curium Coordination Complexes

Series of curium, americium, europium, and neodymium complexes with five different polyoxometalate (POM) ligands were structurally and spectroscopically characterized: Na x Cs y H z [M(XW 11 O 39 ) 2 ]· n H 2 O, where M = Nd 3+ , Eu 3+ , Am 3+ , or Cm 3+ and X = B 3+ , Ga 3+ , Si 4+ , Ge 4+ , or P 5+ . This first “serial approach” on transplutonium chemistry was allowed by minimizing the amount of f-element needed per synthesis down to ∼10 μg. This offers a unique opportunity to contrast structural and spectroscopic properties of trivalent actinide with those of lanthanide compounds under identical conditions. The results showcase that although actinide(III) and lanthanide(III) can make isostructural complexes, their solid-state coordination chemistries deviate significantly. The curium and americium compounds were found to be more uniform than their europium and neodymium counterparts, respectively. The particular symmetry of the Cm-POMs also enables the observation of rarely seen emission bands, which hints at vibronic coupling or new emissive pathways. Additionally, the Cm-POM and Eu-POM complexes were characterized via steady-state and time-resolved excitation and emission spectroscopy in solutions containing five different counterions (Li + , Na + , K + , Rb + , or Cs + ), representing 50 combinations of alkali/[M III (XW 11 ) 2 ] n − complexes. The Cm-POMs respond to the alkali counterions by improving symmetry from Li + to Cs + , while the Eu-POMs do the opposite. The luminescence lifetimes of the Cm-POMs also depart from the ideal Kimura equation. Finally, this study provides a rare, and perhaps the most comprehensive, experimental data set comparing heavy actinides and lanthanides via single crystal XRD plus, solution-state, solid-state, and time-resolved emission spectroscopy and demonstrates that lanthanides are only a coarse approximation for actinide elements.

chemical structure↗
DOE OSTI2025

Determination of Intermolecular Distances in Dilute Solutions of Macroions and Their Direct Correlations with Self-Assembly and Macrophase Transitions

Here, this work demonstrates new attempts to determine the intermolecular distances between charged macroionic solutes in their dilute solutions, which regulate the solute’s microphase (self-assembly) and macrophase transitions. Small-angle X-ray scattering (SAXS) and analytical ultracentrifuge (AUC) techniques were applied to determine the intermolecular distances for the 2.42 nm-sized, spherical uranyl peroxide molecular cluster {U 60 } in their self-assembled states in dilute aqueous solution and concentrated phases, respectively. The counterion-mediated attraction among {U 60 } leads to characteristic, inter-{U 60 } distances in solutions, which increase gradually with decreasing the strength of introduced counterions (e.g., lower valency). The gradual increment of inter-{U 60 } distance demonstrates a nice correlation with the macroscopic phase transitions of {U 60 }, from single crystals to concentrated fluids containing rigid 2-D sheets, then to dilute solutions containing a small amount of standalone, floating 2-D sheets of {U 60 }, and finally to dilute solutions containing a limited amount of self-assembled single-layered, spherical blackberry structures of different sizes from the bending of more flexible 2-D sheets.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Establishing Self-Dopant Design Principles from Structure–Function Relationships in Self-n-Doped Perylene Diimide Organic Semiconductors

Self-doping is a particular doping method that has been applied to a wide range of organic semiconductors. However, there is a lack of understanding regarding the relationship between dopant structure and function. A structurally diverse series of self-n-doped perylene diimides (PDIs) is investigated to study the impact of steric encumbrance, counterion selection, and dopant/PDI tether distance on functional parameters such as doping, stability, morphology, and charge-carrier mobility. The studies show that self-n-doping is best enabled by the use of sterically encumbered ammoniums with short tethers and Lewis basic counterions. Additionally, water is found to inhibit doping, which concludes that thermal degradation is merely a phenomenological feature of certain dopants, and that residual solvent evaporation is the primary driver of thermally activated doping. In situ grazing-incidence wide-angle X-ray scattering studies show that sample annealing increases the π–π stacking distance and shrinks grain boundaries for improved long-range ordering. These features are then correlated to contactless carrier-mobility measurements with time-resolved microwave conductivity before and after thermal annealing. In conclusion, the collective relationships between structural features and functionality are finally used to establish explicit self-n-dopant design principles for the future design of materials with improved functionality.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Exploring Phosphonium‐Based Anion Exchange Polymers for Moisture Swing Direct Air Capture of Carbon Dioxide

This study explores the performance and stability of ammonium and phosphonium-based polymeric ionic liquids (PILs) with methyl and butyl substituents in moisture-swing direct air capture of CO 2 . The polymers are synthesized with chloride counterions, followed by ion exchange to the bicarbonate ion, and tests for CO 2 capture capacity and stability under cyclic wet–dry conditions. The phosphonium polymer with methyl substituents [PVBT-MeP] demonstrates the highest CO 2 capture capacity at ≈510 µmol g⁻¹, attributed to minimal steric hindrance and stronger ion pairing with bicarbonate. However, oxidative degradation is detected by 31 P NMR spectroscopy after the moisture swing experiment, with the appearance of a phosphine oxide peak at 61.28 ppm, which indicates phosphorus oxidation as the primary degradation pathway. In contrast, the ammonium polymer with butyl substituents [PVBT-BuN] exhibits the highest stability, showing no degradation over five moisture swing cycles. Additional stability experiments in 0.5 m KHCO 3 solutions reveal no degradation for any PIL, suggesting that oxidative degradation is driven by dynamic acid-base reactions during the moisture swing cycles in the air. Furthermore, these findings reveal the potential of phosphonium-based PILs for moisture-swing direct air capture, achieving high capacity while highlighting the need for optimized stability through counterion and structural design.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Enhanced mobility of cations and anions in the redox state: The polaronium mechanism

Hugely enhanced slow-ion diffusivity has been widely observed under extreme redox conditions and for unclear reasons. Aided by first-principles calculations on model systems of ZrO 2 , CeO 2 , BaTiO 3 and Li 4/3 Mn 2/3 O 2 , here we successfully explained the intriguing phenomenon by a polaronium mechanism. We found a polaronium, defined as a transitory complex of a polaronic electron or hole and a migrating counterion, becomes highly mobile when the counterion comes from a d 0 or f 0 cation (e.g., Zr 4+ , Ti 4+ and Ce 4 ) or a p 6 anion (e.g., O 2- ) in the host compound. Upon a redox reaction, the complex attains a d 1 /f 1 or p 5 configuration, which spontaneously forms because it is favored by an electron-phonon interaction (manifest as the Jahn-Teller effect in high symmetry systems) that enables local relaxation and lowers the system energy. Here our calculations found such interaction reaching its peak at the saddle point where the local environment is softest, so soft that it allows a reorientation of the anisotropic d/f/p orbital to minimize the electron repulsion locally. Since the complex may dissolve after a successful ion-migration event, the redox electron/hole can be recycled to form another free-radical-like polaronium elsewhere, thereby enhancing ion migration repeatedly. The proposed polaronium mechanism, which also operates in ceramics doped with mixed-valence cations, is most relevant under dynamic and extreme thermal/field/irradiation conditions where extra electrons/holes are abundantly generated by non-equilibrium redox reactions. For such operations, some with emerging applications, our diffusion-enhancing mechanism may provide new theoretical insight to help understand their material/microstructure stability and performance.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Mechanism of ion pairing–induced co-ion penetration into the stern layer

In classical electric double layer (EDL) theory, co-ions—ions carrying the same charge as the electrode—are assumed to be excluded from the Stern layer because of electrostatic repulsion. However, ion pairing with adsorbed counterions may enable co-ions to access the EDL, a phenomenon that remains underexplored. Key questions therefore remain unresolved, including how deeply co-ions penetrate, how ion-specific properties govern them, and how ion pairing within the EDL differs from bulk solution. Here, in this study, we employ molecular dynamics simulations combined with well-tempered metadynamics to quantify co-ion penetration and ion-pairing thermodynamics at positively charged Mg-Al layered double hydroxide interfaces (σ = +20.6 μC/cm 2 ) in aqueous chloride electrolytes across 0.005–3 M range. We reveal a concentration-driven crossover: at dilute conditions, counterion–co-ion contact ion pairing (CIP) is thermodynamically unfavorable or only weakly favorable (ΔG CIP = +0.12 kcal/mol for Li + to −0.42 kcal/mol for K + ), leading to co-ion exclusion. As concentration increases (≥0.5 M), enhanced screening reduces surface repulsion and stabilizes CIP thermodynamics (ΔG CIP = −0.89 kcal/mol for Li + to −1.05 kcal/mol for K + at 3 M), enabling monovalent co-ions to penetrate the Stern layer following the hydration hierarchy K + > Na + > Li + , while divalent co-ions (Mg 2+ , Ca 2+ ) remain in the diffuse layer, forming only solvent-separated pairs. Comparison with bulk solution shows interfacial ion pairing is suppressed at low concentrations but converges to bulk-like thermodynamics at high concentrations, indicating surface effects vanish under strong screening. These results establish a quantitative framework linking ion-pairing thermodynamics to co-ion penetration and providing strategies to design electrochemical interfaces.

Co-ion penetration↗
DOE OSTI2025

Coarse-grained simulation of colloidal self-assembly, cation exchange, and rheology in Na/Ca smectite clay gels

Knowledge Gap: The aggregation of clay minerals—layered silicate nanoparticles—strongly impacts fluid flow, solute migration, and solid mechanics in soils, sediments, and sedimentary rocks. Experimental and computational characterization of clay aggregation is inhibited by the delicate water-mediated nature of clay colloidal interactions and by the range of spatial scales involved, from 1 nm thick platelets to flocs with dimensions up to micrometers or more. Simulations: Using a new coarse-grained molecular dynamics (CGMD) approach, we predicted the microstructure, dynamics, and rheology of hydrated smectite (more precisely, montmorillonite) clay gels containing up to 2,000 clay platelets on length scales up to 0.1 μm. Further, simulations investigated the impact of simulation time, platelet diameters (6 to 25nm), and the ratio of Na to Ca exchangeable cations on the assembly of tactoids (i.e., stacks of parallel clay platelets) and larger aggregates (i.e., assemblages of tactoids). We analyzed structural features including tactoid size and size distribution, basal spacing, counterion distribution in the electrical double layer, clay association modes, and the rheological properties of smectite gels. Findings: Our results demonstrate new potential to characterize and understand clay aggregation in dilute suspensions and gels on a scale of thousands of particles with explicit representation of counterion clouds and with accuracy approaching that of all-atom molecular dynamics (MD) simulations. For example, our simulations predict the strong impact of Na/Ca ratio on clay tactoid formation and the shear-thinning rheology of clay gels.

42 ENGINEERING↗