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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A family of ionic supersalts with covalent-like directionality and unconventional multiferroicity

Ionic crystals composed of elemental ions such as NaCl are non-polar due to directionless ionic bonding interactions. Here, we show that these can develop polarity by changing their building blocks from elemental ions to superalkali and superhalogen cluster-ions, which mimic the chemistry of alkali and halogen atoms, respectively. Due to the non-spherical geometries of these cluster ions, corresponding supersalts form anisotropic polar structures with ionic bonding, yet covalent-like directionality, akin to sp 3 hybridized systems. Using density functional theory and extensive structure searches, we predict a series of stable ferroelectric/ferroelastic supersalts, PnH 4 MX 4 (Pn = N, P; M = B, Al, Fe; X = Cl, Br) composed of superalkali PnH 4 and superhalogen MX 4 ions. Unlike traditional ferroelectric/ferroelastic materials, the cluster-ion based supersalts possess ultra-low switching barrier and can endure large ion displacements and reversible strain. In particular, PH 4 FeBr 4 exhibits triferroic coupling of ferroelectricity, ferroelasticity, and antiferromagnetism with controllable spin directions via either ferroelastic or 90-degree ferroelectric switching.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Accurate numerical simulations of open quantum systems using spectral tensor trains

Decoherence between qubits is a major bottleneck in quantum computations. Decoherence results from intrinsic quantum and thermal fluctuations as well as noise in the external fields that perform the measurement and preparation processes. With prescribed colored noise spectra for intrinsic and extrinsic noise, we present a numerical method, Quantum Accelerated Stochastic Propagator Evaluation (Q-ASPEN), to solve the time-dependent noise-averaged reduced density matrix in the presence of intrinsic and extrinsic noise. Q-ASPEN is arbitrarily accurate and can be applied to provide estimates for the resources needed to error-correct quantum computations. We employ spectral tensor trains, which combine the advantages of tensor networks and pseudospectral methods, as a variational ansatz to the quantum relaxation problem and optimize the ansatz using methods typically used to train neural networks. Here, the spectral tensor trains in Q-ASPEN make accurate calculations with tens of quantum levels feasible. We present benchmarks for Q-ASPEN on the spin-boson model in the presence of intrinsic noise and on a quantum chain of up to 32 sites in the presence of extrinsic noise. In our benchmark, the memory cost of Q-ASPEN scales as a low-order polynomial in the size of the system once the number of system states surpasses the number of basis functions used in the spectral expansion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Monolayer CrCl 3 as an Ideal Test Bed for the Universality Classes of 2D Magnetism

The monolayer halides CrX 3 (X=Cl, Br, I) attract significant attention for realizing 2D magnets with genuine long-range order (LRO), challenging the Mermin-Wagner theorem. In this work, we show that monolayer CrCl 3 has the unique benefit of exhibiting tunable magnetic anisotropy upon applying a compressive strain. This opens the possibility to use CrCl 3 for producing and studying both ferromagnetic and antiferromagnetic 2D Ising-type LRO as well as the Berezinskii-Kosterlitz-Thouless (BKT) regime of 2D magnetism with quasi-LRO. Using state-of-the-art density functional theory, we explain how realistic compressive strain could be used to tune the monolayer's magnetic properties so that it could exhibit any of these phases. Building on large-scale quantum Monte Carlo simulations, we compute the phase diagram of strained CrCl 3 , as well as the magnon spectrum with spin-wave theory. Our results highlight the eminent suitability of monolayer CrCl 3 to achieve very high BKT transition temperatures, around 50 K, due to their singular dependence on the weak easy-plane anisotropy of the material.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Spectral-density estimation with the Gaussian integral transform

The spectral-density operator $\hat{ρ}(ω) = δ(ω–\hat{H})$ plays a central role in linear response theory as its expectation value, the dynamical response function, can be used to compute scattering cross sections. In this work, we describe a near optimal quantum algorithm providing an approximation to the spectral density with energy resolution $\Delta$ and error $\epsilon$ using $O(\sqrt{\text{log}_2 (1/ε)[\text{log}_2 (1 / Δ) + \text{log}_2 (1/ε)]/ Δ)}$ operations. This is achieved without using expensive approximations to the time-evolution operator, but instead exploiting qubitization to implement an approximate Gaussian integral transform of the spectral density. Finally, we also describe appropriate error metrics to assess the quality of the spectral function approximations more generally.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Unconventional surface state pairs in a high-symmetry lattice with anti-ferromagnetic band-folding

