Search NASA⌕ Search

SEARCH · Search NASA

Results for “Classical density functional theory”

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

Deep Understanding of Strong Metal Interface Confinement: A Journey of Pd/FeO x Catalysts

Tuning the atomic interface configuration of noble metals (NMs) and transition-metal oxides is an effective straightforward yet challenging strategy to modulate the activity and stability of heterogeneous catalysts. Herein, Pd supported on mesoporous Fe 2 O 3 with a high specific surface area was rationally designed and chosen to construct the Pd/iron oxide interface. As a versatile model, the physicochemical environments of Pd nanoparticles (NPs) could be precisely controlled by taming the reduction temperature. The experimental and density functional theory calculation results unveiled that the catalyst in the support–metal interface confinement (SMIC) state showed significantly enhanced catalytic activity and sintering resistance for CO oxidation. The constructed Fe sites at the interfaces between FeO x overlayers and Pd NPs not only provided additional coordinative unsaturated ferrous sites for the adsorption and activation of O 2 , thereby facilitating the activation efficiency of O 2 , but also impressively changed the reaction pathway of CO oxidation. As a result, the catalyst followed the Pd/Fe dual-site mechanism instead of the classical Mars–van Krevelen mechanism. For the catalyst in the strong metal–support interaction (SMSI) state, its catalytic activity was seriously suppressed because of the excessive encapsulation of the active Pd sites by FeO x overlayers. Thus, the present study therefore provides detailed insights into the SMIC and SMSI in ferric oxide-supported Pd catalysts, which could guide the preparation of highly efficient supported catalysts for practical applications.

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

Classical and quantum computing of shear viscosity for ( 2 + 1 ) D SU(2) gauge theory

We perform a nonperturbative calculation of the shear viscosity for ( 2 + 1 )-dimensional SU(2) gauge theory by using the lattice Hamiltonian formulation. The retarded Green’s function of the stress-energy tensor is calculated from real time evolution via exact diagonalization of the lattice Hamiltonian with a local Hilbert space truncation, and the shear viscosity is obtained via the Kubo formula. When taking the continuum limit, we account for the renormalization group flow of the coupling but no additional operator renormalization. We find the ratio of the shear viscosity and the entropy density η s is consistent with a well-known holographic result 1 4 π at several temperatures on a 4 × 4 honeycomb lattice with the local electric representation truncated at j max = 1 2 . We also find the ratio of the spectral function and frequency ρ x y ( ω ) ω exhibits a peak structure when the frequency is small. Both the exact diagonalization method and simple matrix product state classical simulation method beyond j max = 1 2 on bigger lattices require exponentially growing resources. So we develop a quantum computing method to calculate the retarded Green’s function and analyze various systematics of the calculation including j max truncation and finite size effects, Trotter errors and the thermal state preparation efficiency. Our thermal state preparation method still requires resources that grow exponentially with the lattice size, but with a very small prefactor at high temperature. We test our quantum circuit on both the Quantinuum emulator and the IBM simulator for a small lattice and obtain results consistent with the classical computing ones. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Salt-flux synthesis, crystal structure and theoretical characterization of Rb 0.74 Ga 6.62 Ti 0·38 O 11

Here single crystals of Rb 0.74 Ga 6.62 Ti 0·38 O 11 (RGTO) were grown from a mixed RbCl–RbF flux at 850 °C. The compound crystallizes in the RbGa 7 O 11 structure type, which is reminiscent of the hollandite and β-Ga 2 O 3 structure types. RGTO crystallizes in the monoclinic space group P2/m with lattice parameters a = 8.3355 (8) Å, b = 3.0286 (3) Å, c = 9.5028 (9) Å, and β = 114.620 (3)°. The crystal structure of RGTO is comprised of GaO 6 and mixed (Ga/Ti)O 6 octahedra and GaO 4 tetrahedra connected in a complex three-dimensional, anionic framework exhibiting eight-sided channels that are occupied by disordered Rb cations required for charge balance. First-principles calculations in the form of density functional theory were performed, which indicated the complex to be a charge transfer semiconductor.

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

