Search NASA⌕ Search

SEARCH · Search NASA

Results for “potential energy”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

How good are polarizable and flexible models for water: Insights from a many-body perspective

We present a systematic analysis of state-of-the-art polarizable and flexible water models from a many-body perspective, with a specific focus on their ability to represent the Born–Oppenheimer potential energy surface of water from the gas to the liquid phase. Using coupled cluster data in the completed basis set limit as a reference, we examine the accuracy of the polarizable models in reproducing individual many-body contributions to interaction energies and harmonic frequencies of water clusters and compare their performance with that of MB-pol, an explicit many-body model that has been shown to correctly predict the properties of water across the entire phase diagram. Based on these comparisons, we use MB-pol as a reference to analyze the ability of the polarizable models to reproduce the energy landscape of liquid water under ambient conditions. We find that, while correctly reproducing the energetics of minimum-energy structures, the polarizable models examined in this study suffer from inadequate representations of many-body effects for distorted configurations. To investigate the role played by geometry-dependent representations of 1-body charge distributions in reproducing coupled cluster data for both interaction and many-body energies, we introduce a simplified version of MB-pol that adopts fixed atomic charges and demonstrate that the new model retains the same accuracy as the original MB-pol model. Based on the analyses presented in this study, we believe that future developments of both polarizable and explicit many-body models should continue in parallel and would benefit from synergistic efforts aimed at integrating the best aspects of the two theoretical/computational frameworks.

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

Mechanism of $\mathrm{E}$-bridge formation by various $\mathrm{PAH}$ molecules: A theoretical study

Potential energy surfaces (PES) for dimerization of acepyrene-pyrene, aceanthracene-anthracene, and acenaphthalene-naphthalene were computed by the G3(MP2,CC) and ONIOM2(G3(MP2,CC):B3LYP/6-31G(d)) methods. Close agreement of relative energies of corresponding species along reaction paths for different reactions was found suggesting that the PES for E-bridge formation depends weakly on the number of rings in the interacting monomers. Our results indicate that the reaction rate is controlled by the dynamics of the initial collision, which may depend on the monomer size. The ONIOM approach is proposed as an accurate tool for calculations of PES in large PAH systems being more reliable than various DFT methods.

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

Quantum reactive scattering calculations for the cold and ultracold Li + LiNa → Li 2 + Na reaction

A first-principles based quantum dynamics study of the Li + LiNa(ν = 0, j = 0) → Li 2 (ν', j') + Na reaction is reported for collision energies spanning the ultracold (1 nK) to cold (1 K) regimes. A full-dimensional ab initio potential energy surface for the ground electronic state of Li 2 Na is utilized that includes an accurate treatment of the long-range interactions. The Li + LiNa reaction is barrierless and exoergic and exhibits a deep attractive potential well that supports complex formation. Thus, significant reactivity occurs even for collision temperatures approaching absolute zero. The reactive scattering calculations are based on a numerically exact time-independent quantum dynamics methodology in hyperspherical coordinates. Total and rotationally resolved rate coefficients are reported at 56 collision energies and include all contributing partial waves. Several shape resonances are observed in many of the rotationally resolved rate coefficients and a small resonance feature is also reported in the total rate coefficient near 50 mK. Furthermore, of particular interest, the angular distributions or differential cross sections are reported as a function of both the collision energy and scattering angle. Unique quantum fingerprints (bumps, channels, and ripples) are observed in the angular distributions for each product rotational state due to quantum interference and shape resonance contributions. The Li + LiNa reaction is under active experimental investigation so that these intriguing features could be verified experimentally when sufficient product state resolution becomes feasible for collision energies below 1 K.

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

Formation and detection of metastable formic acid in a supersonic expansion: High resolution infrared spectroscopy of the jet-cooled cis -HCOOH conformer

