Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density functional 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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

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↗
DOE OSTI2022

Competing electronic states emerging on polar surfaces

Excess charge on polar surfaces of ionic compounds is commonly described by the two-dimensional electron gas (2DEG) model, a homogeneous distribution of charge, spatially-confined in a few atomic layers. Here, by combining scanning probe microscopy with density functional theory calculations, we show that excess charge on the polar TaO 2 termination of KTaO 3 (001) forms more complex electronic states with different degrees of spatial and electronic localization: charge density waves (CDW) coexist with strongly-localized electron polarons and bipolarons. These surface electronic reconstructions, originating from the combined action of electron-lattice interaction and electronic correlation, are energetically more favorable than the 2DEG solution. They exhibit distinct spectroscopy signals and impact on the surface properties, as manifested by a local suppression of ferroelectric distortions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Atom Equivalent Energies for the Rapid Estimation of the Heat of Formation of Explosive Molecules from Density Functional Tight Binding Theory

Atom equivalent energies have been derived from which the gas-phase heat of formation of explosive molecules can be estimated from fast, semiempirical density functional tight binding total energy calculations. The root-mean-square deviation and maximum deviation of the heats of formation from the experimental values for the set of 45 energetic molecules compiled by Byrd and Rice [ J. Phys. Chem. A, 2006, 110, 1005–1013] are 10.4 and 25.5 kcal/mol, respectively, using 4 atom equivalent energies and 7.4 and 15.0 kcal/mol, respectively, using 7 atom equivalent energies. These errors are around a factor of 2–3 larger than those obtained from density functional theory calculations but are smaller than those obtained from other semiempirical electronic structure methods. Heats of formation calculated with density functional tight binding theory using the 4 and 7 atom equivalent energies, the Byrd and Rice scheme, and the atom pair contribution method for a new set of 531 energetic molecules that contain only carbon, hydrogen, nitrogen, and oxygen are provided.

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

The Electron Thermal Conductivity of Pu and Zr Substituted $\mathcal{γ}$-U

Uranium alloys are attractive recycled nuclear fuels because of their high thermal conductivity (𝑘) and fissile density. Limited experimental studies of the 𝑘 of U-Pu-Zr alloys in the range of 15 to 20 wt% Pu and 6 to 15 wt% Zr indicate that increasing the content of either Zr or Pu tends to lower 𝑘. However, which element has the greater effect on 𝑘, and the associated mechanisms, remains unclear. Here, in this study, the electron thermal conductivity (𝑘 𝑒 ) of U-Pu-Zr compositions are calculated using density functional theory. The electronic structure is evaluated to understand the effects of plutonium (Pu) and zirconium (Zr) substitution on the 𝑘 𝑒 of 𝛾-U. Alloys of up to 37.5 at. % Pu and 37.5 at. % Zr are examined. Two methods are applied to calculate 𝑘 𝑒 ; we find that the accuracy of each method depends on the electronic and mass similarities between the solute and solvent atoms. Specifically, when the solute atom is similar in electronic structure and mass, the more accurate method is that which employs the electron relaxation time of 𝛾-U, while if the elements are dissimilar, a mixed method that mixes several parameters associated with JNW_S⁢3033426825100132 from each element in the alloy is best. The introduction of all alloying elements decreases 𝑘 𝑒 ; however, in binary compounds, Pu and Zr have different effects. Pu flattens the electronic bands but compensates for this deleterious effect by increasing electron density near the Fermi level. Zr flattens the electronic bands more severely without adding electron density near the Fermi level. Therefore, Zr decreases 𝑘 𝑒 more than Pu in binary compounds. In ternary compounds, the difference between Pu and Zr is minimal due to the phononic change from the large mass change of Zr substitution, even at 12.5 at. %. Thus, we predict that higher loadings of Pu, and potentially other actinides, can be added to U-Pu-Zr compositions for faster recycling of spent fuel without sacrificing 𝑘. We also note that these 𝑘 𝑒 calculation methods can be applied to non-fuel alloys that require 𝑘 𝑒 predictions, such as cladding, heat exchanger, and structural materials.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Unconventional pressure-induced magnetic transitions and mechanical properties in non-magnetic LaFeSi

Here, this work demonstrates that two consecutive magnetic transitions occur in non-magnetic LaFeSi compound under negative hydrostatic pressure (small volume expansion). Electronic structure properties of LaFeSi were calculated using density functional theory (DFT) to understand the origin of these magnetic transitions. Mechanical properties of LaFeSi, especially elastic properties (shear and Young’s modulus, and Poisson’s ratio), were determined by DFT distortion calculations to elucidate lattice anisotropy and stability. The antiferromagnetic transition at –7.52 GPa, predicted in our work for LaFeSi, is consistent with the reported experimentally observed antiferromagnetic ground states in CeMnSi and LaMnSi compounds having similar crystal geometry. This finding of unexpected magnetic ordering predicts that layered non-magnetic materials, like LaFeSi, may become magnetically active with lattice expansion.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Selectivity in Yttrium Manganese Oxide Synthesis via Local Chemical Potentials in Hyperdimensional Phase Space

