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

First-principles wave-vector- and frequency-dependent exchange-correlation kernel for jellium at all densities

Here we propose a spatially and temporally nonlocal exchange correlation (XC) kernel for the spin-unpolarized fluid phase of ground-state jellium for use in time-dependent density functional and linear response calculations. The kernel is constructed to satisfy known properties of the exact XC kernel to accurately describe the correlation energies of bulk jellium and to satisfy frequency-moment sum rules at a wide range of bulk jellium densities, including those low densities that display strong correlation and symmetry breaking. These effects are easier to understand in the simple jellium model than in real systems. All exact constraints satisfied by the recent MCP07 kernel are maintained in the revised MCP07 (rMCP07) kernel, while others are added. The revision $f^{rMCP07}_{XC}$ (q, ω) differs from MCP07 only for nonzero frequencies ω. Only at densities much lower than those of real bulk metals is the frequency dependence of the kernel important for the correlation energy of jellium. As the wave vector q tends to zero, the kernel has a -4$πα(ω)/q^2$ divergence whose frequency-dependent ultranonlocality coefficient $α(ω)$ vanishes in jellium, and is predicted by rMCP07 to be extremely small for the real metals Al and Na.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Revisiting W–ZrC interfaces: A first principles study

We investigate W–ZrC interfaces using first-principles calculations based on the density functional theory. There have been theoretical and experimental studies exploring W–ZrC interfaces, however, the debate regarding the most stable interface continues to persist. In this study, we systematically simulated various W–ZrC interfaces merging W and ZrC surfaces with different orientations. Subsequently, we evaluated their stabilities and explained the corresponding stabilities in terms of the nature of bonding and charge-transfer processes at the interface. We find ZrC(111)–W(110) is the most stable interface with higher adhesive energy than the other interfaces. The additional stability associated with the ZrC(111)–W(110) results from significant interface reconstruction. Three layers of W and ZrC adjacent to the interface are involved in the charge-transfer process leading to stronger ionic bonds in ZrC(111)–W(110) as compared to the other potential candidate: ZrC(100)–W(100). The C and W atoms are found to be displaced from their symmetric position during the reconstruction process at the interface to facilitate stronger bonds with shorter W–C and W–Zr bonds in ZrC(111)–W(110) as compared to ZrC(100)–W(100). This leads to stronger covalent bonds in ZrC(111)–W(110) than that in ZrC(100)–W(100). Therefore, we conclude that the stronger covalent and ionic forces in ZrC(111)–W(110) than those in ZrC(100)–W(100) are responsible for making ZrC(111)–W(110) to be the most stable interface. In conclusion, this study addresses the long-standing question of the most stable W–ZrC interface and derives a number of implications for other W-transition metal carbide interfaces which are potential candidates for improving the mechanical properties of plasma facing materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

DeePKS + ABACUS as a Bridge between Expensive Quantum Mechanical Models and Machine Learning Potentials

Recently, the development of machine learning (ML) potentials has made it possible to perform large-scale and long-time molecular simulations with the accuracy of quantum mechanical (QM) models. However, for different levels of QM methods, such as density functional theory (DFT) at the meta-GGA level and/or with exact exchange, quantum Monte Carlo, etc., generating a sufficient amount of data for training an ML potential has remained computationally challenging due to their high cost. In this work, we demonstrate that this issue can be largely alleviated with Deep Kohn–Sham (DeePKS), an ML-based DFT model. DeePKS employs a computationally efficient neural network-based functional model to construct a correction term added upon a cheap DFT model. Upon training, DeePKS offers closely matched energies and forces compared with high-level QM method, but the number of training data required is orders of magnitude less than that required for training a reliable ML potential. As such, DeePKS can serve as a bridge between expensive QM models and ML potentials: one can generate a decent amount of high-accuracy QM data to train a DeePKS model and then use the DeePKS model to label a much larger amount of configurations to train an ML potential. Further, this scheme for periodic systems is implemented in a DFT package ABACUS, which is open source and ready for use in various applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Huzinaga projection embedding for efficient and accurate energies of systems with localized spin-densities

