Search NASA⌕ Search

SEARCH · Search NASA

Results for “multi-component 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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

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↗
DOE OSTI2023

A General Framework for Progressive Data Compression and Retrieval

In scientific simulations, observations, and experiments, the transfer of data to and from disk and across networks has become a major bottleneck for data analysis and visualization. Compression techniques have been employed to tackle this challenge, but traditional lossy methods often demand conservative error tolerances to meet the numerical accuracy requirements of both anticipated and unknown data analysis tasks. Progressive data compression and retrieval has emerged as a promising solution, where each analysis task dictates its own accuracy needs. However, few analysis algorithms inherently support progressive data processing, and adapting compression techniques, file formats, client/server frameworks, and APIs to support progressivity can be challenging. Here, this paper presents a framework that enables progressive-precision data queries for any data compressor or numerical representation. Our strategy hinges on a multi-component representation that successively reduces the error between the original and compressed field, allowing each field in the progressive sequence to be expressed as a partial sum of components. We have implemented this approach with four established scientific data compressors and assessed its effectiveness using real-world data sets from the SDRBench collection. The results show that our framework competes in accuracy with the standalone compressors it is based upon. Additionally, (de)compression time is proportional to the number of components requested by the user. Finally, our framework allows for fully lossless compression using lossy compressors when a sufficient number of components are employed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Magnetic contributions to corundum-eskolaite and corundum-hematite phase equilibria: A DFT cluster expansion study

Magnetic contributions have the potential to significantly influence predicted phase stability within alloy and mineral mixing phase diagrams, yet have been historically challenging to incorporate due to a significant increase to phase space sampling. Here, in this work, we employ a computational protocol that includes spin orientation as an additional configurational component within multi-component cluster expansions between magnetic and non-magnetic metal oxide alloys [calculated using density functional theory (DFT) and the generalized gradient approximation]. This approach was used to determine the effect of magnetic contributions to corundum-eskolaite and corundum-hematite phase equilibria from first principles. Two-component cluster expansions of the magnetic components of eskolaite and hematite were first performed showing the ability of this method to properly calculate their respective magnetic properties. Two-component cluster expansions were then performed for non-magnetic Al(III) and ferromagnetic Cr(III) and Fe(III), and phase diagrams were calculated for later comparison. Finally, a non-magnetic Al(III) and “up” and “down” magnetic configurations for anti-ferromagnetic Cr(III) and Fe(III) were performed. Magnetic contributions to the calculated phase diagram for the corundumeskolaite system were shown to be inconsequential but are vital for accurate determination of the corundum-hematite solvus.

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

A generalized and adaptable tensor-contraction-based cluster expansion formalism for multicomponent solids

Density functional theory (DFT)-based simulations of materials have first-principles accuracy, but are very computationally expensive. For simulating various properties of multi-component alloys, the cluster expansion (CE) technique has served as the standard workaround to improve computational efficiency. However, the standard CE technique is difficult to extend to exotic and/or low-symmetry lattices, often implemented via iteration over particular cluster types, which must be enumerated per lattice structure. In this work, we introduce the tensor cluster expansion (TCE), implemented in the open-source code tce-lib, which maps correlation functions to mixed tensor contractions, eliminating the need to iterate over cluster types and additionally making the calculation of correlation functions well-suited for massively parallel architectures like GPUs. We show that local interaction energies are an immediate consequence of the TCE formalism, yielding nearly $\mathcal{O}$(1) energy difference calculations. We then use this formalism to fit CE models for the TaW and CoNiCrFeMn systems, and use these models to respectively compute the enthalpy of mixing curve and Cowley short-range order parameters, showing excellent agreement with ground truth data.

Cluster expansion↗
DOE OSTI2025

Accelerating the discovery of low-energy structure configurations: A computational approach that integrates first-principles calculations, Monte Carlo sampling, and Machine Learning

Finding Minimum Energy Configurations (MECs) is essential in fields such as physics, chemistry, and materials science, as they represent the most stable states of the systems. In particular, identifying such MECs in multi-component alloys considered candidate PFMs is key because it determines the most stable arrangement of atoms within the alloy, directly influencing its phase stability, structural integrity, and thermo-mechanical properties. However, since the search space grows exponentially with the number of atoms considered, obtaining such MECs using computationally expensive first-principles DFT calculations often results in a cumbersome task. To escape the above compromise between physical fidelity and computational efficiency, we have developed a novel physics-based data-driven approach that combines Monte Carlo sampling, first-principles DFT calculations, and Machine Learning to accelerate the discovery of MECs in multi-component alloys. More specifically, we have leveraged well-established Cluster Expansion (CE) techniques with Local Outlier Factor models to establish strategies that enhance the reliability of the CE method. In this work, we demonstrated the capabilities of the proposed approach for the particular case of a tungsten-based quaternary high-entropy alloy. However, the method is applicable to other types of alloys and enables a wide range of applications.

