Search NASA⌕ Search

SEARCH · Search NASA

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

Influence of Heat Treatments in H 2 and Ar on the E 1 center in β -Gs 2 O 3

In this report the influence of heat treating n-type bulk β-Ga 2 O 3 in hydrogen (H2) and argon (Ar) gases on the presence of the defect level commonly labeled as E 1 was studied. Fourier transform-infrared spectroscopy confirms that hydrogen (H) is incorporated into β-Ga 2 O 3 during H 2 annealing at 900 °C. Deep-level transient spectroscopy measurements reveal that the concentration of the E 1 level is promoted by the introduction of H, in contrast to what is observed in samples heat-treated in an Ar flow. We further find the E 1 level to be stable against heat treatments at 650 K, both with and without an applied reverse-bias voltage. Potential candidates for the defect origin of E 1 are investigated using hybrid-functional calculations, and three types of defect complexes involving H are found to exhibit charge-state transition levels compatible with E 1 , including substitutional H at one of the threefold coordinated O sites, Ga-substitutional shallow donor impurities passivated by H, and certain configurations of singly hydrogenated Ga–O divacancies. Among these types, only the latter exhibit H binding energies that are consistent with the observed thermal stability of E 1 .

36 MATERIALS SCIENCE↗
DOE OSTI2021

Computational catalyst discovery: Active classification through myopic multiscale sampling

We report the recent boom in computational chemistry has enabled several projects aimed at discovering useful materials or catalysts. We acknowledge and address two recurring issues in the field of computational catalyst discovery. First, calculating macro-scale catalyst properties is not straightforward when using ensembles of atomic-scale calculations [e.g., density functional theory (DFT)]. We attempt to address this issue by creating a multi-scale model that estimates bulk catalyst activity using adsorption energy predictions from both DFT and machine learning models. The second issue is that many catalyst discovery efforts seek to optimize catalyst properties, but optimization is an inherently exploitative objective that is in tension with the explorative nature of early-stage discovery projects. In other words, why invest so much time finding a “best” catalyst when it is likely to fail for some other, unforeseen problem? We address this issue by relaxing the catalyst discovery goal into a classification problem: “What is the set of catalysts that is worth testing experimentally?” Here, we present a catalyst discovery method called myopic multiscale sampling, which combines multiscale modeling with automated selection of DFT calculations. It is an active classification strategy that seeks to classify catalysts as “worth investigating” or “not worth investigating” experimentally. Our results show an ~7–16 times speedup in catalyst classification relative to random sampling. These results were based on offline simulations of our algorithm on two different datasets: a larger, synthesized dataset and a smaller, real dataset.

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

Ultra-large-scale ab initio quantum chemical computation of bio-molecular systems: The case of spike protein of SARS-CoV-2 virus

The COVID-19 pandemic poses a severe threat to human health with an unprecedented social and economic disruption. Spike (S) glycoprotein of the SARS-CoV-2 virus is pivotal in understanding the virus anatomy, since it initiates the first contact with the ACE2 receptor in the human cell. We report results of ab initio computation of the spike protein, the largest ab initio quantum chemical computation to date on any bio-molecular system, using a divide and conquer strategy by focusing on individual structural domains. In this approach we divided the S-protein into seven structural domains: N-terminal domain (NTD), receptor binding domain (RBD), subdomain 1 (SD1), subdomain 2 (SD2), fusion peptide (FP), heptad repeat 1 with central helix (HR1-CH) and connector domain (CD). The entire Chain A has 14,488 atoms including the hydrogen atoms but excluding the amino acids with missing coordinates based on the PDB data (ID: 6VSB). The results include structural refinement, ab initio calculation of intramolecular bonding mechanism, 3- dimensional non-local inter-amino acid interaction with implications for the inter-domain interaction. Details of the electronic structure, interatomic bonding, partial charge distribution and the role played by hydrogen bond network are discussed. In the interaction among structural domains, we present new insights for crucial hinge-like movement and fusion process. Extension of such calculation to the interface between the S-protein binding domain and ACE2 receptor can provide a pathway for computational understanding of mutations and the design of therapeutic drugs to combat the COVID-19 pandemic.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Mechanistic Insights and Rational Design of Ca-Doped CeO 2 Catalyst for Acetic Acid Ketonization

