Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuum modeling”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Ionic Interdiffusion at Cathode|Solid-Electrolyte Interface: A Machine Learning–Assisted Multiscale Investigation and Mitigation Strategies

Future lithium batteries are expected to use solid electrolytes to achieve higher energy density and fast charge capabilities. However, most solid electrolytes are thermodynamically unstable against layered oxide cathodes. In this study, the stability of LiCoO2 (LCO) cathode with Li10GeP2S12 (LGPS) solid electrolyte is investigated using ab initio molecular dynamics (AIMD) and machine learning molecular dynamics (MLMD). The propensity of ionic interdiffusion, formation of a passivating interphase layer, and corresponding decay in cell performance is addressed using a continuum model. Large-scale MLMD simulations confirm that the LCO|LGPS interface permits interdiffusion of cobalt (Co) and other ionic species, leading to the formation and growth of a resistive interphase and to dramatic capacity fade even in the first cycle. We examine the literature evidence that incorporating a thin layer of LiNb0.5Ta0.5O3 (LNTO) between LCO and LGPS prevents the interdiffusion of ions. Atomistic simulations suggest that substituting lithium (Li) in LNTO with Co is thermodynamically unfavorable, thereby inhibiting ionic interdiffusion. The stable Nb5+/Ta5+ states form a rigid metal-oxide framework, which consequently also prevents the substitution of niobium (Nb) or tantalum (Ta). However, continuum-level analysis suggests that the higher mechanical stiffness of LNTO can lead to interfacial delamination between the LCO and LNTO. This phenomenon reduces the effectiveness of the protective layer. This paper, therefore, highlights the need to develop novel interlayers that balance low ionic interdiffusion with low mechanical stiffness.

Ncube, Musawenkosi K.↗
DOE OSTI2024

Multiscale modeling-enabled design of multifunctional composites

This study aims to create a comprehensive model that considers multiple scales and physics for predicting the electromechanical behavior of fiber-reinforced composites enhanced with barium titanate (BaTiO3). In our earlier work, we have demonstrated that depositing BaTiO3 microparticles of 200-nm-diameter, on fiber surfaces during fiber-reinforced composite fabrication enhances mechanical strength, passive self-sensing, and energy harvesting properties. The key is to carefully control the microparticle concentration to prevent agglomeration. Since the particles are micron-sized, understanding how agglomeration affects the composites' electromechanical properties is crucial for guiding such multifunctional materials’ design. This study introduces a micromechanics-based approach to explore the impact of microparticle dispersion on the bulk composites' electromechanical properties. Insights gained from this investigation are applied in experiments, enabling accurate predictions of mechanical and self-sensing responses in BaTiO3-enhanced fiber-reinforced composites. Micro-level findings from this computational approach can be integrated into larger continuum models to comprehensively capture the electromechanical behavior of the composite structures at bulk scale. The proposed model is validated by comparing predictions with experimental results, accounting for the nonlinear mechanical and electromechanical behaviors of constituent materials. Consequently, this computational model serves as a digital platform for efficiently designing multifunctional composites.

Gupta, Sumit↗
DOE OSTI2022

Magnus induced diode effect for skyrmions in channels with periodic potentials

Using a particle based model, we investigate the skyrmion dynamical behavior in a channel where the upper wall contains divots of one depth and the lower wall contains divots of a different depth. Under an applied driving force, skyrmions in the channels move with a finite skyrmion Hall angle that deflects them toward the upper wall for -x direction driving and the lower wall for +x direction driving. When the upper divots have zero height, the skyrmions are deflected against the flat upper wall for -x direction driving and the skyrmion velocity depends linearly on the drive. For +x direction driving, the skyrmions are pushed against the lower divots and become trapped, giving reduced velocities and a nonlinear velocity-force response. When there are shallow divots on the upper wall and deep divots on the lower wall, skyrmions get trapped for both driving directions; however, due to the divot depth difference, skyrmions move more easily under -x direction driving, and become strongly trapped for +x direction driving. Further, the preferred -x direction motion produces what we call a Magnus diode effect since it vanishes in the limit of zero Magnus force, unlike the diode effects observed for asymmetric sawtooth potentials. We show that the transport curves can exhibit a series of jumps or dips, negative differential conductivity, and reentrant pinning due to collective trapping events. We also discuss how our results relate to recent continuum modeling on a similar skyrmion diode system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Revisiting stellar properties of star-forming galaxies with stellar and nebular spectral modelling