36 MATERIALS SCIENCE↗
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 OSTI2016

Efficient Ab initio Modeling of Random Multicomponent Alloys

Here, we present in this Letter a novel small set of ordered structures (SSOS) method that allows extremely efficient ab initio modeling of random multi-component alloys. Using inverse II-III spinel oxides and equiatomic quinary bcc (so-called high entropy) alloys as examples, we also demonstrate that a SSOS can achieve the same accuracy as a large supercell or a well-converged cluster expansion, but with significantly reduced computational cost. In particular, because of this efficiency, a large number of quinary alloy compositions can be quickly screened, leading to the identification of several new possible high entropy alloy chemistries. Furthermore, the SSOS method developed here can be broadly useful for the rapid computational design of multi-component materials, especially those with a large number of alloying elements, a challenging problem for other approaches.

36 MATERIALS SCIENCE↗
DOE OSTI2017

The Electromagnetic Counterpart of the Binary Neutron Star Merger LIGO/Virgo GW170817. II. UV, Optical, and Near-infrared Light Curves and Comparison to Kilonova Models

We present UV, optical, and near-infrared (NIR) photometry of the first electromagnetic counterpart to a gravitational wave source from Advanced Laser Interferometer Gravitational-wave Observatory (LIGO)/Virgo, the binary neutron star merger GW170817. Our data set extends from the discovery of the optical counterpart at 0.47–18.5 days post-merger, and includes observations with the Dark Energy Camera (DECam), Gemini-South/FLAMINGOS-2 (GS/F2), and the Hubble Space Telescope (HST). The spectral energy distribution (SED) inferred from this photometry at 0.6 days is well described by a blackbody model with T≈8300 K, a radius of R≈4.5×10{sup 14} cm (corresponding to an expansion velocity of v≈0.3c), and a bolometric luminosity of L{sub bol}≈5×10{sup 41} erg s{sup −1}. At 1.5 days we find a multi-component SED across the optical and NIR, and subsequently we observe rapid fading in the UV and blue optical bands and significant reddening of the optical/NIR colors. Modeling the entire data set, we find that models with heating from radioactive decay of {sup 56}Ni, or those with only a single component of opacity from r-process elements, fail to capture the rapid optical decline and red optical/NIR colors. Instead, models with two components consistent with lanthanide-poor and lanthanide-rich ejecta provide a good fit to the data; the resulting “blue” component has M{sub ej}{sup blue}≈0.01 M{sub ⊙} and v{sub ej}{sup blue}≈0.3 c, and the “red” component has M{sub ej}{sup red}≈0.04 M{sub ⊙} and v{sub ej}{sup red}≈0.1 c. These ejecta masses are broadly consistent with the estimated r-process production rate required to explain the Milky Way r-process abundances, providing the first evidence that binary neutron star (BNS) mergers can be a dominant site of r-process enrichment.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Ab initio property predictions of quinary solid solutions using small binary cells

The Set of Small Ordered Structures (SSOS) approach is an ab initio technique for modelling random solid solutions in which many small structures are averaged so that their correlation functions match those of a desired composition. SSOS has been shown to be effective in reducing the cost of density functional theory calculations relative to other well-known techniques such as cluster expansions and special quasirandom structures for modelling solid solutions. Here in this work, we demonstrate that SSOS’s can be constructed using cells with only a subset of elements while still accurately modelling multi-component systems. Specifically, we show that small binary cells can effectively model two quinary high entropy alloys – NbTaTiHfZr and MoNbTaVW – accurately capturing properties such as formation energy, lattice parameters, elastic constants, and root-mean-square atomic displacements. Overall, this insight is useful for those looking to construct databases of such small structures for predicting the properties of multi-component solid solutions, as it greatly decreases the number of structures that needs to be considered.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Intrinsic fluctuations of dust grain charge in multi-component plasmas

A master equation is formulated to model the states of the grain charge in a general multi-component plasma, where there are electrons and various kinds of positive or negative ions that are singly or multiply charged. A Fokker-Planck equation is developed from the master equation through the system-size expansion method. The Fokker-Planck equation has a Gaussian solution with a mean and variance governed by two initial-value differential equations involving the rates of the attachment of ions and electrons to the dust grain. Also, a Langevin equation and a discrete stochastic method are developed to model the time variation of the grain charge. Grain charging in a plasma containing electrons, protons, and alpha particles with Maxwellian distributions is considered as an example problem. The Gaussian solution is in very good agreement with the master equation solution numerically obtained for this problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Reheating signature in the gravitational wave spectrum from self-ordering scalar fields

