Showing 1–39 of 39 results for author: Arroyave, R

Visualizing High Entropy Alloy Spaces: Methods and Best Practices

Authors: Brent Vela, Trevor Hastings, Raymundo Arróyave

Abstract: Multi-Principal Element Alloys (MPEAs) have emerged as an exciting area of research in materials science in the 2020s, owing to the vast potential for discovering alloys with unique and tailored properties enabled by the combinations of elements. However, the chemical complexity of MPEAs poses a significant challenge in visualizing composition-property relationships in high-dimensional design spac… ▽ More Multi-Principal Element Alloys (MPEAs) have emerged as an exciting area of research in materials science in the 2020s, owing to the vast potential for discovering alloys with unique and tailored properties enabled by the combinations of elements. However, the chemical complexity of MPEAs poses a significant challenge in visualizing composition-property relationships in high-dimensional design spaces. Without effective visualization techniques, designing chemically complex alloys is practically impossible. In this methods/protocols article, we present a `toolbox' of visualization techniques that allow for meaningful and insightful visualizations of MPEA composition spaces and property spaces. Our contribution to this toolbox are UMAP projections of entire alloy spaces. We deploy this visualization tool-kit on the following MPEA case studies: 1) Visualizing literature reviews, 2) constraint-satisfaction alloy design scheme, 3) Bayesian optimization alloy design campaigns. Furthermore, we show how this method can be applied to any barycentric design space. While there is no one-size-fits-all visualization technique, our toolbox offers a range of methods and best practices that can be tailored to specific MPEA research needs. This article is intended for materials scientists interested in performing research on multi-principal element alloys, chemically complex alloys, or high entropy alloys, and is expected to facilitate the discovery of novel and tailored properties in MPEAs. △ Less

Submitted 14 August, 2024; originally announced August 2024.

Comments: 24 pages, 17 figures, for associated data availability, see https://doi.org/10.24433/CO.7775216.v1

Illustrating an Effective Workflow for Accelerated Materials Discovery

Authors: Mrinalini Mulukutla, A. Nicole Person, Sven Voigt, Lindsey Kuettner, Branden Kappes, Danial Khatamsaz, Robert Robinson, Daniel Salas, Wenle Xu, Daniel Lewis, Hongkyu Eoh, Kailu Xiao, Haoren Wang, Jaskaran Singh Saini, Raj Mahat, Trevor Hastings, Matthew Skokan, Vahid Attari, Michael Elverud, James D. Paramore, Brady Butler, Kenneth Vecchio, Surya R. Kalidindi, Douglas Allaire, Ibrahim Karaman , et al. (4 additional authors not shown)

Abstract: Algorithmic materials discovery is a multi-disciplinary domain that integrates insights from specialists in alloy design, synthesis, characterization, experimental methodologies, computational modeling, and optimization. Central to this effort is a robust data management system paired with an interactive work platform. This platform should empower users to not only access others data but also inte… ▽ More Algorithmic materials discovery is a multi-disciplinary domain that integrates insights from specialists in alloy design, synthesis, characterization, experimental methodologies, computational modeling, and optimization. Central to this effort is a robust data management system paired with an interactive work platform. This platform should empower users to not only access others data but also integrate their analyses, paving the way for sophisticated data pipelines. To realize this vision, there is a need for an integrative collaboration platform, streamlined data sharing and analysis tools, and efficient communication channels. Such a collaborative mechanism should transcend geographical barriers, facilitating remote interaction and fostering a challenge-response dynamic. In this paper, we present our ongoing efforts in addressing the critical challenges related to an accelerated Materials Discovery Framework as a part of the High-Throughput Materials Discovery for Extreme Conditions Initiative. Our BIRDSHOT Center has successfully harnessed various tools and strategies, including the utilization of cloud-based storage, a standardized sample naming convention, a structured file system, the implementation of sample travelers, a robust sample tracking method, and the incorporation of knowledge graphs for efficient data management. Additionally, we present the development of a data collection platform, reinforcing seamless collaboration among our team members. In summary, this paper provides an illustration and insight into the various elements of an efficient and effective workflow within an accelerated materials discovery framework while highlighting the dynamic and adaptable nature of the data management tools and sharing platforms. △ Less

Submitted 21 May, 2024; originally announced May 2024.

Comments: 28 pages, 9 figures, 2 tables, with appendix that has 8 pages, accepted for publication at IMMI

An Interoperable Multi Objective Batch Bayesian Optimization Framework for High Throughput Materials Discovery

Authors: Trevor Hastings, Mrinalini Mulukutla, Danial Khatamsaz, Daniel Salas, Wenle Xu, Daniel Lewis, Nicole Person, Matthew Skokan, Braden Miller, James Paramore, Brady Butler, Douglas Allaire, Ibrahim Karaman, George Pharr, Ankit Srivastava, Raymundo Arroyave

Abstract: In this study, we introduce a groundbreaking framework for materials discovery, we efficiently navigate a vast phase space of material compositions by leveraging Batch Bayesian statistics in order to achieve specific performance objectives. This approach addresses the challenge of identifying optimal materials from an untenably large array of possibilities in a reasonable timeframe with high confi… ▽ More In this study, we introduce a groundbreaking framework for materials discovery, we efficiently navigate a vast phase space of material compositions by leveraging Batch Bayesian statistics in order to achieve specific performance objectives. This approach addresses the challenge of identifying optimal materials from an untenably large array of possibilities in a reasonable timeframe with high confidence. Crucially, our batchwise methods align seamlessly with existing material processing infrastructure for synthesizing and characterizing materials. By applying this framework to a specific high entropy alloy system, we demonstrate its versatility and robustness in optimizing properties like strain hardening, hardness, and strain rate sensitivity. The fact that the Bayesian model is adept in refining and expanding the property Pareto front highlights its broad applicability across various materials, including steels, shape memory alloys, ceramics, and composites. This study advances the field of materials science and sets a new benchmark for material discovery methodologies. By proving the effectiveness of Bayesian optimization, we showcase its potential to redefine the landscape of materials discovery. △ Less

Submitted 14 May, 2024; originally announced May 2024.

Comments: 12 pages, 6 figures, with Supplementary Appendix that has 17 pages, 9 figures

Two-Shot Optimization of Compositionally Complex Refractory Alloys

Authors: James D. Paramore, Brady G. Butler, Michael T. Hurst, Trevor Hastings, Daniel O. Lewis, Eli Norris, Benjamin Barkai, Joshua Cline, Braden Miller, Jose Cortes, Ibrahim Karaman, George M. Pharr, Raymundo Arroyave

Abstract: In this paper, a synergistic computational/experimental approach is presented for the rapid discovery and characterization of novel alloys within the compositionally complex (i.e., "medium/high entropy") refractory alloy space of Ti-V-Nb-Mo-Hf-Ta-W. This was demonstrated via a material design cycle aimed at simultaneously maximizing the objective properties of high specific hardness (hardness norm… ▽ More In this paper, a synergistic computational/experimental approach is presented for the rapid discovery and characterization of novel alloys within the compositionally complex (i.e., "medium/high entropy") refractory alloy space of Ti-V-Nb-Mo-Hf-Ta-W. This was demonstrated via a material design cycle aimed at simultaneously maximizing the objective properties of high specific hardness (hardness normalized by density) and high specific elastic modulus (elastic modulus normalized by density). This framework utilizes high-throughput computational thermodynamics and intelligent filtering to first reduce the untenably large alloy space to a feasible size, followed by an iterative design cycle comprised of high-throughput synthesis, processing, and characterization in batch sizes of 24 alloys. After the first iteration, Bayesian optimization was utilized to inform selection of the next batch of 24 alloys. This paper demonstrates the benefit of using batch Bayesian optimization (BBO) in material design, as significant gains in the objective properties were observed after only two iterations or "shots" of the design cycle without using any prior knowledge or physical models of how the objective properties relate to the design inputs (i.e., composition). Specifically, the hypervolume of the Pareto front increased by 54% between the first and second iterations. Furthermore, 10 of the 24 alloys in the second iteration dominated all alloys from the first iteration. △ Less

Submitted 11 May, 2024; originally announced May 2024.

Deciphering Chemical Ordering in High Entropy Materials: A Machine Learning-Accelerated High-throughput Cluster Expansion Approach

Authors: Guillermo Vazquez, Daniel Sauceda, Raymundo Arróyave

Abstract: The Cluster Expansion (CE) Method encounters significant computational challenges in multicomponent systems due to the computational expense of generating training data through density functional theory (DFT) calculations. This work aims to refine the cluster and structure selection processes to mitigate these challenges. We introduce a novel method that significantly reduces the computational loa… ▽ More The Cluster Expansion (CE) Method encounters significant computational challenges in multicomponent systems due to the computational expense of generating training data through density functional theory (DFT) calculations. This work aims to refine the cluster and structure selection processes to mitigate these challenges. We introduce a novel method that significantly reduces the computational load associated with the calculation of fitting parameters. This method employs a Graph Neural Network (GNN) model, leveraging the M3GNet network, which is trained using a select subset of DFT calculations at each ionic step. The trained surrogate model excels in predicting the volume and energy of the final structure for a relaxation run. By employing this model, we sample thousands of structures and fit a CE model to the energies of these GNN-relaxed structures. This approach, utilizing a large training dataset, effectively reduces the risk of overfitting, yielding a CE model with a root-mean-square error (RMSE) below 10 meV/atom. We validate our method's effectiveness in two test cases: the (Cr,Hf,Mo,Ta,Ti,Zr)B$_2$ diboride system and the Refractory High-Entropy Alloy (HEA) AlHfNbTaTiZr system. Our findings demonstrate the significant advantages of integrating a GNN model, specifically the M3GNet network, with CE methods for the efficient predictive analysis of chemical ordering in High Entropy Materials. The accelerating capabilities of the hybrid ML-CE approach to investigate the evolution of Short Range Ordering (SRO) in a large number of stoichiometric systems. Finally, we show how it is possible to correlate the strength of chemical ordering to easily accessible alloy parameters. △ Less

Submitted 27 March, 2024; originally announced March 2024.

Comments: 20 pages, 11 figures, presented at TMS 2024

Data-driven study of composition-dependent phase compatibility in NiTi shape memory alloys

Authors: Sina Hossein Zadeh, Cem Cakirhan, Danial Khatamsaz, John Broucek, Timothy D. Brown, Xiaoning Qian, Ibrahim Karaman, Raymundo Arroyave

Abstract: The martensitic transformation in NiTi-based Shape Memory Alloys (SMAs) provides a basis for shape memory effect and superelasticity, thereby enabling applications requiring solid-state actuation and large recoverable shape changes upon mechanical load cycling. In order to tailor the transformation to a particular application, the compositional dependence of properties in NiTi-based SMAs, such as… ▽ More The martensitic transformation in NiTi-based Shape Memory Alloys (SMAs) provides a basis for shape memory effect and superelasticity, thereby enabling applications requiring solid-state actuation and large recoverable shape changes upon mechanical load cycling. In order to tailor the transformation to a particular application, the compositional dependence of properties in NiTi-based SMAs, such as martensitic transformation temperatures and hysteresis, has been exploited. However, the compositional design space is large and complex, and experimental studies are expensive. In this work, we develop an interpretable piecewise linear regression model that predicts the $λ_2$ parameter, a measure of compatibility between austenite and martensite phases, and an (indirect) factor that is well-correlated with martensitic transformation hysteresis, based on the chemical features derived from the alloy composition. The model is capable of predicting, for the first time, the type of martensitic transformation for a given alloy chemistry. The proposed model is validated by experimental data from the literature as well as in-house measurements. The results show that the model can effectively distinguish between $B19$ and $B19^{\prime}$ regions for any given composition in NiTi-based SMAs and accurately estimate the $λ_2$ parameter. Our analysis also reveals that the weighted average of the quotient of the first ionization energy and the Voronoi coordination number is a key compositional characteristic that correlates with the $λ_2$ parameter and thermodynamic responses, including the transformation hysteresis, martensite start temperature, and critical temperature. The work herein demonstrates the potential of data-driven methodologies for understanding and designing NiTi-based SMAs with desired transformation characteristics. △ Less

Submitted 19 February, 2024; originally announced February 2024.

Phase-Field Model of Silicon Carbide Growth During Isothermal Condition

Authors: Elias J. Munoz, Vahid Attari, Marco Martinez, Matthew B. Dickerson, Miladin Radovic, Raymundo Arroyave

Abstract: Silicon carbide (SiC) emerges as a promising ceramic material for high-temperature structural applications, especially within the aerospace sector. The utilization of SiC-based ceramic matrix composites (CMCs) instead of superalloys in components like engine shrouds, combustors, and nozzles offers notable advantages, including a 25% improvement in fuel efficiency, over 10% enhanced thrust, and the… ▽ More Silicon carbide (SiC) emerges as a promising ceramic material for high-temperature structural applications, especially within the aerospace sector. The utilization of SiC-based ceramic matrix composites (CMCs) instead of superalloys in components like engine shrouds, combustors, and nozzles offers notable advantages, including a 25% improvement in fuel efficiency, over 10% enhanced thrust, and the capability to withstand up to 500$^{\circ}$C higher operating temperatures. Employing a CALPHAD-reinforced multi-phase-field model, our study delves into the evolution of the SiC layer under isothermal solidification conditions. By modeling the growth of SiC between liquid Si and solid C at 1450$^{\circ}$C, we compared results with experimental microstructures and quantitatively examined the evolution of SiC thickness over time. Efficient sampling across the entire model space mitigated uncertainty in high-temperature kinetic parameters, allowing us to predict a range of growth rates and morphologies for the SiC layer. The model accounts for parameter uncertainty stemming from limited experimental knowledge and successfully predicts relevant morphologies for the system. Experimental results validated the kinetic parameters of the simulations, offering valuable insights and potential constraints on the reaction kinetics. We further explored the significance of multi-phase-field model parameters on two key outputs, and found that the diffusion coefficient of liquid Si emerges as the most crucial parameter significantly impacting the SiC average layer thickness and grain count over time. This study provides valuable insights into the microstructure evolution of the Si-C binary system, offering pertinent information for the engineering of CMCs in industrial applications. △ Less

Submitted 9 November, 2023; originally announced November 2023.

High-throughput Alloy and Process Design for Metal Additive Manufacturing

Authors: Sofia Sheikh, Brent Vela, Pejman Honarmandi, Peter Morcos, David Shoukr, Abdelrahman Mostafa Kotb, Raymundo Arroyave, Ibrahim Karaman, Alaa Elwany

Abstract: Designing alloys for additive manufacturing (AM) presents significant opportunities. Still, the chemical composition and processing conditions required for printability (ie., their suitability for fabrication via AM) are challenging to explore using solely experimental means. In this work, we develop a high-throughput (HTP) computational framework to guide the search for highly printable alloys an… ▽ More Designing alloys for additive manufacturing (AM) presents significant opportunities. Still, the chemical composition and processing conditions required for printability (ie., their suitability for fabrication via AM) are challenging to explore using solely experimental means. In this work, we develop a high-throughput (HTP) computational framework to guide the search for highly printable alloys and appropriate processing parameters. The framework uses material properties from state-of-the-art databases, processing parameters, and simulated melt pool profiles to predict process-induced defects, such as lack-of-fusion, keyholing, and balling. We accelerate the printability assessment using a deep learning surrogate for a thermal model, enabling a 1,000-fold acceleration in assessing the printability of a given alloy at no loss in accuracy when compared with conventional physics-based thermal models. We verify and validate the framework by constructing printability maps for the CoCrFeMnNi Cantor alloy system and comparing our predictions to an exhaustive 'in-house' database. The framework enables the systematic investigation of the printability of a wide range of alloys in the broader Co-Cr-Fe-Mn-Ni HEA system. We identified the most promising alloys that were suitable for high-temperature applications and had the narrowest solidification ranges, and that was the least susceptible to balling, hot-cracking, and the formation of macroscopic printing defects. A new metric for the global printability of an alloy is constructed and is further used for the ranking of candidate alloys. The proposed framework is expected to be integrated into ICME approaches to accelerate the discovery and optimization of novel high-performance, printable alloys. △ Less

Submitted 8 April, 2023; originally announced April 2023.

Comments: 46 pages, 20 figures

An Automated Fully-Computational Framework to Construct Printability Maps for Additively Manufactured Metal Alloys

Authors: Sofia Sheikh, Meelad Ranaiefar, Pejman Honarmandi, Brent Vela, Peter Morcos, David Shoukr, Raymundo Arroyave, Ibrahim Karaman, Alaa Elwany

Abstract: In additive manufacturing, the optimal processing conditions need to be determined to fabricate porosity-free parts. For this purpose, the design space for an arbitrary alloy needs to be scoped and analyzed to identify the areas of defects for different laser power-scan speed combinations and can be visualized using a printability map. Constructing printability maps is typically a costly process d… ▽ More In additive manufacturing, the optimal processing conditions need to be determined to fabricate porosity-free parts. For this purpose, the design space for an arbitrary alloy needs to be scoped and analyzed to identify the areas of defects for different laser power-scan speed combinations and can be visualized using a printability map. Constructing printability maps is typically a costly process due to the involvement of experiments, which restricts their application in high-throughput product design. To reduce the cost and effort of constructing printability maps, a fully computational framework is introduced in this work. The framework combines CALPHAD models and a reduced-order model to predict material properties. THen, an analytical thermal model, known as the Eagar-Tsai model, utilizes some of these materials' properties to calculate the melt pool geometry during the AM processes. In the end, printability maps are constructed using material properties, melt pool dimensions, and commonly used criteria for lack of fusion, balling, and keyholing defects. To validate the framework and its general application to laser powder-bed fusion alloys, five common additive manufacturing alloys are analyzed. Furthermore, NiTi-based alloys at three different compositions are evaluated to show the further extension of the framework to alloy systems at different compositions. The defect regions in these printability maps are validated with corresponding experimental observations to compare and benchmark the defect criteria and find the optimal criterion set with the maximum accuracy for each unique material composition. Furthermore, printability maps for NiTi that are obtained from our framework are used in conjunction with process maps resulting from a multi-model framework to guide the fabrication of defect-free additive manufactured parts with tailorable properties and performance. △ Less

Submitted 8 April, 2023; originally announced April 2023.

Comments: 18 Figures, 35 pages

An Interpretable Boosting-based Predictive Model for Transformation Temperatures of Shape Memory Alloys

Authors: Sina Hossein Zadeh, Amir Behbahanian, John Broucek, Mingzhou Fan, Guillermo Vazquez Tovar, Mohammad Noroozi, William Trehern, Xiaoning Qian, Ibrahim Karaman, Raymundo Arroyave

Abstract: In this study, we demonstrate how the incorporation of appropriate feature engineering together with the selection of a Machine Learning (ML) algorithm that best suits the available dataset, leads to the development of a predictive model for transformation temperatures that can be applied to a wide range of shape memory alloys. We develop a gradient boosting ML surrogate model capable of predictin… ▽ More In this study, we demonstrate how the incorporation of appropriate feature engineering together with the selection of a Machine Learning (ML) algorithm that best suits the available dataset, leads to the development of a predictive model for transformation temperatures that can be applied to a wide range of shape memory alloys. We develop a gradient boosting ML surrogate model capable of predicting Martensite Start, Martensite Finish, Austenite Start, and Austenite Finish transformation temperatures with an average accuracy of more than 95% by explicitly taking care of potential distribution changes when modeling different alloy systems. We included heat treatment, rolling, extrusion processing parameters, and alloy system categorical features in the model input features to achieve more accurate and realistic results. In addition, using Shapley values, which are calculated based on the average marginal contribution of features to all possible coalitions, this study was able to gain insights into the governing features and their effect on predicted transformation temperatures, providing a unique opportunity to examine the critical parameters and features in martensite transformation temperatures. △ Less

Submitted 4 February, 2023; originally announced February 2023.

A ductility metric for refractory-based multi-principal-element alloys

Authors: Prashant Singh, Brent Vela, Gaoyuan Ouyang, Nicolas Argibay, Jun Cui, Raymundo Arroyave, Duane D. Johnson

Abstract: We propose a quantum-mechanical dimensionless metric, the local$-$lattice distortion (LLD), as a reliable predictor of ductility in refractory multi-principal-element alloys (RMPEAs). The LLD metric is based on electronegativity differences in localized chemical environments and combines atomic$-$scale displacements due to local lattice distortions with a weighted average of valence$-$electron cou… ▽ More We propose a quantum-mechanical dimensionless metric, the local$-$lattice distortion (LLD), as a reliable predictor of ductility in refractory multi-principal-element alloys (RMPEAs). The LLD metric is based on electronegativity differences in localized chemical environments and combines atomic$-$scale displacements due to local lattice distortions with a weighted average of valence$-$electron count. To evaluate the effectiveness of this metric, we examined body$-$centered cubic (bcc) refractory alloys that exhibit ductile$-$to$-$brittle behavior. Our findings demonstrate that local$-$charge behavior can be tuned via composition to enhance ductility in RMPEAs. With finite$-$sized cell effects eliminated, the LLD metric accurately predicted the ductility of arbitrary alloys based on tensile$-$elongation experiments. To validate further, we qualitatively evaluated the ductility of two refractory RMPEAs, i.e., NbTaMoW and Mo$_{72}$W$_{13}Ta$_{10}Ti$_{2.5}Zr$_{2.5}, through the observation of crack formation under indentation, again showing excellent agreement with LLD predictions. A comparative study of three refractory alloys provides further insights into the electronic-structure origin of ductility in refractory RMPEAs. This proposed metric enables rapid and accurate assessment of ductility behavior in the vast RMPEA composition space. △ Less

Submitted 26 June, 2023; v1 submitted 28 November, 2022; originally announced November 2022.

Comments: 36 pages, 12 figures, 5 Table

On the Effect of Nucleation Undercooling on Phase Transformation Kinetics

Authors: José Mancias, Vahid Attari, Raymundo Arróyave, Damien Tourret

Abstract: We carry out an extensive comparison between Johnson-Mehl-Avrami-Kolmogorov (JMAK) theory of first-order phase transformation kinetics and phase-field (PF) results of a benchmark problem on nucleation. To address the stochasticity of the problem, several hundreds of simulations are performed to establish a comprehensive, statistically-significant analysis of the coincidences and discrepancies betw… ▽ More We carry out an extensive comparison between Johnson-Mehl-Avrami-Kolmogorov (JMAK) theory of first-order phase transformation kinetics and phase-field (PF) results of a benchmark problem on nucleation. To address the stochasticity of the problem, several hundreds of simulations are performed to establish a comprehensive, statistically-significant analysis of the coincidences and discrepancies between PF and JMAK transformation kinetics. We find that PF predictions are in excellent agreement with both classical nucleation theory and JMAK theory, as long as the original assumptions of the latter are appropriately reproduced - in particular, the constant nucleation and growth rates in an infinite domain. When deviating from these assumptions, PF results are at odds with JMAK theory. In particular, we observe that the size of the initial particle radius $r_0$ relative to the critical nucleation radius $r^*$ has a significant effect on the rate of transformation. While PF and JMAK agree when $r_0$ is sufficiently higher than $r^*$, the duration of initial transient growth stage of a particle, before it reaches a steady growth velocity, increases as $r_0/r^*\to 1$. This incubation time has a significant effect on the overall kinetics, e.g. on the Avrami exponent of the multi-particle simulations. In contrast, for the considered conditions and parameters, the effect of interface curvature upon transformation kinetics - in particular negative curvature regions appearing during particle impingement, present in PF but absent in JMAK theory - appears to be minor compared to that of $r_0/r^*$. We argue that rigorous benchmarking of phase-field models of stochastic processes (e.g. nucleation) need sufficient statistical data in order to make rigorous comparisons against ground truth theories. In these problems, analysis of probability distributions is clearly preferable to a deterministic approach. △ Less

Submitted 30 October, 2022; originally announced October 2022.

Journal ref: Integrating Materials and Manufacturing Innovation 11 (2022) 628-636

Efficient machine-learning model for fast assessment of elastic properties of high-entropy alloys

Authors: Guillermo Vazquez, Prashant Singh, Daniel Sauceda, Richard Couperthwaite, Nicholas Britt, Khaled Youssef, Duane D. Johnson, Raymundo Arróyave

Abstract: We combined descriptor-based analytical models for stiffness-matrix and elastic-moduli with mean-field methods to accelerate assessment of technologically useful properties of high-entropy alloys, such as strength and ductility. Model training for elastic properties uses Sure-Independence Screening (SIS) and Sparsifying Operator (SO) method yielding an optimal analytical model, constructed with me… ▽ More We combined descriptor-based analytical models for stiffness-matrix and elastic-moduli with mean-field methods to accelerate assessment of technologically useful properties of high-entropy alloys, such as strength and ductility. Model training for elastic properties uses Sure-Independence Screening (SIS) and Sparsifying Operator (SO) method yielding an optimal analytical model, constructed with meaningful atomic features to predict target properties. Computationally inexpensive analytical descriptors were trained using a database of the elastic properties determined from density functional theory for binary and ternary subsets of Nb-Mo-Ta-W-V refractory alloys. The optimal Elastic-SISSO models, extracted from an exponentially large feature space, give an extremely accurate prediction of target properties, similar to or better than other models, with some verified from existing experiments. We also show that electronegativity variance and elastic-moduli can directly predict trends in ductility and yield strength of refractory HEAs, and reveals promising alloy concentration regions. △ Less

Submitted 4 April, 2022; originally announced April 2022.

Comments: 17 Pages; 10 Figures; 3 Tables

Machine-learning enabled thermodynamic model for the design of new rare-earth compounds

Authors: Prashant Singh, Tyler Del Rose, Guillermo Vazquez, Raymundo Arroyave, Yaroslav Mudryk

Abstract: We employ a descriptor based machine-learning approach to assess the effect of chemical alloying on formation-enthalpy of rare-earth intermetallics. Application of machine-learning approaches in rare-earth intermetallic design have been sparse due to limited availability of reliable datasets. In this work, we developed an `in-house' rare-earth database with more than 600$+$ compounds, each entry w… ▽ More We employ a descriptor based machine-learning approach to assess the effect of chemical alloying on formation-enthalpy of rare-earth intermetallics. Application of machine-learning approaches in rare-earth intermetallic design have been sparse due to limited availability of reliable datasets. In this work, we developed an `in-house' rare-earth database with more than 600$+$ compounds, each entry was populated with formation enthalpy and related atomic features using high-throughput density-functional theory (DFT). The SISSO (sure independence screening and sparsifying operator) based machine-learning method with meaningful atomic features was used for training and testing the formation enthalpies of rare earth compounds. The complex lattice function coupled with the machine-learning model was used to explore the effect of transition metal alloying on the energy stability of Ce based cubic Laves phases (MgCu$_{2}$ type). The SISSO predictions show good agreement with high-fidelity DFT calculations and X$-$ray powder diffraction measurements. Our study provides quantitative guidance for compositional considerations within a machine-learning model and discovering new metastable materials. The electronic-structure of Ce$-$Fe$-$Cu based compound was also analyzed in$-$depth to understand the electronic origin of phase stability. The interpretable analytical models in combination with density$-$functional theory and experiments provide a fast and reliable design guide for discovering technologically useful materials. △ Less