Many complex magnetic structures in a high-symmetry lattice can arise from a superposition of well-defined magnetic wave vectors. These “multi-q” structures have garnered much attention because of interesting real-space spin textures such as skyrmions. However, the role multi-q structures play in the topology of electronic bands in momentum space has remained rather elusive. Here we show that the type-I anti-ferromagnetic 1q, 2q and 3q structures in an face-centered cubic sublattice with band inversion, such as NdBi, can induce unconventional surface state pairs inside the band-folding hybridization bulk gap. Our density functional theory calculations match well with the recent experimental observation of unconventional surface states with hole Fermi arc-like features and electron pockets below the Neel temperature. We further show that these multi-q structures have Dirac and Weyl nodes. Our work reveals the special role that band-folding from anti-ferromagnetism and multi-q structures can play in developing new types of surface states.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Toward Addressing the Challenge to Predict the Heat Capacities of RDX and HMX Energetic Materials

Availability of heat capacity as function of pressure and temperature is an essential prerequisite for development of a computational multiscale strategy capable to address the evolution of microstructure and energy release in advanced high energy density materials. In the case of 1,3,5-trinitro-1,3,5-triazinane (RDX) and 1,3,5,7-tetranitro-1,3,5,7-tetrazocane (HMX) systems as two of the most studied energetic materials, there are substantial gaps in experimental data, with available heat capacities values distributed only in a region close to standard ambient conditions. In this study we demonstrate how these major experimental limitations can be addressed in the case of the RDX and HMX systems based on the combined use of classical and quantum mechanical calculations. We show that by considering ideal gas properties evaluated using quantum mechanical methods, and residual properties obtained from molecular simulations using fully flexible atomistic force field models, excellent agreement can be obtained for the predicted heat capacities to the most recent experimental values. An important advantage of the current computational methodology is that it allows evaluation of both constant-volume and constant-pressure heat capacities for a broad interval of temperatures and pressures, which encompasses solid and liquid phases conditions. In the case of the solid α and γ phases of RDX and the β phase of HMX, the predicted results follow closely both the available experimental data at standard ambient conditions and the results obtained using density functional theory calculations at high pressures, a regime where experimental data are not available. A perspective to expand the current methodology is also discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Impact of classical statistics on thermal conductivity predictions of BAs and diamond using machine learning molecular dynamics

Machine learning interatomic potentials (MLIPs) have greatly enhanced molecular dynamics (MD) simulations, achieving near-first-principles accuracy in thermal conductivity studies. In this work, we reveal that this accuracy, observed in BAs and diamond at sub-Debye temperatures, stems from an accidental error cancelation: classical statistics overestimates specific heat while underestimating phonon lifetimes, balancing out in thermal conductivity predictions. However, this balance is disrupted when isotopes are introduced, leading MLIP-based MD to significantly underpredict thermal conductivity compared to experiments and quantum statistics-based Boltzmann transport equation. This discrepancy arises not from classical statistics affecting phonon–isotope scattering rates but from its impact on the interplay between phonon–isotope and phonon–phonon scattering in the normal scattering-dominated BAs and diamond. In conclusion, this work underscores the limitations of MLIP-based MD for thermal conductivity studies at sub-Debye temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Tuning Active Site Flexibility by Defect Engineering of Graphene Ribbon Edge–hosted Fe–N3 Sites

Nitrogen-doped, carbon-supported transition metal catalysts are excellent for several reactions. Structural engineering of metal-Nx sites to boost catalytic activity is rarely studied. Here, we demonstrate that the structural flexibility of Fe-N3 site is vital for tuning the electronic structure of Fe atoms and regulating the catalytic transfer hydrogenation (CTH) activity. By introducing carbon defects, we construct Fe-N3 sites with varying Fe-N bond lengths distinguishable by X-ray absorption spectroscopy. We investigate the CTH activity by density functional theory and microkinetic calculations and reveal that the vertical displacement of the Fe atom out of the plane of the support, induced by the Fe-N3 distortion, raises the Fe [[EQUATION]]orbital and strengthens binding. We propose that the activity is controlled by the relaxation of the reconstructed site, which is further affected by Fe-N bond length, an excellent activity descriptor. Furthermore, we elucidate the origin of the CTH activity and principles for high-performing Fe-N-C catalysts by defect engineering.

