## Uncertainty Propagation in a Multiscale CALPHAD-Reinforced Elastochemical Phase-field Model

Vahid Attari a ; GLYPH&lt;3&gt; , Pejman Honarmandi , Thien Duong , Daniel J Sauceda a a a , Douglas Allaire b , Raymundo Arroyave a b ;

a Materials Science and Engineering Department, Texas A&amp;M University, College Station, TX, USA 77840

b Mechanical Engineering Department, Texas A&amp;M University, College Station, TX, USA 77840

## Abstract

ICME approaches provide decision support for materials design by establishing quantitative process-structure-property relations. Confidence in the decision support, however, must be achieved by establishing uncertainty bounds in ICME model chains. The quantification and propagation of uncertainty in computational materials science, however, remains a rather unexplored aspect of computational materials science approaches. Moreover, traditional uncertainty propagation frameworks tend to be limited in cases with computationally expensive simulations. A rather common and important model chain is that of CALPHAD-based thermodynamic models of phase stability coupled to phase field models for microstructure evolution. Propagation of uncertainty in these cases is challenging not only due to the sheer computational cost of the simulations but also because of the high dimensionality of the input space. In this work, we present a framework for the quantification and propagation of uncertainty in a CALPHADbased elasto-chemical phase field model. We motivate our work by investigating the microstructure evolution in Mg (Si Sn 2 x 1 GLYPH&lt;0&gt; x ) thermoelectric materials. We first carry out a Markov Chain Monte Carlo-based inference of the CALPHAD model parameters for this pseudobinary system and then use advanced sampling schemes to propagate uncertainties across a high-dimensional simulation input space. Through high-throughput phase field simulations we generate 200,000 time series of synthetic microstructures and use machine learning approaches to understand the effects of propagated uncertainties on the microstructure landscape of the system under study. The microstructure dataset has been curated in the Open Phase-field Microstructure Database (OPMD), available at http://microstructures.net.

Keywords: Phase-field modeling, Uncertainty propagation, Uncertainty quantification, Thermoelectrics, Microstructure, Mass scattering, Phonon scattering

## 1. Introduction

Uncertainty Quantification (UQ) has a long and successful history of application to very diverse areas such as climate change [1], structural engineering [2], aerospace engineering and design [3], and medicine [4], to name a few. In the field of materials science, however, notions of UQ remain relatively unexplored even though it is of critical importance as the field progresses towards more quantitative/predictive approaches to materials development. Indeed, uncertainty quantification (UQ) and its propagation (UP) across model/simulation chains are key elements of decision-based [5-8] materials design in the framework of Integrated Computational Materials Engineering (ICME) [9]. The latter prescribes the integration of databases, multi-scale modeling and experiments with the aim to reduce the time and effort of the materials development cycle [10]. Given the complexity and computational cost of most materials simulation frameworks, it is necessary to have a systematic approach to quantify uncertainties in the parameters/variables in any system of interest and to propagate these uncertainties to the respective responses of individual or multi-scale systems.

materials and their performance. In that work, UQ of the parameters and UP across the multi-scale constitutive models (i.e. UP from structure to property to performance) were performed through a Bayesian stochastic method and a stochastic projection technique, respectively. Some works for UP across the multiscale modeling for the probabilistic predictions of plastic flow behavior in poly-crystalline materials have been described thoroughly in [14-16]. Recently, Honarmandi et al . [17] demonstrated the use of UQ approaches to the parameterization of thermodynamically rigorous models for the response of NiTi-based shape memory alloys, followed by uncertainty propagation over the model parameter space.

Despite the importance of UQ/UP in multi-scale modeling [11, 12], there are very few works in the literature dealing with uncertainty quantification (UQ) and/or propagation (UP) across multi-scale models in the field of materials science and engineering. Liu et al. [13], for example, focused on the probabilistic prediction of the effective properties in heterogeneous composite

GLYPH&lt;3&gt; Corresponding author email:attari.v@tamu.edu

In computational materials science, thermodynamic assessments using the CALculation of PHAse Diagrams (CALPHAD) method [18] constitutes the basis for a broad range of approaches to materials simulations, including microstructure evolution through phase field modeling [19]. Given the foundational nature of CALPHAD-based descriptions of phases' free energies in any attempt to predict processing-(micro)structure relationships, UQ/UP in CALPHAD [20] play a very important role, although sparse examples in the literature address this. Honarmandi et al . [21] used a Bayesian framework to quantitate and propagate uncertainty in the context of CALculation of PHAseDiagrams(CALPHAD)thermodynamicassessmentsand showed how information fusion approaches [22] can be used to fuse propagated uncertainties from different competing models. Other groups have also demonstrated different frameworks for the quantification and propagation of uncertainty in CALPHAD models [23-25].

Similarly to the case of CALPHAD-based thermodynamic as-

sessments, the application of UQ/UP frameworks to phase-field modeling remains relatively unexplored [26-31]. Koslowski et al . [15] characterized how uncertainties propagate across spatial and temporal scales in a physics-based model of nanocrystalline plasticity of fcc metals, combining molecular dynamics (MD)withphasefielddislocation dynamics (PFDD) simulations. Wang et al . [32] carried out an asymptotic and uncertainty analysis of void formation during irradiation. Leon et al . [33] used subset selection and active subspace techniques to identify dominant parameters in a continuum phase-field poly-domain model for ferroelectric materials. These earlier approaches focused on UQ/UP over a single modeling framework, but B¨ ottger [34] recently demonstrated the propagation of uncertainty across an entire ICME-based model chain.

Across different fields [35], UP is practically implemented through many different approaches, including Monte-Carlo-, local expansion-,functional expansion-, and numerical integrationbased methods. By far, the most basic and common approach to propagating uncertainty through computational models is via Monte Carlo simulation [36]. For expensive computational models, however, the use of sample-based approaches are often computationally prohibitive. While one can use surrogate models to efficiently sample the input/output relationships [37], numerical efficiency often comes at the expense of fidelity. Moreover, such approaches tend to fail in cases in which the model output changes qualitatively (not only quantitatively) in different regions of the input space. Phase-field simulations belong to the class of computational problems where these conventional approaches to UP tend to be ineffective. The challenges associated with the computational expense and the non-regular nature of the output of phase field simulations are augmented by the large dimensionality of the input space. Finally, the complexity of the output of phase-field simulations makes the analysis of UP frameworks challenging.

The major objective of the current work is to show the propagation of statistically quantified uncertainties of thermodynamic parameters to the Gibbs free energy of phases, and equilibrium phase diagram, in combination with the uncertainties of microelastic and kinetic parameters. Consequently, uncertainties in the Gibbs free energies upto microstructural characteristics are studied by using high throughput analyses across a chain of models that includes a CALPHAD, a microelasticity, and a phase field model, consecutively. The uncertainty of thermodynamic parameters in the CALPHAD model is quantified through a Markov Chain Monte Carlo (MCMC) sampling technique in the context of Bayesian statistics, while the uncertainty of microelastic and kinetic parameters is determined through prior knowledge, calculations and/or estimations.

Figure 1 demonstrates the proposed strategy and steps toward developing a framework for propagating the uncertainty across CALPHAD-phase field model chains. The steps are based on determining 1) type of processing conditions, 2) interacting physics (sub-models), 3) Quantities of Interest (QoIs), and 4) microscopic/macroscopic properties of interest. We propose to demonstrate the framework in the investigation of the microstructure evolution of nanostructured Mg (Si Sn 2 x 1 GLYPH&lt;0&gt; x ) thermoelectric (TE) materials[38, 39]. This pseudo-binary system is characterized by a miscibility gap between two isomorphous cubic phases [39]. The tendency of this system to phase-separate or homogenize is taken into account by employing a fully parameterized elasto- chemical phase-field model that accounts for the effect of process conditions on the resulting microstructure configurations.

Since the performance of multi-phase thermoelectric materials is greatly determined by their (multi-scale) microstructure, it is expected that changes in processing schemes may have considerable impact on performance. In fact, for the case of Mg (Si Sn 2 x 1 GLYPH&lt;0&gt; x )-based thermoelectrics, it has been already shown that equilibrium and non-equilibrium processing can lead to dramatic changes to the TE figure of merit, ZT [39-41]. Understanding of this TE system is further hindered by the considerable uncertainty in the location of the miscibility gap, with different experimental phase boundary estimations disagreeing by several tens of atomic percent [39].

In this work, we attempt to propagate uncertainty through phase field simulations, accounting for the uncertainty not only in the parameters directly associated to the phase field model, but also arising from uncertainty in the CALPHAD parameterization. This paper is structured as follows: in section 2, we motivate the present work through the application of ICME-based frameworks on the design of nanostructured TE materials. We note that the framework put forward is generalizable to a wide range of materials problems. The details of the models (CALPHAD and phase-field) and uncertainty propagation strategy is provided in section 3. Section 4 discussed the prior and the resulting uncertainty in the phase diagram, the subsequent uncertainty in the microstructures under elasto-chemical simulations, and the methods for data interpretation/classification. Furthermore, we present a summary on our findings and draw our conclusions in section 5.

## 2. Nanostructured Thermoelectric Materials

As mentioned above, the motivating example is the propagation of uncertainty in microstructure evolution simulations in nanostructured composite thermoelectrics (TE) [42-45] via direct coupling of CALPHAD thermodynamic assessments [46] with multi-physics phase field models (PFM) [28, 39] that account for both chemical and elastic driving forces for structure formation. The example is motivated by recent work by some of the present authors [39] on the dramatic effect that processing has on the microstructure (and TE performance) in Mg Si-Mg Sn alloys, but 2 2 has a much broader applicability as CALPHAD/PFM-basedmicrostructure simulations are pervasive in ICME-based frameworks for microstructure-sensitive materials design [47-51], and properly accounting for uncertainty is necessary to make design choices with proper confidence bounds.

Current interest in the thermoelectric (TE) effect originates from the ever increasing demand for energy and the associated detrimental effects on global climate. Current TE materials, unfortunately, do not have the efficiency-described by the figure of merit z T = GLYPH&lt;11&gt; 2 sT =GLYPH&lt;11&gt; , where GLYPH&lt;11&gt; is the Seebeck coefficient, s is the electrical conductivity, T is the absolute temperature and GLYPH&lt;20&gt; is the thermal conductivity-that would turn TE-based devices into competitive power-generators [42, 52]. An ideal TE material would have a large Seebeck coefficient, while being electrically conductive and thermally insulating [42, 52]. These properties, however, are coupled and their individual tuning is thus challenging. Over the past decade, a sophisticated arsenal of strategies for the rational design of TE materials has emerged [53-55],

Figure 1: The process-structure-property uncertainty propagation framework deployed in composition, strain and temperature space based on the process definition and natural uncertainties in input parameters.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00638v1-figure-f9974ec3-47da-42d5-a780-d22ac2976b11.jpg)

including the exploitation of spontaneous self-assembly or nonequilibrium processing of nanostructures to enhance phononscattering [39, 42].

Among the hundreds of TE systems investigated to date, environmentally-benign Mg (Si Sn 2 x 1 GLYPH&lt;0&gt; x ) alloys [56, 57] have attracted considerable attention due to their relatively high figure of merit ( zT &gt; 1 ) [58], comparable with the intermediate temperature TE materials such as PbTe and filled skutterudites [39, 57, 59, 60]. The Mg (Si Sn 2 x 1 GLYPH&lt;0&gt; x ) pseudo-binary system exhibits a miscibility gap [61-64] and this has been exploited to realize nanostructures with optimal TE performance [56, 58, 65]. Experimental determination of ZT even in a single alloy exhibits considerable variance, perhaps due to changes in the way these materials are synthesized and processed. Recently, the present authors and collaborators investigated the effect of non-equilibrium processing on the microstructure evolution (and transport properties) in the Mg 2 Si 0 7 : Sn 0 3 : system and found that instead of phase-separating, the system tended to form a solidsolution with superior TE performance, contrary to expectations and prior works [56, 65]. This was ascribed to (elastic) coherency effects and was verified via quantitative multi-physics phase field simulations [39].

