Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “EV”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Exploring impacts of electricity tariff on charging infrastructure planning: An activity-based approach

In the past decade, electric vehicles (EVs) have gained popularity for their efficiency and environmental benefits. Advances in battery technology and charging equipment have yielded long-range EVs and fast-charging. However, many major cities lack adequate charging infrastructure for daily EV use. This study addresses this gap by integrating activity-based modeling, charging behavior simulation, and charging infrastructure optimization. The research utilizes the POLARIS agent-based transportation model to accurately capture user activities, trip patterns, and traffic flows. Additionally, the study investigates the impact of fixed and spatiotemporal electricity rate distributions on optimal charging infrastructure deployment. The framework is applied to the Chicago regional area network and analyzed under various EV ownership scenarios. Further, the results reveal significant impacts of the charging pricing strategy on user decision-making and charging demand distribution. There is also a need for consistent pricing policies in charging infrastructure planning and operational phases to avoid drops in service quality.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

Investigation of Off-Stoichiometry in Ternary Nitrides by EDS and HRTEM

We demonstrate an inverted metamorphic multijunction (IMM) photovoltaic cell comprising lattice-mismatched 1.2 eV AlGaInAs and 1.0 eV GaInAs junctions optimized for high-temperature thermophotovoltaic (TPV) applications. This device differs from traditional IMM solar cells because the mismatched junctions are grown at a single lattice constant. This architecture enables removal of the compositionally graded buffer which otherwise filters light from the junctions below and absorbs sub-bandgap light via free-carrier absorption. Sub-bandgap absorption dramatically reduces the efficiency of TPV systems using high reflectivity cells to enable bandedge spectrum filtering. Three components required development to enable this device: 1) a lattice-mismatched 1.2 eV AlGaInAs junction, 2) a metamorphic contact layer grown after the graded buffer, and 3) a transparent tunnel junction that sits in front of the 1.0 eV GaInAs junction. Growth conditions that minimize oxygen defect incorporation maximize AlGaInAs cell quality, enabling a 0.41 V bandgap open circuit voltage offset at 22 mA/cm2 under AM1.5D. A mismatched GaInAs:Se layer is developed as a low resistance contact. Lastly we develop a GaAsSb:C/GaInP:Se tunnel junction suitable for high-power densities with more transparency than the GaAsSb:C/GaInP:Se structure used in past IMM cells. We characterize the tandem device under a high-intensity spectrum that approximates the emission from a 2150 °C blackbody radiator and deduce a projected ideal TPV efficiency of 39.9% at ~30% of the blackbody irradiance, and 36% ideal TPV efficiency under the full 118 W/cm2 irradiance. Improvements to the back-surface reflectivity and series resistance are expected to increase the ideal TPV efficiency well above 40%.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Structure−Property Relationships in a Series of Hybrid Halopyridinium Ruthenium(IV) Halides

Reported is the synthesis and characterization of a family of hybrid ruthenium halides consisting of haloruthenium octahedra (X = Cl, Br) charge balanced with halopyridinium (XPy; X= H, Cl, Br, I) organic cations which assemble via noncovalent interactions between ion pairs. Diffuse reflectance spectroscopy showed that compounds containing RuBr 6 2- octahedra displayed a lower bandgap (1.05 eV< x < 1.08 eV) compared to compounds with RuCl 6 2- octahedra (1.22 eV < x < 1.43 eV). Additionally, computational density functional theory (DFT) based natural bonding orbital (NBO) analysis and density of state (DOS) methods were used to characterize second-sphere non-covalent interaction strengths and elucidate molecular orbital perturbations to elucidate their influence on Ru-X orbital constructs and in turn rationalize band gap trends. Analyses showed ligand to metal charge transfer is responsible for the small bandgap energies and are unperturbed by second-sphere interactions. Here, this report serves as a platform for probing the relationship between the structural and photophysical properties within the haloruthenium family of low dimensional transition metal halide perovskites.

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

Local Structural Disorder in Metavanadates MV 2 O 6 (M = Zn and Cu) Synthesized by the Deep Eutectic Solvent Route: Photoactive Oxides with Oxygen Vacancies

