Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “polymorph”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

From Layered Antiferromagnet to 3D Ferromagnet: LiMnBi-to-MnBi Magneto-Structural Transformation

Here, the intermetallic compound LiMnBi was synthesized by the two-step solid-state reaction from the elements. The synthesis temperature of 850 K was selected based on in-situ high-temperature powder X-ray diffraction data. LiMnBi crystalizes in the layered-like PbClF structure type (a = 4.3131(7) Å, c = 7.096(1) Å at 100 K, P4/nmm space group, Z = 2). LiMnBi structure is built of the alternating [MnBi] and Li layers, as determined from single-crystal X-ray diffraction data. Magnetic properties measurements and solid-state 7 Li Nuclear Magnetic Resonance data collected for polycrystalline LiMnBi samples indicate the long-range antiferromagnetic ordering of Mn sublattice at ~340 K, with no superconductivity down to 5 K detected. LiMnBi is air- and water-sensitive. In aerobic conditions, Li can be extracted from LiMnBi structure to form Li 2 O/LiOH and MnBi (NiAs structure type, P6 3 /mmc). The obtained MnBi polymorph was previously reported to be one of the strongest rare-earth-free ferromagnets, yet its bulk synthesis in powder form is cumbersome. The proposed magneto-structural transformation from ternary LiMnBi to ferromagnetic MnBi involves condensation of the MnBi4 tetrahedra upon Li deintercalation and is exclusive to LiMnBi. In contrast, ferromagnetic MnBi cannot be obtained from either isostructural NaMnBi and KMnBi, or from the structurally related CaMn 2 Bi 2 . Such a distinctive transformation in the case of LiMnBi is presumed to be due to its fitting reactivity to yield MnBi and favorable interlayer distance between [MnBi] layers, while the interlayer distance in NaMnBi and KMnBi structural analogs is unfavorably long. The studies of delithiation from the layered-like LiMnBi under different chemical environments indicate that the yield of the MnBi depends on the type of solvent used and the kinetics of the reaction. A slow rate and mild reaction media lead to a high fraction of the MnBi product. The saturation magnetization of the “as-prepared” MnBi is ~50 % of the expected value of 81.3 emu/g. Overall, this study adds a missing member to the family of ternary pnictides and illustrates how soft-chemistry methods can be used to obtain “difficult-to-synthesize” compounds.

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

Toggling Stereochemical Activity through Interstitial Positioning of Cations between 2D V 2 O 5 Double Layers

The 5/6s 2 lone-pair electrons of p-block cations in their lower oxidation states are a versatile electronic and geometric structure motif that can underpin lattice anharmonicity and often engender electronic and structural instabilities that underpin the function of active elements in nonlinear optics, thermochromics, thermoelectrics, neuromorphic computing, and photocatalysis. In contrast to periodic solids where lone-pair-bearing cations are part of the structural framework, installing lone-pair-bearing cations in the interstitial sites of intercalation hosts provides a means of a systematically modulating electronic structure through the choice of the group and the period of the inserted cation while preserving the overall framework connectivity. The extent of stereochemical activity and the energy positioning of lone-pair-derived mid-gap states depend on the cation identity, stoichiometry, and strength of anion hybridization. V 2 O 5 polymorphs are versatile insertion hosts that can accommodate a broad range of s-, p-, and d-block cations. However, the insertion of lone-pair-bearing cations remains largely underexplored. In this article, we examine the implications of varying the 6s 2 cations situated in interlayer sites between condensed [V 4 O 10 ] n double layers. Systematic modulations of lattice distortions, electronic structure, and magnetic ordering are observed with increasing strength of stereochemical activity from group 12 to group 14 cations. We compare and contrast p-block-layered M x V 2 O 5 (M = Hg, Tl, and Pb) compounds and map the significance of local off-centering arising from the stereochemical activity of lone-pair cations to the emergence of filled antibonding lone-pair 6s 2 –O 2p-hybridized mid-gap states mediated by second-order Jahn–Teller distortions. Crystallographic studies of cation coordination environments and the resulting modulation of V–V interactions have been used in conjunction with variable-energy hard X-ray photoelectron spectroscopy measurements, first-principles electronic structure calculations, and crystal orbital Hamilton population analyses to decipher the origins of stereochemical activity. Magnetic susceptibility measurements reveal antiferromagnetic signatures for all the three compounds. However, the differences in V–V interactions significantly affect the energy balance of the superexchange interactions, resulting in an ordering temperature of 160 and 260 K for Hg 0.5 V 2 O 5 and δ-Tl 0.5 V 2 O 5 , respectively, as compared to 7 K for δ-Pb 0.5 V 2 O 5 . In δ-Pb 0.5 V 2 O 5 , the strong stereochemical activity of electron lone pairs and the resulting electrostatic repulsions enforce superlattice ordering, which strongly modifies the electronic localization patterns along the [V 4 O 10 ] slabs, resulting in disrupted magnetic ordering and an anomalously low ordering temperature. The results demonstrate a versatile strategy for toggling the stereochemical activity of electron lone pairs to modify the electronic structure near the Fermi level and to mediate superexchange interactions.

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

