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 91 records · Page 5

Sparse expansions of multicomponent oxide configuration energy using coherency and redundancy

We report that compressed sensing has become a widely accepted paradigm to construct high dimensional cluster expansion models used for statistical mechanical studies of atomic configuration in complex multicomponent crystalline materials. However, strict sampling requirements necessary to obtain minimal coherence measurements for compressed sensing to guarantee accurate estimation of model parameters are difficult and in some cases impossible to satisfy due to the inability of physical systems to access certain configurations. Nevertheless, the dependence of energy on atomic configuration can still be adequately learned without these strict requirements by using compressed sensing by way of coherent measurements using redundant function sets known as frames. We develop a particular frame constructed from the union of all occupancy-based cluster expansion basis sets. We illustrate how using this highly redundant frame yields sparse expansions of the configuration energy of complex oxide materials that are competitive and often surpass the prediction accuracy and sparsity of models obtained from standard cluster expansions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Construction of ground-state preserving sparse lattice models for predictive materials simulations

First-principles based cluster expansion models are the dominant approach in ab initio thermodynamics of crystalline mixtures enabling the prediction of phase diagrams and novel ground states. However, despite recent advances, the construction of accurate models still requires a careful and time-consuming manual parameter tuning process for ground-state preservation, since this property is not guaranteed by default. In this paper, we present a systematic and mathematically sound method to obtain cluster expansion models that are guaranteed to preserve the ground states of their reference data. The method builds on the recently introduced compressive sensing paradigm for cluster expansion and employs quadratic programming to impose constraints on the model parameters. The robustness of our methodology is illustrated for two lithium transition metal oxides with relevance for Li-ion battery cathodes, i.e., Li 2x Fe 2(1-x) O 2 and Li 2x Ti 2(1-x) O 2 , for which the construction of cluster expansion models with compressive sensing alone has proven to be challenging. We demonstrate that our method not only guarantees ground-state preservation on the set of reference structures used for the model construction, but also show that out-of-sample ground-state preservation up to relatively large supercell size is achievable through a rapidly converging iterative refinement. This method provides a general tool for building robust, compressed and constrained physical models with predictive power.

36 MATERIALS SCIENCE↗

Diagnostics of plasma produced by femtosecond laser pulse impact upon a target with an internal nanostructure

X-ray diagnostics of the interaction of femtosecond laser pulses with intensities of 10{sup 16}-10{sup 18} W/cm{sup 2} with CO{sub 2} clusters and frozen nanosize water particles is carried out. The stage of cluster expansion and the formation of a plasma channel, which governs the parameters of the formed X-ray radiation source and accelerated ion flows, is studied. The measurements are based on recording spatially resolved X-ray spectra of H- and He-like oxygen ions. Utilization of Rydberg transitions for spectra diagnostics makes it possible to determine plasma parameters on a time scale of t {approx} 10 ps after the beginning of a femtosecond pulse. The role of the rear edge of the laser pulse in sustaining the plasma temperature at a level of {approx}100 eV in the stage of a nonadiabatic cluster expansion is shown. The analysis of the profiles and relative intensities of spectral lines allows one to determine the temperature and density of plasma electrons and distinguish the populations of 'thermal' ions and ions that are accelerated up to energies of a few tens of kiloelectronvolts. It is shown that the use of solid clusters made of frozen nanoscale water droplets as targets leads to a substantial increase in the number of fast He-like ions. In this case, however, the efficiency of acceleration of H-like ions does not increase, because the time of their ionization in plasma exceeds the time of cluster expansion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Parallel kinetic Monte Carlo simulation framework incorporating accurate models of adsorbate lateral interactions

