Search NASA⌕ Search

SEARCH · Search NASA

Results for “First Principles”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

2D in-Plane Ordered MXene Nanosheets Derived from (Mo 2/3 Er 1/3 ) 2 AlC Rare-Earth i-MAX for Energy Storage Applications

MXenes have become one of the most versatile families of two-dimensional (2D) materials due to their high conductivity, hydrophilicity, and remarkable electrochemical performance. This has stimulated intense efforts to design and synthesize MXenes, including structurally unique in-plane ordered 2D MXenes called i-MXenes. Here, we have synthesized the quaternary rare earth (RE)-based i-MAX phase (Mo 2/3 Er 1/3 ) 2 AlC using an arc melting method, and the corresponding 2D i-MXene was then obtained through a LiF/HCl soft etching process. Literature studies have shown that Al and the RE element are etched out during the etching process, leading to the formation of pure vacancy-ordered Mo1.33C 2D i-MXene. However, our investigation reveals that upon exposure to a fluorine solution, the i-MAX phase forms RE fluoride impurities, which are challenging to remove through HCl−DI water washing and persist in the final product, resulting in impure Mo 1.33 C@Er i-MXene. These results were confirmed by various characterizations such as X-ray diffraction, Raman spectroscopy, X-ray photoelectron spectroscopy, and scanning transmission electron microscopy. Although the Mo 1.33 C@Er electrode showed a 24-fold increase in specific capacitance compared to its parent i-MAX phase, it still exhibited a high charge-transfer resistance arising from the insulating nature of RE fluoride byproducts, which adversely influence the overall capacitance behavior of the synthesized 2D Mo 1.33 C@Er i-MXenes. This study contributes to identifying pathways for the preparation of pure 2D i-MXenes from RE-based i-MAX phases and developing improved synthesis methods. With additional process optimization, the 2D i-MXene holds a strong potential for electrochemical energy storage applications. Additionally, the electronic structures of Mo 1.33 C were theoretically studied using first-principles density functional theory calculations, which revealed that pristine Mo 1.33 C is metallic, and this metallic nature is preserved even with −O, −F, and mixed functionalization.

chemical structure↗
DOE OSTI2026

Pseudocapacitive Titanium Oxynitride Nanowires for Ultrahigh Capacitance Supercapacitors

High-quality, multifunctional two-dimensional (2D) titanium oxynitide (TiNO) thin films and one-dimensional (1D) TiNO nanowires have been synthesized using a pulsed laser deposition, a simple, fast, and congruent evaporation method. First-principles calculations as a function of surface orientation and termination indicate that surface oxidation of TiNO nanowires can stabilize the (110) orientation observed experimentally. The specific capacitance value for the TiNO nanowire samples (2725 mF/cm 2 ) has been found to be nearly six times more than that of the TiNO thin film samples (400 mF/cm 2 ), which is attributed to the high packing density of TiNO nanowires over a given area. The nanowire samples have also been found to exhibit a significantly higher energy density (1.35 μWh/cm 2 ) than the TiNO thin-film samples (0.33 μWh/cm 2 ). Thus, the TiNO material system in thin-film and nanowire forms has been demonstrated to be a promising candidate for use as an electrode material in supercapacitors and other charge-storage applications.

capacitors↗
DOE OSTI2026

Synthetic Tuning of Exciton–Phonon Coupling in Janus WS 2(1-x) Se 2x Monolayers Revealed by Resonant Raman Excitation Spectroscopy for Optoelectronic Applications