09 BIOMASS FUELS↗
DOE OSTI2023

Grey Galaxies’ as an endpoint of the Kerr-AdS superradiant instability

Kerr-AdS d+1 black holes for d ≥ 3 suffer from classical superradiant instabilities over a range of masses above extremality. We conjecture that these instabilities settle down into Grey Galaxies (GGs) — a new class of coarse-grained solutions to Einstein’s equations which we construct in d = 3. Grey Galaxies are made up of a black hole with critical angular velocity ω = 1 in the ‘centre’ of AdS, surrounded by a large flat disk of thermal bulk gas that revolves around the centre of AdS at the speed of light. The gas carries a finite fraction of the total energy, as its parametrically low energy density and large radius are inversely related. GGs exist at masses that extend all the way down to the unitarity bound. Their thermodynamics is that of a weakly interacting mix of Kerr-AdS black holes and the bulk gas. Their boundary stress tensor is the sum of a smooth ‘black hole’ contribution and a peaked gas contribution that is delta function localized around the equator of the boundary sphere in the large N limit. We also construct another class of solutions with the same charges; ‘Revolving Black Holes (RBHs)’. RBHs are macroscopically charged SO(d, 2) descendants of AdS-Kerr solutions, and consist of ω = 1 black holes revolving around the centre of AdS at a fixed radial location but in a quantum wave function in the angular directions. RBH solutions are marginally entropically subdominant to GG solutions and do not constitute the endpoint of the superradiant instability. Nonetheless, we argue that supersymmetric versions of these solutions have interesting implications for the spectrum of supersymmetric states in, e.g. N = 4 Yang-Mills theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Cosmological decoherence from thermal gravitons

A bstract We study the effects of gravitationally-driven decoherence on tunneling processes associated with false vacuum decays, such as the Coleman-De Luccia instanton. We compute the thermal graviton-induced decoherence rate for a wave function describing a perfect fluid of nonzero energy density in a finite region. When the effective cosmological constant is positive, the thermal graviton background sourced by a de Sitter horizon provides an unavoidable decoherence effect, which may have important consequences for tunneling processes in cosmological history. We discuss generalizations and consequences of this effect and comment on its observability and applications to black hole physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

CuAu binary alloy with 32 atoms - LSMS-3 data

This dataset contains the estimate of atomic charge density, atomic magnetic moment and total energy for 32,000 configurations of the iron-platinum (CuAu) binary alloy with face-centered cubic (FCC) structure. The configurations span all the compositions from 0%Cu - 100% Au through 100%Cu - 0% Au. The results have been produced running ab-initio density functional theory (DFT) calculations with the LSMS-3 code on OLCF supercomputer Titan. LSMS-3 GitHub repository: https://github.com/mstsuite/lsms Deep Learning research papers published with results based on this dataset: Fast and stable deep-learning predictions of material properties for solid solution alloys Massimiliano Lupo Pasini, Ying Wai Li, Junqi Yin, Jiaxin Zhang, Kipton Barros and Markus Eisenbach Published 14 December 2020. 2020 IOP Publishing Ltd Journal of Physics: Condensed Matter, Volume 33, Number 8 Citation Massimiliano Lupo Pasini et al 2021 J. Phys.: Condens. Matter 33 084005 https://iopscience.iop.org/article/10.1088/1361-648X/abcb10

36 MATERIALS SCIENCE↗
DOE OSTI2022

Toward a QUBO-Based Density Matrix Electronic Structure Method

