Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Cluster expansion”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Phase Stability and Ordering in Rock Salt-Based Thermoelectrics: NaSbX 2 , AgSbX 2 , and Their Alloys with PbX and SnX (X = S, Se, Te)

Using first-principles density functional calculations, we study the phase stability and cation ordering behavior in ABX 2 , PbX-ABX 2 , and SnTe-ABTe 2 (A = Na, Ag; B = Sb; X = S, Se, Te) thermoelectrics that crystallize in rock salt-based lattices. We construct, separately for each ABX 2 system, cluster expansions fitted to T = o K energies of cation-ordered arrangements in order to identify the respective ground-state structures. We calculate the mixing energetics of A/B and {Pb, Sn}/AB cation-disordered solid solutions in ABX 2 and {Pb, Sn}X-ABX 2 systems, respectively, using special quasirandom structures. We find that (i) L11, a rock salt-based superlattice with A 1 B 1 stacking of cations along [111], is often the most favored cation ordering type across ABX 2 systems due to the ability of this superlattice structure to accommodate large size mismatches between A and B; (ii) A/B cation ordering is only weakly preferred over disorder in all ABX 2 systems, by <34 meV/cation, consistent with experimental observations of disordered rock salt phases at high temperatures; (iii) Na-containing ABX 2 systems strongly favor the formation of the ternary compound over the corresponding ground-state mixture of A 2 X + B 2 X 3 binary compounds, while Ag-containing ABX 2 systems do not; (iv) the experimentally reported noncubic phases of NaSbS 2 and AgSbS 2 are in fact distorted cation-ordered rock salt structures; and (v) all {Pb, Sn}X-ABX 2 systems show miscibility gaps in their phase diagrams, and this is partly due to the unfavorable electrostatics of ternary cation mixing relative to Pb or Sn + AB cation phase separation; however, estimation of the miscibility gap temperature using mean-field configurational entropy indicates that all systems except PbTe-AgSbTe 2 will readily mix at relatively low temperatures. Furthermore, our study of the structural properties, alloying behavior, and solubility limits in these promising materials is a crucial step forward in improving their thermoelectric performance.

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

Phase Diagrams of Alloys and Their Hydrides via On-Lattice Graph Neural Networks and Limited Training Data

Efficient prediction of sampling-intensive thermodynamic properties is needed to evaluate material performance and permit high-throughput materials modeling for a diverse array of technology applications. To alleviate the prohibitive computational expense of high-throughput configurational sampling with density functional theory (DFT), surrogate modeling strategies like cluster expansion are many orders of magnitude more efficient but can be difficult to construct in systems with high compositional complexity. We therefore employ minimal-complexity graph neural network models that accurately predict and can even extrapolate to out-of-train distribution formation energies of DFT-relaxed structures from an ideal (unrelaxed) crystallographic representation. This enables the large-scale sampling necessary for various thermodynamic property predictions that may otherwise be intractable and can be achieved with small training data sets. Two exemplars, optimizing the thermodynamic stability of low-density high-entropy alloys and modulating the plateau pressure of hydrogen in metal alloys, demonstrate the power of this approach, which can be extended to a variety of materials discovery and modeling problems.

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

Mysterious SiB 3 : Identifying the Relation between α- and β-SiB 3

