MONTE-CARLO STOCHASTIC-DYNAMICS STUDY OF DIELECTRIC RESPONSE AND NONERGODICITY IN PROTON GLASS

Dynamics of the O-H ... O bond proton glass of the type M(1-x)(NW4)(x) W(2)AO(4) (M=Rb or K, W=H or D, A=P or As) has been simulated using th e Monte Carlo stochastic-dynamics method that allows one to simulate r eal time dynamics. The simulation is based both on microscopic interac tions of protons and on interaction with an external static electric f ield. The polarization decay and response to step field has been compa red with the Kohlrausch-Williams-Watts stretched exponential form and with predictions of a microscopic ''bound charge carrier'' model. Stud ying the proton dynamics by a field cooling simulation has revealed no nergodic behavior at low temperatures.