A Monte Carlo simulation of two-dimensional patch-domain formation and domain coarsening in antiferromagnetically coupled compensated multilayers with fourfold in-plane anisotropy is presented. The simulation accounts for both the emergence of small patch domains on unsaturation and domain coarsening on spin dispartial in hard direction. The simulated domain patterns are in good agreement with published Kerr microscopic images.1. Introduction Antiferromagnetically (AF) coupled metallic multilayers (ML) show giant magnetoresistance (GMR) effect [1]. Recently a domain tailoring mechanism was found on an AF-coupled Fe/Cr ML with fourfold in-plane anisotropy by off-specular synchrotron Mö ssbauer reflectometry (SMR) and verified by polarised neutron reflectometry (PNR) [2]. Small (mm size) domains form when the external magnetic field H ext is decreased from the saturation to remanence. Most of these domains transform into at least one order of magnitude bigger ones (domain coarsening) on the spin flop [3][4][5]. A simple phenomenological model [2] and Monte Carlo simulations [6] of domain patterns were presented to reproduce the domain nucleation on 'unsaturation', i.e., on releasing the saturation magnetic field applied in an easy direction. The same algorithm with similar 'flipping rules' was applied to simulate the domain coarsening observed on spin flop [7].In a strongly AF-coupled metallic ML of fourfold crystalline anisotropy the domain structure of the individual ferromagnetic layers is correlated through the ML stack from the substrate to the surface allowing for a two-dimensional representation of the AF domains. Easy-axis and hard-axis scenarios are qualitatively different. In our model [2] we associate the domain formation and coarsening with the effective correlation length of the saturation field H sat and that of the spin-flop field H sf , respectively. Unsaturating the sample along an easy axis down to remanence leads to a single stable configuration with the magnetisations perpendicular to the field. Magnetic field and anisotropy act in the same direction. The effective correlation length of H sat will determine the average 'primary domain' [2] size in remanence, the anisotropy playing a minor role in the formation. On passing the spin flop in an easy-axis scenario, however, the much larger effective correlation length of the spin-flop field will dominantly determine the resulting 'secondary-domain' size and a remarkable domain coarsening effect is observed [2,6,7]. In a hard-direction-unsaturation scenario, however, these two phases of