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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Rational Design of Superconducting Metal Hydrides via Chemical Pressure Tuning**

Abstract The high critical superconducting temperatures ( T c s) of metal hydride phases with clathrate‐like hydrogen networks have generated great interest. Herein, we employ the Density Functional Theory‐Chemical Pressure (DFT‐CP) method to explain why certain electropositive elements adopt these structure types, whereas others distort the hydrogenic lattice, thereby decreasing the T c . The progressive opening of the H 24 polyhedra in MH 6 phases is shown to arise from internal pressures exerted by large metal atoms, some of which may favor an even higher hydrogen content that loosens the metal atom coordination environments. The stability of the LaH 10 and LaBH 8 phases is tied to stuffing of their shared hydrogen network with either additional hydrogen or boron atoms. The predictive capabilities of DFT‐CP are finally applied to the Y−X−H system to identify possible ternary additions yielding a superconducting phase stable to low pressures.

Hilleke, Katerina P.↗
DOE OSTI2022

Excitonic Effects in X-ray Absorption Spectra of Fluoride Salts and Their Surfaces

Given their natural abundance and thermodynamic stability, fluoride salts may appear as evolving components of electrochemical interfaces in Li-ion batteries and emergent multivalent ion cells. This is due to the practice of employing electrolytes with fluorine-containing species (salt, solvent, or additives) that electrochemically decompose and deposit on the electrodes. Operando X-ray absorption spectroscopy (XAS) can probe the electrode–electrolyte interface with a single-digit nanometer depth resolution and offers a wealth of insights into the evolution and Coulombic efficiency or degradation of prototype cells, provided that the spectra can be reliably interpreted in terms of local oxidation state, atomic coordination, and electronic structure about the excited atoms. Here we explore fluorine K-edge XAS of mono- (Li, Na, and K) and di-valent (Mg, Ca, and Zn) fluoride salts from a theoretical standpoint and discover a surprising level of detailed electronic structure information about these materials despite the relatively predictable oxidation state and ionicity of the fluoride anion and the metal cation. Utilizing a recently developed many-body approach based on the ΔSCF method, we calculate the XAS using density functional theory and experimental spectral profiles are well reproduced despite some experimental discrepancies in energy alignment within the literature, which we can correct for in our simulations. We outline a general methodology to explain shifts in the main XAS peak energies in terms of a simple exciton model and explain line-shape differences resulting from the mixing of core-excited states with metal d character (for K and Ca specifically). Given ultimate applications to evolving interfaces, some understanding of the role of surfaces and their terminations in defining new spectral features is provided to indicate the sensitivity of such measurements to changes in interfacial chemistry.

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

Iron(III)–Oxo Cluster Chemistry with Dimethylarsinate Ligands: Structures, Magnetic Properties, and Computational Studies

A program has been initiated to develop Fe III /oxo cluster chemistry with the ‘pseudo-carboxylate’ ligand dimethylarsinate (Me 2 AsO 2 - ) for comparison with the well investigated Fe III /oxo/carboxylate cluster area. The synthesis and characterization of three polynuclear Fe III complexes are reported, [Fe 12 O 4 (O 2 C t Bu) 8 (O 2 AsMe 2 ) 17 (H 2 O) 3 ]Cl 3 (1), Na 2 [Fe 12 Na 2 O 4 (O 2 AsMe 2 ) 20 (NO 3 ) 6 (Me 2 AsO 2 H) 2 (H 2 O) 4 ](NO 3 ) 6 (2) and [Fe 3 (O 2 AsMe 2 ) 6 (Me 2 AsO 2 H) 2 (hqn) 2 ](NO 3 ) (3), where hqn is 8-hydroxyquinoline. The Fe 12 core of 1 is a type never previously encountered in Fe III carboxylate chemistry, consisting of two Fe 6 units each of which comprises two {Fe 3 (μ 3 -O 2- )} units bridged by three Me 2 AsO 2 - groups and linked into an Fe 12 loop structure by two anti-anti η 1 :η 1 :μ Me 2 AsO 2 - groups, a bridging mode extremely rare with carboxylates. 2 also consists of two Fe 6 units, differing in their ligation from those in 1, and this time linked together into a linear structure by a central {Na 2 (NO 3 ) 2 } bridging unit. 3 is a linear Fe 3 complex with no monoatomic bridges between Fe III ions, a very rare situation in Fe III chemistry with any ligands, and unprecedented in Fe carboxylate chemistry. The distinct differences observed in arsinate vs carboxylate ligation modes are rationalized largely based on the greater basicity of the former vs the latter. Variable-temperature dc and ac magnetic susceptibility data reveal all Fe 2 pairwise interactions to be antiferromagnetic. For 1 and 2, the different J ij couplings were estimated by use of a magnetostructural correlation for high nuclearity Fe III -oxo clusters and by density functional theory calculations using broken symmetry methods, allowing identification of their relative spin vector alignments and thus rationalization of their S = 0 ground states. The J ij values were then used as input values to give excellent fits of the experimental χM T vs T data. For 3, the fits of the experimental χM T vs T data to the Van Vleck equation or with PHI gave a very weak J 12 = -0.8(1) cm -1 (H = –2JŠ i ·Š j convention) between adjacent Fe III ions, and an S = 5/2 ground state. Furthermore, these initial Fe III arsinate complexes also provide structural parameters that help validate literature assignments of arsinate binding modes to iron oxide/hydroxide minerals as part of environmental concerns of using arsenic-containing herbicides in agriculture.

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