Low-frequency Raman signature of Ag-intercalated few-layer MoS 2

Density functional theory based calculations and experimental analysis on a limited number of real samples are performed to study how the presence of silver intercalated in the van der Waals gap of few-layer MoS 2 affects the low-frequency Raman active modes of this material. Silver is found to predominantly affect the breathing-like and shear-like vibrational modes of MoS 2 . These modes correspond to quasi-rigid movements of each individual layer with a restoring force (and, in turn, frequency) that is determined by modulations in the weak interlayer interactions. Noticeable red-shifts with increasing Ag concentration are found for all low-frequency modes. This finding indicates the potential for low-frequency vibrations as useful gauges for practical determination of silver concentration using low-frequency Raman spectroscopy. Finally, this work also describes a semi-classical linear chain model that allows to extrapolate results to a large number of layers. Further, first-principles calculations show how Raman spectroscopy can be used to characterize the quality of the two-dimensional interface between MoS 2 and a silver substrate.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Application of the locally self-consistent embedding approach to the Anderson model with non-uniform random distributions

Highlights: • Typical Medium Theory (TMT) for the Anderson Localization. • Locally Self-Consistent Multiple Scattering Method (LSMS) for Random Disordered Systems. • Linear Scaling Computational Method for Random Disordered Systems. We apply the recently developed embedding scheme for the locally self-consistent method to random disorder electrons systems. The method is based on the locally self-consistent multiple scattering theory and the typical medium theory. The locally self-consistent multiple scattering theory divides a system into many small designated local interaction zones. The subsystem within each local interaction zone is embedded in a self-consistent field from the typical medium theory. This approximation allows the study of random systems with large numbers of sites. We present results for the three dimensional Anderson model with different random disorder potential distributions. Using the typical density of states as an indicator of Anderson localization, we find that the method can capture the localization for commonly studied disorder potentials. These include the uniform distribution, the Gaussian distribution, and even the unbounded Cauchy distribution.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Energy Frontier Research Centers: Center for the Computational Design of Functional Layered Materials (CCDM) August 1, 2014 - July 31, 2018; Center for Complex Materials from First Principles (CCM) August 1, 2018 - July 31, 2021 (Final Report)

The mission of the DOE Energy Frontier Research Centers CCDM (2014-2018) and CCM (2018-2021) was to theoretically develop, computationally apply, and experimentally validate electronic structure methods for all materials, with a focus on the complex materials, especially layered and two-dimensional materials, strongly-correlated materials, and liquid water. This was achieved by over 200 published journal articles authored by about 17 senior investigators from physics and chemistry and from theory, computation, and experiment, plus their collaborators. In particular, the Centers confirmed the predictive power of the SCAN (strongly constrained and appropriately normed) density functional, which was constructed to satisfy 17 known exact constraints and several appropriate norms. Without being fitted to real bonded systems, and at a modest computational cost, SCAN correctly predicted covalent, ionic, metallic, hydrogen, and van der Waals bonds in many challenging materials. SCAN gave an improved description of defects in semiconductors, surface properties of metals, seven phases of ice, liquid water, liquid and supercooled silicon, subtle structural distortions in ferroelectrics, formation energies and structural predictions for solids, and critical pressures for structural phase transitions. Perhaps most remarkably, SCAN correctly described some strongly-correlated materials that were previously believed to be beyond the reach of density-functional approximations. SCAN is the only density functional that correctly predicts the band gap closing under chemical doping of the cuprate high-temperature superconducting materials. SCAN also predicts a landscape of competing stripe and magnetic phases in the cuprates. For some materials with some codes, SCAN has convergence problems that are greatly reduced by the CCM-developed r 2 SCAN, without loss of accuracy or rigor. SCAN and r 2 SCAN still make some self-interaction error, which is greatly reduced by the CCDM/ CCM-developed local orbital scaling correction (LOSC). These Centers further proved that the fundamental energy gaps of a solid from an orbital energy difference and from total energy differences are the same for a large class of generalized Kohn-Sham (GKS) functionals, including SCAN and standard hybrid functionals, and that symmetry breaking arises when a dynamic density fluctuation drops to zero frequency. The Centers identified new mechanisms for catalysis in layered materials with ions intercalated between the layers, investigated charge density waves both in model systems and in real layered materials, studied changes of band gap with the number of layers, and explored topological ultrathin films, bent nanoribbons, and defects.

