Engineering PapersSearch

SEARCH · Engineering Papers

Results for “density functional theory calculations”

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 55 records · Page 3

Using Density Functional Theory (DFT) for the Calculation of Atomization Energies

The calculation of atomization energies using density functional theory (DFT), using the B3LYP hybrid functional, is reported. The sensitivity of the atomization energy to basis set is studied and compared with the coupled cluster singles and doubles approach with a perturbational estimate of the triples (CCSD(T)). Merging the B3LYP results with the G2(MP2) approach is also considered. It is found that replacing the geometry optimization and calculation of the zero-point energy by the analogous quantities computed using the B3LYP approach reduces the maximum error in the G2(MP2) approach. In addition to the 55 G2 atomization energies, some results for transition metal containing systems will also be presented.

Bauschlicher, Charles W., Jr.

Selective Hydrogenation of Furfural Acetone over a Cu Catalyst: Combined Theoretical and Experimental Study

The selective hydrogenation of biomass-derived hydroxymethyl furfural (HMF)-acetone-HMF (HAH) presents an alternative route to producing highly functionalized polyesters and polyurethanes. While HAH undergoes furan ring hydrogenation over Pd, Ru, and Ni catalysts, furan ring hydrogenation is not observed over Cu. Herein, we combine reaction kinetics experiments and density functional theory calculations to elucidate the selective hydrogenation behavior of HAH over Cu catalysts. We identified furfural acetone (FAc) as a suitable surrogate for modeling HAH hydrogenation over Cu surfaces and performed reaction kinetics experiments between temperatures of 313–393 K and a H 2 partial pressure of 55 bar. Similar to the behavior of HAH, hydrogenation of FAc follows a consecutive two-step hydrogenation pathway over Cu and does not undergo furan ring hydrogenation. The apparent activation energy barriers for hydrogenation of the aliphatic double bond (0.58 eV) and carbonyl (0.43 eV) of FAc measured in a continuous flow reactor setup are consistent with those reported in prior studies for batch HAH hydrogenation. Reaction orders with respect to each reactant, including H 2 and FAc, were determined to be nearly one. Density functional theory (DFT; GGA-PBE-D3) calculations on Cu(111) showed that the hydrogenation of the aliphatic double bond of FAc is more facile than the hydrogenation of the furan ring, which displays weak interactions with the Cu surface. We determined an apparent activation energy barrier for FAc hydrogenation (0.58 eV) that agreed with our DFT predictions (highest barrier for FAc hydrogenation of 0.57 eV). Finally, our DFT calculations further show that weak interactions between the furan ring and Cu surface are responsible for the selective hydrogenation behavior.

Biomass

Planar Defect Layers Template a High-Pressure InBi Polymorph

The short- and long-range order of III–V materials under high pressure has long been the subject of debate, with advancements in structural characterization leading to significant revisions to the accepted structural models. Despite these revisions, previous high-pressure structural assignments in the In–Bi system include the site-disordered β-Sn structure type, a structure type demonstrated to be nonexistent in analogous III–V systems. While X-ray diffraction is consistent with site disordering in InBi at high pressure, cluster expansion calculations indicate that disordering requires temperatures above 3000 K. Here, we propose InBi as a model material for studying unique high-pressure planar defects due to its highly anisotropic stress-dependent properties and structure. Specifically, we identify two sets of planar defects that mimic the diffraction pattern of a site disordered β-Sn structure type and are compatible with the calculated disorder barrier. We derive these defects by symmetry relations over crystallographic transitions. Density functional theory calculations of the proposed defects suggest that these defects are stabilized by diminishing interlayer separations with pressure. Further, we find that one of the proposed defects closely resembles a bulk high-pressure phase of InBi, InBi-ϵ, and we assert that the proposed defects order upon heating, acting as a template for InBi-ϵ growth. The proposed defects and their electronic structure provide a basis for the trend of superconducting critical temperature with increasing pressure. These methods for identifying defects are generalizable to other materials with reports of site disorder at high pressure, prompting a broader search for related high-pressure defects.

36 MATERIALS SCIENCE

Defluorination Mechanisms and Real-Time Dynamics of Per- and Polyfluoroalkyl Substances on Electrified Surfaces