Submitted 4 March, 2022; originally announced March 2022.

Comments: 19 pages, 9 figures

Towards Stacking Fault Energy Engineering in FCC High Entropy Alloys

Authors: Tasneem Khan, Tanner Kirk, Guillermo Vazquez, Prashant Singh, A V Smirnov, Duane D Johnson, Khaleed Youssef, Raymundo Arroyave

Abstract: Stacking Fault Energy (SFE) is an intrinsic alloy property that governs much of the plastic deformation mechanisms observed in fcc alloys. While SFE has been recognized for many years as a key intrinsic mechanical property, its inference via experimental observations or prediction using, for example, computationally intensive first-principles methods is challenging. This difficulty precludes the e… ▽ More Stacking Fault Energy (SFE) is an intrinsic alloy property that governs much of the plastic deformation mechanisms observed in fcc alloys. While SFE has been recognized for many years as a key intrinsic mechanical property, its inference via experimental observations or prediction using, for example, computationally intensive first-principles methods is challenging. This difficulty precludes the explicit use of SFE as an alloy design parameter. In this work, we combine DFT calculations (with necessary configurational averaging), machine-learning (ML) and physics-based models to predict the SFE in the fcc CoCrFeMnNiV-Al high-entropy alloy space. The best-performing ML model is capable of accurately predicting the SFE of arbitrary compositions within this 7-element system. This efficient model along with a recently developed model to estimate intrinsic strength of fcc HEAs is used to explore the strength-SFE Pareto front, predicting new-candidate alloys with particularly interesting mechanical behavior. △ Less