08 HYDROGEN↗
DOE OSTI2020

Interlayer magnetism in Fe 3-x GeTe 2

Fe 3 - x GeTe 2 is a layered van der Waals magnetic material with a relatively high ordering temperature and large anisotropy. While most studies have concluded the interlayer ordering to be ferromagnetic, there have also been reports of interlayer antiferromagnetism in Fe 3 - x GeTe 2 . In this work, we investigate the interlayer magnetic ordering by neutron diffraction experiments, scanning tunneling microscopy (STM) and spin-polarized STM measurements, density functional theory plus U calculations, and STM simulations. We conclude that the layers of Fe 3 - x GeTe 2 are coupled ferromagnetically and that in order to capture the magnetic and electronic properties of Fe 3 - x GeTe 2 within density functional theory, Hubbard U corrections need to be taken into account.

2-dimensional systems↗
DOE OSTI2020

NWChem: Past, present, and future

Specialized computational chemistry packages have permanently reshaped the landscape of chemical sciences by providing tools to support and guide the experimental effort and for prediction of chemical and materials properties. In this regard, a special role has been played by electronic structure packages where complex chemical and materials processes can be modeled using first-principle-driven methodologies. Over the last few decades, the rapid development of computing technologies and a tremendous increase in computational power has offered a unique chance to study complex chemical transformations using sophisticated and predictive many-body techniques to describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, a critical role has been played by novel parallel algorithms capable of taking advantage of computational resources to address polynomial scaling of electronic structure methods. NWChem was among the first electronic structure codes that focused on delivering scalable performance for electronic structure simulations. Herein, we briefly review the NWChem suite of computational codes including its history, design principles, parallel tools, current capabilities, outreach and outlook.

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

Spatiotemporal mapping of mesoscopic liquid dynamics

The study of liquid dynamics at mesoscopic scales is still strewn with difficulty due to limitations in theory and experiment. Historically, significant attention has been given to the analysis of space-time correlation functions and their frequency-Fourier transforms at a few discrete wave numbers. The massive computing power afforded by modern high performance computing clusters and the advent of a wide-angle neutron spin-echo spectrometer, however, have unlocked a more intuitive and fruitful approach to this problem. Using molecular dynamics simulations, here we demonstrate the benefits of spatiotemporally mapping intermediate scattering functions on a dense grid of correlation times and wave numbers. Four model systems are investigated: a Lennard-Jones liquid, a coarse-grained bead-spring polymer, a molten sodium chloride, and a poly(ethylene oxide) melt. We show that the spatiotemporal mapping approach is particularly useful for elucidating the mesoscopic dynamics in these liquids, where several underlying mechanisms, such as molecular relaxations, hydrodynamic modes, and nonhydrodynamic excitations, are potentially at play. Therefore, compared to the traditional method, direct visualization of density space-time correlation functions on two-dimensional color maps permits appraisals of complicated dynamical behavior at mesoscales in a global manner. For example, the scaling relations between space and time for different types of molecular motions can be straightforwardly identified on these plots, without any model-dependent analysis. Additionally, we show how theoretical ideas regarding collective mesoscopic dynamics, such as the classical hydrodynamic theory, the convolution approximation, and a recently proposed phenomenological model, can be discussed in terms of the global features of spatiotemporal maps of intermediate scattering functions. The new perspective offered by the spatiotemporal mapping method should prove useful for the study of liquid dynamics in general.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Accretion of matter and spectra of binary X-ray sources in massive gravity

