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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Two-Dimensional Moiré Polaronic Electron Crystals

Two-dimensional moiré materials have emerged as the most versatile platform for realizing quantum phases of electrons. Here, we explore the stability origins of correlated states in WSe 2 /WS 2 moiré superlattices. We find that ultrafast electronic excitation leads to partial melting of the Mott states on timescales 5 times longer than predictions from the charge hopping integrals and that the melting rates are thermally activated, with activation energies of 1⁢8±3 and 1⁢3±2 meV for the one- and two-hole Mott states, respectively, suggesting significant electron-phonon coupling. A density functional theory calculation of the one-hole Mott state confirms polaron formation and yields a hole-polaron binding energy of 16 meV. These findings reveal a close interplay of electron-electron and electron-phonon interactions in stabilizing the polaronic Mott insulators at transition metal dichalcogenide moiré interfaces.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Exact-two-component block-localized wave function: A simple scheme for the automatic computation of relativistic ΔSCF

Block-localized wave function is a useful method for optimizing constrained determinants. In this article, we extend the generalized block-localized wave function technique to relativistic two-component framework. Optimization of excited state determinants for two-component wave functions presents a unique challenge because the excited state manifold is often quite dense with degenerate states. Furthermore, we test the degree to which certain symmetries result naturally from the ΔSCF optimization such as time reversal symmetry and symmetry with respect to the total angular momentum operator on a series of atomic systems. Variational optimizations may often break symmetry in order to lower the overall energy, just as unrestricted Hartree-Fock breaks spin symmetry. Overall, we demonstrate that time reversal symmetry is roughly maintained when using Hartree-Fock, but less so when using Kohn-Sham density functional theory. Additionally, maintaining total angular momentum symmetry appears to be system dependent and not guaranteed. Finally, we were able to trace the breaking of total angular momentum symmetry to the relaxation of core electrons.

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

A Comparison of Three Ways to Measure Time-Dependent Densities With Quantum Simulators

Quantum algorithms are touted as a way around some classically intractable problems such as the simulation of quantum mechanics. At the end of all quantum algorithms is a quantum measurement whereby classical data is extracted and utilized. In fact, many of the modern hybrid-classical approaches are essentially quantum measurements of states with short quantum circuit descriptions. Here, we compare and examine three methods of extracting the time-dependent one-particle probability density from a quantum simulation: direct Z -measurement, Bayesian phase estimation, and harmonic inversion. We have tested these methods in the context of the potential inversion problem of time-dependent density functional theory. Our test results suggest that direct measurement is the preferable method. We also highlight areas where the other two methods may be useful and report on tests using Rigetti's quantum virtual device. This study provides a starting point for imminent applications of quantum computing.

Yang, Jun↗
DOE OSTI2020

Effects of cluster expansion on the locations of phase transition boundary as a first step to quantify uncertainty in first principles statistical mechanics framework

Predicting phase diagrams from first principle calculations eliminates the need of tedious experimental trials and errors. Fully automating first principle phase diagram calculations without any sort of human intervention has been a long daunting task and troubling scientific communities for decades. This grand problem remains not fully resolved, largely due to the vastly high-dimensional parameter space associated with density functional theory, cluster expansion, lattice Monte Carlo, and the substantial uncertainty propagating through a set of complex simulations. As a first step to tackle this grand problem, we reported a first demonstration of how sensitive phase boundary locations can be to various cluster expansion fittings and input DFT training data. To the best knowledge of the authors, this study reported the first ever attempt to quantify uncertainty in first principles statistical mechanics framework. In addition, a semi-automated phase transition detection algorithm has been devised in this paper to deal with the associated statistical errors and uncertainties from Monte Carlo method and its predecessor DFT calculations and cluster expansions. This algorithm has been applied in a classical cluster expansion Mg-Cd binary alloy system to detect phase transitions at various locations of phase diagram using different cluster expansions to demonstrate its predictive power and quantify uncertainties. The results suggested that using Chebyshev basis function shifted the transition locations toward the dilute solid solution phase and expanded the phase stability of concentrated ordered phases to wider composition range. The addition of perturbed defect configurations into training data set lowered the order-disorder transition temperature to be away from true transition temperature, suggesting the transition nature is indeed configurational disorder dominated rather than defect assisted. Finally, we found that the weighting has negligible effect on the transition locations, except for the case of using Chebyshev basis function to fit all non-weighted configurations that can be susceptible to Monte Carlo sampling hysteresis.