Submitted 5 November, 2021; originally announced November 2021.

Comments: 27 pages, 16 figures, 1 Table

Pseudoelastic deformation in Mo-based refractory multi-principal element alloys

Authors: Aayush Sharma, Prashant Singh, Tanner Kirk, Velary I. Levitas, Peter K. Liaw, Ganesh Balasubramanian, Raymundo Arroyave, Duane D Johnson

Abstract: Phase diagrams supported by density functional theory methods can be crucial for designing high-entropy alloys that are subset of multi-principal$-$element alloys. We present phase and property analysis of quinary (MoW)$_{x}$Zr$_{y}$(TaTi)$_{1-x-y}$ refractory high-entropy alloys from combined Calculation of Phase Diagram (CALPHAD) and density-functional theory results, supplemented by molecular d… ▽ More Phase diagrams supported by density functional theory methods can be crucial for designing high-entropy alloys that are subset of multi-principal$-$element alloys. We present phase and property analysis of quinary (MoW)$_{x}$Zr$_{y}$(TaTi)$_{1-x-y}$ refractory high-entropy alloys from combined Calculation of Phase Diagram (CALPHAD) and density-functional theory results, supplemented by molecular dynamics simulations. Both CALPHAD and density-functional theory analysis of phase stability indicates a Mo-W-rich region of this quinary has a stable single-phase body-centered-cubic structure. We report first quinary composition from Mo$-$W$-$Ta$-$Ti$-$Zr family of alloy with pseudo-elastic behavior, i.e., hysteresis in stress$-$strain. Our analysis shows that only Mo$-$W$-$rich compositions of Mo$-$W$-$Ta$-$Ti$-$Zr, i.e., Mo$+$W$\ge$ 85 at.%, show reproducible hysteresis in stress-strain responsible for pseudo-elastic behavior. The (MoW)$_{85}$Zr$_{7.5}$(TaTi)$_{7.5}$ was down-selected based on temperature-dependent phase diagram analysis and molecular dynamics simulations predicted elastic behavior that reveals twinning assisted pseudoelastic behavior. While mostly unexplored in body-centered-cubic crystals, twinning is a fundamental deformation mechanism that competes against dislocation slip in crystalline solids. This alloy shows identical cyclic deformation characteristics during uniaxial $\lt$100$\gt$ loading, i.e., the pseudoelasticity is isotropic in loading direction. Additionally, a temperature increase from 77 to 1500 K enhances the elastic strain recovery in load-unload cycles, offering possibly control to tune the pseudoelastic behavior. △ Less

