Search NASA⌕ Search

SEARCH · Search NASA

Results for “Density functional theory simulations”

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

Ab initio studies of the impact of the Debye-Waller factor on the structural and dynamical properties of amorphous semiconductors: The case of a-Si

This paper presents a first-principles study of the Debye-Waller factor and the Debye temperature for amorphous silicon (a-Si) from lattice-dynamical calculations and direct molecular-dynamics simulations using density-functional theory (DFT). The effects of temperature and structural disorder on the intensity of the diffraction maxima and the vibrational mean-square displacement (MSD) of Si atoms are studied in the harmonic approximation, with particular emphasis on the bond-length disorder, the presence of coordination defects, and microvoids in a-Si networks. It has been observed that the MSDs associated with tetrahedrally-bonded Si atoms are considerably lower than their dangling-bond counterparts – originating from isolated and vacancy-induced clustered defects – and those on the surface of microvoids, leading to an asymmetric non-gaussian tail in the distribution of atomic displacements. An examination of the effect of anharmonicity on the MSD at high temperatures using direct ab initio molecular-dynamics simulations (without the harmonic approximation) suggests that the vibrational motion in a-Si is practically unaffected by anharmonic effects at temperatures below 400 K, as far as the present DFT calculations are concerned. The Debye temperature of a-Si is found to be in the range of 488–541 K from specific-heat and MSD calculations using first-principles lattice-dynamical calculations in the harmonic approximation, which matches closely with the experimental value of 487–528 K obtained from specific-heat measurements of a-Si at low temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Evaluation of Thermal Scattering Law for U 3 Si 2 Using Temperature-Dependent Phonon Density of States

The thermal scattering law for uranium silicide, a candidate material for high-density nuclear fuel, was evaluated from the temperature-dependent phonon density of states obtained using the quasi-harmonic approximation. For uranium silicide model systems of various volumes near the equilibrium volume, Hellmann-Feynman forces were calculated by performing ab initio lattice dynamics simulations based on density functional theory. Subsequently, the quasi-harmonic approximation was applied by interpolating the results of the harmonic approximation for each volume system. As a sensitivity test, the phonon density of states was calculated for the structures with temperature-dependent volumes at 298 K and 1000 K, obtained by the quasi-harmonic approximation, and the corresponding thermal neutron scattering cross section was evaluated. The temperature-dependent volume and heat capacity of uranium silicide obtained by the quasi-harmonic approximation showed reasonable agreement with experimental data while lattice thermal expansion showed deviation from experimental observations. Thermal scattering cross sections evaluated from 0 K and temperature-dependent phonon density of states exhibited differences of up to 4.1% at 1000 K. This study examines the impact of temperature on the thermal scattering law, employing ab initio lattice dynamics simulations through the quasi-harmonic approximation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Intermetallics Based on Sodium Chalcogenides Promote Stable Electrodeposition–Electrodissolution of Sodium Metal Anodes

Sodiophilic micro-composite films of sodium-chalcogenide intermetallics (Na 2 Te and Na 2 S) and Cu particles are fabricated onto commercial copper foam current collectors (Na 2 Te@CF and Na 2 S@CF). For the first time a controllable capacity thermal infusion process is demonstrated. Enhanced wetting by the metal electrodeposition leads to state-of-the-art electrochemical performance. For example, Na 2 Te@CF-based half-cells demonstrate stable cycling at 6 mA cm -2 and 6 mAh cm -2 , corresponding to 54 µm of Na electrodeposited/electrodissolved by geometric area. Sodium metal batteries with Na 3 V 2 (PO 4 ) 3 cathodes are stable at 30C (7 mA cm -2 ) and for 10 000 cycles at 5C and 10C. Cross-sectional cryogenic focused ion beam (cryo-FIB) microscopy details deposited and remnant dissolved microstructures. Sodium metal electrodeposition onto Na 2 Te@CF is dense, smooth, and free of dendrites or pores. On unmodified copper foam, sodium grows in a filament-like manner, not requiring cycling to achieve this geometry. Substrate–metal interaction critically affects the metal–electrolyte interface, namely the thickness and morphology of the solid electrolyte interphase. Density functional theory and mesoscale simulations provide insight into support-adatom energetics, nucleation response, and early-stage morphological evolution. On Na 2 Te sodium atomic dispersion is thermodynamically more stable than isolated clusters, leading to conformal adatom coverage of the surface.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Propane Dehydrogenation on Pt x Zn y Active Sites in Silicalite‐1

