Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atomistic simulation”

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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Composition and Configuration Dependence of Glass-Transition Temperature in Binary Copolymers and Blends of Polyhydroxyalkanoate Biopolymers

Polyhydroxyalkanoates (PHAs), a promising class of biomaterials, have gained considerable attention to replace petroleum-based plastics owing to their excellent biocompatibility and biodegradability. Homopolymers of PHA suffer from poor tunability in thermal and mechanical properties. Going from homopolymers to copolymers, the design space can be substantially enhanced by combining two or more monomers in different compositions (i.e., relative ratios of the different monomers) and configurations (i.e., relative positions of the different monomers in the polymer backbone) leading to a substantially large chemical space where application-specific optimization for the targeted functionality can be performed. However, this composition and configuration dependence of properties in the vast PHA copolymer chemical space remains largely unexplored. In this contribution, further building on our past work with PHA homopolymers, we systematically explore these chemical trends for glass-transition temperature (Tg) in PHA copolymers and blends. Overall, our molecular dynamics simulations, utilizing a previously validated force field for PHAs, suggest that these trends are largely governed not only by the homopolymer Tg values but also configuration-dependent interchain interactions in the copolymer system. In particular, our results indicate that the configuration-dependent variation in the target property at a fixed composition can be significant in the presence of hydrogen-bond-forming monomers. These qualitative observations are further rationalized by quantitatively analyzing various closely related atomic level descriptors of copolymers and blends such as monomer mobility, number of hydrogen bonds, and pair correlation functions. The findings presented in this work help to develop a deeper atomistic-level understanding of thermomechanical behavior of PHA-based copolymers and can potentially guide the rational design of biopolymers as environmental-friendly functional materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2009

A Continuum-Atomistic Analysis of Transgranular Crack Propagation in Aluminum

A concurrent multiscale modeling methodology that embeds a molecular dynamics (MD) region within a finite element (FEM) domain is used to study plastic processes at a crack tip in a single crystal of aluminum. The case of mode I loading is studied. A transition from deformation twinning to full dislocation emission from the crack tip is found when the crack plane is rotated around the [111] crystallographic axis. When the crack plane normal coincides with the [112] twinning direction, the crack propagates through a twinning mechanism. When the crack plane normal coincides with the [011] slip direction, the crack propagates through the emission of full dislocations. In intermediate orientations, a transition from full dislocation emission to twinning is found to occur with an increase in the stress intensity at the crack tip. This finding confirms the suggestion that the very high strain rates, inherently present in MD simulations, which produce higher stress intensities at the crack tip, over-predict the tendency for deformation twinning compared to experiments. The present study, therefore, aims to develop a more realistic and accurate predictive modeling of fracture processes.

Yamakov, V.↗
DOE OSTI2025

Universal Electronic‐Structure Relationship Governing Intrinsic Magnetic Properties in Permanent Magnets

An electronic-structure-centered perspective is presented on permanent-magnet (PM) design, highlighting two key levers, that is, saturation magnetization (M s ), governed by 3d-band filling and exchange physics, and magnetocrystalline anisotropy energy (MAE), arising from spin-orbit coupling (SOC) on anisotropic orbital populations. Reviewing current practices, including DFT-based MAE/J ij extraction, atomistic-spin and micromagnetic modeling, and high-throughput machine learning (ML) pipelines, three bottlenecks limiting predictive discovery is identified that is i) electronic-structure accuracy for small MAE (sensitive to functional choice, Hubbard U, and many-body effects), ii) finite-temperature and kinetic realism (phonon/magnon renormalization, ordering kinetics), and iii) descriptor and multiscale decoupling (lack of SOC-weighted and orbital-resolved fingerprints). Deep dives into the electronic-structure of Nd─Fe─B and Fe─N show how these fingerprints govern magnetic performance, motivating DFT- and quantum-mechanics-based descriptors for discovery. Unbiased, structure-driven exploration, coupled with high-throughput simulations, ML, generative AI, and reasoning models, accelerates candidate identification and propagates insights across scales. Addressing supply-chain risks, on future needs of designing “critical-element-free” magnets with tailored microstructure and high energy products is emphasized. By integrating electronic fingerprints, AI reasoning, and multiscale modeling, a practical roadmap is provided for rare-earth-lean or rare-earth free, high-performance, sustainable PMs.