Ab initio kinetic Monte Carlo (KMC) simulations have been successfully applied for over two decades to elucidate the underlying physico-chemical phenomena on the surfaces of heterogeneous catalysts. These simulations necessitate detailed knowledge of the kinetics of elementary reactions constituting the reaction mechanism, and the energetics of the species participating in the chemistry. The information about the energetics is encoded in the formation energies of gas and surface-bound species, and the lateral interactions between adsorbates on the catalytic surface, which can be modeled at different levels of detail. The majority of previous works accounted for only pairwise-additive first nearest-neighbor interactions. More recently, cluster-expansion Hamiltonians incorporating long-range interactions and many-body terms have been used for detailed estimations of catalytic rate [C. Wu, D. J. Schmidt, C. Wolverton, and W. F. Schneider, J. Catal. 286, 88 (2012)]. In view of the increasing interest in accurate predictions of catalytic performance, there is a need for general-purpose KMC approaches incorporating detailed cluster expansion models for the adlayer energetics. We have addressed this need by building on the previously introduced graph-theoretical KMC framework, and we have developed Zacros, a FORTRAN2003 KMC package for simulating catalytic chemistries. To tackle the high computational cost in the presence of long-range interactions we introduce parallelization with OpenMP. We further benchmark our framework by simulating a KMC analogue of the NO oxidation system established by Schneider and co-workers [J. Catal. 286, 88 (2012)]. We show that taking into account only first nearest-neighbor interactions may lead to large errors in the prediction of the catalytic rate, whereas for accurate estimates thereof, one needs to include long-range terms in the cluster expansion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dual chain perturbation theory: A new equation of state for polyatomic molecules

In the development of equations of state for polyatomic molecules, thermodynamic perturbation theory (TPT) is widely used to calculate the change in free energy due to chain formation. TPT is a simplification of a more general and exact multi-density cluster expansion for associating fluids. In TPT, all contributions to the cluster expansion which contain chain–chain interactions are neglected. That is, all inter-chain interactions are treated at the reference fluid level. This allows for the summation of the cluster theory in terms of reference system correlation functions only. The resulting theory has been shown to be accurate and has been widely employed as the basis of many engineering equations of state. While highly successful, TPT has many handicaps which result from the neglect of chain–chain contributions. The subject of this document is to move beyond the limitations of TPT and include chain–chain contributions to the equation of state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

First Principles Prediction of Structure, Structure Selectivity, and Thermodynamic Stability under Realistic Conditions

Novel materials are often the enabler for new energy technologies. In ab-initio computational materials science, method are developed to predict the behavior of materials starting from the laws of physics, so that properties can be predicted before compounds have to be synthesized and tested. As such, a virtual materials laboratory can be constructed, saving time and money. The objectives of this program were to develop first-principles theory to predict the structure and thermodynamic stability of materials. Since its inception the program focused on the development of the cluster expansion to deal with the increased complexity of complex oxides. This research led to the incorporation of vibrational degrees of freedom in ab-initio thermodynamics, developed methods for multi-component cluster expansions, included the explicit configurational degrees of freedom of localized electrons, developed the formalism for stability in aqueous environments, and culminated in the first ever approach to produce exact ground state predictions of the cluster expansion. Many of these methods have been disseminated to the larger theory community through the Materials Project, pymatgen software, or individual codes. We summarize three of the main accomplishments.

36 MATERIALS SCIENCE↗

Ab initio based examination of the kinetics and thermodynamics of oxygen in Fe-Cr alloys

The behavior of oxygen in Fe-Cr binary alloys is important to understand for a variety of applications including fuel cell interconnects and nuclear energy systems. In this work, we performed an ab initio investigation of the binding of oxygen in the entire compositional range of Fe-Cr at minima and saddles. This database was subsequently used to parametrize concentration-dependent local cluster expansions which were utilized to examine the kinetics of oxygen transport in the alloy and the thermodynamics of the system. The behavior of oxygen in the alloy system was investigated using our cluster expansion model and a convex hull diagram was recorded. The most favorable composition for the incorporation of oxygen was found to occur at 70% Cr. The dependence of the oxygen solution energy on the Cr occupation of the nearest neighbor shells was carefully characterized. It was determined that the third nearest neighbor Cr occupation was unfavorable and that, generally speaking, the first and second nearest neighbor shell Cr occupation was favored with a few exceptions. Kinetic Monte Carlo models at dilute oxygen levels were performed using two models [a kinetically resolved activation (kRA) model and a cluster expansion model for the saddle points]. The models qualitatively agreed with the observed trends in the literature which reported a decrease in diffusivity at dilute Cr levels, but the models predicted different behaviors beyond the dilute Cr limit. Finally, using grand canonical Monte Carlo, we further examined the thermodynamics of oxygen incorporation into the alloy. Furthermore, our results offer new insight into the behavior of oxygen in Fe-Cr alloys that have ramifications on the early stages of the corrosive behavior of the material.

