Results for “Complex systems theory” — Page 11

DOE OSTI2025

Toward improved property prediction of 2D materials using many-body quantum Monte Carlo methods

The field of 2D materials has grown dramatically in the past two decades. 2D materials can be utilized for a variety of next-generation optoelectronic, spintronic, clean energy, and quantum computing applications. These 2D structures, which are often exfoliated from layered van der Waals materials, possess highly inhomogeneous electron densities and can possess short- and long-range electron correlations. The complexities of 2D materials make them challenging to study with standard mean-field electronic structure methods such as density functional theory (DFT), which relies on approximations for the unknown exchange-correlation functional. To overcome the limitations of DFT, highly accurate many-body electronic structure approaches such as diffusion Monte Carlo (DMC) can be utilized. In the past decade, DMC has been used to calculate accurate magnetic, electronic, excitonic, and topological properties in addition to accurately capturing interlayer interactions and cohesion and adsorption energetics of 2D materials. Here, this approach has been applied to 2D systems of wide interest, including graphene, phosphorene, MoS 2 , CrI 3 , VSe 2 , GaSe, GeSe, borophene, and several others. In this review article, we highlight some successful recent applications of DMC to 2D systems for improved property predictions beyond standard DFT.

2D materials↗

DOE OSTI2025

Two-dimensional ring polymer molecular dynamics determination of the MnO + + H 2 /D 2 reaction rates on a Δ-machine learned potential energy surface

In this work, we investigate the impact of nuclear quantum effects in the kinetics of the MnO + + H 2 reaction, a prototypical system for gas-phase H 2 activation by transition metal oxide ions. The DFT based potential energy surfaces (PESs) for the lowest-lying quintet and septet spin states reported in our previous work are improved by 2953 newly calculated CCSD(T)/AVDZ points using a delta-machine learning (Δ-ML) method. To examine nuclear quantum effects, the rate coefficients are computed using ring-polymer molecular dynamics (RPMD). Due to the mechanistic complexity of the reaction, two reaction coordinates are necessary to map out the free-energy surface and an extended RPMD rate theory is developed. The calculated RPMD rate coefficients on the new PES are in better agreement with experimental data. The calculated kinetic isotope effects (KIEs) range from 1.6 to 1.8, also in good agreement with the experimental values. Additionally, the converged RPMD rate coefficients are 2.3 to 3.1 times higher than their classical counterparts across the studied temperature range, suggesting the presence of moderate nuclear quantum effects in this reaction.