Density matrix electronic structure theory is used in many quantum chemistry methods to “alleviate” the computational cost that arises from directly using wave functions. Although density matrix based methods are computationally more efficient than wave function based methods, significant computational effort is involved. Because the Schrödinger equation needs to be solved as an eigenvalue problem, the time-to-solution scales cubically with the system size in mean-field type approaches such as Hartree–Fock and density functional theory and is solved as many times in order to reach charge or field self-consistency. We hereby propose and study a method to compute the density matrix by using a quadratic unconstrained binary optimization (QUBO) solver. This method could be useful to solve the problem with quantum computers and, more specifically, quantum annealers. Our proposed approach is based on a direct construction of the density matrix using a QUBO eigensolver. We explore the main parameters of the algorithm focusing on precision and efficiency. We show that, while direct construction of the density matrix using a QUBO formulation is possible, the efficiency and precision have room for improvement. Moreover, calculations performed with quantum annealing on D-Wave’s new Advantage quantum computer are compared with results obtained with classical simulated annealing, further highlighting some problems of the proposed method. Finally, we also suggest alternative methods that could lead to a more efficient QUBO-based density matrix construction.

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

Predicting Hugoniot equation of state in erythritol with ab initio and reactive molecular dynamics

Erythritol has been proposed as an inert surrogate for developing theoretical and computational models to study aging in energetic materials. In this work, we present a comparison of mechanical and shock properties of erythritol computed using the ReaxFF reactive force field and from ab initio calculations employing density functional theory (DFT). We screened eight different ReaxFF parameterizations, of which the CHO parameters developed for hydrocarbon oxidation provide the most accurate predictions of mechanical properties and the crystal structure of erythritol. Further validation of the applicability of this ReaxFF parameterization for modeling erythritol is demonstrated by comparing predictions of the elastic constants, crystal structure, vibrational density of states, and Hugoniot curves against DFT calculations. The ReaxFF predictions are in close agreement with the DFT simulations for the elastic constants and shock Hugoniot when the crystal is loaded along its c axis but show as much as 30% disagreement in the elastic constants in the ab plane and 12% difference in shock pressures when shocked along the a or b crystal axes. Last, we compare thermomechanical properties predicted from classical molecular dynamics with those calculated using the quasi-harmonic approximation and show that quantum mechanical effects produce large discrepancies in the computed values of heat capacity and thermal expansion coefficients compared with classical assumptions. Combining classical molecular dynamics predictions of mechanical behavior with phonon-based calculations of thermal behaviors, we show that predicted shock-induced temperatures for pressures up to 6.5 GPa do not exceed the pressure-dependent melting point of erythritol.

Hu, Jing↗
DOE OSTI2025

Spin-informed universal graph neural networks for simulating magnetic ordering

The screening and discovery of magnetic materials are hindered by the computational cost of first-principles density-functional theory (DFT) calculations required to find the ground state magnetic ordering. Although universal machine-learning interatomic potentials (uMLIPs), also known as atomistic foundation models, offer high-fidelity models of many atomistic systems with significant speedup, they currently lack the inputs required for predicting magnetic ordering. In this work, we present a data-efficient, spin-informed graph neural network framework that incorporates spin degrees of freedom as inputs and preserves physical symmetries, extending the functionality of uMLIPs to simulate magnetic orderings. This framework speeds up DFT calculations through better initial guesses for magnetic moments, determines the ground-state ordering of bulk materials and even generalizes to magnetic ordering in surfaces. Furthermore, we implement a closed-loop anomaly detection approach that effectively addresses the classic "chicken-and-egg" problem of creating a high-quality dataset while developing a uMLIP, unearthing anomalies in large benchmark datasets and boosting model accuracy.

Xu, Wenbin↗
DOE OSTI2021

$\mathrm{Perturbo}$: A software package for ab initio electron–phonon interactions, charge transport and ultrafast dynamics

We report Perturbo is a software package for first-principles calculations of charge transport and ultrafast carrier dynamics in materials. The current version focuses on electron–phonon interactions and can compute phonon-limited transport properties such as the conductivity, carrier mobility and Seebeck coefficient. It can also simulate the ultrafast nonequilibrium electron dynamics in the presence of electron–phonon scattering. Perturbo uses results from density functional theory and density functional perturbation theory calculations as input, and employs Wannier interpolation to reduce the computational cost. It supports norm-conserving and ultrasoft pseudopotentials, spin–orbit coupling, and polar electron–phonon interactions for bulk and 2D materials. Hybrid MPI plus OpenMP parallelization is implemented to enable efficient calculations on large systems (up to at least 50 atoms) using high-performance computing. Taken together, Perturbo provides efficient and broadly applicable ab initio tools to investigate electron–phonon interactions and carrier dynamics quantitatively in metals, semiconductors, insulators, and 2D materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Coexisting structural disorder and robust spin-polarization in half-metallic FeMnVAl

