We have developed an order-N method to calculate the total energy and
atomic forces for a nonorthogonal tight-binding Hamiltonian. We have f
irst confirmed the validity of our approach by comparing our results w
ith the exact calculation for a small silicon cluster. The efficiency
of our approach for a large system is demonstrated by determining the
equilibrium configuration of a Si-1000 cluster using quenching and ann
ealing techniques based on molecular dynamics. [S0163-1829(98)02908-7]
.