Per- and polyfluoroalkyl substances (PFAS) are persistent environmental contaminants found in groundwater sources and a wide variety of consumer products. In recent years, electrochemical approaches for the degradation of these harmful contaminants have garnered a significant amount of attention due to their efficiency and chemical-free modular nature. However, these electrochemical processes occur in open, highly non-equilibrium systems, and a detailed understanding of PFAS degradation mechanisms in these promising technologies is still in its infancy. To shed mechanistic insight into these complex processes, we present the first constant-electrode potential (CEP) quantum calculations of PFAS degradation on electrified surfaces. These advanced CEP calculations provide new mechanistic details about the intricate electronic processes that occur during PFAS degradation in the presence of an electrochemical bias, which cannot be gleaned from conventional density functional theory calculations. We complement our CEP calculations with large-scale ab initio molecular dynamics simulations in the presence of an electrochemical bias to provide time scales for PFAS degradation on electrified surfaces. Taken together, our CEP-based quantum calculations provide critical reaction mechanisms for PFAS degradation in open electrochemical systems, which can be used to prescreen candidate material surfaces and optimal electrochemical conditions for remediating PFAS and other environmental contaminants.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Structural Evidence of Interanionic Hydrogen Bonding in Phosphoric Acid Solutions

Interanionic hydrogen bonding (IAHB) is a noncovalent interaction between like-charged ions that challenges conventional electrostatic understanding. This study provides direct structural evidence of IAHB in concentrated aqueous phosphoric acid (PA) solutions, which exhibit >60% dissociation under these conditions. Oxygen K-edge X-ray absorption fine structure spectroscopy, combined with electron affinity time-dependent density functional theory calculations, reveals the formation of stable, cyclic phosphate-phosphate IAHB dimers at PA concentrations ≥7 M. Extended X-ray absorption fine structure data show distinct long-range ordering consistent with these dimers, and near-edge X-ray absorption fine structure spectra confirm a concentration-dependent transition from monomeric to dimeric species. Energy decomposition analysis through density functional theory shows that the formation of solution-phase IAHB is energetically favored and is attributed to polarization of and the charge transfer between the two fragments driven by the surrounding solvent molecules, in addition to the permanent electrostatics. These findings offers crucial structural insights into the H-bonded networks in concentrated PA, highlighting the critical role of solvent in facilitating anion–anion association.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Electronic Structural Optimization of the Air Electrodes for Reversible Protonic Ceramic Electrochemical Cells via IIIA Cation Doping

Reversible protonic ceramic electrochemical cells (R-PCECs) have emerged as a novel technology for clean and efficient energy generation and storage. Improving the oxygen and proton conduction characteristics and stability of air electrodes at intermediate temperatures is crucial for achieving a high performance. Herein, we optimize the electronic structure of a state-of-the-art PrBa 0.8 Ca 0.2 Co 2 O 5+δ (PBCC) air electrode via doping IIIA cations (Al 3+ , Ga 3+ , and In 3+ ). Also, it is shown that PrBa 0.8 Ca 0.2 Co 1.9 Ga 0.10 O 5+δ (PBCCGa 0.10 ) exhibits improved oxygen reaction activity and hydration capability. Density functional theory calculations confirm that Ga doping provides the most favorable electronic structure. An R-PCEC with PBCCGa 0.10 achieves a peak power density of 2.21 W cm -2 and a current density of −4.64 A cm -2 at 1.3 V at 650 °C. Additionally, the PBCCGa 0.10 electrode performs good operational stability in FC mode (for about 100 h), EC mode (for about 100 h), and reversible cyclic testing (over 200 h) at 600 °C.

30 DIRECT ENERGY CONVERSION

Machine Learning the COSMO Model for Predicting Thermodynamics of Electrolyte Mixtures

