A Convergent ADMM Algorithm for Grain Boundary Energy Minimization

TL;DR

This paper introduces a modified ADMM algorithm with convergence guarantees for a nonconvex constrained minimization problem in materials science, specifically for modeling grain boundary energy minimization with dislocation structures.

Contribution

The paper develops a convergent ADMM variant tailored for nonconvex problems with linear constraints, extending existing convergence analysis to new settings.

Findings

The modified ADMM converges to a stationary point under certain conditions.

Numerical experiments validate the convergence and effectiveness of the proposed method.

The approach outperforms traditional penalty and augmented Lagrangian methods in the tested scenarios.

Abstract

In this paper, we study a constrained minimization problem that arise from materials science to determine the dislocation (line defect) structure of grain boundaries. The problems aims to minimize the energy of the grain boundary with dislocation structure subject to the constraint of Frank's formula. In this constrained minimization problem, the objective function, i.e., the grain boundary energy, is nonconvex and separable, and the constraints are linear. To solve this constrained minimization problem, we modify the alternating direction method of multipliers (ADMM) with an increasing penalty parameter. We provide a convergence analysis of the modified ADMM in this nonconvex minimization problem, with settings not considered by the existing ADMM convergence studies. Specifically, in the linear constraints, the coefficient matrix of each subvariable block is of full column rank. This property makes each subvariable minimization strongly convex if the penalty parameter is large enough, and contributes to the convergence of ADMM without any convex assumption on the entire objective function. We prove that the limit of the sequence from the modified ADMM is primal feasible and is the stationary point of the augmented Lagrangian function. Furthermore, we obtain sufficient conditions to show that the objective function is quasi-convex and thus it has a unique minimum over the given domain. Numerical examples are presented to validate the convergence of the algorithm, and results of the penalty method, the augmented Lagrangian method, and the modified ADMM are compared.