Spectral synthesis is a powerful tool for interpreting the physical properties of galaxies by decomposing their spectral energy distributions (SEDs) into the main luminosity contributors (e.g. stellar populations of distinct age and metallicity or ionised gas). However, the impact nebular emission has on the inferred properties of star-forming (SF) galaxies has been largely overlooked over the years, with unknown ramifications to the current understanding of galaxy evolution. The objective of this work is to estimate the relations between stellar properties (e.g. total mass, mean age, and mean metallicity) of SF galaxies by simultaneously fitting the stellar and nebular continua and comparing them to the results derived through the more common purely stellar spectral synthesis approach. The main galaxy sample from SDSS DR7 was analysed with two distinct population synthesis codes: FADO, which estimates self-consistently both the stellar and nebular contributions to the SED, and the original version of STARLIGHT, as representative of purely stellar population synthesis codes. Differences between codes regarding average mass, mean age and mean metallicity values can go as high as ~0.06 dex for the overall population of galaxies and ~0.12 dex for SF galaxies (galaxies with EW(Hα) > 3 Å), with the most prominent difference between both codes in the two populations being in the light-weighted mean stellar age. FADO presents a broader range of mean stellar ages and metallicities for SF galaxies than STARLIGHT, with the latter code preferring metallicity solutions around the solar value (Z ⊙ = 0.02). A closer look into the average light- and mass-weighted star formation histories of intensively SF galaxies (EW(Hα) > 75 Å) reveals that the light contributions of simple stellar populations (SSPs) younger than ≤10 7 (10 9 ) years in STARLIGHT are higher by ~5.41% (9.11%) compared to FADO. Moreover, FADO presents higher light contributions from SSPs with metallicity ≤ Z ⊙ /200 (Z ⊙ /50) of around 8.05% (13.51%) when compared with STARLIGHT. This suggests that STARLIGHT is underestimating the average light-weighted age of intensively SF galaxies by up to ~0.17 dex and overestimating the light-weighted metallicity by up to ~0.13 dex compared to FADO (or vice versa). The comparison between the average stellar properties of passive, SF and intensively SF galaxy samples also reveals that differences between codes increase with increasing EW(Hα) and decreasing total stellar mass. Moreover, comparing SF results from FADO in a purely stellar mode with the previous results qualitatively suggests that differences between codes are primarily due to mathematical and statistical differences and secondarily due to the impact of the nebular continuum modelling approach (or lack thereof). However, it is challenging to adequately quantify the relative role of each factor since they are likely interconnected. This work finds indirect evidence that a purely stellar population synthesis approach negatively impacts the inferred stellar properties (e.g. mean age and mean metallicity) of galaxies with relatively high star formation rates (e.g. dwarf spirals, ‘green peas’, and starburst galaxies). In turn, this can bias interpretations of fundamental relations such as the mass-age or mass-metallicity, which are factors worth bearing in mind in light of future high-resolution spectroscopic surveys at higher redshifts (e.g. MOONS and 4MOST-4HS).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Chemo‐Mechanical Coupling in Hydrogels: Dynamics in the Diffusion‐Limited Regime

Hydrogels are characterized by substantial volume changes in response to external stimuli, making them promising candidates for developing smart materials with enhanced adaptability and responsiveness. By integrating chemical reactions, hydrogels acquire dynamic and tunable responsiveness to external stimuli through chemo‐mechanical coupling, expanding their potential in emerging applications. However, capturing their transient behavior remains challenging due to the complex interplay of chemical reactions, solvent transport, and polymer network deformation. Classical theories capture equilibrium swelling but fail to describe time‐dependent phenomena. To address this, a time‐dependent continuum model is developed that explicitly couples these processes. Volume phase transition in hydrogels with homogeneous chemical reactions is investigated, then the effects of reaction kinetics on these transitions are analyzed. The coupling of these mechanisms is further explored through a study of transient mechanical instabilities. To illustrate the impact of distinct reaction and diffusion timescales, a photo‐active gripper is studied for robotic applications. Finally, a photo‐active microswimmer that exhibits non‐reciprocal motion is proposed to highlight the solvent diffusion for locomotion at the micro‐ and nano‐ scales. The work establishes that transient dynamics of chemo‐mechanical hydrogels generate functions not accessible by steady state models and provides a predictive platform for designing adaptive materials in emerging applications.

