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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Adaptive cluster expansion approach for predicting the structure evolution of graphene oxide

An adaptive cluster expansion (CE) method is used to explore surface adsorption and growth processes. Unlike the traditional CE method, suitable effective cluster interaction (ECI) parameters are determined, and then the selected fixed number of ECIs is continually optimized to predict the stable configurations with gradual increase of adatom coverage. Comparing with traditional CE method, the efficiency of the adaptive CE method could be greatly enhanced. As an application, the adsorption and growth of oxygen atoms on one side of pristine graphene was carefully investigated using this method in combination with first-principles calculations. The calculated results successfully uncover the structural evolution of graphene oxide for the different numbers of oxygen adatoms on graphene. The aggregation behavior of the stable configurations for different oxygen adatom coverages is revealed for increasing coverages of oxygen atoms. As a targeted method, adaptive CE can also be applied to understand the evolution of other surface adsorption and growth processes.

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

Ternary alloy material prediction using genetic algorithm and cluster expansion

This thesis summarizes our study on the crystal structures prediction of Fe-V-Si system using genetic algorithm and cluster expansion. Our goal is to explore and look for new stable compounds. We started from the current ten known experimental phases, and calculated formation energies of those compounds using density functional theory (DFT) package, namely, VASP. The convex hull was generated based on the DFT calculations of the experimental known phases. Then we did random search on some metal rich (Fe and V) compositions and found that the lowest energy structures were body centered cube (bcc) underlying lattice, under which we did our computational systematic searches using genetic algorithm and cluster expansion. Among hundreds of the searched compositions, thirteen were selected and DFT formation energies were obtained by VASP. The stability checking of those thirteen compounds was done in reference to the experimental convex hull. We found that the composition, 24-8-16, i.e., Fe 3 VSi 2 is a new stable phase and it can be very inspiring to the future experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Prediction of the Al-rich part of the Al-Cu phase diagram using cluster expansion and statistical mechanics