Bottom-up design of electrolyte mixtures for battery systems requires predicting macro thermodynamic properties from molecular constituents. For instance, molten salt electrolyte batteries require conditions far above room temperature to operate. Therefore, discovering mixtures with increasingly lower eutectic melting points is desirable. A model that can approximate chemical activity is a valuable tool to search through the vast compositional design space. Machine learning can predict properties of materials such as vibrational free energies, electronic energy gaps, and thermal conductivities. Moreover, they can learn physical models such as interatomic potentials. The COSMO-SAC model uses theory and empirical parameterization to predict liquid-vapor and liquid-solid properties using first-principles calculations. However, obtaining activity coefficients required for parameterizing the COSMO-SAC model is costly and limited to a select chemical space. In this work, we explored if machine learning methods could improve the COSMO-SAC model and bridge density functional theory calculations to liquid phase thermodynamic properties. Our data-driven approach uses existing databases for sigma-profiles of organic solvents and reconciles their methodological differences via ensemble averaging. First, an optimal machine learning model is constructed for each dataset. Our machine learning algorithms use the sigma-profile as an input feature to predict binary mixtures' activity coefficients using multi-output regression. Each dataset uses different choices of functionals, methods, and basis sets. Therefore, our ensemble model attempts to predict corrected activity coefficients given the combination of all the model outputs. The activity coefficients used for training are generated using the COSMO-SAC model. This approach enables the extraction of meaningful information from the existing datasets to improve the COSMO-SAC model for obtaining thermodynamic properties of electrolyte mixtures. With the liquid phase activities, we can identify electrolyte mixtures that meet desired phase equilibria conditions.

Thermodynamics

Predicting phase transitions in PbTi⁢O 3 using zentropy through quasiharmonic phonon calculations

According to x-ray diffraction (XRD) measurements, PbTi⁢O 3 undergoes a phase transition from a tetragonal ferroelectric (FE) phase to a cubic paraelectric phase at 763 K. However, x-ray absorption fine-structure (XAFS) measurements indicate that PbTi⁢O 3 is locally tetragonal even after the phase transition. The difference in these results is because XAFS measurements can probe local features of a structure, while XRD averages over such local features. For both measurements to be consistent, PbTi⁢O 3 is macroscopically cubic but locally tetragonal after the phase transition. Despite this, most models, such as the Landau-Ginsburg-Devonshire theory and effective Hamiltonians, are still unable to explain this phenomenon. Moreover, these methods involve model parameters fitted to experimental or theoretical data and do not consider other tetragonal configurations, such as domain walls, to predict the phase transition. In our previous study, we used our zentropy approach to predict the phase transition by considering the tetragonal FE ground-state configuration and the tetragonal 90° and 180° domain wall configurations with their total energies at 0 K. Here, in this paper, the Helmholtz energies of the three configurations are obtained from density functional theory calculations through energy-volume curves and phonon calculations. The predicted phase transition temperature using the meta-GGA 𝑟 2⁢ SCAN and revised multiplicities of configurations is 716 K, showing good agreement with the experimental value of 763 K.

36 MATERIALS SCIENCE

Unveiling the Role of Lithium Iodide in Stabilizing Solid Interfaces in All-Solid-State Li Metal Batteries

A critical challenge in all-solid-state lithium metal batteries (ASSLMBs) is achieving a stable interface between the lithium metal anode and the solid electrolyte. Leveraging its success in Li/I 2 batteries, lithium iodide has garnered significant attentions for its potential to enhance interfacial stability and overall cell performance in ASSLMBs. Here, we elucidate the role of lithium iodide in stabilizing the solid interface in all-solid-state Li metal batteries with a Li argyrodite electrolyte, particularly focusing on its influence on lithium deposition behavior and interfacial evolution. Through in situ optical imaging, we demonstrate more uniform lithium deposition on an iodide-contained argyrodite electrolyte compared to a chloride-based counterpart. Complementary density functional theory calculations attribute improved lithium plating behavior to the enhanced lithiophilicity and better ionic conductivity of lithium iodide at the solid interface, effectively reducing localized current density. In conclusion, these findings provide useful insights into the mechanisms through which lithium iodide enhances the interfacial stability in ASSLMBs.

Anions

Size-dependent attraction of Cu solutes to clusters formed at Ag grain boundaries