Abstract The improvement of Pt‐based catalysts for propane dehydrogenation (PDH) has progressed by recent investigations that have identified Zn as a promising promoter for Pt subnanometer catalysts. It is desirable to gain insights into the structure, stability, and activity of such active sites and the factors that influence them, such as Zn : Pt ratio, Pt coordination and nuclearity. Here, we employ density functional theory and microkinetic simulations to investigate the stability of Pt x Zn y ( x =1–3, y=0–3) active sites grafted on silanols of Silicalite‐1 and the PDH activity of Pt. We find that the coordination of a Pt atom to a nest of grafted Zn(II) atoms increases the stability of the Pt 1 Zn y sites, whose activity is similar for y=0–2 and drops dramatically for y>2. We further demonstrate, via linear scaling relations and microkinetic simulations, that the turnover frequency obeys a volcano law as a function of propylene binding strength. The Pt 2 Zn 1 and Pt 3 Zn 1 sites are stable and exhibit activity similar to Pt 1 Zn 2 , but only Pt 1 Zn 2 manifests reaction kinetics consistent with experimental data, strongly suggesting the active site composition in the synthesized catalyst samples. The methodology presented here suggests a general strategy for deducing active site information such as composition through simple kinetic experiments.

Liu, Yilang↗
DOE OSTI2024

Propane Dehydrogenation on Pt x Zn y Active Sites in Silicalite‐1

Abstract The improvement of Pt‐based catalysts for propane dehydrogenation (PDH) has progressed by recent investigations that have identified Zn as a promising promoter for Pt subnanometer catalysts. It is desirable to gain insights into the structure, stability, and activity of such active sites and the factors that influence them, such as Zn : Pt ratio, Pt coordination and nuclearity. Here, we employ density functional theory and microkinetic simulations to investigate the stability of Pt x Zn y ( x =1–3, y=0–3) active sites grafted on silanols of Silicalite‐1 and the PDH activity of Pt. We find that the coordination of a Pt atom to a nest of grafted Zn(II) atoms increases the stability of the Pt 1 Zn y sites, whose activity is similar for y=0–2 and drops dramatically for y>2. We further demonstrate, via linear scaling relations and microkinetic simulations, that the turnover frequency obeys a volcano law as a function of propylene binding strength. The Pt 2 Zn 1 and Pt 3 Zn 1 sites are stable and exhibit activity similar to Pt 1 Zn 2 , but only Pt 1 Zn 2 manifests reaction kinetics consistent with experimental data, strongly suggesting the active site composition in the synthesized catalyst samples. The methodology presented here suggests a general strategy for deducing active site information such as composition through simple kinetic experiments.

Liu, Yilang↗
DOE OSTI2021

Solid‐state nuclear magnetic resonance of spin‐9/2 nuclei 115 In and 209 Bi in functional inorganic complex oxides

Abstract Indium and bismuth are technologically important elements, in particular as oxides for optoelectronic applications. 115 In and 209 Bi are both I = 9/2 nuclei with high natural abundances and moderately high frequencies but large nuclear electric quadrupole moments. Leveraging the quadrupolar interaction as a measure of local symmetry and polyhedral distortions for these nuclei could provide powerful insights on a range of applied materials. However, the absence of reported nuclear magnetic resonance (NMR) parameters on these nuclei, particularly in oxides, hinders their use by the broader materials community. In this contribution, solid‐state 115 In and 209 Bi NMR of three recently discovered quaternary bismuth or indium oxides are reported, supported by density functional theory calculations, numerical simulations, diffraction and additional multinuclear ( 27 Al, 69,71 Ga, and 121 Sb) solid‐state NMR measurements. The compounds LiIn 2 SbO 6 , BiAlTeO 6 , and BiGaTeO 6 are measured without special equipment at 9.4 T, demonstrating that wideline techniques such as the QCPMG pulse sequence and frequency‐stepped acquisition can enable straightforward extraction of quadrupolar tensor information in I = 9/2 115 In and 209 Bi even in sites with large quadrupolar coupling constants. Relationships are described between the NMR observables and local site symmetry. These are amongst the first reports of the NMR parameters of 115 In, 121 Sb, and 209 Bi in oxides.