chemo-mechanical coupling↗
DOE OSTI2020

Molecular Statics Analyses of Thermodynamics and Kinetics of Hydrogen Cottrell Atmosphere Formation Around Edge Dislocations in Aluminum

Aluminum alloys are being explored as lightweight structural materials for use in hydrogen-containing environments. In this paper, to understand hydrogen effects on deformation, we perform molecular statics studies of the hydrogen Cottrell atmosphere around edge dislocations in aluminum. First, we calculate the hydrogen binding energies at all interstitial sites in a periodic aluminum crystal containing an edge dislocation dipole. This allows us to use the Boltzmann equation to quantify the hydrogen Cottrell atmosphere. Based on these binding energies, we then construct a continuum model to study the kinetics of the hydrogen Cottrell atmosphere formation. Finally, we compare our results with existing theories and discuss the effects of hydrogen on deformation of aluminum-based alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The role of deviatoric stress and dislocations on the α to ω phase transformation in Ti

Under extreme conditions, α-Ti becomes unstable and transforms either into β-Ti at high temperature or into ω-Ti at high pressure. In what concerns the α to ω phase transformation (PT), there has been a wide range of experimentally reported transition pressures from approximately 2 to 15 GPa at room temperature. Deviatoric stresses and internal defects are often assumed to be the root cause of this variation. Here, in this study, these postulates are revisited using both continuum mechanics and molecular dynamics (MD) simulations. First, a simple continuum model, assuming linear elasticity and isotropic plasticity, is developed to describe the effects of applied stress and dislocations on the stability of an ω nucleus in an infinite α domain. Second, a new MD simulation method is developed to generate an ω nucleus in the α domain utilizing the displacement field identified from the topological analysis. Results from MD simulations show that despite the fact that phase diagrams typically delineate the limits between two phases in terms of only P and T, deviatoric stress promotes the α to ω phase transformation by reducing the critical radius above which an ω nucleus is stable. Furthermore, the required deviatoric stress to nucleate and stabilize a nanoscale ω nucleus is likely emanating from the internal stress of defects such as dislocations. The MD-informed micromechanics models are used to identify favorable configurations where dislocations help favor the α to ω transformation. These configurations show that the interaction with a basal or prismatic dislocation reduces the critical radius of a ω nucleus by about 10 or 16 %, respectively. In addition, prismatic edge dislocations are found to promote the growth of ω nucleus when interacting with the ($\bar{1}$$\bar{1}$20) α //(0001) ω interface. Importantly, a simple model of the arrival of dislocation at an ω nucleus suggests that PT does not necessarily require a pile-up to be present but could alternatively be mediated by a constant rapid flow of dislocations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Bridging knowledge gaps in liquid- and vapor-fed CO 2 electrolysis through active electrode area

Increased use of gas diffusion electrodes for CO 2 electroreduction widens the experimental phase space that was previously inaccessible using foil electrodes, raising fundamental questions over the impacts of key variables that translate between liquid- and vapor-fed CO 2 electrolysis systems. This work focuses on studying the interplay of current-potential profiles and electrochemically active surface area (ECSA) by implementing a Cu nanoflower catalyst morphology. The results show decreased overpotentials as much as 460 and 174 mV for foil and gas diffusion electrodes, respectively, while maintaining or improving multi-carbon product current density. Furthermore, these overpotential shifts and product activities normalized by ECSA lead to current-potential relationships akin to those of the Tafel description, which are found through a continuum model to be useful for describing the roughness dependence for both liquid- and vapor-fed systems. This analysis establishes a holistic approach for establishing catalyst design criteria to improve materials development for CO 2 .

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Molecular scale insight of pore morphology relation with mechanical properties of amorphous silica using ReaxFF

