Based on the stress fields, a global optimization method to search the critical slip surface in slopes is proposed by parallel recombination simulated annealing algorithm.