Griffith, Kent J.↗
DOE OSTI2025

A generalized and adaptable tensor-contraction-based cluster expansion formalism for multicomponent solids

Density functional theory (DFT)-based simulations of materials have first-principles accuracy, but are very computationally expensive. For simulating various properties of multi-component alloys, the cluster expansion (CE) technique has served as the standard workaround to improve computational efficiency. However, the standard CE technique is difficult to extend to exotic and/or low-symmetry lattices, often implemented via iteration over particular cluster types, which must be enumerated per lattice structure. In this work, we introduce the tensor cluster expansion (TCE), implemented in the open-source code tce-lib, which maps correlation functions to mixed tensor contractions, eliminating the need to iterate over cluster types and additionally making the calculation of correlation functions well-suited for massively parallel architectures like GPUs. We show that local interaction energies are an immediate consequence of the TCE formalism, yielding nearly $\mathcal{O}$(1) energy difference calculations. We then use this formalism to fit CE models for the TaW and CoNiCrFeMn systems, and use these models to respectively compute the enthalpy of mixing curve and Cowley short-range order parameters, showing excellent agreement with ground truth data.

Cluster expansion↗
DOE OSTI2023

Pressure-induced evolution of the lattice dynamics for selected UO 3 polymorphs

Uranium trioxide (UO 3 ) is a stable chemical form of uranium oxide with multiple polymorphs found throughout the nuclear fuel cycle. The pressure-induced changes in the structure and lattice dynamics of four of these polymorphs are simulated with density functional perturbation theory and analyzed. Two phases, α- and δ-UO 3 are found to exhibit an isotropic response to pressure and do not undergo any changes in coordination geometry up to ~40 GPa. In contrast, the other two phases investigated, β- and γ-UO 3 , exhibit an anisotropic response to pressure. Decomposition of the phonon eigenvectors allows us to assign specific pressure-induced structural changes to individual phonon modes. Here, this analysis has been performed on a per atom basis for the relatively simple α- and δ-UO 3 structures, which have one symmetrically unique uranium site, and on a per coordination environment basis for β- and γ-UO 3 , which have multiple U sites.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Stability and mobility of tungsten clusters on tungsten (110) surface: Ab initio and atomistic simulations

Quantifying the surface transport properties of tungsten (W) is of prime importance to understand the formation of nano-fuzz in fusion plasma-facing conditions. The stability and mobility of W adatom clusters (W n , n = 2-9) on the W(110) surface has been investigated by computer simulations, including ab initio calculations using density functional theory (DFT) and molecular statics (MS) simulations with multiple W interatomic potentials. The DFT results demonstrate that the sequential binding energy generally increases with number of W adatoms, except for the 5 th and 7 th W adatoms. The most common elemental migration steps of W n (n>2) clusters are observed to consist of monomer and dimer hops, while larger W n clusters can also diffuse by dissociation and recombination of smaller clusters. The threshold migration energy of W 9 is the highest, then followed by W 8 , W 4 , and W 6 , while W 3 , W 5 , and W 7 have similar migration energies. Compared to DFT, each interatomic potential evaluated overestimates the binding energies of W n clusters. Additionally, the embedded-atom potential developed by Juslin and Wirth adequately predicts the threshold migration energy of W n (n>2) clusters, although it predicts different underlying migration mechanisms. The results show that interaction mechanism between W adatoms controls the stability and mobility of W n clusters on the W(110) surface.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Net Electronic Charge as an Effective Electronic Descriptor for Oxygen Release and Transport Properties of SrFeO 3 -based Oxygen Sorbents

Perovskite oxides, as oxygen sorbents, exhibit excellent potential in thermochemical redox applications such as chemical looping air separation (CLAS), resulting from their excellent redox properties and high tunability. The structural and compositional flexibility of perovskites, while making these mixed oxides highly versatile, also poses challenges for their design and optimization. Moreover, the elevated operating temperature and dynamic oxygen partial pressure swings during redox processes add additional complications for oxygen sorbent design. From a fundamental standpoint, finding simple yet effective descriptors for the effect of cation substitutions on oxygen sorbents’ redox and oxygen transport properties are of significant importance for rational optimization of perovskite-based oxygen sorbents. In the present study, a series of A/B-site-doped SrFeO 3-δ perovskites are investigated using density functional theory (DFT) + U simulations. From these simulations, an effective electronic structural descriptor, the net electronic charge (Δe) of oxygen anions, is extracted based on its excellent correlations with vacancy formation energy, experimentally derived oxygen capacity, and oxygen vacancy migration barrier. Overall, these findings provide an effective tool to correlate the oxygen sorbents’ redox performances with a simple yet physically meaningful descriptor for rational design and optimization of perovskites in the context of CLAS.

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