36 MATERIALS SCIENCE↗
DOE OSTI2026

An Ab Initio Molecular Dynamics Study of Key Thermodynamic Input Parameters for Computer Simulation of U-6Nb Solidification

The key to metallic fuel development is the fabrication of uranium metal and alloys into fuel forms. U-Nb alloys are one of the best candidates for a metallic fuel alloy with high-temperature strength sufficient to support the core, acceptable nuclear properties, good fabricability, and compatibility with usable coolant media. Melt processing has been a key component of the metallic fuel cycle, and process models require thermophysical parameters at elevated temperatures, particularly above the melting temperatures, regarding which experimental data are scarce, for accurate simulations and process development. By means of ab initio density-functional theory (DFT) quantum molecular dynamics (QMD), we have calculated the main thermophysical parameters—the density, thermal expansion coefficient, specific heat, thermal conductivity, melting temperature, latent heat of fusion, and viscosity—used in the modeling of the U-6 wt.% Nb alloy casting. The melting temperature of the U-6 wt.% Nb alloy at ambient pressure is obtained by means of QMD simulations using the Z-method. The ambient volume change and latent heat of melting of U-6 wt.% Nb are also derived from QMD simulations in conjunction with analytical fitting for the energy and pressure. The thermal conductivity for the solid U-Nb alloy is calculated from the semi-classical Boltzmann transport equation combined with an estimate of the electron relaxation time obtained from DFT simulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Investigating the non-classical M-H 2 bonding in OsClH 3 (PPh 3 ) 3

Dihydrogen ligation to metal centers results in a wide range of metal-hydrogen interactions. Compounds containing non-classical dihydrogen ligands in which the H-H contacts are significantly longer than the 0.74 A distance found in solid H 2 are classified as Kubas-like compounds, stretched dihydrogen ligands, or compressed dihydrides, depending on the H-H contact length. Compounds containing stretched dihydrogen and compressed dihydride ligands are somewhat rare. Developing a more complete understanding of the metal-hydrogen and hydrogen-hydrogen interactions in these species may provide insights into the mechanisms for the oxidative addition reaction and hydrogen storage. Here, we use diffraction and inelastic neutron scattering (INS) measurements, paired with both molecular and solid-state density functional theory (DFT), to fully characterize the nature of the metal-hydrogen and hydrogen-hydrogen interactions in a previously reported compound, OsClH 3 (PPh 3 ) 3 . In that report, the H-H distance was measured as 1.48(2) A by single-crystal neutron diffraction and was described as a "stretched" dihydrogen. The INS-generated phonon density of states is well described by the DFT calculations, including those normal modes dominated by the vibrational motions of the three Os-ligated hydrogen atoms. Additionally, the resulting calculated electron charge density landscape indicates that there is no remnant bonding character between the Os-ligated hydrogen atoms in the compound. Based on these findings, we update the understanding of the metal-hydrogen bonding in this complex and reclassify it as a compressed dihydride. As such, OsClH 3 (PPh 3 ) 3 represents a rare example of a non-classical compressed dihydride species.

08 HYDROGEN↗
DOE OSTI2025

Multiscale Modeling Framework for Lithium Nucleation in 3D Porous Carbon Anodes

