SEARCH · Engineering Papers

Results for “atomic dynamics”

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.

SourceResearch connectionFrom yearThrough year

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

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2025

A Finite Element Method for Compressible and Turbulent Multiphase Flow Instabilities with Heat Transfer

We present a new finite element framework for modeling compressible, turbulent multiphase flows with heat transfer. For two-fluid systems with a free surface, the Volume of Fluid (VOF) method is implemented without the need for interface reconstruction, while turbulence is resolved using a dynamic Vreman large eddy simulation (LES) model. Unlike most two-phase VOF studies, which neglect heat transfer, the present approach incorporates energy transport equations within the VOF formulation to account for heat exchange, an effect particularly important in turbulent flows. Conjugate heat transfer is often challenging in finite volume methods, which require explicit specification of heat fluxes at the solid–fluid interface, limiting accuracy and predictive capability. By contrast, the finite element formulation does not require heat flux inputs, allowing more accurate and robust simulation of heat transfer between solids and fluids. The method is demonstrated through three representative cases. First, a two-fluid instability with a single-mode perturbation is simulated and validated against analytical growth rates. Second, conjugate heat transfer is examined in a high-temperature flow over a cold metal cylinder, with validation performed both quantitatively—via pressure coefficient comparisons with experimental data—and qualitatively using vector field topology. Finally, compressible spray injection and breakup are modeled, demonstrating the ability of the framework to capture interfacial dynamics and atomization under turbulent, high-speed conditions. In the compressible spray injection and breakup case, the results indicate that the finite element formulation achieved higher predictive accuracy and robustness than the finite-volume method. With the same mesh resolution, the FEM reduced the root mean square error (RMSE) and mean absolute percentage error (MAPE) from 6.96 mm and 26.0% (for the FVM) to 4.85 mm and 12.7%, respectively, demonstrating improved accuracy and robustness in capturing interfacial dynamics and heat transfer. The study also introduced vector field topology to visualize and interpret coherent flow structures and instabilities, offering insights beyond conventional scalar-field analyses.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1988

Design and demonstration of a system for the deposition of atomic-oxygen durable coatings for reflective solar dynamic power system concentrators

A system for the vacuum deposition of atomic-oxygen durable coatings for reflective solar dynamic power systems (SDPS) concentrators was designed and demonstrated. The design issues pertinent to SDPS were developed by the Government Aerospace Systems Division of the Harris Corporation and are described in NASA-CR-179489. Both design and demonstration phases have been completed. At the time of this report the deposition system was ready for coating of facets for SDPS concentrators. The materials issue relevant to the coating work were not entirely resolved. These issues can only be resolved when substrates which are comparable to those which will be used in flight hardware are available. The substrates available during the contract period were deficient in the areas of surface roughness and contamination. These issues are discussed more thoroughly in the body of the report.

Mcclure, Donald J.↗

DOE OSTI2026

Quantifying Membrane Structure and Dynamics during Bioproduct Production in Zymomonas mobilis by Molecular Simulation

The conversion of lignocellulosic biomass into biofuels and bioproducts by microbial biorefineries is central to a sustainable chemical industry. Zymomonas mobilis is one such biorefinery chassis and is resistant to ethanol stress, leading to its use in biomass conversion to biofuels and bioproducts. However, Z. mobilis growth is often inhibited by organic acids, aldehydes, alcohols, ketones, and amides found in biomass hydrolysate. The resulting slow growth inhibits production and as a result drives up the price for the resulting products. One hypothesis is that these molecules interact with or disrupt the bacterial membrane, triggering stress responses and hindering growth. To test this hypothesis at the molecular level, we employ all-atom molecular dynamics (MD) simulations to investigate lignocellulose-derived small molecules and their impact on a biologically relevant Z. mobilis membrane model. Simulations were conducted across a range of inhibitor concentrations from 0 to 2.5 mol %, analyzing key membrane properties such as area per lipid (APL), membrane thickness, lipid-order parameter (−S CH ), lateral diffusion coefficient (D xy ), and permeability coefficient (Pm). From simulation, we observed altered membrane structure and dynamics at these modest small molecule concentrations commonly found in hydrolysates. Generally, the membranes become thinner, with a higher area per lipid and lower-order parameter as the small molecule concentration increases. These trends are stronger for more hydrophobic molecules with greater hydrophobic bulk, as isobutanol, propanol, and propanoic acid showed greater membrane perturbations as the concentration increased compared to other small molecules. Tracking small molecule distributions directly in our equilibrium simulations allows us to determine concentration-dependent free energy profiles for these molecules. While the trends are noisy, generally the barriers to crossing the membrane decrease as the concentration increases, indicating that the membranes become leakier as small molecule concentrations rise. Comparing between native Z. mobilis membranes with hopanoids and membranes sharing the same phospholipid composition but without hopanoids, hopanoids stabilize and order the membrane for smaller molecules to maintain membrane structure but appear insufficient for larger hydrophobic molecules like isobutanol. These findings provide a mechanistic understanding of how small molecules found in biomass degradation streams interact with the Z. mobilis membrane, offering valuable insights for future strain engineering efforts to optimize biofuel and bioproduct synthesis from biomass feedstocks by highlighting limits to small molecule tolerance. This knowledge can guide the modification of membrane composition to develop more robust microbes, thereby improving microbial survival and yields in industrial contexts.