Abnormally Low Activation Energy in Cubic Na3SbS4 Superionic Conductor

Inorganic Na-ion superionic conductors play a vital role in all-solid-state Na batteries that operate at room temperature. Sodium thioantimonate (Na3SbS4), a popular sulfide-based solid electrolyte, has attracted serious attention due to its advantages of high ionic conductivity at room temperature and impressive chemical stability under ambient conditions. Much research detailing Na3SbS4 focused on its synthetic approaches and interfacial stability against Na metal, yet, there is limited information elucidating a fundamental understanding of the Na- ion diffusion mechanisms in Na3SbS4 with different crystal structures (e.g., tetragonal and cubic). Herein, we combine real-time electrochemical impedance measurements with theoretical simulations based on density functional theory and in situ quasi-elastic neutron scattering to study the Na-ion conductive properties of Na3SbS4 during its phase transition from a tetragonal to cubic structure. Although there is a slight change in the lattice parameters, the energy barrier for Na-ion diffusion in the tetragonal structure was determined to be much larger (5-10 times) than that in the cubic structure from both theoretical and experimental perspectives. The high degree of symmetry in cubic Na3SbS4 leads to less interatomic correlations between Na and S(Sb) atoms, a shorter jump distance (2.85 angstrom), and a larger diffusion coefficient. This research provides insight into understanding the Na-ion diffusion in solid electrolytes with phase transitions and provides fundamental guidance for designing novel solid-state Na-ion conductors.

Zhang, Qian↗
DOE OSTI2022

Stabilization of Dinuclear Rhodium and Iridium Clusters on Layered Titanate and Niobate Supports

Atomically dispersed organometallic clusters can provide well-defined nuclearity of active sites for both fundamental studies as well as new regimes of activity and selectivity in chemical transformations. More recently, dinuclear clusters adsorbed onto solid surfaces have shown novel catalytic properties resulting from the synergistic effect of two metal centers to anchor different reactant species. Difficulty in synthesizing, stabilizing, and characterizing isolated atoms and clusters without agglomeration challenges allocating catalytic performance to atomic structure. Here, we explore the stability of dinuclear rhodium and iridium clusters adsorbed onto layered titanate and niobate supports using molecular precursors. Both systems maintain their nuclearity when Statistical analysis of HAADF-STEM images revealed that rhodium and iridium dimers had mean cluster-to-cluster distances very similar to what is expected from a random distribution of atoms over a large area, indicating that they are dispersed without aggregation. The stability of dinuclear rhodium clusters supported on titanate nanosheets was also investigated by X-ray absorption fine structure (EXAFS), DRIFTS, and first-principles calculations. Both X-ray absorption spectroscopy and HAADF-STEM simulations, guided by density functional theory (DFT)-optimized structure models, suggested that rhodium dimers adsorb onto the nanosheets in an end-on binding mode that is stable up to 100 degrees C under reducing conditions. Finally, this study highlights that crystalline nanosheets derived from layered metal oxides can be used as model supports to selectively stabilize dinuclear clusters, which could have implications for heterogeneous catalysis.

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

Photophysical Properties and Electronic Structure of Zinc(II) Porphyrins Bearing 0–4 meso -Phenyl Substituents: Zinc Porphine to Zinc Tetraphenylporphyrin (ZnTPP)