Porous carbon scaffolds offer a promising route for mitigating non-uniform lithium (Li) plating to enhance the safety and longevity of Li metal batteries. However, the influence of microstructural morphology on Li nucleation is not well understood. Here, we present a multiscale modeling framework to investigate how the porous microstructure of carbon materials affects Li nucleation behavior. Ab initio molecular dynamics simulations quantify the nucleation energy barriers of Li on graphene as a function of Li content, surface curvature, and applied potential, providing key parameters for a classical nucleation theory (CNT) model. From macroscale half-cell simulations, we obtained Li concentration and electrical potential profiles to define boundary conditions for mesoscopic simulations. At the mesoscale, three distinct synthetic 3D microstructures with different porosities and characteristic feature sizes are generated to resolve local distributions of Li flux, current density, and mechanical stress. These outputs are integrated into the CNT model to map spatial variation in nucleation rates. Our findings reveal trade-offs between suppressing nucleation rates and achieving spatial uniformity, offering design guidelines for optimizing porous carbon anodes to balance nucleation control and mechanical integrity.

Materials science↗
DOE OSTI2020

Crystal and magnetic structures of magnetic topological insulators MnBi 2 Te 4 and MnBi 4 Te 7

Using single-crystal neutron diffraction, we present a systematic investigation of the crystal structure and magnetism of the van der Waals topological insulators MnBi 2 Te 4 and MnBi 4 Te 7 , where rich topological quantum states have been recently predicted and observed. Structural refinements reveal that considerable Bi atoms occupied on the Mn sites in both materials, distinct from the previously reported antisite disorder. We show unambiguously that MnBi 2 Te 4 orders antiferromagnetically below 24 K, featured by a magnetic symmetry R I - 3 c , while MnBi 4 Te 7 is antiferromagnetic below 13 K with a magnetic space group P c - 3 c 1 . They both present antiferromagnetically coupled ferromagnetic layers with spins along the c axis. We put forward a stacking rule for the crystal structure of an infinitely adaptive series MnBi 2 n Te 3 n + 1 ( n ≥ 1 ) with a building unit of [ Bi 2 Te 3 ] . Overall, by comparing the magnetic properties between MnBi 2 Te 4 and MnBi 4 Te 7 , together with recent density-functional theory calculations, we concluded that a two-dimensional magnetism limit might be realized in the derivatives. Our work may promote theoretical studies of topological magnetic states in the series of MnBi 2 n Te 3 n + 1 .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Crystal structure of valbenazine, C 24 H 38 N 2 O 4