We study low-mass binary X-ray sources involving stellar mass black holes within massive gravity. Regarding the accretion disk, we adopt the standard model for an optically thick, cool, and geometrically thin disk by Shakura–Sunyaev. For the gravitational field generated by the black hole, we consider the analogue of the Schwarzschild–de Sitter space–time of Einstein’s theory in massive gravity, for which we found an additional term linear in the radial coordinate. Then, we compute the radial velocity, the energy density and the pressure as a function of the radial coordinate, and the X-ray emission’s soft spectral component expected from the disk. We also investigated in detail the impact of this new geometry. Our result indicates that by using observed spectra from confirmed X-ray binaries involving astrophysical black holes, we can put strong constraints on alternative theories of gravity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Methanol Partial Oxidation on Cu(111) and PtCu(111) Single-Atom Alloy Surfaces: Effect of Surface Oxygen Coverage on Selectivity

The selective oxidation of methanol to formaldehyde on Cu surfaces is an important and well-studied reaction. However, a systematic analysis of product selectivity as a function of oxygen coverage on Cu(111) and Cu-based single-atom alloys (SAAs) has not been previously reported. In this work, we present a comprehensive investigation of deuterated methanol (CD 3 OH) partial oxidation on Cu(111) and 1% PtCu(111) SAA surfaces as a function of preadsorbed oxygen coverage. Temperature-programmed desorption (TPD) and X-ray photoelectron spectroscopy (XPS) reveal that isolated Pt atoms inhibit the initial surface oxidation of Cu(111) under low oxygen exposures. Despite this difference in oxidation kinetics, the product selectivity remains largely unaffected: on both Cu(111) and PtCu(111), formaldehyde (CD 2 O) is the predominant partial oxidation product over a broad range of oxygen coverages. The selectivity toward formaldehyde peaks at intermediate oxygen coverages (∼0.3 monolayers, ML), indicating the existence of an optimal oxygen loading for partial oxidation. Notably, the similar product selectivities on Cu(111) and PtCu(111) over a range of surface oxygen coverage indicate that Pt single atoms do not significantly alter the reaction pathway or shift the optimal oxygen coverage for formaldehyde formation. Control experiments confirm that Cu(111) is unreactive toward methanol in the absence of oxygen, while PtCu(111) surfaces produce a small amount of formaldehyde even when oxygen is not preadsorbed, indicating that isolated Pt atoms facilitate O−H activation at below 150 K, leading to H 2 desorption, followed by C−D activation at higher temperatures (∼350 K). Density functional theory (DFT)-based calculations show that Pt atoms increase the O 2 dissociation barrier relative to Cu(111), consistent with the observed inhibition of oxidation at low exposures. Overall, this work provides the first detailed selectivity map for methanol oxidation on oxidized Cu(111) and PtCu(111) SAA surfaces. By linking classical mechanistic insights such as methoxy- and formate-mediated pathways with single-atom alloy catalyst design, this work demonstrates that while Pt substitution modulates the oxidation kinetics and oxygen binding, the overall selectivity toward formaldehyde is governed primarily by oxygen coverage. These findings underscore the potential of isolated dopants to tune surface oxidation behavior without compromising the intrinsic partial oxidation selectivity of copper-based catalysts.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Diagonalizing the Born–Oppenheimer Hamiltonian via Moyal perturbation theory, nonadiabatic corrections, and translational degrees of freedom

