Track papers across labs, topics, and trends.
Search papers, labs, and topics across Lattice.
To overcome severe convergence bottlenecks in computing non-dilute multicomponent alloy diffusion coefficients, the authors formulate two analytical frameworks that eliminate the need to sample rare vacancy composition fluctuations in Monte Carlo simulations. Across binary to quinary Hf–Mo–Nb–Ti–Zr alloys, a single-vacancy truncation of the semi-grand canonical partition function reproduces full diffusion matrix eigenspectra with a 90% reduction in sampling effort. The formulation analytically proves that the dominant eigenvalue of the diffusion tensor is universally the vacancy tracer diffusion coefficient, establishing an exact scaling relation for non-equilibrium transport.
Sampling rare vacancy fluctuations is completely unnecessary for multicomponent alloy diffusion: a single-vacancy partition function truncation slashes Monte Carlo compute by 90% while analytically isolating the true vacancy tracer diffusion coefficient.
Computing non-dilute diffusion coefficients from Onsager transport coefficients of a multicomponent alloy requires the thermodynamic factor, which measures the curvature of the free energy with respect to composition. As diffusion in metals is mediated by vacancies, the curvature that matters is not that of the vacancy-free alloy, which is straightforward to compute, but that of an alloy carrying a dilute concentration of vacancies. The resulting matrix is nearly singular, and composition fluctuations from which it is obtained converge slowly in Monte Carlo simulations. Here we develop two routes that avoid sampling vacancy composition fluctuations. The first constructs the thermodynamic factor from two quantities, the free-energy curvature of the vacancy-free alloy and the vacancy concentration. The second truncates the semi-grand canonical partition function at a single vacancy and recovers the full thermodynamic factor of the vacancy-containing alloy. The same separation shows that the largest eigenvalue of the diffusion matrix is the vacancy tracer diffusion coefficient for any number of components and any degree of non-ideality, while the remaining eigenvalues scale with the vacancy concentration. Across binary through quinary alloys of the Hf--Mo--Nb--Ti--Zr system, both routes reproduce eigenvalues and eigenvectors of diffusion matrices obtained from fully converged Monte Carlo simulations. The single-vacancy expansion reaches this accuracy with roughly one tenth of the sampling effort. A scaling relation extends both routes to local vacancy chemical potentials away from equilibrium, without additional simulations. These results allow thermodynamic descriptions from atomistic models or CALPHAD assessments to be used directly in mesoscale simulations of mass transport.