Order/disorder phase transition in cordierite and its possible relationship to the development of symplectite reaction textures in granulites

Based on a consistent set of empirical interatomic potentials, static structure energy calculations of various Al/Si configurations in the supercell of Mg-cordierite and Monte Carlo simulations the phase transition between the orthorhombic and hexagonal modifications of cordierite (Crd) is predicted at 1623 K. The temperature dependences of the enthalpy, entropy, and free energy of the Al/Si disorder were calculated using the method of thermodynamic integration. The simulations suggest that the commonly observed crystallization of cordierite in the disordered hexagonal form could be related to a tendency of Al to occupy T1 site, which is driven by local charge balance. The increase in the Al fraction in the T1 site over the ratio of 2/3(T1): 1/3(T2), that characterizes the ordered state, precludes formation of the domains of the orthorhombic phase. This intrinsic tendency to the crystallization of the metastable hexagonal phase could have significantly postponed the formation of the association of orthorhombic cordierite and orthopyroxene over the association of quartz and garnet in metapelites subjected to granulite facies metamorphism. The textures of local metasomatic replacement (the formation of Crd + Opx or Spr + Crd symplectites between the grains of garnet and quartz) indicate the thermodynamic instability of the association of Qtz + Grt at the moment of the metasomatic reaction. This instability could have been caused by the difficulty of equilibrium nucleation of orthorhombic cordierite