Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Level density”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Optimized Auxiliary Functions for Robust Mitigation of Finite-Size Errors in Periodic Hybrid Density Functional Theory

When calculating properties of periodic systems at the thermodynamic limit (TDL), the dominant source of finite size error (FSE) arises from the long-range Coulomb interaction, and can manifest as a slowly converging quadrature error when approximating an integral in the reciprocal space by a finite sum. The singularity subtraction (SS) method offers a systematic approach for reducing this quadrature error and thus the FSE. Here, in this work, we first investigate the performance of the SS method in the simplest setting, aiming at reducing the FSE in exact exchange calculations by subtracting the Coulomb contribution with a single, adjustable Gaussian auxiliary function. We demonstrate that a simple fitting method can robustly estimate the optimal Gaussian width and leads to rapid convergence toward the TDL. Furthermore, we suggest new forms of the auxiliary function, whose optimal parameters could also be determined through least-squares fitting. For a range of semiconductors and insulators, the proposed auxiliary functions achieve robust, millihartree-level accuracy in hybrid density functional theory calculations, including cases with sparse k-meshes and large basis sets.

Quiton, Stephen Jon [University of California, Ber↗
DOE OSTI2024

Electromagnetic Energy-Assisted Thermal Conversion of Fossil-Based Hydrocarbons to Low-Cost Hydrogen

The goal of this project was to develop and optimize catalysts for methane decomposition, particularly focusing on regeneration via an electromagnetic energy-assisted mechanism, to produce hydrogen more cost-effectively compared to electrolysis routes. To achieve this goal, the project pursued several key objectives. The project began with the preparation and testing of various catalysts. A nickel-silica based catalyst was identified as the most promising material for the pyrolysis of methane into carbon and hydrogen. Kinetic parameters for methane decomposition were determined, aiding in computational modeling efforts. Structured catalysts were investigated, highlighting the need for frequent cleaning or regeneration to maintain performance, with methane conversion rates exceeding 70% in tube furnace tests. Computational fluid dynamics modeling was employed to optimize reactor designs and electrode angles, leading the project team to propose a multi-compartment thermal conversion system for larger setups. This modeling work was important in understanding reaction characteristics, carbon deposition rates, and temperature profiles under various conditions. A bench-scale reactor system was assembled to evaluate catalyst regeneration using electromagnetic energy-assisted mechanisms. Experiments demonstrated the potential for carbon removal, though further optimization is needed. The carbon produced from the methane conversion process was evaluated for potential use in lithium-ion battery electrodes. The carbon exhibited properties similar to commercially available high-purity multi-walled carbon nanotubes and nanofibers, with a carbon content greater than 95%. Coin cell batteries assembled with this carbon showed that lower replacement levels (10% to 33%) outperformed the control group, improving specific capacity density and stability. However, higher replacement levels (100%) demonstrated poorer performance, suggesting that excessive carbon substitution negatively impacts battery performance. These findings indicate the potential marketability of the produced carbon as a component in lithium-ion batteries, though further testing is necessary to confirm long-term advantages and disadvantages associated with the use of the carbon product. These results however justified further technoeconomic assessments to determine if the process can provide low-cost hydrogen. The economic feasibility and technical performance of methane decomposition for hydrogen production were assessed, focusing on three plant configurations: 100E (electrically heated), 100C (combustion heated using produced hydrogen), and PE-Hybrid (a combination of pyrolysis (indicating decomposition) and electrolysis). The Levelized Cost of Hydrogen (LCOH) for the pyrolysis configurations was found to be approximately 25% lower than that of electrolysis. The 100E configuration had the lowest LCOH at $\$$3.12/kg. Including carbon product sales significantly improved the economics, with the 100C configuration achieving a negative LCOH of -$\$$0.35/kg. The PE-Hybrid configuration was not economically advantageous compared to pure pyrolysis plants due to its complexity and additional equipment requirements. Ultimately, methane pyrolysis presents a viable method for near carbon dioxide-free hydrogen production, with significant economic advantages over electrolysis, especially when considering the sale of carbon byproducts. The 100E and 100C configurations showed the most promise, with the choice between them ultimately depending on the prices of power and natural gas. In conclusion, this technology has the potential to lower hydrogen production costs by leveraging the methane decomposition process with the sale of valuable carbon byproducts. By optimizing catalyst performance and integrating electromagnetic energy-assisted regeneration, the process can achieve higher efficiency and economic viability, making it a competitive alternative to traditional hydrogen production methods.

08 HYDROGEN↗
DOE OSTI2022

Strain-retardant coherent perovskite phase stabilized Ni-rich cathode

The use of state-of-the-art Ni-rich layered oxides (LiNi x Co y Mn 1-x-y O 2 , x > 0.5) as the cathode material for lithium-ion batteries can push the energy and power density to a higher level than is currently available. However, volume variation associated with anisotropic lattice strain and stress that is being developed during lithium (de) intercalation induces severe structural instability and electrochemical decay of the cathode materials, which is amplified further when the battery is operating at a high voltage (above 4.5 V), which is essential for unlocking its high energy. Even after much effort by the research community, an intrinsic strain-retardant method for directly alleviating the continuous accumulation of lattice strain remains elusive. By introducing a coherent perovskite phase into the layered structure functioning as a ‘rivet’, we significantly mitigate the pernicious structural evolutions by a pinning effect. The lattice strain evolution in every single cycle is markedly reduced by nearly 70% when compared with conventional materials, which significantly enhances morphological integrity leading to a notable improvement in battery cyclability. This strain-retardant approach broadens the perspective for lattice engineering to release the strain raised from lithium (de)intercalation and paves the way for the development of high-energy-density cathodes with long durability.

25 ENERGY STORAGE↗
DOE OSTI2022

The impacts of charge transfer, localization, and metallicity on hydrogen retention and transport capacity

Solid state hydrides such as early transition metal hydrides are of inestimable importance for the future of hydrogen energy and are actively being investigated for energy conversion and storage applications such as fuel cells, solid-state batteries and neutron moderators. The retention and transport behavior of hydrogen in these hydrides has a huge role on the extended performance of components. While early transition-metal-based compounds exhibit many peculiar properties due to their unique correlated electronic signatures arising from -orbital electrons, the fundamental chemistry and transport behavior of hydrogen in such hydrides is not well understood. In the present work, using density functional theory, a highly intricate bonding feature is revealed through the theoretical investigation of the electronic structure of early transition metal hydrides YH 2 and ZrH 2 . In particular, a pronounced charge transfer from the transition element to H, results in localized electron densities at deep energy levels. The interplay between intrinsic charge transfer, charge localization, and metallicity in YH 2 and ZrH 2 leads to strong chemical bonding between metal and hydrogen atoms and large energy barriers for the migration of hydrogen vacancies. Specifically, hydrogen vacancies are found to be stable in the neutral state due to electron screening effects, accompanied by substantially high migration barriers between 0.8–1.2 eV along different crystallographic directions. In contrast, recent literature shows the migration barrier for charged H vacancies in insulating s-block metal hydrides lie between 0.1–0.4 eV, which is suitable for fast conduction applications. This pivotal electron structure difference exploited between early transition metal hydrides and alkali/alkaline earth metal hydrides determines extended hydrogen retention in these early transition metal hydrides. Finally, this work explains fundamental differences between the electronic structure of s-block and d-block metal hydrides, and its impact on the mobility of hydrogen vacancies.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Structural and mechanical properties of monolayer amorphous carbon and boron nitride

Amorphous materials exhibit various characteristics that are not featured by crystals and can sometimes be tuned by their degree of disorder (DOD). Here, we report results on the mechanical properties of monolayer amorphous carbon (MAC) and monolayer amorphous boron nitride (maBN) with different DOD. The pertinent structures are obtained by kinetic-Monte-Carlo (kMC) simulations using machine-learning potentials (MLP) with density-functional-theory (DFT)-level accuracy. An intuitive order parameter, namely the areal fraction F x occupied by crystallites within the continuous random network, is proposed to describe the DOD. We find that F x captures the essence of the DOD: Samples with the same F x but different sizes and arrangements of crystallites, obtained using two distinct kMC procedures, have virtually identical radial distributions functions as well as bond-length and bond-angle distributions. Furthermore, by simulating the fracture process with molecular dynamics, we found that the mechanical responses of MAC and maBN before fracture are mainly determined by F x and are insensitive to the sizes and specific arrangements and to some extent the numbers and area distributions of the crystallites. The behavior of cracks in the two materials is analyzed and found to mainly propagate in meandering paths in the CRN region and to be influenced by crystallites in distinct ways that toughen the material. Furthermore, the present results reveal the relation between structure and mechanical properties in amorphous monolayers and may provide a universal toughening strategy for 2D materials.

2-dimensional systems↗
DOE OSTI2020

Significantly increased Raman enhancement on defect-rich O-incorporated 1T-MoS 2 nanosheets

Improving charge transfer is the key to the performance of non-noble metal semiconductor-based surface enhanced Raman scattering (SERS) substrates. In this work, the O-incorporated 1T-MoS 2 nanosheets with rich sulfur defects (ID-MoS 2 ) are obtained by simple calcination of 1T-MoS 2 nanosheets in air atmosphere. Using rhodamine 6G (R6G) as typical probe molecules, ID-MoS 2 nanosheets show ultrahigh Raman enhancement effects with an enhancement factor of 1.24 × 10 7 due to sulfur defects and O incorporation in the 1T-MoS 2 lattice. First-principle density functional theory calculations suggest that the existence of sulfur defects and O incorporation significantly increase the Fermi energy level (E f ) and electronic density of states of ID-MoS 2 . Furthermore, O incorporation can enhance the interactions between the substrate and the adsorbed molecules through electrostatic and hydrogen bonding. All these improve the charge transfer resonance and result in the remarkable SERS activity of ID-MoS 2 nanosheets. This is the first study on the increasing SERS performance of semiconductor substrates by simultaneously employing defect and dopant incorporation. This study provides an approach to optimize the performance of semiconductor-based SERS substrates.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Rotational excitations in rare-earth nuclei: A comparative study within three cranking models with different mean fields and treatments of pairing correlations

High-spin rotational bands in rare-earth Er (Z = 68), Tm (Z = 69), and Yb (Z = 70) isotopes are investigated by three different nuclear models. These are (i) the cranked relativistic Hartree-Bogoliubov (CRHB) approach with approximate particle number projection by means of the Lipkin-Nogami (LN) method, (ii) the cranking covariant density functional theory (CDFT) with pairing correlations treated by a shell-model-like approach (SLAP) or the so-called particle-number conserving (PNC) method, and (iii) cranked shell model (CSM) based on the Nilsson potential with pairing correlations treated by the PNC method. A detailed comparison between these three models in the description of the ground-state rotational bands of even-even Er and Yb isotopes is performed. The similarities and differences between these models in the description of the moments of inertia, the features of band crossings, and equilibrium deformations and pairing energies of the even-even nuclei under study are discussed. These quantities are considered as a function of rotational frequency and proton and neutron numbers. The changes in the properties of the first band crossings with increasing neutron number in this mass region are investigated. On average, a comparable accuracy of the description of available experimental data is achieved in these models. However, the differences between model predictions become larger above the first band crossings. Because of the time-consuming nature of numerical calculations in the CDFT-based models, a systematic study of the rotational properties of both ground-state and excited-state bands in odd-mass Tm nuclei is carried out only by the PNC-SCM. With few exceptions, the rotational properties of experimental 1-quasiparticle and 3-quasiparticle bands in 165,167,169,171 Tm are reproduced reasonably well. As a result, the appearance of backbendings or upbendings in these nuclei is well understood from the analysis of the variations of the occupation probabilities of the single-particle states and their contributions to total angular momentum alignment with rotational frequency.

150 ≤ A ≤ 189↗
DOE OSTI2024

Macroscopic Traffic Modeling Using Probe Vehicle Data: A Machine Learning Approach

Abstract The macroscopic fundamental diagram (MFD) captures an orderly relationship among traffic flow, density, and speed at the network level. It is a simple yet powerful tool for modeling traffic dynamics in large urban networks with broad application in traffic control and management. However, empirically derived MFDs in urban regions require high-resolution traffic data from the network. Having the network flow and vehicular density estimated at the (granular) census tract level using vehicle probe data, we apply machine learning methods to predict the MFDs across U.S. urban areas and capture the impacts of location-specific input features on the network flow–density relationships at a large scale. The results show that, among the four tested machine learning approaches (Random Forest, XGBoost, Support Vector Machine, and Neural Network), XGBoost delivers the best performance in predicting network traffic flow based on vehicular density and location attributes. Using interaction Shapley Additive explanation (SHAP) values and partial correlation analysis, we examine the factors influencing MFD shapes across different locations. Our empirical findings reveal that across U.S. urban areas, network topology, transportation infrastructure, and land use are primary factors shaping MFD curves, while demand and trip-related factors play a lesser role. Specifically, higher ranking roads, centrality, and development levels correlate positively with network capacity and critical density, whereas negative associations are observed for network connectivity, mixed-use development, and road roughness levels.

Jin, Ling↗
DOE OSTI2022

Understanding the role of segmentation on process-structure–property predictions made via machine learning

Here, the present study investigated the effect of porosity surface determination methods on performance of machine learning models used to predict the tensile properties of AlSi10Mg processed by laser powder bed fusion from micro-computed tomography data. Machine learning models applied in this work include support vector machines, neural networks, decision trees, and Bayesian classifiers. The effects of isosurface thresholding and local gradient approaches for porosity segmentation, as well as image filtering schemes, on model precision were evaluated for samples produced under differing levels of global energy density.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Effects of Zr addition on lattice strains and electronic structures of NbTaTiV high-entropy alloy

We report the room-temperature (RT) deformation behavior for two single-phase body-centered-cubic (BCC) refractory high-entropy alloys (RHEAs), NbTaTiV and NbTaTiVZr, has been comprehensively investigated via in-situ neutron-diffraction experiments. Our work shows that the addition of Zr leads to the transition of mechanical response from ductile to brittle behavior. The results of lattice-strain evolutions obtained from in-situ neutron diffraction for the ductile NbTaTiV RHEA exhibit atypical plastic-deformation behavior, i.e., the reduced plasticanisotropic deformation, leading to an even distribution of the applied stress amongst the grains with different orientations rather than forming stress concentrations in {200}-oriented grains during plastic-deformation. Density functional theory (DFT) analysis shows that NbTaTiVZr has a lower electron density at the Fermi level, larger lattice distortion, and stronger charge transfer, as compared to NbTaTiV, suggesting higher strength and lower ductility in NbTaTiVZr, which are consistent with the current experimental results.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Accelerated test protocols to predict service life and durability of solid oxide fuel cells

Reliable accelerated test protocols are needed for solid oxide fuel cell research to facilitate rapid learning on key durability issues, identify potential modes of failure expeditiously, and eventually predict the calendar lifetime of an electrochemical cell. In this work, solid oxide fuel cells operated at a constant current density were compared to cells undergoing accelerated measurements, which are composed of intermittent current injection to the cell. A general accelerated test profile was developed by cycling a solid oxide fuel cell from open circuit to a predetermined operating current density that is the same as the current density during a steady-state operation, to accelerate the local redox environment. The following parameters were studied: current density, operation temperature, moist level, sintering temperature, cycling current, cycling frequency, and operation time. Up to 1,320,000 cycles were generated in this work. The cell degradation was accelerated by nearly 10 times, suggesting the feasibility of using this protocol for acceleration test to predict life performance and durability of solid oxide fuel cells.

08 HYDROGEN↗
DOE OSTI2023

Quantum Vibronic Effects on the Electronic Properties of Molecular Crystals

We present a study of molecular crystals, focused on the effect of nuclear quantum motion and anharmonicity on their electronic properties. We consider a system composed of relatively rigid molecules, a diamondoid crystal, and one composed of floppier molecules, NAI-DMAC, a thermally activated delayed fluorescence compound. We compute fundamental electronic gaps at the density functional theory (DFT) level of theory, with the Perdew–Burke–Erzenhof (PBE) and strongly constrained and approximately normed (SCAN) functionals, by coupling first-principles molecular dynamics with a nuclear quantum thermostat. We find a sizable zero-point renormalization (ZPR) of the band gaps, which is much larger in the case of diamondoids (0.6 eV) than for NAI-DMAC (0.22 eV). We show that the frozen phonon (FP) approximation, which neglects intermolecular anharmonic effects, leads to a large error (~50%) in the calculation of the band gap ZPR. Instead, when using a stochastic method, we obtain results in good agreement with those of our quantum simulations for the diamondoid crystal. However, the agreement is worse for NAI-DMAC where intramolecular anharmonicities contribute to the ZPR. Our results highlight the importance of accurately including nuclear and anharmonic quantum effects to predict the electronic properties of molecular crystals.

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

Efficient 3D Molecular Design with an E(3) Invariant Transformer VAE

This work introduces a three-dimensional (3D) invariant graph-to-string transformer variational autoencoders (VAE) (Vagrant) for generating molecules with accurate density functional theory (DFT)-level properties. Vagrant learns to model the joint probability distribution of a 3D molecular structure and its properties by encoding molecular structures into a 3D-aware latent space. Directed navigation through this latent space implicitly optimizes the 3D structure of a molecule, and the latent embedding can be used to condition a generative transformer to predict the candidate structure as a one-dimensional (1D) sequence. Additionally, we introduce two novel sampling methods that exploit the latent characteristics of a VAE to improve performance. Here, we show that our method outperforms comparable 3D autoregressive and diffusion methods for predicting quantum chemical property values of novel molecules in terms of both sample quality and computational efficiency.

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

Elucidation of Ultrafast Decay, Vibrational Beating, and Slow Decay Processes of Excited Azulene

Azulene’s nonradiative decay dynamics and kinetics from its singlet excited-states were studied using mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) combined with trajectory surface-hopping nonadiabatic molecular dynamics (NAMD) and Fermi’s golden rule (FGR) rate theory. The NAMD dynamics reproduce experimental observations that the S 1 → S 0 decay is accelerated and exhibits a crossover from mono- to biexponential kinetics with increasing excess vibrational energy. Minimum-energy-path analyses reveal a continuous S 1 /S 0 conical-intersection seam slightly above the S 1 minimum, providing readily accessible funnels. Time-resolved normal-mode projections reveal selective energy funneling into C–C stretching modes at 1272, 1528, and 1692 cm –1 . Constructive combinations of the latter two modes appear to promote rapid access to the conical-intersection seam, whereas their beating at ∼ 160 fs imposes a natural limit on the decrease of the decay time upon a further increase of excess energies, suggesting that their interference delays the decay. Here, the FGR rate calculation data for the S 1 → S 0 transition reaffirm that the nonradiative decay proceeds mainly near or through the conical intersection, rather than via simple nonadiabatic derivative coupling around the minimum of S 1 . On the other hand, the FGR rates indicate that S 2 decays predominantly to S 0 , as the S 2 –S 1 vibronic coupling is exceptionally weak, which serves as a primary cause for azulene’s characteristic anti-Kasha emission.

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

