Search NASA⌕ Search

SEARCH · Search NASA

Results for “coordination complexes”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Unraveling Local Structure of Molten Salts via X-ray Scattering, Raman Spectroscopy, and Ab Initio Molecular Dynamics

In this work, we resolve a long-standing issue concerning the local structure of molten MgCl 2 by employing a multimodal approach, including X-ray scattering and Raman spectroscopy, along with the theoretical modeling of the experimental spectra based on ab initio molecular dynamics (AIMD) simulations utilizing several density functional theory (DFT) methods. We demonstrate the reliability of AIMD simulations in achieving excellent agreement between the experimental and simulated spectra for MgCl 2 and 50 mol % MgCl 2 + 50 mol % KCl, and ZnCl 2 , thus allowing structural insights not directly available from experiment alone. A thorough computational analysis using five DFT methods provides a convergent view that octahedrally coordinated magnesium in pure MgCl 2 upon melting preferentially coordinates with five chloride anions to form distorted square pyramidal polyhedra that are connected via corners and to a lesser degree via edges. This is contrasted with the results for ZnCl 2 , which does not change its tetrahedral coordination on melting. Although the five-coordinate MgCl 5 3– complex was not considered in the early literature, together with an increasing tendency to form a tetrahedrally coordinated complex with decreasing the MgCl 2 content in the mixture with alkali metal chloride systems, current work reconciles the results of most previous seemingly contradictory experimental studies.

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

Second Coordination Sphere Effects in an Evolved Ru Complex Based on Highly Adaptable Ligand Results in Rapid Water Oxidation Catalysis