We demonstrate the accuracy and efficiency of the restricted open-shell and unrestricted formulation of the absolutely localized Huzinaga projection operator embedding method. Restricted open-shell and unrestricted Huzinaga projection embedding in the full system basis is formally exact to restricted open-shell and unrestricted Kohn–Sham density functional theory, respectively. By utilizing the absolutely localized basis, we significantly improve the efficiency of the method while maintaining high accuracy. Furthermore, the absolutely localized basis allows for high accuracy open-shell wave function methods to be embedded into a closed-shell density functional theory environment. The open-shell embedding method is shown to calculate electronic energies of a variety of systems to within 1 kcal/mol accuracy of the full system wave function result. For certain highly localized reactions, such as spin transition energies on transition metals, we find that very few atoms are necessary to include in the wave function region in order to achieve the desired accuracy. This extension further broadens the applicability of our absolutely localized Huzinaga level-shift projection operator method to include open-shell species. Here, we apply our method to several representative examples, such as spin splitting energies, catalysis on transition metals, and radical reactions.

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

Role of Intermolecular Interactions in Deep Eutectic Solvents for CO 2 Capture: Vibrational Spectroscopy and Quantum Chemical Studies

Recent research and reviews on CO 2 capture methods, along with advancements in industry, have highlighted high costs and energy-intensive nature as the primary limitations of conventional direct air capture and storage (DACS) methods. In response to these challenges, deep eutectic solvents (DESs) have emerged as promising absorbents due to their scalability, selectivity, and lower environmental impact compared to other absorbents. However, the molecular origins of their enhanced thermal stability and selectivity for DAC applications have not been explored before. Therefore, the current study focuses on a comprehensive investigation into the molecular interactions within an alkaline DES composed of potassium hydroxide (KOH) and ethylene glycol (EG). Combining Fourier transform infrared (FT-IR) and quantum chemical calculations, the study reports structural changes and intermolecular interactions induced in EG upon addition of KOH and its implications on CO 2 capture. Experimental and computational spectroscopic studies confirm the presence of noncovalent interactions (hydrogen bonds) within both EG and the KOH-EG system and point to the aggregation of ions at higher KOH concentrations. Additionally, molecular electrostatic potential (MESP) surface analysis, natural bond orbital (NBO) analysis, quantum theory of atoms-in-molecules (QTAIM) analysis, and reduced density gradient-noncovalent interaction (RDG-NCI) plot analysis elucidate changes in polarizability, charge distribution, hydrogen bond types, noncovalent interactions, and interaction strengths, respectively. Evaluation of explicit and hybrid models assesses their effectiveness in representing intermolecular interactions. This research enhances our understanding of molecular interactions in the KOH-EG system, which are essential for both the absorption and desorption of CO 2 . The study also aids in predicting and selecting DES components, optimizing their ratios with salts, and fine-tuning the properties of similar solvents and salts for enhanced CO 2 capture efficiency.

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

Quantum-classical hybrid algorithm for the simulation of all-electron correlation

While chemical systems containing hundreds to thousands of electrons remain beyond the reach of quantum devices, hybrid quantum-classical algorithms present a promising pathway toward a quantum advantage. Hybrid algorithms treat the exponentially scaling part of the calculation-the static correlation-on the quantum computer and the non-exponentially scaling part-the dynamic correlation-on the classical computer. While a variety of algorithms have been proposed, the dependence of many methods on the total wave function limits the development of easy-to-use classical post-processing implementations. Here, we present a novel combination of quantum and classical algorithms, which computes the all-electron energy of a strongly correlated molecular system on the classical computer from the 2-electron reduced density matrix (2-RDM) evaluated on the quantum device. Significantly, we circumvent the wave function in the all-electron calculations by using density matrix methods that only require input of the statically correlated 2-RDM. Although the algorithm is completely general, we test it with two classical density matrix methods, the anti-Hermitian contracted Schrödinger equation (ACSE) and multiconfiguration pair-density functional theories, using the recently developed quantum ACSE method for simulating the statically correlated 2-RDM. Furthermore, we obtain experimental accuracy for the relative energies of all three benzyne isomers and thereby demonstrate the ability of the developed algorithm to achieve chemically relevant and accurate results on noisy intermediate-scale quantum devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1984

Effect of microcracks on the thermal expansion of composite laminates

A finite element analysis was used to quantitatively predict the effect of matrix microcracks in the 90 deg plies of graphite/epoxy composites on the coefficient of thermal expansion in the 0 deg direction, alpha (y) (perpendicular to the cracks). Results were generated for (0m/90n)s, (0/+ or 45/90)s and (0/90/ + or - 45)s graphite/epoxy laminate configurations. Analytical predictions were compared with experimental results for the two quasi-isotropic laminate configurations. Both analytical and experimental results showed that microcracks reduced the effective stiffness of the 90 deg plies, thus causing the laminates, thermal response to be more like that of a (0) laminate. The change in alpha(y) was a function of lamina material properties, microcrack density, fiber orientation, and stacking sequence. A combination of classical lamination theory and finite element analysis was used to predict the effect of microcracks in both the 90 deg and 0 deg plies. Analytical results showed that the addition of microcracks in the 0 deg plies do affect alpha(y), but to a lesser extent than those in the 90 deg plies.

