Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “strongly correlated systems 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 19 records

Electronic structure of strongly correlated systems: recent developments in multiconfiguration pair-density functional theory and multiconfiguration nonclassical-energy functional theory

Strong electron correlation plays an important role in transition-metal and heavy-metal chemistry, magnetic molecules, bond breaking, biradicals, excited states, and many functional materials, but it provides a significant challenge for modern electronic structure theory. The treatment of strongly correlated systems usually requires a multireference method to adequately describe spin densities and near-degeneracy correlation. However, quantitative computation of dynamic correlation with multireference wave functions is often difficult or impractical. Multiconfiguration pair-density functional theory (MC-PDFT) provides a way to blend multiconfiguration wave function theory and density functional theory to quantitatively treat both near-degeneracy correlation and dynamic correlation in strongly correlated systems; it is more affordable than multireference perturbation theory, multireference configuration interaction, or multireference coupled cluster theory and more accurate for many properties than Kohn–Sham density functional theory. This perspective article provides a brief introduction to strongly correlated systems and previously reviewed progress on MC-PDFT followed by a discussion of several recent developments and applications of MC-PDFT and related methods, including localized-active-space MC-PDFT, generalized active-space MC-PDFT, density-matrix-renormalization-group MC-PDFT, hybrid MC-PDFT, multistate MC-PDFT, spin–orbit coupling, analytic gradients, and dipole moments. We also review the more recently introduced multiconfiguration nonclassical-energy functional theory (MC-NEFT), which is like MC-PDFT but allows for other ingredients in the nonclassical-energy functional. We discuss two new kinds of MC-NEFT methods, namely multiconfiguration density coherence functional theory and machine-learned functionals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Interwoven atypical quantum states in CeLiBi 2

Here we report the discovery of CeLiBi 2 , the first example of a material in the tetragonal CeTX 2 ( T = transition metal; X = pnictogen) family wherein an alkali cation replaces the typical transition metal. Magnetic susceptibility and neutron powder diffraction measurements are consistent with a crystal-field Γ 6 ground-state Kramers doublet that orders antiferromagnetically below T N = 3.4 K with an incommensurate propagation wave vector k = ( 0, 0.0724(4), 0.5) that generates a nanometric modulation of the magnetic structure. The best model of the ordered state is an elliptical cycloid with Ce moments primarily residing in the ab plane. This is highly unusual, as all other Γ 6 CeTX 2 members order ferromagnetically. Further, we observe an atypical hard-axis metamagnetic transition at 2 T in magnetostriction, magnetization, and resistivity measurements. CeLiBi 2 is a rare example of a highly conductive material with dominant skew scattering leading to a large anomalous Hall effect. Quantum oscillations with five frequencies arise in magnetostriction and magnetic susceptibility data to T = 30 K and μ 0 H = 55 T, which indicate small Fermi pockets of light carriers with effective masses as low as 0.07 m e . Density functional theory calculations indicate that square-net Dirac-like Bi-p bands are responsible for these ultralight carriers. Together, our results show that CeLiBi 2 enables multiple atypical magnetic and electronic properties in a single clean material.

36 MATERIALS SCIENCE↗

How Well Does Kohn–Sham Regularizer Work for Weakly Correlated Systems?

Kohn–Sham regularizer (KSR) is a differentiable machine learning approach to finding the exchange-correlation functional in Kohn–Sham density functional theory that works for strongly correlated systems. Here we test KSR for a weak correlation. We propose spin-adapted KSR (sKSR) with trainable local, semilocal, and nonlocal approximations found by minimizing density and total energy loss. We assess the atoms-to-molecules generalizability by training on one-dimensional (1D) H, He, Li, Be, and Be 2+ and testing on 1D hydrogen chains, LiH, BeH 2 , and helium hydride complexes. The generalization error from our semilocal approximation is comparable to other differentiable approaches, but our nonlocal functional outperforms any existing machine learning functionals, predicting ground-state energies of test systems with a mean absolute error of 2.7 mH.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Benchmarking Density Functional Theory Methods for Efficient Calculations of a Strongly Correlated Li 1– x Ni 1– y O 2−δ System