This article describes a method for calculating higher order or nonadiabatic corrections in Born–Oppenheimer theory and its interaction with the translational degrees of freedom. The method uses the Wigner–Weyl correspondence to map nuclear operators into functions on the classical phase space and the Moyal star product to represent operator multiplication on those functions. These are explained in the body of the paper. The result is a power series in κ 2 , where κ = (m/M)1/4 is the usual Born–Oppenheimer parameter. The lowest order term is the usual Born–Oppenheimer approximation, while higher order terms are nonadiabatic corrections. These are needed in calculations of electronic currents, momenta, and densities. The separation of nuclear and electronic degrees of freedom takes place in the context of the exact symmetries (for an isolated molecule) of translations and rotations, and these, especially translations, are explicitly incorporated into our discussion. This article presents an independent derivation of the Moyal expansion in molecular Born–Oppenheimer theory. Here, we show how electronic currents and momenta can be calculated within the framework of Moyal perturbation theory; we derive the transformation laws of the electronic Hamiltonian, the electronic eigenstates, and the derivative couplings under translations; we discuss in detail the rectilinear motion of the molecular center of mass in the Born–Oppenheimer representation; and we show how the elimination of the translational components of the derivative couplings leads to a unitary transformation that has the effect of exactly separating the translational degrees of freedom.

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

Integrated Quantum-Classical Protocol for the Realistic Description of Solvated Multinuclear Mixed-Valence Transition-Metal Complexes and Their Solvatochromic Properties

Linear cyanide-bridged polymetallic complexes, which undergo photoinduced metal-to-metal charge transfer, represent prototypical systems for studying long-range electron-transfer reactions and understanding the role played by specific solute–solvent interactions in modulating the excited-state dynamics. Here, to tackle this problem, while achieving a statistically meaningful description of the solvent and of its relaxation, one needs a computational approach capable of handling large polynuclear transition-metal complexes, both in their ground and excited states, as well as the ability to follow their dynamics in several environments up to nanosecond time scales. Here, we present a mixed quantum classical approach, which combines large-scale molecular dynamics (MD) simulations based on an accurate quantum mechanically derived force field (QMD-FF) and self-consistent QMD polarized point charges, with IR and UV–vis spectral calculations to model the solvation dynamics and optical properties of a cyano-bridged trinuclear mixed-valence compound (trans-[(NC) 5 Fe III (μ-CN)Ru II (pyridine) 4 (μ-NC)Fe III (CN) 5 ] 4– ). We demonstrate the reliability of the QMD-FF/MD approach in sampling the solute conformational space and capturing the local solute–solvent interactions by comparing the results with higher-level quantum mechanics/molecular mechanics (QM/MM) MD reference data. The IR spectra calculated along the classical MD trajectories in different solvents correctly predict the red shift of the CN stretching band in the aprotic medium (acetonitrile) and the subtle differences measured in water and methanol, respectively. By explicitly including the solvent molecules around the cyanide ligands and calculating the thermal averaged absorption spectra using time-dependent density functional theory calculations within the Tamm–Dancoff approximation, the experimental solvatochromic shift is quantitatively reproduced going from water to methanol, while it is overestimated for acetonitrile. This discrepancy can likely be traced back to the lack of important dispersion interactions between the solvent cyano groups and the pyridine substituents in our micro solvation model. The proposed protocol is applied to the ground state in water, methanol, and acetonitrile and can be flexibly generalized to study excited-state nonequilibrium solvation dynamics.

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

Plasmon-induced excitation energy transfer in silver nanoparticle dimers: A real-time TDDFTB investigation

Here, using real-time quantum dynamics calculations, we perform theoretical investigations of light-induced interactions and electronic excitation transfer in a silver nanoparticle dimer. Real-time time-dependent density functional tight-binding (RT-TDDFTB) calculations provide details of the quantum dynamical processes at an electronic/atomistic level with attosecond resolution. The computational efficiency of RT-TDDFTB allows us to examine electronic dynamics up to picosecond time scales. With time scales varying over six orders of magnitude, we provide insight into interactions between the nanoparticle and laser and between nanoparticles. Our results show that the coupling between nanoparticle monomers is dependent on the separation distance between the nanoparticles in the dimer. As the interparticle distance is varied, the dipole–dipole interactions and electronic excitation transfer mechanisms are markedly different. At large distances (from 50 to 20 Å), the energy transfer from NP1 to NP2 becomes more efficient as the interparticle distance decreases. The total dipole moment of the Ag 14 nanoparticle dimer increases linearly at an interparticle distance of 20 Å and reaches its maximum after 1.2 ps. The electronic excitation transfer is also the most efficient at 20 Å. At short distances, back-transfer effects reduce the ability of the dimer and NP1 to accept energy from the incident electric field. We attribute the distance-dependent features of the nanoparticle dimer to the beating between the laser acting on NP1 and the back transfer from NP2 to NP1.

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

