Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DENSITY FUNCTIONAL METHOD”

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

Calculated Reduction Potentials of Electrolyte Species in Lithium–Sulfur Batteries

In this work, reduction potentials of electrolyte molecules in the lithium–sulfur (Li/S) battery and their variations in several solvent environments are studied using the density functional theory method with Dunning’s triple-ζ correlation consistent basis set. Reliable reduction potential values are key for electrolyte additive design needed for suppressing polysulfide dissolution shuttle mechanism, resulting in poor cycle performance and severe self-discharge of the Li/S battery. Although isolated electrolyte molecules have reduction potentials outside the operating voltage range of the Li/S battery, complexation with other electrolyte species enables the electrolyte molecules to be reduced within the operating voltage range. Among the electrolyte species considered in this study, bis(fluorosulfonyl)imide (FSI – ) and fluoroethylene carbonate (FEC) yield reduction potentials within the expected range, suggesting the development of fluorine-containing additives as a promising line of research.

25 ENERGY STORAGE↗
DOE OSTI2023

Ethylene Hydrogenation Molecular Mechanism on MoC y Nanoparticles

Ethylene hydrogenation catalyzed by MoC y nanoparticles has been studied by means of density functional theory methods and several models. These include MetCar (Mo 8 C 12 ), Nanocube (Mo 14 C 13 ), and Mo 12 C 12 nanoparticles as representatives of experimental MoC y nanostructures. The effect of hydrogen coverage has been studied in detail by considering low-, intermediate-, and high-hydrogen regimes. The calculated enthalpy and energy barriers show that ethylene hydrogenation is feasible on the MetCar, Mo 12 C 12 , and Nanocube but at low, medium, and high hydrogen coverages, respectively. An additional step, related to the H* migration from a Mo to a C site in the nanoparticle, has been found to be the key to establishing the best hydrogenation system. In most cases, the reactions are exothermic, featuring low hydrogenation energy barriers, especially for the Nanocube at high hydrogen coverage. In addition, the calculated adsorption Gibbs free energy shows that, for this system, the C 2 H 4 adsorption is feasible in the 300–400 K temperature range and pressures from 10 –10 to 2 atm. For the hydrogenation steps, calculated transition state theory rates show that the overall process is limited by the first hydrogenation step (C 2 H 4 → C 2 H 5 ) at temperatures of 330–400 K. However, at the lower temperatures of 300–320 K, the reaction rates are comparable for the two steps. The present results indicate that the Mo 14 C 13 Nanocube models of MoC y nanoparticles exhibit appropriate thermodynamic and kinetic features to catalyze ethylene hydrogenation at a high-hydrogen-coverage regime. The present findings provide a basis for understanding the chemistry of active MoC y catalysts, suggest appropriate working conditions for the reaction to proceed, and provide a basis for future experimental studies.

03 NATURAL GAS↗
DOE OSTI2021

Optimizing accuracy and efficacy in data-driven materials discovery for the solar production of hydrogen

The production of hydrogen fuels, via water splitting, is of practical relevance for meeting global energy needs and mitigating the environmental consequences of fossil-fuel-based transportation. Water photoelectrolysis has been proposed as a viable approach for generating hydrogen, provided that stable and inexpensive photocatalysts with conversion efficiencies over 10% can be discovered, synthesized at scale, and successfully deployed. While a number of first-principles studies have focused on the data-driven discovery of photocatalysts, in the absence of systematic experimental validation, the success rate of these predictions may be limited. We address this problem by developing a screening procedure with co-validation between experiment and theory to expedite the synthesis, characterization, and testing of the computationally predicted, most desirable materials. Here, starting with 70 150 compounds in the Materials Project database, the proposed protocol yielded 71 candidate photocatalysts, 11 of which were synthesized as single-phase materials. Experiments confirmed hydrogen generation and favorable band alignment for 6 of the 11 compounds, with the most promising ones belonging to the families of alkali and alkaline-earth indates and orthoplumbates. This study shows the accuracy of a nonempirical, Hubbard-corrected density-functional theory method to predict band gaps and band offsets at a fraction of the computational cost of hybrid functionals, and outlines an effective strategy to identify photocatalysts for solar hydrogen generation.

