Magnetization plateaus in the antiferromagnetic Ising chain with single-ion anisotropy and quenched disorder

We have studied the presence of plateaus on the low-temperature magnetization of an antiferromagnetic spin-1 chain, as an external uniform magnetic field is varied. A crystal-field interaction is present in the model and the exchange constants follow a random quenched (Bernoulli or Gaussian) distribution. Using a transfer-matrix technique we calculate the largest Lyapunov exponent and, from it, the magnetization at low temperatures as a function of the magnetic field, for different values of the crystal field and the width of the distributions. For the Bernoulli distribution, the number of plateaus increases, with respect to the uniform case [Litaiff et al., Solid State Commun. 147, 494 (2008)] and their presence can be linked to different ground states, when the magnetic field is varied. For the Gaussian distributions, the uniform scenario is maintained, for small widths, but the plateaus structure disappears as the width increases.