Proponemos el desarrollo teórico e implementación computacional de una metodología que nos permita estudiar teóricamente el efecto en los puentes de hidrógeno de conglomerados de agua (H2O_n, n=2,5) de sustituciones isotópicas graduales de núcleos de protio por deuterio y tritio. El modelo teórico que se propone desarrollar es una extensión del método de orbitales moleculares nucleares y electrónicos (OMNE) para mezclas de núcleos cuanticos. El método propuesto está basado en la teoría de perturbaciones Møller-Plesset de segundo orden (MP2), esperamos que este nos provea una mejor descripción de los puentes de hidrogeno que su contraparte Hartree-Fock al incluir parte de la correlación electrónica, nuclear y núcleo -electronica. La implementación computacional se realizará en FORTRAN y C++, siguiendo una programación orientada a objetos. Se espera que el estudio permita observar el efecto de la sustitución isotópica parcial sobre las energias de enlace, las geometrías moleculares, las energía de punto cero y los niveles vibracionales de los diferentes complejos. Los resultados de esta investigacion serán claves en futuras investigaciones sobre el tunelamiento nuclear y procesos dinámicos en conglomerados de agua y especies solvatadas.