Janus monolayers, such as WSSe, have broken out-of-plane symmetry and an intrinsic dipole moment, impacting exciton transport, lifetime, and phonon interactions while imbuing piezoelectric, photocatalytic, and Rashba spin-splitting properties to transition metal dichalcogenides (TMDs). The new properties of this atomically thin material can be used for optoelectronic device applications. As TMDs are converted into Janus monolayers, e.g., top selenization of WS2 to WSSe, the bandgap and structure smoothly evolve, impacting not only the formation of excitons but also their complex interactions with different phonon modes. Resonant Raman excitation profiles (REPs) are uniquely well-suited to reveal both excitonic transitions and exciton–phonon coupling. Here, the resonant REPs of $A^{'}_{1}$ WS 2 and A 1 WSSe modes are measured to understand the strength of their coupling with the A, B, and C excitonic bands of a WS2 monolayer throughout its stepwise transformation into Janus WSSe by pulsed laser deposition (PLD) of energetic selenium species. In situ Raman spectroscopy during deposition is used to controllably prepare stable intermediate Janus structures, WS 2(1-x) Se 2x (0 ≤ x ≤ 0.5), for ex situ measurement of their resonant REPs. As x increases, REPs reveal not only pronounced excitonic bands that gradually shift toward lower photon energies but also strong, mode-selective exciton–phonon coupling. First-principles resonant Raman simulations independently predict this spectral behavior and are shown capable of matching the spectrally broadened, experimentally observed REP profiles in this model system, indicating their strong predictive capability for future experiments. The combination of controlled synthesis, REP characterization, and predictive theory employed here demonstrates a powerful pathway to understand and ultimately tune exciton–phonon interactions for future quantum optical devices.

Janus monolayers↗
DOE OSTI2024

Unveiling the Electrocatalytic Hydrogen Evolution Reaction Pathway on RuP 2 through Ab Initio Grand Canonical Monte Carlo

In this study, the high catalytic reactivity of ruthenium phosphide (RuP 2 ) has been identified by first-principles density functional theory (DFT) calculations for the electrocatalytic hydrogen evolution reaction (HER). Complex surface reconstructions are considered by applying the ab initio grand canonical Monte Carlo (ai-GCMC) algorithm, efficiently providing a sufficient phase-space exploration of possible surfaces. Combined with surface-phase Pourbaix diagrams, we are able to identify the actual surfaces that obtained under specific experimental environments, thus leading to a more accurate understanding of the nature of the active sites and the binding strength of adsorbates. Specifically, through hundreds of surface reconstructions and hydrogenation states generated with ai-GCMC, we identify the most favorable surface phases of RuP 2 under aqueous acidic conditions. We discover that the HER activity is determined by multiple surfaces with different stoichiometries within a narrow electrode potential window. Low HER overpotential (η) has been found for each of the identified surfaces, as low as 0.04 V. High H-coverage reconstructed surfaces have been discovered under acidic conditions, and the surface Ru sites introduced by additional Ru adatoms or exposed by P-vacancies serve as the active sites for HER based on their nearly reversible H binding. Furthermore, this work provides atomistic insights into the origin of high HER activity on RuP 2 by exploring the dynamic surface phases of electrocatalysts and features a generalizable method to explore the reconstructed/hydrogenated surface space as a function of experimental conditions.

25 ENERGY STORAGE↗
DOE OSTI2025

Heterogeneous and Framework-Bound Copper Species Contribute to Catalytic Partial Methane Oxidation in Cu–Chabazite Zeolites

The relationship between continuous partial methane oxidation (PMO) rates and Cu site speciation in Cu-CHA zeolite catalysts is explored through differential rate measurements across a series of samples of varying compositions combined with density functional theory, first-principles thermodynamics, and statistical models that characterize Cu speciation. Under continuous PMO conditions (573 K, 0.07 kPa O 2 , 3 kPa H 2 O), Cu ions are shown to anchor to the CHA framework in both monomeric (Z 2 Cu and Z 2 CuH 2 O) and dimeric (O- and OH-bridged Cu) forms that are sensitive to the identity of the local framework anchoring site. Consequently, across the studied compositional range, Cu-CHA catalysts are predicted to contain a mixture of monomeric and dimeric Cu sites. Cu-normalized CH 3 OH formation rates extrapolated to zero conversion reflect contributions from multiple site types. Predicted CH 3 OH formation rates indicate that the Cu site reactivity toward CH 4 is influenced by zeolite composition and is likely limited by the reduction half-cycle.

03 NATURAL GAS↗
DOE OSTI2025

Enhanced Selectivity for C 2 H 4 Production from C 2 H 6 on Partially Chlorinated IrO 2 (110) Surfaces

Modifying metal oxide surfaces to limit their oxidizing activity can provide a means of improving catalytic selectivity toward the partial oxidation of light alkanes. Here, in this study, we investigated the oxidation of C 2 H 6 on Cl-modified IrO 2 (110) surfaces using temperature-programmed reaction spectroscopy (TPRS) and first-principles microkinetic modeling. We find that substituting Cl for O in the IrO 2 (110) surface enhances the selectivity for C 2 H 6 conversion to C 2 H 4 during TPRS by suppressing extensive oxidation to CO x products, while also either enhancing C 2 H 4 production or altering it to a lesser extent, depending on the initial C 2 H 6 coverage. The C 2 H 4 selectivity increased with increasing C 2 H 6 and Cl coverage, but reached a limiting value below 50%. The Cl coverage changed negligibly during C 2 H 6 oxidation, and the surface reactivity decreased only marginally for Cl coverages up to 0.5 ML (monolayer). TPRS simulations using a microkinetic model predict C 2 H 4 and CO x product yields as a function of the Cl coverage that agree closely with the experimental results. According to the simulations, C 2 H 6 conversion to C 2 H 4 occurs on Cl-IrO 2 (110) by the hydrogenation of C 2 H 3 * species adsorbed in blocked states, in which neighboring sites are occupied only by unreactive HO and Cl species. The microkinetic modeling shows that H-hopping away from surface HO groups provides a relatively efficient route for C 2 H 3 * to escape blocked configurations and dehydrogenate, and that this process can limit the C 2 H 4 selectivity on Cl-IrO 2 (110) under the conditions studied. Overall, our results demonstrate that Cl-substitution into IrO 2 (110) enhances the selectivity for C 2 H 4 production from C 2 H 6 and provides insights into the reaction mechanism that can guide strategies to further improve the C 2 H 4 selectivity.

IrO2↗
DOE OSTI2025

Low-Temperature Catalyst Redispersion: A Route to Enhanced Stability of Supported Metal Catalysts?

Sintering poses a significant challenge to achieving the long-term stability of supported metal catalysts under reaction conditions. Here, in this study, we report a low-temperature catalyst redispersion mechanism, in which platinum single atoms, which aggregate into nanoparticles under Reverse Water Gas Shift (RWGS) conditions at elevated temperatures, fragment into atomically dispersed species upon cooling. Using multimodal operando characterization combined with first-principles theoretical modeling, we track the structural evolution of Pt single atoms supported on ceria nanodomes, deposited either on ceria or ceria–titania mixed oxides. We find that fragmentation is more pronounced when cooling occurs under RWGS conditions compared to CO alone, owing to a synergistic interplay of the effects of H 2 , CO 2 , and CO. The support architecture has a strong influence on the extent of redispersion: while CO alone induces fragmentation on ceria, interfacial confinement and vacancy pinning at the ceria–titania interface suppress restructuring. In contrast, RWGS conditions overcome these barriers, enabling redispersion across both supports. These findings point toward a pathway for catalyst stabilization via reaction-induced redispersion under mild conditions.

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

Strain Engineering: Reduction of Microstrain at the Perovskite Surface via Alkali Metal Chloride Treatment Enhances Stability

Degradation of halide perovskites under a humid atmosphere is the major challenge preventing widespread commercial deployment of this material class. Here it is shown that strain engineering via alkali metal chloride treatment at the FAPbI 3 /SnO 2 interface effectively improves moisture-related stability. CsCl and KCl treatments reduce microstrain at the perovskite surface and slow the α- to δ-phase transformation. Alkali metal treatments with LiCl, NaCl, and RbCl led to an increase in microstrain and faster degradation. The compressive strain at the perovskite surface was the smallest for CsCl and was linked to improved stability. First-principles density functional theory calculations confirm the preferential formation of alkali defects at interstitial positions at the perovskite surface. Particularly CsCl and KCl treatments lead to a release of compressive strain at the perovskite surface and local structural distortions that may favor passivation of surface defects. In contrast, the room-temperature dynamics of Li interstitials result in an overall expansion of lattice volume, which may be linked to more facile lattice degradation.

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

Predictive Chemical Kinetic Modeling: Where We Succeed, Where We Struggle, and What Comes Next

Chemical kinetic modeling plays a foundational role in fields ranging from energy to environmental science, pharmaceuticals, and advanced materials. The past two decades have seen remarkable progress, particularly in modeling gas-phase reactions for thermochemical processes, leading to impactful industrial applications such as steam cracking and air quality management. However, new challenges are emerging. The successful development of systematic methodologies for the description of gas-phase kinetics opens the possibility to apply the same approach to the study of more challenging systems. Here, we review recent advances, including ab initio transition state theory-based master equation estimation of elementary rates, automated mechanism generation, machine-learning-assisted kinetics, and uncertainty quantification, and discuss the advances needed to apply the same methodological approach in areas such as heterogeneous catalysis, electrochemistry, liquid-phase and solid-state reactivity, and multiscale model integration. We advocate for the development of targeted tools, especially methods that go beyond empirical tuning toward first-principles-based predictions. We highlight the need for accessible software and AIaugmented workflows to democratize modeling for industry and academia alike. In this perspective, we call attention to not only what has worked but also what remains unsolved, advocating to avoid overemphasizing successes in scientific works at the expense of realism. The next decade should focus on predictive capability, physical accuracy, and community infrastructure (e.g., databases and services) to enable innovation across diverse fields. We argue that kinetic modeling, properly equipped, can accelerate discovery far beyond its traditional domains.

ab initio calculations↗
DOE OSTI2024

Impact of Mg Substitution on the Structure, Stability, and Properties of the Na 2 Fe 2 F 7 Weberite Cathode

Of the few weberite-type Na-ion cathodes explored to date, Na 2 Fe 2 F 7 exhibits the best performance, with capacities up to 184 mAh/g and energy densities up to 550 Wh/kg reported for this material. However, the development of robust structure–property relationships for this material is complicated by its tendency to form as a mixture of metastable polymorphs, and transform to a lower-energy Na y FeF 3 perovskite compound during electrochemical cycling. Our first-principles-guided exploration of Fe-based weberite solid solutions with redox-inactive Mg 2+ and Al 3+ predicts an enhanced thermodynamic stability of Na 2 Mg x Fe 2–x F 7 as the Mg content is increased, and the x = 0.125 composition is selected for further exploration. We demonstrate that the monoclinic polymorph (space group C2/c) of Na 2 Fe 2 F 7 (Mg0) and of a new Mg-substituted weberite composition, Na 2 Mg 0.125 Fe 1.875 F 7 (Mg0.125), can be isolated using an optimized synthesis protocol. The impact of Mg substitution on the stability of the weberite phase during electrochemical cycling, and on the extent and rate of Na (de)intercalation, is examined. Irrespective of the Mg content, we find that the weberite phase is retained when cycling over a narrow voltage window (2.8–4.0 V vs Na/Na+). Over a wider voltage range (1.9–4.0 V), Mg0 shows steady capacity fade due to its transformation to the NayFeF3 perovskite phase, while Mg0.125 displays more reversible cycling and a reduced phase transformation. Yet, Mg incorporation also leads to kinetically limited Na extraction and a reduced overall capacity. These findings highlight the need for the continued compositional optimization of weberite cathodes to improve their structural stability while maximizing their energy density.

25 ENERGY STORAGE↗
DOE OSTI2024

Dimensionality-Dependent Electronic Properties of the Highly Conducting n-Type Polymer, Poly(benzodifurandione)

Poly(benzodifurandione) (n-PBDF) has garnered significant interest as it displays the highest reported n-type electrical conductivity among π-conjugated polymers. Earlier theoretical studies of n-PBDF could not rationalize this high conductivity. Here, we explore the geometric and electronic properties of two-dimensional (2D) and three-dimensional (3D) n-PBDF networks using first-principles calculations and tight-binding models. In 2D networks, a metallic electronic configuration occurs when considering a coplanar geometry with BDF moieties bounded to protons on the same side; however, backbone torsions disrupt the metallic behavior. In contrast, all 3D architectures consistently lead to a metallic nature, which is not impacted by variations in proton positions and stacking patterns. Tight-binding models allowed us to evaluate the respective strengths of intra- and interchain electronic couplings in n-PBDF. Altogether, our investigations provide a comprehensive picture into the electronic properties of n-PBDF and shed light on how they are affected by system dimensionality, proton positions, and stacking patterns.

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

Interplay between Topological States and Rashba States as Manifested on Surface Steps at Room Temperature

The unique spin texture of quantum states in topological materials underpins many proposed spintronic applications. However, realizations of such great potential are stymied by perturbations, such as temperature and local fields imposed by impurities and defects, that can render a promising quantum state uncontrollable. Here, in this study, we report room-temperature scanning tunneling microscopy/spectroscopy observation of interaction between Rashba states and topological surface states, which manifests local electronic structure along step edges controllable by the layer thickness of thin films. The first-principles theoretical calculation elucidates the robust Rashba states coexisting with topological surface states along the surface steps with characteristic spin textures in momentum space. Furthermore, the Rashba edge states can be switched off by reducing the thickness of a topological insulator Bi 2 Se 3 to bolster their interaction with the hybridized topological surface states. The study unveils a manipulating mechanism of the spin textures at room temperature, reinforcing the necessity of thin film technology in controlling the quantum states.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Gold–Thiolate Nanocluster Dynamics and Intercluster Reactions Enabled by a Machine Learned Interatomic Potential

Monolayer protected metal clusters comprise a rich class of molecular systems and are promising candidate materials for a variety of applications. While a growing number of protected nanoclusters have been synthesized and characterized in crystalline forms, their dynamical behavior in solution, including prenucleation cluster formation, is not well understood due to limitations both in characterization and first-principles modeling techniques. Recent advancements in machine-learned interatomic potentials are rapidly enabling the study of complex interactions such as dynamical behavior and reactivity on the nanoscale. Here, we develop an Au-S-C-H atomic cluster expansion (ACE) interatomic potential for efficient and accurate molecular dynamics simulations of thiolate-protected gold nanoclusters (Au n (SCH 3 ) m ). Trained on more than 30,000 density functional theory calculations of gold nanoclusters, the interatomic potential exhibits ab initio level accuracy in energies and forces and replicates nanocluster dynamics including thermal vibration and chiral inversion. Long dynamics simulations (up to 0.1 μs time scale) reveal a mechanism explaining the thermal instability of neutral Au 25 (SR) 18 clusters. Specifically, we observe multiple stages of isomerization of the Au 25 (SR) 18 cluster, including a chiral isomer. Additionally, we simulate coalescence of two Au 25 (SR) 18 clusters and observe series of clusters where the formation mechanisms are critically mediated by ligand exchange in the form of [Au-S] n rings.

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

Kelvin Probe Force Microscopy Imaging of Plasticity in Hydrogenated Perovskite Nickelate Multilevel Neuromorphic Devices

Ion drift in nanoscale electronically inhomogeneous semiconductors is among the most important mechanisms being studied for designing neuromorphic computing hardware. However, nondestructive imaging of the ion drift in operando devices directly responsible for multiresistance states and synaptic memory represents a formidable challenge. Here, we present Kelvin probe force microscopy imaging of hydrogen-doped perovskite nickelate device channels subject to high-speed electric field pulses to directly visualize proton distribution by monitoring surface potential changes spatially, which is also supported with finite element-based electric field distribution studies. First-principles calculations provide mechanistic insights into the origin of surface potential changes as a function of hydrogen donor doping that serves as the contrast mechanism. We demonstrate 128 (7-bit) nonvolatile conductance levels in such devices relevant to in-memory computing applications. The synaptic plasticity measurements are implemented in spiking neural networks and show promising results for classification (SciKit Learn’s Iris and Wine data sets) and control (OpenAI’s CartPole-v1 and BipedalWalker-v3) simulation tasks.

Kelvin probe force microscopy↗
DOE OSTI2025

Connectivity-Dependent Exciton–Phonon Coupling in Cesium Bismuth Halide Quantum Dots

Metal halide octahedra form the fundamental functional building blocks of metal halide perovskites, dictating their structures, optical properties, electronic structures, and dynamics. Here, in this study, we show that the connectivity of bismuth halide octahedra in Cs 3 Bi 2 Br 9 and Cs 3 Bi 2 I 9 quantum dots (QDs) changes with different halide elements. We use first-principles calculations to reveal the key role of the connectivity of bismuth halide octahedra on the wave function symmetry, Huang-Rhys factor, and exciton-phonon interaction strength. Following QD synthesis via a ligand-mediated transport method, the effect of connectivity is verified with transient absorption spectroscopy, where we contrast Cs 3 Bi 2 Br 9 and Cs 3 Bi 2 I 9 QD exciton dynamics. In photoexcited Cs 3 Bi 2 I 9 QDs, phonons related to the vibrational motions of face-sharing [BiI 6 ] 3- bioctahedra couple strongly to the electronic state and drive rapid carrier relaxation. Equivalent signals are not observed for photoexcited Cs 3 Bi 2 Br 9 QDs, implying a lack of phonon involvement in band-edge absorption and subsequent exciton relaxation. Our findings suggest that structural engineering can effectively tune the exciton-phonon coupling and therefore influence exciton relaxation and recombination in perovskite nanomaterials.

TDDFT↗
DOE OSTI2025

Orbital Engineering Band Degeneracy in a Dual-Square Carbon-Oxide Framework

Electron band degeneracies in momentum space give rise to exotic quantum phenomena that have sparked intense interest in condensed matter physics and materials science. Nodal-lines─isolines in k-space formed by the incidental touching of two bands that share the same energy but belong to discrete eigenstates─arise in the presence of symmetries that preclude effective hybridization. Despite recent advances in the design, bottom-up assembly, and engineering of exotic electronic states in graphene nanomaterials, the extension of this approach to access synthetic two-dimensional (2D) quantum materials derived from metal- or covalent-organic frameworks (COFs) has lagged behind. Here we present a molecular orbital engineering approach for designing and fabricating an edge-centered dual square lattice within a π-conjugated 2D-tetraoxa[8]circulene (2D-TOC) COF. First-principles calculations and scanning tunnelling spectroscopy reveal the emergence of Frontier states at the center of a 3 × 3 lattice that give rise to Dirac nodal-lines in 2D-TOC. Our findings not only provide a general guide for the design of conjugated COFs with custom tailored electronic properties from molecular fragments but enable the exploration of emergent topological phenomena in synthetic 2D materials with potential application for high-speed, low-power data processing, transmission, and storage.

Dirac nodal line semimetal↗
DOE OSTI2025

Optomechanical Tuning of Second Harmonic Generation Anisotropy in Janus MoSSe/MoS 2 Heterostructures

Symmetry breaking in van der Waals materials enables the realization of quantum states and advanced device functionalities. Janus transition-metal dichalcogenides (TMDs) exhibit distinctive nonlinear optical properties due to their broken out-of-plane mirror symmetry. However, the dynamic control of second harmonic generation (SHG) anisotropy and resonance behavior via optical excitation remains elusive. Here, in this work, we investigate the SHG response of Janus MoSSe/MoS 2 heterostructures with 2H and 3R stackings. We can tune the SHG response by varying the incident photon wavelength from 800 to 1000 nm, which shows a resonance-dependent enhancement in intensity and a deviation from 6-fold symmetry, indicating wavelength-dependent anisotropy. The ratio between maximum and minimum intensity in the armchair directions, associated with the SHG anisotropy, reaches a value of 1.73 at the excitation wavelength of 1000 nm. Group theory analysis and first-principles calculations reveal that the observed anisotropy arises from optically induced strain. Our findings highlight the role of symmetry breaking and optical resonance contributing to the optomechanical tuning of SHG anisotropy, offering opportunities for developing Janus TMD-based photonic devices for frequency conversion, light generation, and optical switching.

Janus transition metal dichalcogenides↗
DOE OSTI2025

A Machine Learning Framework for Modeling Ensemble Properties of Atomically Disordered Materials

Atomic disorder can strongly influence material properties such as charge transport, optical response, and catalytic activity. However, efficiently modeling these disorder effects remains challenging for first-principles methods due to the cost of sampling large configurational spaces and computing complex physical quantities. Recent advances of machine learning techniques, particularly graph neural networks (GNNs), has enabled the efficient and accurate predictions of complex material properties, offering promising tools for studying disordered systems. In this work, we present a general machine-learning-assisted computational framework that integrates equivariant GNNs with Monte Carlo simulations to compute the thermodynamic and ensemble-averaged functional properties of disordered materials. Using the surface-termination-disordered MXene monolayer Ti 3 C 2 T 2–x as a representative system, we find that electrical conductivity exhibits an emergent peak near the order–disorder phase transition temperature due to the interplay between electron scattering and doping. In contrast, optical conductivity remains largely insensitive to local atomic disorder and reflects the global surface chemical composition. These results highlight the role of atomic disorder in affecting material properties and demonstrate the potential of our approach for statistically modeling disorder effects in a wide range of materials such as high-entropy alloys and spin liquids.

MXene↗