A new modification of Metropolis algorithm, which allows exploring thermodynamic species of seldom thermal fluctuations of small atomic clusters, has been proposed. Based on this algorithm, a new hybrid molecular dynamics with Monte-Carlo method was combined. It allows to simulate thermal decay of small atomic clusters at low temperatures, where conventional molecular dynamics method is unavailable.