Carboxylic acid ketonization has recently gained significant attention to produce biomass-derived hydrocarbon fuels as it not only removes the highly reactive carboxylic functional group but also increases the size of the carbon chain. In this work, Ca-doped CeO 2 -based catalysts were investigated for acetic acid ketonization using a combined experimental and computational approach. Acetic acid conversion was performed across a range of temperatures including higher temperatures relevant to catalytic hot gas filtration (450 °C). Ca addition slightly decreases overall acetic acid ketonization reactivity yet stabilizes the catalyst at the higher temperatures necessary for catalytic hot gas filtration. From density functional theory calculations of the ketonization reaction mechanism, the C–C coupling and water formation steps are identified as two of the most energy-consuming steps on a CeO 2 surface with a proximal oxygen vacancy and the presence of a Ca dopant stabilizes the key intermediates. Calculations predict an optimal structure comprising three Ca ensembles to minimize the reaction free energies for C–C coupling and water formation steps. These findings provide a priori information to guide future experiments for ketonization catalyst design and development.

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

Density functional theory study on the coupling and reactions of diferuloylputrescine as a lignin monomer

Diferuloylputrescine has been found in a variety of plant species, and recent work has provided evidence of its covalent bonding into lignin. Results from nuclear magnetic resonance spectroscopy revealed the presence of bonding patterns consistent with homo-coupling of diferuloylputrescine and the possibility of cross-coupling with lignin. In the present work, density functional theory calculations have been applied to assess the energetics associated with radical coupling, rearomatization, and dehydrogenation for possible homo-coupled dimers of diferuloylputrescine and cross-coupled dimers of diferuloylputrescine and coniferyl alcohol. Furthermore, the values obtained for these reaction energetics are consistent with those reported for monolignols and other novel lignin monomers. As such, this study shows that there would be no thermodynamic impediment to the incorporation of diferuloylputrescine into the lignin polymer and its addition to the growing list of non-canonical lignin monomers.

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

Solution phase high repetition rate laser pump x-ray probe picosecond hard x-ray spectroscopy at the Stanford Synchrotron Radiation Lightsource

We present a dedicated end-station for solution phase high repetition rate (MHz) picosecond hard x-ray spectroscopy at beamline 15-2 of the Stanford Synchrotron Radiation Lightsource. A high-power ultrafast ytterbium-doped fiber laser is used to photoexcite the samples at a repetition rate of 640 kHz, while the data acquisition operates at the 1.28 MHz repetition rate of the storage ring recording data in an alternating on-off mode. The time-resolved x-ray measurements are enabled via gating the x-ray detectors with the 20 mA/70 ps camshaft bunch of SPEAR3, a mode available during the routine operations of the Stanford Synchrotron Radiation Lightsource. As a benchmark study, aiming to demonstrate the advantageous capabilities of this end-station, we have conducted picosecond Fe K-edge x-ray absorption spectroscopy on aqueous [FeII(phen)3]2+, a prototypical spin crossover complex that undergoes light-induced excited spin state trapping forming an electronic excited state with a 0.6–0.7 ns lifetime. In addition, we report transient Fe Kβ main line and valence-to-core x-ray emission spectra, showing a unique detection sensitivity and an excellent agreement with model spectra and density functional theory calculations, respectively. Notably, the achieved signal-to-noise ratio, the overall performance, and the routine availability of the developed end-station have enabled a systematic time-resolved science program using the monochromatic beam at the Stanford Synchrotron Radiation Lightsource.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Nucleation and Initial Stages of Growth during the Atomic Layer Deposition of Titanium Oxide on Mesoporous Silica

A chemical approach to the deposition of thin films on solid surfaces is highly desirable but prone to affect the final properties of the film. To better understand the origin of these complications, the initial stages of the atomic layer deposition of titania films on silica mesoporous materials were characterized. Adsorption-desorption measurements indicated that the films grow in a layer-by-layer fashion, as desired, but initially exhibit surprisingly low densities, about one-quarter of that of bulk titanium oxide. Electron microscopy, X-ray diffraction, UV/visible, and X-ray absorption spectroscopy data pointed to the amorphous nature of the first monolayers, and EXAFS and 29 Si CP/MAS NMR results to an initial growth via the formation of individual tetrahedral Ti-oxide units on isolated Si-OH surface groups with unusually long Ti-O bonds. Here, density functional theory calculations were used to propose a mechanism where the film growth starts at the nucleation centers to form an open 2D structure.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Enhanced CO 2 Reactive Capture and Conversion Using Aminothiolate Ligand–Metal Interface