Tuning Optical and Electrical Properties of Vanadium Oxide with Topochemical Reduction and Substitutional Tin

Vanadium oxides are widely tunable materials, with many thermodynamically stable phases suitable for applications spanning catalysis to neuromorphic computing. The stability of vanadium in a range of oxidation states enables mixed-valence polymorphs of kinetically accessible metastable materials. Low-temperature synthetic routes to, and the properties of, these metastable materials are poorly understood and may unlock new optoelectronic and magnetic functionalities for expanded applications. In this work, we demonstrate topochemical reduction of α-V 2 O 5 to produce metastable vanadium oxide phases with tunable oxygen vacancies (>6%) and simultaneous substitutional tin incorporation (>3.5%). The chemistry is carried out at low temperature (65 °C) with solution-phase SnCl 2 , where Sn 2+ is oxidized to Sn 4+ as V 5+ sites are reduced to V 4+ during oxygen vacancy formation. Despite high oxygen vacancy and tin concentrations, the transformations are topochemical in that the symmetry of the parent crystal remains intact, although the unit cell expands. Band structure calculations show that these vacancies contribute electrons to the lattice, whereas substitutional tin contributes holes, yielding a compensation doping effect and control over the electronic properties. The SnCl 2 redox chemistry is effective on both solution-processed V 2 O 5 nanoparticle inks and mesoporous films cast from untreated inks, enabling versatile routes toward functional films with tunable optical and electronic properties. The electrical conductance rises concomitantly with the SnCl 2 concentration and treatment time, indicating a net increase in density of free electrons in the host lattice. This work provides a valuable demonstration of kinetic tailoring of electronic properties of vanadium–oxygen systems through top-down chemical manipulation from known thermodynamic phases.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Reversible Disorder-to-Order Transition Induced by Aqueous Lithiation in Vanadate Electrode Materials

Vanadium-based oxides are intriguing electrode materials in aqueous electrochemical systems owing to their low cost and high theoretical capacity for alkali storage, especially lithium (Li) ions. However, a sequence of phase transformations and irreversible structure distortion upon Li-ion intercalation causes structural instability and has been a lingering problem for vanadium oxide electrodes. Here, in this work, we investigate lithium vanadate (Li–V 3 O 8 ) for aqueous Li-ion intercalation and deintercalation processes. Unlike its crystalline V 2 O 5 polymorph, Li–V 3 O 8 retains monophasic lithiation, which is attributed to its disordered crystalline nature and large interplanar distance. Importantly, we show a unique and reversible sequence of disorder-to-order structural transition induced by the extent of lithiation, which indicates sequential interlayer and intralayer lithiation process, and vice versa in delithiation process, supported by electrokinetic analysis, in situ X-ray diffraction (XRD), and Debye scattering simulations. The absence of distortive phase transitions and multilithiation pathways facilitates Li-ion diffusion across the vanadate electrode materials to improve storage capacity. This work opens a new dimension for vanadium-based disordered oxides, accelerating the development of low-cost, aqueous electrochemical systems.

36 MATERIALS SCIENCE↗
DOE OSTI2025

High-Pressure Phase Transition of Metastable Wurtzite-Like CuInSe 2 Nanocrystals

Ternary I–III–VI 2 semiconductors, such as CuInSe 2 , exhibit diverse polymorphs with unique structural characteristics and optoelectronic properties. This study investigates the pressure-induced phase transitions of metastable wurtzite-like CuInSe 2 nanocrystals. Using a combination of synchrotron X-ray diffraction, pair distribution function analysis, and density functional theory calculations, we reveal a transition from cation-ordered wurtzite-like (Pmc2 1 ) to cation-disordered NaCl-like (Fm3̅m) structures at 7.7 GPa. The cation-disordered NaCl-like phase persists upon decompression. Bulk modulus calculations highlight size-dependent deviations from bulk material behavior. These findings deepen our understanding of phase stability in colloidal I–III–VI 2 semiconductor nanocrystals, with implications for tailoring functional materials under extreme conditions.

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