Metavanadates MV 2 O 6–δ (M = Zn and Cu) are synthesized by using a deep eutectic solvent (DES), a mixture of hydrogen bond donor and acceptor, as a reaction medium. Dissolution of stable binary metal oxides in a DES followed by a heat treatment yields phase-pure vanadates. According to in situ powder X-ray diffraction, ternary phases (α-Zn 2 V 2 O 7 and β-Cu 2–x V 2 O 7 , x ~ 0.27) are intermediates in the reaction pathway taken. Identifying a polymorphic phase transformation temperature for CuV 2 O 6 as well as the narrow temperature range between formation and decomposition for both metavanadates allows for tackling the challenge of the synthesis of these materials. The oxygen vacancy introduced by the DES route is accompanied by the formation of reduced V 4+ and Cu + in the oxide matrix, based on X-ray photoelectron spectroscopy. These oxygen vacancies modify the vibrational modes in the corresponding Raman spectra and are also responsible for broad optical absorptions in the 1.8–1.1 eV range. The optical band gaps of the materials are found at 1.8 eV (CuV 2 O 6 ) and 2.2 eV (ZnV 2 O 6 ), approximately 0.1–0.3 eV below the values reported in the literature. The reduced band gaps and sub-band gap photon absorption are key features of the oxygen-deficient metavanadates. Surface photovoltage spectroscopy reveals that all the synthesized vanadates are n-type materials with electrons as the majority charge carriers, and photoelectrochemical measurements confirm photoanodic currents for methanol oxidation. Furthermore, these results provide insight into the synthesis and structure–property relationship of the metavanadates for their potential use as photoanodes.

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

Phase Stability, Band Gap Tuning, and Rashba Splitting in Selenium-Alloyed Bournonite: CuPbSb(S 1– x Se x ) 3

Recently, bournonite (CuPbSbS 3 ) has been identified as a potential ferroelectric photovoltaic (PV) material with a 1.3 eV band gap, which falls in an appropriate range for single junction PV devices. Progress in applying bournonite has yielded several studies reporting successful thin-film processing, the most recent of which demonstrated a 2.65% power conversion efficiency for a bournonite-based PV device. In an effort to explore bournonite band gap engineering for PV and other prospective applications (e.g., thermoelectric, spintronic), we here consider the solid-state synthesis of selenium-alloyed bournonite, CuPbSb(S 1-x Se x ) 3 , across the full range of x (0.0 ≤ x ≤ 1.0), and report phase purity for 0.0 ≤ x ≤ 0.5. We characterize the crystal structure and band gap of the samples using X-ray diffraction and diffuse reflectance spectroscopy, determining a band gap decrease for single-phase samples from 1.25 eV at x = 0.0 to 1.06 eV at x = 0.5. Formation energies determined using dispersion corrected hybrid density functional theory (DFT) support experimental findings and are consistent with the stability (instability) of the x = 0.5 (x = 1.0) structure relative to starting materials. The computed band structures show a decreasing trend in the band gap with increasing x, again consistent with experiment, with the states at the conduction (valence) band minima (maxima) being principally derived from Pb (S/Se) states. Spin texture analysis and detailed comparison of spin splitting parameters show that Se-alloying modifies the band gap while maintaining Rashba spin splitting character within the electronic structure. Rashba splitting is most pronounced for the conduction band minimum along the Γ-X path in the Brillouin zone. Energy separation between spin states is maintained at ΔE~0.2 eV for the various Se contents, with σz as the spin direction. Furthermore, a smaller spin splitting occurs along Γ-Z (decreasing value with increasing x), with σ x as the spin direction.

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

Flux-Assisted Boron Chalcogen Mixture (BCM) Method for Synthesizing Mixed Chalcogenide Semiconductors (AkRE 2 Si 2 Se x S 8– x and CaRE 2 Si 2 Se 8 ) ( Ak = Ca and Sr; RE = La, Ce, Pr, Nd, and Sm): Investigation of Their Magnetic and Optical Properties