Metallic catalyst modification by organic ligands is an emerging catalyst design in enhancing the activity and selectivity of electrocatalytic carbon dioxide (CO 2 ) reactive capture and reduction to value-added fuels. However, a lack of fundamental science on how these ligand–metal interfaces interact with CO 2 and key intermediates under working conditions has resulted in a trial-and-error approach for experimental designs. With the aid of density functional theory calculations, we provided a comprehensive mechanism study of CO 2 reduction to multicarbon products over aminothiolate-coated copper (Cu) catalysts. Our results indicate that the CO 2 reduction performance was closely related to the alkyl chain length, ligand coverage, ligand configuration, and Cu facet. The aminothiolate ligand–Cu interface significantly promoted initial CO 2 activation and lowered the activation barrier of carbon–carbon coupling through the organic (nitrogen (N)) and inorganic (Cu) interfacial active sites. Experimentally, the selectivity and partial current density of the multicarbon products over aminothiolate-coated Cu increased by 1.5-fold and 2-fold, respectively, as compared to the pristine Cu at –1.16 V RHE , consistent with our theoretical findings. Furthermore, this work highlights the promising strategy of designing the ligand–metal interface for CO 2 reactive capture and conversion to multicarbon products.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Unraveling the structural dependency of Weyl nodes in Co 2 MnGa

Conventionally, the modulation of the intrinsic Weyl nodes in Weyl semimetals is challenging, due to topological protection. Here we report the structural dependence of the Weyl nodes in a Co2⁢MnGa Heusler thin film via a temperature-dependent tetragonal distortion. The ability to manipulate these Weyl nodes allows for the control of the intrinsic electromagnetic properties. Temperature-dependent x-ray diffraction (XRD) measurements identify a compressive tetragonal distortion with decreasing temperature from 300 to 20 K. The calculated Weyl properties can be directly compared with experimental parameters through the temperature-dependent XRD measurements which show the intrinsic correlation between Weyl properties and important magnetic parameters. Further, the microscopic momentum space properties of Weyl nodes such as the distance (d W ), solid angle (Ω W ), tilt (φ W ), and nodal point energy (E W ) directly affect the macroscopic observable properties such as exchange stiffness (A), magnetization (M), and effective anisotropy field $H^{eff}_{K}$, as shown via structure-dependent density functional theory calculations. These predictions are experimentally observed as large variations in the bulk magnetization and effective anisotropy field as a function of temperature. These results highlight a unique degree of freedom in the control of macroscopic magnetic properties via the modulation of the intrinsic properties of Weyl nodes through structural distortions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Computational descriptor for electrochemical currents of carbon dioxide reduction on Cu facets

Computation screening is crucial for designing efficient electrochemical catalysts for carbon dioxide (CO 2 R) reduction that produce valuable hydrocarbons and oxygenates. In this work, leveraging density functional theory calculations for the CO adsorption energy ΔE CO on seventeen Cu terminations, we discover a strong linear correlation between ΔE CO and the experimentally measured CO 2 R electrochemical currents (ACS Catal. 2022, 12, 11, 6578–6588). Examining ab initio thermodynamics of early critical intermediates CO*, COH*, and CHO*, we find that CO* → CHO* is the thermodynamically controlling step. Beyond the general CO adsorption energy that only shows a linear trend with CO 2 R activity, we show that the reaction free energy of CO* → CHO* is the descriptor for the overall CO 2 R activity for Cu facets, as it displays a volcano relationship with the experimental current. Importantly, we show that high step and kink density of the Cu terminations not only enhances CO adsorption strength but also modulates the CO* → CHO* pathway, as respectively exemplified in the (941) and (741) facets. In addition, we explain that the high activity of (741) is due to its relatively low hydrogen evolution reaction activity compared with the other Cu surfaces.

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