We report a size-dependent solute clustering mechanism at grain boundaries in a sputtered ultrafine-grained Ag-Cu alloy, where large Cu clusters form despite weak individual solute-solute interactions. X-ray diffraction confirms limited Cu solubility in the Ag matrix, while scanning transmission electron microscopy reveals Cu clustering at both ordinary GBs and GB junctions. Density functional theory calculations show that 12-atom Cu clusters are energetically preferred, while smaller three-atom clusters are significantly less stable. Additional calculations demonstrate a marked increase in solute-cluster attraction energy with cluster size. As a result, these findings point to a previously unrecognized pathway for grain-boundary solute clustering in immiscible systems, driven by collective solute-cluster interactions, with implications for segregation behavior and stability in nanocrystalline and ultrafine-grained alloys.

Alloys

Near-Infrared Emission from Sulfur Heteroatom Defects in Single-Walled Carbon Nanotubes

Quantum defects introduced into single-walled carbon nanotubes (SWCNTs) that support near-infrared emission are promising candidates for many applications. Here, we present a type of quantum defect in which sulfur atoms are incorporated into the crystal lattice of SWCNTs through their reactions with a mild reducing agent, sodium dithionite. When radical anions generated by sodium dithionites react with (6,5)-SWCNTs, two new emission features emerge at 1120 and 1250 nm. Comparing these spectral features with density functional theory calculation results, we attribute them to defects with episulfide or thioether structures. Here, by investigating the influences of reaction conditions including reactant concentration, light irradiation, and solvent type, we identify the reaction mechanism between SWCNTs and sodium dithionite and achieve good control over defect densities. Due to the intrinsic reactivity of sulfur species, the sulfur-containing quantum defects discovered here may facilitate the introduction of more complex and functionally diverse moieties into SWCNTs.

Anions

Low-temperature etching of silicon oxide and silicon nitride with hydrogen fluoride

Etching of high aspect ratio features into alternating SiO 2 and SiN layers is an enabling technology for the manufacturing of 3D NAND flash memories. In this paper, we study a low-temperature or cryo plasma etch process, which utilizes HF gas together with other gas additives. Compared with a low-temperature process that uses separate fluorine and hydrogen gases, the etching rate of the SiO 2 /SiN stack doubles. Both materials etch faster with this so-called second generation cryo etch process. Pure HF plasma enhances the SiN etching rate, while SiO 2 requires an additional fluorine source such as PF 3 to etch meaningfully. Further, the insertion of H 2 O plasma steps into the second generation cryo etch process boosts the SiN etching rate by a factor of 2.4, while SiO 2 etches only 1.3 times faster. We observe a rate enhancing effect of H 2 O coadsorption in thermal etching experiments of SiN with HF. Ammonium fluorosilicate (AFS) plays a salient role in etching of SiN with HF with and without plasma. AFS appears weakened in the presence of H 2 O. Density functional theory calculations confirm the reduction of the bonding energy when NH 4 F in AFS is replaced by H 2 O.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Nanoscale modulation of flat bands via controllable charge density wave defects in 4⁢𝐻⁡𝑏−Ta⁢S 2

Electron correlation is a main driver of exotic quantum phases and their interplay. The 4H b-TaS 2 system, possessing an intrinsic heterostructure of 1T- and 1H -TaS 2 monolayers, offers a unique opportunity to control electron correlation by distorting the atomic lattice or tuning interlayer coupling. Here, we investigated intrinsically deformed charge-density waves (CDWs) in the 1T layer of 4H b-TaS 2 to elucidate and control their effects on flat bands using scanning tunneling microscopy and spectroscopy (STM/S) combined with first-principles calculations. We identified two types of CDW defects: Type 1 has structural distortion and locally suppressed flat bands, while Type 2 features an increased flat band filling factor of intact CDW structure. Density functional theory calculations indicate that a sulfur vacancy in the 1T layer distorts the CDW structure and gives rise to a Type 1, whereas a sulfur vacancy in the 1H layer reduces the interlayer charge transfer and leads to a Type 2. Furthermore, we demonstrated creating and erasing individual CDW defects via STM manipulation. Here, our findings provide a pathway to not only tune flat bands but also selectively manipulate the interaction between CDW, the atomic lattice, and interlayer coupling in strongly correlated systems with atomic precision.

Charge density waves