Single Crystals of Vanadium Oxides as a Lens for Understanding Structural and Electronic Phase Transformations, Ion Transport, Chemo-Mechanical Coupling, and Electrothermal Neuronal Emulation

Vanadium oxides cystallize in a diverse array of structures and compositions arising from the redox versatility of vanadium, variable covalency of V−O bonds, and myriad coordination geometries. Their open frameworks present abundant interstitial sites that enable insertion of guest-ions. In such compounds, V3d electron and spin localization and disorder couple strongly to structural preferences. The rich structural diversity manifests as a “rugged” free energy landscape with multiple interconvertible polymorphs. Such a landscape sets up structural, electronic, and magnetic transitions that underpin the promise of these materials as ion-insertion battery electrodes; compact primitives for brain-inspired computing, and heterogeneous catalysts. Here, we examine the structural and compositional diversity, electronic instabilities, defect dynamics, structure transformations, mechanical properties, and surface structure of vanadium oxides using single crystals as a distinctive lens. Single crystals enable the measurement of structure−function correlations without the ensemble and orientational averaging inevitable in polycrystalline materials. Their well-defined surfaces further enable examination of facet-dependent reactivity toward molecular adsorbates, ion fluxes, and lattice (mis)matched solids. We provide a comprehensive account of vanadium-oxide single-crystal studies, from delineation of common structural motifs to single-crystal growth techniques, topochemical modification strategies, mechanisms underpinning electronic instabilities, and implementation as electrothermal neurons and battery electrode materials.