High-resolution direct absorption infrared spectra of metastable cis-formic acid (HCOOH) trapped in a cis-well resonance behind a 15 kcal/mol barrier are reported for the first time, with the energetically unstable conformer produced in a supersonic slit plasma expansion of trans-formic acid/H 2 mixtures. We present a detailed high-resolution rovibrational analysis for cis-formic acid species in the OH stretch (ν 1 ) fundamental, providing first precision vibrational band origin, rotational constants, and term values, which in conjunction with ab initio calculations at the couple-cluster with single, double, and perturbative triple [CCSD(T)]/ANOn (n = 0, 1, 2) level support the experimental assignments and establish critical points on the potential energy surface for internal rotor trans-to-cis isomerization. Relative intensities for a- and b-type transitions observed in the spectra permit the transition dipole moment components to be determined in the body fixed frame and prove to be in good agreement with ab initio CCSD(T) theoretical estimates but in poor agreement with simple bond-dipole predictions. As a result, the observed signal dependence on H 2 in the discharge suggests the presence of a novel H atom radical chemical mechanism for strongly endothermic “up-hill” internal rotor isomerization between trans- and cis-formic acid conformers.

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

Ab initio molecular dynamics on quantum computers

Ab initio molecular dynamics (AIMD) is a valuable technique for studying molecules and materials at finite temperatures where the nuclei evolve on potential energy surfaces obtained from accurate electronic structure calculations. In this work, we present an approach to running AIMD simulations on noisy intermediate-scale quantum (NISQ)-era quantum computers. The electronic energies are calculated on a quantum computer using the variational quantum eigensolver (VQE) method. Algorithms for computation of analytical gradients entirely on a quantum computer require quantum fault-tolerant hardware, which is beyond NISQ-era. Therefore, we compute the energy gradients numerically using finite differences, the Hellmann-Feynman theorem, and a correlated sampling technique. This method only requires additional classical calculations of electron integrals for each degree of freedom without any additional computations on a quantum computer beyond the initial VQE run. As a proof of concept, AIMD simulations are demonstrated for the H-2 molecule on IBM quantum devices. In addition, we demonstrate the validity of the method for larger molecules using full configuration interaction wave functions. As quantum hardware and noise mitigation techniques continue to improve, the method can be utilized for studying larger molecular systems.

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

Phase Diagram of a Deep Potential Water Model

Using the Deep Potential methodology, we construct a model that reproduces accurately the potential energy surface of the SCAN approximation of density functional theory for water, from low temperature and pressure to about 2400 K and 50 GPa, excluding the vapor stability region. The computational efficiency of the model makes it possible to predict its phase diagram using molecular dynamics. Satisfactory overall agreement with experimental results is obtained. Here, the fluid phases, molecular and ionic, and all the stable ice polymorphs, ordered and disordered, are predicted correctly, with the exception of ice III and XV that are stable in experiments, but metastable in the model. The evolution of the atomic dynamics upon heating, as ice VII transforms first into ice VII" and then into an ionic fluid, reveals that molecular dissociation and breaking of the ice rules coexist with strong covalent fluctuations, explaining why only partial ionization was inferred in experiments.

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

Theoretical studies of chemical reactions related to the formation and growth of polycyclic aromatic hydrocarbons (PAH) and molecular properties of their key intermediates (Final Progress Report)