Dielectronic satellite emission from a solid-density Mg plasma: Relationship to models of ionization potential depression

Here, we report on experiments where solid-density Mg plasmas are created by heating with the focused output of the Linac Coherent Light Source x-ray free-electron laser. We study the K-shell emission from the helium- and lithium-like ions using Bragg crystal spectroscopy. Observation of the dielectronic satellites in lithium-like ions confirms that the M-shell electrons appear bound for these high charge states. An analysis of the intensity of these satellites indicates that when modeled with an atomic-kinetics code, the ionization potential depression model employed needs to produce depressions for these ions which lie between those predicted by the well known Stewart-Pyatt and Ecker-Kroll models. These results are largely consistent with recent density functional theory calculations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Revisiting the structures and phase transitions of Ba 2 NaNb 5 O 15

The room-temperature and low-temperature structure(s) of Ba 2 NaNb 5 O 15 (BNN) have been debated since the structure was proposed in the 1960s. This work revisits the structures and phase transitions of BNN, combining high-resolution X-ray and neutron powder diffraction with density functional theory calculations. Temperature-dependent high-resolution X-ray powder diffraction patterns are collected from 4 to 918 K, and sequential batch Rietveld refinement using a symmetry mode approach to describe the structure is used to extract the main structural changes as a function of temperature. The data show that the average structure of BNN is best described by the Ama2 space group, and no other structural phase transitions were observed below the ferroelastic transition. The symmetry mode analysis, combining results from diffraction and density functional theory, shows significant octahedral tilting and corrugations of both the A1 and A2 sites along the c direction. A strong correlation between the spontaneous strain and the octahedral tilting was observed, and a potential connection with emerging microstructure at low temperatures is proposed, all enabled by the symmetry mode approach used in this work.

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

VASP calculations on Chicoma: CPU vs. GPU [Slides]

Density functional theory (DFT) calculations with the VASP code provide insight into materials behavior. In particular, the behavior of Pu still poses challenges to our understanding. Calculations for a 108-atom system that represents delta phase Pu with one substitutional Ga show a significant speedup running on GPUs vs CPUs. The gain in speed (2.7), however, remains well below the ratio of SU per node hour conversions, 540/108 = 5.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Towards Quantum-Chemical Level Calculations of SARS-CoV-2 Spike Protein Variants of Concern by First Principles Density Functional Theory

The spike protein (S-protein) is a crucial part of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), with its many domains responsible for binding, fusion, and host cell entry. In this review we use the density functional theory (DFT) calculations to analyze the atomic-scale interactions and investigate the consequences of mutations in S-protein domains. We specifically describe the key amino acids and functions of each domain, which are essential for structural stability as well as recognition and fusion processes with the host cell; in addition, we speculate on how mutations affect these properties. Such unprecedented large-scale ab initio calculations, with up to 5000 atoms in the system, are based on the novel concept of amino acid–amino acid-bond pair unit (AABPU) that allows for an alternative description of proteins, providing valuable information on partial charge, interatomic bonding and hydrogen bond (HB) formation. In general, our results show that the S-protein mutations for different variants foster an increased positive partial charge, alter the interatomic interactions, and disrupt the HB networks. We conclude by outlining a roadmap for future computational research of biomolecular virus-related systems.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Canted antiferromagnetic order and spin dynamics in the honeycomb-lattice compound Tb 2 Ir 3 Ga 9

Here, single-crystal neutron diffraction, inelastic neutron scattering, bulk magnetization measurements, and first-principles calculations are used to investigate the magnetic properties of the honeycomb lattice Tb 2 Ir 3 Ga 9 . While the R ln 2 magnetic contribution to the low-temperature entropy indicates a J eff = 1 / 2 moment for the lowest-energy crystal-field doublet, the Tb 3 + ions form a canted antiferromagnetic structure below 12.5 K. Due to the Dzyaloshinskii-Moriya interactions, the Tb moments in the a b plane are slightly canted towards b with a canted moment of 1.22 μ B per formula unit. A minimal x x z spin Hamiltonian is used to simultaneously fit the spin-wave frequencies along the high-symmetry directions and the field dependence of the magnetization along the three crystallographic axes. Long-range magnetic interactions for both in-plane and out-of-plane couplings up to the second nearest neighbors are needed to account for the observed static and dynamic properties. The z component of the exchange interactions between Tb moments is larger than the x and y components. This compound also exhibits bond-dependent exchange with negligible nearest-neighbor exchange coupling between moments parallel and perpendicular to the 4 f orbitals. Despite the J eff = 1 / 2 moments, the spin Hamiltonian is denominated by a large in-plane anisotropy K z ~ – 1 meV . DFT calculations confirm the antiferromagnetic ground state and the substantial interplane coupling at larger Tb-Tb distances.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Co(salen)-Catalyzed Oxidation of Lignin Models to Form Benzoquinones and Benzaldehydes: A Computational and Experimental Study