08 HYDROGEN↗
DOE OSTI2021

Correction: Enhanced 5f-δ bonding in [U(C 7 H 7 ) 2 ] - : C K-edge XAS, magnetism, and ab initio calculations

The authors regret that two references to the multi-configurational pair-density functional theory methods used in the wavefunction calculations did not appear in the main article. The references were originally included in the Methods section of the Electronic Supplementary Information as ref. 17a and b, but should also have appeared in the second paragraph on page 9564. The references are listed below as ref. 1 and 2. The Royal Society of Chemistry apologises for these errors and any consequent inconvenience to authors and readers.

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

Catlas: an automated framework for catalyst discovery demonstrated for direct syngas conversion

Catalyst discovery is paramount to support access to energy and key chemical feedstocks in a post fossil fuel era. Exhaustive computational searches of large material design spaces using ab initio methods like density functional theory (DFT) are infeasible. We seek to explore large design spaces at relatively low computational cost by leveraging large, generalized, graph-based machine learning (ML) models, which are pretrained and therefore require no upfront data collection or training. We present Catlas, a framework that distributes and automates the generation of adsorbate-surface configurations and ML inference of DFT energies to achieve this goal. Catlas is open source, making ML assisted catalyst screenings easy and available to all. To demonstrate its efficacy, we use Catlas to explore catalyst candidates for the direct conversion of syngas to multi-carbon oxygenates. For this case study, we explore 947 stable/metastable binary, transition metal intermetallics as possible catalyst candidates. On this subset of materials, we are able to predict the adsorption energy of key descriptors, *CO and *OH, with good accuracy (0.16, 0.14 eV MAE, respectively). Using the projected selectivity towards C 2+ oxygenates from an existing microkinetic model, we identified 144 candidate materials. For 10 promising candidates, DFT calculations reveal a good correlation with our assessment using ML. Finally, among the top elemental combinations were Pt–Ti, Pd–V, Ni–Nb, and Ti–Zn, all of which appear unexplored experimentally.

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

Tricyclononenes and tricyclononadienes as efficient monomers for controlled ROMP: understanding structure–propagation rate relationships and enabling facile post-polymerization modification

Grubbs 3rd-generation (G3) pre-catalyst-initiated ring-opening metathesis polymerization (ROMP) remains an indispensable tool in the polymer chemist's toolbox. Tricyclononenes (TCN) and tricyclononadienes (TCND) represent under-explored classes of monomers for ROMP that have the potential to both advance fundamental knowledge (e.g., structure-polymerization kinetics relationships) and serve as practical tools for the polymer chemist (e.g., post-polymerization functionalization). In this work, a library of TCN and TCND imides, monoesters, and diesters, along with their exo-norbornene counterparts, were synthesized to compare their behaviors in G3-initiated ROMP. Real-time 1 H NMR was used to study their polymerization kinetics; propagation rates ( k p ) were extracted for each monomer. To understand the relationships between monomer structure and ROMP propagation rates, density functional theory methods were used to calculate a variety of electronic and steric parameters for each monomer. While electronic parameters (e.g., HOMO energy levels) correlated positively with the measured k p values, steric parameters generally gave improved correlations, which indicates that monomer size and shape are better predictors for k p than electronic parameters for this data set. Furthermore, the TCND diester—which contains an electron-deficient cyclobutene that is resistant to ROMP—and its polymer p(TCND) are shown to be highly reactive toward DBU-catalyzed conjugate addition reactions with thiols, providing a protecting- and activating-group free strategy for post-polymerization modification.

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

Excited states in variational Monte Carlo using a penalty method

Here, in this article, the authors present a technique using variational Monte Carlo to solve for excited states of electronic systems. This technique is based on enforcing orthogonality to lower energy states, which results in a simple variational principle for the excited states. Energy optimization is then used to solve for the excited states. This technique is applied to the well-characterized benzene molecule, in which ~10 000 parameters are optimized for the first 12 excited states. Agreement within ~0.2 eV is obtained with higher scaling coupled cluster methods; small disagreements with experiment are likely due to vibrational effects.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Atomic isotropic hyperfine properties for first row elements (B–F) revisited