These results are interesting as they exemplify the influence of processing on the microstructural evolution in TE materials [42] and the corresponding change in performance. Most importantly, this constitutes one of the very few examples-to the best of our knowledge-in which phase-field modeling (PFM), in combination with CALPHAD free energies has been used to investigate the microstructure evolution of TE materials. Further investigation of the PFM developed to study the Mg Si 2 0 7 : Sn 0 3 : showed that rather small changes in the strength of the elastic couplingsmediated via lattice parameter differences between Si- and Snrich domains-resulted in qualitatively different microstructures, which in turn could be expected to exhibit different phonon transport behavior.

The ability to quantitatively understand [42] and control the different materials and processing parameters related to mi- crostructural morphology, topology, size and spacing in composite TE materials has already been demonstrated [45, 66-69]. Many of these approaches have been inspired by metallurgy and thus the time is ripe to translate much of what has been learned on ICME-enabling microstructure-sensitive (structural) alloy design to the problem of designing (self-assembled) TE microstructures for optimal performance. While the modeling framework via PFM has shown to result in (semi-)quantitative predictions that compare well with experiments [39], a robust ICME research program on microstructure design of TE materials requires reliable and efficient UQ/ UP frameworks.

Figure 2 illustrates an schematic phase diagram in which the material shows an inherent instability in certain regions of the composition space. This material is uniform at high temperatures and upon reducing the temperature decomposes into distinct phases. Similarly, Mg X Sn,Si 2 f g system has a miscibility where the uncertainties in the Sn-rich and Si-rich boundaries (Refer to Fig. 3 of [39]) has been the center of discussions in several articles [38, 39, 62, 63, 70]. As a consequence, these uncertainties impact the chemical free energy and (the predicted) microstructure of this system.

The extent of instabilities, the uncertainty in the boundaries of the unstable region, and extensive possibilities by means of using different synthesis conditions all suggest that analysis and quantification of uncertainty is necessary in this material to establish the correlation between the thermodynamics and microstructural phenomena. These analyses are essential to estimate how the variance in processing and (epistemic or aleatoric) uncertainties in material parameters can affect the microstructure and ultimately the TE response of this material. This also can help to better realize the reliable thermodynamic conditions for enhancement of the properties of TE materials [71].

## 3. Models and Methodologies

Here, we first define the thermodynamic state variables relevant to the elasto-chemical phase field model formulation and

Figure 2: A general strategy to enhance thermoelectric performance by engineering better phonon scattering via nanoscale precipitates or solid solution.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00638v1-figure-5f75fb29-4131-4865-bdab-4bf3b5b0ba34.jpg)

then quantify the uncertainties in the given material variables and propagate these uncertainties with the aim of quantifying the variations in the micro-structure of the material that ultimately alters the predicted macroscopic response.

## 3.1. Thermodynamic State Variables

For an isothermal and isobaric state, the total free energy functional ( F tot ) for an undeformed material configuration with confined boundaries can be constructed from the sum of all contributing fields over it. We restrict F tot to chemical and elastic contributing effects here. Consequently, the energy is a function of composition ( c ), temperature ( T ), strain ( " ), gradients of composition ( r c ), and other fields, if they were present. We write the total free energy functional, F tot as the sum of contributing fields:

$$F ^ { t o t } ( c, T, \nabla c, \varepsilon ) = f ^ { b u l k } + f ^ { i n t e r f a c i a l } + f ^ { e l a s } \quad \ \ ( 1 )$$

where the bulk free energy density, f bulk , the interfacial free energy density, f inter f acial , and elastic strain energy density, f elas are:

$$f ^ { b u l k } ( c, T ) = \int _ { \Omega } \left ( f ^ { 0 } ( c, T ) \right ) d \Omega \quad \quad ( 2 ) \quad \text{late} \\ \intertext { a i n t e r f o c i a l } = \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \end{array}$$

$$\bar { f } ^ { i n t e r f a c i a l } ( \nabla c ) \stackrel { J _ { \Omega } } { = } \int _ { \lambda } \kappa ( \nabla c ) ^ { 2 } d \Omega \quad \quad ( 3 ) \quad \text{the} \quad.$$

$$\overset { J } { \nu } = \overset { \nu } { \int } _ { \Omega } \overset { J } { \sigma _ { i j } } \varepsilon _ { i j } ^ { e l } d \Omega \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
 \text{$f^{e l a s } } ( \varepsilon ) = \overset { J } { \int } _ { \Omega } \overset { \nu \nu } { \sigma _ { i j } } \overset { \nu \nu } { \varepsilon _ { i j } } d \Omega \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
 \text{$f^{r} } ( \overset { J } { \int } _ { \Omega } \overset { \nu } { \sigma _ { i j j } } \overset { \nu } { \varepsilon _ { i j } } \overset { \nu } { \varepsilon _ { i j } } \overset { \nu } { \varepsilon _ { i j } } \overset { \nu } { \varepsilon _ { i j } } \overset { \nu } { \varepsilon _ { i j } } \overset { \nu } { \varepsilon _ { i j } } \overset { \nu } { \varepsilon _ { i j } } \overset { \nu } { \varepsilon _ { i j } } \overset { \nu } { \varepsilon _ { i j } } \overset { \nu } { \varepsilon _ { i j } } \overset { \nu } { \varepsilon _ { i j } } \overset { \nu \nu } { \varepsilon _ { i j } } \overset { \nu } { \varepsilon _ { i j } } \overset { \nu } { \varepsilon _ { i j } } \overset { \nu } { \varepsilon _ { i j } } \overset { \nu } { \varepsilon _ { i j } } \overset { \nu } { \varepsilon _ { i j } } \overset { \nu } { \varepsilon _ { i j } } \overset { \nu } { \varepsilon _ { i j } } \overset { \nu } { \varepsilon _ { i j } } \overset { \nu } { \varepsilon _ { i j } } \overset { \nu } { \varepsilon_ { i j } } \overset { \nu } { \varepsilon _ { i j } } \overset { \nu } { \varepsilon _ { i j } } \overset { \nu } { \varepsilon _ { i j } } \overset { \nu } { \varepsilon _ { i j } } \overset { \nu } { \varepsilon _ { i j } } \overset { \nu } { \varepsilon _ { i j } } \overset { \nu } { \varepsilon _ { i j } } \overset { \nu } { \varepsilon _ { i j } } \overset { \nu } { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { n } \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \varepsilon _ { i j } } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { n } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu } \overset { \nu } { \nu }$$

where f 0 ( c T ; ) is the free energy of a unit volume of homogeneous material, GLYPH&lt;20&gt; is the gradient energy coefficient, " el i j and GLYPH&lt;27&gt; i j are the local elastic strain and stress in the material, respectively. The chemical free energy is composed of the interfacial and bulk energy contributions, and it determines the compositions and volumefractions of the equilibrium phases. The strain energy affects the equilibrium compositions and volume fractions of the coexisting phases, but also determines the shapes and configurations of the phase domains.

The bulk free energy of the Mg (Si Sn 2 x 1 GLYPH&lt;0&gt; x ) pseudo-binary system is described through the sub-regular solution model as:

$$f ^ { b u l k } ( c _ { \alpha }, T ) = \sum _ { \alpha } c _ { \alpha } G _ { \alpha } ^ { 0 } + R T \sum _ { \alpha } c _ { \alpha } l n ( c _ { \alpha } ) + \sum _ { \alpha } \sum _ { \beta > \alpha } c _ { \alpha } c _ { \beta } \sum _ { \nu } L _ { \alpha \beta } ^ { \nu } ( c _ { \alpha } - c _ { \beta } ) ^ { \nu } \ ( 5 )$$

where c GLYPH&lt;11&gt; is the mole fraction of the constituent GLYPH&lt;11&gt; (either Mg Sn 2 or Mg Si), 2 G 0 GLYPH&lt;11&gt; is the reference energy of the constituent, T is temperature (in K), R is the gas constant, and L GLYPH&lt;23&gt; is given as:

GLYPH&lt;11&gt;GLYPH&lt;12&gt;

$$L _ { \alpha \beta } ^ { \nu } = \nu a _ { \alpha \beta } + \nu b _ { \alpha \beta } T$$

where GLYPH&lt;23&gt; a GLYPH&lt;11&gt;GLYPH&lt;12&gt; and GLYPH&lt;23&gt; b GLYPH&lt;11&gt;GLYPH&lt;12&gt; are model parameters, which describe the interactions between the constituents beyond those of ideal mixing. These parameters can be calibrated against the available data deterministically or probabilistically, as discussed later. In order to build the phase diagram, the total Gibbs free energy of the system is minimized at different temperatures for volume fraction and composition of each existing phase.

## 3.2. Material Modeling Strategy: CALPHAD Reinforced Phasefield Method

From linear kinetic theory, the local mass flux in the presence of a gradient in composition, the diffusion flux, ~ J (in units of mol m GLYPH&lt;0&gt; 2 s GLYPH&lt;0&gt; 1 ) is given by:

$$\vec { J } = - \vec { M } \nabla \mu ^ { t o t }$$

where ~ M is the interface mobility assumed to be constant due to the isotropic nature of the crystal structures of the two phases, and GLYPH&lt;22&gt; tot = GLYPH&lt;14&gt; F tot GLYPH&lt;14&gt; c is the total potential for the kinetic transition. We postulate the following form of the Cahn-Hilliard (C-H) kinetic equation along with other micro-elasticity equations to study the evolution of the Mg (Si Sn 2 x 1 GLYPH&lt;0&gt; x ) microstructure in space:

$$\frac { \partial c } { \partial t } = \nabla \vec { M } \nabla \left ( \frac { \delta F ^ { t o t } } { \delta c } \right )$$

$$\frac { \partial \sigma _ { i j } } { A _ { \nu }. } = 0$$

$$\sigma _ { i j } \stackrel { \bullet } { = } C _ { i j k l } ^ { ^ { \prime } } \varepsilon _ { k l } ^ { e l } \stackrel { \bullet ^ { \prime } } { = } \quad \quad \quad$$

$$\frac { \frac { \dot { \nu } } { \partial r _ { j } } } { \partial r _ { j } } = \frac { 1 } { 2 } \left ( \frac { \partial u _ { i } } { \partial r _ { j } } - \frac { \partial u _ { j } } { \partial r _ { i } } \right ) \quad \ \ ( 1 0 ) \\ \nu \, - \, \tilde { \nu } \dots \tilde { \tilde { \nu } } ^ { e l }$$

$$\varepsilon _ { k l } ^ { e l } = \varepsilon _ { k l } ^ { t o t } - \varepsilon _ { k l } ^ { 0 } \text{ \quad \quad \ \ } ( 1 2 )$$

Equations 8-12 are the C-H equation, mechanical equilibrium condition, kinematics, Hooke's microscopic constitutive law for linear elasticity, and strain relationship, respectively. The dilatational eigenstrain term is given by " 0 kl = " T GLYPH&lt;14&gt; kl h c ( ) , and it is the consequence of lattice strain between the phases. GLYPH&lt;14&gt; is the Kronecker-delta function and h c ( ) is an interpolation function. Cijkl is the composition-dependent fourth order elastic modulus tensor. It is convenient to describe Cijkl using the following expression:

$$C _ { i j k l } ( c ) = C _ { i j k l } ^ { e f f } - g ( c ) \Delta C _ { i j k l } \text{ \quad \ \ } ( 1 3 )$$

where GLYPH&lt;1&gt; Cijkl = C i jkl GLYPH&lt;11&gt; GLYPH&lt;0&gt; C i jkl GLYPH&lt;12&gt; is the difference between the elastic moduli tensor of GLYPH&lt;11&gt; and GLYPH&lt;12&gt; phases. Cijkl is a 6 GLYPH&lt;2&gt; 6 symmetric tensor which can be denoted by the Voigt notation for any linearly elastic medium in the form of:

$$C _ { i j k l } \to C _ { m n } = \left [ \begin{array} { c c c c } C _ { 1 1 } & C _ { 1 2 } & C _ { 1 3 } & C _ { 1 4 } & C _ { 1 5 } & C _ { 1 6 } \\ C _ { 2 2 } & C _ { 2 3 } & C _ { 2 4 } & C _ { 2 5 } & C _ { 2 6 } \\ C _ { 3 3 } & C _ { 3 4 } & C _ { 3 5 } & C _ { 3 6 } \\ & C _ { 4 4 } & C _ { 4 5 } & C _ { 4 6 } \\ \text{Symmetric} & & C _ { 5 5 } & C _ { 5 6 } \\ & & & & C _ { 6 6 } \end{array} \right ] \quad ( 1 4 )$$

Similar to h c ( ) , g c ( ) is a scalar-valued interpolation function such that h c ( GLYPH&lt;11&gt; ) = g c ( GLYPH&lt;11&gt; ) = 0 and h c ( GLYPH&lt;12&gt; ) = g c ( GLYPH&lt;12&gt; ) = 1 . The phasefield microstructural evolution problem (eqn. 8) is solved by utilizing a semi-implicit Fourier spectral approach [72]. The microelasticity problem (eqns. 9-12) is solved by the FFT-based iterative solver described in [29]. In this approach, the following conditions are considered: 1) stress-free transformation strains (SFTS) for each phase and inhomogeneous elastic constants in the domain, 2) GLYPH&lt;1&gt; Cijkl = 0 in the first iteration, 3) strain-control based on stress-control, 4) convergence of the problem when the L 2 normof jj u n + 1 GLYPH&lt;0&gt; u n jj is less than 10 GLYPH&lt;0&gt; 8 , and 5) Periodic boundary conditions in all sides of the domain.