Lignin is a highly abundant polyphenolic polymer that imparts mechanical strength to plant biomass. Transition-metal complexes can catalyze lignin oxidation to produce value-added products, but low catalytic efficiency has hampered their use in industry. Identifying the chemical and structural factors that govern catalytic activity is a prerequisite to rational design of catalysts with improved activity. In this study, we combine computational and experimental approaches to investigate the mechanism of Co(salen)-catalyzed oxidation of the monomeric lignin models syringyl (S), vanillyl (G), and 4-hydroxybenzyl alcohol (H) to produce benzoquinone and benzaldehyde products. Experimentally, S oxidation to form dimethoxybenzoquinone proceeded efficiently with a Co(salen) catalyst coordinated by a pyridine ligand, but G and H did not undergo oxidation. Density functional theory calculations reveal that catalyst regeneration is energetically unfavorable in the presence of H, which prevents oxidation. In contrast, S readily facilitates catalyst regeneration. Formation of methoxybenzoquinone from G was achieved experimentally by adding bulky, noncoordinating bases. These findings provide a fundamental baseline for enhancing the activity of Co-Schiff base catalysts toward lignin-like molecules by adding sterically hindered nitrogenous bases or potentially by including a cocatalyst that promotes catalyst regeneration.

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

Two-dimensional Dirac semimetal based on the alkaline earth metal CaP 3

Using an evolutionary algorithm in combination with first-principles density-functional theory calculations, we identify a two-dimensional (2D) CaP 3 monolayer as a new Dirac semimetal due to inversion and nonsymmorphic spatial symmetries of the structure. This new topological material, composed of light elements, exhibits high structural stability (higher than the phase known in the literature), which is confirmed by thermodynamic and kinetic stability analysis. Moreover, it satisfies the electron filling criteria, so that its Dirac state is located near the Fermi level. The existence of the Dirac state predicted by the theoretical symmetry analysis is also confirmed by first-principles electronic band structure calculations. We find that the energy position of the Dirac state can be tuned by strain, while the Dirac state is unstable against an external electric field since it breaks the spatial inversion symmetry. In conclusion, our findings should be instrumental in the development of 2D Dirac fermions based on light elements for their application in nanoelectronic devices and topological electronics.

2-dimensional systems↗
DOE OSTI2026

First-principles simulation of a shocked H-He mixture along the principal Hugoniot

Recent laser-shock experiments on an H–He mixture containing 11% helium (atomic fraction) have suggested the presence of an immiscibility region inside Jupiter. Reflectivity measurements were used as the primary diagnostic of H–He demixing, with discontinuities in the optical reflectivity proposed as a signature of phase separation under conditions relevant to Jupiter's interior. Here, we investigate shock-compressed H–He using ab initio molecular dynamics simulations with optical properties evaluated within the Kubo–Greenwood formalism. The equation of state and ionic configurations were obtained using the thermal Tr 2 SCANL meta-GGA exchange–correlation (XC) functional, while optical properties were computed using the recently developed RS-KDT0 range-separated thermal hybrid XC, which provides state-of-the-art accuracy for band-gap predictions in the warm dense matter regime. The calculated reflectivity shows overall good agreement with experimental measurements; however, no discontinuity is observed at elevated temperatures. Moreover, the reflectivity predictions for the mixed system are consistent with the experimental measurements in the temperature range where the mixture is inferred to be demixed. Furthermore, these results suggest that reflectivity alone may not provide a unique or sensitive diagnostic of H–He demixing at low helium concentrations under these conditions.

Density functional calculations↗