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

Discovery of single-atom alloy catalysts for CO 2 -to-methanol reaction by density functional theory calculations

The transformations of CO 2 molecules into valuable products are of increasing interest due to the negative impact of anthropogenic CO 2 emissions on global warming. The CO 2 -to-methanol hydrogenation is an economically profitable reaction of carbon fixation, but it still steps away from widespread industrialization because of the lack of efficient and selective catalysts. Recently, single-atom alloy (SAA) catalysts have been developed to work remarkably in CO 2 hydrogenation reactions. Doping isolated single atoms into metallic catalyst can dramatically alter the catalytic performance of the host. Here, we have performed a screening discovery on Ru and 6 RuX (X = Fe, Co, Ni, Cu, Ir and Pt) SAAs using density functional theory (DFT) computations. We considered 13 possible elementary reactions in 4 possible reaction pathways on Ru and all RuX surfaces. In the computed mechanisms, we found that the formation of *H 2 COOH and *HCOO intermediates plays a critical role in determining catalysts’ activities. Doping Co and Pt isolated single atoms into Ru surface can thermodynamically and kinetically facilitate these intermediates formation processes, eventually promoting the production of methanol. The combination of weak binding and enhanced charge redistribution on RuCo and RuPt surfaces give them improved catalytic activities over pure Ru. This work will ultimately facilitate the discovery and development of SAAs for CO 2 to methanol, serving as guidance to experiments and theoreticians alike.

25 ENERGY STORAGE↗
DOE OSTI2021

Magnetization distribution in Cu 0.6 Mn 2.4 Ge 2 ferromagnet from polarized and non-polarized neutron powder diffraction aided by density-functional theory calculations

The crystal structure and magnetic properties of Cu 0.6 Mn 2.4 Ge 2 have been re-investigated by a combination of extensive magnetic measurements and neutron scattering experiments, aided by electronic structure calculations. The material is found to be a soft ferromagnet with the ordering temperature T C = 316 K. The magnetocaloric effect evaluated from field-dependent magnetization isotherms is equal to 1.2 J/(kg·K) and 2.5 J/(kg·K) under the maximum applied magnetic field of 2 T and 5 T, respectively. The compound crystallizes in the hexagonal space group P6 3 /mmc. A complex structural disorder necessitated testing of several disorder models against the results of non-polarized and polarized neutron scattering experiments and magnetization measurements. Simulations at the density-functional theory level were also performed to identify the most robust solution that properly described the data observed. The final magnetic structure model reveals non-equal magnetic moments on the Mn1 and Mn2 atoms (2.29(9) µ B and 2.7(1) µ B , respectively) and the presence of vacancies and minor Cu substitution defects in both Mn sites. The work demonstrates how the non-polarized and polarized neutron scattering methods can be combined with electronic structure calculations to establish the microscopic structure of magnetic materials with complex crystallographic disorder.

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

Local order of orthorhombic weberite-type Y 3 TaO 7 as determined by neutron total scattering and density functional theory calculations

Weberite-type oxides are a family of oxides with A 3 BO 7 stoichiometry that can adopt different space groups depending on chemical composition. There is a discrepancy in previous studies as to whether Y 3 TaO 7 is the orthorhombic C222 1 or the Ccmm space group. Here, we describe the short- and long-range structural properties of weberite-type Y 3 TaO 7 using neutron total scattering data, which has a high sensitivity to the oxygen sublattice. Simultaneous analysis of both short- and long-range structural data via conventional Rietveld and small box modelling demonstrates that Y 3 TaO 7 is best modeled as C222 1 . While the Ccmm describes equally well the long-range structure, pair distribution function analysis revealed that the local atomic configuration can be best modelled as C222 1 . This is corroborated by first-principles calculations that confirm the energetic preference of the C222 1 over Ccmm. Further, neutron total scattering data are reported for the first time for the Y 3 TaO 7 weberite-type ceramics.

36 MATERIALS SCIENCE↗
DOE OSTI2023

LaMnO 3 Dopants for Efficient Thermochemical Water Splitting Identified by Density Functional Theory Calculations

