Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional theory development”

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

Form factors and spectral densities from Lightcone Conformal Truncation

We use the method of Lightcone Conformal Truncation (LCT) to obtain form factors and spectral densities of local operators $\mathcal{O}$ in $\phi^4$ theory in two dimensions. We show how to use the Hamiltonian eigenstates from LCT to obtain form factors that are matrix elements of a local operator $\mathcal{O}$ between single-particle bra and ket states, and we develop methods that significantly reduce errors resulting from the finite truncation of the Hilbert space. We extrapolate these form factors as a function of momentum to the regime where, by crossing symmetry, they are form factors of $\mathcal{O}$ between the vacuum and a two-particle asymptotic scattering state. We also compute the momentum-space time-ordered two-point functions of local operators in LCT. These converge quickly at momenta away from branch cuts, allowing us to indirectly obtain the time-ordered correlator and the spectral density at the branch cuts. We focus on the case where the local operator $\mathcal{O}$ is the trace Θ of the stress tensor.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Strong Magnetocrystalline Anisotropy Arising from Metal–Ligand Covalency in a Metal–Organic Candidate for 2D Magnetic Order

Layered metal–organic frameworks are promising candidates for new two-dimensional magnets, as their synthetic programmability of these materials can provide a route to diverse structural and electronic properties. However, such framework materials typically lack the heavy elements that engender magnetocrystalline anisotropy in the monolayer ferromagnets reported to date. Alternative sources of magnetic anisotropy are therefore needed in these materials. Here, we report the synthesis of single crystals of the framework material (NMe 4 ) 2 [Fe 2 L 3 ] (H 2 L = 3,6-dichloro-2,5-dihydroxybenzoquinone) and evaluate the angular dependence of its magnetic properties. Oriented-crystal magnetization measurements reveal strong uniaxial anisotropy, where the easy axis is aligned with the crystallographic c axis. While the spin carriers of this structure are isotropic $S = ^5/_2$ Fe III metal centers and $S = ^1/_2$ organic linkers, the anisotropy energy of the framework material is comparable to that of reported 2D ferromagnets. Density functional theory calculations indicate that the observed magnetocrystalline anisotropy arises from ligand-to-metal charge transfer that enhances the magnetic anisotropy of the otherwise isotropic Fe centers, suggesting that metal–ligand covalency can be utilized as a general additive for the development of 2D magnets. These results in this study show the possibility for (NMe 4 ) 2 [Fe 2 L 3 ] to retain magnetic order down to the 2D monolayer limit. In addition, the combination of large magnetic anisotropy and semiconducting character in (NMe 4 ) 2 [Fe 2 L 3 ] highlights its potential as a new 2D magnetic semiconductor.

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

Integrating Machine Learning Potential and X-ray Absorption Spectroscopy for Predicting the Chemical Speciation of Disordered Carbon Nitrides

Precise determination of atomic structural information in functional materials holds transformative potential and broad implications for emerging technologies. Spectroscopic techniques, such as X-ray absorption near-edge structure (XANES), have been widely used for material characterization; however, extracting chemical information from experimental probes remains a significant challenge, particularly for disordered materials. We present an integrated approach that combines atomic simulations, data-driven techniques, and experimental measurements to investigate chemical speciation of amorphous carbon nitride systems as a case study. Here, we discuss the development of machine learning potentials that can efficiently explore the vast configuration space of amorphous carbon nitrides. By employing statistical methods, this structural database enables the elucidation of the most representative local structures and how they evolve with chemical compositions and density. Density functional theory simulations are used to establish a correlation between the local structure and spectroscopic signatures, which then serve as the basis for interpreting and extracting chemical content from experimental data. Although our framework is specifically demonstrated for XANES and carbon nitrides, the approach described herein is readily adaptable as applied to other experimental characterization probes and materials classes.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Nuclear–Electronic Orbital QM/MM Approach: Geometry Optimizations and Molecular Dynamics

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods allow simulations of chemical reactions in atomistic solvent and heterogeneous environments such as proteins. Herein, the nuclear–electronic orbital (NEO) QM/MM approach is introduced to enable the quantization of specified nuclei, typically protons, in the QM region using a method such as NEO-density functional theory (NEO-DFT). This approach includes proton delocalization, polarization, anharmonicity, and zero-point energy in geometry optimizations and dynamics. Expressions for the energies and analytical gradients associated with the NEO-QM/MM method, as well as the previously developed polarizable continuum model (NEO-PCM), are provided. Geometry optimizations of small organic molecules hydrogen bonded to water in either dielectric continuum solvent or explicit atomistic solvent illustrate that aqueous solvation can strengthen hydrogen-bonding interactions for the systems studied, as indicated by shorter intermolecular distances at the hydrogen-bond interface. We then performed a real-time direct dynamics simulation of a phenol molecule in explicit water using the NEO-QM/MM method. Furthermore, these developments and initial examples provide the foundation for future studies of nuclear–electronic quantum dynamics in complex chemical and biological environments.

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

