Analytical gradient-based optimization of CALPHAD model parameters

The calibration of CALPHAD (CALculation of PHAse Diagrams) models involves the solution of a very challenging high-dimensional multiobjective optimization problem. Traditional approaches to parameter fitting predominantly rely on gradient-free methods, which while robust, are computationally inefficient and often scale poorly with model complexity. In this work, we introduce and demonstrate a generalizable framework for analytic gradient-based optimization of the parameters of the CALPHAD model enabled by the recently formalized Jansson derivative technique. This method allows for efficient evaluation of gradients of thermodynamic properties at equilibrium with respect to model parameters, even in the presence of arbitrarily complex internal degrees of freedom. Leveraging these semi-analytic gradients, we employ the conjugate gradient (CG) method to optimize thermodynamic model parameters for four binary alloy systems: Cu-Mg, Fe-Ni, Cr-Ni, and Cr-Fe. Across all systems, CG achieves comparable or superior optimality relative to Bayesian ensemble Markov Chain Monte Carlo (MCMC) with improvements in computational efficiency ranging from one to three orders of magnitude. Furthermore, our results establish a new paradigm for CALPHAD assessments in which high fidelity data-rich model calibration becomes tractable using deterministic gradient-informed algorithms.. Authors: Kunselman, Courtney [Texas A&M University, College Station, TX (United States)] (ORCID:000000034903874X); Bocklund, Brandon [Lawrence Livermore National Laboratory (LLNL), Livermore, CA (United States)] (ORCID:0000000232648413); Otis, Richard [Proteus Space Inc., Los Angeles, CA (United States)] (ORCID:0000000211479032); Arróyave, Raymundo [Texas A&M University, College Station, TX (United States)]. DOE Contract: AC52-07NA27344. Subjects: CALPHAD; Equilibrium calculations; Thermodynamics