Understanding Inlet Concentration Effects on the Electrocatalytic Conversion of CO 2 to Formic Acid in Gas-Fed Electrolyzers

The electrochemical CO 2 reduction reaction (CO2RR) to produce value-added products remains a developing technology for utilizing waste CO 2 streams. Most device-level CO2RR studies use pure CO 2 gas feeds; however, the effect of dilute CO 2 on the electrolyzer performance is an important consideration for large-scale electrolyzer operation, single-pass conversion, and real-world CO 2 source utilization. This work investigates the effect that the CO 2 concentration has on the performance of formic acid (HCOOH) producing tin oxide (SnO 2 ) and bismuth oxide (Bi 2 O 3 ) catalysts in an electrolyzer device setting. Surprisingly, SnO2 demonstrated an approximately 20% increase in HCOOH selectivity (Faradaic efficiency) when the CO 2 concentration decreased from 100 to 20%. In contrast, Bi 2 O 3 consistently demonstrated high selectivity toward HCOOH across the same CO 2 concentration range. The effects of the CO 2 concentration on selectivity were further investigated with half-cell experiments and in situ Raman spectroscopy, which revealed dynamic changes in the cathodic overpotential and chemical state of the catalyst that depended on the CO 2 concentration. Density functional theory calculations showed how changes in the surface oxidation state of Sn, varying from fully oxidized SnO 2 to metallic Sn(0), affect the thermodynamic barriers of the three main observed products: HCOOH, CO, and H 2 . Our results indicate that dilute CO 2 concentrations required larger cathodic overpotentials to sustain a fixed current density, which, in turn, pushed the Sn-based catalyst toward a more reduced surface that was favorable to HCOOH formation. On the other hand, the Bi-based catalyst remained in a metallic state at CO2RR-relevant potentials and demonstrated a consistent product selectivity regardless of CO 2 concentration. These findings highlight how varying the CO 2 inlet gas concentrations affects the chemical state of catalysts and the resulting performance metrics.

42 ENGINEERING

Single-crystalline orthorhombic GdAlGe as a rare-earth magnetic Dirac nodal-line metal

Crystal engineering is a method for discovering new quantum materials and phases, which may be achieved using external pressure or strain. Chemical pressure is unique in that it generates internal pressure perpetually to the lattice. As an example, GdAlSi from the rare-earth (𝑅) 𝑅⁢Al⁢𝑋 (𝑋=Si or Ge) family of Weyl semimetals is considered. Replacing Si with the larger isovalent element Ge creates sufficiently large chemical pressure to induce a structural transition from the tetragonal structure of GdAlSi, compatible with a Weyl semimetallic state, to an orthorhombic phase in GdAlGe, resulting in an inversion-symmetry-protected nodal-line metal. We find that GdAlGe hosts an antiferromagnetic ground state with two successive orderings, at 𝑇 N ⁢1=35K and 𝑇 N ⁢2=30K. In-plane isothermal magnetization shows a magnetic field induced metamagnetic transition at 6.2 T for 2 K. Furthermore, electron-hole compensation gives rise to a large magnetoresistance of ∼100% at 2 K and 14 T. Angle-resolved photoemission spectroscopy measurements and density functional theory calculations reveal a Dirac-like linear band dispersion over an exceptionally large energy range of ∼1.5eV with a high Fermi velocity of ∼10 6 m/s, a rare feature not observed in any magnetic topological materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND

Computational Investigation of a CO 2 Conversion Strategy via Diels–Alder Reaction in a Carbon Capture Solvent

Molecular-level insights into reactive separations are crucial for the design of new conversion pathways of carbon dioxide (CO 2 ). This work explores a postulated pathway that directs CO 2 to undergo inverse-electron-demand Diels–Alder reactions to produce heterocycles using the CO 2 chemically fixed on water-lean solvent molecules. Density functional theory calculations are applied to evaluate the lowest unoccupied molecular orbital (LUMO) energies of three types of reactants (1,3-butadiene, 1,3-cyclohexadiene, and 1,2,4,5-tetrazine) with various functional substituents. These calculations also provide a data set (5.8k data) for developing a machine learning model to efficiently predict LUMO energies. A computational screening of LUMO energies for an additional 47k diene and tetrazine candidates is performed, and a list of candidates with lowered LUMO energies by electron-withdrawing substituents is provided. These candidates are further examined by their reaction energy barriers computed from the interatomic potential or density functional theory. Two major energy barriers are identified, one for the proton transfer within the water-lean solvent and the other for the CO 2 transfer from the solvent molecule to the reactant candidate (diene or tetrazine). The functional substituents have a more significant impact on the second barrier but a very slight one on the first barrier. This exploratory work demonstrates a new possibility for guiding experimental efforts toward the chemical conversion of fixated CO 2 to value-added compounds.