In sharp contrast to molecular synthesis, materials synthesis is generally presumed to lack selectivity. The few known methods of designing selectivity in solid-state reactions have limited scope, such as topotactic reactions or strain stabilization. This contribution describes a general approach for searching large chemical spaces to identify selective reactions. This novel approach explains the ability of a nominally “innocent” Na 2 CO 3 precursor to enable the metathesis synthesis of single-phase Y 2 Mn 2 O 7 : an outcome that was previously only accomplished at extreme pressures and which cannot be achieved with closely related precursors of Li 2 CO 3 and K 2 CO 3 under identical conditions. By calculating the required change in chemical potential across all possible reactant-product interfaces in an expanded chemical space including Y, Mn, O, alkali metals, and halogens, using thermodynamic parameters obtained from density functional theory calculations, we identify reactions that minimize the thermodynamic competition from intermediates. In this manner, only the Na-based intermediates minimize the distance in the hyperdimensional chemical potential space to Y 2 Mn 2 O 7 , thus providing selective access to a phase which was previously thought to be metastable. Experimental evidence validating this mechanism for pathway-dependent selectivity is provided by intermediates identified from in situ synchrotron-based crystallographic analysis. Finally, this approach of calculating chemical potential distances in hyperdimensional compositional spaces provides a general method for designing selective solid-state syntheses that will be useful for gaining access to metastable phases and for identifying reaction pathways that can reduce the synthesis temperature, and cost, of technological materials.

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

Extreme sensitivity of higher-order interatomic force constants and thermal conductivity to the energy surface roughness of exchange-correlation functionals

In this Letter, we report that the fourth-order interatomic force constants (4th-IFCs) are significantly sensitive to the energy surface roughness of exchange-correlation (XC) functionals in density functional theory calculations. This sensitivity, which is insignificant for the second- (2nd-) and third-order (3rd-) IFCs, varies for different functionals in different materials and can cause misprediction of thermal conductivity by several times of magnitude. As a result, when calculating the 4th-IFCs using the finite difference method, the atomic displacement needs to be taken large enough to overcome the energy surface roughness, in order to accurately predict phonon lifetime and thermal conductivity. We demonstrate this phenomenon on a benchmark material (Si), a high-thermal conductivity material (BAs), and a low thermal conductivity material (NaCl). For Si, we find that the LDA, PBE, and PBEsol XC functionals are all smooth to the 2nd- and 3rd-IFCs but all rough to the 4th-IFCs. This roughness can lead to a prediction of nearly one order of magnitude lower thermal conductivity. For BAs, all three functionals are smooth to the 2nd- and 3rd-IFCs, and only the PBEsol XC functional is rough for the 4th-IFCs, which leads to a 40% underestimation of thermal conductivity. For NaCl, all functionals are smooth to the 2nd- and 3rd-IFCs but rough to the 4th-IFCs, leading to a 70% underprediction of thermal conductivity at room temperature. In conclusion, with these observations, we provide general guidance on the calculation of 4th-IFCs for an accurate thermal conductivity prediction.

42 ENGINEERING↗
DOE OSTI2021

High-throughput search for magnetic topological materials using spin-orbit spillage, machine learning, and experiments

Magnetic topological insulators and semi-metals have a variety of properties that make them attractive for applications including spintronics and quantum computation. Here, we use systematic high-throughput density functional theory calculations to identify magnetic topological materials from the ≈ 40000 three-dimensional materials in the JARVIS-DFT database. First, we screen materials with net magnetic moment > 0.5 μB and spin-orbit spillage > 0.25, resulting in 25 insulating and 564 metallic candidates. The spillage acts as a signature of spin-orbit induced band-inversion. Then, we carry out calculations of Wannier charge centers, Chern numbers, anomalous Hall conductivities, surface bandstructures, and Fermi-surfaces to determine interesting topological characteristics of the screened compounds. We also train machine learning models for predicting the spillage, bandgaps, and magnetic moments of new compounds, to further accelerate the screening process. We experimentally synthesize and characterize a few candidate materials to support our theoretical predictions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Structure, Bonding, and Vibrational Dynamics of a Triamine High Energy Density Material under Pressure