Twists and Puckers: Tuning Crystal Chemistry in the La(Au x Ge 1– x ) 2 Compositional Series

The physical properties of solid-state materials are closely tied to their crystal structure, yet our understanding of how competing structural arrangements energetically compare is limited. In this work, we explore how small differences in composition affect structure in the La(Au x Ge 1-x ) 2 series of compounds, comprising four unique structure types between LaGe 2 and LaAu 2 . This family includes the previously unknown AlB2-type compound with the stoichiometry La(Au 0.375 Ge 0.625 ) 2 , as well as La(Au 0.25 Ge 0.75 ) 2 , an intergrowth of the AlB 2 and ThSi 2 structure types. We then study the chemical forces driving the structure changes, including using phonon band structure calculations and DFT Chemical Pressure to evaluate atomic size effects. These calculations show that the parent AlB 2 structure type is disfavored in Au-rich compounds due to soft atomic motions along the c axis. The instability of AlB 2 -type LaAuGe is confirmed by the presence of imaginary modes in the phonon band structure that correspond to a ‘puckering’ of the hexagonal AlB 2 -type lattice, resulting in the experimentally observed LiGaGe structure type. The impact of size effects is less clear for Au-poor compositions; instead, ‘twisting’ the AlB 2 structure type to form the ThSi 2 type opens a pseudogap at the Fermi level in the electronic density of states. Here, this investigation demonstrates how crystal structure in solid-state materials can be compositionally tuned based on balancing size and electronics when multiple structure types are in close thermodynamic competition.

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

