This paper addresses a novel optimization technique based on molecular dynamics simulation which has been utilized for physical model simulation at various disciplines. In this study, objective functions are considered to be potential functions, which depict molecular interactions. Comparisons of typical optimization method such as the steepest descent and the present method for several test functions are made. The present method shows applicability and stability in finding a global optimum.