Phosphine Ligand Binding and Catalytic Activity of Group 10–14 Heterobimetallic Complexes

We report heterobimetallic complexes have attracted much interest due to their broad range of structures and reactivities as well as unique catalytic abilities. Additionally, these complexes can be utilized as single-source precursors for the synthesis of binary intermetallic compounds. An example is the family of bis(pyridine-2-thiolato)dichloro-germanium and tin complexes of group 10 metals (Pd and Pt). The reactivity of these heterobimetallic complexes is highly tunable through substitution of the group 14 element and the neutral ligand bound to the transition metal. Here, we study the binding energies of three different phosphorous-based ligands, PR 3 (R = Bu, Ph, and OPh) by density functional theory and restricted Hartree–Fock methods. The PR 3 ligand-binding energies follow the trend of PBu 3 > PPh 3 > P(OPh) 3 , in agreement with their sigma-bonding ability. These results are confirmed by ligand exchange experiments monitored with 31 P NMR spectroscopy, in which a weaker binding PR 3 ligand is replaced with a stronger one. Furthermore, we demonstrate that the heterobimetallic complexes are active catalysts in the Negishi coupling reaction, where stronger binding PR 3 ligands inhibit access to an active site at the metal center. Similar strategies could be applied to other complexes to better understand their ligand-binding energetics and predict their reactivity as both precursors and catalysts.

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

Effect of Linker Structure and Functionalization on Secondary Gas Formation in Metal–Organic Frameworks

Rare-earth terephthalic acid (BDC)-based metal–organic frameworks (MOFs) are promising candidate materials for acid gas separation and adsorption from flue gas streams. However, previous simulations have shown that acid gases (H 2 O, NO 2 , and SO 2 ) react with the hydroxyl on the BDC linkers to form protonated acid gases as a potential degradation mechanism. Herein, gas-phase computational approaches were used to identify the formation energies of these secondary protonated acid gases across multiple BDC linker molecules. Formation energies for secondary protonated acid gases were evaluated using both density functional theory (DFT) and correlated wave function methods for varying BDC–gas reaction mechanisms. Upon validation of DFT to reproduce wave function calculation results, rotated conformational linkers and chemically functionalized BDC linkers with -OH, -NH 2 , and -SH were investigated. Furthermore, the calculations show that the rotational conformation affects the molecule stability. Double-functionalized BDC linkers, where two functional groups are substituted onto BDC, showed varied reaction energies depending on whether the functional groups donate or withdraw electrons from the aromatic system. Based on these results, BDC linker design must balance adsorption performance with degradation via linker dehydrogenation for the design of stable MOFs for acid gas separations.

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

C–C Bond Formation during Electrochemical CO 2 Reduction on Pristine Cu(100) Unlikely to Involve Adsorbed CO at Any Potential

Formation of hydrocarbons containing two or more carbon atoms (C 2+ ) during heterogeneous electrochemical CO and CO 2 reduction (ECOR and ECO 2 R) only occurs, among pure metals, on Cu electrodes. Moreover, the activity and selectivity is facet dependent, with Cu(100) generally preferentially forming ethylene over methane. Previously, we found via quantum-mechanics-based modeling that, unlike standard density functional theory, more accurate correlated wavefunction methods predict that non-electroactive coupling pathways involving two adsorbed COs (*CO) or a *CO and a *COH to form C–C bonds on Cu(100) are kinetically inhibited, with the former also thermodynamically unfavorable. Here, we extend that embedded complete active space second order perturbation theory (ECASPT2) study, further showing that electrochemical coupling of two *COs to form an anionic dimer [OC*–*CO] (1+δ)– , followed by protonation to form [OC*–*COH] δ− , is not kinetically competitive with the reduction of *CO to *COH at relevant ECO/CO 2 R potentials. Our simulations therefore suggest that the ability of Cu(100) to electrochemically synthesize C 2+ molecules from CO and CO 2 is unlikely to be via *CO, at least on pristine Cu(100). Instead, hydrogenated CO species (*COH, *CH x OH, or *CH x ) are most likely to be the key intermediates in C–C bond formation.

