Simulations of cubic-tetragonal ferroelastics

We study domain patterns in cubic-tetragonal ferroelastics by solving numerically equations of motion derived from a Landau model of the phase transition, including dissipative stresses. Our system sizes, of up to 256(3) points, are large enough to reveal many structures observed experimentally. Most patterns found at late stages in the relaxation are multiply banded; all three tetragonal variants appear, but inequivalently. Two of the variants form broad primary bands; the third intrudes into the others to form narrow secondary bands with the hosts. On colliding with walls between the primary variants, the third either terminates or forms a chevron. The multiply banded patterns, with the two domain sizes, the chevrons and the terminations, are seen in the microscopy of zirconia and other cubic-tetragonal ferroelastics. We examine also transient structures obtained much earlier in the relaxation; these show the above features and others also observed in experiment.