The formation mechanisms of polycyclic aromatic hydrocarbons, (PAHs) – organic molecules carrying fused benzene rings – are of great interest to scientists and engineers due to their importance in combustion chemistry and astrochemistry. On Earth, PAHs are largely produced in incomplete combustion of fossil fuel and are considered as critical precursors to unwanted soot particles leading to combustion inefficiency and causing air pollution along with detrimental health effects. Simple PAH molecules initially formed in the gas phase, are further involved in a build-up process in combustion flames leading to larger PAH, bowl-shaped nanostructures, fullerenes, and solid-phase species including carbonaceous dust, graphene particles, and soot. In deep space, PAH and their derivatives are potential key intermediates and nucleation sites leading eventually to carbonaceous nanoparticles (“interstellar grains”). Therefore, the understanding of the key processes in the synthesis of PAHs along with their precursors and their degradation mechanisms in combustion systems and in interstellar, circumstellar, and planetary atmospheric environments will provide critical insights into how complex aromatic structures, carbonaceous nanoparticles, and fullerenes are formed and destroyed. Achieving this understanding is an important step in the development of the efficient combustion processes and of the ecofriendly devices with reduced environmental pollution as well as technological strategies for the production of hydrogen and solid carbon through thermal or plasma-assisted pyrolysis of natural gas and biomass. Also, the understanding of the key processes of PAH and soot growth will help in our comprehension of chemical evolution in the universe. Detailed information on the mechanisms and reliable rate constants of the key elementary chemical reactions involved in PAH formation and destruction processes and in inception of soot particles is often missing, with the main deficiencies being the absence of temperature- and pressure-dependent rate constants for the broad range of conditions occurring in various terrestrial and interstellar processes and the lack of data on the reaction products and their branching ratios. Complementary to experimental studies, these gaps in knowledge can be filled by using quantum chemical calculations of reaction potential energy surfaces providing us with accurate energies of reaction products, intermediates, and transition states, revealing the reaction mechanism, and giving the molecular properties required to compute rate constants for relevant reaction steps and product branching ratios using the RRKM-Master Equation (ME) method. Molecular dynamics (MD) simulations can be used in cases when a reaction rate cannot be properly described by statistical theories. During the terminal renewal project period we employed these ab initio/RRKM-ME and MD approaches to complete our studies on several key reactions relevant to the formation/growth of PAH and inception of soot particles including (1) the reaction mechanism and kinetics of the resonance stabilized fulvenallenyl radical with propargyl and C 3 H 4 isomers; (2) the reaction mechanism and kinetics for the C + indene and C 2 + styrene reactions producing naphthyl or azulenyl radicals in low-temperature environments; (3) the MD study of non-equilibrium dimerization of acepyrene and coronene and its radical. The information derived from our theoretical calculations contributed to a better fundamental understanding of the reaction mechanisms and provide missing critical kinetic data to improve combustion models of hydrocarbon fuels and astrochemical models of the growth of carbonaceous molecules and particles in cold molecular clouds, circumstellar envelopes, and planetary atmospheres.

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

The Reaction of o-Benzyne with Vinylacetylene: An Unexplored Way to Produce Naphthalene

Herein, the mechanisms and kinetics of the reaction of ortho-benzyne with vinylactylene have been studied by ab initio and density functional CCSD(T)-F12/cc-pVTZ-f12//B3LYP/6-311G(d,p) calculations of the pertinent potential energy surface combined with Rice-Ramsperger-Kassel-Marcus - Master Equation calculations of reaction rate constants at various temperatures and pressures. Under prevailing combustion conditions, the reaction has been shown to predominantly proceed by the biradical acetylenic mechanism initiated by the addition of C4H4 to one of the C atoms of the triple bond in ortho-benzyne by the acetylenic end, with a significant contribution of the concerted addition mechanism. Following the initial reaction steps, an extra six-membered ring is produced and the rearrangement of H atoms in this new ring leads to the formation of naphthalene, which can further dissociate to 1- or 2-naphthyl radicals. The o-C 6 H 4 + C 4 H 4 reaction is highly exothermic, by ~143 kcal/mol to form naphthalene and by 31-32 kcal/mol to produce naphthyl radicals + H, but features relatively high entrance barriers of 9-11 kcal/mol. Although the reaction is rather slow, much slower than the reaction of phenyl radical with vinylacetylene, it forms naphthalene and 1- and 2-naphthyl radicals directly, with their relative yields controlled by the temperature and pressure, and thus represents a viable source of the naphthalene core under conditions where ortho-benzyne and vinylacetylene are available.

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

Formation of phenanthrene via H–assisted isomerization of 2–ethynylbiphenyl produced in the reaction of phenyl with phenylacetylene