Benchmark quality isotropic hyperfine properties have been obtained for first row elements (B–F) using a systematic composite approach consisting of a sequence of core/valence correlation consistent basis sets, up through aug-cc-pCV8Z, along with configuration interaction and coupled cluster theory methods. The best nonrelativistic final values (in MHz) are 10.64 (B), 20.22 (C), 10.59 (N), –31.74 (O), and 318.30 (F) and are in very good agreement with available experimental values for these difficult-to-describe properties. Agreement is especially close in the case of N, which has the most accurate experimental value. The spin densities derived from the best composite level of theory were found to closely follow a simple quadratic scaling with the atomic number, Z. Finally, observed convergence rates in the 1-particle and n-particle expansions obtained here may be useful in judging likely accuracy that can be expected in studies of molecular systems.

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

Atomic isotropic hyperfine properties for second row elements (Al–Cl)

Isotropic hyperfine properties have been obtained for the second row elements Al–Cl using a systematic composite approach consisting of a sequence of core/valence correlation consistent basis sets, up through aug-cc-pCV7Z, along with configuration interaction and coupled cluster methods. The best nonrelativistic final values for the atomic ground states (in MHz) are -1.80 27 Al ( 2 P o 1/2 ), -24.31 29 Si ( 3 P 0 ), 63.70 31 P ( 4 So 3/2 ), 20.77 33 S ( 3 P 2 ), and 35.42 35 Cl ( 2 P o 3/2 ). We find a large K shell contribution to the spin density at the nucleus that is almost canceled by the L and M shell contributions. The spin density in atomic units is approximately linear with respect to the atomic number.

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

Examination of probability distribution of mixture fraction in LES/FDF modelling of a turbulent partially premixed jet flame

An accurate prediction of the probability density function (PDF) of the mixture fraction is crucial to the prediction of combustion since mixing plays an important role in turbulent non-premixed and partially premixed flames. This work provides an assessment of the large-eddy simulation (LES)/filtered density function (FDF) method for the prediction of the PDF of the mixture fraction. The advantage of the LES/FDF method is that it provides the full predictions of the statistical distribution of scalars including the mixture fraction. The predictive accuracy of the method for the PDF is yet to be fully validated. Assessing the prediction of the PDF of the mixture fraction, a conserved scalar, is an important starting point. The Sydney/Sandia inhomogeneous inlet jet flame is used as a test case. A quick comparison shows that the LES/FDF predicted PDF shapes of the mixture fraction deviate significantly from the commonly presumed Beta-PDF as well as from the experimental data in the flame. Here, to examine the source of the discrepancy, we clarify the different PDF definitions used in the comparison among the predictions, measurements, and the presumed shape PDFs. The discrepancy observed from the comparison is largely reconciled by clarifying the difference between the PDFs that are examined. The PDF of the resolved mixture fraction is shown to be close to the Beta-PDF in both the measurements and predictions, while the PDF directly deduced from the LES/FDF particles deviates significantly from the Beta-PDF. A multimodal PDF analysis and a pseudo convergence analysis are conducted to provide plausible evidence to support the predicted multimodal PDF shapes. The sub-filter scale FDF is shown to be close to the Beta-PDF too through the construction of a synthesized PDF, which supports the common presumed Beta-PDF assumption used in the presumed PDF methods when combined with LES.

42 ENGINEERING↗
DOE OSTI2025

First principles study of the Fermi surface topology of CeCu 2 ⁢Si 2

Since the discovery of heavy-fermion superconductivity in CeCu 2 ⁢Si 2 , the material has attracted great interest, particularly with regard to the nature of the superconducting pairing and its mechanism. Consequently, it is essential to better understand the electronic Fermi surface topology and its role in strong antiferromagnetic fluctuations. The standard density functional theory method is insufficient to model the interplay of strong on-site Coulomb repulsion in localized 4⁢𝑓 electrons and their hybridization with itinerant ligand-orbital electrons. We have performed electronic ground-state calculations on CeCu 2 ⁢Si 2 using the Gutzwiller wave function approximation. The Gutzwiller approximation captures the quasiparticle band renormalization from the strong on-site Coulomb repulsion. We have performed an analysis of this effect on the electronic structure and the Fermi surface topology by varying the interaction strength and taking into account the crystal-field splitting. Using the de Haas-van Alphen effect, the extremal Fermi surface cross-sectional areas were calculated to quantify the effects of quasiparticle mass renormalization on the Fermi surface. Our results confirm the presence of two Fermi surface sheets corresponding to the heavy (488⁢𝑚 𝑒 ) and light (4.35⁢𝑚 𝑒 ) quasiparticles when the crystal-field splitting is accounted for on equal footing with the electronic correlations. This method gives the best agreement with experimental measurements as well as the renormalized band method.