Binary silicon boride SiB 3 has been reported to occur in two forms, as disordered and nonstoichiometric α-SiB 3–x , which relates to the α-rhombohedral phase of boron, and as strictly ordered and stoichiometric β-SiB 3 . Similar to other boron-rich icosahedral solids, these SiB 3 phases represent potentially interesting refractory materials. However, their thermal stability, formation conditions, and thermodynamic relation are poorly understood. Here, we map the formation conditions of α-SiB 3–x and β-SiB 3 and analyze their relative thermodynamic stabilities. α-SiB 3–x is metastable (with respect to β-SiB 3 and Si), and its formation is kinetically driven. Pure polycrystalline bulk samples may be obtained within hours when heating stoichiometric mixtures of elemental silicon and boron at temperatures 1200–1300 °C. At the same time, α-SiB 3–x decomposes into SiB6 and Si, and optimum time-temperature synthesis conditions represent a trade-off between rates of formation and decomposition. The formation of stable β-SiB 3 was observed after prolonged treatment (days to weeks) of elemental mixtures with ratios Si/B = 1:1–1:4 at temperatures 1175–1200 °C. The application of high pressures greatly improves the kinetics of SiB 3 formation and allows decoupling of SiB 3 formation from decomposition. Quantitative formation of β-SiB 3 was seen at 1100 °C for samples pressurized to 5.5–8 GPa. β-SiB 3 decomposes peritectoidally at temperatures between 1250 and 1300 °C. The highly ordered nature of β-SiB 3 is reflected in its Raman spectrum, which features narrow and distinct lines. In contrast, the Raman spectrum of α-SiB 3–x is characterized by broad bands, which show a clear relation to the vibrational modes of isostructural, ordered B 6 P. The detailed composition and structural properties of disordered α-SiB 3–x were ascertained by a combination of single-crystal X-ray diffraction and 29 Si magic angle spinning NMR experiments. Notably, the compositions of polycrystalline bulk samples (obtained at T ≤ 1200 °C) and single crystal samples (obtained from Si-rich molten Si–B mixtures at T > 1400 °C) are different, SiB 2.93(7) and SiB 2.64(2) , respectively. The incorporation of Si in the polar position of B 12 icosahedra results in highly strained cluster units. This disorder feature was accounted for in the refined crystal structure model by splitting the polar position into three sites. The electron-precise composition of α-SiB 3–x is SiB 2.5 and corresponds to the incorporation of, on average, two Si atoms in each B 12 icosahedron. Accordingly, α-SiB 3–x constitutes a mixture of B 10 Si 2 and B 11 Si clusters. The structural and phase stability of α-SiB 3–x were explored using a first-principles cluster expansion. The most stable composition at 0 K is SiB 2.5 , which however is unstable with respect to the decomposition β-SiB 3 + Si. Modeling of the configurational and vibrational entropies suggests that α-SiB 3–x only becomes more stable than β-SiB 3 at temperatures above its decomposition into SiB 6 and Si. Hence, we conclude that α-SiB 3–x is metastable at all temperatures. Density functional theory electronic structure calculations yield band gaps of similar size for electron-precise α-SiB 2.5 and β-SiB 3 , whereas α-SiB 3 represents a p-type conductor.

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

Correlated Anion Disorder in Heteroanionic Cubic TiOF 2

Resolving anion configurations in heteroanionic materials is crucial for understanding and controlling their properties. For anion-disordered oxyfluorides, conventional Bragg diffraction cannot fully resolve the anionic structure, necessitating alternative structure determination methods. We have investigated the anionic structure of anion-disordered cubic (ReO 3 -type) TiOF 2 using X-ray pair distribution function (PDF), 19 F MAS NMR analysis, density functional theory (DFT), cluster expansion modeling, and genetic-algorithm structure prediction. Our computational data predict short-range anion ordering in TiOF 2 , characterized by predominant cis-[O 2 F 4 ] titanium coordination, resulting in correlated anion disorder at longer ranges. To validate our predictions, we generated partially disordered supercells using genetic-algorithm structure prediction and computed simulated X-ray PDF data and 19 F MAS NMR spectra, which we compared directly to experimental data. To construct our simulated 19 F NMR spectra, we derived new transformation functions for mapping calculated magnetic shieldings to predicted magnetic chemical shifts in titanium (oxy)fluorides, obtained by fitting DFT-calculated magnetic shieldings to previously published experimental chemical shift data for TiF 4 . We find good agreement between our simulated and experimental data, which supports our computationally predicted structural model and demonstrates the effectiveness of complementary experimental and computational techniques in resolving anionic structure in anion-disordered oxyfluorides. From additional DFT calculations, we predict that increasing anion disorder makes lithium intercalation more favorable by, on average, up to 2 eV, highlighting the significant effect of variations in short-range order on the intercalation properties of anion-disordered materials.

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

Planar Defect Layers Template a High-Pressure InBi Polymorph