In this work, the thermodynamic properties of α-Al and other phases (GP zones, θ'', θ' and θ) in the Al-rich part of the Al-Cu system have been obtained by means of the cluster expansion formalism in combination with statistical mechanics. This information was used to build the Al-rich part of the Al-Cu phase-diagram taking into account vibrational entropic contributions for θ', as those of the other phases were negligible. The simulation predictions of the phase boundaries between α-Al and either θ'', θ' or θ phases as a function of temperature are in good agreement with experimental data and extend the phase boundaries to a wider temperature range. The DFT calculations reveal the presence of a number of metastable Guinier-Preston-zone type configurations that may coexist with α-Al and θ'' at low temperatures. They also demonstrate that θ' is the stable phase below 550K but it is replaced by θ above this temperature due to the vibrational entropic contribution to the Gibbs energy of θ'. This work shows how the combination of cluster expansion and statistical mechanics can be used to expand our knowledge of the phase diagram of metallic alloys and to provide Gibbs free energies of different phases that can be used as input in mesoscale simulations of precipitation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A systematic analysis of phase stability in refractory high entropy alloys utilizing linear and non-linear cluster expansion models

We report that obtaining surrogate models for various types of materials has long been a goal of computational material scientists and physicists, as they allow the calculation of thermodynamic quantities and properties of interest far more efficiently than DFT. Many surrogate models, including interatomic potentials generated by machine learning algorthims and cluster expansions, rely on incorporating many-body interactions such as three-body or four-body interactions. When utilizing machine learning algorithms for interatomic potentials, it is common to start with two-body interactions and systematically include three- and four-body clusters to improve the fit, stopping when including the additional clusters does not measurably change the predictive power of the model. In cluster expansions for FCC and related lattices, it is often enough to include mostly two-body interactions and a handful of three-body interactions. However, it remains to be seen whether those will be enough in the case of lattices with a lower packing fraction, such as BCC.

36 MATERIALS SCIENCE↗
DOE OSTI2024

First-principles and cluster expansion study of the effect of magnetism on short-range order in Fe–Ni–Cr austenitic stainless steels

Short-range order (SRO), the regular and predictable arrangement of atoms over short distances, alters the mechanical properties of technologically relevant structural materials such as medium/high entropy alloys and austenitic stainless steels. In this study, we present a generalized spin cluster expansion (CE) model and show that magnetism is a primary factor influencing the level of SRO present in austenitic Fe-Ni-Cr alloys. The spin CE consists of a chemical cluster expansion combined with an Ising model for Fe-Ni-Cr austenitic alloys. It explicitly accounts for local magnetic exchange interactions, thereby capturing the effects of finite temperature magnetism on SRO. Model parameters are obtained by fitting to a first-principles data set comprising both chemically and magnetically diverse FCC configurations. The magnitude of the magnetic exchange interactions are found to be comparable to the chemical interactions. Compared to a conventional implicit magnetism CE built from only magnetic ground state configurations, the spin CE shows improved performance on several experimental benchmarks over a broad spectrum of compositions, particularly at higher temperatures due to the explicit treatment of magnetic disorder. We find that SRO is strongly influenced by alloy Cr content, since Cr atoms prefer to align antiferromagnetically with nearest neighbors but become magnetically frustrated with increasing Cr concentration. Using the spin CE, we predict that increasing the Cr concentration in typical austenitic stainless steels promotes the formation of SRO and increases order-disorder transition temperatures. Furthermore, this study underscores the significance of considering magnetic interactions explicitly when exploring the thermodynamic properties of complex transition metal alloys. It also highlights guidelines for customizing SRO through adjustments of alloy composition.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Atomic cluster expansion potential for large scale simulations of hydrocarbons under shock compression

We present an Atomic Cluster Expansion (ACE) machine learned potential developed for high-fidelity atomistic simulations of hydrocarbons, targeting pressures and temperatures near and above supercritical fluid regimes for molecular fluids. A diverse set of stoichiometries were covered in training, including 1:0 (pure carbon), 1:4 (methane), and 1:1 (benzene), and rich bonding environments sampled at supercritical temperatures, hydrogen rich, reactive mixtures where metastable stoichiometries arise, including 1:2 (ethylene) and 1:3 (ethane). A high-fidelity training database was constructed by performing large-scale quantum molecular dynamic simulations [density functional theory (DFT) MD] of diamond, graphite, methane, and benzene. A novel approach to selecting structures from DFT MD is also presented, which allows for the rapid selection of unique DFT MD frames from complex trajectories. Comparisons to DFT and experimental data demonstrate that the presented ACE potential accurately reproduces isotherms, carbon melting curves, radial distribution functions, and shock Hugoniots for carbon and hydrocarbon systems for pressures up to 100 GPa and temperatures up to 6000 K for hydrocarbon systems and up to 9000 K for pure carbon systems. This work delivers a potential that can be used for accurate, large-scale simulations of shocked hydrocarbons and demonstrates a methodology for fitting and validating machine learning interatomic potentials to complex molecular environments, which can be applied to energetic materials in future works.

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

Cluster expansion by transfer learning for phase stability predictions

Recent progress towards universal machine-learned interatomic potentials holds considerable promise for materials discovery. Yet the accuracy of these potentials for predicting phase stability may still be limited. In contrast, cluster expansions provide accurate phase stability predictions but are computationally demanding to parameterize from first principles, especially for structures of low dimension or with a large number of components, such as interfaces or multimetal catalysts. We overcome this trade-off via transfer learning. Using Bayesian inference, we incorporate prior statistical knowledge from machine-learned and physics-based potentials, enabling us to sample the most informative configurations and to efficiently fit first-principles cluster expansions. Furthermore, this algorithm is tested on Pt:Ni, showing robust convergence of the mixing energies as a function of sample size with reduced statistical fluctuations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Shadow Molecular Dynamics and Atomic Cluster Expansions for Flexible Charge Models

Here, a shadow molecular dynamics scheme for flexible charge models is presented where the shadow Born–Oppenheimer potential is derived from a coarse-grained approximation of range-separated density functional theory. The interatomic potential, including the atomic electronegativities and the charge-independent short-range part of the potential and force terms, is modeled by the linear atomic cluster expansion (ACE), which provides a computationally efficient alternative to many machine learning methods. The shadow molecular dynamics scheme is based on extended Lagrangian (XL) Born–Oppenheimer molecular dynamics (BOMD). XL-BOMD provides stable dynamics while avoiding the costly computational overhead associated with solving an all-to-all system of equations, which normally is required to determine the relaxed electronic ground state prior to each force evaluation. To demonstrate the proposed shadow molecular dynamics scheme for flexible charge models using atomic cluster expansion, we emulate the dynamics generated from self-consistent charge density functional tight-binding (SCC-DFTB) theory using a second-order charge equilibration (QEq) model. The charge-independent potentials and electronegativities of the QEq model are trained for a supercell of uranium oxide (UO 2 ) and a molecular system of liquid water. The combined ACE+XL-QEq molecular dynamics simulations are stable over a wide range of temperatures both for the oxide and for the molecular systems and provide a precise sampling of the Born–Oppenheimer potential energy surfaces. Accurate ground Coulomb energies are produced by the ACE-based electronegativity model during an NVE simulation of UO 2 , predicted to be within 1 meV of those from SCC-DFTB on average during comparable simulations.

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

Shadow molecular dynamics and atomic cluster expansions for flexible charge models

A shadow molecular dynamics scheme for flexible charge models is presented, where the shadow Born-Oppenheimer potential is derived from a coarse-grained approximation of range-separated density functional theory. The interatomic potential, including the atomic electronegativities and the charge-independent short-range part of the potential and force terms, are modeled by the linear atomic cluster expansion (ACE), which provides a computationally efficient alternative to many machine learning methods. The shadow molecular dynamics scheme is based on extended Lagrangian (XL) Born-Oppenheimer molecular dynamics (BOMD) [Eur. Phys. J. B 94, 164 (2021)]. XL-BOMD provides a stable dynamics, while avoiding the costly computational overhead associated with solving an all-to-all system of equations, which normally is required to determine the relaxed electronic ground state prior to each force evaluation. To demonstrate the proposed shadow molecular dynamics scheme for flexible charge models using the atomic cluster expansion, we emulate the dynamics generated from self-consistent charge density functional tight-binding (SCC-DFTB) theory using a second-order charge equilibration (QEq) model. The charge-independent potentials and electronegativities of the QEq model are trained for a supercell of uranium oxide (UO2) and a molecular system of liquid water. The combined ACE + XL-QEq dynamics are stable over a wide range of temperatures both for the oxide and the molecular systems, and provide a precise sampling of the Born-Oppenheimer potential energy surfaces. Accurate ground Coulomb energies are produced by the ACE-based electronegativity model during an NVE simulation of UO 2 , predicted to be within 1 meV of those from SCC-DFTB on average during comparable simulations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2017

Rapidly convergent cluster expansion and application to lithium ion battery materials

The convergence of a cluster expansion for lithium transition-metal (TM) oxides is improved by explicit treatment of TM magnetic moments. The approach is applied to layered LiC o y N i 1 - y O 2 (NC). The ground state and low-lying excited state structures are identified, and the distribution of TM ions and magnetic moment in those structures is investigated to explain the origin of Ni-antisite ions and Jahn-Teller distortion. In this work, the developed model also reveals the mechanisms governing the atomic arrangement of NC, including in-plane Co-Co vs Co-Ni competition, magnetic frustration vs disproportionation competition, and cationic interactions spanning adjacent layers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Predicting the Pseudocapacitive Windows for MXene Electrodes with Voltage-Dependent Cluster Expansion Models

MXene transition-metal carbides and nitrides are of growing interest for energy storage applications. These compounds are especially promising for use as pseudocapacitive electrodes due to their ability to convert energy electrochemically at fast rates. Using voltage-dependent cluster expansion models, we predict the charge storage performance of MXene pseudocapacitors for a range of electrode compositions. M 3 C 2 O 2 electrodes based on group-VI transition metals have up to 80% larger areal energy densities than prototypical titanium-based (e.g., Ti 3 C 2 O 2 ) MXene electrodes. We attribute this high pseudocapacitance to the Faradaic voltage windows of group-VI MXene electrodes, which are predicted to be 1.2 to 1.8 times larger than those of titanium-based MXenes. The size of the pseudocapacitive voltage window increases with the range of oxidation states that are accessible to the MXene transition metals. Here, by similar mechanisms, the presence of multiple ions in the solvent (Li + and H + ) leads to sharp changes in the transition-metal oxidation states and can significantly increase the charge capacity of MXene pseudocapacitors.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Phase diagram of magnetic shape memory alloy Ni 50 Mn $50–x$ In $x$ , 0 < $x$ , 25 from first principles, via spin cluster expansion and phonon vibrational entropies

The metamagnetic shape memory Heusler alloy Ni 50 Mn $50–x$ In $x$ exhibits a rich phase diagram featuring competing magnetic states, coupled magnetic–structural phase transitions, and strong compositional sensitivity. Existing first-principles approaches struggletocapturetheintertwinedchemical, magnetic, andvibrationaleffectsinthesealloys, necessitating a more integrated modeling framework. We develop a spin cluster expansion (spin-CE) framework augmented by a quasi-harmonic phonon model to capture both configurational (chemical and magnetic) and vibrational contributions to the free energy of Ni 50 Mn $50–x$ In $x$ over the full compositional range 0 ≤ x ≤25. The spin-CE includes both chemical clusters and composition-dependent Ising spin interactions, with parameters fit to a first-principles density functional theory (DFT) dataset. Using this approach, we predict the complete magnetostructural phase diagram and transformation temperatures of Ni 50 Mn $50–x$ In $x$ across the composition space. We find that vibrational entropy alone is insufficient to reproduce the martensitic transformation in the magnetic shape memory alloy regime, highlighting the essential role of magnetism. Incorporating both magnetic and vibrational contributions allows us to reproduce all experimentally known phases, including the disappearance of the stable martensite phase at a critical In concentration and the Curie temperature of the austenite phase. The method also captures the transition with increasing In in martensite from antiferromagnetic to ferromagnetic order and predicts re-entrant ferromagnetism, though the latter occurs at higher In content than reported experimentally. We discuss possible sources of this discrepancy and highlight the broader applicability of the method to other magnetostructurally complex materials, where it may offer mechanistic insight and predictive design capabilities.

Cluster expansion↗
DOE OSTI2018

Error controlling of the combined Cluster-Expansion and Wang–Landau Monte-Carlo method and its application to FeCo

Simulating order–disorder phase transitions in magnetic materials requires the accurate treatment of both the atomic and magnetic interactions, which span a vast configuration space. Taking FeCo as a prototype system, we demonstrate that this can be accurately and efficiently addressed by the combined Cluster-Expansion method and Wang–Landau algorithm. All involved key parameters in this method are carefully tested. Based on these tests, error controlling approaches are also developed. Furthermore our method is readily applicable to the study of similar second-order phase transitions in other binary and multi-component systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Multisublattice cluster expansion study of short-range ordering in iron-substituted strontium titanate

Owing to the challenges in obtaining realistic atomic configurations in large chemical phase spaces, it is not straightforward to describe structure–property relations in materials exhibiting configurational disorder. One example is iron-substituted strontium titanate (SrTi 1–x Fe x O 3–d , STF), a promising perovskite-derivative cathode material in solid oxide fuel cells that exhibits full solid solubility 0 ≤ x ≤ 1 and a tendency to exhibit short-range order. Here we demonstrate a multisublattice cluster expansion (CE) framework and apply it to STF across the full composition range. The CE approach is distinct from more traditional CE formulations in that clusters are defined explicitly by the chemical species distributed among multiple sublattices, rather than via cluster functions of occupation variables with decoration. The modified CE approach makes it easy to distinguish meaningful chemical interactions that are harder to extract from conventional CE, since for the latter chemical identity in a cluster is expressed as a product of site occupations. The least absolute shrinkage and selection operator (LASSO) is implemented as a regression analysis tool to select key clusters and avoid overfitting. We demonstrate this formulation on STF, and show that it can accurately predict configurational energies in comparison to conventional CE. From the key clusters, we identify that short-range ordering between substitutional Fe and oxygen vacancies (V O ) results in the formation of Fe–VO strings. In addition, we consider the stability of STF through CE-based Monte Carlo (MC) simulations and confirm the presence of superstructures that were previously observed in transmission electron microscopy. In this work, analysis of atomic configurations from MC samples reveals variations in the oxidation state of Fe atoms, which can be explained by the ordering tendency of Fe and V O . The cluster description and selection formalism described here may be applied to other disordered multisublattice systems for accurate and efficient material modeling.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The Generalized Green’s function Cluster Expansion: A Python package for simulating polarons

We present an efficient implementation of the Generalized Green’s function Cluster Expansion (GGCE), which is a new method for computing the ground-state properties and dynamics of polarons (single electrons coupled to lattice vibrations) in model electron-phonon systems. The GGCE works at arbitrary temperature and is well suited for a variety of electron-phonon couplings, including, but not limited to, site and bond Holstein and Peierls (Su-SchriefferHeeger) couplings, and couplings to multiple phonon modes with different energy scales and coupling strengths. Quick calculations can be performed efficiently on a laptop using solvers from NumPy and SciPy, or in parallel at scale using the PETSc sparse linear solver engine.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Assessing Mg–Sc–(rare earth) ternary phase stability via constituent binary cluster expansions

The disordered Mg–Sc body-centered cubic (bcc) phase is both lightweight and strong; however, the system is impractical for general industrial use due to the high cost of scandium. Here we propose a computationally efficient metric that assesses ternary rare earth element additions that may stabilize the bcc phase at lower Sc concentrations. We find that the bcc phase is stabilized by the ternary addition of Y or Er, but not by La, Ce, or Nd, and we validate these predictions by experimental production and characterization of Mg–Sc–(Y,Er,Nd) alloys. The results suggest a computationally efficient method to anticipate integration of ternary elements into binary systems using cluster expansions of constituent binaries.

36 MATERIALS SCIENCE↗
DOE OSTI2025

An atomic cluster expansion potential for twisted multilayer graphene

Twisted multilayer graphene, characterized by its moiré patterns arising from inter-layer rotational misalignment, serves as a rich platform for exploring quantum phenomena. Machine learning interatomic potentials (MLIPs) are a promising approach to model such systems. Our work develops a method to generate training and test datasets for fitting MLIPs that capture all possible misalignments but remain small-scale to facilitate efficient data generation and parameter estimation. To achieve this, we generate configurations with periodic boundary conditions suitable for density functional theory calculations, and then introduce an internal twist and shift within those supercell structures. Using this technique, supplemented with an active learning workflow, we fit an Atomic Cluster Expansion potential for simulating twisted multilayer graphene and test it for accuracy and robustness on a range of simulation tasks.

2D materials↗
DOE OSTI2019

Tensor optimized Fermi sphere method for nuclear matter—Power series correlated wave function and a cluster expansion

A new formalism, called “tensor optimized Fermi sphere (TOFS) method”, is developed to treat the nuclear matter using a bare interaction among nucleons. In this method, the correlated nuclear matter wave function is taken to be a power series type, Ψ{sub N}=[∑{sub n=0}{sup N}(1∕n!)F{sup n}]Φ{sub 0} and an exponential type, Ψ{sub ex}=exp(F)Φ{sub 0}, with the uncorrelated Fermi-gas wave function Φ{sub 0}, where the correlation operator F can induce central, spin–isospin, tensor, etc. correlations, and Ψ{sub ex} corresponds to a limiting case of Ψ{sub N} (N→∞). In the TOFS formalism based on Hermitian form, it is shown that the energy per particle in nuclear matter with Ψ{sub ex} can be expressed in terms of a linked-cluster expansion. On the basis of these results, we present the formula of the energy per particle in nuclear matter with Ψ{sub N}. We call the Nth-order TOFS calculation for evaluating the energy with Ψ{sub N}, where the correlation functions are optimally determined in the variation of the energy. The TOFS theory is applied for the study of symmetric nuclear matter using a central NN potential with short-range repulsion. The calculated results are fairly consistent to those of other theories such as the Brueckner–Hartree–Fock approach etc.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