Martirez, John Mark P. [Princeton Plasma Physics L↗
DOE OSTI2020

A Drift‐Kinetic Method for Obtaining Gradients in Plasma Properties From Single‐Point Distribution Function Data

Abstract In this paper, we derive a new drift‐kinetic method for estimating gradients in the plasma properties through a velocity space distribution at a single point. The gradients are intrinsically related to agyrotropic features of the distribution function. This method predicts the gradients in the magnetized distribution function and can predict gradients of arbitrary moments of the gyrotropic background distribution function. The method allows for estimates on density and pressure gradients on the scale of a Larmor radius, proving to resolve smaller scales than any method currently available to spacecraft. The model is verified with a set of fully kinetic VPIC particle‐in‐cell simulations.

Wetherton, B. A.↗
DOE OSTI2025

Low-iridium stabilized ruthenium oxide anode catalyst for durable proton-exchange membrane water electrolysis

While mixing iridium (Ir) with ruthenium oxide (RuO 2 ) has proven to be an effective strategy for reducing Ir loading in anode catalysts for proton-exchange membrane (PEM) water electrolysers, achieving industrially relevant long-term stability typically requires an Ir-rich, Ru-lean combination. Here, in this study, by combining density functional theory with Metropolis Monte Carlo methods, we discovered that sufficient stabilization in the RuO 2 lattice could be achieved with less than 50 at.% of Ir, and that Ir in the first subsurface layer plays a critical role. By effectively dispersing Ir dopants within the RuO 2 lattice, we demonstrated an Ir:Ru atomic ratio of only 1:6 that exhibited exceptional stability for over 1,500 h of continuous water electrolysis at 2 A cm −2 . Our Ru 6 IrO x catalyst has the potential to reduce Ir loading by 80% compared with current commercial PEM water electrolysers, and its stability was further validated under industrial testing conditions in a 25-cm 2 PEM electrolyser.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Molecular structure determination: Equilibrium structure of pyrimidine (m-C 4 H 4 N 2 ) from rotational spectroscopy (reSE) and high-level ab initio calculation (re) agree within the uncertainty of experimental measurement

The pure rotational spectrum of pyrimidine (m-C 4 H 4 N 2 ), the meta-substituted dinitrogen analog of benzene, has been studied in the millimeter-wave region from 235 GHz to 360 GHz. The rotational spectrum of the ground vibrational state has been assigned and fit to yield accurate rotational and distortion constants. Over 1700 distinct transitions were identified for the normal isotopologue in its ground vibrational state and least-squares fit to a partial sextic S-reduced Hamiltonian. Transitions for all four singly substituted 13 C and 15 N isotopologues were observed at natural abundance and were likewise fit. Deuterium-enriched samples of pyrimidine were synthesized, giving access to all eleven possible deuterium-substituted isotopologues, ten of which were previously unreported. Experimental values of rotational constants and computed values of vibration–rotation interaction constants and electron-mass corrections were used to determine semi-experimental equilibrium structures (r e SE ) of pyrimidine. The r e SE structure obtained using coupled-cluster with single, double, and perturbative triple excitations [CCSD(T)] corrections shows exceptional agreement with the re structure computed at the CCSD(T)/cc-pCV5Z level (≤0.0002 Å in bond distance and ≤0.03° in bond angle). Of the various computational methods examined, CCSD(T)/cc-pCV5Z is the only method for which the computed value of each geometric parameter lies within the statistical experimental uncertainty (2σ) of the corresponding semi-experimental coordinate. The exceptionally high accuracy and precision of the structure determination is a consequence of the large number of isotopologues measured, the precision and extent of the experimental frequency measurements, and the sophisticated theoretical treatment of the effects of vibration–rotation coupling and electron mass. Taken together, these demanding experimental and computational studies establish the capabilities of modern structural analysis for a prototypical monocyclic aromatic compound.

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

Ab initio calculations in atoms, molecules, and solids, treating spin–orbit coupling and electron interaction on an equal footing

We incorporate explicit, non-perturbative treatment of spin–orbit coupling into ab initio auxiliary-field quantum Monte Carlo (AFQMC) calculations. The approach allows a general computational framework for molecular and bulk systems in which material specificity, electron correlation, and spin–orbit coupling effects can be captured accurately and on an equal footing, with favorable computational scaling vs system size. We adopt relativistic effective-core potentials that have been obtained by fitting to fully relativistic data and that have demonstrated a high degree of reliability and transferability in molecular systems. This results in a two-component spin-coupled Hamiltonian, which is then treated by generalizing the ab initio AFQMC approach. We demonstrate the method by computing the electron affinity in Pb, the bond dissociation energy in Br 2 and I 2 , and solid Bi.

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

Identifying native point defect configurations in α-alumina

Intimately intertwined atomic and electronic structures of point defects govern diffusion-limited corrosion and underpin the operation of optoelectronic devices. For some materials, complex energy landscapes containing metastable defect configurations challenge first-principles modeling efforts. Here, we thoroughly reevaluate native point defect geometries for the illustrative case of α-Al 2 O 3 by comparing three methods for sampling candidate geometries in density functional theory calculations: displacing atoms near a naively placed defect, initializing interstitials at high-symmetry points of a Voronoi decomposition, and Bayesian optimization. We find symmetry-breaking distortions for oxygen vacancies in some charge states, and we identify several distinct oxygen split-interstitial geometries that help explain literature discrepancies involving this defect. We also report a surprising and, to our knowledge, previously unknown trigonal geometry favored by aluminum interstitials in some charge states. Importantly, these new configurations may have transformative impacts on our understanding of defect migration pathways in aluminum-oxide scales protecting metal alloys from corrosion. Overall, the Voronoi scheme appears most effective for sampling candidate interstitial sites because it always succeeded in finding the lowest-energy geometry identified in this study, although no approach found every metastable configuration. Finally, we show that the position of defect levels within the band gap can depend strongly on the defect geometry, underscoring the need to conduct careful searches for ground-state geometries in defect calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

A method to fuse multiphysics waveforms and improve predictive explosion detection: theory, experiment and performance

Natural and human-made sources of transient energy often emit multiple geophysical signatures that include mechanical and electromagnetic waveforms. We present a constructive method to fuse and evaluate statistics that we derive from such multiphysics waveforms that improves our capability to detect small, near-ground explosions over similar methods that consume single signature waveforms. Our method advances Fisher's Combined Probability Test (Fisher's Method) to operate under both hypotheses of a binary test on noisy data and provide researchers with the density functions required to forecast the ability of Fisher's Method to screen fused explosion signatures from noise. We apply this method against 12 d, multisignature explosion and noise records to show (1) that a fused multiphysics waveform statistic that combines radio, acoustic and seismic waveform data can identify explosions roughly 0.8 magnitude units lower than an acoustic emission, STA/LTA detector for the same detection probability and (2) that we can quantitatively predict how this fused, multiphysics statistic performs with Fisher's Method. Our work thereby offers a baseline method for predictive waveform fusion that supports multiphenomenological explosion monitoring (multiPEM) and is applicable to any binary testing problem in observational geophysics.