Bowles, D. E.↗
DOE OSTI2021

Predicted Reaction Mechanisms, Product Speciation, Kinetics, and Detonation Properties of the Insensitive Explosive 2,6-Diamino-3,5-dinitropyrazine-1-oxide (LLM-105)

2,6-Diamino-3,5-dinitropyrazine-1-oxide (LLM-105) is a relatively new and promising insensitive high-explosive (IHE) material that remains only partially characterized. IHEs are of interest for a range of applications and from a fundamental science standpoint, as the root causes behind insensitivity are poorly understood. In this work, we adopt a multitheory approach based on reactive molecular dynamic simulations performed with density functional theory, density functional tight-binding, and reactive force fields to characterize the reaction pathways, product speciation, reaction kinetics, and detonation performance of LLM-105. We compare and contrast these predictions to 1,3,5-triamino-2,4,6-trinitrobenzene (TATB), a prototypical IHE, and 1,3,5,7-tetranitro-1,3,5,7-tetrazoctane (HMX), a more sensitive and higher performance material. The combination of different predictive models allows access to processes operative on progressively longer timescales while providing benchmarks for assessing uncertainties in the predictions. We find that the early reaction pathways of LLM-105 decomposition are extremely similar to TATB; they involve intra- and intermolecular hydrogen transfer. Additionally, the detonation performance of LLM-105 falls between that of TATB and HMX. We find agreement between predictive models for first-step reaction pathways but significant differences in final product formations. Predictions of detonation performance result in a wide range of values, and one-step kinetic parameters show the similar reaction rates at high temperatures for three out of four models considered.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Machine Learned Force Field Modeling of Metal Organic Frameworks for CO2 Direct Air Capture

Direct air capture (DAC) is a method for removing CO2 directly from air. Metal organic frameworks (MOFs) have been studied as DAC sorbent materials because of their structural and chemical diversity. Thermodynamic calculations using classical force fields are often used to evaluate MOFs for their performance in separations such as CO2 capture. Machine-learned force fields (MLFFs) can use machine learning to form quantitative relationships between a material’s chemical structure and the forces and energies predicted by more accurate quantum mechanical calculations, such as dispersion-corrected density functional theory (DFT). These descriptions of forces and energies can be used to improve the accuracy of adsorption calculations. In this work, classical models were used to pre-screen MOFs for CO2 capture. DFT calculations were then used to examine the adsorption mechanism. Next, MLFF models were developed for MOFs to achieve DFT-level accuracy for the forces and energies associated with MOF flexibility and CO2 adsorption. These methods were parametrized based on thousands of DFT calculations of CO2 in flexible MOFs and used to predict MOF structural properties as well as CO2 adsorption properties.

Findley, John↗
DOE OSTI2022

Rotating quantum turbulence in the unitary Fermi gas

Quantized vortices carry the angular momentum in rotating superfluids, and are key to the phenomenon of quantum turbulence. Advances in ultracold-atom technology enable quantum turbulence to be studied in regimes with both experimental and theoretical control, unlike the original contexts of superfluid helium experiments. While much work has been performed with bosonic systems, detailed studies of fermionic quantum turbulence are nascent, despite wide applicability to other contexts such as rotating neutron stars. In this paper, we present a large-scale study of quantum turbulence in rotating fermionic superfluids using an accurate time-dependent density functional theory called the superfluid local density approximation. We identify two different modes of turbulent decay in the dynamical equilibration of a rotating fermionic superfluid, and contrast these results with a computationally simpler description provided by the Gross-Pitaevskii equation, which we find can qualitatively reproduce these decay mechanisms if dissipation is explicitly included. These results demonstrate that dissipation mechanisms intrinsic to fermionic superfluids play a key role in differentiating fermionic from bosonic turbulence, which manifests by enhanced damping of Kelvin waves.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Electronic energy gap closure and metal-insulator transition in dense liquid hydrogen