We report a detailed structural analysis of a series of ten quaternary rare-earth-containing seleno-thiosilicates AkRE 2 Si 2 Se x S 8-x and selenosilicates, CaRE 2 Si 2 Se 8 (Ak = Ca and Sr; RE = La, Ce, Pr, Nd, and Sm). Single crystals were obtained by using the flux-assisted boron chalcogen mixture (BCM) method and single crystal X-ray diffraction was used to determine their structures. All members of the AkRE 2 Si 2 Se x S 8-x and CaRE 2 Si 2 Se 8 series crystallize in the space group R$\bar{3}$c (space group number 167) of the trigonal crystal system. The single-crystal X-ray diffraction analysis revealed a strong preference for Se/S atoms to occupy one vs. the other of the two available sites. Polycrystalline samples were used for magnetic susceptibility and UV–visible diffuse reflectance measurements. Magnetic measurements show that CaCe₂Si₂Se₁.₇₃S₆.₂₇ and CaNd₂Si₂Se₂.₅S₅.₅ are paramagnetic with negative Weiss constants (θ = –60.1 and –26.2). Diffuse reflectance analysis gives optical band gaps of 2.7(1) eV (CaLa₂Si₂Se₂.₃₈S₅.₆₂), 2.2(1) eV (CaCe₂Si₂Se₁.₇₃S₆.₂₇), 2.5(1) eV (CaNd₂Si₂Se₂.₅S₅.₅), and 2.0(1) eV (CaCe₂Si₂Se₈), consistent with density functional theory calculations. By partially or fully replacing S sites with Se it was possible to achieve band gap tuning. Photoluminescence behavior was also investigated for CaCe 2 Si 2 Se 1.73 S 6.27 via irradiation with 375 nm ultraviolet light.

BCM method↗
DOE OSTI2021

Quantitative Account of the Bonding Properties of a Rubredoxin Model Complex [Fe(SCH3)4]q, q = -2, -1, +2, +3

Iron-sulfur clusters play important roles in biology as parts of electron transfer chains and catalytic cofactors. Here, we report a detailed computational analysis of a structural model of the simplest natural iron-sulfur cluster of rubredoxin and its cationic counterparts. Specifically, we report results for the ground and low-lying electronically excited states of the complex [Fe(SCH3)4]2-/1-/2+/3+, using Multi-Reference (CASSCF, MRCISD), and Coupled Cluster [CCSD(T)] methodology in order to provide accurate adiabatic reduction energies, dissociation energies and insights into the bonding analysis. The nature of the Fe-S chemical bond and the magnitude of the ionization potentials in the anionic and cationic [Fe(SCH3)4] complexes offer a physical rationale for the relative stabilization, structure and speciation of these complexes. Anionic and cationic complexes present different types of chemical bonds: prevalently ionic in [Fe(SCH3)4]2-/1- complexes and covalent in [Fe(SCH3)4]2+/3+ complexes. The ionic bonds result in an energy gain for the transition [Fe(SCH3)4]2-®[Fe(SCH3)4]- (i.e., FeII®FeIII) of 1.5 eV, while the covalent bonds result in an energy loss for the transition [Fe(SCH3)4]2+®[Fe(SCH3)4]3+ of 16.6 eV, almost half of the IP of Fe2+. The ionic vs covalent bond character influences the Fe-S bond strength and length, i.e., ionic Fe-S bonds are longer than covalent ones by about 0.2 Å (for FeII) and 0.04 Å (for FeII). Finally, the average Fe-S heterolytic bond strength is 6.7 eV (FeII) and 14.6 (FeIII) eV at the RCCSD(T) level of theory.

Tzeli, Demeter↗
DOE OSTI2026

Low-Lying Excited States of Linear All- Trans Polyenes: Insights from Analytic Gradient and Nonadiabatic Coupling Calculations Based on Multireference Configuration Interaction