High energy density materials have complex intermolecular interactions which influence their stability and performance. We used a combination of synchrotron X-ray diffraction, synchrotron infrared spectroscopy, and Raman vibrational spectroscopy, supplemented by density functional theory calculations, to probe pressure-induced changes in structure and intermolecular interactions of 1H,4'H-[3,3'-bis(1,2,4-triazole)]-4',5,5'-triamine as a model high energy density material up to 40 GPa. We find that compression of the triamine is accompanied by increased intermolecular interactions that give rise to an interesting evolution of the structure and bonding with pressure. Analysis of the equation of state determined from the X-ray diffraction indicates a change in compression mechanism near 19 GPa consistent with changes in vibrational spectra that provide evidence for a structural rearrangement associated with changes in hydrogen bonding near that pressure. As a result, the overall compressional behavior calculated theoretically agrees with that observed experimentally though differences are found that indicate the need for improved treatment of the intermolecular interactions including hydrogen bonding under pressure.

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

Interpreting Tafel behavior of consecutive electrochemical reactions through combined thermodynamic and steady state microkinetic approaches

Assessing the reaction pathway of multi-electron-transfer reactions is an essential yet difficult task for the rational design of electrocatalysts. In this work, we develop a heuristic approach that combines thermodynamic adsorption energetics calculated through density functional theory with microkinetic modeling using the steady state approximation to interpret the potential-dependent Tafel behavior of consecutive electrochemical reactions. In doing so, we introduce a kinetic framework for ab initio calculations that ensures self-consistent adsorption energetics based on kinetically limited adsorbate coverages. The approach is applied to experimental results on CoO x (OH) 2-x single crystal electrocatalyst particles yielding coverage dependent mechanistic information and identification of the rate-limiting step with standard rate constants for the oxygen evolution reaction on the ($11\bar{2}0$) surfaces of the β-Co(OH) 2 , β-CoOOH, and CoO 2 bulk phases. Finally, this generalizable method enables catalyst benchmarking based on determining the active species involved and associated intrinsic reaction rate constants in consecutive multi-electron-transfer reactions.

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

Effect of chemical substitution and external strain on phase stability and ferroelectricity in two dimensional M 2 CT 2 MXenes

Two dimensional ferroelectric materials are gaining increasing attention for use in ultrathin electronic devices owing to the presence of a spontaneous polarization down to one or two monolayers. However, such materials are difficult to identify, especially those with out-of-plane electric polarizations. Previous work predicted that a metastable ferroelectric phase exists in the 2D MXene Sc 2 CO 2 , while further studies have predicted that this phase exists in other MXene chemistries. However, questions remain about the origin of ferroelectricity, the stability of this phase relative to other competing phases, and the effect of external stimuli in these materials. In this work, we use density functional theory calculations to investigate 12 M 2 CT 2 MXenes (M = transition metal, T = surface terminating group) and determine which have the ferroelectric phase as their ground state. We compute these materials’ polarizations, densities of states, phonon band structures, Bader charges, and Born effective charges in the ferroelectric phase to elucidate the reasons for its stabilization. We demonstrate that this ferroelectric phase can be preferentially stabilized in non-ferroelectric MXenes through full chemical substitution of Sc or O, alloying of the Sc sites, or application of epitaxial strain. Finally, we show that these materials have excellent piezoelectric properties as well. This work provides a detailed understanding of ferroelectric MXenes and show how the number of 2D ferroelectric materials can be increased through chemical substitution or application of external stimuli.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Understanding initial zeolite oligomerization steps with first principles calculations

Zeolites are porous aluminosilicate materials that find various applications in the chemical industry in separations, catalysis, ion exchange, and so forth. However, despite their widespread use, the reaction mechanisms occurring during zeolite growth are still unclear. In the study, we use density functional theory calculations to gain insights into the thermodynamics of oligomerization, which constitute the initial steps of zeolite growth. By taking into consideration solvent and temperature effects, our results demonstrate that the growth of aluminosilicate systems is significantly more exothermic than their pure silicate counterparts. Under pH neutral conditions, water prefers to dissociate on the early-growth-stage aluminosilicate complexes rather than desorb, thus generating potential Brønsted acid sites on the oligomers. Additionally, (alumino)silicate growth pathways are evaluated in the presence of Na + cation, as well as the Ca 2+ cation for the pure silicate pathway. The presence of cations increases the exothermicity of growth, with Ca 2+ exhibiting the most energetically favorable growth environment for the silicate systems. Importantly, we demonstrate through reaction extent analysis that the presence of cations modulates the speciation of the formed oligomers, with Na + favoring linear species in addition to the generally preferred cyclic ones. Overall, this work provides a fundamental understanding of the thermodynamics of complex reaction paths that occur during early stages of zeolite growth and suggests that the initial growth steps can have significant impact on the final zeolite structure.

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

Targeting One- and Two-Dimensional Ta–Te Structures via Nanotube Encapsulation