36 MATERIALS SCIENCE↗

Generalized constructive tree weights

The Loop Vertex Expansion (LVE) is a quantum field theory (QFT) method which explicitly computes the Borel sum of Feynman perturbation series. This LVE relies in a crucial way on symmetric tree weights which define a measure on the set of spanning trees of any connected graph. In this paper we generalize this method by defining new tree weights. They depend on the choice of a partition of a set of vertices of the graph, and when the partition is non-trivial, they are no longer symmetric under permutation of vertices. Nevertheless we prove they have the required positivity property to lead to a convergent LVE; in fact we formulate this positivity property precisely for the first time. Our generalized tree weights are inspired by the Brydges-Battle-Federbush work on cluster expansions and could be particularly suited to the computation of connected functions in QFT. Several concrete examples are explicitly given.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Time-resolved x-ray imaging of a laser-induced nanoplasma and its neutral residuals

The evolution of individual, large gas-phase xenon clusters, turned into a nanoplasma by a high power infrared laser pulse, is tracked from femtoseconds up to nanoseconds after laser excitation via coherent diffractive imaging, using ultra-short soft x-ray free electron laser pulses. A decline of scattering signal at high detection angles with increasing time delay indicates a softening of the cluster surface. Here we demonstrate, for the first time a representative speckle pattern of a new stage of cluster expansion for xenon clusters after a nanosecond irradiation. The analysis of the measured average speckle size and the envelope of the intensity distribution reveals a mean cluster size and length scale of internal density fluctuations. Furthermore, the measured diffraction patterns were reproduced by scattering simulations which assumed that the cluster expands with pronounced internal density fluctuations hundreds of picoseconds after excitation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An exact, finite, gauge-invariant, non-perturbative approach to QCD renormalization

A particular choice of renormalization, within the simplifications provided by the non-perturbative property of Effective Locality, leads to a completely finite, non-perturbative approach to renormalized QCD, in which all correlation functions can, in principle, be defined and calculated. In this Model of renormalization, only the Bundle chain-Graphs of the cluster expansion are non-zero. All Bundle graphs connecting to closed quark loops of whatever complexity, and attached to a single quark line, provided no ‘self-energy’ to that quark line, and hence no effective renormalization. However, the exchange of momentum between one quark line and another, involves only the cluster-expansion’s chain graphs, and yields a set of contributions which can be summed and provide a finite color-charge renormalization that can be incorporated into all other QCD processes. An application to High Energy elastic pp scattering is now underway.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A multi-scale study of the adsorption of lanthanum on the (110) surface of tungsten

In this study, we utilize a multi-scale approach to studying lanthanum adsorption on the (110) plane of tungsten. The energy of the system is described from density functional theory calculations within the framework of the cluster expansion method. It is found that including two-body figures up to the sixth nearest neighbor yielded a reasonable agreement with density functional theory calculations as evidenced by the reported cross validation score. The results indicate that the interaction between the adsorbate atoms in the adlayer is important and cannot be ignored. The parameterized cluster expansion expression is used in a lattice gas Monte Carlo simulation in the grand canonical ensemble at 773 K and the adsorption isotherm is recorded. Implications of the obtained results for the pyroprocessing application are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Modeling magnetic evolution and exchange hardening in disordered magnets: The example of Mn 1 - x Fe x Ru 2 Sn Heusler alloys