The parameter space for CALPHAD, micro-elasticity, and phase-field models is listed in Table 1.

## 3.3. Uncertainty Quantification/Propagation

The quantification of uncertainties associated with model parameters and predictions is one of the most important tasks in computational aided materials design. In this regard, quantifying the uncertainties of the model parameters given data is an inverse problem, known as inverse UQ, while obtaining the uncertainties of the predictions of either an individual model or a chain of models through the propagation of the parameter uncertainties is a forward problem, known as forward UP.

Generally, UQ of model parameters can be performed in the context of either frequentist or Bayesian inference. In this work, the latter has been applied to probabilistically calibrate the relevant model parameters. Here, the parameter calibration or uncertainty quantification is performed based on Bayes' rule where the parameter posterior probability given data is proportional to the parameter prior probability times likelihood, i.e. P ( GLYPH&lt;18&gt; j D ) / P ( GLYPH&lt;18&gt; ) P D ( j GLYPH&lt;18&gt; ) assuming GLYPH&lt;18&gt; and D are the parameters and data, respectively. Here, the prior probability distribution of the model parameters-which is defined through the available prior knowledge about the parameters-is updated to their posterior probability distribution given data using the likelihood distributionwhich is a conditional probability of data given the model parameters. Therefore, the Bayesian inference results from the combination of prior beliefs and likelihood rather than the inference just based on likelihood, as is the case in frequentist approaches.

Theresulting posterior distribution represents the plausible optimal value of the parameters and their uncertainties that can be propagated to the model predictions through different analytical or numerical approaches. To find the statistical properties (mean value and covariance matrix) of the posterior probability distribution, some intractable integrals should be solved; however, these integrations are very difficult to carry out through the conventional numerical approaches. For this reason, Monte Carlo (MC)-based approaches are often used to estimate these integrals [73, 74]. In order to obtain well-converged estimates for the relevant integrals, a substantial number of evaluations of the model must be made. Generally, Markov Chain Monte Carlo (MCMC) approaches are the most common used sampling techniques in the context of Bayesian inference. Among these techniques, the Metropolis-Hastings (MH) algorithm has been selected in this work, as described in detail in section 3.3.1.

We note that it is also important to propagate uncertainties from the parameters to the model predictions since these prob- abilistic predictions of the given system are the quantities of interest in materials design.

Three commonly used approaches for UP are the analytical methods, the numerical Monte Carlo (MC) methods, and stochastic finite elements. In analytical methods, the output/prediction uncertainties are represented explicitly as functions of parameter uncertainties. These methods are useful when only small ranges of uncertainty are considered, and cannot always be applied to complex, nonlinear models with large parameter uncertainties. Monte Carlo methods involve the estimation of probability density functions (PDFs) for selected model outputs by performing a sufficiently large number of model runs with randomly sampled parameters. For computationally intensive models, the time and resources required for Monte Carlo methods could be prohibitively expensive. A degree of computational efficiency is accomplished by the use of Modified Monte Carlo (MMC) methods that sample from the parameter distribution in an efficient manner, so that the number of necessary solutions are significantly reduced compared to the simple Monte Carlo method. Such methods include the Fourier Amplitude Sensitivity Test (FAST) [75] and Latin Hyper-cube Sampling [76]. However, even these require a substantial number of model runs to obtain a good approximation of the output PDFs, especially for cases involving a large number of uncertain parameters. Therefore, there is a need for even more computationally efficient ways for propagating uncertainty in complex/expensive models and this is currently being investigated by the present authors.

In this work, a forward UP analysis based on optimal sampling methods has been applied to propagate the uncertainties across the chain of models, which is explained in section 3.3.2.

## 3.3.1. Applied MCMC Approach for the Parameter Uncertainty Quantification

The MCMC-MH toolbox in Matlab has been utilized to calibrate the Gibbs free energy parameters of a CALPHAD model for the pseudo-binary Mg (Si Sn 2 x 1 GLYPH&lt;0&gt; x ) system. Since there is no prior knowledge about the distribution of the parameters and their correlations, a uniform prior probability distribution has been considered for the model parameters in this work. However, the initial values ( GLYPH&lt;18&gt; 0 ) and ranges of the parameters have been determined based on the deterministically calibrated parameters in ThermoCalc software in order to make the parameter convergence faster during MCMC sampling. After defining the prior probability distribution, the parameters (candidates) are randomly sampled from an arbitrary posterior proposal distribution ( q ) iteratively. In this work, the proposal distribution is an adaptive multi-variate Gaussian distribution that is centered at the last accepted parameter vector in the MCMC chain with a covariance matrix adapting during MCMC sampling based on the resulting covariance from the previous parameters in the chain, according to Haario et al [77]. In each iteration, the sampled candidate is accepted or rejected based on the Metropolis-Hastings ratio:

$$M - H = \frac { p ( \theta ^ { c a n d } ) p ( D | \theta ^ { c a n d } ) } { p ( \theta ^ { z - 1 } ) p ( D | \theta ^ { z - 1 } ) } \frac { q ( \theta ^ { z - 1 } | \theta ^ { c a n d } ) } { q ( \theta ^ { c a n d } | \theta ^ { z - 1 } ) } \quad \ ( 1 5 )$$

This ratio compares the posterior probability of the sampled candidate ( GLYPH&lt;18&gt; cand ) with its counterpart for the last accepted parameter vector in the MCMC chain ( GLYPH&lt;18&gt; z GLYPH&lt;0&gt; 1 ) through the metropolis ratio

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00638v1-figure-843f31f6-9429-4b7a-9687-6ea7ac3da1e9.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00638v1-figure-19f71a92-57e5-40bc-a451-817656f7d4f0.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00638v1-figure-4fddeb8a-83ee-4966-851d-25b258319c4c.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00638v1-figure-08e1c469-00d9-4b24-9ef6-c47e3cd6c665.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00638v1-figure-608bb3d7-0c24-48bc-a789-89fe63dad5bb.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00638v1-figure-15534387-8dde-468f-b762-b836e34703cf.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00638v1-figure-f6c9c875-f2ca-4a73-af4c-b4d27737888d.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00638v1-figure-e0ac3a99-c2d7-48ad-afbf-67fbae49f85f.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00638v1-figure-83e4da13-2cf1-4285-a3a5-645e891598a3.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00638v1-figure-200c95a4-9bb4-4e9f-8bc5-3c9d3788d3c3.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00638v1-figure-aac3c630-93a3-4ebf-8152-7015eb6edd65.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00638v1-figure-823b1591-5088-4d15-a8f5-6556bd936e18.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00638v1-figure-6a41944a-f0fd-4ab5-acf8-a4e409ea12d0.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00638v1-figure-3c2c0734-2a05-485f-b45c-8587ce70c329.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00638v1-figure-d1811a80-bd5c-40c9-94a3-14fd85174640.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00638v1-figure-51813f8e-d3f4-4448-99af-0e3409d0f46f.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00638v1-figure-6965ce9a-090a-41c0-b515-66296bd25596.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00638v1-figure-410a3301-0dd3-40a7-9a7a-8c5f685c02df.jpg)

Table 1: The quantification of uncertainty in the CALPHAD model and 1 st order statistics on the outputs of inputs of the phase-field model. The outputs of the CALPHAD model were fed to the phase-field model. The fitted truncated PDFs which correspond to the plausible optimal values and uncertainty of the model parameters. 'ss' and 'liq' denote Mg 2 (Si x Sn 1 GLYPH&lt;0&gt; x ) and liquid phases, respectively. ** SFTS stands as Stress-Free Transformation Strain.

