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 multipy 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.