Polyenes serve as a rigorous test for theoretical models and electronic structure methods, playing a key role in advancing computational and theoretical chemistry. Here, we present a high-level theoretical investigation of linear, all-trans polyenes using energy gradients and nonadiabatic coupling vectors based on an MR-CISD wave function to describe electronic transitions involving the ground state (1 1 A g – ) and three low-lying excited states (2 1 A g – , 1 1 B u + , and 2 1 B u – ) of hexatriene, octatetraene, and decapentaene. This approach enables accurate evaluation of both adiabatic and vertical excitation and emission energies, yielding results in excellent agreement with experiment, as well as locating minima on the crossing seam between adiabatic states. Our results show that vertical excitation energies to the 1 1 B u + state are blue-shifted by 0.2–0.3 eV relative to the experimental absorption maximum, whereas the vertical emission energy from the 2 1 A g – state is red-shifted by ∼0.2 eV relative to the experimental emission maximum. Upon relaxation from the Franck–Condon geometry, the 2 1 A g – state stabilizes by around 1 eV, compared to 0.2–0.3 eV for the 1 1 B u + state. An analysis of the S 1 /S 0 crossing seam in hexatriene shows that its minimum involves asymmetric backbone deformations and provides an efficient channel for ultrafast internal conversion to the ground state, consistent with the absence of detectable fluorescence in this molecule. These results demonstrate the power of analytic gradients and nonadiabatic coupling vectors based on an MR-CISD wave function for accurately characterizing the electronic structure and photophysics of polyenes.

Excited states↗
DOE OSTI2022

Interaction of Th with H 0/–/+ : Combined Experimental and Theoretical Thermodynamic Properties

In this work, high-level electronic structure calculations of the lowlying energy electronic states for ThH, ThH – , and ThH + are reported and compared to experimental measurements. The inclusion of spin–orbit coupling is critical to predict the ground-state ordering as inclusion of spin–orbit switches the coupled-cluster CCSD(T) ordering of the two lowest energy states for ThH and ThH + . At the multireference spin–orbit SO-CASPT2 level, the ground states of ThH, ThH – , and ThH + are predicted to be the 2 Δ 3/2 , 3 Φ 2 , and 3 Δ 1 states, respectively. The adiabatic electron affinity is calculated to be 0.820 eV, and the vertical detachment energy is calculated to be 0.832 eV in comparison to an experimental value of 0.87 ± 0.02 eV. The observed ThH – photoelectron spectrum has many transitions, which approximately correlate with excitations of Th + and/or Th. The adiabatic ionization energy of ThH including spin–orbit corrections is calculated to be 6.181 eV. The natural bond orbital results are consistent with a significant contribution of the Th + H – ionic configuration to the bonding in ThH. The bond dissociation energies for ThH, ThH – , and ThH + using the Feller–Peterson–Dixon approach were calculated to be similar for all three molecules and lie between 259 and 280 kJ/mol.

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

Experimental and Computational Description of the Interaction of H and H – with U

In this work, the results of ab initio correlated molecular orbital theory electronic structure calculations for low-lying electronic states are presented for UH and UH – and compared to photoelectron spectroscopy measurements. The calculations were performed at the CCSD(T)/CBS and multireference CASPT2 including spin–orbit effects by the state interacting approach levels. The ground states of UH and UH – are predicted to be 4 Ι 9/2 and 5 Λ 6 , respectively. The spectroscopic parameters T e , r e , ω e , ω e x e , and Be were obtained, and potential energy curves were calculated for the low energy Ω states of UH. The calculated adiabatic electron affinity is 0.468 eV in excellent agreement with an experimental value of 0.462 ± 0.013 eV. The lowest vertical detachment energy was predicted to be 0.506 eV for the ground state, and the adiabatic ionization energy (IE) is predicted to be 6.116 eV. The bond dissociation energy (BDE) and heat of formation values of UH were obtained using the IE calculated at the Feller–Peterson–Dixon level. For UH, UH – , and UH + , the BDEs were predicted to be 225.5, 197.9, and 235.5 kJ/mol, respectively. The BDE for UH is predicted to be ~20% lower in energy than that for ThH. The analysis of the natural bond orbitals shows a significant U + H – ionic component in the bond of UH.

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

Theoretical and Experimental Study of the Spectroscopy and Thermochemistry of UC +/0/–