The short- and long-range order of III–V materials under high pressure has long been the subject of debate, with advancements in structural characterization leading to significant revisions to the accepted structural models. Despite these revisions, previous high-pressure structural assignments in the In–Bi system include the site-disordered β-Sn structure type, a structure type demonstrated to be nonexistent in analogous III–V systems. While X-ray diffraction is consistent with site disordering in InBi at high pressure, cluster expansion calculations indicate that disordering requires temperatures above 3000 K. Here, we propose InBi as a model material for studying unique high-pressure planar defects due to its highly anisotropic stress-dependent properties and structure. Specifically, we identify two sets of planar defects that mimic the diffraction pattern of a site disordered β-Sn structure type and are compatible with the calculated disorder barrier. We derive these defects by symmetry relations over crystallographic transitions. Density functional theory calculations of the proposed defects suggest that these defects are stabilized by diminishing interlayer separations with pressure. Further, we find that one of the proposed defects closely resembles a bulk high-pressure phase of InBi, InBi-ϵ, and we assert that the proposed defects order upon heating, acting as a template for InBi-ϵ growth. The proposed defects and their electronic structure provide a basis for the trend of superconducting critical temperature with increasing pressure. These methods for identifying defects are generalizable to other materials with reports of site disorder at high pressure, prompting a broader search for related high-pressure defects.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Atomic-scale probing of short-range order and its impact on electrochemical properties in cation-disordered oxide cathodes

Abstract Chemical short-range-order has been widely noticed to dictate the electrochemical properties of Li-excess cation-disordered rocksalt oxides, a class of cathode based on earth abundant elements for next-generation high-energy-density batteries. Existence of short-range-order is normally evidenced by a diffused intensity pattern in reciprocal space, however, derivation of local atomic arrangements of short-range-order in real space is hardly possible. Here, by a combination of aberration-corrected scanning transmission electron microscopy, electron diffraction, and cluster-expansion Monte Carlo simulations, we reveal the short-range-order is a convolution of three basic types: tetrahedron, octahedron, and cube. We discover that short-range-order directly correlates with Li percolation channels, which correspondingly affects Li transport behavior. We further demonstrate that short-range-order can be effectively manipulated by anion doping or post-synthesis thermal treatment, creating new avenues for tailoring the electrochemical properties. Our results provide fundamental insights for decoding the complex relationship between local chemical ordering and properties of crystalline compounds.

25 ENERGY STORAGE↗
DOE OSTI2023

Uncertainty-aware molecular dynamics from Bayesian active learning for phase transformations and thermal transport in SiC

Abstract Machine learning interatomic force fields are promising for combining high computational efficiency and accuracy in modeling quantum interactions and simulating atomistic dynamics. Active learning methods have been recently developed to train force fields efficiently and automatically. Among them, Bayesian active learning utilizes principled uncertainty quantification to make data acquisition decisions. In this work, we present a general Bayesian active learning workflow, where the force field is constructed from a sparse Gaussian process regression model based on atomic cluster expansion descriptors. To circumvent the high computational cost of the sparse Gaussian process uncertainty calculation, we formulate a high-performance approximate mapping of the uncertainty and demonstrate a speedup of several orders of magnitude. We demonstrate the autonomous active learning workflow by training a Bayesian force field model for silicon carbide (SiC) polymorphs in only a few days of computer time and show that pressure-induced phase transformations are accurately captured. The resulting model exhibits close agreement with both ab initio calculations and experimental measurements, and outperforms existing empirical models on vibrational and thermal properties. The active learning workflow readily generalizes to a wide range of material systems and accelerates their computational understanding.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Machine learning potential assisted exploration of complex defect potential energy surfaces

Abstract Atomic-scale defects generated in materials under both equilibrium and irradiation conditions can significantly impact their physical and mechanical properties. Unraveling the energetically most favorable ground-state configurations of these defects is an important step towards the fundamental understanding of their influence on the performance of materials ranging from photovoltaics to advanced nuclear fuels. Here, using fluorite-structured thorium dioxide (ThO 2 ) as an exemplar, we demonstrate how density functional theory and machine learning interatomic potential can be synergistically combined into a powerful tool that enables exhaustive exploration of the large configuration spaces of small point defect clusters. Our study leads to several unexpected discoveries, including defect polymorphism and ground-state structures that defy our physical intuitions. Possible physical origins of these unexpected findings are elucidated using a local cluster expansion model developed in this work.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Complexity of many-body interactions in transition metals via machine-learned force fields from the TM23 data set