While ceria is the standard material for two-step water splitting, perovskites are emerging as viable alternatives. In this work, based on the orthorhombic LaMnO 3 supercell, we substitute Li Na K Rb Mg Ca Sr Ba on the A-sites (La sites) and Al Ga In Mg Zn on the B-sites (Mn sites) at a concentration of 37.5%. The range of temperature and oxygen partial pressure at which each composition is stable is predicted. For compositions that are stable in relevant temperature and pressure ranges, the oxygen vacancy formation energies are determined for all of the oxygen vacancy site positions available in the computational supercell. Mg, Ca, Sr, and Ba A-site-substituted LaMnO 3 and Al and In B-site-substituted LaMnO 3 meet these two criteria for candidates in solar-thermal water splitting applications. Finally, oxygen vacancy formation energy can also be controlled by adjusting the doping strategy.

08 HYDROGEN↗
DOE OSTI2023

Machine-learning-aided density functional theory calculations of stacking fault energies in steel

A combined large-scale first principles approach with machine learning and materials informatics is proposed to quickly sweep the chemistry-composition space of advanced high strength steels (AHSS). AHSS are composed of iron and key alloying elements such as aluminum and manganese. A systematic exploration of the distribution of aluminum and manganese atoms in iron is used to investigate low stacking fault energies configurations using first principles calculations. To overcome the computational cost of exploring the composition space, this process is sped up using an automated machine learning tool: DeepHyper. Here our results predict that it is energetically favorable for Al to stay away from a stacking fault, but Mn atoms do not affect the stacking fault energy and can stay in the vicinity of the fault. The distribution of Al and Mn atoms in systems containing stacking faults and the effects of their interactions on the equilibrium distribution are systematically analyzed.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Elucidating Biomass-Derived Pyrolytic Lignin Structures from Demethylation Reactions through Density Functional Theory Calculations

Pyrolytic lignin is a fraction of pyrolysis oil that contains a wide range of phenolic compounds that can be used as intermediates to produce fuels and chemicals. However, the characteristics of the raw lignin structure make it difficult to establish a pyrolysis mechanism and determine pyrolytic lignin structures. Herein this study proposes dimer, trimer, and tetramer structures based on their relative thermodynamic stability for a hardwood lignin model in pyrolysis. Different configurations of oligomers were evaluated by varying the positions of the guaiacyl (G) and syringyl (S) units and the bonds βO4 and β5 in the hardwood model lignin through electronic structure calculations. The homolytic cleavage of βO4 bonds is assumed to occur and generate two free radical fragments. These can stabilize by taking hydrogen radicals that may be in solution during the intermediate liquid (pathway 1) formation before the thermal ejection. An alternative pathway (pathway 2) could occur when the radicals use intramolecular hydrogen, turning themselves into stable products. Subsequently, a demethylation reaction can take place, thus generating a methane molecule and new oligomeric lignin-derived molecules. The most probable resulting structures were studied. We used FTIR and NMR spectra of selected model compounds to evaluate our calculation approach. Thermophysical properties were calculated using group contribution methods. The results give insights into the lignin oligomer structures and how these molecules are formed. They also provide helpful information for the design of pyrolysis oil separation and upgrading equipment.

09 BIOMASS FUELS↗
DOE OSTI2024

Sparse-Stochastic Fragmented Exchange for Large-Scale Hybrid Time-Dependent Density Functional Theory Calculations

Here we extend our recently developed sparse-stochastic fragmented exchange formalism for ground-state near-gap hybrid DFT to calculate absorption spectra within linear-response time-dependent generalized Kohn-Sham DFT (LR-GKS-TDDFT) for systems consisting of thousands of valence electrons within a grid-based/plane-wave representation. A mixed deterministic/fragmented-stochastic compression of the exchange kernel, here using long-range explicit exchange functionals, provides an efficient method for accurate optical spectra. Both real-time propagation as well as frequency-resolved Casida-equation-type approaches for spectra are presented, and the method is applied to large molecular dyes.

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

Exploring Anomalous Photoelectron Angular Distributions in the Photoelectron Spectra of Gd 3 O 3 – : Study of Gd 3 O 2 – and Gd 3 O 3 – Using Photoelectron Spectroscopy and Density Functional Theory Calculations

