Brownian dynamics simulations of the filming process of a mixed polymer latex in the water evaporation stage were performed in order to explore the effect of surface potential on latex particle packing and distribution at a temperature far below the glass transitions of polymers in bulk. Polymer latex particles are modeled as spheres that interact via DLVO potential with various surface charge densities for emulsifier-free emulsion polymerized particles and dispersion polymerized particles. It is found that the distribution of modeled poly(methyl methacrylate) and polystyrene latex particles in the finally formed film exhibits a noticeable dependence of surface potentials of latex particles. When the difference of the surface potentials between binary mixed latex particles is small, the particles distribute randomly. In contrast, when the difference of the surface potentials between binary mixed latex particles is large, heterocoagulation occurs and the polymer latex in which the repulsive electrostatic potential is weak will form clusters in the film. The results are in agreement with laser confocal fluorescence microscopy observations of fluorescent dye labeled poly(methyl methacrylate) and polystyrene mixed latex films. The correlation between latex particles increases with increasing repulsive electrostatic potential, and the spatial order can be obtained at the end of the water evaporation stage. Copyright 2000 Academic Press.