Six zinc(II) porphyrins bearing 0–4 meso-phenyl substituents have been examined spectroscopically and theoretically. Comparisons with previously examined free base analogues afford a deep understanding of the electronic and photophysical effects of systematic addition of phenyl groups in porphyrins containing a central zinc(II) ion versus two hydrogen atoms. Trends in the wavelengths and relative intensities of the absorption bands are generally consistent with predictions from time-dependent density functional theory calculations and simulations from Gouterman’s four-orbital model. These trends derive from a preferential effect of the meso-phenyl groups to raise the energy of the highest occupied molecular orbital. The calculations reveal additional insights, such as a progressive increase in oscillator strength in the violet-red (B-Q) absorption manifold with increasing number of phenyls. Progressive addition of 0–4 phenyl substituents to the zinc porphyrins in O 2 -free toluene engenders a reduction in the measured lifetime of the lowest singlet excited state (2.5–2.1 ns), an increase in the S 1 → S o fluorescence yield (0.022–0.030), a decrease in the yield of S 1 → T 1 intersystem crossing (0.93–0.88), and an increase in the yield of S 1 → S o internal conversion (0.048–0.090). The derived rate constants for S 1 decay reveal significant differences in the photophysical properties of the zinc chelates versus free base forms. Here, the unexpected finding of a larger rate constant for internal conversion for zinc chelates versus free bases is particularly exemplary. Collectively, the findings afford fundamental insights into the photophysical properties and electronic structure of meso-phenylporphyrins, which are widely used as benchmarks for tetrapyrrole-based architectures in solar energy and life sciences research.

14 SOLAR ENERGY↗
DOE OSTI2020

Passivating Nucleobases Bring Charge Transfer Character to Optically Active Transitions in Small Silver Nanoclusters

DNA-wrapped silver nanoclusters (DNA–AgNCs) are known for their efficient luminescence. However, their emission is highly sensitive to the DNA sequence, the cluster size, and its charge state. To get better insights into photophysics of these hybrid systems, simulations based on density functional theory (DFT) are performed. Our calculations elucidate the effect of the structural conformations, charges, solvent polarity, and passivating bases on optical spectra of DNA–AgNCs containing five and six Ag atoms. It is found that inclusion of water in calculations as a polar solvent media results in stabilization of nonplanar conformations of base-passivated clusters, while their planar conformations are more stable in vacuum, similar to the bare Ag5 and Ag6 clusters. Cytosines and guanines interact with the cluster twice stronger than thymines, due to their larger dipole moments. In addition to the base–cluster interactions, hydrogen bonds between bases notably contribute to the structure stabilization. While the relative intensity, line width, and the energy of absorption peaks are slightly changing depending on the cluster charge, conformations, and base types, the overall spectral shape with five well-resolved bands at 2.5–5.5 eV is consistent for all structures. Independent of the passivating bases and the cluster size and charge, the low energy optical transitions at 2.5–3.5 eV exhibit a metal to ligand charge transfer (MLCT) character with the main contribution emerging from Ag-core to the bases. Cytosines facilitate the MLCT character to a larger degree comparing to the other bases. However, the doublet transitions in clusters with the open shell electronic structure (Ag5 and Ag6 + ) result in appearance of additional red-shifted (<2.5 eV) and optically weak band with negligible MLCT character. The passivated clusters with the closed shell electronic structure (Ag5 + and Ag6) exhibit higher optical intensity of their lowest transitions with much higher MLCT contribution, thus having better potential for emission, than their open shell counterparts.

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

Predicted Reaction Mechanisms, Product Speciation, Kinetics, and Detonation Properties of the Insensitive Explosive 2,6-Diamino-3,5-dinitropyrazine-1-oxide (LLM-105)

2,6-Diamino-3,5-dinitropyrazine-1-oxide (LLM-105) is a relatively new and promising insensitive high-explosive (IHE) material that remains only partially characterized. IHEs are of interest for a range of applications and from a fundamental science standpoint, as the root causes behind insensitivity are poorly understood. In this work, we adopt a multitheory approach based on reactive molecular dynamic simulations performed with density functional theory, density functional tight-binding, and reactive force fields to characterize the reaction pathways, product speciation, reaction kinetics, and detonation performance of LLM-105. We compare and contrast these predictions to 1,3,5-triamino-2,4,6-trinitrobenzene (TATB), a prototypical IHE, and 1,3,5,7-tetranitro-1,3,5,7-tetrazoctane (HMX), a more sensitive and higher performance material. The combination of different predictive models allows access to processes operative on progressively longer timescales while providing benchmarks for assessing uncertainties in the predictions. We find that the early reaction pathways of LLM-105 decomposition are extremely similar to TATB; they involve intra- and intermolecular hydrogen transfer. Additionally, the detonation performance of LLM-105 falls between that of TATB and HMX. We find agreement between predictive models for first-step reaction pathways but significant differences in final product formations. Predictions of detonation performance result in a wide range of values, and one-step kinetic parameters show the similar reaction rates at high temperatures for three out of four models considered.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