Submitted 5 September, 2021; originally announced September 2021.

Comments: 29 pages, 13 Figures (9 main figures, 4 supplement figures)

On the martensitic transformation in Fe$_{x}$Mn$_{80-x}$Co$_{10}$Cr$_{10}$ high-entropy alloy

Authors: Prashant Singh, Sezer Picak, Aayush Sharma, Y. I. Chumlyakov, Raymundo Arroyave, Ibrahim Karaman, Duane D. Johnson

Abstract: High-entropy alloys (HEAs), and even medium-entropy alloys (MEAs), are an intriguing class of materials in that structure and property relations can be controlled via alloying and chemical disorder over wide ranges in the composition space. Employing density-functional theory combined with the coherent-potential approximation to average over all chemical configurations, we tune free energies betwe… ▽ More High-entropy alloys (HEAs), and even medium-entropy alloys (MEAs), are an intriguing class of materials in that structure and property relations can be controlled via alloying and chemical disorder over wide ranges in the composition space. Employing density-functional theory combined with the coherent-potential approximation to average over all chemical configurations, we tune free energies between face-centered-cubic (fcc) and hexagonal-close-packed (hcp) phases in Fe$_{x}$Mn$_{80-x}$Co$_{10}$Cr$_{10}$ systems.~Within Fe-Mn-based alloys, we show that the martensitic transformation and chemical short-range order directly correlate with the fcc-hcp energy difference and stacking-fault energies, which are in quantitative agreement with recent experiments on a $x$=40~at.\% polycrystalline HEA/MEA. Our predictions are further confirmed by single-crystal measurements on a$x$=40at.\% using transmission-electron microscopy, selective-area diffraction, and electron-backscattered-diffraction mapping. The results herein offer an understanding of transformation-induced/twinning-induced plasticity (TRIP/TWIP) in this class of HEAs and a design guide for controlling the physics behind the TRIP effect at the electronic level. △ Less

Submitted 3 August, 2021; originally announced August 2021.

Comments: 6 pages, 5 figures, 1 table

Metric-driven search for structurally stable inorganic compounds

Authors: R. Villarreal, P. Singh, R. Arroyave

Abstract: We report a facile `metric' for the identification of structurally and dynamically (positive definite phonon structure) stable inorganic compounds. The metric considers charge-imbalance within the local substructures in crystalline compounds calculated using first-principles density-functional theory. To exemplify, we chose carbon-based nitrides as it provides a large pool of structurally stable a… ▽ More We report a facile `metric' for the identification of structurally and dynamically (positive definite phonon structure) stable inorganic compounds. The metric considers charge-imbalance within the local substructures in crystalline compounds calculated using first-principles density-functional theory. To exemplify, we chose carbon-based nitrides as it provides a large pool of structurally stable and unstable phases. We showcase how local structural information related to Wyckoff symmetry uniquely identifies four new carbon-nitride phases. The metric predicts three new structurally stable phases of 4 (C)$:$3 (N) stoichiometry, i.e., C$_{4}$N$_{3}$, which is further confirmed by direct phonon calculations. The structurally stable phases also satisfy the thermodynamic stability and mechanical stability criteria. New phases possess extraordinary mechanical properties, similar to diamond, and show insulating bandgap ranging from optimal, 1.45 eV (optimal) to 5.5 eV (large gap). The structural, electronic, and optical properties of Pm(1)-C$_{4}$N$_{3}$, discussed in detail, indicate possible application in optoelectronic and photovoltaic technologies. The way metric is design, it can be used to predict stability of any material with three-dimensional bonding network. The metric was also able to predict structural stability of other nitride polymorphs with 100$\%$ accuracy. We believe that the proposed `Metric' will accelerate the search of unknown and unexplored inorganic compounds by quick filtering of dynamically stable phases without expensive density-functional theory calculations. △ Less