Porous materials are typically heterogeneous and they contain large variations of micro-/nano-pore structures, causing complicated behaviors. In continuum models, most mechanical properties of porous materials are estimated based on porosity, while the variations of micro/nano structures are ignored. That could be problematic as the microscopic heterogeneity may affect the mechanical response of porous materials. Thus, understanding micro/nano heterogeneity impact has been the focus in many scientific and engineering subjects. In the present study, we investigated the effect of nanopore structure (including pore shape and orientation) as well as porosity on mechanical properties of amorphous silica (a-SiO 2 ). The pore sizes in our simulations are comparable to the corresponding ones observed in a-SiO 2 based materials. We found that the existing of nanopores strongly influences Young’s modulus (E) and critical energy release rate (G IC ). These properties decrease with increasing porosity. Importantly, the impact of nanopores was characterized by structural parameters of porous materials. In addition to dependency on porosity, Young’s modulus also was found to vary as a function of potential energy per atom, which highly depends on nanopore shape. Furthermore, critical energy release rate was found to increase with increasing ligament length (also known as pore wall thickness). Finally, the results highlighted the importance of nanopore structures, which must be taken into account when studying fracture mechanisms in porous materials. Based on our findings, it was proposed that mechanical properties of porous materials can be controlled by nano-engineering pore structures.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Pore-scale simulation of multiphase flow and reactive transport processes involved in geologic carbon sequestration

Multiphase flow and reactive transport are two essential physicochemical processes that govern the effectiveness of geological carbon sequestration (GCS). The interaction and feedback among different phases and components during intricate physicochemical processes hold great significance in understanding CO 2 sequestration. Pore-scale simulations can account for multiphase flow and reactive transport processes in porous media and obtain spatial distributions of parameters (density, velocity, concentration, etc.) in the pore space as well as their temporal evolutions. This proves especially valuable considering that experiments can be hindered by constraints in spatial and temporal resolution. The comprehensive insights garnered from pore-scale research can be leveraged for continuum modeling using the representative elementary volume (REV) concept. In this contribution, four sequential mechanisms of CO 2 -brine-rock interaction in three zones delineated by CO 2 saturation are elaborated to elucidate complicated physicochemical processes involved in GCS, which are followed by general descriptions of mathematical equations and pore-scale numerical methods. In addition, as interested and commonly encountered processes, leakage risks during GCS and CO 2 -enhanced oil recovery (CO 2 -EOR) processes are presented. The existing challenges and future directions are discussed for both the performance of the pore-scale models and the current gaps in the field of GCS. Importantly, we expect that this review will prove beneficial for researchers interested in pore-scale simulations, GCS, and related disciplines.

58 GEOSCIENCES↗
DOE OSTI2021

Quantifying the hierarchy of structural and mechanical length scales in granular systems

Continuum modeling of granular media is made possible by the existence of a length scale at and above which grain-resolved properties can be meaningfully homogenized. Progress has been made in identifying such length scales relevant to local structural properties such as porosity. However, a systematic analysis of scales above which different mechanical properties can be homogenized has yet to emerge. Here, X-ray tomography and 3D X-ray diffraction data are examined to identify such length scales. The data was obtained in-situ in compressed granular materials with rigid and flexible confinement. The experimental data are supplemented with validated discrete element simulations which examine different system sizes and different boundary conditions. Overall, our study reveals a hierarchy in the length scales of granular solids, with lengths governing structural variables being the shortest, lengths of stress variables being intermediate, and lengths of energy dissipation being the longest. All structural and mechanical length scales obey a power law based on the theory of Geostatistics, implying that the length scales can be found by analyzing samples significantly smaller than the length scales themselves. The length scales are also found to be sensitive to boundary conditions, implying that they are extrinsic features of granular media.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Breakdown of electroneutrality in nanopores

Ion transport in extremely narrow nanochannels has gained increasing interest in recent years due to unique physical properties at the nanoscale and the technological advances that allow us to study them. It is tempting to approach this confined regime with the theoretical tools and knowledge developed for membranes and microfluidic devices, and naively apply continuum models, such as the Poisson-Nernst-Planck and Navier-Stokes equations. However, it turns out that some of the most basic principles we take for granted in larger systems, such as the complete screening of surface charge by counter-ions, can break down under extreme confinement. Additionally, we show that in a truly one-dimensional system of ions interacting with three-dimensional electrostatic interactions, the screening length is exponentially large, and can easily exceed the macroscopic length of a nanotube. Without screening, electroneutrality breaks down within the nanotube, with fundamental consequences for ion transport and electrokinetic phenomena. In this work, we build a general theoretical framework for electroneutrality breakdown in nanopores, focusing on the most interesting case of a one-dimensional nanotube, and show how it provides an elegant interpretation for the peculiar scaling observed in experimental measurements of ionic conductance in carbon nanotubes.

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

An adaptive scalable fully implicit algorithm based on stabilized finite element for reduced visco-resistive MHD