Here, using quantum Monte Carlo (QMC) calculations, we investigate the insulator-metal transition observed in liquid hydrogen at high pressure. Below the critical temperature of the transition from the molecular to the atomic liquid, the fundamental electronic gap closure occurs abruptly, with a small discontinuity reflecting the weak first-order transition in the thermodynamic equation of state. Above the critical temperature, molecular dissociation sets in while the gap is still open. When the gap closes, the decay of the off-diagonal reduced density matrix shows that the liquid enters a gapless, but localized, phase: there is a crossover between the insulating and the metallic liquids. Compared to different density functional theory (DFT) functionals, our QMC calculations provide larger values for the fundamental gap and the electronic density of states close to the band edges, indicating that optical properties from DFT potentially benefit from error cancellations.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Thermodynamics and its prediction and CALPHAD modeling: Review, state of the art, and perspectives

Thermodynamics is a science concerning the state of a system, whether it is stable, metastable, or unstable, when interacting with its surroundings. The combined law of thermodynamics derived by Gibbs about 150 years ago laid the foundation of thermodynamics. In Gibbs combined law, the entropy production due to internal processes was not included, and the 2nd law was thus practically removed from the Gibbs combined law, so it is only applicable to systems under equilibrium, thus commonly termed as equilibrium or Gibbs thermodynamics. Gibbs further derived the classical statistical thermodynamics in terms of the probability of configurations in a system in the later 1800's and early 1900's. With the quantum mechanics (QM) developed in 1920's, the QM-based statistical thermodynamics was established and connected to classical statistical thermodynamics at the classical limit as shown by Landau in the 1940's. In 1960's the development of density functional theory (DFT) by Kohn and co-workers enabled the QM prediction of properties of the ground state of a system. On the other hand, the entropy production due to internal processes in non-equilibrium systems was studied separately by Onsager in 1930's and Prigogine and co-workers in the 1950's. In 1960's to 1970's the digitization of thermodynamics was developed by Kaufman in the framework of the CALculation of PHAse Diagrams (CALPHAD) modeling of individual phases with internal degrees of freedom. CALPHAD modeling of thermodynamics and atomic transport properties has enabled computational design of complex materials in the last 50 years. Our recently termed zentropy theory integrates DFT and statistical mechanics through the replacement of the internal energy of each individual configuration by its DFT-predicted free energy. The zentropy theory is capable of accurately predicting the free energy of individual phases, transition temperatures and properties of magnetic and ferroelectric materials with free energies of individual configurations solely from DFT-based calculations and without fitting parameters, and is being tested for other phenomena including superconductivity, quantum criticality, and black holes. Those predictions include the singularity at critical points with divergence of physical properties, negative thermal expansion, and the strongly correlated physics. Furthermore, those individual configurations may thus be considered as the genomic building blocks of individual phases in the spirit of the materials genome®. This has the potential to shift the paradigm of CALPHAD modeling from being heavily dependent on experimental inputs to becoming fully predictive with inputs solely from DFT-based calculations and machine learning models built on those calculations and existing experimental data through newly developed and future open-source tools. Furthermore, through the combined law of thermodynamics including the internal entropy production, it is shown that the kinetic coefficient matrix of independent internal processes is diagonal with respect to the conjugate potentials in the combined law, and the cross phenomena that the phenomenological Onsager flux and reciprocal relationships are due to the dependence of the conjugate potential of a molar quantity on nonconjugate molar quantities and other potentials, which can be predicted by the zentropy theory and CALPHAD modeling.

42 ENGINEERING↗
DOE OSTI2022

A 9.2-GHz clock transition in a Lu(II) molecular spin qubit arising from a 3,467-MHz hyperfine interaction

Spins in molecules are particularly attractive targets for next-generation quantum technologies, enabling chemically programmable qubits and potential for scale-up via self-assembly. Here, we report observation of one of the largest hyperfine interactions for a molecular system, A iso = 3467±50 MHz, along with an associated clock transition of unprecedented magnitude. This is achieved through chemical control of the degree of s-orbital mixing into the spin-bearing d-orbital associated with a series of spin-½ La(II) and Lu(II) complexes. Increased s-orbital character reduces spin-orbit coupling and enhances the electron-nuclear Fermi contact interaction. Both outcomes are advantageous for quantum applications: the former reduces spin-lattice relaxation, while the latter maximizes the hyperfine interaction that, in turn, generates a 9 gigahertz clock transition, leading to an increase in phase memory time from 1.0±0.4 to 12±1 microseconds for one of the Lu(II) complexes. Furthermore, these findings suggest strategies for development of molecular quantum technologies, akin to trapped ion systems.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Structural changes upon magnetic ordering in magnetocaloric AlFe 2 B 2