Submitted 17 November, 2020; originally announced November 2020.

Comments: 37 Pages, 9 Figures, 5 Tables

Journal ref: Acta Mater. 202 (2021) 437-447

Impact of Particle Arrays on Phase Separation Composition Patterns

Authors: Supriyo Ghosh, Arnab Mukherjee, Raymundo Arroyave, Jack F. Douglas

Abstract: We examine the symmetry-breaking effect of fixed constellations of particles on the surface-directed spinodal decomposition of binary blends in the presence of particles whose surfaces have a preferential affinity for one of the components. Our phase-field simulations indicate that the phase separation morphology in the presence of particle arrays can be tuned to have a continuous, droplet, lamell… ▽ More We examine the symmetry-breaking effect of fixed constellations of particles on the surface-directed spinodal decomposition of binary blends in the presence of particles whose surfaces have a preferential affinity for one of the components. Our phase-field simulations indicate that the phase separation morphology in the presence of particle arrays can be tuned to have a continuous, droplet, lamellar, or hybrid morphology depending on the interparticle spacing, blend composition, and time. In particular, when the interparticle spacing is large compared to the spinodal wavelength, a transient target pattern composed of alternate rings of preferred and non-preferred phases emerge at early times, tending to adopt the symmetry of the particle configuration. We reveal that such target patterns stabilize for certain characteristic length, time, and composition scales characteristic of the pure phase separating mixture. To illustrate the general range of phenomena exhibited by mixture-particle systems, we simulate the effects of single-particle, multi-particle, and cluster-particle systems having multiple geometrical configurations of the particle characteristic of pattern substrates on phase separation. Our simulations show that tailoring the particle configuration, or substrate pattern configuration, a relative fluid-particle composition should allow the desirable control of the phase separation morphology as in block copolymer materials, but where the scales accessible to this approach of organizing phase-separated fluids usually are significantly larger. Limited experiments confirm the trends observed in our simulations, which should provide some guidance in engineering patterned blend and other mixtures of technological interest. △ Less

Submitted 10 June, 2020; originally announced June 2020.

Journal ref: The Journal of Chemical Physics, 152, 224902 (2020)

Accelerated design of Fe-based soft magnetic materials using machine learning and stochastic optimization

Authors: Yuhao Wang, Yefan Tian, Tanner Kirk, Omar Laris, Joseph H. Ross, Jr., Ronald D. Noebe, Vladimir Keylin, Raymundo Arróyave

Abstract: Machine learning was utilized to efficiently boost the development of soft magnetic materials. The design process includes building a database composed of published experimental results, applying machine learning methods on the database, identifying the trends of magnetic properties in soft magnetic materials, and accelerating the design of next-generation soft magnetic nanocrystalline materials t… ▽ More Machine learning was utilized to efficiently boost the development of soft magnetic materials. The design process includes building a database composed of published experimental results, applying machine learning methods on the database, identifying the trends of magnetic properties in soft magnetic materials, and accelerating the design of next-generation soft magnetic nanocrystalline materials through the use of numerical optimization. Machine learning regression models were trained to predict magnetic saturation ($B_S$), coercivity ($H_C$) and magnetostriction ($λ$), with a stochastic optimization framework being used to further optimize the corresponding magnetic properties. To verify the feasibility of the machine learning model, several optimized soft magnetic materials -- specified in terms of compositions and thermomechanical treatments -- have been predicted and then prepared and tested, showing good agreement between predictions and experiments, proving the reliability of the designed model. Two rounds of optimization-testing iterations were conducted to search for better properties. △ Less

Submitted 4 May, 2020; v1 submitted 4 February, 2020; originally announced February 2020.

Journal ref: Acta Mater. 194, 144 (2020)

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

Authors: Vahid Attari, Pejman Honarmandi, Thien Duong, Daniel J Sauceda, Douglas Allaire, Raymundo Arroyave

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 mater… ▽ More 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}. △ Less

Submitted 1 August, 2019; originally announced August 2019.

The Effect of Chemical Disorder on Defect Formation and Migration in Disordered MAX Phases

Authors: Prashant Singh, Daniel Sauceda, Raymundo Arroyave

