One of the most important design issues for phase change materials is to engineer the phase transition process. The challenge of accurately predicting a phase transition is estimating the true value of transition rate, which is determined by the saddle point with the minimum energy barrier between stable states on the potential energy surface (PES). In this thesis, a new algorithm for searching the minimum energy path (MEP) is presented. The new algorithm is able to locate both the saddle point and local minima simultaneously. Therefore no prior knowledge of the precise positions for the reactant and product on the PES is needed. Unlike existing pathway search methods, the algorithm is able to search multiple transition paths on the PES simultaneously, which gives us a more comprehensive view of the energy landscape than searching individual ones. In this method, a Bézier curve is used to represent each transition path. During the searching process, the reactant and product states are located by minimizing the two end control points of the curve, while the shape of the transition pathway is refined by moving the intermediate control points of the curve in the conjugate directions. A curve subdivision scheme is developed so that multiple transitions paths can be located. The algorithm is demonstrated by examples of LEPS potential, LEPS plus harmonic oscillator potential, and PESs defined by Rastrigin function and Schwefel function.