Anion photoelectron (PE) spectra of lanthanide oxide clusters obtained previously have exhibited anomalous photoelectron angular distributions which were attributed to strong PE–valence electron (PEVE) interactions. Here, to further explore this effect, we have obtained the PE spectra of Gd 3 O 2 – and Gd 3 O 3 – , two clusters that have similarly complex electronic structures but contrasting symmetries. The spectra exhibit manifolds of detachment transitions at similar binding energies in a 0.5 eV window of energy. The electron affinity of Gd 3 O 2 is measured to be 1.29 ± 0.05 eV, and that of Gd 3 O 3 is 1.31 ± 0.05 eV. As seen in previous studies on lanthanide oxide cluster anions in lower than conventional oxidation states, transitions in spectra obtained lower photon energies are more congested than those obtained with higher photon energy, a signature of strong PEVE interactions. While the detachment transitions have predominantly parallel photoelectron angular distributions (PAD), the PAD varies across the manifold of transitions in the PE spectrum of Gd 3 O 3 – in a way that suggests four different subgroups of transitions. Results of calculations on Gd 3 O 2 – suggest kite or V-shape structures with antiferromagnetic coupling between one of the 4f 7 subshells with the two others. Calculations on Gd 3 O 3 – more definitively point to ring structures with a nearly isoenergetic ferromagnetically coupled high spin (24-tet) state and a dectet state in which one of the 4f 7 subshells is antiferromagnetically coupled with the other two. Taking these results as qualitative, we propose that strong mixing between the unperturbed states predicted computationally leads to overlapping transitions with different PADs.

anions↗
DOE OSTI2024

Dynamic and relativistic effects on Pt–Pt indirect spin–spin coupling in aqueous solution studied by ab initio molecular dynamics and two- vs four-component density functional NMR calculations

Treating 195Pt nuclear magnetic resonance parameters in solution remains a considerable challenge from a quantum chemistry point of view, requiring a high level of theory that simultaneously takes into account the relativistic effects, the dynamic treatment of the solvent–solute system, and the dynamic electron correlation. A combination of Car-Parrinello molecular dynamics (CPMD) and relativistic calculations based on two-component zeroth order regular approximation spin–orbit Kohn–Sham (2c-ZKS) and four-component Dirac–Kohn–Sham (4c-DKS) Hamiltonians is performed to address the solvent effect (water) on the conformational changes and JPtPt1 coupling. A series of bridged PtIII dinuclear complexes [L1–Pt2(NH3)4(Am)2–L2]n+ (Am = α–pyrrolidonate and pivalamidate; L = H2O, Cl−, and Br−) are studied. The computed Pt–Pt coupling is strongly dependent on the conformational dynamics of the complexes, which, in turn, is correlated with the trans influence among axial ligands and with the angle N–C–O from the bridging ligands. The J-coupling is decomposed in terms of dynamic contributions. The decomposition reveals that the vibrational and explicit solvation contributions reduce JPtPt1 of diaquo complexes (L1 = L2 = H2O) in comparison to the static gas-phase magnitude, whereas the implicit solvation and bulk contributions correspond to an increase in JPtPt1 in dihalo (L1 = L2 = X−) and aquahalo (L1 = H2O; L2 = X−) complexes. Relativistic treatment combined with CPMD shows that the 2c-ZKS Hamiltonian performs as well as 4c-DKS for the JPtPt1 coupling.

Chemistry↗
DOE OSTI2021

High-pressure phase diagram of beryllium from ab initio free-energy calculations

In this report we use first-principles molecular dynamics simulations coupled with the thermodynamic integration method to study the hexagonal close-packed (hcp) to body-centered cubic (bcc) transition and melting of beryllium up to a pressure of 1600 GPa. We derive the melting line by equating solid and liquid Gibbs free energies and represent it by a Simon-Glatzel fit T m = 1564 K [1 + P/(15.6032 GPa)] 0.383 , which is in good agreement with previous two-phase simulations <6000 K. We also derive the hcp-bcc solid-solid phase boundary and show that the quasiharmonic approximation underestimates the stability of the hcp structure, predicting lower transition pressures between hcp and bcc phases. Our results are consistent with the stability regime predicted by the phonon quasiparticle method. We also predict that the hcp-bcc-liquid triple point is located at 164.7 GPa and 4314 K. In addition, we compute the shock Hugoniot curve and show that it is in good agreement with experiments, intersecting our derived melting curve at ~235 GPa and 4900 K. Finally, we make predictions for future ramp compression experiments. Starting with an isentropic compression of the liquid, we predict the path to intersect the melting line at low pressure and temperature, then to continue along the melting line over a large temperature interval of 7000 K as the sample remains in the mixed solid-liquid state before it enters the solid phase.

36 MATERIALS SCIENCE↗