UiO-66 Metal–Organic Framework as an Anode for a Potassium-Ion Battery: Quantum Mechanical Analysis

The natural abundance of potassium in the earth’s crust is 1000 times higher than that of lithium, so energy technologies built on potassium are more sustainable. Potassium-ion batteries have attracted considerable attention because of their relatively low cost and high operating potential, but questions remain about the best anode material for such batteries. Here, in this paper, we report first-principles computations based on density functional theory to investigate the performance of the UiO-66 metal–organic framework as an anode material for potassium-ion batteries; the goal is to provide a fundamental understanding of metal–organic framework (MOF)-based electrodes to guide the design and development of high-performance potassium-ion batteries. Our study includes the stability and electronic properties of potassiated structures and the mechanisms of potassium intercalation and diffusion in the framework. The results indicate that UiO-66 has a maximum specific capacity of 644 mAh/g as the anode of a potassium-ion battery. During potassiation, we observe charge transfer from potassium to carbon or oxygen of UiO-66 near the intercalated K. During K diffusion, the K migrates along the UiO-66 framework with a maximal migration energy barrier of 0.377 eV in the optimal pathway, which is much lower than the barriers for Li and Na diffusion in UiO-66. The diffusion coefficient of K in the anode is several orders of magnitude larger than those of Li and Na. This favors potassium ions over lithium ions or sodium ions when UiO-66 is the anode.

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

From Micro-environments to Macroscopic Effects: How the Alkaline Hydrogen Evolution Reaction Drives Cu Cathodic Corrosion

Cathodic corrosion of copper (Cu) has posed a significant challenge for over a century, impeding various technological progresses such as electrochemical conversion of CO 2 (eCO 2 RR) into fuels and other value-added carbon products. Here, in this study, employing a combined Density Functional Theory (DFT) and kinetic Monte Carlo (kMC) simulation approach, we delve into the atomistic level mechanism driving this phenomenon in Cu. Our hypothesis posits the pivotal role of alkaline hydrogen evolution reaction (HER) in facilitating cathodic corrosion in Cu. We rigorously develop a pH-dependent hydroxide (OH) adsorption mechanism and calculate the equilibrium OH coverage (𝜃 OH ) at varying pH levels, the thermodynamic stability of subsurface oxygen (O sub ), as well as the Cu-vacancy mediated diffusion of subsurface oxygens (O sub ). Through comprehensive analysis, we establish correlation among various microenvironments, including oxygen diffusion in subsurface layers, pH-dependent OH adsorption, and Cu dissolution into the electrolyte as (Cu-OH) complexes. Furthermore, our investigation explores the correlation between surface coordination environment of active sites and cathodic corrosion of Cu. Finally, by integrating DFT-derived thermodynamic data into a kMC model, we successfully predict the formation of experimentally observed corrosion pits on Cu-surfaces. This combined approach not only advances our fundamental understanding of Cu cathodic corrosion but also offers insights crucial for developing effective corrosion mitigation strategies.

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