| Model       | Sub-model       | Parameter                                 | PDF form   | Unit                            | GLYPH<22>                             | Dispersion ( GLYPH<27> 2 GLYPH<22> )   | Lower Bound                   | Upper Bound                   |
|-------------|-----------------|-------------------------------------------|------------|---------------------------------|---------------------------------------|----------------------------------------|-------------------------------|-------------------------------|
|             |                 | 0 a ss                                    |            | J : mol GLYPH<0> 1              | 12840.4                               | 729.0                                  | 6824.89                       | 20474.69                      |
|             |                 | 0 b ss                                    |            | J : mol GLYPH<0> 1              | 7.20                                  | 0.59                                   | 3.67                          | 11.02                         |
|             | Parameters      | 1 a ss                                    |            | J : mol GLYPH<0> 1              | -3324.3                               | -288.3                                 | -5208.34                      | -1736.11                      |
|             | Model           | 0 a liq                                   |            | J : mol GLYPH<0> 1              | 80635.7                               | 26688.6                                | 43550.19                      | 130650.57                     |
|             |                 | 0 b liq                                   |            | J : mol GLYPH<0> 1              | -61.25                                | 3.98                                   | -86.03                        | -28.68                        |
| Parameters  | CALPHAD         | 1 a liq                                   |            | J : mol GLYPH<0> 1              | 6400.5                                | 529.01                                 | 3314.80                       | 9944.41                       |
| Model       |                 | Alloy composition ( c )                   |            | mol                             | 0.40                                  | 0.0083                                 | 0.3                           | 0.5                           |
|             |                 | Interface mobility ( M )                  |            | m 2 s GLYPH<0> 1 : J GLYPH<0> 1 | 5 : 62 GLYPH<2> 10 GLYPH<0> 19        | 4 : 64 GLYPH<2> 10 GLYPH<0> 25         | 10 GLYPH<0> 18 = ( RT )       | 10 GLYPH<0> 20 = ( RT )       |
| Phase-field | Parameters      | Gradient energy Coefficient ( GLYPH<20> ) |            | J : m GLYPH<0> 2                | 1 : 16 GLYPH<2> 10 GLYPH<0> 24        | 4 : 64 GLYPH<2> 10 GLYPH<0> 25         | 2 : 0 GLYPH<2> 10 GLYPH<0> 26 | 2 : 0 GLYPH<2> 10 GLYPH<0> 24 |
|             | Kinetic         | SFTS ** ( " T )                           |            | -                               | GLYPH<0> 3 : 1 GLYPH<2> 10 GLYPH<0> 5 | -4.28                                  | -0.02                         | +0.02                         |
|             |                 | C 11 Mg 2 Sn                              |            | GPa                             | 76.56                                 | 0.20                                   | 68.30                         | 83.71                         |
|             |                 | C 12 Mg 2 Sn                              |            | GPa                             | 27.75                                 | 1.04                                   | 17.68                         | 39.79                         |
|             | Parameters      | C 44 Mg 2 Sn                              |            | GPa                             | 29.94                                 | 1.26                                   | 16.03                         | 41.94                         |
|             |                 | C 11 Mg 2 Si                              |            | GPa                             | 120.15                                | 0.07                                   | 114.07                        | 126.00                        |
|             | Model           | C 12 Mg 2 Si                              |            | GPa                             | 22.62                                 | 0.10                                   | 19.56                         | 26.00                         |
|             | Microelasticity | C 44 Mg 2 Si                              |            | GPa                             | 46.79                                 | 0.63                                   | 33.32                         | 58.20                         |
|             |                 | Molar volume ( V Mg 2 Si m )              |            | m 3 : mol GLYPH<0> 1            | 5 : 78 GLYPH<2> 10 GLYPH<0> 5         | 2 : 40 GLYPH<2> 10 GLYPH<0> 7          | 4 : 73 GLYPH<2> 10 GLYPH<0> 5 | 6 : 38 GLYPH<2> 10 GLYPH<0> 5 |
|             |                 | Molar volume ( V Mg 2 Sn m )              |            | m 3 : mol GLYPH<0> 1            | 4 : 88 GLYPH<2> 10 GLYPH<0> 5         | 1 : 55 GLYPH<2> 10 GLYPH<0> 7          | 3 : 95 GLYPH<2> 10 GLYPH<0> 5 | 5 : 33 GLYPH<2> 10 GLYPH<0> 5 |

(the first ratio in Equation 15), and also compares the probability of moving from GLYPH&lt;18&gt; z GLYPH&lt;0&gt; 1 to GLYPH&lt;18&gt; cand with the probability of the reverse move through the Hastings ratio (the second ratio in Equation 15) in the case that the posterior proposal distribution is asymmetric. It should be noted that the likelihood function is also a Gaussian distribution centered at the given data ( D ) with a variance that corresponds to the uncertainty of the data. Since the uncertainty of the data used in this work is not clear, an unknown variance has been considered for the likelihood which has been updated as a hyper-parameter with the model parameters during MCMC sampling, which has been explained in more detail by Gelman et al. [73]. Min(M-H,1) indicates the acceptance probability of the candidate in each iteration. GLYPH&lt;18&gt; z equals GLYPH&lt;18&gt; cand in the case of acceptance which updates the mean value of the proposal distribution; while it equals GLYPH&lt;18&gt; z GLYPH&lt;0&gt; 1 if the candidate is rejected. The sampling process proceeds until the convergence of MCMC chain to a stationary distribution, which represents the parameter convergence during MCMCprocess. After discarding the burn-in period which is the initial MCMC samples before the parameter convergence, the mean values of the remaining samples and the square root of the diagonal terms in their covariance matrix represent the plausible optimal values and uncertainties of the parameters, respectively.

## 3.3.2. Applied Uncertainty Propagation Approach

As mentioned earlier, a forward UP analysis has been considered in this work to propagate uncertainties from thermodynamic parameterstoGibbsfreeenergyofphasestomicrostructuralcharacteristics in Mg (Si Sn 2 x 1 GLYPH&lt;0&gt; x ) system. In this regard, a group of sampled parameter vectors has been considered as representative of the multi-variate posterior distribution of the thermodynamic parameters and run through the applied CALPHAD model to find the corresponding responses for Gibbs free energy of phases and phase diagram. To find 95% Bayesian credible intervals (BCI), 2.5% of the samples from both the lower and upper bounds of the response distributions have been removed. The same approach has been used to propagate uncertainties from thermodynamic, microelastic and kinetic parameters to microstructural characteristics throughout a phase field model. However, a Gaussian copula approach has been used to sample a reasonable number of parameter vectors from the distributions of microelastic and kinetic parameters-assuming they are independent-as well as the multi-variate posterior distribution of thermodynamic parameters obtained by the MCMC technique.

## 3.3.3. Sampling Methodology

Here we note that propagating uncertainties from highdimensional input sets has to be carried out in a way that minimizes the number of samples ( a typical naive sampling scheme using MCMC approaches may require O (1 ; 000 000) ; random samples) while at the same time accounting for the statistical correlations among input parameters-for example, parameters in CALPHAD thermodynamic descriptions tend to be highly correlated.

To construct sample sets with correct marginal distributions and preserved pairwise correlations, we instead make use of Gaussian copulas. A copula is a function that relates the joint cumulative distribution function (CDF) of multiple variables to their marginal CDFs and their correlations [78]. To begin, we assume that we have available marginal distributions, f Xi ( xi ) , for each parameter, where Xi denotes the random variable associated with the i th parameter, i 2 f 1 2 ; ; : : : ; K g , and xi is a particular realization of Xi . Wealso have available pairwise correlation coefficients, GLYPH&lt;26&gt; i ; j , where i ; j 2 f 1 2 ; ; : : : ; K g , and

$$\rho _ { i, j } = \frac { \text{Cov} ( X _ { i }, X _ { j } ) } { \sigma _ { X _ { i } } \sigma _ { X _ { j } } }, \quad \text{(1.6)} \quad \text{an} \quad \text{flu} \quad \text{ }$$

where Cov denotes the covariance and GLYPH&lt;27&gt; denotes the standard deviation. These correlation coefficients are stored in a matrix, R 2 [ GLYPH&lt;0&gt; 1 1] ; K GLYPH&lt;2&gt; K , which we use to create the proper correlation structure among the pairwise joint distributions. To do this, we create a set of independent, identically distributed random vectors, G 1 ; G 2 ; : : : ; G n i.i.d. GLYPH&lt;24&gt; N ( 0 ; R ) , thus, each Gj , j 2 f 1 2 ; ; : : : ; n g is a K -dimensional random vector with a zero mean multivariate normal distribution with covariance, R . From this set of random vectors we can create a sample set, g 1 ; g 2 ; : : : ; g n , of n samples from N ( 0 ; R ) , where g j are realizations from each identically distributed G j . Thus, we have n samples from standard normal distributions marginally, and each pairwise joint density has the desired correlation in the sample set. From this set of samples, we generate uniformly distributed samples from f u j = ( GLYPH&lt;8&gt; ( g 1 j ) ; GLYPH&lt;8&gt; ( g 2 j ) ; : : : ; GLYPH&lt;8&gt; ( g K j )) g n j = 1 , where GLYPH&lt;8&gt; is the cumulative distribution function of a standard normal random variable, and g j = ( g 1 j ; g 2 j ; : : : ; g K j ) . This results in n samples from a vector of uniformly distributed marginal random variables with our desired pairwise correlations preserved. The final step makes use of the inverse cumulative marginal distributions of each of our parameters, F Xi GLYPH&lt;0&gt; 1 ( xi ) . From these inverse cumulative marginal distributions and the samples u j , we compute f x j = ( F X GLYPH&lt;0&gt; 1 1 ( u 1 j ) ; F X GLYPH&lt;0&gt; 1 2 ( u 2 j ) ; : : : ; F XK GLYPH&lt;0&gt; 1 ( u K j )) g n j = 1 , where u j = ( u 1 j ; u 2 j ; : : : ; u K j ) , which is a set of n sample vectors sampled from the correct marginal distributions provided at the outset and preserved pairwise correlation information.

Our sampling methodology is demonstrated here for a twodimensional random vector in Figure 3 to make the preceding discussion more concrete. The top left plot of the figure is the original two-dimensional joint distribution with each marginal shown as well. This information is distilled into the marginal distributions, f X 1 ( x 1 ) and f X 2 ( x 2 ) , as well as the correlation coefficient matrix, R . Fromthiscorrelation information, two correlated standard normal marginal distributions are created and shown in the bottom left plot of the figure. Samples from these distributions are passed individually through the standard normal cumulative distribution, GLYPH&lt;8&gt; , which leads to two correlated uniform random variables shown in the bottom right plot. Samples from these distributions are passed individually through the inverse marginal distributions of each parameter, respectively. This results in samples from correct marginal distributions with preserved pairwise correlation. This can be seen by comparing the top left and top right plots. The top left plot is the true joint distribution and the top right plot is the joint distribution generated via this methodology, which is correct in terms of the marginal distributions and the correlation coefficient between X 1 and X 2 .

## 4. Results and Discussion

Solid state reactions are not only influenced by chemical driving forces, but also by the interfacial and elastic energy contributions. The bulk free energy ( f bulk ) greatly dominates the overall phase stability of a system but strong couplings to elastic fields arising from lattice and elastic mismatch as well as anisotropy in the properties play a dominant role in controlling the overall thermodynamic stability of the system as well as the topology and morphology of the resulting microstructures. Given the influence of bulk thermodynamic properties, the uncertainties in the bulk free energy require precise quantification first. This will be described in section 4.1. Then the quantified uncertainties of the bulk free energy, the kinetic and microelastic parameters are propagated to microstructural characteristics in section 4.2 through the elastochemical phase field model.

## 4.1. Uncertainty Quantification of Gibbs free energy parameters and Phase Diagram

A thorough parameter uncertainty analysis of the CALPHAD model for the Mg (Si Sn 2 x 1 GLYPH&lt;0&gt; x ) pseudo-binary system has been performed through the MCMC technique (explained in section 3.3.1) against the calculated composition-temperature data sampled from the phase diagram proposed by Kozlov et al. [62]. Then, the calculated parameter uncertainties have been propagated to the Gibbs free energy of phases and phase diagram.

GLYPH&lt;3&gt;

Figure 3: Two-Dimensional demonstration of our sampling strategy that preserves marginal distributions and pairwise correlations.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00638v1-figure-58575312-fe09-40a7-bbc9-805cd7ef6406.jpg)

As shown in Table 1, three parameters 0 a , 0 b and 1 a are selected for each phase in the CALPHAD model, i.e. Mg X Sn,Si 2 f g and liquid (six parameters in total). As mentioned earlier, the deterministically optimized parameters obtained from ThermoCalc have been utilized as initial parameter values for MCMC sampling process; however, lack of knowledge about the parameter probability distributions resulted in the consideration of a uniform (non-informative) prior distributions for the parameters. In addition, GLYPH&lt;6&gt; 50% of the parameter initial values have been considered as the parameter ranges during this process.

During MCMC calibration, 100,000 samples have been generated to ensure parameter convergence. In this regard, plotting the joint frequency distribution of each pair of the model parameters can graphically show the parameter convergences in the parameter space. For example, one of these plots has been shown in Figure 4. As observed in this figure, the red region with the highest density of parameter samples indicates the convergence region in the pair parameter space. Moreover, these kinds of plots can qualitatively show a linear correlation between the model parameters based on the linearity and direction of the convergence region.

Thedegreeofcorrelation between each parameter pair can also be quantified through the Pearson linear coefficient [79] denoted by GLYPH&lt;26&gt; i ; j in Equation 16. It is worth noting that the linear coefficient is a quantity between -1 and 1. Generally, parameters are uncorrelated if GLYPH&lt;26&gt; is close to 0, but highly correlated if GLYPH&lt;26&gt; is close to either -1 or 1. In addition, the negative and positive signs indicate the correlation direction. In Figure 4, a semi-circular shape

Figure 4: Joint frequency distribution between a selected pair of parameters.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00638v1-figure-d298b2e8-eddb-4dc5-be91-07f7463ed587.jpg)

of convergence region and a very close value of GLYPH&lt;26&gt; to zero imply negligible correlation between 0 ass and 0 bliq .