Abstract This work examines challenges associated with the accuracy of machine-learned force fields (MLFFs) for bulk solid and liquid phases ofd-block elements. In exhaustive detail, we contrast the performance of force, energy, and stress predictions across the transition metals for two leading MLFF models: a kernel-based atomic cluster expansion method implemented using sparse Gaussian processes (FLARE), and an equivariant message-passing neural network (NequIP). Early transition metals present higher relative errors and are more difficult to learn relative to late platinum- and coinage-group elements, and this trend persists across model architectures. Trends in complexity of interatomic interactions for different metals are revealed via comparison of the performance of representations with different many-body order and angular resolution. Using arguments based on perturbation theory on the occupied and unoccupieddstates near the Fermi level, we determine that the large, sharpddensity of states both above and below the Fermi level in early transition metals leads to a more complex, harder-to-learn potential energy surface for these metals. Increasing the fictitious electronic temperature (smearing) modifies the angular sensitivity of forces and makes the early transition metal forces easier to learn. This work illustrates challenges in capturing intricate properties of metallic bonding with current leading MLFFs and provides a reference data set for transition metals, aimed at benchmarking the accuracy and improving the development of emerging machine-learned approximations.

Chemistry↗
DOE OSTI2024

Towards accurate prediction of configurational disorder properties in materials using graph neural networks

Abstract The prediction of configurational disorder properties, such as configurational entropy and order-disorder phase transition temperature, of compound materials relies on efficient and accurate evaluations of configurational energies. Previous cluster expansion methods are not applicable to configurationally-complex material systems, including those with atomic distortions and long-range orders. In this work, we propose to leverage the versatile expressive capabilities of graph neural networks (GNNs) for efficient evaluations of configurational energies and present a workflow combining attention-based GNNs and Monte Carlo simulations to calculate the disorder properties. Using the dataset of face-centered tetragonal gold copper without and with local atomic distortions as an example, we demonstrate that the proposed data-driven framework enables the prediction of phase transition temperatures close to experimental values. We also elucidate that the variance of the energy deviations among configurations controls the prediction accuracy of disorder properties and can be used as the target loss function when training and selecting the GNN models. The work serves as a fundamental step toward a data-driven paradigm for the accelerated design of configurationally-complex functional material systems.

Chemistry↗
DOE OSTI2026

Optimal invariant sets for atomistic machine learning

The representation of atomic configurations for machine learning models has led to numerous sets of descriptors. However, many descriptor sets are incomplete and/or functionally dependent. Incomplete sets cannot faithfully represent atomic environments. Yet complete constructions often suffer from a high degree of functional dependence, where some descriptors are functions of others. These redundant descriptors do not improve discrimination between atomic environments. We employ pattern recognition techniques to remove dependent descriptors to produce the smallest possible set that satisfies completeness. We apply this in two ways: First, we refine an existing description, the atomic cluster expansion. Second, we augment an incomplete construction, yielding a new message-passing neural network architecture that can recognize up to 5-body patterns. This architecture shows strong accuracy on state-of-the-art benchmarks while retaining low computational cost. Our results demonstrate the utility of this strategy to optimize descriptor sets across a range of descriptors and application datasets.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

The critical role of intrinsic defects and many-body interactions on the stability of MnBi2Te4

Intrinsic antisite defects pose a major challenge to understanding and predicting the exotic properties of the layered topological magnetic insulator MnBi2Te4 (MBT). In this work, we study the origin of the abundance of intrinsic defects in MBT, including many-body defect–defect interactions and many-body electronic correlations. Until now, ab initio methods have struggled to explain thermodynamic stability and properties influenced by defect behavior in MBT. We model native Mn–Bi antisite defects in MBT at finite temperatures using a cluster expansion that includes defect–defect interactions. To overcome the limitations of conventional density functional theory (DFT), we introduce a hybrid approach that incorporates high-accuracy quantum Monte Carlo (QMC) calculations, introducing missing correlations. This strategy allows for accurate estimation of defect energetics and finite-temperature properties. We compute the configurational free energy, defect concentration, and configurational heat capacity, revealing a second-order order–disorder phase transition near the experimental synthesis temperature. Our study provides the first theoretical insight into the thermodynamics of intrinsic defects in MBT. The negative free energy relative to pristine MBT at synthesis temperatures indicates that Mn–Bi antisite formation is thermodynamically spontaneous. We also present a broadly applicable general framework for correcting low-level theoretical theories using highly accurate many-body corrections from QMC.