The magnetohydrodynamics (MHD) equations are continuum models used in the study of a wide range of plasma physics systems, including the evolution of complex plasma dynamics in tokamak disruptions. However, efficient numerical solution methods for MHD are extremely challenging due to disparate time and length scales, strong hyperbolic phenomena, and nonlinearity. Additionally, therefore the development of scalable, implicit MHD algorithms and high-resolution adaptive mesh refinement strategies is of considerable importance. In this work, we develop a high-order stabilized finite-element algorithm for the reduced visco-resistive MHD equations based on the MFEM finite element library (mfem.org). The scheme is fully implicit, solved with the Jacobian-free Newton-Krylov (JFNK) method with a physics-based preconditioning strategy. Our preconditioning strategy is a generalization of the physics-based preconditioning methods in Chacón et al. (2002) to adaptive, stabilized finite elements. Algebraic multigrid methods are used to invert sub-block operators to achieve scalability. A parallel adaptive mesh refinement scheme with dynamic load-balancing is implemented to efficiently resolve the multi-scale spatial features of the system. Our implementation uses the MFEM framework, which provides arbitrary-order polynomials and flexible adaptive conforming and non-conforming meshes capabilities. Results demonstrate the accuracy, efficiency, and scalability of the implicit scheme in the presence of large scale disparity. The potential of the AMR approach is demonstrated on an island coalescence problem in the high Lundquist-number regime (≥ 10 7 ) with the successful resolution of plasmoid instabilities and thin current sheets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Dislocation dynamics in heterogeneous nanostructured materials

Crystalline materials can be strengthened by introducing dissimilar phases that impede dislocation glide. At the same time, the changes in microstructure and chemistry usually make the materials less ductile. One way to circumvent the strength–ductility dilemma is to take advantage of heterogeneous nanophases which simultaneously serve as dislocation barriers and sources. Owing to their superior mechanical properties, heterogeneous nanostructured materials (HNMs) have attracted a lot of attention worldwide. Nevertheless, it has been difficult to characterize dislocation dynamics in HNMs using classical continuum models, mainly due to the challenges in describing the elastic and plastic heterogeneity among the phases. Here, in this work, we advance a phase-field dislocation dynamics (PFDD) model to treat multi-phase materials, consisting of phases differing in composition, structural order, and size in the same system. We then apply the advanced PFDD model to exploring two important but divergent materials design problems in HNMs: dislocation/obstacle interactions and dislocation/interface interactions. Results show that the interactions between a dislocation and distribution of obstacles varying in structure and composition cannot be understood by simply interpolating from their individual interactions with a dislocation. It is also found that materials containing interfaces with nanoscale thicknesses and compositional gradients have a much higher dislocation bypass stress than those with sharp interfaces, providing an explanation for the simultaneous high strength and toughness of thick interface-containing nanolaminates as observed in recent experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Dislocation nano-hydrides in nickel: Nucleation, evolution and effects on dislocation behaviors

Nano-hydrides have been predicted to precipitate at the core of edge dislocations in the Ni-H system, a mechanism that may promote hydrogen embrittlement. However, nano-hydride nucleation, growth, and effects on dislocation behavior have seldom been explored. This work combines molecular dynamics grand canonical Monte Carlo (MD-GCMC) simulations and continuum modeling to uncover a wide range of phenomena linked to dislocation nano-hydrides. Simulations reveal that nano-hydrides can be stabilized at dislocation cores with all character angles, including screw segments, due to the hydrostatic stresses around the cores of the Shockley partials. Nano-hydride nucleation takes place in these regions, and growth is dictated by the character angle θ of the perfect dislocation. The equilibrium stacking fault width d eq varies dynamically to increase the local hydrostatic stress field and facilitate the formation of the nano-hydride, forming a constriction-like feature and leading to three distinct behaviors: d eq decreases for θ > 30°, d eq remains unchanged for θ = 30°, and d eq increases for θ < 30°. Remote hydrostatic and Escaig stresses are also shown to influence the nucleation stage, implying stress concentrations such as those ahead of crack tips may facilitate nano-hydride precipitation. Moreover, we identify a new hydrogen-induced 60° dislocation reaction that emits a Shockley partial on a conjugate plane, with potential implications for twin nucleation. Testable predictions from this study are then used to reinterpret previous results from the literature. In conclusion, these findings provide a comprehensive framework to assess nano-hydride formation and evolution at dislocations in nickel and other face-centered cubic metals, with important implications to hydrogen embrittlement.

Dislocations↗
DOE OSTI2021

