Lattice instability during phase transformations under multiaxial stress: Modified transformation work criterion

A continuum/atomistic approach for predicting lattice instability during crystal-crystal phase transformations (PTs) is developed for the general loading with an arbitrary stress tensor and large strains. It is based on a synergistic combination of the generalized Landau-type theory for PTs and molecular dynamics (MD) simulations. The continuum approach describes the entire dissipative transformation process in terms of an order parameter, and the general form of the instability criterion is derived utilizing the second law of thermodynamics. The feedback from MD allowed us to present the instability criterion for both direct and reverse PTs in terms of the critical value of the modified transformation work, which is linear in components of the true stress tensor. It was calibrated by MD simulations for direct and reverse PTs between semiconducting silicon Si I and metallic Si II phases under just two different stress states. Then, it describes hundreds of MD simulations under various combinations of three normal and three shear stresses. In particular, the atomistic simulations show that the effects of all three shear stresses along cubic axes on lattice instability of Si I are negligible, which is in agreement with our criterion.