The crystal structure of valbenazine has been solved and refined using synchrotron X-ray powder diffraction data and optimized using density functional theory techniques. Valbenazine crystallizes in space groupP2 1 2 1 2 1 (#19) witha= 5.260267(17),b= 17.77028(7),c= 26.16427(9) Å,V= 2445.742(11) Å 3 , andZ= 4 at 295 K. The crystal structure consists of discrete molecules and the mean plane of the molecules is approximately (8,−2,15). There are no obvious strong intermolecular interactions. There is only one weak classical hydrogen bond in the structure, from the amino group to the ether oxygen atom. Two intramolecular and one intermolecular C–H⋯O hydrogen bonds also contribute to the lattice energy. The powder pattern has been submitted to ICDD for inclusion in the Powder Diffraction File™ (PDF®)

Materials Science↗
DOE OSTI2022

Atomic Molybdenum for Synthesis of Ammonia with 50% Faradic Efficiency

The electrochemical dinitrogen (N 2 ) reduction reaction (NRR) under ambient conditions has gained significant interest as an environmentally friendly alternative to the traditional Haber–Bosch process for the synthesis of ammonia (NH 3 ). However, up to now, most of the reported NRR electrocatalysts with satisfactory catalytic activities have been hindered by the large overpotential in N 2 activation. The preparation of highly efficient Mo-based NRR electrocatalyst in acidic electrolytes under ambient conditions is demonstrated here, consisting of stabilized single Mo atoms anchored on holey nitrogen-doped graphene synthesized through a convenient potassium-salt-assisted activation method. At -0.05 V versus a reversible hydrogen electrode (RHE), an electrode consisting of the resultant electrocatalyst immobilized on carbon fiber paper can attain an exceptional Faradaic efficiency of 50.2% and a NH3 yield rate of 3.6 µg h -1 mg cat -1 with low overpotentials. Density functional theory calculations further unveil that compared to the original graphene without holes, the edge coordinated Mo atoms and the existence of vacancies on holey graphene lower the overpotential of N 2 reduction, thereby promoting the NRR catalytic activity. Finally, this work could provide new guidelines for future designs in single-atom catalysis that would be beneficial to ambient N 2 fixation, and replacement of classical synthesis processes that are very energy-intensive.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A Transferable Force Field for Simulating Adsorption in Metal–Organic Frameworks with Open Metal Sites Based on the 12–6–4 Lennard-Jones Potential

Metal−organic frameworks (MOFs) that contain coordinatively unsaturated open metal sites (OMSs) provide strong host− guest interactions, making them promising sorbents for low-concentration gas adsorption applications such as direct air capture and atmospheric water harvesting. However, accurately modeling host−guest interactions involving OMSs remains challenging for classical force fields (FFs) based on the 12−6 Lennard−Jones (LJ) potential, as the polarization effect of the guest molecule induced by the positively charged OMS is not considered. Here, we introduce an FF based on the 12−6−4 LJ potential, which incorporates charge−induced dipole interactions and is parametrized against a diverse set of host−guest potential energy surfaces (PESs) obtained from density functional theory (DFT). The resulting FF, trained on a generic trimetallic cluster, performs well in both host−guest binding energetics and gas adsorption isotherms across different OMS-containing MOFs, including MOF-74 series and Cu-BTC. These results highlight the excellent transferability of our approach and its potential to enhance the accuracy and robustness of high-throughput MOF discovery workflows, particularly for gas adsorption and separation in large and diverse MOF databases.

36 MATERIALS SCIENCE↗
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 OSTI2022

Improving the Accuracy of Atomistic Simulations of the Electrochemical Interface

Atomistic simulation of the electrochemical double layer is an ambitious undertaking, requiring quantum mechanical description of electrons, phase space sampling of liquid electrolytes, and equilibration of electrolytes over nanosecond time scales. All models of electrochemistry make different trade-offs in the approximation of electrons and atomic configurations, from the extremes of classical molecular dynamics of a complete interface with point-charge atoms to correlated electronic structure methods of a single electrode configuration with no dynamics or electrolyte. Here, we review the spectrum of simulation techniques suitable for electrochemistry, focusing on the key approximations and accuracy considerations for each technique. Additionally, we discuss promising approaches, such as enhanced sampling techniques for atomic configurations and computationally efficient beyond density functional theory (DFT) electronic methods, that will push electrochemical simulations beyond the present frontier.

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

Parameterizing empirical interatomic potentials for predicting thermophysical properties via an irreducible derivative approach: the case of ThO 2 and UO 2

The accuracy of classical physical property predictions using molecular dynamics simulations is determined by the quality of the interatomic potentials. Here we introduce a training approach for empirical interatomic potentials (EIPs) which is well suited for capturing phonons and phonon-related properties. Our approach is based on direct comparisons of the second- and third-order irreducible derivatives (IDs) between an EIP and the Born–Oppenheimer potential within density functional theory (DFT) calculations. IDs fully exploit space group symmetry and allow for training without redundant information. We demonstrate the fidelity of our approach in the context of ThO 2 and UO 2 , where we optimize parameters of an embedded-atom method potential in addition to core–shell interactions. Our EIPs provide thermophysical properties in good agreement with DFT and outperform widely utilized EIPs for phonon dispersion and thermal conductivity predictions. Reasonable estimates of thermal expansion and formation energies of Frenkel pairs are also obtained.

empirical interatomic potential↗
DOE OSTI2021

Collective excitations and quantum incompressibility in electron-hole bilayers

In this study, we apply quantum continuum mechanics to the calculation of the excitation spectrum of a coupled electron-hole bilayer. The theory expresses excitation energies in terms of ground-state intra- and interlayer pair correlation functions, which are available from quantum Monte Carlo calculations. The final formulas for the collective modes deduced from this approach coincide with the formulas obtained in the “quasilocalized particle approximation” by Kalman et al. [G. Kalman, V. Valtchinov, and K. I. Golden, Phys. Rev. Lett. 82 , 3124 (1999) ], and likewise the theory predicts the existence of gapped excitations in the charged channel, with the gap arising from electron-hole correlation. An immediate consequence of the gap is that the static density-density response function of the charged channel vanishes as q 2 for wave vector q → 0 , rather than linearly in q , as commonly expected. In this sense, the system is incompressible . This feature, which has no analog in the classical electron-hole plasma, is consistent with the existence of an excitonic ground state and implies the existence of a discontinuity in the chemical potential of electrons and holes when the numbers of electrons and holes are equal. It should be experimentally observable by monitoring the densities of electrons and holes in response to potentials that attempt to change these densities in opposite directions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Evolution of the Kondo lattice electronic structure above the transport coherence temperature

Significance The temperature ( T )-dependent evolution of the Kondo lattice electronic structure is a long-standing topic of theoretical and experimental investigation, still lacking a truly microscopic theory that agrees with a full experimental characterization. Here multiple characteristic T scales of the interaction of localized f moments with conduction electrons in the Kondo lattice CeCo I n 5 are identified and investigated using angle-resolved photoemission measurements that substantiate dynamical mean-field theory which newly includes the full realism of crystalline electric-field (CEF) f splittings. Thereby errors in the itinerant versus localized f -state classification from standard density functional theory are corrected, microscopic insight into the broad T -range crossover of f -hybridization effects is gained, and a prediction of CEF degeneracy crossover below the lattice coherence T is made.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Reaction nanoscopy of ion emission from sub-wavelength propanediol droplets

Abstract Droplets provide unique opportunities for the investigation of laser-induced surface chemistry. Chemical reactions on the surface of charged droplets are ubiquitous in nature and can provide critical insight into more efficient processes for industrial chemical production. Here, we demonstrate the application of the reaction nanoscopy technique to strong-field ionized nanodroplets of propanediol (PDO). The technique’s sensitivity to the near-field around the droplet allows for the in-situ characterization of the average droplet size and charge. The use of ultrashort laser pulses enables control of the amount of surface charge by the laser intensity. Moreover, we demonstrate the surface chemical sensitivity of reaction nanoscopy by comparing droplets of the isomers 1,2-PDO and 1,3-PDO in their ion emission and fragmentation channels. Referencing the ion yields to gas-phase data, we find an enhanced production of methyl cations from droplets of the 1,2-PDO isomer. Density functional theory simulations support that this enhancement is due to the alignment of 1,2-PDO molecules on the surface. The results pave the way towards spatio-temporal observations of charge dynamics and surface reactions on droplets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Ionization state and dielectric constant in cold rarefied hydrocarbon plasmas of inertial confinement fusion

A combined approach to study cold rarefied matter is introduced that includes a semianalytical method based on the free-energy minimization and ab initio calculations based on the finite-temperature density-functional theory. The approach is used to calculate the ionization state of hydrocarbon (CH) under the shock-release conditions in inertial confinement fusion. Furthermore, the dielectric constant of CH is calculated using the Kubo-Greenwood formulation and contribution from atomic polarizabilities is found to be as important as the free-electron contribution. Using the ionization state and dielectric constant, the electron density profile in the rarefaction wave of the shock-release plasma is obtained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