Singh, Nitin Kumar [Michigan State Univ., East Lan↗

DOE OSTI2022

Hierarchical Self-Assembly Pathways of Peptoid Helices and Sheets

Peptoids (N-substituted glycines) are a class of tailorable synthetic peptidomic polymers. Amphiphilic diblock peptoids have been engineered to assemble 2D crystalline lattices with applications in catalysis and molecular separations. Assembly is induced in an organic solvent/water mixture by evaporating the organic phase, but the assembly pathways remain uncharacterized. We conduct all-atom molecular dynamics simulations of Nbrpe6Nc6 as a prototypical amphiphilic diblock peptoid comprising an NH 2 -capped block of six hydrophobic N-((4-bromophenyl)ethyl)glycine residues conjugated to a polar NH 3 (CH 2 ) 5 CO tail. We identify a thermodynamically controlled assembly mechanism by which monomers assemble into disordered aggregates that self-order into 1D chiral helical rods then 2D achiral crystalline sheets. We support our computational predictions with experimental observations of 1D rods using small-angle X-ray scattering, circular dichroism, and atomic force microscopy and 2D crystalline sheets using X-ray diffraction and atomic force microscopy. Furthermore, this work establishes a new understanding of hierarchical peptoid assembly and principles for the design of peptoid-based nanomaterials.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Magnetic domains and domain wall pinning in atomically thin CrBr 3 revealed by nanoscale imaging

The emergence of atomically thin van der Waals magnets provides a new platform for the studies of two-dimensional magnetism and its applications. However, the widely used measurement methods in recent studies cannot provide quantitative information of the magnetization nor achieve nanoscale spatial resolution. These capabilities are essential to explore the rich properties of magnetic domains and spin textures. Here, we employ cryogenic scanning magnetometry using a single-electron spin of a nitrogen-vacancy center in a diamond probe to unambiguously prove the existence of magnetic domains and study their dynamics in atomically thin CrBr 3 . By controlling the magnetic domain evolution as a function of magnetic field, we find that the pinning effect is a dominant coercivity mechanism and determine the magnetization of a CrBr 3 bilayer to be about 26 Bohr magnetons per square nanometer. The high spatial resolution of this technique enables imaging of magnetic domains and allows to locate the sites of defects that pin the domain walls and nucleate the reverse domains. Our work highlights scanning nitrogen-vacancy center magnetometry as a quantitative probe to explore nanoscale features in two-dimensional magnets.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Oxidation in laser-generated metal plumes

The temporal evolution of atoms and molecules in a laser-produced plasma was investigated using optical emission spectroscopy for several metal targets (i.e., Al, Ti, Fe, Zr, Nb, and Ta). Plasmas from metal targets were generated by focusing 1064 nm, 6 ns pulses from an Nd:YAG laser. Gas-phase oxidation/plasma chemistry was initiated by adding O2 (partial pressures up to ≈20%) to an N2 environment where the total background pressure was kept at a constant 1 atmosphere. Temporally resolved emission spectral features were used to track the gas-phase oxidation. The dynamics of atomic and molecular species were monitored using space-resolved time-of-flight emission spectroscopy. Our results highlight that the partial pressure of O2 strongly influences spectral features and molecular formation in laser-produced plasmas. Atoms and molecules co-exist in plasmas, although with different temporal histories depending on the target material due to differences in thermo- and plasma chemical reactions occurring in the plume.

Kautz, Elizabeth J. (ORCID:0000000263389223)↗

DOE OSTI2021

Phonon-Assisted Intervalley Scattering Determines Ultrafast Exciton Dynamics in MoSe 2 Bilayers

While valleys (energy extrema) are present in all band structures of solids, their preeminent role in determining exciton resonances and dynamics in atomically thin transition metal dichalcogenides (TMDC) is unique. Using two-dimensional coherent electronic spectroscopy, we find that exciton decoherence occurs on a much faster timescale in MoSe 2 bilayers than that in the monolayers. We further identify two population relaxation channels in the bilayer, a coherent and an incoherent one. Our microscopic model reveals that phonon-emission processes facilitate scattering events from the K valley to other lower-energy Γ and Λ valleys in the bilayer. Our combined experimental and theoretical studies unequivocally establish different microscopic mechanisms that determine exciton quantum dynamics in TMDC monolayers and bilayers. Understanding exciton quantum dynamics provides critical guidance to the manipulation of spin-valley degrees of freedom in TMDC bilayers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Structural insights into the modulation of coronavirus spike tilting and infectivity by hinge glycans

Abstract Coronavirus spike glycoproteins presented on the virion surface mediate receptor binding, and membrane fusion during virus entry and constitute the primary target for vaccine and drug development. How the structure dynamics of the full-length spikes incorporated in viral lipid envelope correlates with the virus infectivity remains poorly understood. Here we present structures and distributions of native spike conformations on vitrified human coronavirus NL63 (HCoV-NL63) virions without chemical fixation by cryogenic electron tomography (cryoET) and subtomogram averaging, along with site-specific glycan composition and occupancy determined by mass spectrometry. The higher oligomannose glycan shield on HCoV-NL63 spikes than on SARS-CoV-2 spikes correlates with stronger immune evasion of HCoV-NL63. Incorporation of cryoET-derived native spike conformations into all-atom molecular dynamic simulations elucidate the conformational landscape of the glycosylated, full-length spike that reveals a role of hinge glycans in modulating spike bending. We show that glycosylation at N1242 at the upper portion of the stalk is responsible for the extensive orientational freedom of the spike crown. Subsequent infectivity assays implicated involvement of N1242-glyan in virus entry. Our results suggest a potential therapeutic target site for HCoV-NL63.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Charge transfer rate constants for the carotenoid-porphyrin-C 60 molecular triad dissolved in tetrahydrofuran: The spin-boson model vs the linearized semiclassical approximation

Charge transfer rate constants were calculated for the carotenoid-porphyrin-C 60 (CPC 60 ) molecular triad dissolved in explicit tetrahydrofuran. The calculation was based on mapping the all-atom anharmonic Hamiltonian of this system onto the spin-boson Hamiltonian. The mapping was based on discretizing the spectral density from the time correlation function of the donor–acceptor potential energy gap, as obtained from all-atom molecular dynamics simulations. Different spin-boson Hamiltonians were constructed for each of the possible transitions between the three excited electronic states in two different triad conformations. The rate constants of three possible transitions were calculated via the quantum-mechanically exact Fermi’s golden rule (FGR), as well as a progression of more approximate expressions that lead to the classical Marcus expression. The advantage of the spin-boson approach is that once the mapping is established, the quantum-mechanically exact FGR and the hierarchy of approximations are known in closed form. The classical Marcus charge transfer rate constants obtained with the spin-boson Hamiltonians were found to reproduce those obtained from all-atom simulations with the linearized semiclassical approximation, thereby confirming the equivalence of the two approaches for this system. Furthermore, within the spin-boson Hamiltonian, we also found that the quantum-mechanically exact FGR rate constants were significantly enhanced compared to the classical Marcus theory rate constants for two out of three transitions in one of the two conformations under consideration. The results confirm that mapping to the spin-boson model can yield accurate predictions for charge transfer rate constants in a system as complex as CPC 60 dissolved in tetrahydrofuran.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Electronic Structure of the Au Nanoparticle-TiO 2 Heterojunction: Influence of Nanoparticle Size, Shape, Oxygen Vacancies, and Temperature

The electronic structure of the gold nanoparticle-titanium dioxide (AuNP-TiO 2 ) heterojunction plays a critical role for charge transfer and recombination dynamics that underpin its photocatalytic function. However, building a representative model to capture the key physics remains a significant challenge. Here, we investigate the influence of the AuNP size and shape, as well as oxygen vacancy (V O ) defects at the anatase-phase TiO 2 (101) surface and the temperature of the heterojunction, on its interfacial electronic properties. Using density functional theory (DFT), we compare the closed-shell Au 20 and open-shell Au 19 clusters interfaced with pristine and V O defect TiO 2 surfaces. We find that the presence of a V O defect transforms pure TiO 2 from a p-type to an n-type semiconductor, reversing the interfacial band bending from downward to upward. For the heterosystem, density of states (DOS) analysis shows that V O minimally affects the open-shell Au 19 -TiO 2 system, but it significantly alters the closed-shell Au 20 -TiO 2 system, converting it from p-type to n-type at the Γ-point at 0 K. Furthermore, ab initio molecular dynamics (AIMD) simulations at 300 K reveal significant thermal fluctuations in AuNP positions relative to the TiO 2 surface. These fluctuations result in dynamic variations in the gap between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) in all systems studied: Au 19 NP-Pristine TiO 2 , Au 20 NP-Pristine TiO 2 , Au 19 NP-V O defect TiO 2 , and Au 20 NP-V O defect TiO 2 . Complementary Bader charge analysis performed at both 0 K and finite-temperature AIMD snapshots supports the emergence of an upward band bending and the formation of a Schottky barrier at the V O containing heterojunctions. Notably, we find that at finite temperature, an Au atom can dynamically passivate the V O , leading to a pronounced widening of the HOMO-LUMO gap in the AuNP-V O defect TiO 2 heterojunctions. Our computational findings underscore the pivotal role of the V O defect and thermal effects in modulating interfacial band alignment, electronic states, and HOMO-LUMO gaps, providing insights for designing Au-TiO 2 heterojunctions with tailored electronic properties.