A combination of high-level ab initio calculations and anion photoelectron detachment (PD) measurements is reported for the UC, UC – , and UC + molecules. To better compare the theoretical values with the experimental photoelectron spectrum (PES), a value of 1.493 eV for the adiabatic electron affinity (AEA) of UC was calculated at the Feller–Peterson–Dixon (FPD) level. The lowest vertical detachment energy (VDE) is predicted to be 1.500 eV compared to the experimental value of 1.487 ± 0.035 eV. A shoulder to lower energy in the experimental PD spectrum with the 355 nm laser can be assigned to a combination of low-lying excited states of UC – and excited vibrational states. The VDEs calculated for the low-lying excited electronic states of UC at the SO-CASPT2 level are consistent with the observed additional electron binding energies at 1.990, 2.112, 2.316, and 3.760 eV. Potential energy curves for the Ω states and the associated spectroscopic properties are also reported. Compared to UN and UN + , the bond dissociation energy (BDE) of UC (411.3 kJ/mol) is predicted to be considerably lower. The natural bond orbitals (NBO) calculations show that the UC 0/+/– molecules have a bond order of 2.5 with their ground-state configuration arising from changes in the oxidation state of the U atom in terms of the 7s orbital occupation: UC (5f 2 7s 1 ), UC – (5f 2 7s 2 ), and UC + (5f 2 7s 0 ). Furthermore, the behavior of the UN and UC sequence of molecules and anions differs from the corresponding sequences for UO and UF.

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

Properties of NoO and YbO

Bond dissociation energies (BDEs) and spectroscopic parameters for the ground states of YbO and NoO were calculated at the Brueckner doubles (BD(T)) level and for NoO at the coupled cluster CCSD(T) level based on the closed-shell configurations Yb 2+ (4f 14 )O 2 – (2p 6 ) and No 2+ (5f 14 )O 2 – (2p 6 ). For YbO, the BD(T) BDE (3.93 eV) and vibrational frequency (686.3 cm –1 ) are consistent with the experiment; the BD(T) bond distance was 1.8492 Å, about 0.04 Å longer than the experiment. For NoO, the calculated BDEs are 3.07 eV (BD(T)) and 3.08 eV (CCSD(T)) and are expected to be good to ±3 kcal/mol (±0.13 eV). Potential energy curves (PECs) for the spin–orbit (Ω) states associated with the f 13 s 1 configurations of YbO and NoO were calculated at the SA-CASSCF/SO-CASPT2/aQ-DK level. PECs considering both the f 13 s 1 and f 14 configurations were calculated at this level, as the ground states are expected to contain both configurations. Due to a bias in the SA-CASSCF to the f 13 s 1 configurations, the f 14 configurations could not be properly described simultaneously with the f 13 s 1 ones. Finally, the natural population analysis based on the natural bond orbitals shows that YbO and NoO are quite ionic, with 5d orbitals being more important in YbO than the 6d orbitals in NoO.

basis sets↗
DOE OSTI2025

Energetic and Electronic Properties of UX +/0/– for X = Li and Be and Comparison of the Properties of the Uranium Atom Binding to 2nd Row Elements Li–F

The bonding and spectroscopic properties of ULi +/0/– and UBe +/0/– to complete the series for UX +/0/– for X = Li to F were investigated by high-level ab initio SO-CASPT2 and CCSD(T) electronic structure calculations. The low-lying spin–orbit states were obtained at the SA-CASPT2/aQ-PP level; bond dissociation energies (BDEs), ionization energies (IEs), adiabatic electronic affinities (AEAs), and vertical detachment energies (VDEs) were calculated at the Feller-Peterson-Dixon (FPD) level. A dense manifold of low-lying states was predicted for ULi +/0/– and UBe +/0/– . Here, the calculated BDEs for ULi (37.7 kJ/mol) and UBe (8.0 kJ/mol) show that UBe is weakly bound. For redox processes, the BDEs increased for ULi + (109.3 kJ/mol), ULi – (47.4 kJ/mol), UBe + (35.6 kJ/mol), and UBe – (72.3 kJ/mol). The IE(ULi) = 4.650 eV is lower than IE(Li); the IE(UBe) = 5.901 eV is close to the IE(U) and to the IEs of UB, UC, UN, UO, and UF. The AEAs of ULi (0.708 eV) and UBe (0.989 eV) are lower than those for UB, UC, UN, and UO but higher than that for EA(UF). Natural bond orbital (NBO) calculations show that ULi has the 5f 3 6d 1 7s 2 configuration for U and 2s 1 for Li, with a small partial negative charge slightly delocalized on U. UBe arises from the U(5f 3 6d 1 7s 2 ) and Be(2s 2 ) electron configurations with no charge separation. The same calculations were made for WX (X = Li, Be, C–F) to enable detailed comparisons of the properties for UX with WX (X = Li–F). For WX, BDE(WX) is higher than that for UX for X = Li to N and lower than BDE(UX) for X = O and F, mostly due to the higher IE of W than U as ionic character becomes more important going from Li to F.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Experimental and Theoretical Study of the Electronic Structures of Lanthanide Indium Perovskites LnInO 3