We demonstrate how exchange hardening can arise in a chemically disordered solid solution from a first-principles statistical mechanics approach. A general mixed-basis chemical and magnetic cluster expansion has been developed, and applied to the Mn 1-x Fe x Ru 2 Sn Heusler alloy system; single-phase solid solutions between antiferromagnetic MnRu 2 Sn and ferromagnetic FeRu 2 Sn with disorder on the Mn/Fe sublattice that exhibit unexpected exchange hardening. Monte Carlo simulations applied to the cluster expansion are able to reproduce the experimentally measured magnetic transition temperatures and the bulk magnetization as a function of composition. The magnetic ordering around a site is shown to be dependent not only on bulk composition, but also on the identity of the site and the local composition around that site. The simulations predict that local antiferromagnetic orderings form inside a bulk ferromagnetic region at intermediate compositions that drives the exchange hardening. Furthermore, the antiferromagnetic regions disorder at a lower temperature than the ferromagnetic regions, providing an atomistic explanation for the experimentally observed decrease in exchange hardening with increasing temperature. These effects occur on a length scale too small to be resolved with previously used characterization techniques.

36 MATERIALS SCIENCE↗

First principles calculation of phase diagrams of V-Nb, V-Ta and Nb-Ta alloys

We report the solid state phase diagram of V-Nb, V-Ta and Nb-Ta alloys computed by combining the density functional theory total energies with the cluster expansion and Monte Carlo techniques. From the computed phase diagrams, we find that V-Nb and Nb-Ta alloys form continuous series of solid solutions and the solid solution phase is stable down to ambient temperatures, consistent with experiments. The calculated bcc V-Ta phase diagram exhibits complete miscibility. Since the current cluster expansion ignore V{sub 2}Ta phase, the chemical interaction due to relatively large electronegativity difference, which cause the ordering of V{sub 2}Ta phase from the bcc solid solution, appears to manifest by making the solid solution phase remain stable for the complete concentration range, down to ambient temperatures, perhaps with some short-range-order. This work further demonstrates the dominant role of constituent strains in the accurate calculation of phase diagram of alloys of constituents with significant size mismatches.

36 MATERIALS SCIENCE↗

Modeling Intercalation Chemistry with Multiredox Reactions by Sparse Lattice Models in Disordered Rocksalt Cathodes

Modern battery materials can contain many elements with substantial site disorder, and their configurational state has been shown to be critical for their performance. The intercalation voltage profile is a critical parameter to evaluate the performance of energy storage. The application of commonly used cluster expansion techniques to model the intercalation thermodynamics of such systems ab initio is challenged by the combinatorial increase in configurational degrees of freedom as the number of species grows. Such challenges necessitate the efficient generation of lattice models without overfitting and proper sampling of the configurational space under the requirement of charge balance in ionic systems. In this work, we introduce a combined approach that addresses these challenges by (1) constructing a robust cluster expansion Hamiltonian using the sparse regression technique, including -norm regularization and structural hierarchy; and (2) implementing semigrand-canonical Monte Carlo to sample charge-balanced ionic configurations using the table-exchange method and an ensemble average approach. These techniques are applied to a disordered rocksalt oxyfluoride (LMNOF) that is part of a family of promising earth-abundant cathode materials. The simulated voltage profile is found to be in good agreement with experimental data and particularly provides a clear demonstration of the and oxygen contributions to the redox potential as a function of content.

25 ENERGY STORAGE↗

Self-organization of S adatoms on Au(111): √3R30° rows at low coverage

Using scanning tunneling microscopy, we observe an adlayer structure that is dominated by short rows of S atoms, on unreconstructed regions of a Au(111) surface. This structure forms upon adsorption of low S coverage (less than 0.1 monolayer) on a fully reconstructed clean surface at 300 K, then cooling to 5 K for observation. The rows adopt one of three orientations that are rotated by 30° from the close-packed directions of the Au(111) substrate, and adjacent S atoms in the rows are separated by √3 times the surface lattice constant, a. Monte Carlo simulations are performed on lattice-gas models, derived using a limited cluster expansion based on density functional theory energetics. Models which include long-range pairwise interactions (extending to 5a), plus selected trio interactions, successfully reproduce the linear rows of S atoms at reasonable temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ab initio investigation of Ti{sub 2}Al(C,N) solid solutions

