El potencial zeta es una propiedad fisicoquímica de gran relevancia en el estudio de nanopartículas y en la industria del petróleo, principalmente por su influencia en la determinación de la estabilidad de sistemas coloidales. En este trabajo de investigación se implementó una metodología para calcular el potencial zeta mediante simulaciones de dinámica molecular de no equilibrio (NEMD). Se planteó como modelo de prueba un sistema de superficie de rutilo (110) que interactúa con una solución acuosa de NaCl, para posteriormente verificar la transferibilidad del método en superficies de sílice con ambientes Q2 a pH 2 y pH 5. Los resultados obtenidos muestran que la metodología propuesta permite calcular valores de potencial zeta acordes a las superficies evaluadas y muestra la capacidad de poder ser usada en diferentes sistemas cristalinos.