We report half-metallic ferromagnets (HMFs) are one of the most promising materials in the field of spintronics due to their unique band structure, consisting of one spin subband having metallic characteristics along with another subband with semiconductor-like behavior. In this work, we report the synthesis of a novel quaternary Heusler alloy—FeMnVAl—and we have studied the structural, magnetic, transport, and electronic properties complemented with first-principles calculations. Among different possible structurally ordered arrangements, the optimal structure is identified by theoretical energy minimization. The corresponding spin-polarized band-structure calculations indicate the presence of a half-metallic ferromagnetic ground state. A detailed and careful investigation of the x-ray diffraction data, Mössbauer spectra, and nuclear magnetic resonance spectra suggests the presence of site disorder between the Fe and Mn atoms in the stable ordered structure of the system. The magnetic susceptibility measurement clearly establishes a ferromagnetic-like transition below ~213 K. The Fe Mössbauer spectrometry measurements suggest that only the Mn spins could be responsible for the magnetic order, which is consistent with our theoretical calculation. Surprisingly, the density-functional-theory calculations reveal that the spin-polarization value is almost immunized (92.4% → 90.4%) from the Mn-Fe structural disorder, even when nonmagnetic Fe and moment-carrying Mn sites are entangled inseparably. The robustness of spin polarization and half-metallicity in the studied FeMnVAl compound comprising structural disorder is thus quite interesting and could provide a new direction to investigate and understand the exact role of disorders on spin polarization in this class of materials, beyond the available knowledge.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Crystal structure of ractopamine hydrochloride, C 18 H 24 NO 3 Cl

The crystal structure of ractopamine hydrochloride has been solved and refined using synchrotron X-ray powder diffraction data, and optimized using density functional theory techniques. Ractopamine hydrochloride crystallizes in space groupPbca(#61) witha= 38.5871(49),b= 10.7691(3),c= 8.4003(2) Å,V= 3490.75(41) Å 3 , andZ= 8. The ractopamine cation contains two chiral centers, and the sample consists of a mixture of the S,S/R,R/S,R and R,S forms. Models for the two diastereomers S,S and S,R were refined, and yielded equivalent residuals, but the S,R form is significantly lower in energy. The crystal structure consists of layers of molecules parallel to thebc-plane. In each structure one of the H atoms on the protonated N atom acts as a donor in a strong discrete N–H⋯Cl hydrogen bond. Hydroxyl groups act as donors in O–H⋯Cl and O–H⋯O hydrogen bonds. Both the classical and C–H⋯Cl and C–H⋯O hydrogen bonds differ between the forms, helping to explain the large microstrain observed for the sample. The powder pattern has been submitted to ICDD® for inclusion in the Powder Diffraction File™ (PDF®).

Materials Science↗
DOE OSTI2021

\( \mathcal{N} \) = 4 supersymmetric Yang-Mills thermodynamics to order λ 2

We calculate the resummed perturbative free energy of \( \mathcal{N} \) = 4 supersymmetric Yang-Mills in four spacetime dimensions (SYM 4 , 4 ) through second order in the ’t Hooft coupling λ at finite temperature and zero chemical potential. Our final result is ultraviolet finite and all infrared divergences generated at three-loop level are canceled by summing over SYM 4 , 4 ring diagrams. Non-analytic terms at \( \mathcal{O} \) ( λ 3 / 2 ) and \( \mathcal{O} \) ( λ 2 log λ ) are generated by dressing the A 0 and scalar propagators. The gauge-field Debye mass m D and the scalar thermal mass M D are determined from their corresponding finite-temperature self-energies. Based on this, we obtain the three-loop thermodynamic functions of SYM 4 , 4 to \( \mathcal{O} \) ( λ 2 ). We compare our final result with prior results obtained in the weak- and strong-coupling limits and construct a generalized Padé approximant that interpolates between the weak-coupling result and the large- N c strong-coupling result. Our results suggest that the \( \mathcal{O} \) ( λ 2 ) weak-coupling result for the scaled entropy density is a quantitatively reliable approximation to the scaled entropy density for 0 ≤ λ ≲ 2.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