Machine-learned impurity level prediction for semiconductors: the example of Cd-based chalcogenides

The ability to predict the likelihood of impurity incorporation and their electronic energy levels in semiconductors is crucial for controlling its conductivity, and thus the semiconductor's performance in solar cells, photodiodes, and optoelectronics. The difficulty and expense of experimental and computational determination of impurity levels makes a data-driven machine learning approach appropriate. In this work, we show that a density functional theory-generated dataset of impurities in Cd-based chalcogenides CdTe, CdSe, and CdS can lead to accurate and generalizable predictive models of defect properties. By converting any semiconductor + impurity system into a set of numerical descriptors, regression models are developed for the impurity formation enthalpy and charge transition levels. These regression models can subsequently predict impurity properties in mixed anion CdX compounds (where X is a combination of Te, Se and S) fairly accurately, proving that although trained only on the end points, they are applicable to intermediate compositions. We make machine-learned predictions of the Fermi-level-dependent formation energies of hundreds of possible impurities in 5 chalcogenide compounds, and we suggest a list of impurities which can shift the equilibrium Fermi level in the semiconductor as determined by the dominant intrinsic defects. Machine learning predictions for the dominating impurities compare well with DFT predictions, revealing the power of machine-learned models in the quick screening of impurities likely to affect the optoelectronic behavior of semiconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Modeling dark photon oscillations in our inhomogeneous Universe

Not provided.A dark photon may kinetically mix with the Standard Model photon, leading to observable cosmological signatures. The mixing is resonantly enhanced when the dark photon mass matches the primordial plasma frequency, which depends sensitively on the underlying spatial distribution of electrons. Crucially, inhomogeneities in this distribution can have a significant impact on the nature of resonant conversions. We develop and describe, for the first time, a general analytic formalism to treat resonant oscillations in the presence of inhomogeneities. Our formalism follows from the theory of level crossings of random fields and only requires knowledge of the one-point probability density function (PDF) of the underlying electron number density fluctuations. We validate our formalism using simulations and illustrate the photon-to-dark photon conversion probability for several different choices of PDFs that are used to characterize the low-redshift Universe.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Influence of Ga doping on magnetic properties, magnetocaloric effect, and electronic structure of pseudobinary GdZn 1 - x Ga x ( x = 0 - 0.1 )

Here, we explore the impact of introducing IIIA-group element Ga in place of IIB-group element Zn in binary intermetallic GdZn on its magnetic and magnetocaloric properties, as well as explicate the modified electronic band structure of the compound. The magnetic transition temperature of the compound decreases with the increase of Ga concentration in GdZn 1-x ⁢Ga x ⁢(x = 0–0.1) while the crystal structure (CsCl prototype) and lattice parameters remain unchanged. Our detailed analysis of magnetization and magnetocaloric data conclusively proves that long-ranged magnetic ordering exists in the sample, despite the magnetic interaction considerably weakening with the increase of Ga. The experimental data are rationalized using both a theoretical machine-learning model and first-principle density functional theory. The electronic band structure of GdZn is manifested with some unusual complex features which gradually diminish with Ga doping and the conventional sinusoidal feature of Ruderman-Kittel-Kasuya-Yosida-type interactions also disappears. A mean-field theory model is developed and can successfully describe the overall magnetocaloric behavior of the GdZn 1-x⁢ Ga x series of samples.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Green synthesis and characterization of nontoxic L-methionine capped silver and gold nanoparticles