Ternary lanthanide indium oxides LnInO 3 (Ln = La, Pr, Nd, Sm) were synthesized by high-temperature solid-state reaction and characterized by Xray powder diffraction. Rietveld refinement of the powder patterns showed the LnInO 3 materials to be orthorhombic perovskites belonging to the space group Pnma, based on almost-regular InO 6 octahedra and highly distorted LnO 12 polyhedra. Experimental structural data were compared with results from density functional theory (DFT) calculations employing a hybrid Hamiltonian. Valence region X-ray photoelectron and K-shell X-ray emission and absorption spectra of the LnInO 3 compounds were simulated with the aid of the DFT calculations. Photoionization of lanthanide 4f orbitals gives rise to a complex final-state multiplet structure in the valence region for the 4f n compounds PrInO 3 , NdInO 3 , and SmInO 3 , and the overall photoemission spectral profiles were shown to be a superposition of final-state 4f n-1 terms onto the cross-section weighted partial densities of states from the other orbitals. The occupied 4f states are stabilized in moving across the series Pr-Nd-Sm. Band gaps were measured using diffuse reflectance spectroscopy. These results demonstrated that the band gap of LaInO 3 is 4.32 eV, in agreement with DFT calculations. This is significantly larger than a band gap of 2.2 eV first proposed in 1967 and based on the idea that In 4d states lie above the top of the O 2p valence band. However, both DFT and X-ray spectroscopy show that In 4d is a shallow core level located well below the bottom of the valence band. Band gaps greater than 4 eV were observed for NdInO 3 and SmInO 3 , but a lower gap of 3.6 eV for PrInO 3 was shown to arise from the occupied Pr 4f states lying above the main O 2p valence band.

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

Nature of Oxygen Adsorption on Defective Carbonaceous Materials

Here, plane-wave density functional theory has been used to study oxygen adsorption on graphene, graphite, and (12,0) zigzag single-walled carbon nanotubes with and without Stone–Wales (SW) and single-vacancy (SV) defects to understand the role of defects on carbonaceous material reactivity. Atomic oxygen adsorption leads to the formation of an epoxide on defect-free graphene and graphite and an ether on the exterior wall of carbon nanotubes and SW-defected materials. O 2 chemisorption is endothermic on defect-free graphene and graphite and slightly exothermic on defect-free nanotubes. O 2 chemisorption energies are predicted to be –1.1 to –1.4 eV on an SW defect and –6.0 to –8.0 eV on an SV defect. An SW defect lowers the energy barriers by 0.90 and 0.50 eV for O 2 chemisorption on graphene and nanotubes, respectively. The formation of a C–O–O–C group is important for O 2 dissociation on defect-free and SW-defected materials. The energy barrier is less than 0.30 eV on an SV defect. The more reactive SW defect toward O adsorption on graphene is mostly due to the strained defective carbon atoms being able to donate more electrons to an O to form an ether. The larger 2s character in the hybrid orbitals in an ether than in an epoxide makes the ether C–O bond stronger. Stronger C–O binding on an SW-defective carbon nanotube than on a defect-free nanotube is in part due to more flexibility of the defect to release the epoxide ring strain to form an ether.

25 ENERGY STORAGE↗
DOE OSTI2022

Structure, Stability, and Electronic Properties of Boron Suboxide: A Density Functional Theory Study