58 GEOSCIENCES↗
DOE OSTI2024

Diffusion quantum Monte Carlo approach to the polaritonic ground state

Making and using polaritonic states (i.e., hybrid electron-photon states) for chemical applications has recently become one of the most prominent and active fields that connects the communities of chemistry and quantum optics. Modeling of such polaritonic phenomena using ab initio approaches calls for new methodologies, leading to the reinvention of many commonly used electronic structure methods, such as Hartree-Fock, density functional, and coupled cluster theories. Here in this work, we explore the formally exact diffusion quantum Monte Carlo approach to obtain numerical solutions to the polaritonic ground state during the dissociation of the H 2 molecular system. Here, we examine various electron-nuclear-photon properties throughout the dissociation, such as changes to the minimum of the cavity Born-Oppenheimer surface, the localization of the electronic wave function, and the average mode occupation. Finally, we directly compare our results to that obtained with state-of-the-art, yet approximate, polaritonic coupled cluster approaches.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Time-dependent orbital-free density functional theory: Background and Pauli kernel approximations

Time-dependent orbital-free density functional theory (DFT) is an efficient method for calculating the dynamic properties of large-scale quantum systems due to the low computational cost compared to standard time-dependent DFT. In this work, we formalize this method by mapping the real system of interacting fermions onto a fictitious system of noninteracting bosons. The dynamic Pauli potential and associated kernel emerge as key ingredients of time-dependent orbital-free DFT. Using the uniform electron gas as a model system, we derive an approximate frequency-dependent Pauli kernel. Pilot calculations suggest that space nonlocality is a key feature for this kernel. Nonlocal terms arise already in the second-order expansion with respect to unitless frequency and reciprocal space variable ($\frac{ω}{qk_F}$ and $\frac{q}{2k_F}$, respectively). Given the encouraging performance of the proposed kernel, we expect it will lead to more accurate orbital-free DFT simulations of nanoscale systems out of equilibrium. Additionally, the proposed path to formulate nonadiabatic Pauli kernels presents several avenues for further improvements which can be exploited in future work to improve the results.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Nuclear quantum effects on the quasiparticle properties of the chloride anion aqueous solution within the GW approximation