The simple green method for synthesis of stable L-Methionine (L-Met) capped silver (Ag@LM NPs) and gold (Au@LM NPs) nanoparticles (NPs) without adding any additional reduction agent or stabilizer was developed. Colloidal dispersions were characterized by UV–Vis spectrophotometry. The size and spherical shape of NPs were evaluated by transmission electron microscopy. Their surface covering was confirmed by atomic force microscopy, Fourier transform infrared spectroscopy, dynamic light scattering, and zeta potential measurements. Density functional theory calculations pointed that the preferential adsorption mode of L-Met on both Ag and Au surfaces was a vertical binding geometry via –NH2 group, while horizontal binding mode via [sbnd]S[sbnd] and –NH2 groups is also possible. The genotoxicity (evaluated by the micronucleus assay) of NPs, as well as their effects on some oxidative stress parameters (catalase activity, malondialdehyde level), were assessed in vitro using human peripheral blood cells as a model system. The influence of NPs on the morphology of lymphocyte cells studied using atomic force microscopy revealed that the membrane of cells remained unaffected after the treatment with NPs. When considering the effects of NPs on catalase activity and malondialdehyde level, neither particle type promoted oxidative stress. However, the treatment of lymphocytes with Ag@LM NPs induced a concentration-dependent enhancement of the micronuclei incidence and suppression of the cell proliferation while Au@LM NPs promoted cell proliferation, with no significant effects on micronuclei formation. The Ag@LM NPs were more prone to induce DNA damage than Au@LM NPs, which makes the latter type more suitable for further studies in nano-medicine. © 2019.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Prediction of Redox Potentials for the Late Actinides Cm to Lr Using Electronic Structure Methods

Our previously developed computational method for calculating the aqueous redox potentials of the early actinides has been extended to the later elements in the actinide series: Cm, Bk, Cf, Es, Fm, Md, No, and Lr in multiple oxidation states. These calculations were performed using density functional theory with small-core pseudopotentials and their associated basis sets. Solvation effects were considered via a supermolecule-continuum approach, with 30 water molecules representing two solvation shells. Both the COSMO and SMD implicit solvation models were utilized. The structural parameters and hydration numbers for Cm(III), Bk(III), Bk(IV), and Cf(III) are in reasonable agreement with the available experimental data. For redox processes involving atomic cations in solution, the B3LYP/COSMO approach predicted redox potentials to within ±0.2 V of experiment for most redox couples, consistent with our prior work. Inclusion of spin-orbit corrections in specific redox pairs, especially those with the later actinides in high oxidation states, yields improved results relative to calculations including only scalar-relativistic corrections. The An +m /An(0) redox potentials were calculated using a Born-Haber cycle incorporating sublimation, ionization, and hydration energies. Due to a lack of experimental data, three sets of ionization energies were used for the Born-Haber cycle. The calculated An(III/0) potentials showed better agreement with experimental data when using the COSMO solvation model and the test set comprising the NIST recommended ionization energies. Furthermore, the Md(II/0) potential was better described with the SMD model, whereas No(II/0) was not well described by all methods. Finally, the computational approach was able to predict redox potentials that for most cases agreed with the current available experimental or estimated data.

Actinides↗
DOE OSTI2020

Coverage-Dependent Adsorption of Hydrogen on Fe(100): Determining Catalytically Relevant Surface Structures via Lattice Gas Models

Hydrogen adatoms are a critical surface species for several reactions catalyzed by Fe surfaces such as Fischer–Tropsch synthesis, ammonia synthesis, and the hydrodeoxygenation of biomass-derived oxygenates. Parameterizing the energetics for H/Fe in terms of both coverage and configuration space can significantly aid in the development of multiscale models as well as provide atomic level insight into the dominant surface structures present under realistic reaction conditions. Here, we construct a lattice gas model for H/Fe(100), where the lateral interactions are determined from first-principles using density functional theory. Using 950 symmetrically unique H/Fe(100) configurations, we generate a cluster expansion with a predictive accuracy in terms of surface energy of 3.8 meV/site over a coverage range from 0 to 3 monolayers. Ten electronic ground state structures are identified from this thorough scan (including the structures at 0 and 3 monolayers), which were subsequently used to generate ab initio phase diagrams under a range of temperatures and pressures. Under reaction conditions typical of Fischer–Tropsch synthesis, ammonia synthesis, and biomass oxygenate hydrodeoxygenation, we find that the 1.0 monolayer structure is dominant. Furthermore, examination of the total H–H lateral interactions for the H/Fe(100) electronic ground state structures shows that H/Fe(100) can be accurately modeled via a mean-field ideal lattice gas model for coverages less than 1.0 monolayers. Altogether, this work enables the incorporation of H–H lateral interactions on Fe(100) into multiscale models, via either mean-field or site-dependent techniques, and provides atomic insight into the catalytically relevant H/Fe(100) structures for a range of heterogeneous reactions.