Boron suboxide (B 6 O) is a boron-rich compound derived from the α-rhombohedral boron lattice with extreme hardness and unusual semiconducting properties. In this work, density functional theory (DFT) was used to show that unit cell volume, mechanical strength, band gaps, and thermodynamic stabilities of B 6 O were influenced by the interstitial elements and point defects at the icosahedral sites. While the hexagonal unit cell volume (HUCV) varies with interstitial occupancy, it is the icosahedral defect that weakens the intrinsic bulk modulus of B 6 O. Using the hybrid HSE functional, we confirmed that the perfect B 6 O bulk is a p-type semiconductor with a direct band gap of 2.8 eV. Furthermore, by screening α-boron compounds systematically, we found that a simple octet rule may offer a consistent explanation for the variations in the computed electronic structures. The formation free energies calculated over a wide range of temperatures (0–2500 K) and pressures (0–80 GPa) predict that formations of interstitial defects become favorable only at higher temperatures (ca. 1800 K) in bulk B 6 O lattices. The nudged elastic band (NEB) method was employed to identify the minimum energy pathways for the diffusions of dislocated B and O atoms. The diffusion of icosahedral B atoms has an energy barrier of 0.16 eV. More complex B diffusion paths involving the reorganization of icosahedral boron atoms incur higher barriers (>1 eV). In contrast, the diffusion of interstitial O atoms is facile with a barrier of 0.4 eV. Lastly, successive O insertions into the α-B lattice were performed using DFT to generate a basic understanding of the oxidation process. Furthermore, these calculations provide fundamental atomistic insights into the growth of B 6 O crystals and control of their point defects.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Resolving the Coverage Dependence of Surface Reaction Kinetics with Machine Learning and Automated Quantum Chemistry Workflows

Microkinetic models for catalytic systems require estimation of many thermodynamic and kinetic parameters that can be calculated for isolated species and transition states using ab initio methods. However, the presence of nearby coadsorbates on the surface can dramatically alter these thermodynamic and kinetic parameters causing them to be dependent on species coverage fractions. As there are combinatorially many coadsorbed configurations on the surface, computing the coverage dependence of these parameters is far less straightforward. We present a framework for generating and applying machine learning models to predict coverage-dependent parameters for microkinetic models. Our toolkit enables automatic calculation and evaluation of coadsorbed configurations allowing us to sample 2,000 coadsorbed adsorbates and transition states (TSs) for a diverse set of 9 reactions on Cu(111), a challenging surface, with four possible coadsorbates. This dataset was then used to train subgraph isomorphic decision trees (SIDTs) to predict the stability and association energy of configurations. We were able to achieve mean absolute errors (MAEs) of 0.106 eV on adsorbates, 0.172 eV on TSs, and due to natural error cancellation in SIDTs for relative properties, 0.130 eV on reaction energies and 0.180 eV on activation barriers. In conclusion, we describe how to use these models to predict coverage-dependent corrections for adsorbates and TSs and demonstrate on H*, HO*, and O* comparing the generated SIDT model with an iteratively refined version.

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

Energy Barriers for H 2 Uptake and Release by Defects in Boron Nitride

Density functional theory simulations (DFT) have been used in combination with the climbing image nudged elastic band method (CI-NEB) to determine the energy barriers required to dissociate H 2 at four types of h-BN defects and evaluate the potential for reacted H species to hop to adjacent available sites. The defects were selected based on previous work in which ab initio thermodynamics was used to identify specific defects with optimal thermodynamics for reversible H 2 storage, that is, the defects would dissociate molecular H 2 while not binding H species too strongly so that hydrogen release is favorable. The four defects investigated are the N and B monovacancies, which only offer either B or N defective sites respectively, as well as larger cluster vacancy defects named 3V(1B2N) and hexagonal 6V(3B3N), which present a mixture of B and N defect sites available for H binding. While homogeneous H 2 dissociation by B and N monovacancies have energy barriers larger than 0.4 eV, we found that mixed BN terminated defects have energy barrier lower than 0.4 eV for H 2 dissociation. In particular, H 2 dissociation by the hexagonal 6V(3B3N) defect, exclusively made of Frustrated Lewis Pairs (FLPs), is almost barrierless. The energy barrier for subsequent H hopping to available adjacent B or N sites of the defect were also calculated. It was found that the hopping of H species generally involves high energy barriers ranging from 0.26 eV to 2.57 eV. In conclusion, this work completes previous thermodynamic theoretical investigations and provides a more comprehensive picture of H 2 reaction with defective h-BN for potential H 2 storage and recovery.

Sassi, Michel J. [Pacific Northwest National Labor↗