After discarding the burn-in period in the beginning of the MCMC sampling, the marginal probability distribution of each model parameter can be plotted as shown with blue color in Table 1 and Figure 5. As observed in these distributions, all parameters show distributions close to truncated normal. Therefore, a truncated normal probability distribution has been fitted to each marginal probability distribution, where the hyperparameters ( GLYPH&lt;22&gt; and GLYPH&lt;27&gt; ) of these truncated distributions have been

estimated through a maximum likelihood method. The estimated GLYPH&lt;22&gt; and GLYPH&lt;27&gt; of each parameter probability distribution has also been listed in Table 1, which represent the plausible optimal values and uncertainties of the model parameters.

Theparameteruncertainties in Table 1 have been propagated to the molar Gibbs free energy of mixing for different phases in the system to the phase diagram. As mentioned in section 3.3.2, the last 5,000 MCMC samples as an ensemble of the whole convergence region have been used in a forward model analysis scheme to identify the variation of the lines in the phase diagram. Then, 95%BCIshavebeendeterminedbydiscarding2.5%oftheresulting samples from the above and below the variation band at any specified Si composition in the Gibbs energy of mixing curves or the phase diagram. In this regard, Figure 5 demonstrates uncertainty propagation across different levels of CALPHAD modeling. Here, just an example of molar Gibbs energy of mixing for solid phases (Mg X Sn,Si ) at 700 2 f g o K has been plotted to show howuncertainties propagate from the Gibbs free energy curves to the phase diagram. This analysis has been repeated for different temperatures in the range from 200 to 1400 o K to construct the whole phase diagram and its uncertainties. Based on Figure 5, it should be noted that there are very high uncertainties in the curves for the Gibbs energy of mixing as well as the resulting phase diagram in most temperatures. These high uncertainties imply very high impacts of the Gibbs free energy parameter variations on the phase equilibria in this system.

## 4.2. Effect of Parameter Uncertainties on Microstructure

A complete phase-field model for the Mg (Si Sn 2 x 1 GLYPH&lt;0&gt; x ) system requires the definition of a large number of parameters-see Table 1. The proper characterization of all relevant parameters has been a major challenge in mesoscale microstructural models, and it is a common practice to define the values of these parameters deterministically and ignore actual variabilities. In turn, this might ignore the evidence variability on the model outputs. This implies that common sources of uncertainty may have little or no impact on the selection of the value of model parameters. The uncertainties in parameters of the Gibbs energy, obtained through the MCMC approach, are combined with the prior distribution of the parameters of other sub-models (microelasticity and kinetic) to study the microstructural evolution using the phase-field model. The prior information of the all parameters are reported in Table 1. Here we elaborate on the source of these prior distributions.

The process of phase separation is influenced by the elastic anisotropy, and the hardness enhancement observed upon the age hardening relies on a shear modulus difference between the formed domains as well as their coherency strain [80]. Thus, it is of primary interest to take into account the uncertainities in the elastic properties of the Mg 2 (Si x Sn 1 GLYPH&lt;0&gt; x ) system to study the contributions in microstructural uncertainty.

The ab-initio -based elastic constants (C 11 , C 12 , and C 44 ) for (c)-Mg Sn and (c)-Mg Si are reported in the literature-we note 2 2 that these quantities have yet to be measured experimentally. These include both 0 K, and high temperature data which are GLYPH&lt;14&gt; calculated by taking the entropic and/or harmonic effects in the structure. The empirical cumulative density functions (ecdf) of C 11 , C 12 , and C 44 for both phases are shown in the left side of the

Figure 5: An illustration of UP from thermodynamic parameters to the curves associated with the molar Gibbs energy of mixing to the phase diagram for Mg X Sn,Si 2 f g quasi-binary system at 700 o K.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00638v1-figure-2d12232b-6e98-4011-a478-b08f0b2ba6ba.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00638v1-figure-17b5e3b2-1d5e-41a6-a4c2-95017dc50b35.jpg)

microstructure palette in Fig. 6(a). Both ab-initio and experimental lattice parameters for the cubic Mg (Si Sn 2 x 1 GLYPH&lt;0&gt; x ) system at room/high temperature for different ranges of compositions are considered. These values are either provided for individual phases (Mg Sn and Mg Si), or for the parent phase as a function 2 2 of composition. The ecdf of the lattice parameter data are shown

in Fig. 6(b). Using these data, the range of SFTS ( " 0 i j ) for Mg Sn 2 and Mg Si is estimated. This range is used to draw samples from 2 a uniform distribution for propagation of uncertainty in the microelasticity model.

The molar volumes ( Vm ) of the phases are also obtained based on both calculated and measured data. The calculated data is obtained from V cryst m = Vcell : Nv = Z where Vcell is the volume of the unit cell, Nv is the Avogadro's number, and Z is the number of atoms in the cell. Z=3 for both phases considering the fact that the crystal structure of both Mg Sn and Mg Si is Fm 2 2 ¯ 3 m. The distribution is shown in Fig. 6(c), which is very similar to the distribution of lattice parameters shown in this figure. Lower bound of this distribution belongs to the data obtained by division of the molar mass and density distributions.

Experimental diffusion data are used to estimate the atomic mobility. The recent combinatorial diffusion couple study by Viv´ es et al . [70], and the indirect calculation of inter-diffusion coefficient using the forward-simulation method provides an insight about the diffusion at 600 GLYPH&lt;14&gt; C. Their calculations suggest an estimated diffusion coefficient of GLYPH&lt;25&gt; 10 GLYPH&lt;0&gt; 15 ms GLYPH&lt;0&gt; 2 . There is no other information regarding diffusivity measurements to the best of our knowledge. Accordingly, the mobility is estimated by the Einstein's rule ( M = D RT ), or ( M = D ( @ 2 f c 2 ) ) [39]. Hence, a uni-

@ form sampling is carried out around this value. Initial composition of the alloy is taken to be uniformly distributed between the xSi = 0 3 : GLYPH&lt;0&gt; 0 5 : to make sure we are under the chemical miscibility region.

In order to propagate the uncertainties of the prior data, it is necessary to carry out high-throughput phase field simulations of microstructure evolution in the Mg (Si Sn 2 x 1 GLYPH&lt;0&gt; x ) system. Using the strategy described in section 3.3.3, we sampled 10,000 combinations of the parameters out of the prior distributions, which have been fed to the phase-field solver (Fortran code). The whole process of 1) data initiation, 2) environment preparation and 3) Fortran code run in the #Terra super-computing cluster at Texas A&amp;M University is automated through a Python wrapper. A square simulation cell with 512 GLYPH&lt;2&gt; 512 grid numbers, where Lx = Ly = 350 nm is used to perform the simulations. The snapshots of the obtained microstructures for 800 samples are shown in the form of a palette in Fig. 6. Sixteen example microstructures (same evolution time) are also extracted from the palette for better effective demonstration and is located in the bottom section of this figure.

Phase-field models tend to be highly nonlinear so that its output can differ qualitatively, depending on the region in the input/parameter space where the sample is taken. Phase-field models are also highly complex in its formulation and are thus not amenable to intrusive approaches of UP. Moreover, they are computationally expensive, with full three-dimensional realizations of the simulations requiring upwards of 10,000 CPU-hours in some of the fastest supercomputers available. Finally, the input/parameter space is high-dimensional, with more than &gt; 20 stochastic input conditions (i.e. temperature-time) and model parameters that dramatically affect the thermodynamic and kinetic state of the system, as shown in Figure 6.

It must be noted that since we are interested in the elastochemical interactions in the microstructure, the effect of variations in elastic constants on the microstructure were considered. The lo- cal long-range interactions (i.e. strain and/or stress fields) of the heterogeneous multi-phase nanostructure are quite sensitive to the selection of the elastic parameters, and must be considered with a great care. Figure 7 demonstrates the elastic strain maps ( " el i j ) for two distinct microstructures with different set of elastic constants (i.e. shear constants), and very similar eigenstrain ( " 0 ii ) values. Both cases can be categorized as cuboid-type microstructures. Though in the first case, (Fig. 7a,b, and c) the morphology is perfectly cuboid with a deviation from unimodal particle-size distribution. In the second case (Fig. 7d,e, and f), the particles are sheared at the corners to evolve towards a cuboidal shape during the coarsening stages as well. While microelastic effects can be considered to be very similar, slight differences in other parameters of the chain model yielded different observations.

## 4.3. Microstructure Quantification

Uncertainty propagation consists of linking the uncertainty in the input conditions and model parameters to the generated uncertainty in the model output. To properly quantify the uncertainty and study the way it propagates across the simulations, it is necessary to define quantity(ies) of interest (QoI) whose distribution can then be interpreted as a measure of variance resulting from the stochastic nature of the input space.

Atraditional way to quantify a microstructure is to exploit several QoI in the (micro)structure. A challenging aspect of the uncertainty propagation effort in this work is the large dimensionality of both the input and the output space. The large dimensionality of the input space originates from the many parameters that are necessary to complete the phase field model. Given the fact that the microstructure space is highly heterogeneous, it is to be expected that many QoIs would be necessary to complete characterize the effect of model inputs on the morphology and topology of the resulting microstructures. In order to handle these very complex spaces, we will resort to machine learning frameworks, as will be described below.

The propagation of uncertainty first requires the quantification of changes in the microstructure. Given the large dimensionality of the output (512x512 pixels), dimensional reduction is necessary. Moreover, in order to enhance the interpretability of the analysis it is desirable for the reduced dimensions to have physical significance. The determination of QoIs typically involves extraction of features or disclosing a bank of descriptors that can be used to train a classifier based on the frequency of observations [81, 82]. A conventional yet very useful approach is to use semantic texton forests [83]. This is specifically useful when the phasefield variable is composition and can also be broadly used in the real-image data. In addition, visual words [84] can be used as powerful discrete image representations for categorization. Another way to tackle this is to utilize one or a combination of Filterbank responses (e.g. Fourier or other sort of wavelets), and one or combination of invariant descriptors (e.g. SIFT [85, 86]).

A series of microstructure analytic tools were used to determine QoIs to evaluate the developed UP framework. We use metrics such as average feature size, area fraction, composition of the phases, aspect ratio as well as increasingly common approaches such as n-point statistics [82], frequency-domain analysis (c.f. Figure 8). In the later case, a general information about the morphology and orientation of the particles are reflected in the frequency space. Many of the commonly used metrics provide in-

Microstructure palette (800 out of 10,000 for a fixed time)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00638v1-figure-9d576f06-f746-4deb-a7d4-cae584bc3d35.jpg)

Figure 6: Center) Microstructure library from 800 phase-field runs (out of 10,000 run) for the Mg 2 X Si,Sn f g system at a fixed time. 10,000 points in a 18-dimensional parameter space were sampled from the prior distributions shown in Table 1. ecdf plots for a) elastic constants, b) lattice parameters, c) molar volumes of the two product phases. d) Example microstructures extracted from the palette.

Figure 7: Effect of the local strain due to inhomogenous elastic effects in the microstructure. a and d) " el 11 , b and e) " el 22 , c and f) " el 12 . All input parameters are different in these two cases. The SFTS value is " 0 ii = GLYPH&lt;0&gt; 0 0176 : , elastic constants are GLYPH&lt;2&gt; Mg Si 2 : C 11 = 120 2 : ; C 12 = 22 1 : ; C 44 = 56 3 : GPa , GLYPH&lt;3&gt; GLYPH&lt;2&gt; Mg Sn 2 : C 11 = 76 9 : ; C 12 = 19 2 : ; C 44 = 38 4 : GPa GLYPH&lt;3&gt; for the case in the first row, and " 0 ii = GLYPH&lt;0&gt; 0 0191 : , GLYPH&lt;2&gt; Mg Si 2 : C 11 = 120 5 : ; C 12 = 22 2 : ; C 44 = 48 8 : GPa , GLYPH&lt;3&gt; GLYPH&lt;2&gt; Mg Sn 2 : C 11 = 79 14 : ; C 12 = 38 2 : ; C 44 = 21 5 : GPa GLYPH&lt;3&gt; for the case in the second row.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00638v1-figure-bdb4c8de-7a55-4eef-b780-548f6562ef46.jpg)