09 BIOMASS FUELS↗
DOE OSTI2020

Proton-coupled electron transfer across benzimidazole bridges in bioinspired proton wires

Designing molecular platforms for controlling proton and electron movement in artificial photosynthetic systems is crucial to efficient catalysis and solar energy conversion. The transfer of both protons and electrons during a reaction is known as proton-coupled electron transfer (PCET) and is used by nature in myriad ways to provide low overpotential pathways for redox reactions and redox leveling, as well as to generate bioenergetic proton currents. Herein, we describe theoretical and electrochemical studies of a series of bioinspired benzimidazole-phenol (BIP) derivatives and a series of dibenzimidazole-phenol (BI 2 P) analogs with each series bearing the same set of terminal proton-accepting (TPA) groups. The set of TPAs spans more than 6 pK a units. These compounds have been designed to explore the role of the bridging benzimidazole(s) in a one-electron oxidation process coupled to intramolecular proton translocation across either two (the BIP series) or three (the BI 2 P series) acid/base sites. These molecular constructs feature an electrochemically active phenol connected to the TPA group through a benzimidazole-based bridge, which together with the phenol and TPA group form a covalent framework supporting a Grotthuss-type hydrogen-bonded network. Infrared spectroelectrochemistry demonstrates that upon oxidation of the phenol, protons translocate across this well-defined hydrogen-bonded network to a TPA group. The experimental data show the benzimidazole bridges are non-innocent participants in the PCET process in that the addition of each benzimidazole unit lowers the redox potential of the phenoxyl radical/phenol couple by 60 mV, regardless of the nature of the TPA group. Using a series of hypothetical thermodynamic steps, density functional theory calculations correctly predicted the dependence of the redox potential of the phenoxyl radical/phenol couple on the nature of the final protonated species and provided insight into the thermodynamic role of dibenzimidazole units in the PCET process. This information is crucial for developing molecular “dry proton wires” with these moieties, which can transfer protons via a Grotthuss-type mechanism over long distances without the intervention of water molecules.

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

Higher loadings of Pt single atoms and clusters over reducible metal oxides: application to C–O bond activation

We develop higher loadings of isolated noble metal atoms and clusters on a mildly reducible metal oxide. Here we demonstrate the approach for Pt supported on TiO 2 and confirmed it by XRD, AC-HAADF-STEM, CO-FTIR, XAS, and XPS. Density functional theory calculations rationalize the experimental stability and the IR shifts using mixtures of CH 3 I and CO. The redispersed catalysts are thermally stable in inert gas or H 2 and afford enhanced selectivity and activity in hydrodeoxygenation reactions compared to metal nanoparticles by creating surface oxygen vacancies that promote C–O cleavage without side reactions. Higher metal loadings, e.g., 1%Pt/TiO 2 , on the oxide surface profoundly increase the activity of the bare oxide catalyst tenfold compared to ultra-low loadings typically used for single atom catalysis.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Real-space Kohn–Sham density functional theory for complex energy applications

Real-space Kohn-Sham density functional theory (real-space KS-DFT) enables large-scale electronic structure simulations that is particularly well-suited for the modern high-performance computing (HPC) architectures. This feature article reviews its theoretical foundations, highlights the algorithmic advances and recent developments, and showcases applications in complex nano systems. We aim to provide a perspective on the trajectory of real-space KS-DFT as an emerging tool for computational chemistry and materials science in the exascale era.

Zhang, Zeyi↗
DOE OSTI2025