Chemical reactions

Spin qubit properties of the boron-vacancy/carbon defect in the two-dimensional hexagonal boron nitride

Spin qubit defects in two-dimensional materials have a number of advantages over those in three-dimensional hosts including simpler technologies for defect creation and control, as well as qubit accessibility. In this work, we select the V B C B defect in the hexagonal boron nitride (hBN) as a possible optically controllable spin qubit and explain its triplet ground state and neutrality. In this defect a boron vacancy is combined with a carbon dopant substituting the closest boron atom to the vacancy. Our density-functional-theory calculations confirmed that the system has dynamically stable spin triplet and singlet ground states. As revealed from our linear response GW calculations, the spin-sensitive electronic states are localized around the three undercoordinated N atoms and make local peaks in the density of electronic states within the bandgap. Using the triplet and singlet ground state energies, as well as the energies of the optically excited states, obtained from solution to the Bethe–Salpeter equation, we construct the spin-polarization cycle, which is found to be favorable for the spin qubit initialization. The calculated zero-field splitting parameters ensure that the splitting energy between the spin projections in the triplet ground state is comparable to that of the known spin qubits. We thus propose the V B C B defect in hBN as a promising spin qubit.

2D BN

Hollow-structured Ni-N-C catalysts for highly selective CO 2 electroreduction

Atomically dispersed single-atom catalysts have emerged as promising non-precious catalyst alternatives to expensive Ag and Au catalysts for electrochemical CO 2 reduction reaction (CO 2 RR). In particular, nickel-nitrogen-carbon (Ni-N-C) catalysts have demonstrated a high faradaic efficiency (FE) toward CO formation at low overpotentials. Nonetheless, the exact nature of Ni active sites under CO 2 RR remains elusive and conventional Ni-N-C catalysts are limited by microporosity and low density of Ni single atoms, hindering performance in CO 2 electrolyzers. Here, we report the synthesis of hollow-structured Ni-N-C ( hs -Ni-N-C) catalysts via a post-synthesis modification (PSM) strategy using partial ligand exchange of 2-methylimidazole with 3-amino-1,2,4-triazole. This approach enables the formation of a hollow structure, resulting in more than a twofold increase in Ni atom density compared to regular Ni-N-C (r-Ni-N-C). In a zero-gap CO 2 electrolyzer, the optimized hs -Ni-N-C allows for achieving an FE CO of 97% at a current density of > 100 mA cm⁻ 2 , while maintaining high CO selectivity with stable performance over 100 h at 2.5 V. hs-Ni-N-C shows a more than sevenfold increase in the CO partial current density relative to r-Ni-N-C resulting from the combined effects of a higher density of Ni single-atom sites, improved kinetics, and lower transport resistance under the operating conditions, as indicated by electrochemical impedance spectra and distribution of relaxation times analysis. Operando high energy-resolution X-ray absorption spectroscopy (XAS) reveals that atop-bonded CO on Ni single sites induces dynamic transformations of the Ni–N coordination environment, leading to a symmetric coordination structure of hs -Ni-N-C. Under CO 2 RR, the catalysts undergo a more pronounced structural change and form a minor fraction of Ni nanoparticles. Density functional theory calculations are consistent with the XAS results and provide molecular insights showing that the interplay between protonation and CO adsorption leads to adsorbate-induced restructuring of the Ni single atom. This work demonstrates the synergistic role of hollow structure and high-density Ni atoms in governing CO 2 RR selectivity and provides mechanistic insights into the structural dynamics of single-atom catalysts under operating conditions.

36 MATERIALS SCIENCE