With a Curie temperature just above room temperature, AlFe 2 B 2 is a useful magnetocaloric material composed of earth-abundant elements. Here, we employ temperature-dependent high-resolution synchrotron X-ray diffraction to establish with high certainty that the paramagnetic to ferromagnetic transition in AlFe 2 B 2 is of second order, showing no discontinuity in lattice parameters or cell volume. Nevertheless, the lattice parameters undergo anisotropic changes across the transition with distinct differences in the thermal expansion coefficients. While the $a$ and $b$ lattice parameters show a positive thermal expansion, $c$ shows a negative thermal expansion. We link these changes to the respective interatomic distances to determine the contribution of magnetism to the anisotropic structural evolution. The work underpins the possible role of magnetostructural coupling in driving the magnetocaloric effect in AlFe 2 B 2 .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Observation of spatially resolved Rashba states on the surface of CH 3 NH 3 PbBr 3 single crystals

Hybrid organic-inorganic perovskites (HOIPs) are prime candidates for studying Rashba effects due to the heavy metal and halogen atoms in their crystal structure coupled with predicted inversion symmetry breaking. Nevertheless, observation of the Rashba effect in cubic CH 3 NH 3 PbBr 3 single crystals that possess bulk inversion symmetry is the subject of extensive debate due to the lack of conclusive experiments and theoretical explanations. In this work, we provide experimental evidence that Rashba state in cubic CH 3 NH 3 PbBr 3 single crystals at room temperature occurs exclusively on the crystal surface and depends on specific surface termination that results in local symmetry breaking. We demonstrate this using a suite of spatially resolved and depth-sensitive techniques, including circular photogalvanic effect, inverse spin Hall effect, and multiphoton microscopy, that are supported by first principle calculations. Our work suggests using surface Rashba states in these materials for spintronic applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Charge density wave activated excitons in TiSe 2 –MoSe 2 heterostructures

Layered materials enable the assembly of a new class of heterostructures where lattice-matching is no longer a requirement. Interfaces in these heterostructures therefore become a fertile ground for unexplored physics as dissimilar phenomena can be coupled via proximity effects. In this article, we identify an unexpected photoluminescence (PL) peak when MoSe 2 interacts with TiSe 2 . A series of temperature-dependent and spatially resolved PL measurements reveal that this peak is unique to the TiSe 2 –MoSe 2 interface, is higher in energy compared to the neutral exciton, and exhibits exciton-like characteristics. The feature disappears at the TiSe 2 charge density wave transition, suggesting that the density wave plays an important role in the formation of this new exciton. We present several plausible scenarios regarding the origin of this peak that individually capture some aspects of our observations but cannot fully explain this feature. These results therefore represent a fresh challenge for the theoretical community and provide a fascinating way to engineer excitons through interactions with charge density waves.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Photophysical properties of zero-dimensional perovskites studied by PBE0 and GW+BSE methods

Due to the exceptional photophysical properties and great potential applications in solid state lighting, zero-dimensional halide perovskites with broadband emitting have recently attracted much attention. In this work, to reveal the origin of the broadband absorption and emission spectra of 0D (C 4 N 2 H 14 Br) 4 SnBr x I 6-x (x = 3), the ground state and exciton dynamics properties have been studied by hybrid PBE0 functional and the GW+BSE method. Different self-trapped excitons (STEs) are obtained, and their emission energies by PBE0 are consistent well with experimental data. By considering the vibrational broadening of emission, our calculated full width at half maximum of emission accords well with the experiment. The absorption of free exciton and STEs is also calculated by considering the quasiparticle and excitonic effect. Our results show that these light-induced trapped states can be observed by experiments, as confirmed in other similar systems by the femtosecond transient absorption spectra.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Improved loss functions for machine-learned atomic potentials

Machine learning (ML) has become an invaluable tool across a wide array of domains in science as researchers find new ways to leverage its predictive power. This is especially true in chemistry, where ML is used to fit chemical properties or desirable attributes to the local structure of molecules and materials. In the pursuit of greater accuracy, it is relatively simple to increase the size or complexity of such models, although this often requires simultaneously seeking larger datasets in order to both fit and interpret the larger number of parameters. However, it is equally important to assess the quality and relative importance of the data and how these factors impact the training process. We, therefore, investigate the impact of using different loss functions for training neural network potentials (NNPs), as the loss function defines the error and parameter gradients used to train the NNP. In particular, we test the mean-squared error and Huber loss functions and, using insight from these functions, derive a new loss function based on the Asinh function, which yields significant improvement in the accuracy and generality of NNPs. We show that by discounting/minimizing errors and anomalies in the optimization process, both the Huber and Asinh loss functions improve the training of NNPs, leading to a final potential with a greater effective dimensionality.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