Transition metal oxides (TMOs), such as LiNiO 2 , are promising candidates for energy storage and electronic devices due to their unique electronic properties, exceptional physical and chemical characteristics, and ability to adopt multiple oxidation states. However, accurately predicting their properties using mean-field density functional theory (DFT) is challenging due to the presence of strongly correlated d-electrons and the complex interplay between their structural, electronic, and magnetic responses. These challenges are further exacerbated by the need to model defects, surfaces, and interfaces, which require computationally efficient, large-scale simulations. To address these issues, we carry out a benchmark study on the Li 1–x NiO 2 system, evaluating the performance of several popular functionals. Our findings demonstrate that combining SCAN functional relaxation with single-step HSE calculations provides a practical and scalable computational strategy. This approach balances accuracy and efficiency, enabling high-throughput simulations of strongly correlated TMOs and improved predictive modeling capability of TMOs for practical applications.

25 ENERGY STORAGE↗

Analytic gradients for compressed multistate pair-density functional theory

Photochemical reactions often involve states that are closely coupled due to near degeneracies, for example by proximity to conical intersections. Therefore, a multistate method is used to accurately describe these states; for example, ordinary perturbation theory is replaced by quasidegenerate perturbation theory. Multiconfiguration pair-density functional theory (MC-PDFT) provides an efficient way to approximate the full dynamical correlation energy of strongly correlated systems, and we recently proposed compressed multistate pair-density functional theory (CMS-PDFT) to treat closely coupled states. In the present paper, we report the implementation of analytic gradients for CMS-PDFT in both OpenMolcas and PySCF, and here we illustrate the use of these gradients by applying the method to the excited states of formaldehyde and phenol.

74 ATOMIC AND MOLECULAR PHYSICS↗

Novel phenomena in transition-metal oxide thin films and heterostructures with strong correlations and spin-orbit coupling

Transition-metal oxides have been a central subject of condensed matter physics for decades. In addition to novel electronic states driven by the influence of strong correlation, relativistic spin-orbit coupling effects have recently attracted much attention for their potential to explore topological phenomena. Here, in this article, we review various experimental and theoretical studies on transition-metal oxides with focus on thin films and heterostructures where their physics is much influenced by correlation effects and spin-orbit coupling. The combination of the heterostructure geometry together with correlation and topology leads to a variety of novel states here reviewed. We also discuss perspectives for future research in this broad promising area.

36 MATERIALS SCIENCE↗

Interpretations of ground-state symmetry breaking and strong correlation in wavefunction and density functional theories

Significance The ground state of a quantum mechanical system is the lowest-energy eigenstate of the Hamiltonian. In isolation, it persists unchanged forever, with symmetries dictated by those of the Hamiltonian. But near-eigenstates of broken symmetry can persist for long times, even on the scale of human measurement. The appearance of broken symmetries of the electron density or spin density in a density functional calculation can reveal strong correlations among the electrons that are present in a symmetry-unbroken wavefunction. Symmetry breaking can arise when a wave-like fluctuation drops to zero frequency. The presented examples are the stretched hydrogen molecule, antiferromagnetism in solids, and the static charge-density wave in a low-density jellium, which is shown quantitatively to be a zero-frequency plasma wave.

42 ENGINEERING↗

Weighted active space protocol for multireference machine-learned potentials

Multireference methods such as multiconfiguration pair-density functional theory accurately capture electronic correlation in systems with strong multiconfigurational character, but their cost precludes direct use in molecular dynamics. Combining these methods with machine-learned interatomic potentials (MLPs) can extend their reach. However, the sensitivity of multireference calculations to the choice of the active space complicates the consistent evaluation of energies and gradients across structurally diverse nuclear configurations. To overcome this limitation, we introduce the weighted active space protocol (WASP), a systematic approach to assign a consistent active space for a given system across uncorrelated configurations. By integrating WASP with MLPs and enhanced sampling techniques, we propose a data-efficient active learning cycle that enables the training of an MLP on multireference data. We demonstrated the approach on the TiC + -catalyzed C–H activation of methane, a reaction that poses challenges for Kohn–Sham density functional theory due to its significant multireference character. This framework enables accurate and efficient modeling of catalytic dynamics, establishing a paradigm for simulating complex reactive processes beyond the limits of conventional electronic-structure methods.

enhanced sampling↗

Engineering 2D Square Lattice Hubbard Models in 90° Twisted GeX/SnX (X =S, Se) Moiré Superlattices