We investigate the imprint of reheating on the gravitational wave spectrum produced by self-ordering of multi-component scalar fields after a global phase transition. The equation of state of the Universe during reheating, which usually has different behaviour from that of a radiation-dominated Universe, affects the evolution of gravitational waves through the Hubble expansion term in the equations of motion. This gives rise to a different power-law behavior of frequency in the gravitational wave spectrum. The reheating history is therefore imprinted in the shape of the spectrum. We perform 512{sup 3} lattice simulations to investigate how the ordering scalar field reacts to the change of the Hubble expansion and how the reheating effect arises in the spectrum. We also compare the result with inflation-produced gravitational waves, which has a similar spectral shape, and discuss whether it is possible to distinguish the origin between inflation and global phase transition by detecting the shape with future direct detection gravitational wave experiments such as DECIGO.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Numerical study of multi-component flow and mixing in a scaled fission product venting system

Numerical simulation of isothermal multi-component gas flow was conducted to study fluid flow and mixing in a scaled low-pressure, low-temperature test facility for the fission product venting system (FPVS) of a gas cooled fast reactor (GFR). The geometry of FPVS test facility was an open loop including gradual expansion coupling and two 90°pipe elbows. First, the validation of large eddy simulation (LES) wall-adapting local eddy-viscosity (WALE) turbulent model was performed for transitional flow in a circular pipe with gradual expansion. The reactingFOAM solver in OpenFOAM v8 was employed. By introducing appropriate turbulence disturbance at the flow inlet, such as turbulence intensity, the numerical results showed good agreement in (1) the velocity profile downstream of the pipe expansion measured by particle image velocimetry (PIV) in this study, and (2) the reattachment length reported in the literature. Using this validated model, the numerical simulation of flow and component mixing was performed for the FPVS with an inlet flowrate corresponding to a Reynolds number of 2,400 to investigate the flow behavior and component mixing. Standard deviation of the mass fraction of component, referred to as absolute mixing index (AMI), was calculated to quantify the mixing of components, showing that the component mixing length determined by AMI is related to the reattachment length downstream of the expansion. Additionally, this study was able to identify the best location in the FPVS test facility to measure the velocity and concentration profies where multicomponent flow is well mixed.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Thin film combinatorial sputtering of Al-Ce alloys: Investigating the phase separation of as-deposited solid solutions and determining the coefficient of thermal expansion

Al x Ce 100–x thin films with a composition range of ~75.0 < x < 99.5 at% (36.5 < x < 97.5 wt%) were synthesized via combinatorial co-sputtering from an Al and an Al 50 Ce 50 target. The crystal structure, phase fraction, film morphology, electrical resistivity, and temperature-dependent coefficients of thermal expansion (CTE) are all correlated to the Al x Ce 100–x composition. The as-deposited films form a metastable solid-solution, and annealing leads to the formation of the thermodynamically stable two-phase system of Al and the α-Al 11 Ce 3 intermetallic. Temperature dependent x-ray diffraction (XRD) reveals that the two phases expand independently of one another, and the thin film Al temperature-dependent CTE is similar to bulk Al. The thin film Al 11 Ce 3 intermetallic phase has a nearly constant CTE of ~1.5 × 10 –5 /°C within the temperature range studied (25–550 °C). To confirm the thin film Al 11 Ce 3 results, bulk stoichiometric Al 11 Ce 3 and +/- 1 wt% Ce samples were prepared and the CTE of each was measured with the same conditions. A Rietveld analysis of the bulk data enabled an estimation of the CTE in each of the 3 orthorhombic lattice parameters, which displayed anisotropic behavior. The thin film and bulk CTE measurements were in very good agreement. Estimations of the temperature dependent CTE of the two-phase alloys are made via the Reuss and Voigt models. Finally, by demonstrating the efficacy of the approach, more complex multi-component rapid materials discovery of low CTE Al-alloys can be pursued via the combinatorial thin film synthesis and XRD measurement.

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 OSTI2017

High-temperature crystallization of nanocrystals into three-dimensional superlattices