Fine control over material synthesis on the nanoscale can facilitate the stabilization of competing crystalline structures. In this report we demonstrate how carbon nanotube reaction vessels can be used to selectively create one-dimensional TaTe 3 chains or two-dimensional TaTe 2 nanoribbons with exquisite control of the chain number or nanoribbon thickness and width. Transmission electron microscopy and scanning transmission electron microscopy reveal the detailed atomic structure of the encapsulated materials. Complex superstructures such as multichain spiraling and apparent multilayer moirés are observed. The rare 2H phase of TaTe 2 (1H in monolayer) is found to be abundant as an encapsulated nanoribbon inside carbon nanotubes. The experimental results are complemented by density functional theory calculations for the atomic and electronic structure, which uncovers the prevalence of 2H-TaTe 2 due to nanotube-to-nanoribbon charge transfer and size confinement. Calculations also reveal new 1T' type charge density wave phases in TaTe 2 that could be observed in experimental studies.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Fast and stable deep-learning predictions of material properties for solid solution alloys

We present a novel deep learning (DL) approach to produce highly accurate predictions of macroscopic physical properties of solid solution binary alloys and magnetic systems. The major idea is to make use of the correlations between different physical properties in alloy systems to improve the prediction accuracy of neural network (NN) models. We use multitasking NN models to simultaneously predict the total energy, charge density and magnetic moment. These physical properties mutually serve as constraints during the training of the multitasking NN, resulting in more reliable DL models because multiple physics properties are correctly learned by a single model. Two binary alloys, copper–gold (CuAu) and iron–platinum (FePt), were studied. Our results show that once the multitasking NN's are trained, they can estimate the material properties for a specific configuration hundreds of times faster than first-principles density functional theory calculations while retaining comparable accuracy. We used a simple measure based on the root-mean-squared errors to quantify the quality of the NN models, and found that the inclusion of charge density and magnetic moment as physical constraints leads to more stable models that exhibit improved accuracy and reduced uncertainty for the energy predictions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

ReaxFF Parameter Set for Boron Clusters and Icosahedral Boron Crystals: Comparison with Density Functional Theory and Machine-Learning Potentials

Icosahedral boron materials, which include regular icosahedra of 12 boron atoms have gained increasing attention due to their potential applications as superhard materials, semiconductors, and energy storage media. However, the synthesis of high quality crystals of these materials has been a major barrier to the development of these applications. To enable computational prediction of synthesis conditions yielding high-quality icosahedral boron crystals, herein we tested and refined a set of ReaxFF parameters for the nucleation and growth of such crystals. We focused on matching the relative energies of small boron clusters obtained by density functional theory since such small clusters and similar motifs are likely present in crystal nuclei and at the interface of growing crystals. Using a training set of B 80 clusters, including a low-energy core–shell structure containing a B 12 icosahedron core and a high-energy single-shell structure produced in preliminary ReaxFF simulations, the ReaxFF parameter set was refined to better reproduce energies calculated by density functional theory (DFT). Among existing ReaxFF parameter sets and the machine-learning interatomic potentials MACE-MP-0, MACE-MP-0b3, MACE-MPA-0, PFP v7.0.0, and SevenNet-MF-ompa, only our new parameter set and PFP v7.0.0 correctly ranked these B 80 clusters. This refinement led to improved agreement with DFT for a test set of 58 clusters consisting of 8–103 boron atoms. Furthermore, our refined parameter set yielded greater local icosahedral structure than the previously existing ReaxFF parameter set for larger scale simulations of crystallization from supercooled liquid boron. Additionally, simulations of solid boron in contact with molten nickel using our refined ReaxFF parameters yielded a boron solubility value that agrees moderately well with experimental expectations, while the previous boron parameters gave a value that was much too low.

boron↗
DOE OSTI2025

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↗
DOE OSTI2025

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↗
DOE OSTI2020

Rapid generation of optimal generalized Monkhorst-Pack grids

In this report computational modeling of the properties of crystalline materials has become an increasingly important aspect of materials research, consuming hundreds of millions of CPU-hours at scientific computing centers around the world each year, if not more. A routine operation in such calculations is the evaluation of integrals over the Brillouin zone. We have previously demonstrated that performing such integrals using generalized Monkhorst-Pack k-point grids can roughly double the speed of these calculations relative to the widely-used traditional Monkhorst-Pack grids. However the generation of optimal generalized Monkhorst-Pack grids is not implemented in most software packages due to the computational cost and difficulty of identifying the best grids. To address this problem, we present new algorithms that allow rapid generation of optimal generalized Monkhorst-Pack grids on the fly. We demonstrate that the grids generated by these algorithms are on average significantly more efficient than those generated using existing algorithms across a range of grid densities. For grids that correspond to a real-space supercell with at least 50 Å between lattice points, which is sufficient to converge density functional theory calculations within 1 meV/atom for nearly all materials, our algorithm finds optimized grids in an average of 0.19 s on a single processing core. To facilitate the widespread adoption of this approach, we present new open-source tools including a library designed for integration with third-party software packages.

36 MATERIALS SCIENCE↗