Because of the large-period superlattices emerging in moiré two-dimensional (2D) materials, electronic states in such systems exhibit low energy flat bands that can be used to simulate strongly correlated physics in a highly tunable setup. While many investigations have thus far focused on moiré flat bands and emergent correlated electron physics in triangular, honeycomb, and quasi-one-dimensional lattices, tunable moiré realizations of square lattices subject to strong correlations remain elusive. Here, in this work, we propose a feasible scheme to construct moiré square lattice systems by twisting two or more layers of 2D materials in a rectangular lattice by 90°. We demonstrate the concept with twisted GeX/SnX (X =S, Se) moiré superlattices and calculate their electronic structures from first principles. We show that the lowest conduction flat band in these systems can be described by a square lattice Hubbard model with parameters which can be controlled by varying the choice of host materials, number of layers, and external electric fields. In particular, twisted double bilayer GeSe realizes a square lattice Hubbard model with strong frustration due to the next-nearest-neighbor hopping that could host unconventional superconductivity, in close analogy to the Hubbard model for copper-oxygen planes of cuprate high-temperature superconductors. The presented scheme uses 90° twisted 2D materials with rectangular unit cells as a promising platform for realizing the physical phenomena of square lattice Hubbard models, establishing a new route for studying its rich phase diagram of magnetism, charge order, and unconventional superconductivity in a highly tunable setting.

2-dimensional systems↗

Hubbard-corrected oxide formation enthalpies without adjustable parameters

A density functional theory (DFT) approach to computing transition metal oxide heat of formation without adjustable parameters is presented. Different degrees of d-electron localization in oxides are treated within the DFT+U approach with site-dependent, first-principles Hubbard U-parameters obtained from linear response theory, and delocalized states in the metallic phases are treated without Hubbard corrections. Comparison of relative stabilities of these differently treated phases is enabled by a local d-electron density matrix-dependent model, which was found by genetic programming against experimental reference formation enthalpies. This mathematically simple model does not explicitly depend on the Hubbard-corrected ionic species and is shown to reproduce the heats of formation of the Mott insulators Ca 2 RuO 4 and Y 2 Ru 2 O 7 within ~3% of experimental results, where the experimental training data did not contain Ru oxides. This newly developed method thus absolves from the need for element-specific corrections fitted to experiments in existing Hubbard-corrected approaches to the prediction of reaction energies of transition metal oxides and metals. The absence of fitting parameters opens up here the possibility to predict relative thermodynamic stabilities and reaction energies involving d-states of varying degree of localization at transition metal oxide interfaces and defects, where site-dependent U-parameters will be particularly important and devising a fitting scheme against experimental data with predictive power would be exceedingly difficult.

transition metal oxides↗

Weyl semimetal in the rare-earth hexaboride family supporting a pseudonodal surface and a giant anomalous Hall effect

Rare-earth hexaborides offer a rich tapestry for exploring the interplay of topological orders, magnetism, and electron-correlation effects associated with protected quantum phenomena. Using first-principles modeling, we identify in cerium hexaboride (CeB 6 ) a topological state that harbors a pseudonodal surface (PNS) along with double-Weyl fermions. An analysis of the electronic states in the vicinity of the Fermi level reveals the presence of M z -mirror-symmetry-protected band crossings that involve itinerant (d) and localized (f) states of Ce. These band crossings are found to remain nearly gapless as one goes away from the mirror plane, driving the formation of PNS states. We also demonstrate the presence of coexisting double-Weyl fermions with chiral charges of ±2 on the C 4z rotational axis. By analyzing the Berry curvature field associated with the topological states, we predict that CeB 6 supports a giant anomalous Hall conductivity of -1840Ω -1 cm -1 in the slightly electron-doped regime, which would be larger than the values reported to date in magnetic Weyl semimetals. Our study thus indicates that CeB 6 , which combines ferromagnetic order with topological physics in a heavy-fermion matrix, would provide a promising playground for investigating topological magnetic states with giant anomalous Hall responses.

36 MATERIALS SCIENCE↗

Hourglass-type bulk Ni 3 d band and Ce 4 f Kondo resonance states in the potential topological Kondo semimetal CeNiSn via angle-resolved photoemission spectroscopy