Supercurrent decay in ballistic magnetic Josephson junctions

We investigate transport properties of ballistic magnetic Josephson junctions and establish that suppression of supercurrent is an intrinsic property of the junctions, even in absence of disorder. By studying the role of ferromagnet thickness, magnetization, and crystal orientation we show how the supercurrent decays exponentially with thickness and identify two mechanisms responsible for the effect: (i) large exchange splitting may gap out minority or majority carriers leading to the suppression of Andreev reflection in the junction, (ii) loss of synchronization between different modes due to the significant dispersion of the quasiparticle velocity with the transverse momentum. Our results for Nb/Ni/Nb junctions are in good agreement with recent experimental studies. Our approach combines density functional theory and the Bogoliubov-de Gennes model and opens a path for material composition optimization in magnetic Josephson junctions and superconducting magnetic spin valves.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Molecular insights into CO 2 -to-bicarbonate transformation in functionalized anion exchange ionomers for electrochemical separations

Bipolar membrane (BPM) electrochemical processes are a promising platform for carbon dioxide (CO 2 ) separations, but the molecular level thermodynamic and kinetic understanding of CO 2 -to-bicarbonate (HCO 3 − ) transformation remain poorly understood. This study employs a multiscale computational approach to systematically explore the adsorption and reactive transformation of CO 2 in five anion exchange ionomer systems. Classical molecular dynamics (MD) simulation results demonstrate that polymers with imidazolium groups significantly reduce CO 2 diffusion and enhance (OH − )–CO 2 interactions due to stronger electrostatic and π-interactions. Compared to the commonly used quaternary ammonium ionomers, imidazolium-functionalized ionomers show improved CO 2 proximity and interaction strength. Ab initio MD and density functional theory (DFT) calculations reveal that the benzyl-substituted imidazolium (IM-Ben) substantially reduces the energy barrier for HCO 3 − formation (∼72 meV lower) compared to the alkyl-substituted IM-nBu, while also mitigating imidazolium deprotonation under moderate hydration conditions. Transition state analysis shows IM-Ben forms more extensive hydrogen-bonding networks, which stabilize the transition state structure and contribute to a lower energy barrier for bicarbonate formation. These findings highlight the advantage of the adjacent benzyl moiety in enabling efficient CO 2 -to-bicarbonate transformation via hydrated hydroxide ion counterions, offering mechanistic insights and clear molecular design principles for optimizing anion exchange ionomers at bipolar membrane interfaces for electrochemical CO 2 separation applications.

Bipolar membranes, Reactive transformation of CO2,↗
DOE OSTI2021

Influence of solution phase environmental heterogeneity and fluctuations on vibronic spectra: Perylene diimide molecular chromophore complexes in solution

Ensembles of ab initio parameterized Frenkel-exciton model Hamiltonians for different perylene diimide dimer systems are used, together with various dissipative quantum dynamics approaches, to study the influence of the solvation environment and fluctuations in chromophore relative orientation and packing on the vibronic spectra of two different dimer systems: a π-stacked dimer in aqueous solution in which the relative chromophore geometry is strongly confined by a phosphate bridge and a side-by-side dimer in dichloromethane involving a more flexible alkyne bridge that allows quasi-free rotation of the chromophores relative to one another. These entirely first-principles calculations are found to accurately reproduce the main features of the experimental absorption spectra, providing a detailed mechanistic understanding of how the structural fluctuations and environmental interactions influence the vibronic dynamics and spectroscopy of solutions of these multi-chromophore complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