Este artículo presenta un modelo bimodal basado en la distribución epsilon-skew-normal. Esta extensión genera distribuciones bimodales similares a las producidas por la mezcla de distribuciones normales. Estudiamos las propiedades básicas de esta nueva familia. Aplicamos estimadores de máxima verosimilitud, calculamos la matriz de información y presentamos un estudio de simulación para evaluar la recuperación de parámetros. Finalmente, ilustramos los resultados en tres conjuntos de datos reales, sugiriendo esta nueva distribución como una alternativa plausible para modelar datos bimodales.