Photoelectron spectroscopy experiments in ionic solutions reveal important electronic structure information, in which the interaction between hydrated ions and water solvent can be inferred. Based on many-body perturbation theory with GW approximation, we theoretically compute the quasiparticle electronic structure of chloride anion solution, which is modeled by path-integral abinitio molecular dynamics simulation by taking account the nuclear quantum effects (NQEs). The electronic levels of hydrated anion as well as water are determined and compared with the recent experimental photoelectron spectra. It is found that NQEs improve the agreement between theoretical prediction and experiment because NQEs effectively weaken the hybridization of the between the Cl – anion and water. Furthermore, our study indicates that NQEs plays a small but non-negligible role in predicting the electronic structure of the aqueous solvation of ions of the Hofmeister series.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Pressure-induced superconductivity in the hydrogen-rich pseudobinary CaB - H n compounds

Here, the crystal structures of CaB – Hn compounds with n = 1 – 12 in a pressure range of 50–300 GPa were studied using the genetic algorithm method and first-principles density-function theory calculations. Stable structures with stoichiometry of CaB H 6 and CaB H 7 were predicted in different pressure range. BH 4 , BH 5 , and BH 6 units were found to be the main motifs in these compounds. Moreover, metastable Imm2 CaB H 7 is dynamically stable above 180 GPa, with the formation of tetrahedral BH 4 unit surrounded by Ca atom and H 3 unit. Electron-phonon coupling (EPC) calculations reveal that the superconducting properties are closely related to the strong hydrogen-boron bonding of the BH 4 unit in Imm2, and T c can reach ~200 K at 200 GPa. As the major units, BH 4 , BH 5 , and BH 6 units exist in the Imm2, P 2 , and P2 1 /m phases of CaB H 7 at 300 GPa, respectively; the corresponding EPC parameter λ decreases with the increase of hydrogen content in CaB H 7 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Electronic and optical properties of crystalline nitrogen versus black phosphorus: A comparative first-principles study

We report crystalline black nitrogen (BN) is an allotrope of nitrogen with the black phosphorus (BP) structure recently synthesized at high pressure by two independent research groups. Here, we present a systematic study of the electronic and optical properties of BN focusing on its comparison with BP. To this end, we use the state-of-the-art quasiparticle self-consistent GW approach with vertex corrections in both the electronic and optical channels. Despite many similarities, the properties of BN are found to be considerably different. Unlike BP, BN exhibits a larger optical gap (2.5 vs 0.26 eV), making BN transparent in the visible spectral region with a highly anisotropic optical response. This difference can be primarily attributed to a considerably reduced dielectric screening in BN, leading to enhancement of the effective Coulomb interaction. Despite relatively strong Coulomb interaction, exciton formation is largely suppressed in both materials. Our analysis of the elastic properties shows exceptionally high stiffness of BN, comparable to that of diamond.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Double and Charge-Transfer Excitations in Time-Dependent Density Functional Theory

Time-dependent density functional theory has emerged as a method of choice for calculations of spectra and response properties in physics, chemistry, and biology, with its system-size scaling enabling computations on systems much larger than otherwise possible. While increasingly complex and interesting systems have been successfully tackled with relatively simple functional approximations, there has also been increasing awareness that these functionals tend to fail for certain classes of approximations. Here I review the fundamental challenges the approximate functionals have in describing double excitations and charge-transfer excitations, which are two of the most common impediments for the theory to be applied in a black-box way. At the same time, I describe the progress made in recent decades in developing functional approximations that give useful predictions for these excitations.

Chemistry↗