Ponis, John [Texas A&M University, College Station↗
DOE OSTI2021

Formation Pathway of Wurtzite-like Cu 2 ZnSnSe 4 Nanocrystals

Cu 2 ZnSnSe 4 is a direct band gap semiconductor composed of Earth-abundant elements, making it an attractive material for thin film photovoltaic technologies. Cu 2 ZnSnSe 4 crystallizes in the kesterite structure type as a bulk material, but it can also crystallize in a metastable wurtzite-like crystal structure when synthesized on the nanoscale. The wurtzite-like polymorph introduces unique and useful properties to Cu 2 ZnSnSe 4 materials, including widely tunable band gaps and superior composi-tional flexibility as compared to kesterite Cu 2 ZnSnSe 4 . Here, we investigate the formation pathway of colloidally prepared wurtzite-like Cu 2 ZnSnSe 4 nanocrystals. We show that this quaternary material forms through a chain of reactions, starting with binary Cu 3 Se 2 nanocrystals that, due to both kinetic and thermodynamic reasons, preferentially react with tin to yield hexagonal copper tin selenide intermediates. Furthermore, these ternary intermediates then react with zinc to form the resulting wurtzite-like Cu 2 ZnSnSe 4 nanocrystals. Based on this formation pathway, we suggest synthetic methods that may prevent the for-mation of unwanted impurity phases that are known to hamper the efficiency of Cu 2 ZnSnSe 4 -based optoelectronic devices.

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

A Low-Temperature Structural Transition in Canfieldite, Ag 8 SnS 6 , Single Crystals

Canfieldite, Ag 8 SnS 6 , is a semiconducting mineral notable for its high ionic conductivity, photosensitivity, and low thermal conductivity. In this paper, we report the solution growth of large single crystals of Ag 8 SnS 6 of mass up to 1 g from a ternary Ag–Sn–S melt. On cooling from high temperature, Ag 8 SnS 6 undergoes a known cubic ($F\bar{4}3m$) to orthorhombic ($Pna2_1$) phase transition at ≈460 K. By studying the magnetization and thermal expansion between 5–300 K, we discover a second structural transition at ≈120 K. Single crystal X-ray diffraction reveals the low-temperature phase adopts a different orthorhombic structure with space group $Pmn2_1$ ($\textit{a}$ = 7.662 9(5) Å, $\textit{b}$ = 7.539 6(5) Å, $\textit{c}$ = 10.630 0(5) Å, Z = 2 at 90 K) that is isostructural to the room-temperature forms of the related Se-based compounds Ag 8 SnSe 6 and Ag 8 GeSe 6 . The 120 K transition is first-order and has a large thermal hysteresis. On the basis of the magnetization and thermal expansion data, the room-temperature polymorph can be kinetically arrested into a metastable state by rapidly cooling to temperatures below 40 K. We last compare the room- and low-temperature forms of Ag8SnS6 with its argyrodite analogues, Ag 8 TQ 6 ($\textit{T}$ = Si, Ge, Sn; $\textit{Q}$ = S, Se), and identify a trend relating the preferred structures to the unit cell volume, suggesting smaller phase volume favors the $Pna2_1$ arrangement. We support this picture by showing that the transition to the $Pmn2_1$ phase is avoided in Ge alloyed Ag 8 Sn 1–x Ge x S 6 samples as well as in pure Ag 8 GeS 6 .

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

Anomalous Ferromagnetic Behavior in Orthorhombic Li 3 Co 2 SbO 6

Monoclinic Li 3 Co 2 SbO 6 has been proposed as a Kitaev spin liquid candidate and investigated intensively, whereas the properties of its polymorph, the orthorhombic phase, are less known. In this work, we report the magnetic properties of orthorhombic Li 3 Co 2 SbO 6 as revealed by dc and ac magnetic susceptibility, muon spin relaxation (μSR), and neutron diffraction measurements. Successive magnetic transitions at 115, 89, and 71 K were observed in the low-field dc susceptibility measurements. The transitions below T N (115 K) are suppressed at higher applied fields. However, zero-field ac susceptibility measurements reveal distinct frequency-independent transitions at about 114, 107, 97, 79, and 71 K. A long-range magnetic ordered state was confirmed by specific heat, μSR, and neutron diffraction measurements, all indicating a single transition at about 115 K. The discrepancy between different measurements is attributed to possible stacking faults and/or local disorders of the ferromagnetic zigzag chains, resulting in ferromagnetic boundaries within the overall antiferromagnetic matrix.

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

Two Neptunium(III) Mellitate Coordination Polymers: Completing the Series Np–Cf of Trans-Uranic An(III) Mellitates

For this work, two neptunium(III) mellitates, 237 Np 2 (mell)(H 2 O) 9 ·1.5H 2 O (Np-1α) and 237 Np 2 (mell)(H 2 O) 8 ·2H 2 O (Np-1β), have been synthesized from 237 NpCl 4 (dme) 2 by reduction with KC 8 and subsequent reaction with an aqueous solution of mellitic acid (H 6 mell). Characterization by single-crystal X-ray crystallography and UV–vis–NIR spectroscopy confirms that the neptunium is in its +3 oxidation state and both polymorphs are isostructural to the previously reported plutonium mellitates. Of the two morphologies, Np-1α is indefinitely stable in air, while Np-1β slowly oxidizes over several months. This is due to the change in the energy of the metal-ligand charge-transfer absorption exhibited by these compounds attributed to differing numbers of carboxylate bonds to Np(III), where in Np-1β the energy is low enough to result in spontaneous oxidation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Phase Discovery and Selected Synthesis of Subvalent Niobium Tellurides Using a Polytelluride Flux Strategy

Transition metal subchalcogenides involve electron-rich metals and can facilitate an in-depth understanding of the relationships among quantum properties such as superconductivity, charge density wave, and topological band structures. However, effective experimental routes toward synthesizing transition metal subchalcogenides are still lacking, hindering the development of new quantum materials. Herein, we propose a eutectic polytelluride flux strategy as an excellent solution to address phase discovery and crystal growth in transition metal subtelluride systems. We report new phases easily and selectively synthesized using a eutectic “K 3 Te 4 ” polytelluride flux upon adjusting the ratio of Nb metal to flux in the starting materials (K/Nb/Te = 3:x:4). Using a high Nb content in the solvent (x = 2 and 1), crystals of KNb 3 Te 3 O 0.38 and K 0.9 Nb 3 Te 4 are obtained. Both subtellurides exhibit diverse Nb clusters, including face-sharing and edge-sharing Nb 6 octahedral columns and zig-zag Nb chains. Reducing the Nb content to x = 0.33 leads to the formation of a layered compound, K 1.06 NbTe 2 . This compound comprises a NbTe 6 trigonal prism with K intercalated between the layers. Single crystals of known binary Nb tellurides can also be grown using another eutectic flux “KTe 3 . 2 ”, and the obtained NbTe 2 exhibits a new polymorphism with extra trimerization along the b-axis in the Nb–Nb bonded double zig-zag cluster. Finally, precise control over the structural dimensionality and oxidation state, combined with the facile crystal growth process, makes our synthetic strategy an efficient route to explore quantum materials in transition metal subchalcogenides.

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

In situ characterization of metastable Pb3O5 and Pb2O3 phases during thermal decomposition of PbO2 to PbO

Nonstoichiometric lead oxides play a key role in the formation and cycling of the positive electrodes in a lead acid battery. These phases have been linked to the underutilization of the positive active material but also play a key role in the battery’s cycle life, providing inter-particle adhesion and the connection to the underlying lead grid. Similar phases have previously been identified by mass loss or color change during thermal annealing of PbO2 to PbO, suggesting that at least two intermediate PbOx phases exist. Using multiple in situ analysis techniques (powder diffraction, x-ray absorption, x-ray photoelectron spectroscopy) and ex situ nuclear magnetic resonance measurements, the structural conversion and changes in the lead oxidation states were identified during this process. Isolation of the PbOx phases enabled confirmation of Pb3O5 and Pb2O3 by diffraction and the first 207Pb NMR measurement of these intermediates.

Kinnibrugh, Tiffany L.↗
DOE OSTI2026

Pressure-Induced Reduction of Dicyanamide by Samarium(II) in a Coordination Polymer

A samarium(II) coordination polymer, [Sm(2.2.2- cryptand)(dca)]I (dca− = dicyanamide), has been prepared from the reaction of SmI2 with tetrabutylammonium dicyanamide and 2.2.2- cryptand. The structure consists of [Sm(2.2.2-cryptand)] 2+ cations bridged by dicyanamide to form corrugated 1D chains. Single crystals of this compound have been studied spectroscopically using UV−vis−NIR and Raman spectroscopy as a function of pressure to reveal a pressureinduced two-electron reduction of dicyanamide by Sm 2+ to cyanide and cyanamide and Sm 3+ between 5 and 8 GPa. High-pressure single-crystal diffraction studies reveal a phase change at 2.6 ± 0.1 GPa, where a polymorph exhibiting a smaller Sm 2+ ···N C bond angle and shorter Sm 2+ ···C distance was observed. These observations provide mechanistic insight into the reduction of dicyanamide into cyanide and cyanamide by pressurized samarium(II).

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

New Nonreactive Force Field for Accurate Molecular Dynamics Simulations of TATB at Extreme Conditions

Insensitive high explosives based on TATB (1,3,5-triamino-2,4,6-trinitrobenzene) are needed in applications when safety is of paramount importance, but the basic material properties that give rise to its insensitivity are not fully understood. Molecular dynamics modeling using empirical force fields (FFs) has been the main route to characterize many complicated dynamical properties of TATB single crystal, but these FFs have not been comprehensively tested at extreme conditions typical of detonation. We collect a benchmark data set of (quasi)static TATB physical properties as determined by experiments and electronic structure calculations and apply this data set to validate four existing TATB FFs along with a new TATB FF that we develop here and denote as the CEA-LLNL-Missouri (CLM) FF. Benchmark data include vibrational spectra, the TATB crystal temperature–pressure–volume equation of state and lattice parameters, properties of TATB crystal polymorphs and transitions to the gaseous and liquid states, dimer energy landscapes, the pressure-dependent elastic tensor, and the energy landscape for inelastic deformation via sliding of TATB crystal layers. As a general assessment, we find that the two existing nonreactive FFs are more accurate in describing TATB’s physical properties compared to the two variants of the ReaxFF reactive FF considered. The new CLM FF is found to consistently yield similar or better agreement with experiments and electronic structure theory than any of the existing FF models, and it presents a distinct improvement in accurately modeling TATB elasticity and equation of state. So this work is expected to help improve the accuracy of FF-based modeling of complicated dynamic responses that ultimately govern the safety and performance characteristics of this material.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Ice Phase Classification Made Easy with Score-Based Denoising

Accurate identification of ice phases is essential for understanding various physicochemical phenomena. However, such classification for structures simulated with molecular dynamics is complicated by the complex symmetries of ice polymorphs and thermal fluctuations. For this purpose, both traditional order parameters and data-driven machine learning approaches have been employed, but they often rely on expert intuition, specific geometric information, or large training data sets. In this work, we present an unsupervised phase classification framework that combines a score-based denoiser model with a subsequent model-free classification method to accurately identify ice phases. Further, the denoiser model is trained on perturbed synthetic data of ideal reference structures, eliminating the need for large data sets and labeling efforts. The classification step utilizes the smooth overlap of atomic position (SOAP) descriptors as the atomic fingerprint, ensuring Euclidean symmetries and transferability to various structural systems. Our approach achieves a remarkable 100% accuracy in distinguishing ice phases of test trajectories using only seven ideal reference structures of ice phases as model inputs. This demonstrates the generalizability of the score-based denoiser model in facilitating phase identification for complex molecular systems. The proposed classification strategy can be broadly applied to investigate structural evolution and phase identification for a wide range of materials, offering new insights into the fundamental understanding of water and other complex systems.

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

Phase Equilibrium of Water with Hexagonal and Cubic Ice Using the SCAN Functional

Machine learning models are rapidly becoming widely used to simulate complex physicochemical phenomena with ab initio accuracy. Here, we use one such model as well as direct density functional theory (DFT) calculations to investigate the phase equilibrium of water, hexagonal ice (Ih), and cubic ice (Ic), with an eye toward studying ice nucleation. The machine learning model is based on deep neural networks and has been trained on DFT data obtained using the SCAN exchange and correlation functional. We use this model to drive enhanced sampling simulations aimed at calculating a number of complex properties that are out of reach of DFT-driven simulations and then employ an appropriate reweighting procedure to compute the corresponding properties for the SCAN functional. This approach allows us to calculate the melting temperature of both ice polymorphs, the driving force for nucleation, the heat of fusion, the densities at the melting temperature, the relative stability of ices Ih and Ic, and other properties. We find a correct qualitative prediction of all properties of interest. In some cases, quantitative agreement with experiment is better than for state-of-the-art semiempirical potentials for water. Furthermore, our results also show that SCAN correctly predicts that ice Ih is more stable than ice Ic.

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

Polarizable Water Potential Derived from a Model Electron Density

A new empirical potential for efficient, large scale molecular dynamics simulation of water is presented. The HIPPO (Hydrogen-like Intermolecular Polarizable POtential) force field is based upon the model electron density of a hydrogen-like atom. This framework is used to derive and parametrize individual terms describing charge penetration damped permanent electrostatics, damped polarization, charge transfer, anisotropic Pauli repulsion, and damped dispersion interactions. Initial parameter values were fit to Symmetry Adapted Perturbation Theory (SAPT) energy components for ten water dimer configurations, as well as the radial and angular dependence of the canonical dimer. The SAPT-based parameters were then systematically refined to extend the treatment to water bulk phases. The final HIPPO water model provides a balanced representation of a wide variety of properties of gas phase clusters, liquid water, and ice polymorphs, across a range of temperatures and pressures. This water potential yields a rationalization of water structure, dynamics, and thermodynamics explicitly correlated with an ab initio energy decomposition, while providing a level of accuracy comparable or superior to previous polarizable atomic multipole force fields. The HIPPO water model serves as a cornerstone around which similarly detailed physics-based models can be developed for additional molecular species.

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

Hierarchical Reinforcement Learning of a Short-Range Bond-Order Potential for Silica: Analytic Embedding of Coordination with Classical Efficiency

Reinforcement learning (RL) has recently emerged as a data-efficient strategy to parametrize short-range interatomic potentials. Building on our past RL optimization of pairwise silica models, we extend the framework to a bond-order (Tersoff-type) potential that provides an analytic embedding of local coordination through a three-body term. A hierarchical RL workflow combining continuous-action Monte Carlo Tree Search and property-based rewards efficiently explores the 26-dimensional parameter space, sequentially optimizing lattice parameters, densities, angles, and cohesive energies of 21 silica polymorphs. The resulting models, Q-Tersoff and ML-Tersoff, reproduce the energetic ordering of low-energy phases and capture the angular correlations and amorphous structure factors of silica with improved fidelity over pairwise force fields, while remaining orders of magnitude faster than high-dimensional machine-learned potentials. Both models underperform for elastic constants and high-energy frameworks, delineating the limits of the current analytic form. The approach establishes a general and interpretable route to angle-aware, short-range potentials that bridge physics-based and machine-learned descriptions of silicate materials.

36 MATERIALS SCIENCE↗