Abstract: MAX phases have attracted increased attention due to their unique combination of ceramic and metallic properties. Point-defects are known to play a vital role in the structural, electronic and transport properties of alloys in general and this system in particular. As some MAX phases have been shown to be stable in non-stoichiometric compositions, it is likely that such alloying effects will affec… ▽ More MAX phases have attracted increased attention due to their unique combination of ceramic and metallic properties. Point-defects are known to play a vital role in the structural, electronic and transport properties of alloys in general and this system in particular. As some MAX phases have been shown to be stable in non-stoichiometric compositions, it is likely that such alloying effects will affect the behavior of lattice point defects. This problem, however, remains relatively unexplored. In this work, we investigate the alloying effect on the structural-stability, energy-stability, electronic-structure, and diffusion barrier of point defects in MAX phase alloys within a first-principles density functional theory framework. The vacancy (V$_{M}$, V$_{A}$, V$_{X}$) and antisite (M-A; M-X) defects are considered with M and A site disorder in (Zr-M)$_{2}$(AA${'}$)C, where M=Cr,Nb,Ti and AA${'}$=Al, Al-Sn, Pb-Bi. Our calculations suggest that the chemical disorder helps lower the V$_{A}$ formation energies compared to V$_{M}$ and V$_{X}$. The V$_{A}$ diffusion barrier is also significantly reduced for M-site disorder compared to their ordered counterpart. This is very important finding because reduced barrier height will ease the Al diffusion at high-operating temperatures, which will help the formation of passivating oxide layer (i.e., Al$_{2}$O$_{3}$ in aluminum-based MAX phases) and will slow down or stop the material degradation. We believe that our study will provide a fundamental understanding and an approach to tailor the key properties that can lead to the discovery of new MAX phases. △ Less

Submitted 24 October, 2019; v1 submitted 31 July, 2019; originally announced August 2019.

Comments: 26 pages, 18 figures, 2 Tables

Probing discontinuous precipitation in U-Nb

Authors: Thien Duong, Robert E. Hackenberg, Vahid Attari, Alex Landa, Patrice E. A. Turchi, Raymundo Arroyave

Abstract: U-Nb's discontinuous precipitation, $γ^{bcc}_{matrix} \rightarrow α^{orth}_{cellular} + γ'^{bcc}_{cellular}$, is intriguing in the sense that it allows formation and growth of the metastable $γ'$ phase during the course of its occurrence. Previous attempts to explain the thermodynamic origin of U-Nb's discontinuous precipitation hypothesized that the energy of $α$ forms an intermediate common tang… ▽ More U-Nb's discontinuous precipitation, $γ^{bcc}_{matrix} \rightarrow α^{orth}_{cellular} + γ'^{bcc}_{cellular}$, is intriguing in the sense that it allows formation and growth of the metastable $γ'$ phase during the course of its occurrence. Previous attempts to explain the thermodynamic origin of U-Nb's discontinuous precipitation hypothesized that the energy of $α$ forms an intermediate common tangent with the first potential of the double-well energy of $γ$ at the $γ'$ composition. While this hypothesis is eligible and consistent with the experimental observation of gradual increase in $γ'$ composition at increasing temperature, it is in conflict with recent experiments whose results indicated a distribution of $γ'$ compositions in the vicinity of 50 at\%Nb. To shed some light onto this issue, the current work investigates the origin of U-Nb's discontinuous precipitation in view of fundamental thermodynamics and kinetics, taken from the perspective of phase-field theory. It has been showed that local misfit strain tends to play a crucial role in the formation and growth the discontinuous precipitation. Depending on the magnitude of strain developed at grain boundaries, either an increasing $γ'$ composition or a random distribution of $γ'$ composition around the equiatomic value with respect to increasing temperature could be expected. △ Less

Submitted 1 July, 2019; originally announced July 2019.

Finite Interface Dissipation Phase Field Modeling of Ni-Nb under Additive Manufacturing Conditions

Authors: Kubra Karayagiz, Luke Johnson, Raiyan Seede, Vahid Attari, Bing Zhang, Xueqin Huang, Supriyo Ghosh, Thien Duong, Ibrahim Karaman, Alaa Elwany, Raymundo Arroyave

Abstract: During the laser powder bed fusion (L-PBF) process, the built part undergoes multiple rapid heating-cooling cycles, leading to complex microstructures with nonuniform properties. In the present work, a computational framework, which weakly couples a finite element thermal model to a non-equilibrium PF model was developed to investigate the rapid solidification microstructure of a Ni-Nb alloy durin… ▽ More During the laser powder bed fusion (L-PBF) process, the built part undergoes multiple rapid heating-cooling cycles, leading to complex microstructures with nonuniform properties. In the present work, a computational framework, which weakly couples a finite element thermal model to a non-equilibrium PF model was developed to investigate the rapid solidification microstructure of a Ni-Nb alloy during L-PBF. The framework is utilized to predict the spatial variation of the morphology and size of cellular segregation structure as well as the microsegregation in single-track melt pool microstructures obtained under different process conditions. A solidification map demonstrating the variation of microstructural features as a function of the temperature gradient and growth rate is presented. A planar to cellular transition is predicted in the majority of keyhole mode melt pools, while a planar interface is predominant in conduction mode melt pools. The predicted morphology and size of the cellular segregation structure agrees well with experimental measurements. △ Less

Submitted 19 February, 2020; v1 submitted 24 June, 2019; originally announced June 2019.

Exploration of the Microstructure Space in TiAlZrN Ultra-Hard Nanostructured Coatings

Authors: Vahid Attari, Aitor Cruzado, Raymundo Arroyave

Abstract: Ti$_{1-x-y}$Al$_{x}$Zr$_{y}$N cubic alloys within the 25-70\% Al composition range have high age-hardening capabilities due to metastable phase transition pathways at high temperatures. They are thus ideal candidates for ultra-hard nano-coating materials. There is growing evidence that this effect is associated with the elasto-chemical field-induced phase separation into compositionally-segregated… ▽ More Ti$_{1-x-y}$Al$_{x}$Zr$_{y}$N cubic alloys within the 25-70\% Al composition range have high age-hardening capabilities due to metastable phase transition pathways at high temperatures. They are thus ideal candidates for ultra-hard nano-coating materials. There is growing evidence that this effect is associated with the elasto-chemical field-induced phase separation into compositionally-segregated nanocrystaline nitride phases. Here, we studied the microstructural evolution in this pseudo-ternary system within spinodal regions at 1200 $^\circ$C by using an elasto-chemical phase field model. Our simulations indicate that elastic interactions between nitride nano-domains greatly affect not only the morphology of the microstructure but also the local chemical phase equilibria. In Al-rich regions of the composition space we further observe the onset of the transformation of AlN-rich phases into their equilibrium wurtzite crystal structure. This work points to a wide palette of microstructures potentially accessible to these nitride systems and their tailoring is likely to result in significant improvements in the performance of transition metal nitride-based coating materials. △ Less

Submitted 12 November, 2018; originally announced November 2018.

Phase stability and elastic properties in the $Al_{1-x-y}Cr_{x}Ti_{y}N$ system from first principles

Authors: Erik Gutiérrez-Valladares, Rurick Santos-Fragoso, Guillermo Vázquez-Tovar, Andrés Manuel Garay-Tapia, Diego Germán Espinosa-Arbeláez, Raymundo Arróyave, Jesús González-Hernández, Juan Manuel Alvarado-Orozco

Abstract: Multicomponent nitrides are a hot research topic in the search of hard coatings. The effect of substitutions on the phase stabilities, magnetic, and elastic properties of $Al_{1-x-y}Cr_{x}Ti_{y}N$ $(0\leq x,y\leq1)$ was studied using first principles calculations based on the density functional theory. These calculations provide the lattice parameter, formation energy, mixing enthalpy and elastic… ▽ More Multicomponent nitrides are a hot research topic in the search of hard coatings. The effect of substitutions on the phase stabilities, magnetic, and elastic properties of $Al_{1-x-y}Cr_{x}Ti_{y}N$ $(0\leq x,y\leq1)$ was studied using first principles calculations based on the density functional theory. These calculations provide the lattice parameter, formation energy, mixing enthalpy and elastic constants. The calculated values are in good agreement with experiments and compare well with other theoretical results. A magnetic transition from antiferromagnetic to ferromagnetic state occurs at concentrations of B1-TiN higher than 60$\%$. The quaternary zone has a lower aluminum solubility than the constituent ternary systems. The Poisson's ratio, Shear and Young modulus were used as predictors of the hardness, indicating that the higher hardness values are found on the transition line from B1 to B4. The obtained results enable the design of new $Al_{1-x-y}Cr_{x}Ti_{y}N$-based materials for coating applications. △ Less

Submitted 12 October, 2018; originally announced October 2018.

Multi-Information Source Fusion and Optimization to Realize ICME: Application to Dual Phase Materials

Authors: Seyede Fatemeh Ghoreishi, Abhilash Molkeri, Ankit Srivastava, Raymundo Arroyave, Douglas Allaire

Abstract: Integrated Computational Materials Engineering (ICME) calls for the integration of computational tools into the materials and parts development cycle, while the Materials Genome Initiative (MGI) calls for the acceleration of the materials development cycle through the combination of experiments, simulation, and data. As they stand, both ICME and MGI do not prescribe how to achieve the necessary to… ▽ More Integrated Computational Materials Engineering (ICME) calls for the integration of computational tools into the materials and parts development cycle, while the Materials Genome Initiative (MGI) calls for the acceleration of the materials development cycle through the combination of experiments, simulation, and data. As they stand, both ICME and MGI do not prescribe how to achieve the necessary tool integration or how to efficiently exploit the computational tools, in combination with experiments, to accelerate the development of new materials and materials systems. This paper addresses the first issue by putting forward a framework for the fusion of information that exploits correlations among sources/models and between the sources and `ground truth'. The second issue is addressed through a multi-information source optimization framework that identifies, given current knowledge, the next best information source to query and where in the input space to query it via a novel value-gradient policy. The querying decision takes into account the ability to learn correlations between information sources, the resource cost of querying an information source, and what a query is expected to provide in terms of improvement over the current state. The framework is demonstrated on the optimization of a dual-phase steel to maximize its strength-normalized strain hardening rate. The ground truth is represented by a microstructure-based finite element model while three low fidelity information sources---i.e. reduced order models---based on different homogenization assumptions---isostrain, isostress and isowork---are used to efficiently and optimally query the materials design space. △ Less

Submitted 17 July, 2018; originally announced September 2018.

Comments: 19 pages, 11 figures, 5 tables

Parametric Analysis of a Phenomenological Constitutive Model for Thermally Induced Phase Transformation in Ni-Ti Shape Memory Alloys

Authors: Pejman Honarmandi, Alex Solomou, Raymundo Arroyave, Dimitris Lagoudas

Abstract: In this work, a thermo-mechanical model that predicts the actuation response of shape memory alloys is probabilistically calibrated against three experimental data sets simultaneously. Before calibration, a design of experiments (DOE) has been performed in order to identify the parameters most influential on the actuation response of the system and thus reduce the dimensionality of the problem. Su… ▽ More In this work, a thermo-mechanical model that predicts the actuation response of shape memory alloys is probabilistically calibrated against three experimental data sets simultaneously. Before calibration, a design of experiments (DOE) has been performed in order to identify the parameters most influential on the actuation response of the system and thus reduce the dimensionality of the problem. Subsequently, uncertainty quantification (UQ) of the influential parameters was carried out through Bayesian Markov Chain Monte Carlo (MCMC). The assessed uncertainties in the model parameters were then propagated to the transformation strain-temperature hysteresis curves (the model output) using first an approximate approach based on the variance-covariance matrix of the MCMC-calibrated model parameters and then an explicit propagation of uncertainty through MCMC-based sampling. Results show good agreement between model and experimental hysteresis loops after probabilistic MCMC calibration such that the experimental data are situated within 95% Bayesian confidence intervals. The application of the MCMC-based UQ/UP approach in decision making for experimental design has also been shown by comparing the information that can be gained from running replicas around a single new experimental condition versus running experiments in different regions of the experimental space. △ Less

Submitted 10 August, 2018; originally announced August 2018.

Comments: 26 pages, 5 figures, 4 tables

Multi-Objective Bayesian Materials Discovery: Application on the Discovery of Precipitation Strengthened NiTi Shape Memory Alloys through Micromechanical Modeling

Authors: Alexandros Solomou, Guang Zhao, Shahin Boluki, Jobin K. Joy, Xiaoning Qian, Ibrahim Karaman, Raymundo Arróyave, Dimitris C. Lagoudas

Abstract: In this study, a framework for the multi-objective materials discovery based on Bayesian approaches is developed. The capabilities of the framework are demonstrated on an example case related to the discovery of precipitation strengthened NiTi shape memory alloys with up to three desired properties. In the presented case the framework is used to carry out an efficient search of the shape memory al… ▽ More In this study, a framework for the multi-objective materials discovery based on Bayesian approaches is developed. The capabilities of the framework are demonstrated on an example case related to the discovery of precipitation strengthened NiTi shape memory alloys with up to three desired properties. In the presented case the framework is used to carry out an efficient search of the shape memory alloys with desired properties while minimizing the required number of computational experiments. The developed scheme features a Bayesian optimal experimental design process that operates in a closed loop. A Gaussian process regression model is utilized in the framework to emulate the response and uncertainty of the physical/computational data while the sequential exploration of the materials design space is carried out by using an optimal policy based on the expected hyper-volume improvement acquisition function. This scalar metric provides a measure of the utility of querying the materials design space at different locations, irrespective of the number of objectives in the performed task. The framework is deployed for the determination of the composition and microstructure of precipitation-strengthened NiTi shape memory alloys with desired properties, while the materials response as a function of microstructure is determined through a thermodynamically-consistent micromechanical model. △ Less

Submitted 18 July, 2018; originally announced July 2018.

Bayesian Uncertainty Quantification and Information Fusion in CALPHAD-based Thermodynamic Modeling

Authors: Pejman Honarmandi, Thien Chi Duong, Seyede Fatemeh Ghoreishi, Douglas Allaire, Raymundo Arroyave

Abstract: Calculation of phase diagrams is one of the fundamental tools in alloy design---more specifically under the framework of Integrated Computational Materials Engineering. Uncertainty quantification of phase diagrams is the first step required to provide confidence for decision making in property- or performance-based design. As a manner of illustration, a thorough probabilistic assessment of the CAL… ▽ More Calculation of phase diagrams is one of the fundamental tools in alloy design---more specifically under the framework of Integrated Computational Materials Engineering. Uncertainty quantification of phase diagrams is the first step required to provide confidence for decision making in property- or performance-based design. As a manner of illustration, a thorough probabilistic assessment of the CALPHAD model parameters is performed against the available data for a Hf-Si binary case study using a Markov Chain Monte Carlo sampling approach. The plausible optimum values and uncertainties of the parameters are thus obtained, which can be propagated to the resulting phase diagram. Using the parameter values obtained from deterministic optimization in a computational thermodynamic assessment tool (in this case Thermo-Calc) as the prior information for the parameter values and ranges in the sampling process is often necessary to achieve a reasonable cost for uncertainty quantification. This brings up the problem of finding an appropriate CALPHAD model with high-level of confidence which is a very hard and costly task that requires considerable expert skill. A Bayesian hypothesis testing based on Bayes' factors is proposed to fulfill the need of model selection in this case, which is applied to compare four recommended models for the Hf-Si system. However, it is demonstrated that information fusion approaches, i.e., Bayesian model averaging and an error correlation-based model fusion, can be used to combine the useful information existing in all the given models rather than just using the best selected model, which may lack some information about the system being modelled. △ Less

Submitted 18 July, 2018; v1 submitted 12 June, 2018; originally announced June 2018.

Comments: 22 pages, 8 Figures

On the Interfacial Phase Growth and Vacancy Evolution during Accelerated Electromigration in Cu/Sn/Cu Microjoints

Authors: Vahid Attari, Supriyo Ghosh, Thien Duong, Raymundo Arroyave

Abstract: In this work, we integrate different computational tools based on multi-phase-field simulations to account for the evolution of morphologies and crystallographic defects of Cu/Sn/Cu sandwich interconnect structures that are widely used in three dimensional integrated circuits (3DICs). Specifically, this work accounts for diffusion-driven formation and disappearance of multiple intermetallic phases… ▽ More In this work, we integrate different computational tools based on multi-phase-field simulations to account for the evolution of morphologies and crystallographic defects of Cu/Sn/Cu sandwich interconnect structures that are widely used in three dimensional integrated circuits (3DICs). Specifically, this work accounts for diffusion-driven formation and disappearance of multiple intermetallic phases during accelerated electromigration and takes into account the non-equilibrium formation of vacancies due to electromigration. The work compares nucleation, growth, and coalescence of intermetallic layers during transient liquid phase bonding and virtual joint structure evolution subjected to accelerated electromigration conditions at different temperatures. The changes in the rate of dissolution of Cu from intermetallics and the differences in the evolution of intermetallic layers depending on whether they act as cathodes or anodes are accounted for and are compared favorably with experiments. The model considers non-equilibrium evolution of vacancies that form due to differences in couplings between diffusing atoms and electron flows. This work is significant as the point defect evolution in 3DIC solder joints during electromigration has deep implications to the formation and coalescence of voids that ultimately compromise the structural and functional integrity of the joints. △ Less

Submitted 23 October, 2018; v1 submitted 19 March, 2018; originally announced March 2018.

Comments: 15 pages, 7 figures

Journal ref: Attari, Vahid, Supriyo Ghosh, Thien Duong, and Raymundo Arroyave. "On the Interfacial Phase Growth and Vacancy Evolution during Accelerated Electromigration in Cu/Sn/Cu Microjoints." Acta Materialia (2018), 160, pp. 185-198

Autonomous Efficient Experiment Design for Materials Discovery with Bayesian Model Averaging

Authors: Anjana Talapatra, Shahin Boluki, Thien Duong, Xiaoning Qian, Edward Dougherty, Raymundo Arróyave

Abstract: The accelerated exploration of the materials space in order to identify configurations with optimal properties is an ongoing challenge. Current paradigms are typically centered around the idea of performing this exploration through high-throughput experimentation/computation. Such approaches, however, do not account fo the always present constraints in resources available. Recently, this problem h… ▽ More The accelerated exploration of the materials space in order to identify configurations with optimal properties is an ongoing challenge. Current paradigms are typically centered around the idea of performing this exploration through high-throughput experimentation/computation. Such approaches, however, do not account fo the always present constraints in resources available. Recently, this problem has been addressed by framing materials discovery as an optimal experiment design. This work augments earlier efforts by putting forward a framework that efficiently explores the materials design space not only accounting for resource constraints but also incorporating the notion of model uncertainty. The resulting approach combines Bayesian Model Averaging within Bayesian Optimization in order to realize a system capable of autonomously and adaptively learning not only the most promising regions in the materials space but also the models that most efficiently guide such exploration. The framework is demonstrated by efficiently exploring the MAX ternary carbide/nitride space through Density Functional Theory (DFT) calculations. △ Less

Submitted 30 October, 2018; v1 submitted 14 March, 2018; originally announced March 2018.

Journal ref: Phys. Rev. Materials 2, 113803 (2018)

Exploration of the High Entropy Alloy Space as a Constraint Satisfaction Problem

Authors: Anas Abu-Odeh, Edgar Galvan, Tanner Kirk, Huahai Mao, Qing Chen, Paul Mason, Richard Malak, Raymundo Arroyave

Abstract: High Entropy Alloys (HEAs), Multi-principal Component Alloys (MCA), or Compositionally Complex Alloys (CCAs) are alloys that contain multiple principal alloying elements. While many HEAs have been shown to have unique properties, their discovery has been largely done through costly and time-consuming trial-and-error approaches, with only an infinitesimally small fraction of the entire possible com… ▽ More High Entropy Alloys (HEAs), Multi-principal Component Alloys (MCA), or Compositionally Complex Alloys (CCAs) are alloys that contain multiple principal alloying elements. While many HEAs have been shown to have unique properties, their discovery has been largely done through costly and time-consuming trial-and-error approaches, with only an infinitesimally small fraction of the entire possible composition space having been explored. In this work, the exploration of the HEA composition space is framed as a Continuous Constraint Satisfaction Problem (CCSP) and solved using a novel Constraint Satisfaction Algorithm (CSA) for the rapid and robust exploration of alloy thermodynamic spaces. The algorithm is used to discover regions in the HEA Composition-Temperature space that satisfy desired phase constitution requirements. The algorithm is demonstrated against a new (TCHEA1) CALPHAD HEA thermodynamic database. The database is first validated by comparing phase stability predictions against experiments and then the CSA is deployed and tested against design tasks consisting of identifying not only single phase solid solution regions in ternary, quaternary and quinary composition spaces but also the identification of regions that are likely to yield precipitation-strengthened HEAs. △ Less

Submitted 25 February, 2018; v1 submitted 6 December, 2017; originally announced December 2017.

Comments: 14 pages, 13 figures

On the stochastic phase stability of Ti2AlC-Cr2AlC

Authors: Thien C. Duong, Anjana Talapatra, Woongrak Son, Miladin Radovic, Raymundo Arroyave

Abstract: The quest towards expansion of the MAX design space has been accelerated with the recent discovery of several solid solution and ordered phases involving at least two MAX end members. Going beyond the nominal MAX compounds enables not only fine tuning of existing properties but also entirely new functionality. This search, however, has been mostly done through painstaking experiments as knowledge… ▽ More The quest towards expansion of the MAX design space has been accelerated with the recent discovery of several solid solution and ordered phases involving at least two MAX end members. Going beyond the nominal MAX compounds enables not only fine tuning of existing properties but also entirely new functionality. This search, however, has been mostly done through painstaking experiments as knowledge of the phase stability of the relevant systems is rather scarce. In this work, we report the first attempt to evaluate the finite-temperature pseudo-binary phase diagram of the Ti2AlC-Cr2AlC via first-principles-guided Bayesian CALPHAD framework that accounts for uncertainties not only in ab initio calculations and thermodynamic models but also in synthesis conditions in reported experiments. The phase stability analyses are shown to have good agreement with previous experiments. The work points towards a promising way of investigating phase stability in other MAX Phase systems providing the knowledge necessary to elucidate possible synthesis routes for MAX systems with unprecedented properties. △ Less

Submitted 8 April, 2017; originally announced April 2017.

Comments: 16 pages, 4 figures

Lattice vibrations boost demagnetization entropy in shape memory alloy

Authors: Paul Stonaha, Mike Manley, Nick Bruno, Ibrahim Karaman, Raymundo Arroyave, Navdeep Singh, Douglas Abernathy, Songxue Chi

Abstract: Magnetocaloric (MC) materials present an avenue for chemical-free, solid state refrigeration through cooling via adiabatic demagnetization. We have used inelastic neutron scattering to measure the lattice dynamics in the MC material Ni45Co5Mn36.6In13.4. Upon heating across the Curie Temperature (TC), the material exhibits an anomalous increase in phonon entropy of 0.22 +/- 0.04 kB/atom, which is t… ▽ More Magnetocaloric (MC) materials present an avenue for chemical-free, solid state refrigeration through cooling via adiabatic demagnetization. We have used inelastic neutron scattering to measure the lattice dynamics in the MC material Ni45Co5Mn36.6In13.4. Upon heating across the Curie Temperature (TC), the material exhibits an anomalous increase in phonon entropy of 0.22 +/- 0.04 kB/atom, which is ten times larger than expected from conventional thermal expansion. This transition is accompanied by an abrupt softening of the transverse optic phonon. We present first-principle calculations showing a strong coupling between lattice distortions and magnetic excitations. △ Less

Submitted 12 August, 2015; v1 submitted 31 July, 2015; originally announced July 2015.

First-principles calculation of the instability leading to giant inverse magnetocaloric efects

Authors: Denis Comtesse, Markus E. Gruner, Vladimir V. Sokolovskiy, Vasiliy D. Buchelnikov, Anna Grünebohm, Raymundo Arroyave, Navdeep Singh, Tino Gottschall, Oliver Gutfleisch, Volodymyr Chernenko, Franca Albertini, Sebastian Fähler, Peter Entel

Abstract: The structural and magnetic properties of functional Ni-Mn-Z (Z = Ga, In, Sn) Heusler alloys are studied by first-principles and Monte Carlo methods. The \textit{ab initio} calculations give a basic understanding of the underlying physics which is associated with the strong competition of ferro- and antiferromagnetic interactions with increasing chemical disorder. The resulting $d$-electron orbita… ▽ More The structural and magnetic properties of functional Ni-Mn-Z (Z = Ga, In, Sn) Heusler alloys are studied by first-principles and Monte Carlo methods. The \textit{ab initio} calculations give a basic understanding of the underlying physics which is associated with the strong competition of ferro- and antiferromagnetic interactions with increasing chemical disorder. The resulting $d$-electron orbital dependent magnetic ordering is the driving mechanism of magnetostructural instability which is accompanied by a drop of magnetization governing the size of the magnetocaloric effect. The thermodynamic properties are calculated by using the \textit{ab initio} magnetic exchange coupling constants in finite-temperature Monte Carlo simulations, which are used to accurately reproduce the experimental entropy and adiabatic temperature changes across the magnetostructural transition. △ Less

Submitted 31 January, 2014; originally announced January 2014.

Comments: First-principles calculations of magnetocaloric effects

Thermodynamic modeling of the Hf-Si-O system

Authors: Dongwon Shin, Raymundo Arróyave, Zi-Kui Liu

Abstract: The Hf-O system has been modeled by combining existing experimental data and first-principles calculations results through the CALPHAD approach. Special quasirandom structures of $α$ and $β$ hafnium were generated to calculate the mixing behavior of oxygen and vacancies. For the total energy of oxygen, vibrational, rotational and translational degrees of freedom were considered. The Hf-O system… ▽ More The Hf-O system has been modeled by combining existing experimental data and first-principles calculations results through the CALPHAD approach. Special quasirandom structures of $α$ and $β$ hafnium were generated to calculate the mixing behavior of oxygen and vacancies. For the total energy of oxygen, vibrational, rotational and translational degrees of freedom were considered. The Hf-O system was combined with previously modeled Hf-Si and Si-O systems, and the ternary compound in the Hf-Si-O system, HfSiO$_4$ has been introduced to calculate the stability diagrams pertinent to the thin film processing. △ Less

Submitted 30 August, 2007; originally announced August 2007.

Comments: 13 pages, 11 figures

Journal ref: CALPHAD, 30(4), 375-386 (2006)

Thermodynamic properties of binary HCP solution phases from special quasirandom structures

Authors: Dongwon Shin, Raymundo Arróyave, Zi-Kui Liu, Axel van de Walle

Abstract: Three different special quasirandom structures (SQS) of the substitutional hcp $A_{1-x}B_x$ binary random solutions ($x=0.25$, 0.5, and 0.75) are presented. These structures are able to mimic the most important pair and multi-site correlation functions corresponding to perfectly random hcp solutions at those compositions. Due to the relatively small size of the generated structures, they can be… ▽ More Three different special quasirandom structures (SQS) of the substitutional hcp $A_{1-x}B_x$ binary random solutions ($x=0.25$, 0.5, and 0.75) are presented. These structures are able to mimic the most important pair and multi-site correlation functions corresponding to perfectly random hcp solutions at those compositions. Due to the relatively small size of the generated structures, they can be used to calculate the properties of random hcp alloys via first-principles methods. The structures are relaxed in order to find their lowest energy configurations at each composition. In some cases, it was found that full relaxation resulted in complete loss of their parental symmetry as hcp so geometry optimizations in which no local relaxations are allowed were also performed. In general, the first-principles results for the seven binary systems (Cd-Mg, Mg-Zr, Al-Mg, Mo-Ru, Hf-Ti, Hf-Zr, and Ti-Zr) show good agreement with both formation enthalpy and lattice parameters measurements from experiments. It is concluded that the SQS's presented in this work can be widely used to study the behavior of random hcp solutions. △ Less

Submitted 29 August, 2007; originally announced August 2007.

Comments: 15 pages, 8 figures

Journal ref: Phys. Rev. B 74, 024204 (2006)

Structural and transport properties of epitaxial NaxCoO2 thin films

Authors: A. Venimadhav, A. Soukiassian, D. A. Tenne, Qi Li, X. X. Xi, D. G. Schlom, R. Arroyave, Z. K. Liu, H. P. Sun, Xiaoqing Pan, Minhyea Lee, N. P. Ong

Abstract: We have studied structural and transport properties of epitaxial NaxCoO2 thin films on (0001) sapphire substrate prepared by topotaxially converting an epitaxial Co3O4 film to NaxCoO2 with annealing in Na vapor. The films are c axis oriented and in-plane aligned with [10 1 0] NaxCoO2 rotated by 30 degrees from [10 1 0] sapphire. Different Na vapor pressures during the annealing resulted in films… ▽ More We have studied structural and transport properties of epitaxial NaxCoO2 thin films on (0001) sapphire substrate prepared by topotaxially converting an epitaxial Co3O4 film to NaxCoO2 with annealing in Na vapor. The films are c axis oriented and in-plane aligned with [10 1 0] NaxCoO2 rotated by 30 degrees from [10 1 0] sapphire. Different Na vapor pressures during the annealing resulted in films with different Na concentrations, which showed distinct transport properties. △ Less

Submitted 10 May, 2005; originally announced May 2005.