M{sub n+1}AX{sub n} phases (M: early transition metal, A: IIIA- or IVA-group element, X: carbon or nitrogen) are layered ternary compounds that possess both metal- and ceramic-like properties with numerous potential applications in bulk and thin film forms, particularly under high-temperature conditions. In this work, we use the cluster expansion formalism to investigate the energetics of C-N interactions across the entire Ti{sub 2}AlC-Ti{sub 2}AlN composition range. It is shown that there is a definite tendency for ordering in the C,N sublattice. However, the molar volume and bulk modulus of the ordered structures found along the Ti{sub 2}AlC-Ti{sub 2}AlN composition range show small deviations from the (linear) rule of mixing, indicating that despite the ordering tendencies, the C-N interactions are not strong and the solution becomes disordered at relatively low temperatures. Random solid solutions of Ti{sub 2}AlC{sub 1-x}N{sub x} are simulated using special quasirandom structures (SQS) with x=0.25, 0.50, and 0.75. The thermodynamic properties of these structures are compared to those of the structures found to belong to the ground state through the cluster expansion approach. It is found that the structural properties of these approximations to random alloys do not deviate significantly from Vegard's law. The trend in the structural parameters of these SQS are found to agree well with available experimental data and the predictions of the bulk modulus suggest a very weak alloying effect--with respect to Vegard's law--on the elastic properties of Ti{sub 2}AlC{sub 1-x}N{sub x}.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Finite-Temperature Structures of Supported Subnanometer Catalysts Inferred via Statistical Learning and Genetic Algorithm-Based Optimization

Single-atom catalysts (SACs) minimize noble metal utilization and can alter the activity and selectivity of supported metal nanoparticles. However, the morphology of active centers, including single atoms and subnanometer clusters of a few atoms, remains elusive due to experimental challenges. The computational cost to describe numerous cluster shapes and sizes makes direct first-principles calculations impractical. We present a computational framework to enable structure determination for single-atom and subnanometer cluster catalysts. As a case study, we obtained the low energy structures of Pd n (n = 1-21) clusters supported on CeO 2 (111), which are critical components of automobile three-way catalysts. Trained on density functional theory data, a three-dimensional cluster expansion is established using statistical learning to describe the Hamiltonian and predict energies of supported Pdn clusters of any structure. Low energy stable and metastable structures are identified using a Metropolis Monte Carlo-based genetic algorithm in the canonical ensemble at 300 K. We observe that supported single atoms sinter to form bilayer clusters and large cluster isomers share similarities in both shape and energy, and elucidate the significance of the support and microstructure on cluster stability. We discovered a simple surrogate structure-energy model, where the energy per atom scales with the square root of the average first coordination number, which can be used to estimate energies and compare the stability of clusters. Our framework, applicable to any metal/support system, fills an important methodological gap to predict the stability of supported metal catalysts in the subnanometer regime.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quaternary phase diagrams of spinel Li$_y\square_{1–y}$Mn$_x$Ni$_{2–x}$O$_4$ and composite cathode voltages for concentration gradient materials

Core-shell coating structures and concentration gradient materials may enhance Li-ion battery performance by integrating advantages of core and shell components without introducing unfavorable problems associated with general coatings. The fundamental thermodynamic properties of concentration gradient composite materials are complex due to the multicomponent nature of the problem. Here, we systematically study the thermodynamics of ordering and phase separation in the quaternary spinel Li$_y\square_{1–y}$Mn$_x$Ni$_{2–x}$O$_4$ ($\square$ means vacancy) system by density functional theory calculations, together with the coupled cluster expansion method with interactions within and between (Li/$\square$) and (Mn/Ni) sublattices. On the basis of coupled cluster expansion interactions and Monte Carlo simulations, we calculate quaternary phase diagrams as a function of temperature as well as voltage profiles of single ordered phases and multiphase composite structures. The phase diagram and voltage results are in good agreement with available experimental observations. We also predict a stable high-voltage ordered compound LiMnNiO 4 , with a very high delithiation voltage of 4.76 V. For the composite (Mn-rich + Ni-rich) cathode materials, the voltage profiles show combinations of plateaus from each component compound. The computational strategy of combining quaternary phase diagrams with voltage calculations provides a pathway to understand and design concentration gradient materials.

36 MATERIALS SCIENCE↗