exaPD: A highly parallelizable workflow for multi-element phase diagram (PD) construction

Phase diagrams (PDs) illustrate the relative stability of competing phases under varying conditions, serving as critical tools for synthesizing complex materials. Reliable phase diagrams rely on precise free energy calculations, which are computationally intensive. We introduce exaPD, a user-friendly workflow that enables simultaneous sampling of multiple phases across a fine mesh of temperature and composition for free energy calculations. The package employs standard molecular dynamics (MD) and Monte Carlo (MC) sampling techniques, as implemented in the LAMMPS package. Various interatomic potentials are supported, including the neural network potentials with near {\it ab initio} accuracy. A global controller, built with Parsl, manages the MD/MC jobs to achieve massive parallelization with near ideal scalability. The resulting free energies of both liquid and solid phases, including solid solutions, are integrated into CALPHAD modeling using the PYCALPHAD package for constructing the phase diagram.