Simulacion Monte Carlo usando el campo ´ de fuerza polarizable amoeba un programa en fortran90

Un simple algoritmo de simulación Monte Carlo basado en este campo de fuerzas fue implementado. Este algoritmo consiste en 6 pasos: generar un esqueleto de carbonó de manera aleatoria; agregar los ´átomos de hidrógeno; convertir la matriz Z a las coordenadas cartesianas con parámetros apropiados; evaluar la energía del campo de fuerzas MMF; probar la condición de Metropolis. Para probar este algoritmo, las simulaciones de tres geometrías iniciales de la molécula de pentano fueron realizadas. Para cada geometría se evaluaron cinco tipos de energía sumando sus contribuciones obteniendo así la energía total inicial de cada una de estas, posteriormente se hizo un cambio aleatorio al ´ángulo diedro modificando su geometría, con esta nueva geometría tenemos una segunda energía total y podemos aplicar el método de Monte Carlo y el algoritmo de metrópolis usando la condición de Boltzmann para encontrar la energía mínima. Se presenta una breve discusión de la implementación y los resultados de las simulaciones.