Ghaffar, Abdul [ORNL] (ORCID:0000000241190168)↗
DOE OSTI2020

Prediction of structure and cation ordering in an ordered normal-inverse double spinel

Abstract Spinels represent an important class of technologically relevant materials, used in diverse applications ranging from dielectrics, sensors and energy materials. While solid solutions combining two “single spinels” have been explored in a number of past studies, no ordered “double” spinels have been reported. Based on our first principles computations, here we predict the existence of such a double spinel compound MgAlGaO 4 , formed by an equimolar mixing of MgAl 2 O 4 normal and MgGa 2 O 4 inverse spinels. After studying the details of its atomic and electronic structure, we use a cluster expansion based effective Hamiltonian approach with Monte Carlo simulations to study the thermodynamic behavior and cation distribution as a function of temperature. Our simulations provide strong evidence for short-ranged cation order in the double spinel structure, even at significantly elevated temperatures. Finally, an attempt was made to synthesize the predicted double spinel compound. Energy Dispersive X-ray Spectrometry and X-ray diffraction Rietveld refinements were performed to characterize the single-phase chemical composition and local configurational environments, which showed a favorable agreement with the theoretical predictions. These findings suggest that a much larger number of compounds can potentially be realized within this chemical space, opening new avenues for the design of spinel-structured materials with tailored functionality.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Chemical short-range order in derivative Cr–Ta–Ti–V–W high entropy alloys from the first-principles thermodynamic study

The development of high-entropy alloys (HEAs) focuses on exploring compositional regions in multi-component systems with all alloy elements in equal or near-equal atomic concentrations. Initially it was based on the main idea that high mixing configurational entropy contributions to the alloy free energy could promote the formation of a single solid solution phase. By using the ab-initio based Cluster Expansion (CE) Hamiltonian model constructed for the quinary bcc Cr–Ta–Ti–V–W system in combination with Monte Carlo (MC) simulations, we show that the phase stability and chemical short-range order (SRO) of the equiatomic quinary and five sub-quaternary systems, as well as their derivative alloys, can dramatically change the order–disorder transition temperatures (ODTT) as a function of alloy compositions. In particular, it has been found, that the equiatomic quaternary Ta–Ti–V–W and Cr–Ta–Ti–W alloys had the lowest order–disorder transition temperature (500 K) among all the analysed equiatomic compositions. In all investigated alloy systems, the strongest chemical ordering has been observed between Cr and V, which led to the conclusion that decreasing the concentration of either Cr or V might be beneficial in terms of decreasing the ODTT. It also predicts that increasing concentration of Ti significantly decreases the ODTT. Our analysis of chemical SRO as a function of alloy composition allows to understand the microstructure evolution of HEAs as a function of temperature in excellent agreement with available experimental observations. Importantly, our free energy of mixing and SRO calculations predict that the origin of precipitates formed by Cr- and V-rich in the sub-quaternary Cr–Ta–V–W system is driven by the thermodynamics. The modelling results are in an excellent agreement with experimental observation of Cr and V segregation in the W 0.38 Ta 0.36 Cr 0.15 V 0.11 alloy which in turns shows an exceptional radiation resistance.

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

Ab initio study of short-range ordering in vanadium-based disordered rocksalt structures

Disordered rocksalt Li-excess (DRX) compounds are attractive new cathode materials for Li-ion batteries as they contain resource-abundant metals and do not require the use of cobalt or nickel. Understanding the delithiation process and cation short-range ordering (SRO) in DRX compounds is essential to improve these promising cathode materials. Herein, we use first-principles calculations along with the cluster-expansion approach to model the disorder in DRX Li 2-x VO 3 , 0 < x < 1. We discuss the SRO of Li in tetrahedral and octahedral sites, and the order in which Li delithiates and V oxidizes with respect to local environments. We reveal that the number of nearest-neighbor V dictates the order of delithiation from octahedral sites and that V is oxidized in a manner that minimizes the electrostatic interactions among V. In conclusion, our results provide valuable insight for tailoring the performance of V-based DRX cathode materials in general by controlling the SRO features that reduce energy density.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Nature of molybdenum carbide surfaces for catalytic hydrogen dissociation using machine-learned potentials: an ensemble-averaged perspective

