Molecular-dynamics simulations of intergranular films (IGF) containing Si, O, N, and Ca in contact with Si3N4 surfaces containing different levels of interface mixing of the species from the IGF with the crystal surfaces were performed using a multibody interatomic potential. This mixing is equivalent to the formation of a roughened silicon oxynitride crystal surface. With significant interphase mixing at the crystal surfaces, less ordering into the IGF caused by the compositionally modified oxynitride interfaces is observed. Such results are in contrast to our earlier data that showed significant ordering into the IGF induced by the ideally terminated crystal surfaces with no interphase mixing. In all cases, the central position of the first peak in the Si-O pair distribution function (PDF) at the interface ranges from 1.62 to 1.64 A, consistent with recent experimental findings. The central position of the first peak in the Si-N PDF ranges from 1.72 to 1.73 A, consistent with experimental results. With increased interphase mixing, the intensity as well as the area of the first peak of the Si-O and Si-N PDFs for Si attached to the crystal decreases, indicating the decrease of coordination number of O or N with these silicon. Such combined decrease in coordination indicates a significant remnant of vacancies in the crystal surfaces due to the exchange process used here. The results imply a significant effect of interface roughness on the extent of ordering in the amorphous IGF induced by the crystal surface. (C) 2005 American Institute of Physics.