What Controls Crystal Diversity and Microphysical Variability in Cirrus Clouds?

Abstract Variability of ice microphysical properties like crystal size and density in cirrus clouds is important for climate through its impact on radiative forcing, but challenging to represent in models. For the first time, recent laboratory experiments of particle growth (tied to crystal morphology via deposition density) are combined with a state‐of‐the‐art Lagrangian particle‐based microphysics model in large‐eddy simulations to examine sources of microphysical variability in cirrus. Simulated particle size distributions compare well against balloon‐borne observations. Overall, microphysical variability is dominated by variability in the particles' thermodynamic histories. However, diversity in crystal morphology notably increases spatial variability of mean particle size and density, especially at mid‐levels in the cloud. Little correlation between instantaneous crystal properties and supersaturation occurs even though the modeled particle morphology is directly tied to supersaturation based on laboratory measurements. Thus, the individual thermodynamic paths of each particle, not the instantaneous conditions, control the evolution of particle properties.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Mechanistic insight of KBiQ 2 (Q = S, Se) using panoramic synthesis towards synthesis-by-design

Solid-state synthesis has historically focused on reactants and end products; however, knowledge of reaction pathways, intermediate phases and their formation may provide mechanistic insight of solid-state reactions. With an increased understanding of reaction progressions, design principles can be deduced, affording more predictive power in materials synthesis. In pursuit of this goal, in situ powder X-ray diffraction is employed to observe crystalline phase evolution over the course of the reaction, thereby constructing a “panoramic” view of the reaction from beginning to end. We conducted in situ diffraction studies in the K–Bi–Q (Q = S, Se) system to understand the formation of phases occurring in this system in the course of their reactions. Powder mixtures of K 2 Q to Bi 2 Q 3 in 1:1 and 1.5:1 ratios were heated to 800 °C or 650 °C, while simultaneously collecting diffraction data. Three new phases, K 3 BiS 3 , β-KBiS 2 , and β-KBiSe 2 , were discovered. Panoramic synthesis showed that K 3 BiQ 3 serves an important mechanistic role as a structural intermediate in both chalcogen systems (Q = S, Se) in the path to form the KBiQ 2 structure. Thermal analysis and calculations at the density functional theory (DFT) level show that the cation-ordered β-KBiQ 2 polymorphs are the thermodynamically stable phase in this compositional space, while Pair Distribution Function (PDF) analysis shows that all α-KBiQ 2 structures have local disorder due to stereochemically active lone pair expression of the bismuth atoms. The formation of the β-KBiQ 2 structures, both of which crystallize in the α-NaFeO 2 structure type, show a boundary where the structure can be disordered or ordered with regards to the alkali metal and bismuth. A cation radius tolerance for six-coordinate cation site sharing of ($\frac{r^+}{r^{3+}}$) ~ 1.3 is proposed. The mechanistic insight the panoramic synthesis technique provides in the K–Bi–Q system is progress towards the overarching goal of synthesis-by-design.

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