36 MATERIALS SCIENCE↗
DOE OSTI2024

First-principles residual resistivity using a locally self-consistent multiple scattering method

The locally self-consistent multiple scattering (LSMS) method can perform efficient first-principles calculations of systems with a large number of atoms. Here, in this paper, we combine the Kubo-Greenwood equation with LSMS, enabling us to calculate the first-principles residual resistivity of large systems. This has been implemented in the open-source code LSMS . We apply this method to selected pure elements and binary random alloys. The results compare well with experiment, and with values obtained from a first-principles effective medium technique (the Korringa-Kohn-Rostoker coherent potential approximation). We discuss future applications of this method to complex systems where other methods are not applicable.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Tuning superconductivity and spin-vortex instabilities in CaKFe 4 As 4 through in-plane antisymmetric strains

Lattice strains of appropriate symmetry have served as an excellent tool to explore the interaction of superconductivity in the iron-based superconductors with orthorhombic-nematic and stripe spin-density-wave (SSDW) order. In this Letter, we contribute to a broader understanding of the coupling of strain to superconductivity and competing normal-state orders by studying CaKFe 4 As 4 under large, in-plane strains of B 1 g and B 2 g symmetry. In contrast to the majority of iron-based superconductors, pure CaKFe 4 As 4 exhibits superconductivity with a relatively high transition temperature of T c ∼ 35 K in proximity of a noncollinear, tetragonal, hedgehog spin-vortex crystal (SVC) order. Through experiments and calculations, we demonstrate an anisotropic in-plane strain response of T c and the favored SVC configuration in CaKFe 4 As 4 . This supports a scenario, in which the change in spin fluctuations dominates the strain response of superconducting T c . Overall, by suggesting moderate B 2 g strains as an effective parameter to change the stability of SVC and SSDW, we outline a pathway to a unified phase diagram of iron-based superconductivity. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Coupling between collective modes in the deformed 98 Zr nucleus: Insights from consistent HFB + QRPA calculations with the Gogny interaction

The zirconium isotopes exhibit structural properties that present multiple challenges to nuclear theory. Investigations of the coupling present within isoscalar modes and within isovector modes are scarce but important for advancing our understanding of the microscopic picture of nuclei. To explore some of these underlying coupling features, and to test the predictive power of a state-of-the-art nuclear structure approach, we provide a detailed analysis of the properties of 90,96,98 Zr . This region includes a benchmarking case and offers insights into nuclear deformation phenomena. Here, to investigate the coupling between collective modes in deformed nuclei, we focused our analysis on the ground and excited-state properties of these isotopes, employing a consistent approach with the axially symmetric deformed Hartree-Fock-Bogoliubov (HFB) and the quasiparticle random phase approximation (QRPA) framework, both using the Gogny D1M force. This approach effectively describes both low-lying and giant-resonance states. We devoted special attention to the deformed 98 Zr nucleus, where we confirm the existence of coupling between monopole and quadrupole excitations through the 𝐾 𝜋 = 0 + QRPA components and demonstrate an analogous dipole-octupole coupling through the 𝐾 𝜋 = 0 − and 𝐾 𝜋 = 1 − components. Intrinsic transition densities and associated radial projections illustrate the coupling. Our work complements and extends earlier studies carried out using density-functional-based methods and notably, we included the complete Coulomb interaction also in the pairing fields, i.e., we treat terms exactly that are approximated in typical calculations that use the Gogny D1 and D2 interaction families.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Possible Kitaev Quantum Spin Liquid State in 2D Materials with S=3/2