Model chemistry G3(MP2,CC)//B3LYP/6–311G(d,p) calculations of the potential energy surface for the reaction of phenyl radical (C 6 H 5 ) with phenylacetylene (C 8 H 6 ) have been carried out and combined with Rice–Ramsperger–Kassel–Marcus/Master Equation calculations of temperature– and pressure–dependent rate constants. The results showed that the reaction can serve as a viable source for the formation of phenanthrene via an indirect route involving a primary reaction of phenyl addition to the ortho carbon in the ring of phenylacetylene and H elimination producing 2–ethynylbiphenyl followed by secondary H–assisted isomerization of 2–ethynylbiphenyl to phenanthrene. In the secondary reaction, the H atom adds to the α carbon of the ethynyl side chain, then a six–member ring closure takes place followed by aromatization via an H loss. The channel of H addition to the side chain of 2–ethynylbiphenyl appears to be much faster than H addition to the ortho carbon in the ethynyl–substituted ring leading back to the initial C 6 H 5 + C 8 H 6 reactants. Rate constants for the primary C 6 H 5 + C 8 H 6 ⇌2–ethynylbiphenyl (p1) + H and secondary p1 + H⇌phenanthrene (p2) + H reactions have been computed in the temperature range of 500–2500 K at pressures of 30 Torr, 1, 10, and 100 atm and fitted to modified Arrhenius expressions. Here, the suggested kinetic scheme and rate constants are proposed as a prototype for the modeling of the growth of polycyclic aromatic hydrocarbons via the phenyl addition–dehydrocyclization (PAC) mechanism involving an addition of a PAH radical to an ethynyl–substituted PAH molecule.

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

Conversion of acenaphthalene to phenalene via methylation: A theoretical study

Ab initio calculations of the C 13 H 10 and C 13 H 9 potential energy surfaces related to the reaction of 1-acenaphthyl and methyl radicals and secondary isomerization of C 13 H 9 primary radical products have been performed at the chemically accurate G3(MP2,CC)//B3LYP/6-311G** level of theory to unravel the mechanism of conversion of acenaphthalene to phenalene or phenalenyl radical + H. The computed energetics and molecular parameters were utilized in Rice-Ramsperger-Kassel-Marcus Master Equation (RRKM-ME) calculations of reaction rate constants and relative product yields. The 1-acenaphthyl + CH 3 reaction is predicted to proceed by a fast radical-radical recombination mechanism and to predominantly produce collisionally stabilized 1-methylacenaphthalene B1 or a C 13 H 9 benzylic radical A1 + H with exothermicities of 106.8 and 25.2 kcal/mol, respectively. The A1 + H channel is preferable at higher temperatures, whereas the stabilization of 1-methylacenaphthalene is favored at higher pressures. The radical A1 can nearly irreversibly interconvert to phenalenyl radical AP via a 22.7 kcal/mol exothermic isomerization process involving formal insertion of the CH 2 group into a C-C bond of the five-member ring leading to the expansion of this ring to a six-member ring. The rate constants for the H addition reaction to phenalenyl radical to form phenalene and the reverse H loss from phenalene were also computed and the strength of the weakest C-H bond in the CH 2 group of phenalene is evaluated as 62.2 kcal/mol. The analysis of the reaction kinetics allowed us to deduce a mechanism for the conversion of 1-acenaphthyl radical to phenalene or phenalenyl radical + H via methylation involving the formation of A1 via a well-skipping channel or stabilization and dissociation of (or H abstraction from) B1 followed by isomerization of A1 to phenalenyl AP, which can add an H atom producing phenalene. As a result, rate constants for the significant elementary reactions are fitted to modified Arrhenius expressions and are proposed for kinetic modeling of the expansion of a five-member ring on a zigzag edge of PAH to a six-member ring by methylation.

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

Metal dissolution from first principles: Potential-dependent kinetics and charge transfer

Metal dissolution and its inverse process are integral to both corrosion and electrodeposition; however, many mechanistic details regarding the dissolution process are challenging to decipher. These include how ion dissolution kinetics and charge transfer are influenced by the competition between metal and solvent interactions under an electrode potential. Here, in this work, we introduce a computational framework based on density functional theory with grand-canonical treatment of electrons to directly predict the potential energy landscape for metal dissolution at a constant potential. Using aluminum as an example, we demonstrate that dissolution kinetics is governed by competing kinetics between two physical processes associated with metal–metal bond breaking and ion-migration within the electrical double layer, respectively. We identify a kinetic transition between regimes dominated by each of these processes and show that this transition depends on the operating electrode potential, among other key factors. It is further found that kinetics and thermodynamics of these processes can be described with a simple, one-parameter Marcus-theory-type model. Beyond offering new understanding of charge transfer during dissolution, our simulation protocol provides a recipe for directly predicting other important quantities in electrochemical reactions from first principles that are difficult to measure, such as the symmetry factor.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Analysis of Energy Savings and CO2 Emission Reduction Contribution for Industrial Facilities in USA

