m_epsilon(t) = -2(m_epsilon(t) + tanh(-epsilon β rho_epsilon(t)))