Crystallization of colloidal nanocrystals into superlattices represents a practical bottom-up process with which to create ordered metamaterials with emergent functionalities. With precise control over the size, shape and composition of individual nanocrystals, various single-and multi-component nanocrystal superlattices have been produced, the lattice structures and chemical compositions of which can be accurately engineered. Nanocrystal superlattices are typically prepared by carefully controlling the assembly process through solvent evaporation or destabilization or through DNA-guided crystallization. Slow solvent evaporation or cooling of nanocrystal solutions (over hours or days) is the key element for successful crystallization processes. Here we report the rapid growth (seconds) of micrometre-sized, face-centred-cubic, three-dimensional nanocrystal superlattices during colloidal synthesis at high temperatures (more than 230 degrees Celsius). Using in situ small-angle X-ray scattering, we observe continuous growth of individual nanocrystals within the lattices, which results in simultaneous lattice expansion and fine nanocrystal size control due to the superlattice templates. Thermodynamic models demonstrate that balanced attractive and repulsive interparticle interactions dictated by the ligand coverage on nanocrystal surfaces and nanocrystal core size are responsible for the crystallization process. The interparticle interactions can also be controlled to form different superlattice structures, such as hexagonal close-packed lattices. In conclusion, the rational assembly of various nanocrystal systems into novel materials is thus facilitated for both fundamental research and for practical applications in the fields of magnetics, electronics and catalysis.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Accurate calculation of many-body energies in water clusters using a classical geometry-dependent induction model

We incorporate geometry-dependent distributed multipole and polarizability surfaces into an induction model that is used to describe the 3- and 4-body terms of the interaction between water molecules. The expansion is carried out up to hexadecapole with the multipoles distributed on the atom sites. Dipole-dipole, dipole-quadrupole, and quadrupole-quadrupole distributed polarizabilities are used to represent the response of the multipoles to an electric field. We compare the model against two large databases consisting of 43,844 3-body terms and 3,603 4-body terms obtained from high level ab initio calculations previously used to fit the MB-pol and q-AQUA interaction potentials. The classical induction model with no adjustable parameters reproduces the ab-initio 3- and 4-body terms contained in these teo Databases with a Root-Mean-Square-Error (RMSE) of 0.104/0.058 and a Mean-Absolute-Error (MAE) of 0.054/0.026 kcal/mol, respectively, results that are on a par with those obtained 1 by fitting the same data using tens of thousands of Permutationally Invariant Polynomials (PIPs). This demonstrates the accuracy of this physically motivated model in describing the 3- and 4-body terms in the interactions between water molecules with no adjustable parameters. The triple-dipole-dispersion energy was included in the 3-body energy and was found to be small but not quite negligible. The model represents a practical, efficient and transferable approach for obtaining accurate non-additive interactions for multi-component systems without the need of performing tens of thousands of high level electronic structure calculations and fitting them with tens of thousands of PIPs.

molecular analysis water, aqueous↗
DOE OSTI2020

The Novel Hybrid Ab Initio Model of High-Performance Structural Alloys Design for Fossil Energy Power Plants

The current project developed a novel HT-CALPHAD/DFT approach, which can quickly design new high-performance structural alloys for the application of FE power plants. The PI will mainly take charge of high-throughput DFT simulations and computational thermodynamics of the selected multicomponent alloy systems for the FE power plant applications. At the end of the project, a novel hybrid model based on high throughput CALPHAD/DFT simulations and computational thermodynamics will be developed to provide guidance on how to identify multi-component new high-performance structural alloys with much less computational effort needed. It will address the extensive computation time needed for DFT on the new alloys design. In addition, it will also address the well-known headache of DFT, i.e. how to make the accurate prediction of the high-temperature equilibria. The novel hybrid model the PI proposed will not only be applied to the design of high-performance structural alloys in FE power plants but in many different applications, such as nuclear reactors. This hybrid modeling approach includes four sections: 1) HT-CALPHAD modeling of Al-Co-Cr-Ni-Fe system with FCC and BCC phase. In this section, 3561 non-equiatomic compositions were randomly generated in order to investigate the phase stability of single FCC and BCC phases. Meanwhile, we proposed a data screening procedure to screen out the good candidates within these compositions, considering the temperature range, average density, and melting temperature, etc. 2) Investigation of FCC-Cr lattice stability in Fe-Cr system. We systematically assessed the reliability of FCC-Cr lattice stability derived by DFT and CALPHAD approaches. Meanwhile, the Fe-Cr binary system was chosen as a case study to verify the Cr lattice stability obtained by both approaches. 3) High-throughput DFT modeling on elastic properties of Al-Co-Cr-Ni-Fe systems. We predicted and established the FCC quinary elastic constant database of the Al-Co-Cr-Fe-Ni systems at 0K by using special quasi-random structure (SQS) approach. The predictions will start with pure elements of Al-Co-Cr-Fe-Ni system and will be continued with binaries, ternaries, quaternaries, and finally the quinary compositions. 3) Modeling of temperature-dependent elastic properties in Al-Co-Cr-Ni-Fe systems. In this part, we predict the thermal expansion coefficient and elastic stiffness coefficient as a function of temperature by applying quasiharmonic approximation. With this approach, the elastic properties of HEAs at elevate temperature can be estimated.

01 COAL, LIGNITE, AND PEAT↗