CO 2 Capture Characteristics of Hyperbranched Poly(alkylene imine): A Molecular Dynamics Simulation Approach

This study explores the CO 2 capture characteristics of hyperbranched poly(ethylenimine) (HB-PEI) and poly- (propyleneimine) (HB-PPI) through molecular dynamics simulations using density functional theory-calibrated force fields. Key features such as density, free volume, glass transition temperature, CO 2 /H 2 O distribution, and molecular diffusion are systematically investigated to elucidate structure−function relationships under dry and hydrated conditions. HB-PEI demonstrates a slightly higher density and lower free volume compared to HB-PPI yet shows superior CO 2 capture due to the high amine concentration. Glass transition analysis indicates a higher thermal mobility in HBPEI, enhancing the CO 2 diffusivity. Pair correlation and coordination analyses confirm a stronger affinity of CO 2 with primary and secondary amines, particularly in hydrated environments where water competes with CO 2 for binding sites. Despite its more compact structure, HB-PEI outperformed HB-PPI in CO 2 and H 2 O transport, as confirmed by higher diffusion coefficients across all hydration levels. These findings highlight a critical balance among polymer architecture, amine accessibility, and hydration in designing next-generation solid amine sorbents for efficient direct air capture applications.

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

Insights into the Electrochemical Oxidation and Reduction of Nickel Oxide Surfaces

Surface oxidation/reduction processes, driven by varying electrochemical potentials, can substantially impact catalyst effectiveness and, consequently, electrolyzer performance. Here, this study combines theoretical and experimental approaches to explore the surface redox behavior of nickel oxides, which are cost-effective and efficient catalysts for many electrochemical reactions. Surface Pourbaix diagrams for three different phases of nickel oxides, i.e., nickel hydroxide (Ni(OH) 2 ), nickel oxyhydroxide (NiOOH), and nickel dioxide (NiO 2 ), were constructed using density functional theory-based simulations. Various experimental methods, including cyclic voltammetry, in situ Raman spectroscopy, and electrochemical titration, were employed to probe the surface redox processes of nickel oxide thin films. Our findings indicate that the ABAB stacking sequence of Ni(OH) 2 lacks stability under oxidizing conditions to host the surface oxidation (deprotonation) events, while the AABBCC stacking sequence of NiOOH is energetically favorable due to the presence of interlayer hydrogen bonding. Rapid charge transfer facilitated by interlayer hydrogen bonding accounts for the higher reactivity of partially oxidized/reduced NiOOH (001) surfaces compared to Ni(OH) 2 (001) and NiO 2 (001) surfaces with the same stoichiometry, where interlayer hydrogen bonding is absent. Insights presented in this work can offer guidelines for optimizing operational conditions and tailoring the surface structures and oxidation states of nickel oxides to enhance performance in applications such as electrocatalysis and supercapacitors.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Effective Storage of Electrons in Water by the Formation of Highly Reduced Polyoxometalate Clusters

Aqueous solutions of polyoxometalates (POMs) have been shown to have potential as high-capacity energy storage materials due to their potential for multi-electron redox processes, yet the mechanism of reduction and practical limits are currently unknown. Herein, we explore the mechanism of multi-electron redox processes that allow the highly reduced POM clusters of the form {MO 3 } y to absorb y electrons in aqueous solution, focusing mechanistically on the Wells-Dawson structure X 6 [P 2 W 18 O 62 ], which comprises 18 metal centers and can uptake up to 18 electrons reversibly (y = 18) per cluster in aqueous solution when the countercations are lithium. This unconventional redox activity is rationalized by density functional theory, molecular dynamics simulations, UV-vis, electron paramagnetic resonance spectroscopy, and small-angle X-ray scattering spectra. These data point to a new phenomenon showing that cluster protonation and aggregation allow the formation of highly electron-rich meta-stable systems in aqueous solution, which produce H 2 when the solution is diluted. Finally, we show that this understanding is transferrable to other salts of [P 5 W 30 O 110 ] 15- and [P 8 W 48 O 184 ] 40- anions, which can be charged to 23 and 27 electrons per cluster, respectively.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