Low-Scaling Perturbative Spin-Orbit Coupling for Closed-Shell Time-Dependent DFT
Explore the source record for details and available documents.
Liao, Can Hua↗
SEARCH · Search NASA
Results for “DFT”
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.
Dynamical structure factors of warm dense matter from time-dependent orbital-free and mixed-stochastic-deterministic density functional theory
Abstract We present the first calculations of the inelastic part of the dynamical structure factor (DSF) for warm dense matter (WDM) using time-dependent orbital-free density functional theory (TD-OF-DFT) and mixed-stochastic-deterministic (mixed) Kohn Sham TD-DFT (KS TD-DFT). WDM is an intermediate phase of matter found in planetary cores and laser-driven experiments, where the accurate calculation of the DSF is critical for interpreting x-ray Thomson scattering measurements. Traditional TD-DFT methods, while highly accurate, are computationally expensive, motivating the exploration of TD-OF-DFT and mixed TD-KS-DFT as more efficient alternatives. We applied these methods to experimentally measured WDM systems, including solid-density aluminum and beryllium, compressed beryllium, and carbon–hydrogen mixtures. Our results show that TD-OF-DFT requires a dynamical kinetic energy potential in order to qualitatively capture the plasmon response. Additionally, it struggles with capturing bound electron contributions. In contrast, mixed TD-KS-DFT offers greater accuracy in distinguishing bound and free electron effects, aligning well with experimental data, though at a higher computational cost. This study highlights the trade-offs between computational efficiency and accuracy, demonstrating that TD-OF-DFT remains a valuable tool for rapid scans of parameter space, while mixed TD-KS-DFT should be preferred for high-fidelity simulations. Our findings provide insight into the future development of DFT methods for WDM and suggest potential improvements for TD-OF-DFT.
36 MATERIALS SCIENCE↗
Large-Scale Materials Modeling at Quantum Accuracy: Ab Initio Simulations of Quasicrystals and Interacting Extended Defects in Metallic Alloys
Ab initio electronic-structure has remained dichotomous between achievable accuracy and length-scale. Quantum many-body (QMB) methods realize quantum accuracy but fail to scale. Density functional theory (DFT) scales favorably but remains far from quantum accuracy. We present a framework that breaks this dichotomy by use of three interconnected modules: (i) invDFT: a methodological advance in inverse DFT linking QMB methods to DFT; (ii) MLXC: a machine-learned density functional trained with invDFT data, commensurate with quantum accuracy; (iii) DFT-FE-MLXC: an adaptive higher-order spectral finite-element (FE) based DFT implementation that integrates MLXC with efficient solver strategies and HPC innovations in FE-specific dense linear algebra, mixed-precision algorithms, and asynchronous compute-communication. Furthermore, we demonstrate a paradigm shift in DFT that not only provides an accuracy commensurate with QMB methods in ground-state energies, but also attains an unprecedented performance of 659.7 PFLOPS (43.1% peak FP64 performance) on 619,124 electrons using 8,000 GPU nodes of Frontier supercomputer.
density functional theory↗
Taming the fixed-node error in diffusion Monte Carlo via range separation
By combining density-functional theory (DFT) and wave function theory via the range separation (RS) of the interelectronic Coulomb operator, we obtain accurate fixed-node diffusion Monte Carlo (FN-DMC) energies with compact multi-determinant trial wave functions. In particular, we combine here short-range exchange-correlation functionals with a flavor of selected configuration interaction known as configuration interaction using a perturbative selection made iteratively (CIPSI), a scheme that we label RS-DFT-CIPSI. One of the take-home messages of the present study is that RS-DFT-CIPSI trial wave functions yield lower fixed-node energies with more compact multi-determinant expansions than CIPSI, especially for small basis sets. Indeed, as the CIPSI component of RS-DFT-CIPSI is relieved from describing the short-range part of the correlation hole around the electron-electron coalescence points, the number of determinants in the trial wave function required to reach a given accuracy is significantly reduced as compared to a conventional CIPSI calculation. Importantly, by performing various numerical experiments, we evidence that the RS-DFT scheme essentially plays the role of a simple Jastrow factor by mimicking short-range correlation effects, hence avoiding the burden of performing a stochastic optimization. Considering the 55 atomization energies of the Gaussian-1 benchmark set of molecules, we show that using a fixed value of mu = 0.5 bohr(-1) provides effective error cancellations as well as compact trial wave functions, making the present method a good candidate for the accurate description of large chemical systems.
Scemama, Anthony↗
Moving closer to experimental level materials property prediction using AI
Abstract While experiments and DFT-computations have been the primary means for understanding the chemical and physical properties of crystalline materials, experiments are expensive and DFT-computations are time-consuming and have significant discrepancies against experiments. Currently, predictive modeling based on DFT-computations have provided a rapid screening method for materials candidates for further DFT-computations and experiments; however, such models inherit the large discrepancies from the DFT-based training data. Here, we demonstrate how AI can be leveraged together with DFT to compute materials properties more accurately than DFT itself by focusing on the critical materials science task of predicting “formation energy of a material given its structure and composition”. On an experimental hold-out test set containing 137 entries, AI can predict formation energy from materials structure and composition with a mean absolute error (MAE) of 0.064 eV/atom; comparing this against DFT-computations, we find that AI can significantly outperform DFT computations for the same task (discrepancies of $$>0.076$$ > 0.076 eV/atom) for the first time.
36 MATERIALS SCIENCE↗
Discrete Fourier Transform Analysis in a Complex Vector Space
Alternative computational strategies for the Discrete Fourier Transform (DFT) have been developed using analysis of geometric manifolds. This approach provides a general framework for performing DFT calculations, and suggests a more efficient implementation of the DFT for applications using iterative transform methods, particularly phase retrieval. The DFT can thus be implemented using fewer operations when compared to the usual DFT counterpart. The software decreases the run time of the DFT in certain applications such as phase retrieval that iteratively call the DFT function. The algorithm exploits a special computational approach based on analysis of the DFT as a transformation in a complex vector space. As such, this approach has the potential to realize a DFT computation that approaches N operations versus Nlog(N) operations for the equivalent Fast Fourier Transform (FFT) calculation.
Dean, Bruce H.↗
Electrochemical CO2 Reduction over Metal-/Nitrogen-Doped Graphene Single-Atom Catalysts Modeled Using the Grand-Canonical Density Functional Theory
Renewably driven, electrochemical conversion of carbon dioxide into value-added products is expected to be a critical tool in global decarbonization. However, theoretical studies based on the computational hydrogen electrode largely ignore the nonlinear effects of the applied potential on the calculated results, leading to inaccurate predictions of catalytic behavior or mechanistic pathways. Here, we use grand canonical density functional theory (GC-DFT) to model electrochemical CO2 reduction (CO2R) over metal- and nitrogen-doped graphene catalysts (MNCs) and explicitly include the effects of the applied potential. We used GC-DFT to compute the CO2 to CO reaction intermediate energies at -0.3, -0.7, and -1.2 VSHE catalyzed by MNCs each doped with 1 of the 10 3d block metals coordinated by four pyridinic nitrogen atoms. Our results predict that Sc-, Ti-, Co-, Cu-, and Zn-N4Cs effectively catalyze CO2R at moderate to large reducing potentials (-0.7 to -1.2 VSHE). ZnN4C is a particularly promising electrocatalyst for CO2R to CO both at low and moderate applied potentials based on our thermodynamic analysis. Our findings also explain the observed pH independence of CO production over FeN4C and predict that the rate-determining step of CO2R over FeN4C is not *CO2- formation but rather *CO desorption. Additionally, the GC-DFT-computed density of states analysis illustrates how the electronic states of MNCs and adsorbates change non-uniformly with applied potential, resulting in a significantly increased *CO2- stability relative to other intermediates and demonstrating that the formation of the adsorbed *CO2- anion is critical to CO2R activation. This work demonstrates how GC-DFT paves the way for physically realistic and accurate theoretical simulations of reacting electrochemical systems.
CO2 reduction↗
Predicting magnetic properties of single-molecule magnets from self-interaction-free density-functional theory (Final Report)
In this project we investigated electronic structure of an intermediate-sized copper-based molecule and magnetic and hyperfine properties of several small non-magnetic and magnetic molecules including transition-metal elements by applying self-interaction corrections to density functional theory (DFT). In a sufficient number of cases, DFT-calculated electronic structure and magnetic properties of single-molecule magnets qualitatively differ from corresponding experimental data. This is partly due to self-interacting electrons within the DFT formalism. Recently, an efficient method to correct the self-interactions was proposed, i.e., Fermi-Lowdin prbital (FLO) based self-interaction corrected (SIC) methodology, within DFT. Henceforth, this method is referred to as FLO-SIC method which exists in FLOSIC code. We used this FLO-SIC method for our studies of electronic structure and magnetic and hyperfine properties of small magnetic molecules and non-magnetic molecules. Our study will provide insight into predictions of magnetic properties of single-molecule magnets where self-interaction corrections play a critical role. There are two components of this project. In the first work, we studied the electronic structure of a planar mononuclear Cu-based molecule in two oxidation states, using DFT with the FLO-SIC method. We chose this system because it is small enough and it includes a transition metal element. We found that the standard FLO-SIC method takes too much compute time even for the small transition-metal molecule and so we slightly modified the method in order to expedite the process. In the dianionic state, we found that the FLO-SIC spin density agrees quantitatively with accurate quantum chemistry methods, while DFT spin density without self-interaction corrections are severely deviated from the quantum chemistry methods. We also showed that the energy gap between the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) of the dianionic state is larger than that of the monoanionic state. This result is consistent with experimental data. In the second work, we investigated how the interaction between the electron spin and nuclear spin is affected by electron self-interactions within small non-magnetic and magnetic molecules. Such an interaction is called hyperfine interaction. For molecules without significant orbital angular momentum, the hyperfine interaction consists of Fermi contact and dipolar interaction terms. Since the Fermi contact term depends on electron spin density at the nuclear site, it would be highly affected by self-interaction corrections. Therefore, we calculated the hyperfine interaction for the small molecules using DFT with the slightly modified expedited FLO-SIC method which was obtained in the first work, and compared the results to experimental data and DFT calculations without self-interaction corrections. We found significant improvement of the Fermi contact term computed using the FLO-SIC method for small magnetic molecules. Overall, the first and second work provided positive outlook of application of the FLO-SIC method to magnetic molecules and systems including transition-metal elements.
36 MATERIALS SCIENCE↗
Orbital-Free Quantum Simulation Methods for Application to Warm Dense Matter (Final Technical Report)
Predictive simulations for prediction of condensed system behavior in state conditions far from ambient is increasingly crucial to DOE priorities. Warm dense matter (WDM) is the paradigm: temperature T > 1-15 eV, pressures P to 1 Mbar or greater. Experiments under such state conditions are difficult and costly. We summarize work driven by the need and opportunity to make free-energy density functional theory (DFT) as powerful a tool for ab initio simulation of matter under such extreme conditions as ground state DFT is for ordinary matter Advancing orbital-free DFT (OF-DFT) to eliminate the Kohn-Sham (KS) scaling bottleneck in such simulations is the other priority. The concurrent challenge for both goals is the intrinsic complexity of WDM. We summarize 15 years of successes and major progress on (1) free energy exchange-correlation functionals; (2) non-interacting free energy functionals (counterpart to T=0 Kohn-Sham kinetic energy density functionals); (3) rigorous results and constraints for free-energy DFT; (4) software for free energy DFT calculations in both conventional Kohn-Sham and OF-DFT form; (5) de-orbitalization of advanced orbital-dependent ground state functionals for use in OF-DFT; (6) demonstration calculations; (7) ancillary achievements (e.g. major review articles, secondary explorations motivated by primary goals).
36 MATERIALS SCIENCE↗
Polyphase-discrete Fourier transform spectrum analysis for the Search for Extraterrestrial Intelligence sky survey
The sensitivity of a matched filter-detection system to a finite-duration continuous wave (CW) tone is compared with the sensitivities of a windowed discrete Fourier transform (DFT) system and an ideal bandpass filter-bank system. These comparisons are made in the context of the NASA Search for Extraterrestrial Intelligence (SETI) microwave observing project (MOP) sky survey. A review of the theory of polyphase-DFT filter banks and its relationship to the well-known windowed-DFT process is presented. The polyphase-DFT system approximates the ideal bandpass filter bank by using as few as eight filter taps per polyphase branch. An improvement in sensitivity of approx. 3 dB over a windowed-DFT system can be obtained by using the polyphase-DFT approach. Sidelobe rejection of the polyphase-DFT system is vastly superior to the windowed-DFT system, thereby improving its performance in the presence of radio frequency interference (RFI).
Zimmerman, G. A.↗
Machine Learning Approaches for Rare-Earth Silicate Environmental Barrier Coating Thermochemical and Thermomechanical Property Predictions
Environmental barrier coatings (EBCs) are a necessary enabling technology for the transition from superalloys to silicon carbide (SiC) ceramic matrix composites (CMCs) in gas turbine engines for increased efficiency and decreased fuel costs. SiC-based CMCs are prone to oxidation-based degradation in the engine hot section, and rare-earth (RE) silicates are promising candidates for EBCs due to their close thermal expansion match to the composite substrate and oxidation resistance. However, the design of EBCs is hindered by the large chemical space of candidate materials and the difficulty in obtaining material properties for engineering optimization. This is especially difficult as research continues into mixed-cation or “high-entropy” RE silicates. First-principles computational methods such as density functional theory (DFT) are highly effective at calculating material properties to guide coating design but are limited by their computational cost. Atomistic simulations have the potential to both accelerate property calculations and expand the properties able to be calculated due to their lower computational compared to DFT. However, they require interatomic potentials (IAPs) specific to the material system of interest, and, to our knowledge, there are no suitable IAPs for RE silicates. Machine learning (ML) is a promising technique to accelerate material property predictions indirectly by generating IAPs for atomistic simulations or via direct prediction. In this work, we present two ML approaches to accelerate the calculation of RE silicate properties relevant to EBC design: 1) a ML-derived interatomic potential (IAP) for atomistic simulations of yttrium disilicate (Y2Si2O7) from DFT training data, and 2) a neural network (NN) model to directly predict thermochemical properties of RE silicates and oxides directly from easily obtainable unit cell parameters. Classical MD simulations using the IAP yield lattice properties and bond lengths in good agreement with both DFT and experimental results from x-ray diffraction. Thermodynamic properties calculated using the finite-displacement phonon method and quasi-harmonic approximation were orders of magnitude faster than DFT with good agreement to the DFT results. The IAP was also used to calculate properties such as coefficient of thermal expansion (CTE) that require large simulation supercells and are therefore difficult with DFT. The IAP correctly predicted the anisotropic nature of the CTE in three different phases of Y2Si2O7. The NN model predicts constant pressure heat capacity, Cp, orders of magnitude faster than DFT calculations, which can enable its use as a surrogate model for multiscale simulations. The two methods presented in this work demonstrate the utility of ML for accelerating the prediction of RE silicate properties, which can in turn accelerate EBC design and optimization.
machine learning↗
Density functional theory of material design: fundamentals and applications—II
Abstract This is the second and the final part of the review on density functional theory (DFT), referred to as DFT-II. In the first review, DFT-I, we have discussed wavefunction-based methods, their complexity, and basics of density functional theory. In DFT-II, we focus on fundamentals of DFT and their implications for the betterment of the theory. We start our presentation with the exact DFT results followed by the concept of exchange-correlation (xc) or Fermi-Coulomb hole and its relationship with xc energy functional. We also provide the exact conditions for the xc-hole, xc-energy and xc-potential along with their physical interpretation. Next, we describe the extension of DFT for non-integer number of electrons, the piecewise linearity of total energy and discontinuity of chemical potential at integer particle numbers, and derivative discontinuity of the xc potential, which has consequences on fundamental gap of solids. After that, we present how one obtains more accurate xc energy functionals by going beyond the LDA. We discuss the gradient expansion approximation (GEA), generalized gradient approximation (GGA), and hybrid functional approaches to designing better xc energy functionals that give accurate total energies. However, these functionals fail to predict properties like the ionization potential and the band gap. Thus, we next describe different methods of modelling these potentials and results of their application for calculation of the band gaps of different solids to highlight accuracy of different xc potentials. Finally, we conclude with a glimpse on orbital-free density functional theory and the machine learning approach.
36 MATERIALS SCIENCE↗
The Zintl Concept Applied to Intergrowth Structures: Electron‐Hole Matching, Stacking Preferences, and Chemical Pressures in Pd 5 InAs
Abstract Enumerating the potential stacking sequences of layers is a fundamental way to account for the structure diversity of solid state compounds. In many cases, these stacking variations represent polymorphs with only small energetic differences. Here, we examine a compound for which the preferred stacking pattern instead reveals key aspects about its chemical bonding: Pd 5 InAs. Its structure is based on the intergrowth of slabs of the AuCu 3 and PtHg 2 (or alternatively, fluorite) structure types. Two basic stacking arrangements are available to this compound, represented by the Pd 5 TlAs and HoCoGa 5 structure types. DFT total energy calculations reveal that the former outcompetes the latter by a staggering 0.65 eV/formula unit. Through a combination of DFT‐reversed approximation Molecular Orbital (DFT‐raMO) and DFT‐Chemical Pressure (DFT‐CP) analysis we trace this preference to two factors. First, with DFT‐raMO analysis, we derive a Zintl‐like bonding scheme of Pd 5 InAs. This scheme, along with the inspection of selected crystal orbitals, is then connected to preferred stacking through the coordination environments of the Pd atoms at the interface between the Pd−In and Pd−As layers. In the hypothetical HoCoGa 5 ‐type and observed Pd 5 InAs‐type structures, different Pd coordination environments arise at the interfaces. The hypothetical structure features square planar PdIn 2 As 2 units, in each of which the same 4d‐orbital serves in the Pd sublattice's role as both Lewis acid (for interactions with the As) and Lewis base (for interactions with the In). In the observed structure, tetrahedral PdIn 2 As 2 units occur instead, so that these contradictory roles are distributed to separate 4d‐orbitals, leading to more effective bonding. DFT‐CP analysis illustrates that this driving force for the Pd 5 TlAs‐type arrangement is supplemented by a favorable alignment of the packing tensions in the parent structures. Altogether, the resulting picture demonstrates how the reaction of simple intermetallic structures to form intergrowths can be guided by recognizable chemical interactions.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Revisiting Competing Paths in Electrochemical CO 2 Reduction on Copper via Embedded Correlated Wavefunction Theory
We re-evaluate two key steps in the mechanism of CO 2 reduction on copper at a higher level of theory capable of correcting inherent errors in density functional theory (DFT) approximations, namely, embedded correlated wavefunction (ECW) theory. Here, we consider the CO reduction step on Cu(111), which is critical to understanding reaction selectivity. Here, we optimize embedding potentials at the periodic plane-wave DFT level using density functional embedding theory (DFET). All possible adsorption sites (adsites) for each adsorbate then are screened with ECW theory at the catalytically active site to refine the local electronic structure. Unsurprisingly, DFT and ECW theory predict different adsite preferences, largely because of DFT’s inability to properly situate the CO 2π* level. Differing preferred adsites suggest that different reaction pathways could emerge from DFT versus ECW theory. Starting from these preferred ECW theory adsites, we then obtain reaction pathways at the plane-wave DFT level using the climbing-image nudged elastic band method to determine minimum energy paths. Thereafter, we perform ECW calculations at the catalytically active site to correct the energetics at each interpolated structure (image) along the reaction pathways. Via this approach, we confirm that the first step in CO reduction via hydrogen transfer on Cu(111) is to form hydroxymethylidyne (*COH) instead of formyl (*CHO). Although the prediction to preferentially form *COH is consistent with that of DFT, the two theories predict quite different structural and mechanistic behaviors, suggesting that verification is needed for other parts of the mechanism of CO 2 reduction, which is the subject of ongoing work.
30 DIRECT ENERGY CONVERSION↗
Defect modeling in semiconductors: the role of first principles simulations and machine learning
Abstract Point defects in semiconductors dictate their electronic and optical properties. Vacancies, interstitials, substitutional defects, and defect complexes can form in the semiconductor lattice and significantly impact its performance in applications such as solar absorption, light emission, electronics, and catalysis. Understanding the nature and energetics of point defects is essential for the design and optimization of next-generation semiconductor technologies. Here, we provide a comprehensive overview of the current state of research on point defects in semiconductors, focusing on the application of density functional theory (DFT) and machine learning (ML) in accelerating the prediction and understanding of defect properties. DFT has been instrumental in accurately calculating defect formation energies, charge transition levels, and other defect-related properties such as carrier recombination rates and lifetimes, and ion migration barriers. ML techniques, particularly neural networks, have emerged as powerful tools for enabling rapid prediction of defect properties at DFT-accuracy in order to overcome the expense of using large supercells and advanced functionals. We begin this article with a discussion of different types of point defects and complexes, their impact on semiconductor properties, and the experimental and DFT approaches typically used for their characterization. Through multiple case studies, we explore how DFT has been successfully applied to understand defect behavior across a variety of semiconductors, and how ML approaches integrated with DFT can efficiently predict defect properties and facilitate the discovery of new materials with tailored defect behavior. Overall, the advent of ‘DFT+ML’ promises to drive advancements in semiconductor technology, catalysis, and renewable energy applications, paving the way for the development of high-performance semiconductors which are defect-tolerant or have desirable dopability.
Rahman, Md Habibur (ORCID:000000027705984X)↗
Electronic structure of LaNiO 2 and CaCuO 2 from a self-consistent vertex-corrected GW approach
In this work, the electronic structure of one of the nickelates (LaNiO 2 ) and one of the cuprates (CaCuO 2 ) is studied with three self-consistent GW-based methods: scGW, sc(GW+vertex), and quasiparticle self-consistent GW. Low-energy features obtained in our study are in many respects similar to the features reported in previous density functional theory plus dynamical mean-field theory (DFT+DMFT) studies. Consistent with the DFT+DMFT conclusion, we find LaNiO 2 to be more correlated than CaCuO 2 . However, correlation effects included in our study change the DFT Fermi surface near the Γ point differently from that reported in DMFT studies. Features that are a few electronvolts away from the Fermi level are broader in our calculations than in the DFT+DMFT, which reflects the differences between the DFT and the GW methods. Our results are in qualitative agreement with previous G0W0 results, but the self-consistency results in quantitative differences. Generally, correlation effects are found to be sufficiently weak in both materials, which allows one to use totally ab initio diagrammatic approaches such as sc(GW+Vertex) and to avoid the methods with adjustable parameters (DFT+U or DFT+DMFT). However, the possibility of some strong correlations at low energy that cannot be captured by perturbative methods cannot be completely excluded. For instance, differences in the Fermi surface should be resolved, thus experimental studies are necessary.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
matsim-agents v1.0
matsim-agents is a multi-agent AI framework for atomistic materials simulation and discovery. It orchestrates large language models (LLMs), machine-learned interatomic potentials (MLIPs), and DFT codes into a single agentic loop running on laptops and DOE leadership-class supercomputers. MULTI-AGENT ORCHESTRATION A LangGraph state machine with three nodes: a Planner that converts a natural-language research objective into structured tasks; an Executor that dispatches atomistic tools and loops until the queue is empty; and an Analyst that summarizes results into a human-readable report. State is checkpointed after every step and human-in-the-loop gates can be inserted at any edge. HYPOTHESIS-DRIVEN DISCOVERY CHAT An interactive REPL (matsim-agents chat) that couples LLM dialogue with atomistic simulation. Chemical formulas are automatically detected in conversation turns and trigger a full crystal-phase exploration: structure generation → relaxation → stability scoring → result injection back into the conversation, creating a closed hypothesis-refinement loop. CRYSTAL PHASE ENUMERATION Given a composition, the phase explorer enumerates prototypes by stoichiometry: elemental (fcc/bcc/hcp/sc/diamond), binary 1:1 (rocksalt/CsCl/zincblende/ wurtzite/fluorite/rutile), ternary 1:1:3 (cubic perovskite), ternary 1:2:4 (perovskite + spinel), quaternary 1:1:2:6 (Fm-3m double perovskite). 2-D prototypes (graphene, h-BN, MoS2 2H/1T) and multilayer stacking are also supported via --include-2d and --num-layers. SUPERCELL GENERATION AND SITE DECORATION Auto-tiling to a minimum atom count (--min-atoms), explicit NxNxN tiling (--supercell), symmetry-distinct site decorations (--n-orderings), and isotropic lattice-scale sweeps (--lattice-scales) for volume bracketing. MLFF RELAXATION AND STABILITY SCORING HydraGNN (multi-headed GNN) drives structure relaxation via ASE with FIRE, BFGS, or BFGSLineSearch. Stability output: delta-E/atom ranking across phases and a max-residual-force dynamical-stability proxy. Other MLIPs (MACE, NequIP, Orb) can be plugged in through the same interface. DFT BACKENDS Quantum ESPRESSO pw.x and VASP 6.6 are first-class labellers. Both have validated GPU builds and SLURM/PBS launchers for three DOE platforms: Frontier (AMD MI250X, ROCm), Aurora (Intel PVC, oneAPI), Perlmutter (NVIDIA A100, CUDA). QE produces ~100 binaries (pw.x, ph.x, epw.x, ...). VASP supports scf, relax, vc-relax, and vc-relax-shape run types. ACTIVE-LEARNING LOOP matsim-agents al run CONFIG.yaml drives an iterative HydraGNN-DFT loop: MD generates candidates → ensemble/MC-dropout uncertainty selects the most informative → DFT labels them in parallel inside one allocation → dataset grows → HydraGNN retrains → repeat. DFT backend is a single YAML toggle (dft.backend: vasp | qe). LLM-generated seed structures are supported (no curated POSCAR library needed). Config uses ${VAR}, ${VAR:-default}, ${VAR:?msg} shell-style substitution for cross-user/cross-site portability. LLM BACKENDS Ollama (local, default), vLLM (HPC multi-GPU serving), OpenAI, Anthropic, HuggingFace Transformers+Accelerate. Selected at runtime via flag or env var with no code changes. HPC PORTABILITY Same Python entry points run on Frontier (ROCm 7.2), Aurora (oneAPI), and Perlmutter (CUDA 12). DFT and ML stacks are never co-loaded in the same shell; they couple through the scheduler and filesystem. Advanced multi-node launchers (serve, discovery-chat, single-relaxation, active-learning, QE warm-start) are provided for all three platforms. CODABENCH COMPETITION BUNDLE A self-contained benchmark: 159 atomistic test structures across 11 material classes, 5 tasks (formation energy, forces, ML relaxation, AI-DFT relaxation, phase stability ranking), public/private leaderboard split (30/70), and four ready-to-run baselines: MACE-MP-0, HydraGNN, UMA, AllScAIP.
Lupo Pasini, Massimiliano [Oak Ridge National Labo↗
Electronic Structures and Magnetism of Zr-, Th-, and U-based Metal-Organic Frameworks (MOFs) by Density Functional Theory
Metal-organic frameworks (MOFs) have recently gained wide interest as candidate materials for nuclear waste immobilization. While the fundamental thermodynamic properties, such as the substitution energies determine the favorability of radionuclide sequestration by utilization of a MOF matrix, the studies of MOF electronic structure reveal the role of d-, and/or f-electrons on changes in physical properties of actinide-containing materials. We use density functional theory (DFT) calculations to investigate the electronic structures of Zr-, Th-, and U-MOFs, including their electronic band structures and, where appropriate, their magnetic properties. We employ various DFT methods including DFT+U, collinear spin-polarization, spin-orbit coupling, and different flavors of exchange-correlation functionals to assess the robustness to the specific exchange-correlation functional. Unlike the Zr-, and Th-MOFs, the U-MOF is found to be sensitive to electron localization and spin; hence we explore the magnetic structure of the U-MOF in further detail.
Metal-Organic Framework, MOF, Actinides, DFT↗