Mora Perez, Carlos [Lawrence Berkeley National Lab↗

DOE OSTI2024

Noncentrosymmetric Triangular Magnet CaMnTeO 6 : Strong Quantum Fluctuations and Role of s 0 versus s 2 Electronic States in Competing Exchange Interactions

Abstract Noncentrosymmetric triangular magnets offer a unique platform for realizing strong quantum fluctuations. However, designing these quantum materials remains an open challenge attributable to a knowledge gap in the tunability of competing exchange interactions at the atomic level. Here, a new noncentrosymmetric triangular S = 3/2 magnet CaMnTeO 6 is created based on careful chemical and physical considerations. The model material displays competing magnetic interactions and features nonlinear optical responses with the capability of generating coherent photons. The incommensurate magnetic ground state of CaMnTeO 6 with an unusually large spin rotation angle of 127°(1) indicates that the anisotropic interlayer exchange is strong and competing with the isotropic interlayer Heisenberg interaction. The moment of 1.39(1) µB , extracted from low‐temperature heat capacity and neutron diffraction measurements, is only 46% of the expected value of the static moment 3 µB . This reduction indicates the presence of strong quantum fluctuations in the half‐integer spin S = 3/2 CaMnTeO 6 magnet, which is rare. By comparing the spin‐polarized band structure, chemical bonding, and physical properties of AMnTeO 6 ( A = Ca, Sr, Pb), how quantum‐chemical interpretation can illuminate insights into the fundamentals of magnetic exchange interactions, providing a powerful tool for modulating spin dynamics with atomically precise control is demonstrated.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Dynamics and lipid membrane coupling of the RAS-RAF complex revealed via multiscale simulations

To gain molecular and mechanistic insights into initiation of the RAS-RAF signaling cascade, we developed and used a combination of multiscale simulation and experimental approaches. The influence and impact of the membrane on RAS and RAF proteins is a factor we are just beginning to understand and appreciate in more detail. Molecular simulation is an ideal methodology to further study this complicated relationship between the membrane and associated proteins. Our previous work using Multiscale Machine-learned Modeling Infrastructure investigated different lipid compositions solely around the KRAS4b protein and the interplay between protein behavior and these membrane environments. Multiscale Machine-learned Modeling Infrastructure uses machine learning to couple adjacent simulation scales and has been efficiently scaled across some of the world’s largest high-performance computers. Recently, we have expanded this multiresolution framework to include the all-atom simulation scale and to incorporate the RAF RBDCRD domains. Here, we present the overall analysis results from this new simulation campaign comprising a mixture of RAS and RAF RBDCRD proteins. Approximately 35,000 coarse-grained and 10,000 all-atom molecular dynamics simulations were completed, sampled from a variety of protein/lipid composition configurations that were generated from a micron-scale continuum simulation containing hundreds of copies of the proteins. Our studies suggest that orientations of the RAS-RBDCRD complex on the membrane occupy distinct configurational states, and the spatial patterns of lipid arrangements around these different protein states are unique to each state. The extent and size of lipid “fingerprints” imposed on the membrane by the RAS-RBDCRD protein complex are significantly larger than observed for just the RAS protein on its own. These protein complexes strongly associate, but we do not observe statistically significant preferred protein-protein orientations. These observations indicate that spatial colocalization of RAS-RBDCRD proteins in the same vicinity may be assisted by specific membrane environments, acting to increase the probability of signaling complex formation.

Carpenter, Timothy S. [Lawrence Livermore National↗

DOE OSTI2021

Effect of surface steps on chemical ordering in the subsurface of Cu(Au) solid solutions

Atomic steps are typically assumed to only influence surface phenomena of solids. In contrast, in this paper we show, using in-situ atomic-resolution electron microscopy observations and atomistic modeling, the pronounced effect from surface stepsin inducing compositional and structural evolution in the subsurface of a Cu(Au) solid solution. We find that Au surface segregation results in a stacked sequence of Cu-Au ordered phases in the subsurface. The presence of a monatomic step at the surface induces the formation of an anti-phase boundary that extends from the surface step to deeper layers by maintaining the same composition profile associated with each terrace. The bunching of surface steps induces chemical disordering of the Cu-Au ordered phases in the subsurface region of the bunched steps. These results demonstrate the instant propagation of surface dynamics of atomic steps into the subsurface region and can find broader applicability in utilizing surface defects to tune the composition and structure in the subsurface of alloys.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Macroscale properties and atomic-scale mechanisms of ash removal in low-temperature hydrothermal carbonization

Biogenic ash is a significant impediment to the utilization of agricultural residues in biofuel production. Such challenge can be addressed by various treatments, as demonstrated in this study on the experimental and computational mechanisms involved in the hydrothermal treatment (HT) of wheat straw. A combination of classical (all-atom) molecular dynamics simulations of cellulose carrying silica and calcium species, along with first principles quantum chemical calculations, indicates the dissociation of inorganics from the cellulose with increased HT temperature. Here, this observation is confirmed by experimental evidence of effective ash removal by HT, showing at least 50% removal of sulfur, chlorine, potassium, and calcium, and 12.5% of silica, leading to a reduced total ash content (from 6.7% to 4.2%).Changes in structural features upon HT, such as surface cellular structure and porosity, were revealed, accompanied by an increased specific surface area (from 1.17 to 6.34 m 2 /g). Our simulations suggest that silica binds tightly to the hydrophobic face of cellulose at room temperature, but HT significantly reduces the binding free energy of association with both hydrophobic and hydrophilic surfaces. Most significantly, ash removal leads to an increased calorific value, rising from approximately 16 MJ/kg to about 19 MJ/kg, along with improved thermal behavior. The improved integration combustion index parameter S indicates that the combustion properties improve with ash removal efficiency. The proposed atomic-level mechanism for the observed removal of inorganics during mild HT underscores the potential of such treatment in producing energy-dense wheat straw, a widely available agricultural residue.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Neural-network quantum states for the nuclear many-body problem

A long-standing goal of nuclear theory is to explain how the structure and dynamics of atomic nuclei and neutron-star matter emerge from the underlying interactions among protons and neutrons. Achieving this goal requires solving the nuclear quantum many-body problem with high accuracy across a wide range of length scales and density regimes. In this review, we discuss how artificial neural network representations of the nuclear many-body wave function have significantly extended the capabilities of continuum quantum Monte Carlo methods. In particular, neural network quantum states enable calculations of larger systems than were previously accessible and provide a flexible framework for capturing phenomena that challenge conventional approaches, including the emergence of nuclear clusters and superfluid phases in dense matter. We highlight recent applications to finite nuclei, infinite nuclear and neutron matter, and dynamical processes relevant to lepton-nucleus and nucleus-nucleus scattering. We also discuss conceptual and methodological connections with condensed matter physics, emphasizing developments in neural network quantum states that bridge strongly correlated systems across disciplines. Together, these developments demonstrate how neural-network methods open new avenues toward unified and accurate descriptions of nuclear structure, matter, and reactions.

Lovato, Alessandro [Argonne; TIFPA-INFN, Trento; V↗

NASA NTRS2003

Diagnostics of Carbon Nanotube Formation in a Laser Produced Plume: An Investigation of the Metal Catalyst by Laser Ablation Atomic Fluorescence Spectroscopy

Carbon nanotubes, elongated molecular tubes with diameters of nanometers and lengths in microns, hold great promise for material science. Hopes for super strong light-weight material to be used in spacecraft design is the driving force behind nanotube work at JSC. The molecular nature of these materials requires the appropriate tools for investigation of their structure, properties, and formation. The mechanism of nanotube formation is of particular interest because it may hold keys to controlling the formation of different types of nanotubes and allow them to be produced in much greater quantities at less cost than is currently available. This summer's work involved the interpretation of data taken last summer and analyzed over the academic year. The work involved diagnostic studies of carbon nanotube formation processes occurring in a laser-produced plume. Laser ablation of metal doped graphite to produce a plasma plume in which carbon nanotubes self assemble is one method of making carbon nanotube. The laser ablation method is amenable to applying the techniques of laser spectroscopy, a powerful tool for probing the energies and dynamics of atomic and molecular species. The experimental work performed last summer involved probing one of the metal catalysts, nickel, by laser induced fluorescence. The nickel atom was studied as a function of oven temperature, probe laser wavelength, time after ablation, and position in the laser produced plume. This data along with previously obtained data on carbon was analyzed over the academic year. Interpretations of the data were developed this summer along with discussions of future work. The temperature of the oven in which the target is ablated greatly influences the amount of material ablated and the propagation of the plume. The ablation conditions and the time scale of atomic and molecular lifetimes suggest that initial ablation of the metal doped carbon target results in atomic and small molecular species. The metal atoms survive for several milliseconds while the gaseous carbon atoms and small molecules nucleate more rapidly. Additional experiments and the development of in situ methods for carbon nanotube detection would allow these results to be interpreted from the perspective of carbon nanotube formation.

deBoer, Gary↗

NASA NTRS2006

Hybrid MD-Nernst Planck Model of Alpha-hemolysin Conductance Properties

Motivated by experiments in which an applied electric field translocates polynucleotides through an alpha-hemolysin protein channel causing ionic current transient blockade, a hybrid simulation model is proposed to predict the conductance properties of the open channel. Time scales corresponding to ion permeation processes are reached using the Poisson-Nemst-Planck (PNP) electro-diffusion model in which both solvent and local ion concentrations are represented as a continuum. The diffusion coefficients of the ions (K(+) and Cl(-)) input in the PNP model are, however, calculated from all-atom molecular dynamics (MD). In the MD simulations, a reduced representation of the channel is used. The channel is solvated in a 1 M KCI solution, and an external electric field is applied. The pore specific diffusion coefficients for both ionic species are reduced 5-7 times in comparison to bulk values. Significant statistical variations (17-45%) of the pore-ions diffusivities are observed. Within the statistics, the ionic diffusivities remain invariable for a range of external applied voltages between 30 and 240mV. In the 2D-PNP calculations, the pore stem is approximated by a smooth cylinder of radius approx. 9A with two constriction blocks where the radius is reduced to approx. 6A. The electrostatic potential includes the contribution from the atomistic charges. The MD-PNP model shows that the atomic charges are responsible for the rectifying behaviour and for the slight anion selectivity of the a-hemolysin pore. Independent of the hierarchy between the anion and cation diffusivities, the anionic contribution to the total ionic current will dominate. The predictions of the MD-PNP model are in good agreement with experimental data and give confidence in the present approach of bridging time scales by combining a microscopic and macroscopic model.

Cozmuta, Ioana↗

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↗