Abstract Energy audits can identify energy consumptions, energy costs of the facility and evolve to develop measures to maximize efficiency, optimize supply energy, and eliminate waste. This paper investigates the potential energy savings at 20 different industrial sectors with 152 assessments for various facilities in Wisconsin, USA. On average, eight energy recommendations were suggested and applied in each facility. This paper provides a detailed guideline for each industry in terms of eight different energy categories: heating, ventilation, and air conditioning (HVAC) systems, heat recovery systems, electrical demand management, and utility bills, motors, compressors, waste management, and productivity enhancement, lighting, besides building envelope. In total, the energy savings were as follows: 98 million kWh in the shape of electricity, 561 billion British thermal units (BTUs) gas savings, 44 million gallons water savings, and 2-million-pound solid waste savings. Based on these savings, a 100-thousand-ton reduction in carbon dioxide emissions was obtained.

Energy & Fuels↗
DOE OSTI2020

Accelerating Finite-temperature Kohn-Sham Density Functional Theory\ with Deep Neural Networks

We present a numerical modeling workflow based on machine learning (ML) which reproduces the the total energies produced by Kohn-Sham density functional theory (DFT) at finite electronic temperature to within chemical accuracy at negligible computational cost. Based on deep neural networks, our workflow yields the local density of states (LDOS) for a given atomic configuration. From the LDOS, spatially-resolved, energy-resolved, and integrated quantities can be calculated, including the DFT total free energy, which serves as the Born-Oppenheimer potential energy surface for the atoms. We demonstrate the efficacy of this approach for both solid and liquid metals and compare results between independent and unified machine-learning models for solid and liquid aluminum. Our machine-learning density functional theory framework opens up the path towards multiscale materials modeling for matter under ambient and extreme conditions at a computational scale and cost that is unattainable with current algorithms.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Accelerating Finite-Temperature Kohn-Sham Density Functional Theory with Deep Neural Networks

We present a numerical modeling workflow based on machine learning (ML) which reproduces the total energies produced by Kohn-Sham density functional theory (DFT) at finite electronic temperature to within chemical accuracy at negligible computational cost. Based on deep neural networks, our workflow yields the local density of states (LDOS) for a given atomic configuration. From the LDOS, spatially-resolved, energy-resolved, and integrated quantities can be calculated, including the DFT total free energy, which serves as the Born-Oppenheimer potential energy surface for the atoms. We demonstrate the efficacy of this approach for both solid and liquid metals and compare results between independent and unified machine-learning models for solid and liquid aluminum. Our machine-learning density functional theory framework opens up the path towards multiscale materials modeling for matter under ambient and extreme conditions at a computational scale and cost that is unattainable with current algorithms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Repurposing Inactive Oil and Gas Wells for Energy Storage: Maximizing the Potential Via Optimal Drivetrain Control: Preprint

In recent years, there has been a growing emphasis on utilizing energy storage to enhance grid resilience against disruptive events. While renewable energy supply continues to expand, gravity-based solutions like pumped hydro remain dominant in the commercial space. However, their geographical limitations constrain availability, scalability, and increase costs for co-locating solar and wind energy. An alternative approach proposes repurposing idle oil and gas wells located closer to existing grid infrastructure, offering a promising and cost-effective solution. This paper addresses the optimization and control of a regenerative drive system coupled to an interior permanent magnet synchronous machine in a 300-meter well with a 100-Newton weight. The study employs a dynamic MATLAB/Simulink model to simulate the operation of the electric drivetrain system during storage and discharge operations. The results demonstrate an initial round-trip efficiency of 85.9% for the electrical system alone and identify crucial factors for maximizing efficiency. The optimized operation and control of the electromechanical drivetrain system hold great potential for minimizing the levelized cost of storage while maximizing efficiency and revenue generation.

ENERGY STORAGE,ENGINEERING↗
DOE OSTI2023

Calorimetric Energies and Chemical Potentials of Metal Atoms in Catalytic Nanoparticles on Oxide and Carbon Supports: Improved Size Dependencies and Adhesion Energies