Prediction of the structures and heats of formation of MO 2 , MO 3 , and M 2 O 5 for M = V, Nb, Ta, Pa

Structures for the mono-, di-, and tri-bridge isomers of M 2 O 5 as well as those for the MO 2 and MO 3 fragments for M = V, Nb, Ta, and Pa were optimized at the density functional theory (DFT) level. Single point CCSD(T) calculations extrapolated to the complete basis set (CBS) limit at the DFT geometries were used to predict the energetics. Here, the lowest energy dimer isomer was the di-bridge for M = V and Nb and the tri-bridge for M = Ta and Pa. The di-bridge isomers were predicted to be composed of MO 2 + and MO 3 - fragments, whereas the mono- and tri-bridge are two MO 2 + fragments linked by an O 2 - . The heats of formation of M 2 O 5 dimers, as well as MO 2 and MO 3 neutral and ionic species were predicted using the Feller–Peterson–Dixon (FPD) approach. The heats of formation of the MF 5 species were calculated to provide additional benchmarks. Dimerization energies to form the M 2 O 5 dimers are predicted to become more negative going down group 5 and range from -29 to -45 kcal mol -1 . The ionization energies (IEs) for VO 2 and TaO 2 are essentially the same at 8.75 eV whereas the IEs for NbO 2 and PaO 2 are 8.10 and 6.25 eV, respectively. The predicted adiabatic electron affinities (AEAs) range from 3.75 eV to 4.45 eV for the MO 3 species and vertical detachment energies from 4.21 to 4.59 eV for MO 3 - . The calculated M[double bond, length as m-dash]O bond dissociation energies increase from 143 kcal mol -1 for M = V to ~170 kcal mol -1 for M = Nb and Ta to ~200 kcal mol -1 for M = Pa. The M–O bond dissociation energies are all similar ranging from 97 to 107 kcal mol -1 . Natural bond analysis provided insights into the types of chemical bonds in terms of their ionic character. Pa 2 O 5 is predicted to behave like an actinyl species dominated by the interactions of approximately linear PaO 2 + groups.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