The results of a Monte Carlo study of oxygenation of 1,2,3 high temperature superconductors was reported. The diffusion of oxygen with locally anisotropic next-nearest-neighbor interactions into an initially empty crystal was considered. The oxygen kinetics obtained had an optimal temperature. The diffusion profile had a sharp front and consisted of two parts. The chain formation drove the diffusion and led to its non-Markovian behavior. The connection between these calculations and experiment was considered.

Monte Carlo Simulation of Oxygen Kinetics in YBa2Cu3O6+x. V.V.Alekseenko, I.I.Titov: Physica C, 1991, 185-189[3], 1715-6