A new Ru complex containing the deprotonated 2,2':6',2"-terpyridine-6,6"-diphosphonic acid (H 4 tPa) and pyridine (py) of general formula [Ru II (H 3 tPa-κ-N 3 O)(py) 2 ] + , 2 + , has been prepared and thoroughly characterized by means of spectroscopic and electrochemical techniques, X-ray diffraction analysis, and density functional theory (DFT) calculations. Complex 2 + presents a dynamic behavior in the solution that involves the synchronous coordination and the decoordination of the dangling phosphonic groups of the tPa 4– ligand. Yet, at oxidation state IV, complex 2 + becomes seven coordinated with the two phosphonic groups now bonded to the metal center. Further, at this oxidation state at neutral and basic pH, the Ru complex undergoes the coordination of an exogenous OH – group from the solvent that leads to an intramolecular aromatic O atom insertion into the CH bond of one of the pyridyl groups, forming the corresponding phenoxo-phosphonate Ru complex [Ru III (tPaO-κ-N 2 O P O C )(py) 2 ] 2– , 4 2– , where tPaO 5– is the 3-(hydroxo-[2,2':6',2"-terpyridine]-6,6"-diyl)bis(phosphonate) ligand. This new in situ generated Ru complex, 4 2– , has been isolated and spectroscopically and electrochemically characterized. In addition, a crystal structure has been also obtained using single-crystal X-ray diffraction techniques. Complex 4 2– turns out to be an exceptional water oxidation catalyst achieving record maximum turnover frequencies (TOF max ) on the order of 16 000 s –1 . A mechanistic analysis complemented with DFT calculations has further been carried out, showing the critical role of intramolecular second coordination sphere effects exerted by the phosphonate groups in lowering the activation energy at the rate-determining step.

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

Iron(II) Complexes of an Anionic Bis(ylide)diphenylborate Ligand

Double deprotonation of the salt [Ph 2 B(PMe 3 ) 2 ][OTf] provides access to a bis(ylide)diphenylborate ligand that is readily transferred in situ to iron(II). Depending on the reaction stoichiometry, both the “ate” complex [Ph 2 B(Me 2 PCH 2 ) 2 Fe(μ-Cl) 2 Li(THF) 2 ] and the homoleptic complex [Ph 2 B(Me 2 PCH 2 ) 2 ] 2 Fe can be prepared from FeCl 2 (THF) 1.5 . Further reaction of with FeCl 2 (THF) 1.5 produces the chloride-bridged dimer [Ph 2 B(Me 2 PCH 2 ) 2 Fe(μ-Cl) 2 Fe(CH 2 PMe 2 ) 2 BPh 2 ]. Attempts to reduce or alkylate provide as the only isolable product, likely a consequence of the low steric hindrance of the bis(ylide)diphenylborate ligand. On the other hand, reaction of with the strong field ligand CN t Bu provides the six-coordinate, diamagnetic complex [Ph 2 B(Me 2 PCH 2 ) 2 Fe(CN t Bu) 4 ][Cl]. Here, the electronic structure calculations for the bis(ylide)diphenylborate ligand and homoleptic complex suggest that the C(ylide) atoms are strong σ-donors with little π-bonding character. These initial results suggest the potential for this bis(ylide)diphenylborate ligand in coordination chemistry.

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

Synthesis and spectroscopic characterization of actinyl(VI) tetrahalide coordination compounds containing 2, 2'-bipyridine

A series of three new actinide tetrahalide coordination compounds containing 2,2'-bipyridinium (H 2 bpy) have been synthesized by room temperature evaporation and characterized via single crystal X-ray diffraction, Raman spectroscopy, and infrared spectroscopy. The solid-state compounds (Np1, U1, and U2) contain the AnO 2 2+ cation coordinated to four halide (chloride or bromide) atoms to form the anionic tetrahalide molecular unit [AnO 2 X 4 ] 2– (X = Cl, Br) and are charge balanced by the 2,2'-bipyridinium cations that engage in electrostatic interactions with the actinyl complex. Np1 (N 2 C 10 H 10 )[NpO 2 Cl 4 ] contains 2,2'-bipyridine in the cis conformation where both the nitrogen atoms are directed towards the [NpO 2 Cl 4 ] 2– monomer. Within U1 (N 2 C 10 H 9 )[UO 2 Cl 4 ], the 2,2'-bipyridine is in a trans conformation, and found in the same orientation in the isomorphous structure U2 (N 2 C 10 H 9 )[UO 2 Br 4 ]. The compounds were further characterized by vibrational spectroscopy with specific interest in the actinyl cation. Raman spectroscopy provided information on the actinyl symmetric stretch (ν 1 ) whereas IR spectroscopy is used to determine the location of the asymmetric stretch (ν 3 ). Here, spectral similarities are observed between the two isomorphous structures, but slight variations across the series of structures are found to be associated with bipyridine bands in the various conformations. Additional analysis of these vibrational features in U1 and U2 provided information on the actinyl force constants, allowing insight into the relative uranyl bond strength that was compared to similar tetrahalide compounds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

Preliminary Application of Formal Verification to An Autonomy Architecture for Unmanned Aircraft

There is a desire to design autonomous systems in such a way that capabilities can be easily added or re- combined to produce new behaviors while preserving their safety properties. ICAROUS, a prototype software architecture for building safety-centric autonomous unmanned aircraft applications, is designed to support this type of extensibility and re-configurability. In ICAROUS, core capabilities are implemented as individual soft- ware services, so that enabling access to new capabilities simply requires adding new services. To make use of these capabilities, ICAROUS includes a specialized service that provides a general framework for config- uring the relative priorities, conditions, and rules that govern how different modules should be engaged and disengaged during flight. The inherent complexity of coordinating multiple modules under changing conditions makes it difficult to determine whether a particular configuration could have erroneous behaviors in certain circumstances. A robust set of integration tests can help discover errors, but testing can only realistically cover a relatively small proportion of total system behaviors. Developing good tests and interpreting the results to pinpoint the cause of errors when they arise can also be very time-consuming. To supplement testing, formal methods can be used to model and analyze complex systems, achieving better coverage and simplifying the process of finding, understanding, and fixing errors. To demonstrate these benefits, this paper explores the ap- plication of formal methods to ICAROUS. In particular, the Spin model checker is used to specify requirements for and model portions of the system, then verify whether the model satisfies the requirements and find and fix errors when it does not.

Formal Methods↗
NASA NTRSDate not supplied

Preliminary Application of Formal Verification to An Autonomy Architecture for Unmanned Aircraft

There is a desire to design autonomous systems in such a way that capabilities can be easily added or re-combined to produce new behaviors while preserving their safety properties. ICAROUS, a prototype software architecture for building safety-centric autonomous unmanned aircraft applications, is designed to support this type of extensibility and re-configurability. In ICAROUS, core capabilities are implemented as individual soft- ware services, so that enabling access to new capabilities simply requires adding new services. To make use of these capabilities, ICAROUS includes a specialized service that provides a general framework for config- uring the relative priorities, conditions, and rules that govern how different modules should be engaged and disengaged during flight. The inherent complexity of coordinating multiple modules under changing conditions makes it difficult to determine whether a particular configuration could have erroneous behaviors in certain circumstances. A robust set of integration tests can help discover errors, but testing can only realistically cover a relatively small proportion of total system behaviors. Developing good tests and interpreting the results to pinpoint the cause of errors when they arise can also be very time-consuming. To supplement testing, formal methods can be used to model and analyze complex systems, achieving better coverage and simplifying the process of finding, understanding, and fixing errors. To demonstrate these benefits, this paper explores the ap- plication of formal methods to ICAROUS. In particular, the Spin model checker is used to specify requirements for and model portions of the system, then verify whether the model satisfies the requirements and find and fix errors when it does not.

Formal Methods↗
DOE OSTI2024

Optimal operation of multi-plant steam district heating systems for enhanced efficiency and sustainability

Despite their crucial role in supplying heat and power to universities, industries, and healthcare facilities, many steam-based district heating systems rely on outdated control methods. Among these, multi-central plant districts are particularly challenging due to the complexities of coordinating multiple plants, optimizing load distributions, and managing system downtime. In response, new operational strategies are developed to enhance the efficiency and sustainability of steam districts while utilizing existing resources. These strategies include reducing plant operational pressure without compromising the reliable supply to buildings and optimizing load allocation across multiple plants. The load allocation considers boiler part-load efficiency, runtime, network losses, and building pressure set points, and is compared with traditional multi-boiler controls. To support this exploration, new dynamic Modelica models are developed. In addition, methods to reduce modeling complexities are incorporated, enhancing their suitability for practical applications. A holistic district-wide analysis using a real university case study demonstrates a 4.7% fuel savings by lowering boiler operational pressure from 900 kPa to 600 kPa, along with a 13.3% reduction in condensation losses across the distribution network. Furthermore, the load allocation approach results in a 13.1% reduction in fuel consumption during peak winter periods and 15.3% during shoulder periods, with corresponding decreases in carbon emissions and fuel costs. This approach can also save maintenance costs by reducing the boiler runtime by 49.6%. In conclusion, this research underscores the benefits of retrofitting aging steam district heating systems, offering immediate operational improvements by enhancing efficiency, meeting regulatory compliance, and extending infrastructure lifespans while delaying costly overhauls.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Force-Modulated Selectivity of the Rhodium-Catalyzed Hydroformylation of 1-Alkenes

Here we quantify the regio- chemo-, and enantioselectivity of the hydroformylation of 1-alkenes catalyzed by Rh(I) complexes containing chiral macrocyclic biaryl bis(phosphine) ligands as a function of mechanical force applied to the biaryl backbone of these ligands. Extension forces increase the regioselectivity of the hydroformylation of 1-octene from l:b = 1.2 to 1.9 and decrease the chemoselectivity for C 9 aldehyde relative to 2-octene from 96% to 20% across a ~230 pN change in applied force. Extension forces decrease the enantioselectivity of the hydroformylation of styrene from 28% to 18% ee across a similar force range. Variable temperature 1 H and 31 P NMR analysis of five-coordinate rhodium complexes of the form (P–P)Rh(CO) 2 H established structures in which the bisphosphine ligand occupies one equatorial and one apical position in the trigonal bipyramid with rapid interconversion of the equatorial and apical phosphorus atoms ($ΔG^‡_{223 \text{K}}$ = 9.6–9.9 kcal/mol). Neither the isomeric composition nor fluxional behavior of these complexes was detectably perturbed across a ~300 pN change in applied force.

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

An improved guess for the variational calculation of charge-transfer excitations in large systems

Ab initio quantum-chemical methods that perform well for computing the electronic ground state are not straightforwardly transferable to electronically excited states, particularly in large molecular systems. Wave function theory offers high accuracy, but is often prohibitively expensive. Methods based on time-dependent density functional theory (TD-DFT) are crucially sensitive to the chosen exchange-correlation functional (XCF) parameterization, and system-specific tuning protocols were therefore proposed to address the method's robustness. Methods based on the variational relaxation of the excited-state electron density showcased promising results for the calculation of charge-transfer excitations, but the complex shape of the electronic hypersurface makes convergence to a specific excited state much more difficult than for the ground state when standard variational techniques are applied. We address the latter aspect by providing suitable initial guesses, which we obtain by two separate constrained algorithms. Combined with the squared-gradient minimization algorithm for all-electrons relaxation in a freeze-and-release scheme (FRZ-SGM), we demonstrate that orbital-optimized density functional theory (OO-DFT) calculations can reliably converge to the charge-transfer states of interest even for large molecular systems. We test the FRZ-SGM method on a phenothiazine-anthraquinone CT excitation in a supramolecular Pd(II) coordination cage complex as a function of the cage conformation. This compound has been studied experimentally prior to our work. We compare this freeze-and-release scheme to two XCF reparameterizations, which were recently proposed as low-cost TD-DFT-based alternatives to variational methods. Two dye-semiconductor complexes, which were previously investigated in the context of photovoltaic applications, serve as a second example to investigate the convergence and stability of the FRZ-SGM approach. Our results demonstrate that FRZ-SGM provides reliable convergence for charge-transfer excited states and avoids variational collapse to lower-lying electronic states, whereas time-dependent DFT calculations with an adequate tuning procedure for the range-separation parameter provide a computationally efficient initial estimate of the corresponding energies, with a computational cost comparable to that of configuration-interaction singles (CIS) calculations.

Bogo, Nicola↗
DOE OSTI2025

Inverse Trans Influence and Uranium-Arene σ-Bonding Drive Molecular Geometry: Ligand Modification from Hard to Soft Flips the Oxide from Axial to Equatorial

A rare example of an equatorially bound terminal uranium(V) oxo complex in a chelating sulfur-based ligand environment, namely [(mes( Me,Ad ArS) 3 )U V (O eq )(THF)] (2), is presented. Octahedrally coordinated 2 is obtained by reaction of the mesitylene-anchored tris-thiophenolate-coordinated uranium(III) complex [U III ((SAr Ad,Me ) 3 mes)] (1) with the oxygen-atom transfer reagent N 2 O. The observed, equatorially bound oxo ligand in 2 is in stark contrast to its known tris-aryloxide analog, [(mes( Me,Ad ArO) 3 )U V (O ax )(THF)] (A), where the oxo ligand occupies the typically observed axial coordination site. Complexes 1 and 2 are characterized by single-crystal X-ray diffraction analyses and spectroscopic and magnetochemical methods, including 1 H NMR, UV/vis/NIR electronic absorption, as well as EPR spectroscopy and SQUID magnetometry, thus confirming the C S symmetry and the pentavalent oxidation state of 2. Encompassing quantum chemical calculations (DFT and CASPT2) on 2 and its tris-phenolate analog A, support and rationalize the structural and electronic differences. The molecular orbital pictures show that a stabilizing σ-bonding interaction arising from the U–O eq inverse trans influence (ITI) is present in 2 but missing in A. In 2, the sulfur 3p orbitals are closer in energy to the uranium 5f manifold than the arene π-system, leading to an ITI, while U–arene σ- or δ-bonding is not observed. Although the arene orbitals remain separated from the uranium 5f orbitals in A, the absence of an ITI allows the arene a 2u orbital to engage in a σ-type interaction with the metal. Thus, incorporating a tris-thiophenolate to an arene anchor introduces a new design concept in molecular f-element chemistry. This approach stabilizes an equatorially bound U(V) oxo center, contrasting with its tris-phenolate counterpart, where oxo coordination is axial. The observed geometric divergence, driven by competing ITI and U–arene interactions, not only tunes electronic structure but also leads to differentiated reactivity: only the phenolate analogs activate H 2 O, while the thiolates do not.

Hydrocarbons↗
DOE OSTI2026

Covalency of M–N Bonds in Isomorphous Lanthanide and Actinide 5-(2-Pyridyl)-1H-tetrazolate Complexes

Experimental and computational analyses of [M(pdtz) 3 (H 2 O) 3 ]·3.5H 2 O (M 3+ = Pu 3+ −Cm 3+ , La 3+ −Nd 3+ , and Sm 3+ −Ho 3+ , pdtz− = 5-(2-pyridyl)-1H-tetrazolate) were conducted to understand potential differences in bonding between lanthanide and actinide complexes with a N-donor ligand. Structural analyses show that the An−N bond distances in the Pu 3+ , Am 3+ , and Cm 3+ complexes are within error of one another. Whereas in the lanthanide series, there is a nearly linear decrease in the Ln−N bond lengths from La 3+ to Ho 3+ (excluding Pm 3+ ). The An−N bond lengths are ∼0.015 Å shorter than their similarly-sized lanthanide analogs, in agreement with computational results that suggest greater covalent character in these bonds versus those with lanthanides. QTAIM analysis indicates that the An−N orbital mixing remains essentially unchanged from Pu 3+ to Cm 3+ , consistent with the nearly identical An−N bond lengths. However, upon deconvolution of the NLMOs into orbital compositions, the metal orbital contributions to An−N bonding decreases slightly overall wherein the 6d involvement remains constant, 7s involvement slightly increases, and 5f participation decreases. The molecular orbital energy diagram indicates that energy degeneracy between the 5f metal and 2p ligand orbitals increases from Pu 3+ to Cm 3+ and counteracts the contraction of the 5f orbtials. Together with prior reports of decreasing energy degeneracy between 5f and 3p orbitals from Np 3+ to Cf 3+ , these observations provide guidance on understanding how chemical bonding evolves in the actinide series.

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

Quantitative Decoupling of Oxygen-Redox and Manganese-Redox Voltage Hysteresis in a Cation-Disordered Rock Salt Cathode

Pronounced voltage hysteresis in Li-excess cathode materials is commonly thought to be associated with oxygen redox. However, these materials often possess overlapping oxygen and transition-metal redox, whose contributions to hysteresis between charge and discharge are challenging to distinguish. In this work, a two-step aqueous redox titration is developed with the aid of mass spectrometry (MS) to quantify oxidized lattice oxygen and Mn 3+/4+ redox in a representative Li-excess cation-disordered rock salt—Li 1.2 Mn 0.4 Ti 0.4 O 2 (LMTO). Two MS-countable gas molecules evolve from two separate titrant-analyte reactions, thereby allowing Mn and O redox capacities to be decoupled. The decoupled O and Mn redox coulombic efficiencies are close to 100% for the LMTO cathode, indicating high charge-compensation reversibility. As incremental Mn and O redox capacities are quantitatively decoupled, each redox voltage hysteresis is further evaluated. Overall, LMTO voltage hysteresis arises not only from an intrinsic charge-discharge voltage mismatch related to O redox, but also from asymmetric Mn-redox overvoltages. The results reveal that O and Mn redox both contribute substantially to voltage hysteresis. In conclusion, this work further shows the potential of designing new analytical workflows to experimentally quantify key properties, even in a disordered material having complex local coordination environments.

25 ENERGY STORAGE↗
DOE OSTI2024

A multilayer network analysis of Alzheimer's disease pathogenesis: Roles for p‐tau, synaptic peptides, and physical activity

INTRODUCTION: In the aging brain, cognitive abilities emerge from the coordination of complex pathways arising from a balance between protective lifestyle and environmental factors and accumulation of neuropathologies. METHODS: As part of the Rush Memory and Aging Project (n = 440), we measured accelerometer-based actigraphy, cognitive performance, and after brain autopsy, selected reaction monitoring mass spectrometry. Multilevel network analysis was used to examine the relationships among the molecular machinery of vesicular neurotransmission, Alzheimer's disease (AD) neuropathology, cognition, and late-life physical activity. RESULTS: Synaptic peptides involved in neuronal secretory function were the most influential contributors to the multilayer network, reflecting the complex interdependencies among AD pathology, synaptic processes, and late-life cognition. Older adults with lower physical activity evidenced stronger adverse relationships among phosphorylated tau peptides, markers of synaptic integrity, and tangle pathology. DISCUSSION: Network-based approaches simultaneously model interdependent biological processes and advance understanding of the role of physical activity in age-associated cognitive impairment. Highlights: Network-based approaches simultaneously model interdependent biological processes. Secretory synaptic peptides were influential contributors to the multilayer network. Older adults with lower physical activity had adverse relationships among pathology. There was interdependence among phosphorylated tau, synaptic integrity, and tangles. Network methods elucidate the role of physical activity in cognitive impairment.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Sulfonate Group Improves the Solubility and Electrocatalytic Performance of Ru‐Based bda‐ and pda‐Type Water Oxidation Catalysts under Neutral Conditions

Four ruthenium water oxidation catalysts that bear carboxylate and sulfonate groups in the active site have been synthesized and analyzed for their catalytic activity. Here, the sulfonate‐containing catalysts show higher electrochemical activity in pH 7 phosphate buffer with 4 times larger catalytic current, improved durability with sacrificial oxidant, and increased solubility compared to similar species containing two carboxylate groups. Density functional theory calculations suggest that the sulfonate group provides a more favorable geometry for water nucleophilic attack, which is both the energetically most favorable mechanism calculated and the experimentally predicted mechanism under electrochemical conditions. Further experimental studies have been performed to show that under certain conditions catalysts can perform well electrochemically under pH conditions as low as 1.6 and that various structural components can greatly change solubility and catalytic operation.

14 SOLAR ENERGY↗
DOE OSTI2022

Multiscale simulation of integrated energy system and electricity market interactions

Accelerating the deep decarbonization of the world's electric grids requires the coordination of complex energy systems and infrastructures across timescales from seconds to decades. Here, we present a new multiscale simulation framework that integrates process- and grid-centric modeling paradigms to better design, operate, and control integrated energy systems (IESs), which combine multiple technologies, in wholesale energy markets. Traditionally, IESs are analyzed with a process-centric paradigm such as levelized cost of electricity (LCOE) or annualized net revenue, ignoring important interactions with electricity markets. This framework explicitly models the complex interactions between an IES's bidding, scheduling, and control decisions and the energy market's clearing and settlement processes, while incorporating operational uncertainties. Through two case studies, we show the importance of understanding and quantifying complex resource-grid interactions. In case study 1, we demonstrate that optimized bidding from one resource shifts the profit distribution for all energy systems in the market. This result suggests new and more flexible IES technologies can disrupt the economics of all market participants, possibly leading to accelerated retirements of less flexible resources. Interestingly, the optimized bidding has little impact on grid-level aggregate statistics, such as total generation costs and renewable penetration rate. While aggregate modeling strategies may remain valid under some IES adoption scenarios for analysis focused on regional outcomes, direct comparisons of IES technologies at specific locations without considering these interactions may lead to misleading or incorrect conclusions. In case study 2, we consider the design and flexible operation of IESs that hybridize conventional generators with energy storage. Through a sensitivity analysis, we find that as the size of the storage system increases, the total number of start-ups for coal- and natural gas-based IESs reduced by 25% and 33.6%, and the total thermal generator ramping (i.e., mileage) reduced by 86.5% and 62.5%, respectively. This shows the primary benefit of storage may not be reduced operational costs (which do not change significantly) but fewer start-ups and less ramping, which may greatly simplify the design, operation, and control of carbon capture systems. The new modeling and optimization capabilities from this work enable the coupling of rigorous, dynamic process models with grid-level production cost models to quantitatively identify the nuanced interdependencies across these vast timescales that must be addressed to realize clean, safe, and secure energy production. Moreover, the proposed general multiscale simulation framework is applicable to all IES technologies and can be easily extended to consider other energy carriers (e.g., hydrogen, ammonia) and energy infrastructures (e.g., natural gas pipelines).

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

A comparative study of the effectiveness of empirical potentials for molecular dynamics simulations of borosilicate glasses

Despite their practical importance, atomistic modeling of B 2 O 3 containing glasses have been challenging due to the lack of reliable empirical potentials. Fortunately, a few recent developments have shown promises to simulate these glasses where the boron coordination has complex and non-linear dependence on glass composition. This work aims to provide an evaluation of the effectiveness of three recently developed potentials by a systematic study of a series (~20) of sodium borosilicate glasses with constant K (ratio of [SiO 2 ]/[B 2 O 3 ]) =2 and varying R (ratio of [Na 2 O]/[B 2 O 3 ] ranging from 0.1 to 4) values and several sodium boroaluminosilicate glasses. A comparative assessment was established on the basis of the short- and medium-range structure features, such as boron N 4 values, total correlation functions, bond angle distribution, oxygen speciation, and mechanical properties using experimental or well-established models as criteria. Furthermore, this work provides insights on the choices of empirical potentials for MD simulations of borosilicate glasses and gives directions of future potential development and refining.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