Computational design of microarchitected porous electrodes for redox flow batteries

Porous electrodes are used as the core reactive component across electrochemical technologies. In flowing systems, controlling the fluid distribution, species transport, and reactive environment is critical to attaining high performance. However, conventional electrode materials like felts and papers provide few opportunities for precise engineering of the electrode and its microstructure. To address these limitations, architected electrodes composed of unit cells with spatially varying geometry determined via computational optimization are proposed. Resolved simulation is employed to develop a homogenized description of the constituent unit cells. These effective properties serve as inputs to a continuum model for the electrode when used in the negative half-cell of a vanadium redox flow battery. Porosity distributions minimizing power loss are then determined via computational design optimization to generate architected porosity electrodes. The architected electrodes are compared to bulk, uniform porosity electrodes and found to lead to increased power efficiency across operating flow rates and currents. The design methodology is further used to generate a scaled-up electrode with comparable power efficiency to the bench-scale systems. Finally, the variable porosity architecture and computational design methodology presented here thus offers a novel pathway for automatically generating spatially engineered electrode structures with improved power performance.

25 ENERGY STORAGE↗
DOE OSTI2026

Hydration size-dependent transport of ions across nanoporous graphene membranes

Nanoporous atomically-thin membranes are promising candidates for metal ion separations due to their chemical stability and high permeance and selectivity, but experimental evidence of the mechanisms responsible for ion-ion selectivity is sparse. Here, in this work, we measured the simultaneous diffusion of a dilute mixture of ten different Group I, Group II, and rare earth cations in a salt background across nanoporous graphene (NPG) membranes with sub-nanometer pores. The membranes exhibited ion-ion selectivity, including between similarly-charged ions. Cation transport was governed primarily by the hydrated ion size that was consistent with continuum models of ion diffusion, with additional influence of the ionic charge. Selectivity enhancement was achieved by modifying ion sizes using an ion-selective complexing agent. Our study provides evidence of the importance of size-sieving and electrostatic mechanisms governing ion transport across NPG, and allows for quantitative prediction of transport rates to guide future development of ion-selective atomically-thin membranes.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Hydrated Anions: From Clusters to Bulk Solution with Quasi-Chemical Theory

The interactions of hydrated ions with molecular and macromolecular solution and interface partners are strong on a chemical energy scale. Here we recount the foremost ab initio theory for the evaluation of the hydration free energies of ions, namely, quasi-chemical theory (QCT). We focus on anions, particularly halides but also the hydroxide anion, because they have been outstanding challenges for all theories. For example, this work supports understanding the high selectivity for F – over Cl – in fluoride-selective ion channels despite the identical charge and the size similarity of these ions. QCT is built by the identification of inner-shell clusters, separate treatment of those clusters, and then the integration of those results into the broader-scale solution environment. Recent work has focused on a close comparison with mass-spectrometric measurements of ion-hydration equilibria. We delineate how ab initio molecular dynamics (AIMD) calculations on ion-hydration clusters, elementary statistical thermodynamics, and electronic structure calculations on cluster structures sampled from the AIMD calculations obtain just the free energies extracted from the cluster experiments. That theory–experiment comparison has not been attempted before the work discussed here, but the agreement is excellent with moderate computational effort. This agreement reinforces both theory and experiment and provides a numerically accurate inner-shell contribution to QCT. The inner-shell complexes involving heavier halides display strikingly asymmetric hydration clusters. Asymmetric hydration structures can be problematic for the evaluation of the QCT outer-shell contribution with the polarizable continuum model (PCM). Nevertheless, QCT provides a favorable setting for the exploitation of PCM when the inner-shell material shields the ion from the outer solution environment. For the more asymmetrically hydrated, and thus less effectively shielded, heavier halide ions clustered with waters, the PCM is less satisfactory. We therefore investigate an inverse procedure in which the inner-shell structures are sampled from readily available AIMD calculations on the bulk solutions. This inverse procedure is a remarkable improvement; our final results are in close agreement with a standard tabulation of hydration free energies, and the final composite results are independent of the coordination number on the chemical energy scale of relevance, as they should be. Finally, a comparison of anion hydration structure in clusters and bulk solutions from AIMD simulations emphasize some differences: the asymmetries of bulk solution inner-shell structures are moderated compared with clusters but are still present, and inner hydration shells fill to slightly higher average coordination numbers in bulk solution than in clusters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