Liu, Yang [University of New Mexico, Albuquerque, ↗

DOE OSTI2024

Glassy Word Problems: Ultraslow Relaxation, Hilbert Space Jamming, and Computational Complexity

We introduce a family of local models of dynamics based on “word problems” from computer science and group theory, for which we can place rigorous lower bounds on relaxation timescales. These models can be regarded either as random circuit or local Hamiltonian dynamics and include many familiar examples of constrained dynamics as special cases. The configuration space of these models splits into dynamically disconnected sectors, and for initial states to relax, they must “work out” the other states in the sector to which they belong. When this problem has a high time complexity, relaxation is slow. In some of the cases we study, this problem also has high space complexity. When the space complexity is larger than the system size, an unconventional type of jamming transition can occur, whereby a system of a fixed size is not ergodic but can be made ergodic by appending a large reservoir of sites in a trivial product state. This finding manifests itself in a new type of Hilbert space fragmentation that we call fragile fragmentation. We present explicit examples where slow relaxation and jamming strongly modify the hydrodynamics of conserved densities. In one example, density modulations of wave vector q exhibit almost no relaxation until times O ( exp ( 1 / q ) ) , at which point they abruptly collapse. We also comment on extensions of our results to higher dimensions. Published by the American Physical Society 2024

Physics↗

DOE OSTI2024

Library of Advanced Materials for Engineering (LAMÉ) 5.20

Accurate and efficient constitutive modeling remains a cornerstone issue for solid mechanics analysis. Over the years, the LAMÉ advanced material model library has grown to address this challenge by implementing models capable of describing material systems spanning soft polymers to stiff ceramics including both isotropic and anisotropic responses. Inelastic behaviors including (visco)plasticity, damage, and fracture have all incorporated for use in various analyses. This multitude of options and flexibility, however, comes at the cost of many capabilities, features, and responses and the ensuing complexity in the resulting implementation. Therefore, to enhance confidence and enable the utilization of the LAMÉ library in application, this effort seeks to document and verify the various models in the LAMÉ library. Specifically, the broader strategy, organization, and interface of the library itself is first presented. The physical theory, numerical implementation, and user guide for a large set of models is then discussed. Importantly, a number of verification tests are performed with each model to not only have confidence in the model itself but also highlight some important response characteristics and features that may be of interest to end-users. Finally, in looking ahead to the future, approaches to add material models to this library and further expand the capabilities are presented.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Library of Advanced Materials for Engineering (LAMÉ) 5.22

Accurate and efficient constitutive modeling remains a cornerstone issue for solid mechanics analysis. Over the years, the LAMÉ advanced material model library has grown to address this challenge by implementing models capable of describing material systems spanning soft polymers to stiff ceramics including both isotropic and anisotropic responses. Inelastic behaviors including (visco)plasticity, damage, and fracture have all incorporated for use in various analyses. This multitude of options and flexibility, however, comes at the cost of many capabilities, features, and responses and the ensuing complexity in the resulting implementation. Therefore, to enhance confidence and enable the utilization of the LAMÉ library in application, this effort seeks to document and verify the various models in the LAMÉ library. Specifically, the broader strategy, organization, and interface of the library itself is first presented. The physical theory, numerical implementation, and user guide for a large set of models is then discussed. Importantly, a number of verification tests are performed with each model to not only have confidence in the model itself but also highlight some important response characteristics and features that may be of interest to end-users. Finally, in looking ahead to the future, approaches to add material models to this library and further expand the capabilities are presented.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Library of Advanced Materials for Engineering (LAMÉ) 5.28

Accurate and efficient constitutive modeling remains a cornerstone issue for solid mechanics analysis. Over the years, the LAMÉ advanced material model library has grown to address this challenge by implementing models capable of describing material systems spanning soft polymers to stiff ceramics including both isotropic and anisotropic responses. Inelastic behaviors including (visco)plasticity, damage, and fracture have all incorporated for use in various analyses. This multitude of options and flexibility, however, comes at the cost of many capabilities, features, and responses and the ensuing complexity in the resulting implementation. Therefore, to enhance confidence and enable the utilization of the LAMÉ library in application, this effort seeks to document and verify the various models in the LAMÉ library. Specifically, the broader strategy, organization, and interface of the library itself is first presented. The physical theory, numerical implementation, and user guide for a large set of models is then discussed. Importantly, a number of verification tests are performed with each model to not only have confidence in the model itself but also highlight some important response characteristics and features that may be of interest to end-users. Finally, in looking ahead to the future, approaches to add material models to this library and further expand the capabilities are presented.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Library of Advanced Materials for Engineering (LAMÉ) 5.30

Accurate and efficient constitutive modeling remains a cornerstone issue for solid mechanics analysis. Over the years, the LAMÉ advanced material model library has grown to address this challenge by implementing models capable of describing material systems spanning soft polymers to stiff ceramics including both isotropic and anisotropic responses. Inelastic behaviors including (visco)plasticity, damage, and fracture have all incorporated for use in various analyses. This multitude of options and flexibility, however, comes at the cost of many capabilities, features, and responses and the ensuing complexity in the resulting implementation. Therefore, to enhance confidence and enable the utilization of the LAMÉ library in application, this effort seeks to document and verify the various models in the LAMÉ library. Specifically, the broader strategy, organization, and interface of the library itself is first presented. The physical theory, numerical implementation, and user guide for a large set of models is then discussed. Importantly, a number of verification tests are performed with each model to not only have confidence in the model itself but also highlight some important response characteristics and features that may be of interest to end-users. Finally, in looking ahead to the future, approaches to add material models to this library and further expand the capabilities are presented.

36 MATERIALS SCIENCE↗

DOE OSTI2025

A review of thermo-hydro-mechanical modeling of coupled processes in fractured rock: From continuum to discontinuum perspective

Coupled thermo-hydro-mechanical (THM) processes in fractured rock are playing a crucial role in geoscience and geoengineering applications. Diverse and conceptually distinct approaches have emerged over the past decades in both continuum and discontinuum perspectives leading to significant progress in their comprehending and modeling. This review paper offers an integrated perspective on existing modeling methodologies providing guidance for model selection based on the initial and boundary conditions. By comparing various models, one can better assess the uncertainties in predictions, particularly those related to the conceptual models. The review explores how these methodologies have significantly enhanced the fundamental understanding of how fractures respond to fluid injection and production, and improved predictive capabilities pertaining to coupled processes within fractured systems. It emphasizes the importance of utilizing advanced computational technologies and thoroughly considering fundamental theories and principles established through past experimental evidence and practical experience. The selection and calibration of model parameters should be based on typical ranges and applied to the specific conditions of applications. The challenges arising from inherent heterogeneity and uncertainties, nonlinear THM coupled processes, scale dependence, and computational limitations in representing field scale fractures are discussed. Realizing potential advances on computational capacity calls for methodical conceptualization, mathematical modeling, selection of numerical solution strategies, implementation, and calibration to foster simulation outcomes that intricately reflect the nuanced complexities of geological phenomena. Future research efforts should focus on innovative approaches to tackle the hurdles and advance the state-of-the-art in this critical field of study.

Coupling scheme↗

DOE OSTI2025

Enhanced Electrocatalytic and Selective CO 2 -to-CO Reduction by a Rhenium(I) Complex Bearing 6,6′-Substituted 2,2′-Bipyridines

The electrochemical reduction of CO 2 (CO 2 RR) into value-added chemicals offers a promising route toward a circular carbon economy and reduced reliance on fossil fuels. A detailed understanding of the structural and electronic factors governing the performance of molecular CO 2 RR electrocatalysts is essential for the design of efficient, tunable systems. Here, in this study, we report a series of rhenium(I) complexes, fac-[Re I (6,6′-(R) 2 -bpy)(CO) 3 Cl] (bpy = 2,2′-bipyridine; R = mesityl (mes), 2,4,6-triisopropylphenyl (trip), or isophthalic acid (phth)) and evaluate their electrocatalytic activity. Among these, fac-[Re I (6,6′-(mes) 2 -bpy)(CO) 3 Cl] exhibited the highest performance, enabling selective CO 2 -to-CO conversion for 1 hour with Faradaic efficiency (FE) > 97%, representing an unprecedented activity level for a Re-bpy catalysts. Single-crystal X-ray diffraction and density functional theory (DFT) calculations indicated that favorable CO 2 binding could be promoted by the tilting of the 6,6′-(mes)2-bpy ligand (from the Re-CO coordination plane), providing mechanistic insight into the observed enhancement. The study consequently demonstrates a rational correlation between the CO 2 electrocatalytic performance of Re-bpy catalysts and their structural variations, as derived from X-ray data and corroborated by computational modeling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Tree tensor network hierarchical equations of motion based on time-dependent variational principle for efficient open quantum dynamics in structured thermal environments

In this work, we introduce an efficient method, TTN-HEOM, for exactly calculating the open quantum dynamics for driven quantum systems interacting with highly structured bosonic baths by combining the tree tensor network (TTN) decomposition scheme with the bexcitonic generalization of the numerically exact hierarchical equations of motion (HEOM). The method yields a series of quantum master equations for all core tensors in the TTN that efficiently and accurately capture the open quantum dynamics for non-Markovian environments to all orders in the system–bath interaction. These master equations are constructed based on the time-dependent Dirac–Frenkel variational principle, which isolates the optimal dynamics for the core tensors given the TTN ansatz. The dynamics converges to the HEOM when increasing the rank of the core tensors, a limit in which the TTN ansatz becomes exact. We introduce TENSO, tensor equations for non-Markovian structured open systems, as a general-purpose Python code to propagate the TTN-HEOM dynamics. We implement three general propagators for the coupled master equations: two fixed-rank methods that require a constant memory footprint during the dynamics and one adaptive-rank method with a variable memory footprint controlled by the target level of computational error. We exemplify the utility of these methods by simulating a two-level system coupled to a structured bath containing one Drude–Lorentz component and eight Brownian oscillators, which is beyond what can presently be computed using the standard HEOM. Our results show that the TTN-HEOM is capable of simulating both dephasing and relaxation dynamics of driven quantum systems interacting with structured baths, even those of chemical complexity, with an affordable computational cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Insulating ground state and 2−𝑘 magnetic structure of candidate Weyl hydrogen-atom K2⁢Mn3⁢(AsO4)3

The ideal Weyl "hydrogen-atom" semimetal exhibits only a single pair of Weyl nodes and no other trivial states at the Fermi energy. Such a material would be a panacea in the study of Weyl quasi-particles, allowing direct unambiguous observation of their topological properties. The alluaudite-like K2⁢Mn3⁢(AsO4)3 compound was recently proposed as such a material. Here, we use comprehensive experimental work and first-principle calculations to assess this prediction. We find K2⁢Mn3⁢(AsO4)3 crystallizes in the 𝐶⁢2/𝑐 symmetry with a quasi-one-dimensional Mn sublattice, growing as small needle-like crystals. Bulk property measurements reveal magnetic transitions at ≈8 and ≈4 K, which neutron scattering experiments show correspond to two distinct magnetic orders, first a partially ordered ferrimagnetic 𝐤𝟏=(0,0,0) structure at 8 K and a second transition of 𝐤𝟐=(1,0,0) at 4 K to a fully ordered state. Below the second transition, both ordering vectors are necessary to describe the complex magnetic structure with modulated spin magnitudes. Both of the best-fit magnetic structures in this work are found to break the symmetry necessary for the generation of Weyl nodes, though one of the magnetic structures allowed by 𝐤𝟏 does preserve this symmetry. However, the crystals are optically transparent and ellipsometry measurements reveal a large band gap, undermining expectations of semimetallic behavior. Density functional theory calculations predict an insulating antiferromagnetic ground state, in contrast to previous reports, and suggest potential frustration on the magnetic sublattice. Given the wide tunability of the alluaudite structure, we consider ways to push the system closer to a semimetallic state.

Taddei, Keith M [Argonne National Laboratory]↗

DOE OSTI2025

Parallel-in-time quantum simulation via Page and Wootters quantum time

In the past few decades, researchers have created a veritable zoo of quantum algorithms by drawing inspiration from classical computing, information theory, and even from physical phenomena. Here, we present quantum algorithms for parallel-in-time simulations that are inspired by the Page and Wootters formalism. In this framework, and thus in our algorithms, the classical time variable of quantum mechanics is promoted to the quantum realm by introducing a Hilbert space of “clock” qubits that are then entangled with the “system” qubits. We show that our algorithms can compute temporal properties over 𝑁 different times of many-body systems by only using log⁡(𝑁) clock qubits. As such, we achieve an exponential trade-off between time and spatial complexities. In addition, we rigorously prove that the entanglement created between the system qubits and the clock qubits has operational meaning, as it encodes valuable information about the system’s dynamics. We also provide a circuit depth estimation of all the protocols, showing a running time advantage in computation times over traditional sequential-in-time algorithms. In particular, for the case when the dynamics are determined by the Aubry-Andre model, we present a hybrid method for which our algorithms have a depth that only scales as 𝒪⁡(log⁡(𝑁)⁢𝑛). As a by-product, we can relate the previous schemes to the problem of equilibration of an isolated quantum system, thus indicating that our framework enables a new dimension for studying dynamical properties of many-body systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Theoretical Insights into Reaction-Induced Transformation and Tuning of Catalytic Behavior in Heterogenous Catalysis

Reaction-induced transformations in heterogenous catalysis represent diverse phenomena that challenge traditional views of static catalyst surfaces. From surface adsorbate dynamics, atomic rearrangements, to composition and phase transitions, these processes reveal the profound differences between idealized model systems under ultrahigh vacuum and the complex, evolving interfaces that govern real catalytic behaviors under reaction conditions. Here, this perspective reviews recent theoretical efforts to provide atomic-level mechanistic insights into significant reaction-induced transformations and their impact on catalytic activity and selectivity. It underscores the need for an integrated framework that combines predictive simulations with operando characterization to uncover active sites and mechanisms under realistic operating conditions. Achieving this requires accelerating existing simulations to fully capture diverse reaction-induced surface dynamics, enabling scalable and accurate modeling of catalysts as condition-dependent, dynamically evolving systems. Such approaches are critical to bridge the gap between theory and practice, offering a pathway to more impactful and predictive catalyst design.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2024

Lithium Storage Mechanisms and Electrochemical Behavior of a Molybdenum Disulfide Nanoparticle Anode

This study investigates the electrochemical behavior of molybdenum disulfide (MoS 2 ) as an anode in Li-ion batteries, focusing on the extra capacity phenomenon. Employing advanced characterization methods such as in situ and ex situ X-ray diffraction, Raman spectroscopy, X-ray photoelectron spectroscopy, and transmission electron microscopy, the research unravels the complex structural and chemical evolution of MoS 2 throughout its cycling. A key discovery is the identification of a unique Li intercalation mechanism in MoS 2 , leading to the formation of reversible Li x MoS 2 phases that contribute to the extra capacity of the MoS 2 electrode. Density function theory calculations suggest the potential for overlithiation in MoS 2 , predicting Li 5 MoS 2 as the most energetically favorable phase within the lithiation–delithiation process. Additionally, the formation of a Li-rich phase on the surface of Li 4 MoS 2 is considered energetically advantageous. After the first discharge, the battery system engages in two main reactions. One involves operation as a Li-sulfur battery within the carbonate electrolyte, and the other is the reversible intercalation and deintercalation of Li in Li x MoS 2 . The latter reaction contributes to the extra capacity of the battery. The incorporation of reduced graphene oxide as a conductive additive in MoS 2 electrodes notably improves their rate capability and cycling stability.

Lithium Rich Phase↗

DOE OSTI2025

QM Investigation of Rare Earth Ion Interactions with First Hydration Shell Waters and Protein-Based Coordination Models

Here, conventional methods for extracting rare earth metals (REMs) from mined mineral ores are inefficient, expensive, and environmentally damaging. Recent discovery of lanmodulin (LanM), a protein that coordinates REMs with high-affinity and selectivity over competing ions, provides inspiration for new REM refinement methods. Here, we used quantum mechanical (QM) methods to investigate trivalent lanthanide cation (Ln 3+ ) interactions with coordination systems representing bulk solvent water and protein binding sites. Energy decomposition analysis (EDA) showed differences in the energetic components of Ln 3+ interaction with representatives of solvent (water, H 2 O) and protein binding sites (acetate, CH 3 COO – ), highlighting the importance of accurate description of electrostatics and polarization in computational modeling of REM interactions with biological and bioinspired molecules. Relative binding free energies were obtained for Ln 3+ with coordination complexes originating from binding sites in PDB structures of a lanthanum binding peptide (PDB entry 7CCO) and LanM, with explicit consideration of the first hydration shell waters, according to quasi-chemical theory (QCT). Beyond the first shell, the bulk solvent environment was represented with an implicit continuum model. Ln 3+ interactions with (H 2 O) 9 and both binding site models became more favorable, moving down the periodic series. This trend was more pronounced with the protein binding site models than with water, resulting in affinity increasing with periodic number, except for the last REM, Lu 3+ , which bound less favorably than the preceding element, Yb 3+ . Using the truncated 7CCO binding site model, the magnitude and trend of the experimental Ln 3+ relative binding free energies for the whole 7CCO peptide were reproduced. Conversely, the previously reported experimental data for LanM show a preference for the earlier lanthanides; this is likely due to longer-range interactions and cooperative effects, which are not represented by the reduced models. Using the truncated 7CCO binding site model, the magnitude and trend of the experimental Ln 3+ relative binding free energies for the whole 7CCO peptide were reproduced. In contrast to the previously reported experimental data for LanM, the peptide preferentially binds the earlier lanthanides. This difference likely arises due to longer-range interactions and cooperative effects not represented by the peptide. Further investigation of Ln 3+ interactions with whole proteins using polarizable molecular mechanics models with explicit solvent is warranted to understand the influence of longer-ranged interactions, cooperativity, and bulk solvent. Nevertheless, the present work provides new insights into Ln 3+ interactions with biomolecules and presents an effective computational platform for designing specific single-site REM binding peptides more efficiently.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

High Pressure Synthesis of Rubidium Superhydrides

Through laser-heated diamond anvil cell experiments, we synthesize a series of rubidium superhydrides and explore their properties with synchrotron x-ray powder diffraction and Raman spectroscopy measurements, combined with density functional theory calculations. Upon heating rubidium monohydride embedded in H 2 at a pressure of 18 GPa, we form RbH 9 − I , which is stable upon decompression down to 8.7 GPa, the lowest stability pressure of any known superhydride. At 22 GPa, another polymorph, RbH 9 − II is synthesised at high temperature. Unique to the Rb-H system among binary metal hydrides is that further compression does not promote the formation of polyhydrides with higher hydrogen content. Instead, heating above 87 GPa yields RbH 5 , which exhibits two polymorphs ( RbH 5 − I and RbH 5 − II ). All of the crystal structures comprise a complex network of quasimolecular H 2 units and H − anions, with RbH 5 providing the first experimental evidence of linear H 3 − anions. Published by the American Physical Society 2025

Kuzovnikov, Mikhail A.↗

DOE OSTI2024

Computational Advances in Ionic Liquid Applications for Green Chemistry: A Critical Review of Lignin Processing and Machine Learning Approaches

The valorization and dissolution of lignin using ionic liquids (ILs) is critical for developing sustainable biorefineries and a circular bioeconomy. This review aims to critically assess the current state of computational and machine learning methods for understanding and optimizing IL-based lignin dissolution and valorization processes reported since 2022. The paper examines various computational approaches, from quantum chemistry to machine learning, highlighting their strengths, limitations, and recent advances in predicting and optimizing lignin-IL interactions. Key themes include the challenges in accurately modeling lignin’s complex structure, the development of efficient screening methodologies for ionic liquids to enhance lignin dissolution and valorization processes, and the integration of machine learning with quantum calculations. These computational advances will drive progress in IL-based lignin valorization by providing deeper molecular-level insights and facilitating the rapid screening of novel IL-lignin systems.

09 BIOMASS FUELS↗

DOE OSTI2025

Data for: Subsurface Interface Structure Controlling Local Electronic Properties of Epitaxial Graphene on SiC(0001)

Recently realized high-mobility semiconducting epitaxial graphene on silicon carbide, provided an important step towards integration of the graphene-based system into active components in post-silicon micro- and nano-electronics. However, the exact atomic-scale structure and the complex bonding configurations of the first epitaxial graphene carbon layer remain an open problem. Our recent report has shed new light on understanding this interface, where the external transverse electric field-dependent dynamic switching behavior of the Cbuffer-SiC bonds was observed. Here, using scanning tunneling microscopy and spectroscopy (STM and STS), we present the direct evidence of silicon (Si) vacancies at the interface and provide their distribution at the topmost reconstructed SiC(0001) layer. Experimental STM and density functional theory modeling data were used in the preparation of figures in a published article in the Journal of Physical Chemistry Letters. Files related to the figures and supplementary materials in the article are present in this dataset in .txt format.

Condensed matter imaging↗