Molecular dynamics study of structural changes in berlinite

Abstract

Fractional coordinates and anisotropic temperature factors of atoms in berlinite, AlPO₄ with the quartz topology, were successfully simulated in a molecular dynamics simulation (MDS) at high temperatures. Time-dependent and time-averaged atomic order parameters were analyzed in detail with the aid of spectral densities calculated from trajectory data. These parameters show characteristic behavior of the order-disorder regime for a structure change, where atoms spend most of the time oscillating around the ₁-sites (or ₂-sites) in the low temperature α-phase, but oscillate over both sites in the higher temperature α-phase and the β-phase. In the spectral density functions calculated for atom order parameters, a nearly zero-frequency excitation, which is accompanied by the emergence of large-scale ₁ and ₂ clusters, appears at the Γ point of the Brillouin zone below the transition point T _o, and increases in intensity up to T _o. A low-lying optic branch along Γ-M, which is strongly temperature dependent in the small q-region, is another characteristic of the spectral density functions for the β phase. The spectrum at Γ continuously reduces its frequency from 0.6 THz at temperatures far above T _o to nearly 0 THz at temperatures approaching T _o from above. The dynamical behavior in β berlinite rapidly but continuously changes from that in less oscillatory clusters in the vicinity of T _o to that in the typical β phase at temperatures departing from T _o.