Abstract
Lattice energy minimisation calculations with empirical shell model potentials have been performed on feldspar solid solutions as a function of composition and tetrahedral Al/Si order. The monoclinic-triclinic phase transitions have been simulated across the KxNa1-xAlSi3O8 and SrxCa1-xA12Si2O8 solid solutions. In both cases it is found that the transitions are driven by an elastic softening without critical softening of an optic phonon. The results for the ferroelastic spontaneous strains are consistent with experimental observation, and explain the nature of the coupling between the ferroelastic instability and Al/Si ordering.