The chemical potential of metal atoms, μM, in supported metal nanoparticles is an important descriptor related to both the catalytic activity and the stability of the nanoparticles. Here, we derive an expression relating μM to the radius of the particle’s contact area with the support and the adhesion energy at the metal/support interface (E adh ) that assumes the particles have the shape of spherical caps but of arbitrary contact angle with the support (θ c ) and includes an empirical correction for the increase in metal surface energy and adhesion energy with decreasing radius of curvature. We then show that, at any assumed contact angle, we can simultaneously fit previously reported measurements of both calorimetric μ M (from heats of metal vapor adsorption during nanoparticle growth by vapor deposition) versus metal coverage data and the He + low-energy ion scattering (LEIS) intensities for the metal and/or support versus metal coverage (using our recently developed spherical cap model for quantitative LEIS intensities), to determine the particle size versus coverage and E adh . Only one choice of contact angle gives a pair of values for contact angle and E adh , which is consistent with the Young–Dupré equation for the equilibrium shape of a spherical particle. At this equilibrium shape, we then applied this spherical cap model (SCM) to reanalyze microcalorimetric metal chemical potentials and LEIS signals versus coverage data for nine metal/support combinations that were previously analyzed by assuming that the particles had the shape of hemispherical caps, i.e., with a contact angle of 90°. We show that this revised approach gives close agreement with the calorimetric and LEIS data; the best-fit contact angles vary from 64 to 84°, correcting the earlier assumption of 90°. Furthermore, these results provide significant accuracy improvements in particle size versus coverage, metal chemical potential versus size and coverage, metal/support adhesion energies and contact angles for Cu, Ag and Au on CeO 2 (111), Ni on MgO(100), Ag on Fe 3 O 4 (111) and TiO 2 (100), and Ag, Ni and Pd on Ni-supported graphene. This revised approach is much more broadly applicable than the earlier hemispherical cap model (HCM).

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2022

Preliminary Assessment of Potential for Wind Energy Technology on the Turtle Mountain Band of Chippewa Reservation

Wind energy can provide renewable, sustainable electricity to rural Native homes and power schools and businesses. It can even provide tribes with a source of income and economic development. The purpose of this research is to determine the potential for deploying community and utility-scale wind renewable technologies on Turtle Mountain Band of Chippewa tribal lands. Ideal areas for wind technology development were investigated, based on wind resources, terrain, land usage, and other factors. This was done using tools like the National Renewable Energy Laboratory Wind Prospector, in addition to consulting tribal members and experts in the field. The result was a preliminary assessment of wind energy potential on Turtle Mountain lands, which can be used to justify further investigation and investment into determining the feasibility of future wind technology projects.

17 WIND ENERGY↗
DOE OSTI2025

Basin-Size Mapping: Prediction of Metastable Polymorph Synthesizability Across TaC–TaN Alloys

The sizes of the basins of attraction on the potential energy surface are helpful indicators in determining the experimental synthesizability of metastable phases. In principle, these basins can be controlled with changes in thermodynamic conditions such as composition, pressure, and surface energy. Herein, we use random structure sampling to computationally study how alloying smoothly perturbs basin of attraction sizes. The TaC 1-x N x pseudobinary is an ideal test system given the structural and polymorphic contrast of its parent compounds and their technological relevance as epitaxial substrates for Al 1-x Ga x N. While we find limited thermodynamic stability across all computationally observed phases, random structure sampling shows a significant composition region where the rocksalt basin dominates. As such, we predict the potential for the nonequilibrium synthesis of metastable rocksalt TaC 1-x N x alloys as substrates for Al 1-x Ga x N. At higher nitrogen concentrations, other low-energy metastable polymorphs emerge that continue to retain the hexagonal close packing suitable for III-N growth. Confidence in these trends was established through uncertainty quantification of the basin sizes and energy distributions; such analysis utilized the Beta and Dirichlet distributions. In conclusion, we also find (a) polymorph basin sizes can be rationalized in terms of energetic preferences for different coordination environments; and (b) basin sizes universally shrink with increasing nitrogen content, making the system more prone to amorphous growth.

36 MATERIALS SCIENCE↗