The electronic structure of CeNiSn, which is a potential topological Kondo insulator and a Dirac nodal-loop semimetal, has been investigated by employing temperature (T) dependent angle-resolved photoemission spectroscopy (ARPES). The Fermi surfaces (FSs) and the band structures of CeNiSn for three orthogonal crystallographic planes are measured, in which both the very dispersive bands and the flat bands are observed, having mainly the Ni 3d character and the Ce 4f character, respectively. The measured FSs and ARPES bands agree reasonably well with the density functional theory (DFT) calculations. The Fermi-edge (E F ) photon energy (hν) map along k b (=k (010) ) shows that the metallic E F -crossing states on the (010) surface have the three-dimensional character, suggesting that the observed E F -crossing metallic states do not correspond to the topological surface states of the two-dimensional character. On the other hand, albeit weak, the features of the hourglass-type bulk band crossings are observed along SXS, with the energies and the slopes being similar to those predicted by the DFT calculations, supporting the Dirac semimetallic nature of CeNiSn. In T-dependent ARPES, the Ce 4f Kondo resonance states are clearly revealed at low T, which become much suppressed above ~80 K. This feature is consistent with the Kondo temperature of CeNiSn, estimated from its ρ(T) data. This work demonstrates the importance of the coherent Kondo states in determining the topological properties of CeNiSn.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Localized Active Space Pair-Density Functional Theory

Accurate quantum chemical methods for the prediction of spin-state energy gaps for strongly correlated systems are computationally expensive and scale poorly with the size of the system. This makes calculations for many experimentally interesting molecules impractical even with abundant computational resources. Previous work has shown that the localized active space (LAS) self-consistent field (SCF) method can be an efficient way to obtain multiconfiguration SCF wave functions of comparable quality to the corresponding complete active space (CAS) ones. To obtain quantitative results, a post-SCF method is needed to estimate the complete correlation energy. One such method is multiconfiguration pair-density functional theory (PDFT), which calculates the energy based on the density and on-top pair density obtained from a multiconfiguration wave function. In this work, we introduce localized-active-space PDFT, which uses a LAS wave function for subsequent PDFT calculations. The method is tested by computing spin-state energies and gaps in conjugated organic molecules and a bimetallic compound and comparing to the corresponding CAS-PDFT values.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Advances and Challenges of SCAN and r 2 SCAN Density Functionals in Transition‐Metal Compounds

ABSTRACT Transition‐metal compounds (TMCs) with open‐shell d ‐electrons are characterized by a complex interplay of lattice, charge, orbital, and spin degrees of freedom, giving rise to various fascinating applications. Often exhibiting exotic properties, these compounds are commonly classified as correlated systems due to strong inter‐electronic interactions called Hubbard U . This inherent complexity presents significant challenges to Kohn‐Sham density functional theory (KS‐DFT), the most widely used electronic structure method in condensed matter physics and materials science. While KS‐DFT is, in principle, exact for the ground‐state total energy, its exchange‐correlation energy must be approximated in practice. The mean‐field nature of KS implementations, combined with the limitations of current exchange‐correlation density functional approximations, has led to the perception that DFT is inadequate for correlated systems, particularly TMCs. Consequently, a common workaround involves augmenting DFT with an on‐site Hubbard‐like U correction. In recent years, the strongly constrained and appropriately normed (SCAN) density functional, along with its refined variant r 2 SCAN, has achieved remarkable progress in accurately describing the structural, energetic, electronic, magnetic, and vibrational properties of TMCs, challenging the traditional perception of DFT's limitations. This review explores the design principles of SCAN and r 2 SCAN, highlights their key advancements in studying TMCs, explains the mechanisms driving these improvements, and addresses the remaining challenges in this evolving field.