Quantum spin liquids (QSLs) form an extremely unusual magnetic state in which the spins are highly correlated and fluctuate coherently down to the lowest temperatures, but without symmetry breaking and without the formation of any static long-range-ordered magnetism. Such intriguing phenomena are not only of great fundamental relevance in themselves, but also hold promise for quantum computing and quantum information. Further, among different types of QSLs, the exactly solvable Kitaev model is attracting much attention, with most proposed candidate materials, e.g., RuCl 3 and Na 2 IrO 3 , having an effective $\textit{S}$ = 1/2 spin value. Here, via extensive first-principles-based simulations, we report the investigation of the Kitaev physics and possible Kitaev QSL state in epitaxially strained Cr-based monolayers, such as CrSiTe 3 , that rather possess a $\textit{S}$ = 3/2 spin value. Our study thus extends the playground of Kitaev physics and QSLs to 3$\textit{d}$ transition metal compounds.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Structure of the (110) Ln ScO 3 ( Ln =Gd,Tb,Dy) surfaces

The surface reconstruction on the (110) surface of LnScO 3 (Ln=Gd,Tb,Dy) has been investigated using x-ray photoelectron spectroscopy, atomic force microscopy, and transmission electron diffraction coupled with higher-level density-functional theory methods. The experimental techniques generate constraints on the surface chemistry and structure that are used to inform the theoretical comparison of a number of potential surface structures. The resulting structure is then compared to experimental aberration-corrected transmission electron microscopy results using multislice simulation. The surfaces of both single-crystalline substrates and hydrosauna synthesized nanoparticles exhibit a Sc-rich double layer with a Sc 3 O 4 termination.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Nature of native atomic defects in ZrTe 5 and their impact on the low-energy electronic structure

Over the past decades, investigations of the anomalous low-energy electronic properties of ZrTe 5 have reached a wide array of conclusions. An open question is the growth method's impact on the stoichiometry of ZrTe 5 samples, especially given the very small density of states near its chemical potential. In this work, we report on high-resolution scanning tunneling microscopy and spectroscopy measurements performed on samples grown via different methods. Using density functional theory calculations, we identify the most prevalent types of atomic defects on the surface of ZrTe 5 , namely, Te vacancies and intercalated Zr atoms. Finally, we precisely quantify their density and outline their role as ionized defects in the anomalous resistivity of this material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Insulating band gaps both below and above the Néel temperature in d -electron LaTi O 3 , LaV O 3 , SrMn O 3 , and LaMn O 3 perovskites as a symmetry-breaking event

Metal d-electron oxides having an odd number of electrons per cell should exhibit band degeneracy at the Fermi energy, making them, in band theory, formally metallic. In many cases, however, these are “false metals”, as evidenced by the observation that many ABO 3 oxide perovskites with a magnetic 3d B-atom are observed to be insulators both below and above the Néel temperature. These inconsistencies between experimental observation and expectation have historically been resolved by invoking degeneracy-breaking physics, based largely on -pure electron effects, such as strong interelectronic correlation for d-electron compounds (the Mott mechanism). Such explanations generally consider microscopic lattice or magnetic degrees of freedom (m-DOF) as largely passive spectators, not a cause of the formally metallic being an insulator. Yet, it has long been known that ABO 3 perovskites can manifest an arrangement of m-DOFs in the form of octahedral tilting, bond dimerization, Jahn-Teller distortions, and ordering of local magnetic moments. It appears reasonable that such structural and magnetic local degrees of freedom need to be allowed to compete with purely electronic strong correlation. To answer this question, we explored a range of d-electron oxide perovskites exemplified by the archetypes LaTiO 3 , LaVO 3 , SrMnO 3 , and LaMnO 3 with 1, 2, 3, or 4 d -electrons, respectively. Using a mean-field-like electronic structure method (here, density functional theory), we find that a combination of magnetic symmetry breaking with structural distortions can account for insulating band gaps in this series while at the same time correctly predicting for the control case an intrinsic paramagnetic metal in SrVO 3 as symmetry-breaking is insufficiently strong to remove the degeneracy. Finally, this indicates that calculating quantitatively local magnetic and positional symmetry-breaking motifs in unit cells that avoid averaging at the outset over the low symmetry motifs can provide consistent trends in a Mott transition without Mott U.

36 MATERIALS SCIENCE↗