formation about the average state of microstructure spaces, but in many cases they do not provide information about their topology.

In the case of transport behavior (such as phonon conductivity), topology may play an important role and metrics quantifying this microstructure feature are necessary.

Figure 8: a,c) Example microstructures, and the corresponding b,d) Power spectrum maps obtained by Fourier transformation of the microstructures generated from the elastochemical phase field model.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00638v1-figure-c51365b2-dff5-4724-a6d5-dfb0d06caa3f.jpg)

Table 2 summarizes the list of QoIs and the associated statistics for phase-field model outputs. It contains the posterior, means, index of dispersion denoted by variance to mean ratio (VMR) of the obtained QoIs. VMR is an index to quantify the dispersion of a probability and measures the clusterability and variability in the data. For VMR &gt; 1 , and 0 &lt; VMR &lt; 1 the data are over-dispersed and under-dispersed, respectively. VMR = 0 is associated to a random data-set and VMR = 1 corresponds to a Poisson distribution. Except the last three, all other QoIs demonstrate multimodal distributions. We wish to note that this multi-modality would render simplified uncertainty analysis frameworks-such as those based on sampling of min, max and mean values of input

parameters-ineffective.

It should be noted that C GLYPH&lt;11&gt; and and C GLYPH&lt;12&gt; are considered to be the equilibrum compositions of Mg Sn and Mg Si phases, respec2 2 tively. Both of these demonstrate a trimodal distribution: I) the first peak belong to the case where the process of phase separation has entirely finished and the elastic interactions are weak. In this case, we obtain the equilibrium values corresponding to the sampled phase-diagram. II) In the second case, the composition between elastic and chemical driving forces are great, and morphology of the particles are similar to those reported in elastochemical interactions. III) In the third case, either kinetics is slow or the system tends to dissolve rather than phase separate. In such cases, the elastic driving forces are significantly stronger than the chemical ones. This is believed to be achieved by the non-equilibrium synthesis conditions (refer to [39] for further information).

Table 2: List of extracted QoIs and their posterior distributions, mean, standard deviation, minimum and maximum from the entire phase-field runs.

| Target variables   | Posterior Distribution   | GLYPH<22>                   | Dispersion ( GLYPH<27> 2 )   |
|--------------------|--------------------------|-----------------------------|------------------------------|
| C GLYPH<11>        |                          | 0.57                        | 0.056                        |
| C GLYPH<12>        |                          | 0.34                        | 0.026                        |
| Char. Length       |                          | 8.62 GLYPH<2> 10 GLYPH<0> 8 | 1.08 GLYPH<2> 10 GLYPH<0> 7  |
| Area fraction      |                          | 0.17                        | 0.30                         |
| Roundness          |                          | 0.99                        | 0.70                         |
| Diagonal = width   |                          | 1.18                        | 0.29                         |
| Diagonal = height  |                          | 0.98                        | 0.05                         |
| ¯ GLYPH<22> chem   |                          | -156.06                     | -1012.1                      |
| ¯ GLYPH<22> elas   |                          | -27.83                      | -214.75                      |
| ¯ GLYPH<22> int    |                          | 1.59 GLYPH<2> 10 GLYPH<0> 9 | -6.01 GLYPH<2> 10 GLYPH<0> 6 |

The characteristic length scale for each microstructure is obtained by calculating the radial basis power spectrum of the Fourier transform of the microstructure (2D composition space). The distribution is bimodal and its peaks belong to either the not-decomposed, or the decomposed structures with certain morphologies. The calculated area fractions uniformly ascends from zero fraction for dissolved or not-decomposed microstructures to the decomposed microstructures. The probability density of observation of not-decomposed or dissolved microstructures are higher than its counterpart for the decomposed ones. Roundness is extracted based on the shape of the morphology that appears in the power spectrum domain. Figure 8 demonstrates examples of different power spectrums that are obtained for the shown microstructures. For circular morphology the power spectrum is symmetric and appears to be a circle with a diameter based on the coarsening stage. The power spectrum appears to be in the form of rectangle/square for a perfectly oriented morphology.

Quantities Diagonal = width and Diagonal = height are the relative ratio of the diagonals of the power spectrum centroid with respect to its width and height. The closer the value is to p 2 : height 2 and/or p 2 : width 2 , the morphologies tend to shear more. Moreover, the higher these values get, the anisotropy becomes more evident. On the other hand, the mean value of each of the contributing energy fields are also obtained for the current sampling. The mean of bulk and elastic driving forces in the domain demonstrate a unimodal probability distributions. However, the mean of interfacial driving forces demonstrates multimodal isolated peaks.

Table 2 represents the total variance in the microstructure space assuming the CDF of the input parameters and their statistical correlation. The propagated uncertainty is represented through different QoIs. However, these results are aggregated and further analysis is necessary in order to examine the effect of individual parameters (or sets of parameters) on the resulting microstructures. The large dimensionality of the input and output spaces, as we discussed above, makes it necessary to rely on machine learning approaches that facilitate the analysis of the microstructure space and their relation to model inputs.

## 4.4. Application of the Materials Informatic Techniques in Microstructural Evaluation

Aprimary purpose of data-mining techniques is to help in determination of the possible patterns for better prediction abilities, which form the foundations for understanding materials behavior [87]. Materials informatics is an interdisciplinary blend of statistics, machine learning, artificial intelligence, pattern recognition and materials science. Here we employ a few core tasks (e.g. Cluster analysis, Anomaly detection, etc.) on the posterior data generated in this study. This collective integration of statistical learning tools with experimental and/or computational materials science allows for an informatics driven strategy for materials design and development under the framework of ICME.

A preliminary simple hierarchical clustering using different distance metrics are performed to elucidate the structure of the available microstructural data. The results are demonstrated as a corresponding dendrogram in Fig. 9. This calculation is based on the the 18 material parameter as the inputs and the eight QoI as the target variables. In addition, several distance metrics (i.e. Euclidean, Manhattan, Mahalanobis, Spearman and Pearson) are

Figure 9: Hierarchical clustering using Mahalanobis distance matrix. The horizontal line indicates 78% height ratio. The images are clustered into two main subcategory of 'decomposed' versus 'not-decomposed' microstructures. Each of these main categories are also divided into three sub-levels with an acceptable error.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00638v1-figure-89e89d61-8f5f-4d91-87f9-752dffc1f332.jpg)

used. The result in Fig. 9 is based on Mahalanobis distance metric where the elements are taken in a pairwise fashion between the elements of a given set using dM ( x y ; ) = p ( x GLYPH&lt;0&gt; y ) T S GLYPH&lt;0&gt; 1 ( x GLYPH&lt;0&gt; y ) where S is the Covariance matrix. The dendrogram shown in Fig. 9 is created by the Ward's linkage criterion (increase in variance for the cluster being merged) and aims to indicate the similarity/dissimilarity among annotation categories. The five sub-clusters of the dendrogram shown in different colors are annotated with selected representing microstructures. As seen, this dendrogram clusters the microstructures into two main categories, i.e. 'decomposed' and 'not-decomposed' classes, with an acceptable error.

A comparison of several classifiers is performed on the extracted QoIs data-set. Classification methods are suited for cases where the class label is discrete. Hence, we simply use the 'notdecomposed' versus 'decomposed' class labels. Here we return to the connection with performance in TE materials by noting that microstructures that correspond to non decomposed states could be associated with alloying/mass phonon scattering, while decomposed microstructures corresponds to interfacial phonon scattering. The length scales of different scattering mechanisms are different and are thus expected to change the phonon transport characteristics and the corresponding thermoelectric performance of the Mg 2 f Sn,Si g system [56].

Figure 10 illustrates the nature of such decision boundaries in the 2D input parameter regions which can be used to determine the desired regions (alloy/mass vs interface scattering) in the material parameter space. Several classifiers are tested on all of the parameter pairs and the selected classifiers to report are the Nearest neighbors, Gaussian process, Radial Basis Function (RBF) kernel Support Vector Machine (SVM), Neural Network. Each subfigure in Fig. 10 illustrates the nature of decision boundaries for a given pair of the input parameter using different classifiers. The training points in these plots are shown in solid colors and the testing points are defined by the semi-transparent contour. We first randomly extracted 10% of the original data obtained from a set of microstructures frozen at a fixed time. 60% of this data is used as the training set and 40% as the test set. Figure 10a demon- strates the Nearest neighbors classifier in the (alloy composition ( XSi ), a 0 SS ) space with a transitioning boundary where the points are mixed in the boundary of the two class. Figure 10b demonstrates the Gaussian process classifier [88] result in the ( " T , a 0 SS ) space with a smooth, and continuous separation boundary between the decomposed (blue) and not-decomposed (red) clusters. Figure 10c shows the classification result in (C 11 :Mg Si, 2 a 0 SS ) space using RBF SVM classifier that reveals the transitioning bounderies between the two clusters. The red cluster preserves a larger area while certain points of this cluster are mixed with the blue cluster. The result associated with the trained neural network classifier in ( GLYPH&lt;20&gt; , a 0 SS ) space is shown in Fig. 10d, which again illustrates a smooth, and continuous separation boundary between the two classes. However, it should be noted that no clear classification boundaries are obtained in most pair-parameter spaces which canbeattributed to the lack of recognition through 2D projections of the high dimensional parameter space. We note that in these cases, the application of a dimensional reduction technique, e.g. principal component analysis (PCA), might be helpful for more rigorous classification study.

Whilefurtherinvestigation of the microstructure space through machine learning approaches may be warranted in order to fully characterize the connections between input parameters and resulting microstructures, this is well beyond the scope of the present work. It is in fact our intent to further explore the use of different machine learning frameworks into the obtained microstructure dataset.

## 5. Summary and Conclusion

As with any computational analysis, uncertainty quantification/propagation plays a major role in predicting the outcome of multi-scale models in materials science. One of the most important tasks in the materials design under ICME framework is propagation of uncertainties of parameters across the multi-scale models that connect process-structure-propertyperformance. Hence, we addressed an existing challenge in materials science, i.e., UP from the thermodynamic parameters to the microstructural features through a chain of CALPHAD, microelasticity, and phase field models. This is a multi-step uncertainty propagation with the techniques used in each step motivated by the type of problem that is being dealt with. Technically, thermodynamic parameters in the CALPHAD model and their underlying imposed uncertainties obtained from an MCMC sampling approach has been propagated to Gibbs free energy of phases and equilibrium phase diagram through a forward analysis of an ensemble of these samples. Then, these uncertainties as well as the uncertainties of microelastic and kinetic parameters have been propagated to the microstructural features using a Gaussian copula sampling approach.

While the present study has multiple aims, one of these is to rectify the common belief towards the deterministic assumption about the parameter values in the phase field models as well as the experimental data used for the calibration process.

The propagation of uncertainty in the prior parameter space using model chains resulted in a massive microstructure dataset ( GLYPH&lt;24&gt; 50TB). The quantities of interest in the microstructures are identified and extracted using automated frameworks from the large amounts of data that is generated by HT-phase-field

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00638v1-figure-78ae8f12-bee2-47dd-8434-fe3148dc2a7b.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00638v1-figure-bb638f60-7317-4b1f-8072-cd6907409409.jpg)