Molybdenum carbides with an electronic structure similar to noble metals have gained attention as a promising low-cost catalyst for biomass valorization and the hydrogen evolution reaction. However, our fundamental understanding of the catalyst surface and how different phases of these catalysts behave at varying reaction conditions is limited to ground state density functional theory calculations as ab initio molecular dynamics (AIMD) is computationally prohibitive at relevant length and time scales. Here, in this work, we train a multi-atomic cluster expansion (MACE) machine-learned interatomic potentials (MLIP) to study hydrogen dissociation and dynamics over Mo, δ-MoC, α-Mo 2 C, and β-Mo 2 C surfaces at varying temperatures and hydrogen partial pressures. Our simulations identify unique and different molecular and atomic hydrogen adsorption sites on different surfaces that do not depend on the temperature. At low hydrogen pressures, the surface coverage is monolayer, which transitions to two-layer adsorption at higher pressures. We find that atomic hydrogen diffusion and recombinations are preferred over molybdenum atom hollow sites, while the diffusion over carbon-terminated facets was negligible, signifying particularly strong C–H interactions. In contrast, molecular hydrogen adsorption occurs mostly atop Mo or the bridging sites. At a comparable hydrogen loading, β-Mo 2 C (001) is the most active surface for hydrogen dissociation reaction. This work provides insights into the dynamic nature of the hydrogen dissociation chemistry and the diversity of hydrogen adsorption sites on molybdenum carbides.

08 HYDROGEN↗
DOE OSTI2025

Structural origin of disorder-induced ion conduction in NaFePO 4 cathode materials

Diffusion in NaFePO 4 can be enhanced through amorphization. Based on computations using DFT and machine learning potentials, we ascribe this phenomenon to the formation of less constrained Na-ion environments upon disordering. Most modern battery technologies depend on solid-state crystalline cathode materials. However, some of these materials are constrained by the low ionic conductivity of their most stable phases. An example of this is maricite (NaFePO 4 ). Interestingly, experiments have shown that maricite can improve its rate capability through disordering (amorphization). However, experimental characterization of amorphous cathode materials remains a major challenge, hindering a clear understanding of the structural origin of the disorder-induced improvement in sodium-ion mobility. To address this, we here employ molecular dynamics simulations by first training a machine learning potential for NaFePO 4 based on the atomic cluster expansion approach and a batch active learning potential parameterization scheme. This potential is then applied to explore the structural and dynamical properties of NaFePO 4 glasses as cathode materials. Specifically, we investigate the effect of glass structure on sodium-ion diffusion, revealing the relative influences of short-range and medium-range order features. We find significant heterogeneity in sodium-ion diffusivity in the glass, with fast-conducting ions residing in less constrained atomic environments with fewer P and Fe neighbors. These more mobile ions are also surrounded by larger ring-type structures. Overall, the results and developed approach present promising avenues for developing high-performance glassy cathodes for next-generation batteries.

Christensen, Rasmus↗
DOE OSTI2020

Structure of chalcogen overlayers on Au(111): Density functional theory and lattice-gas modeling

Ordering of different chalcogens, S, Se, and Te, on Au(111) exhibit broad similarities but also some distinct features, which must reflect subtle differences in relative values of the long-range pair and many-body lateral interactions between adatoms. We develop lattice-gas (LG) models within a cluster expansion framework, which includes about 50 interaction parameters. These LG models are developed based on density functional theory (DFT) analysis of the energetics of key adlayer configurations in combination with the Monte Carlo (MC) simulation of the LG models to identify statistically relevant adlayer motifs, i.e., model development is based entirely on theoretical considerations. The MC simulation guides additional DFT analysis and iterative model refinement. Given their complexity, development of optimal models is also aided by strategies from supervised machine learning. The model for S successfully captures ordering motifs over a broader range of coverage than achieved by previous models, and models for Se and Te capture the features of ordering, which are distinct from those for S. More specifically, the modeling for all three chalcogens successfully explains the linear adatom rows (also subtle differences between them) observed at low coverages of ~0.1 monolayer. The model for S also leads to a new possible explanation for the experimentally observed phase with a (5 × 5)-type low energy electron diffraction (LEED) pattern at 0.28 ML and to predictions for LEED patterns that would be observed with Se and Te at this coverage.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