Singh, Prashant [Ames Laboratory, and Iowa State U↗
DOE OSTI2022

Data-driven electron-diffraction approach reveals local short-range ordering in CrCoNi with ordering effects

Abstract The exceptional mechanical strength of medium/high-entropy alloys has been attributed to hardening in random solid solutions. Here, we evidence non-random chemical mixing in a CrCoNi alloy, resulting from short-range ordering. A data-mining approach of electron nanodiffraction enabled the study, which is assisted by neutron scattering, atom probe tomography, and diffraction simulation using first-principles theory models. Two samples, one homogenized and one heat-treated, are observed. In both samples, results reveal two types of short-range-order inside nanoclusters that minimize the Cr–Cr nearest neighbors (L1 2 ) or segregate Cr on alternating close-packed planes (L1 1 ). The L1 1 is predominant in the homogenized sample, while the L1 2 formation is promoted by heat-treatment, with the latter being accompanied by a dramatic change in dislocation-slip behavior. These findings uncover short-range order and the resulted chemical heterogeneities behind the mechanical strength in CrCoNi, providing general opportunities for atomistic-structure study in concentrated alloys for the design of strong and ductile materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Status Report on Design of In-situ Thermomechanical Testing at LANSCE

Nuclear fuel encounters severe thermomechanical environments in which its mechanical response is determined by its microstructure, temperature and stress level histories. Simulating the response of such microstructures is crucial for predicting both performance and transient fuel mechanical responses and experimental verification of such predictions is therefore of great interest. While most of the deformation in a nuclear fuel rod occurs in the cladding, deformation of the fuel itself is still of interest with deformation mechanisms at operating temperature and above including creep, swelling, cracking as well as pellet-clad interaction. Characterization of these properties and understanding of the underlying deformation phenomena at operating or excursion temperatures is therefore of great importance for development and ultimately licensing of improved and novel nuclear fuel forms. Diffraction techniques offer unique insight on the atomistic (e.g. crystal structure) and microstructure (e.g. phase transformations, texture, defects) length scales and have a long history of providing unique data to inform relevant deformation models that enable the required predictive capabilities. For example, dislocations lead to diffraction peak broadening that can be characterized to estimate the dislocation density and study the role of dislocations on the deformation while measuring lattice strains allows to studie load sharing in two phase materials. In this report the requirements for a sample environment for high temperature deformation of nuclear fuels are defined. The HIPPO neutron time-of-flight diffractometer at LANSCE will host this sample environment and is also described. This instrument covers diffraction angles from 140° to 40° and is also equipped with an event-mode neutron imaging detector system, enabling energy-resolved neutron imaging in parallel with the diffraction that could measure sample temperature from Doppler broadening of neutron absorption resonances or measure pore densities from changes in the attenuation. Designs of devices to characterize thermomechanical properties of nuclear fuel without diffraction are also considered to guide the design. While this report is focused on applications for nuclear fuels, the device can also characterize cladding, moderator or structural materials and therefore contribute to other fields of research and development for advanced reactors. The temperatures planned to be reached are above 2000℃, thus enabling characterization of LWR reactor fuels under accident scenarios but also reaching temperatures of fuels developed for nuclear thermal propulsion and providing opportunities to characterize those. In conjunction with the energy-resolved neutron imaging detector, this setup would allow to measure neutron cross-sections at high temperatures, filling a gap towards development of reactors operating at high temperatures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Accelerating multiscale electronic stopping power predictions with time-dependent density functional theory and machine learning

Knowing the rate at which particle radiation releases energy in a material, the “stopping power,” is key to designing nuclear reactors, medical treatments, semiconductor and quantum materials, and many other technologies. While the nuclear contribution to stopping power, i.e., elastic scattering between atoms, is well understood in the literature, the route for gathering data on the electronic contribution has for decades remained costly and reliant on many simplifying assumptions, including that materials are isotropic. We establish a method that combines time-dependent density functional theory (TDDFT) and machine learning to reduce the time to assess new materials to hours on a supercomputer and provide valuable data on how atomic details influence electronic stopping. Our approach uses TDDFT to compute the electronic stopping from first principles in several directions and then machine learning to interpolate to other directions at a cost of 10 million times fewer core-hours. We demonstrate the combined approach in a study of proton irradiation in aluminum and employ it to predict how the depth of maximum energy deposition, the “Bragg Peak,” varies depending on the incident angle—a quantity otherwise inaccessible to modelers and far outside the scales of quantum mechanical simulations. The lack of any experimental information requirement makes our method applicable to most materials, and its speed makes it a prime candidate for enabling quantum-to-continuum models of radiation damage. The prospect of reusing valuable TDDFT data for training the model makes our approach appealing for applications in the age of materials data science.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Developing reliable machine learning interatomic potential for Fe–Cr–Ni austenitic alloys

Gaining atomistic understanding of mechanical behavior of heat-resistant structural materials such as Fe–Cr–Ni-based alloys requires an approach with an accuracy close to density functional theory (DFT) that considers the intrinsic properties of the bulk lattice and important defects such as stacking faults, grain boundaries, and surfaces. This work aims to develop reliable machine learning interatomic potential (MLIAP) at cross-scale for Fe–Cr–Ni ternary alloys with a focus on the face-centered-cubic (fcc) solid solution structure. Leveraging the advantages of moment tensor potentials, which typically necessitate a relatively small training dataset and enable rapid calculations using the large-scale atomic/molecular massively parallel simulator package, we ensure the stability and accuracy of the trained potentials. Important defects such as stacking faults, grain boundaries, and surfaces for wide-range compositions are investigated. Structural, thermal, elastic, and defect properties are determined from molecular dynamics simulations comprising several thousand atoms, generated via canonical Monte Carlo simulations guided by the trained potential. The trained potential allows efficient atomic simulations of structural, thermal, and mechanical properties of fcc Fe–Cr–Ni solid solution alloys as a function of composition and temperature. Therefore, the MLIAP approach represents a major advancement from DFT calculations that are limited to small simulation sizes and traditional molecular dynamics simulations using relatively low accuracy potentials. Furthermore, this work outlines a practical foundation for further investigating the structural evolution and mechanical behavior of austenitic stainless steel and nickel-based alloys in a wide array of applications in extreme environments.

Crystal structure↗
DOE OSTI2024

A transferable classical force field to describe glyme based lithium solvate ionic liquids

A non-polarizable force field for lithium (Li + ) and bis(trifluoromethanesulfonyl)imide (TFSI – ) ions solvated in diglyme at around 0.2 mol fraction salt concentration was developed based on ab initio molecular dynamics (AIMD) simulations and a modified polymer consistent force field model. A force–torque matching based scheme, in conjunction with a genetic algorithm, was used to determine the Lennard-Jones (LJ) parameters of the ion–ion and ion–solvent interactions. This force field includes a partial charge scaling factor and a scaling factor for the 1–4 interactions. The resulting force field successfully reproduces the radial distribution function of the AIMD simulations and shows better agreement compared to the unmodified force field. The new force field was then used to simulate salt solutions with glymes of increasing chain lengths and different salt concentrations. The comparison of the MD simulations, using the new force field, with experimental data at different salt concentrations and AIMD simulations on equimolar concentrations of the triglyme system demonstrates the transferability of the force field parameters to longer glymes and higher salt concentrations. Furthermore, the force field appears to reproduce the features of the experimental x-ray structure factors, suggesting accuracy beyond the first solvation shell, for equimolar salt solutions using both triglyme and tetraglyme as the solvent. Altogether, the new force field was found to accurately reproduce the molecular descriptions of LiTFSI-glyme systems not only at various salt concentrations but also with glymes of different chain lengths. Thus, the new force field provides a useful and accurate tool to perform in silico studies of this family of systems at the atomistic level.

25 ENERGY STORAGE↗
DOE OSTI2023

Initial decomposition mechanisms of 2,4,6-triamino-1,3,5-trinitrobenzene (TATB) and their kinetic isotope effect

2,4,6-triamino-1,3,5-trinitrobenzene (TATB) is an insensitive High Explosive (HE) that is widely studied to better understand the physical properties of safety and sensitivity of HE. A dominant initial decomposition mechanism of TATB is believed to be a dehydration reaction that forms mono- and di-furazans, although other mechanisms have been reported. In this work, seven initial decomposition mechanisms were modeled with ab initio simulations to calculate its free energy barriers, decomposition rates, and kinetic isotope effects. The energy barrier for mono-benzofurazan mechanisms was found to be high, >61 kcal/mol in the gas phase; however, the reaction energy can decrease significantly in a disordered condensed state. The predicted kinetic isotope effect ratio of the furazan mechanism was found to be k h /k d ≈ 1.41 at 60K, in agreement with the experiment. The NO 2 scission mechanism was found to be an entropy-driven mechanism because the free energy barrier decreased significantly with temperature, making it the most energetically favorable mechanism at high temperatures in the gas phase. So the results provide a better understanding of the atomistic decomposition mechanisms of TATB and may be useful for improving models of safety and sensitivity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2022

Discovery of multimechanisms of screw dislocation interaction in bcc iron from open-ended saddle point searches

Dislocation motion and interactions determine mechanical properties in body-centered cubic (bcc) metallic materials. However, studying mechanisms for the screw dislocation interaction is fundamentally challenging since many underlying processes involve mesotimescales and atomistic resolution, currently inaccessible by either experimental techniques or continuum theoretical methods. Here, in this paper, we develop a computational capability based on self-evolving atomistic kinetic Monte Carlo (SEAKMC) to sample the critical events and saddle point energies related to screw dislocations and their junctions. The method is first validated by calculating the stress dependence of Peierls barriers and formation energies of kink pairs and cross-slip kink pairs on a single screw dislocation in bcc iron. Then the method is applied to a binary junction of a pair of intersecting screw dislocations, the structure of which is crucial for low-temperature plastic deformation. We identify three important mechanisms: coplanar cross-slipping, jog-pinning, and a previously unknown unzipping mechanism during the evolution of the binary junction. The mechanisms are then further validated using classical molecular dynamics simulations. The computational capability developed in this paper provides an effective tool to evaluate screw dislocation related thermally activated events in complex stress conditions. The mechanisms discovered in this paper provide critical insights into temperature dependence of the anomalous slip, a breakdown of the Schmidt law, during the plastic deformation in bcc iron and can be generalized to other bcc metals.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Artificial neural network potentials for mechanics and fracture dynamics of two-dimensional crystals **

Understanding the mechanics and failure of materials at the nanoscale is critical for their engineering and applications. The accurate atomistic modeling of brittle failure with crack propagation in covalent crystals requires a quantum mechanics-based description of individual bond-breaking events. Artificial neural network potentials (NNPs) have emerged to overcome the traditional, physics-based modeling tradeoff between accuracy and accessible time and length scales. Previous studies have shown successful applications of NNPs for describing the structure and dynamics of molecular systems and amorphous or liquid phases of materials. However, their application to deformation and failure processes in materials is still uncommon. In this study, we discuss the apparent limitations of NNPs for the description of deformation and fracture under loadings and propose a way to generate and select training data for their employment in simulations of deformation and fracture simulations of crystals. We applied the proposed approach to 2D crystalline graphene, utilizing the density-functional tight-binding method for more efficient and extensive data generation in place of density functional theory. Then, we explored how the data selection affects the accuracy of the developed artificial NNPs. It revealed that NNP’s reliability should not only be measured based on the total energy and atomic force comparisons for reference structures but also utilize comparisons for physical properties, e.g. stress–strain curves and geometric deformation. In sharp contrast to popular reactive bond order potentials, our optimized NNP predicts straight crack propagation in graphene along both armchair and zigzag (ZZ) lattice directions, as well as higher fracture toughness of ZZ edge direction. Our study provides significant insight into crack propagation mechanisms on atomic scales and highlights strategies for NNP developments of broader materials.

2D materials↗
DOE OSTI2023

Energy localization efficiency in 1,3,5-trinitro-2,4,6-triaminobenzene pore collapse mechanisms

Atomistic and continuum scale modeling efforts have shown that the shock-induced collapse of porosity can occur via a wide range of mechanisms dependent on pore morphology, the shockwave pressure, and material properties. The mechanisms that occur under weaker shocks tend to be more efficient at localizing thermal energy but do not result in high, absolute temperatures or spatially large localizations compared to mechanisms found under strong shock conditions. However, the energetic material 1,3,5-trinitro-2,4,6-triaminobenzene (TATB) undergoes a wide range of collapse mechanisms that are not typical of similar materials, leaving the collapse mechanisms and the resultant energy localization from the collapse, i.e., hotspots, relatively uncharacterized. Therefore, we present the pore collapse simulations of cylindrical pores in TATB for a wide range of pore sizes and shock strengths that trigger viscoplastic collapses that occur almost entirely perpendicular to the shock direction for weak shocks and hydrodynamic-like collapses for strong shocks that do not break the strong hydrogen bonds of the TATB basal planes. The resulting hotspot temperature fields from these mechanisms follow trends that differ considerably from other energetic materials; hence, we compare them under normalized temperature values to assess the relative efficiency of each mechanism to localize energy. The local intra-molecular strain energy of the hotspots is also assessed to better understand the physical mechanisms behind the phenomena that lead to a latent potential energy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Machine learning interatomic potential for silicon-nitride (Si 3 N 4 ) by active learning

Silicon nitride (Si 3 N 4 ) is an extensively used material in the automotive, aerospace, and semiconductor industries. However, its widespread use is in contrast to the scarce availability of reliable interatomic potentials that can be employed to study various aspects of this material on an atomistic scale, particularly its amorphous phase. In this work, we developed a machine learning interatomic potential, using an efficient active learning technique, combined with the Gaussian approximation potential (GAP) method. Our strategy is based on using an inexpensive empirical potential to generate an initial dataset of atomic configurations, for which energies and forces were recalculated with density functional theory (DFT); thereafter, a GAP was trained on these data and an iterative re-training algorithm was used to improve it by learning on-the-fly. When compared to DFT, our potential yielded a mean absolute error of 8 meV/atom in energy calculations for a variety of liquid and amorphous structures and a speed-up of molecular dynamics simulations by 3–4 orders of magnitude, while achieving a first-rate agreement with experimental results. Our potential is publicly available in an open-access repository.

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

Charge Transport in Solvated Donor–Acceptor Functionalized Peptoids: Molecular Dynamics and Rate Theory

Scalable solar-energy conversion requires photoactive materials that combine the efficiency of natural photosynthetic systems with the stability and processability needed for practical applications. Achieving reliable charge transport in soft, self-assembled organic materials remains challenging, as structural fluctuations and environmental effects strongly influence charge-transfer (CT) rates. Here, we present a broadly applicable computational framework for evaluating CT rates in the condensed phase, combining Fermi’s golden rule rate theory with inputs from all-atom molecular dynamics (MD) simulations and first-principles electronic-structure calculations. The approach does not rely on system-specific parametrization and is applicable to a wide range of soft and disordered materials. We demonstrate the applicability and usefulness of the framework on redox-active peptoids functionalized with iron–porphyrin (Fe–P) complexes, a bioinspired platform with programmable donor–acceptor units and tunable three-dimensional organization. The calculated CT rates exhibit strong sensitivity to molecular conformation, with variations spanning several orders of magnitude. This dependence is shown to arise from the pronounced variation in diabatic electronic coupling with the relative orientations and separations of the Fe–P complexes across the conformational ensemble. The framework provides a consistent route for connecting atomistic structure to CT kinetics in the condensed phase and enables analysis of structure–rate relations in organic semiconducting systems.

Charge transfer↗
DOE OSTI2022

Ultrafast visualization of incipient plasticity in dynamically compressed matter

Plasticity is ubiquitous and plays a critical role in material deformation and damage; it inherently involves the atomistic length scale and picosecond time scale. A fundamental understanding of the elastic-plastic deformation transition, in particular, incipient plasticity, has been a grand challenge in high-pressure and high-strain-rate environments, impeded largely by experimental limitations on spatial and temporal resolution. Here, we report femtosecond MeV electron diffraction measurements visualizing the three-dimensional (3D) response of single-crystal aluminum to the ultrafast laser-induced compression. We capture lattice transitioning from a purely elastic to a plastically relaxed state within 5 ps, after reaching an elastic limit of ~25 GPa. Our results allow the direct determination of dislocation nucleation and transport that constitute the underlying defect kinetics of incipient plasticity. Large-scale molecular dynamics simulations show good agreement with the experiment and provide an atomic-level description of the dislocation-mediated plasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Localized High-Concentration Electrolytes Get More Localized

Liquid electrolytes in batteries are typically treated as macroscopically homogeneous ionic transport media, despite having complex chemical composition and atomistic solvation structures, thus leaving a knowledge gap of microstructural characteristics. A promising electrolyte system design has recently been developed with the title of a localized high-concentration electrolyte (LHCE). These electrolytes are a summation of three components – a charge-carrying salt, a solvent that solvates the salt, and a diluent that is sparingly interactive with the salt in terms of ion solvation. This creates the perception of a high-concentration electrolyte (HCE) with high anion-cation interactions that supports a salt-derived, stable solid-electrolyte interphase, while an LHCE improves bulk electrolyte properties, such as ionic conductivity and viscosity, when compared to an HCE. Much trial by error design of this LHCE concept has been supported by simulations to confirm ionic interactions, cycling data to confirm cell stability, and post-mortem analysis to provide evidence of mechanisms. However, a systematic design to extend the capabilities of LHCEs is still to be desired. Here, we analyzed a specific LHCE mixture to dictate the best viable mixture that extends the localization of solvation structures while maintaining viable cell-level capabilities. We seek to observe and understand solvent-to-salt interactions that impact electrolyte performance. These interactions range from single-ion solvation, multidentate interactions (i.e., solvation-driven clustering), to more widespread network-level interactions. With the support of both experimental and computational analysis of different electrolyte features, we have begun a methodology of formulating the most ideal LHCE mixtures for improving the cyclability of high-capacity battery chemistries.

25 ENERGY STORAGE↗
DOE OSTI2020

Molecular Origins of the Barriers to Proton Transport in Acidic Aqueous Solutions

The self-consistent iterative multistate empirical valence bond (SCI-MS-EVB) method is used to analyze the structure, thermodynamics, and dynamics of hydrochloric acid solutions. The reorientation timescales of irreversible proton transport are elucidated by simulating 0.43 M, 0.85 M, 1.68 M, and 3.26 M HCl solutions at 270 K, 285 K, 300 K, 315 K, and 330 K. The results indicate increased counter-ion pairing with increasing concentration, which manifests itself via a reduced hydronium oxygen–chloride (O*-Cl) structuring in the radial distribution functions. Increasing ionic concentration also reduces the diffusion of the hydrated excess protons, principally by reducing the contribution of the Grotthuss proton hopping (shuttling) mechanism to the overall diffusion process. In agreement with prior experimental findings, a decrease in the activation energy of reorientation timescales was also observed, which is explicitly explained using activated rate theory and an energy-entropy decomposition of the state-averaged radial distribution functions. Furthermore, these results provide atomistic verification of suggestions from recent two-dimensional infrared spectroscopy experiments that chloride anions (as opposed to hydrated excess protons) create entropic barriers to proton transport.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