Figure 10: Classification results for a) Nearest neighbors classifier in (alloy composition ( XSi ), 0 aSS ) space, b) Gaussian process classifier in ( " T , 0 aSS ) space, c) Radial basis function support vector machine (RBF SVM) classifier in ( C 11 :Mg Si, 2 0 aSS ) space d) Neural network classifier in ( GLYPH&lt;20&gt; , 0 aSS ) space. For interpretation of the colors the reader is referred to the online version. stands for not-decomposed or dissolved microstructures, while stands for decomposed microstructures with different morphologies.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00638v1-figure-9177469f-a46e-4f0f-9697-137990406403.jpg)

runs. We used eight QoIs to map the obtained probability distributions of parameters into probability distributions of the extracted quantities. Then, data-mining techniques are employed to find patterns in the parameter space that can contribute to better understanding of process-microstructure relations. The results show that the data is clustered into two main categories of 'decomposed' and 'not-decomposed' microstructures. These results can be very useful for engineering material behaviour in favor of specific phonon scattering mechanism and/or better thermoelectric response. The proposed framework is generalizable to applications to other materials problems and microstructuresensitive properties. As of this writing further exploration of the input parameter-microstructure space is in progress and will be reported in our future work.

Developing general QoIs of the microstructural space will enable better quantification of the uncertainty propagated through the models establishing inverse maps to connect regions in the microstructure space to the corresponding regions in the input space. While standard microstructure analysis approaches can be used to represent microstructure spaces in a compact way. To establish PSP connections, it is very important to predict properties that depend on the specific features of such microstructures.

## Acknowledgment

The authors would like to acknowledge the #Terra supercomputing facility of the Texas A&amp;M University, for providing computing resources useful in conducting the research reported in this paper. The data is currently hosted by this facility. This research was supported by the National Science Foundation under NSF Grant No. CMMI-1462255 and NSF-CMMI-1663130. RA and DA also acknowledge the support of ARL through [grant No. W911NF-132-0018].

## References

## References

- [1] Y. Qian, C. Jackson, F. Giorgi, B. Booth, Q. Duan, C. Forest, D. Higdon, Z. J. Hou, G. Huerta, Uncertainty quantifi-
- cation in climate modeling and projection, Bulletin of the American Meteorological Society 97 (5) (2016) 821-824.
- [2] S. Gupta, D. Ghosh, Uncertainty quantification in structural engineering: Current status and computational challenges, in: Uncertainty Quantification in Computational Science: TheoryandApplicationinFluidsandStructuralMechanics, World Scientific, 2017, pp. 119-149.
- [3] S. Ghoreishi, D. Allaire, Adaptive uncertainty propagation for coupled multidisciplinary systems, AIAA Journal (2017) 3940-3950.
- [4] E. Begoli, T. Bhattacharya, D. Kusnezov, The need for uncertainty quantification in machine-assisted medical decision making, Nature Machine Intelligence 1 (1) (2019) 20.
- [5] J. H. Panchal, S. R. Kalidindi, D. L. McDowell, Key computational modeling issues in integrated computational materials engineering, Computer-Aided Design 45 (1) (2013) 4-25.
- [6] D. L. McDowell, S. R. Kalidindi, The materials innovation ecosystem: a key enabler for the materials genome initiative, Mrs Bulletin 41 (4) (2016) 326-337.
- [7] R. Arroyave, S. Shields, C.-N. Chang, D. Fowler, R. Malak, D. Allaire, Interdisciplinary research on designing engineering material systems: Results from a national science foundation workshop, Journal of Mechanical Design.
- [8] R. Arr´ oyave, D. L. McDowell, Systems approaches to materials design: Past, present, and future, Annual Review of Materials Research 49.
- [9] J. Allison, Integrated computational materials engineering: A perspective on progress and future steps, Jom 63 (4) (2011) 15-18.
- [10] X. Liu, D. Furrer, J. Kosters, J. Holmes, Vision 2040: A roadmapforintegrated, multiscale modeling and simulation of materials and systems (2018).
- [11] E. Weinan, Principles of multiscale modeling, Cambridge University Press, 2011.

- [12] M. Sanghvi, P. Honarmandi, V. Attari, T. Duong, R. Arroyave, D. L. Allaire, Uncertainty propagation via probability measure optimized importance weights with application to parametric materials models, in: AIAA Scitech 2019 Forum, 2019, p. 0967.
- [13] W.K.Liu, L. Siad, R. Tian, S. Lee, D. Lee, X. Yin, W. Chen, S. Chan, G. B. Olson, L.-E. Lindgen, et al., Complexity science of multiscale materials via stochastic computations, International Journal for Numerical Methods in Engineering 80 (6-7) (2009) 932-978.
- [14] B. Kouchmeshky, N. Zabaras, Microstructure model reduction and uncertainty quantification in multiscale deformation processes, Computational Materials Science 48 (2) (2010) 213-227.
- [15] M. Koslowski, A. Strachan, Uncertainty propagation in a multiscale model of nanocrystalline plasticity, Reliability Engineering &amp; System Safety 96 (9) (2011) 1161-1170.
- [16] S. Salehghaffari, M. Rais-Rohani, E. Marin, D. Bammann, A new approach for determination of material constants of internal state variable based plasticity models and their uncertainty quantification, Computational Materials Science 55 (2012) 237-244.
- [17] P. Honarmandi, A. Solomou, R. Arroyave, D. C. Lagoudas, Uncertainty quantification of the parameters and predictions of a phenomenological constitutive model for thermally induced phase transformation in ni-ti shape memory alloys, Modelling and Simulation in Materials Science and Engineering.
- [18] N. Saunders, A. P. Miodownik, CALPHAD (calculation of phase diagrams): a comprehensive guide, Vol. 1, Elsevier, 1998.
- [19] I. Steinbach, B. B¨ttger, J. Eiken, N. Warnken, S. Fries, Calo phad and phase-field modeling: a successful liaison, Journal of phase equilibria and diffusion 28 (1) (2007) 101-106.
- [20] M. Stan, B. J. Reardon, A Bayesian approach to evaluating the uncertainty of thermodynamic data and phase diagrams, Calphad 27 (3) (2003) 319-323.
- [21] P. Honarmandi, T. C. Duong, S. F. Ghoreishi, D. Allaire, R. Arroyave, Bayesian uncertainty quantification and information fusion in calphad-based thermodynamic modeling, Acta Materialia 164 (2019) 636-647.
- [22] W. D. Thomison, D. L. Allaire, A model reification approach to fusing information from multifidelity information sources, in: 19th AIAA non-deterministic approaches conference, 2017, p. 1949.
- [23] P. Honarmandi, N. Paulson, R. Arroyave, M. Stan, Uncertainty quantification and propagation in calphad modelling, Modelling and Simulation in Materials Science and Engineering.
- [24] R. A. Otis, Z.-K. Liu, High-throughput thermodynamic modeling and uncertainty quantification for icme, JOM 69 (5) (2017) 886-892.
- [25] B. Bocklund, R. Otis, A. Egorov, A. Obaied, I. Roslyakova, Z.-K. Liu, Espei for efficient thermodynamic database development, modification, and uncertainty quantification: application to cu-mg, arXiv preprint arXiv:1902.01269.
- [26] S. G. Kim, W. T. Kim, T. Suzuki, Phase-field model for binary alloys, Physical review e 60 (6) (1999) 7186.
- [27] I. Steinbach, L. Zhang, M. Plapp, Phase-field model with finite interface dissipation, Acta Materialia 60 (6-7) (2012) 2689-2701.
- [28] V. Attari, S. Ghosh, T. Duong, R. Arroyave, On the interfacial phase growth and vacancy evolution during accelerated electromigration in cu/sn/cu microjoints, Acta Materialia 160 (2018) 185-198.
- [29] V. Attari, A. Cruzado, R. Arroyave, Exploration of the microstructure space in tialzrn ultra-hard nanostructured coatings, Acta Materialia 174 (2019) 459 - 476.
- [30] K. Karayagiz, L. Johnson, R. Seede, V. Attari, B. Zhang, X. Huang, S. Ghosh, T. Duong, I. Karaman, A. Elwany, et al., Finite interface dissipation phase field modeling of ninb under additive manufacturing conditions, arXiv preprint arXiv:1906.10200.
- [31] T. Duong, R. E. Hackenberg, V. Attari, A. Landa, P. E. Turchi, R. Arroyave, Probing discontinuous precipitation in u-nb, arXiv preprint arXiv:1907.00918.
- [32] N. Wang, S. Rokkam, T. Hochrainer, M. Pernice, A. ElAzab, Asymptotic and uncertainty analyses of a phase field model for void formation under irradiation, Computational materials science 89 (2014) 165-175.
- [33] L. S. Leon, R. C. Smith, W. S. Oates, P. Miles, Identifiability and active subspace analysis for a polydomain ferroelectric phase field model, in: ASME 2017 Conference on Smart Materials, Adaptive Structures and Intelligent Systems, American Society of Mechanical Engineers, 2017, pp. V002T03A022-V002T03A022.
- [34] B. B¨ ottger, R. Altenfeld, G. Laschet, G. Schmitz, B. St¨hr, o B. Burbaum, An icme process chain for diffusion brazing of alloy 247, Integrating Materials and Manufacturing Innovation 7 (2) (2018) 70-85.
- [35] S. H. Lee, W. Chen, A comparative study of uncertainty propagation methods for black-box-type problems, Structural and Multidisciplinary Optimization 37 (3) (2009) 239.
- [36] P. Honarmandi, A. Solomou, R. Arroyave, D. Lagoudas, Parametric analysis of a phenomenological constitutive model for thermally induced phase transformation in ni-ti shape memory alloys, arXiv preprint arXiv:1808.07377.
- [37] G.Tapia, L. Johnson, B. Franco, K. Karayagiz, J. Ma, R. Arroyave, I. Karaman, A. Elwany, Bayesian calibration and uncertainty quantification for a physics-based precipitation model of nickel-titanium shape-memory alloys, Journal of Manufacturing Science and Engineering 139 (7) (2017) 071002.

- [38] S. Wang, N. Mingo, Improved thermoelectric properties of mg si ge sn 2 x y 1 GLYPH&lt;0&gt; x GLYPH&lt;0&gt; y nanoparticle-in-alloy materials, Applied Physics Letters 94 (20) (2009) 203109.
- [39] S. Yi, V. Attari, M. Jeong, J. Jian, S. Xue, H. Wang, R. Arroyave, C. Yu, Strain-induced suppression of the miscibility gap in nanostructured mg si-mg sn solid solutions, Journal 2 2 of Materials Chemistry A 6 (36) (2018) 17559-17570.
- [40] T. Aizawa, R. Song, Mechanically induced reaction for solid-state synthesis of mg si and mg sn, 2 2 Intermetallics 14 (4) (2006) 382-391.
- [41] Y. Noda, H. Kon, Y. Furukawa, N. Otsuka, I. A. Nishida, K. Masumoto, Preparation and thermoelectric properties of mg si 2 1 GLYPH&lt;0&gt; x ge x (x= 0.0 GLYPH&lt;24&gt; 0.4) solid solution semiconductors, Materials Transactions, JIM 33 (9) (1992) 845-850.
- [42] N. A. Heinz, T. Ikeda, Y. Pei, G. J. Snyder, Applying quantitative microstructure control in advanced functional composites, Advanced Functional Materials 24 (15) (2014) 2135-2153.
- [43] H. Balout, P. Boulet, M.-C. Record, Thermoelectric properties of sn-containing mg2si nanostructures, The Journal of Physical Chemistry C 119 (31) (2015) 17515-17521.
- [44] S. Gorsse, S. Viv` es, P. Bellanger, D. Riou, L. Laversenne, S. Miraglia, D. R. Clarke, Multi-scale architectured thermoelectric materials in the mg2 (si, sn) system, Materials Letters 166 (2016) 140-144.
- [45] S. Gorsse, P. Bauer Pereira, R. Decourt, E. Sellier, Microstructure engineering design for thermoelectric materials: an approach to minimize thermal diffusivity, Chemistry of Materials 22 (3) (2009) 988-993.
- [46] W. Wang, L. Yang, N. Wang, H. Zhang, Y. Jia, Design of a thermoelectric material using the calphad technique: Thermodynamic reassessment of the al-sb-zn system, Journal of Electronic Materials 47 (1) (2018) 261-271.
- [47] D. U. Furrer, Application of phase-field modeling to industrial materials and manufacturing processes, Current Opinion in Solid State and Materials Science 15 (3) (2011) 134140.
- [48] G. J. Schmitz, U. Prahl, Toward a virtual platform for materials processing, JOM 61 (5) (2009) 19-23.
- [49] Y. C. Yabansu, P. Steinmetz, J. H¨ otzer, S. R. Kalidindi, B. Nestler, Extraction of reduced-order process-structure linkages from phase-field simulations, Acta Materialia 124 (2017) 182-194.
- [50] D. Wang, R. Shi, Y. Zheng, R. Banerjee, H. L. Fraser, Y. Wang, Integrated computational materials engineering (icme) approach to design of novel microstructures for tialloys, Jom 66 (7) (2014) 1287-1298.
- [51] G. J. Schmitz, B. B¨ ottger, M. Apel, Microstructure modeling in icme settings, in: Proceedings of the 3rd World Congress on Integrated Computational Materials Engineering (ICME 2015), Springer, 2015, pp. 165-172.
- [52] G. J. Snyder, E. S. Toberer, Complex thermoelectric materials, in: Materials For Sustainable Energy: A Collection of Peer-Reviewed Research and Review Articles from Nature Publishing Group, World Scientific, 2011, pp. 101-110.
- [53] J. Yang, H.-L. Yip, A. K.-Y. Jen, Rational design of advanced thermoelectric materials, Advanced Energy Materials 3 (5) (2013) 549-565.
- [54] G. Tan, L.-D. Zhao, M. G. Kanatzidis, Rationally designing high-performance bulk thermoelectric materials, Chemical reviews 116 (19) (2016) 12123-12149.
- [55] P. Gorai, D. Gao, B. Ortiz, S. Miller, S. A. Barnett, T. Mason, Q. Lv, V. Stevanovi´ c, E. S. Toberer, Te design lab: A virtual laboratory for thermoelectric material design, Computational Materials Science 112 (2016) 368-376.
- [56] A. S. Tazebay, S.-I. Yi, J. K. Lee, H. Kim, J.-H. Bahk, S. L. Kim, S.-D. Park, H. S. Lee, A. Shakouri, C. Yu, Thermal transport driven by extraneous nanoparticles and phase segregation in nanostructured mg2 (si, sn) and estimation of optimum thermoelectric performance, ACS applied materials &amp; interfaces 8 (11) (2016) 7003-7012.
- [57] W. Liu, X. Tan, K. Yin, H. Liu, X. Tang, J. Shi, Q. Zhang, C.Uher, Convergenceofconductionbandsasameansofenhancing thermoelectric performance of n-type mg si 2 1 GLYPH&lt;0&gt; x sn x solid solutions, Physical Review Letters 108 (16) (2012) 166601.
- [58] Q. Zhang, J. He, T. Zhu, S. Zhang, X. Zhao, T. Tritt, High figures of merit and natural nanostructures in mg 2 si 0.4 sn 0.6 based thermoelectric materials, Applied Physics Letters 93 (10) (2008) 102109.
- [59] V. Zaitsev, M. Fedorov, E. Gurieva, I. Eremin, P. Konstantinov, A. Y. Samunin, M. Vedernikov, Highly effective mg si 2 1 GLYPH&lt;0&gt; x sn x thermoelectrics, Physical Review B 74 (4) (2006) 045207.
- [60] G. Nolas, D. Wang, M. Beekman, Transport properties of polycrystalline mg si 2 1 GLYPH&lt;0&gt; y sb y (0 GLYPH&lt;20&gt; y &lt; 0 4) : , Physical Review B 76 (23) (2007) 235204.
- [61] S. Vives, P. Bellanger, S. Gorsse, C. Wei, Q. Zhang, J.-C. Zhao, Combinatorial approach based on interdiffusion experiments for the design of thermoelectrics: Application to the mg2 (si, sn) alloys, Chemistry of Materials 26 (15) (2014) 4334-4337.
- [62] A. Kozlov, J. Gr¨bner, R. Schmid-Fetzer, Phase formation o in mg-sn-si and mg-sn-si-ca alloys, Journal of Alloys and Compounds 509 (7) (2011) 3326-3337.
- [63] R. Viennois, C. Colinet, P. Jund, J.-C. T´ edenac, Phase stability of ternary antifluorite type compounds in the quasibinary systems mg2x-mg2y (x, y= si, ge, sn) via ab-initio calculations, Intermetallics 31 (2012) 145-151.
- [64] E. Nikitin, V. Bazanov, V. Tarasov, Thermoelectric properties of mg2si-mg2sn solid solutions, Sov. Phys. Solid State 3 (1961) 2648-2652.

| [65]   | G. Polymeris, N. Vlachos, A. Khan, E. Hatzikraniotis, C. B. Lioutas, A. Delimitis, E. Pavlidou, K. Paraskevopoulos, T. Kyratsi, Nanostructure and doping stimulated phase sep- aration in high-zt mg2si0.55sn0.4ge0.05compounds, Acta Materialia 83 (2015) 285-293.   |
|--------|-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [66]   | D. Rowe, V. Shukla, N. Savvides, Phonon scattering at grain boundaries in heavily doped fine-grained silicon- germanium alloys, Nature 290 (5809) (1981) 765.                                                                                                         |
| [67]   | D. Medlin, G. Snyder, Interfaces in bulk thermoelectric ma- terials: a review for current opinion in colloid and interface science, Current Opinion in Colloid & Interface Science 14 (4) (2009) 226-235.                                                             |
| [68]   | T. Ikeda, L. A. Collins, V. A. Ravi, F. S. Gascoin, S. M. Haile, G. J. Snyder, Self-assembled nanometer lamellae of thermoelectric pbte and sb2te3 with epitaxy-like interfaces, Chemistry of Materials 19 (4) (2007) 763-767.                                        |
| [69]   | Y. Pei, N. A. Heinz, A. LaLonde, G. J. Snyder, Combina- tion of large nanostructures and complex band structure for high performance thermoelectric lead telluride, Energy & Environmental Science 4 (9) (2011) 3640-3645.                                            |
| [70]   | S. Viv´s, e P. Bellanger, S. Gorsse, C. Wei, Q. Zhang, J.-C. Zhao, Combinatorial approach based on interdiffusion ex- periments for the design of thermoelectrics: Application to the mg2 (si, sn) alloys, Chemistry of Materials 26 (15) (2014) 4334-4337.           |
| [71]   | P. Gorai, V. Stevanovi´, c E. S. Toberer, Computationally guided discovery of thermoelectric materials, Nature Re- views Materials 2 (9) (2017) 17053.                                                                                                                |
| [72]   | L. Q. Chen, J. Shen, Applications of semi-implicit fourier- spectral method to phase field equations, Computer Physics Communications 108 (2-3) (1998) 147-158.                                                                                                       |
| [73]   | A. Gelman, H. S. Stern, J. B. Carlin, D. B. Dunson, A. Ve- htari, D. B. Rubin, Bayesian data analysis, Chapman and Hall/CRC, 2013.                                                                                                                                    |
| [74]   | D. Foreman-Mackey, D. W. Hogg, D. Lang, J. Goodman, emcee: the mcmc hammer, Publications of the Astronomi- cal Society of the Pacific 125 (925) (2013) 306.                                                                                                           |
| [75]   | G. J. McRae, J. W. Tilden, J. H. Seinfeld, Global sensitiv- ity analysis-a computational implementation of the fourier amplitude sensitivity test (fast), Computers&ChemicalEn- gineering 6 (1) (1982) 15-25.                                                         |
| [76]   | D. Huntington, C. Lyrintzis, Improvements to and limita- tions of latin hypercubesampling, Probabilistic engineering mechanics 13 (4) (1998) 245-253.                                                                                                                 |
| [77]   | H. Haario, E. Saksman, J. Tamminen, et al., An adaptive metropolis algorithm, Bernoulli 7 (2) (2001) 223-242.                                                                                                                                                         |
| [78]   | K. Choi, Y. Noh, I. Lee, Reliability-based design optimiza- tion with confidence level for problems with correlated in- put distributions, in: 6th China-Japan-Korea Joint Sympo- sium on Optimization of Structural and Mechanical Sys- tems, 2010.                  |

- [79] D. J. Battle, P. Gerstoft, W. S. Hodgkiss, W. Kuperman, P. L. Nielsen, Bayesian model selection applied to self-noise geoacoustic inversion, J. Acoust. Soc. Am. 116 (4) (2004) 2043-2056.
- [80] E. Fried, M. E. Gurtin, Coherent solid-state phase transitions with atomic diffusion: a thermomechanical treatment, Journal of Statistical Physics 95 (5-6) (1999) 1361-1427.
- [81] M. V. Karsanina, K. M. Gerke, E. B. Skvortsova, D. Mallants, Universal spatial correlation functions for describing and reconstructing soil microstructure, PloS one 10 (5) (2015) e0126515.
- [82] D. T. Fullwood, S. R. Niezgoda, S. R. Kalidindi, Microstructure reconstructions from 2-point statistics using phase-recovery algorithms, Acta Materialia 56 (5) (2008) 942-948.
- [83] J. Shotton, M. Johnson, R. Cipolla, Semantic texton forests for image categorization and segmentation, in: 2008 IEEE Conference on Computer Vision and Pattern Recognition, IEEE, 2008, pp. 1-8.
- [84] G. Csurka, C. Dance, L. Fan, J. Willamowski, C. Bray, Visual categorization with bags of keypoints, in: Workshop on statistical learning in computer vision, ECCV, Vol. 1, Prague, 2004, pp. 1-2.
- [85] D. G. Lowe, Distinctive image features from scale-invariant keypoints, International journal of computer vision 60 (2) (2004) 91-110.
- [86] Y. Ham, M. Kamari, Automated content-based filtering for enhanced vision-based documentation in construction toward exploiting big visual data from drones, Automation in Construction 105 (2019) 102831.
- [87] K. Rajan, Informatics for materials science and engineering: data-driven discovery for accelerated experimentation and application, Butterworth-Heinemann, 2013.
- [88] C. K. Williams, C. E. Rasmussen, Gaussian processes for machinelearning, Vol. 2, MIT Press Cambridge, MA, 2006.
Vahid Attari Pejman Honarmandi Thien Duong Daniel J Sauceda Douglas Allaire Raymundo Arroyave
Uncertainty Propagation in a Multiscale CALPHAD-Reinforced Elastochemical Phase-field Model
ICME approaches provide decision support for materials design by establishing
quantitative process-structure-property relations. Confidence in the decision
support, however, must be achieved by establishing uncertainty bounds in ICME
model chains. The quantification and propagation of uncertainty in
computational materials science, however, remains a rather unexplored aspect of
computational materials science approaches. Moreover, traditional uncertainty
propagation frameworks tend to be limited in cases with computationally
expensive simulations. A rather common and important model chain is that of
CALPHAD-based thermodynamic models of phase stability coupled to phase field
models for microstructure evolution. Propagation of uncertainty in these cases
is challenging not only due to the sheer computational cost of the simulations
but also because of the high dimensionality of the input space. In this work,
we present a framework for the quantification and propagation of uncertainty in
a CALPHAD-based elasto-chemical phase field model. We motivate our work by
investigating the microstructure evolution in Mg$_2$(Si$_x$Sn$_{1-x}$)
thermoelectric materials. We first carry out a Markov Chain Monte Carlo-based
inference of the CALPHAD model parameters for this pseudobinary system and then
use advanced sampling schemes to propagate uncertainties across a
high-dimensional simulation input space. Through high-throughput phase field
simulations we generate 200,000 time series of synthetic microstructures and
use machine learning approaches to understand the effects of propagated
uncertainties on the microstructure landscape of the system under study. The
microstructure dataset has been curated in the Open Phase-field Microstructure
Database (OPMD), available at
\href{http://microstructures.net}{http://microstructures.net}.