An atomic cluster expansion potential for twisted multilayer graphene

Twisted multilayer graphene, characterized by its moiré patterns arising from inter-layer rotational misalignment, serves as a rich platform for exploring quantum phenomena. Machine learning interatomic potentials (MLIPs) are a promising approach to model such systems. Our work develops a method to generate training and test datasets for fitting MLIPs that capture all possible misalignments but remain small-scale to facilitate efficient data generation and parameter estimation. To achieve this, we generate configurations with periodic boundary conditions suitable for density functional theory calculations, and then introduce an internal twist and shift within those supercell structures. Using this technique, supplemented with an active learning workflow, we fit an Atomic Cluster Expansion potential for simulating twisted multilayer graphene and test it for accuracy and robustness on a range of simulation tasks.

2D materials↗
DOE OSTI2026

Excitation energy of fission fragments within nuclear time-dependent density functional theory

The number and properties of the neutrons and photons emitted in nuclear fission are directly related to the excitation energy of the fission fragments when they are formed at scission. Though not observable experimentally because of the extremely short timescales, the excitation energy of fission fragments can be predicted by microscopic theory based on time-dependent density functional theory (TDDFT). Initial results on the value of the total kinetic energy of fission reactions were very promising, but could not probe all possible fragmentations. Here, in this work, we perform large-scale TDDFT calculations in 240 Pu enabled by the development of a new TDDFT solver. We obtain TDDFT trajectories covering nearly all possible fragmentations. We find that the total kinetic energy is close to experimental values only for the most likely fission while it is severely underestimated at both small and large asymmetries. This conclusion seems rather independent of the parametrization of the energy functional, both in its particle-hole and particle-particle channels.

Physics - Nuclear physics and radiation physics↗
NASA NTRS1994

A time-dependent gyro-kinetic model of thermal ion upflows in the high-latitude F region

Ample evidence supports the significance of the high-latitude ionospheric contribution to magnetospheric plasma. Assuming flux conservation along a flux tube, the upward field-aligned ion flows observed in the magnetosphere require high-latitude ionospheric field-aligned ion upflows of the order of 10(exp 8) to 10(exp 9)/sq cm/s. Since radar and satellite observations of high-latitude F region flows at times exceed this flux requirement by an order of magnitude, the thermal ionospheric upflows are not simply the ionospheric response to a magnetospheric flux requirement. Several ionospheric ion upflow mechanisms have been proposed, but simulations based on fluid theory do not reproduce all the observed features of ionospheric ion upflows. Certain asymmetries in the statistical morphology of high-latitude F region ion upflows suggest that the ion upflows may be generated by ion-neutral frictional heating. We developed a single-component (O(+)), time-dependent gyro-kinetic model of the high-latitude F region response to frictional heating in which the neutral exobase is a discontinuous boundary between fully collisional and collisionless plasmas. The concept of a discontinuous neutreal exobase and the assumption of a constant and uniform polarization electric field reduce the ion velocity distribution function, from which we can compute the ion density, parallel velocity, parallel and perpendicular temperature, and parallel flux. Using our model, we simulated the response of a convecting flux tube between 500 km and 2500 km to various frictional heating inputs; the results were both qualitatively and quantitatively different from fluid model results, which may indicate an inadequacy of the fluid theory approach. The gyro-kinetic frictional heating model responses to the various simulations were qualitatively similar: (1) initial perturbations of all the modeled parameters propagated rapidly up the flux tube, (2) transient values of the ion parallel velocity, temperature, and flux exceeded 3 km/s, 2 x 10(exp 4) K, and 10(exp 9)/sq cm/s, respectively, (3) a second transient regime developed wherein the parallel temperature drops to very low values (a few hundred Kelvins), and (4) well after heating ceased, large parallel temperatures and large downward parallel velocities and fluxes developed as the flux tube slowly returned to diffusive equilibrium. The ion velocity distributions during the simulation are often non-Maxwellian and are sometimes composed of two distinct ion populations.

Loranc, M.↗