Zhang, Yubo [Minjiang Collaborative Center for The↗

Shortcomings of meta-GGA functionals when describing magnetism

Several recent studies have shown that SCAN, a functional belonging to the meta-generalized gradient approximation (MGGA) family, leads to significantly overestimated magnetic moments in itinerant ferromagnetic metals. However, this behavior is not inherent to the MGGA level of approximation since TPSS, for instance, does not lead to such severe overestimations. In order to provide a broader view of the accuracy of MGGAfunctionals for magnetism, we extend the assessment to more functionals but also to antiferromagnetic solids. The results show that to describe magnetism there is overall no real advantage in using a MGGA functional compared to GGAs. For both types of approximation, an improvement in ferromagnetic metals is necessarily accompanied by a deterioration (underestimation) in antiferromagnetic insulators, and vice versa. Furthermore, we also provide some analysis in order to understand in more detail the relation between the mathematical form of the functionals and the results.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Dipole Moment Calculations Using Multiconfiguration Pair-Density Functional Theory and Hybrid Multiconfiguration Pair-Density Functional Theory

The dipole moment is the molecular property that most directly indicates molecular polarity. The accuracy of computed dipole moments depends strongly on the quality of the calculated electron density, and the breakdown of single-reference methods for strongly correlated systems can lead to poor predictions of the dipole moments in those cases. Here, we derive the analytical expression for obtaining the electric dipole moment by multiconfiguration pair density functional theory (MCPDFT), and we assess the accuracy of MC-PDFT for predicting dipole moments at equilibrium and nonequilibrium geometries. We show that MC-PDFT dipole moment curves have reasonable behavior even for stretched geometries, and they significantly improve upon the CASSCF results by capturing more electron correlation. Herein, the analysis of a dataset consisting of 18 first-row transition metal diatomics and 6 main-group polyatomic molecules with multireference character suggests that MC-PDFT and its hybrid extension (HMC-PDFT) perform comparably to CASPT2 and MRCISD+Q methods and have a mean unsigned deviation of 0.2–0.3 D with respect to the best available dipole moment reference values. We explored the dependence of the predicted dipole moments upon the choice of the on-top density functional and active space, and we recommend the tPBE and hybrid tPBE0 on-top choices for the functionals combined with the moderate correlated participating orbital scheme for selecting the active space. With these choices, the mean unsigned deviations (in debyes) of the calculated equilibrium dipole moments from the best estimates are 0.77 for CASSCF, 0.29 for MC-PDFT, 0.24 for HMCPDFT, 0.28 for CASPT2, and 0.25 for MRCISD+Q. These results are encouraging because the computational cost of MC-PDFT or HMC-PDFT is largely reduced compared to the CASPT2 and MRCISD+Q methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Orbital-selective Peierls phase in the metallic dimerized chain MoOCl 2

Using ab initio density functional theory, here we systematically study the monolayer MoOCl 2 with a 4 d 2 electronic configuration. Our main result is that an orbital-selective Peierls phase (OSPP) develops in MoOCl 2 , resulting in the dimerization of the Mo chain along the b axis. Specifically, the Mo- d x y orbitals form robust molecular-orbital states inducing localized d x y singlet dimers, while the Mo- d x z / y z orbitals remain delocalized and itinerant. Our study shows that MoOCl 2 is globally metallic, with the Mo- d x y orbital bonding-antibonding splittings opening a gap and the Mo- d x z / y z orbitals contributing to the metallic conductivity. Overall, the results resemble the recently much discussed orbital-selective Mott phase but with the localized band induced by a Peierls distortion instead of Hubbard interactions. Finally, we also qualitatively discuss the possibility of OSPP in the 3 d 2 configuration, as in CrOCl 2 .

1-dimensional systems↗

Role of magnetic and structural symmetry breaking in forming the Mott insulating gap in Nb 3 ⁢Cl 8

The α-phase of the gapped insulator Nb 3 Cl 8 has recently emerged as the long-sought critical testing bed for examining the importance of strong interelectronic correlation vs symmetry breaking in understanding insulation of such Mott compounds. Structural symmetry breaking detected by density functional theory (DFT) energy lowering (such as dimer formation, disproportionation, or Jahn-Teller distortions) explains insulation in both d-electron Mott-like systems and in non-d-electron cases without recourse to strong correlation. Yet, in Nb 3 Cl 8 , structural symmetry breaking alone (viz. formation of Nb trimers) fails to explain insulation, leading instead to a partially occupied metallic flat band, in contrast with experimental observations. We examine the role of magnetic symmetry breaking, noting that Nb 3 Cl 8 is an observed paramagnet (not an antiferromagnet), thus potentially carrying also short-range ordered magnetic moments. Describing the latter as a polymorphous distribution of nonzero local moments with total zero net magnetization is demonstrated to lower the DFT total energy, while gapping the system without recourse to strong correlation or long-range magnetic order. This suggests that degeneracy removal by symmetry breaking in mean-field-like approaches—either structural, or magnetic, or both—can reduce or eliminate the need for strong correlation, allowing the use of DFT for such Mott systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