New technologies such as cluster or molecular ion implantation and flash lamp annealing (FLA) seem to be applied as advanced shallow junction technologies. In this paper, the Molecular Dynamics (MD) model based simulation on BB 10 H 14 and B 18 H 22 implantation is performed. The spike annealing of cluster implantation is simulated by atomistic model. The inactivation and clustering of B implanted at 0.5keV and annealed at 900 C~1200 C are correctly simulated by atomistic model. The simulation on activation ratio of B in FLA is presented. The discussion on cluster evolution in annealing is performed. o o 1. Introduction Boron cluster ion implantation is a potential technology for shallow junction formation in integrated circuits manufacture [1][2]. With the technology becoming more practicable, capability of accurately simulating it for shallow junction formation becomes necessary. As far as we know, by the full MD method, it is still difficult to simulate the implantation at high dose, e.g. 10 15 B/cm 2 , because of the limitation of computation capability. Kinetic Monte Carlo (KMC) method simulation has been successfully applied in the simulation of annealing of B monomer implantation. It is still to be extended to the annealing of B cluster implantation. FLA technology can significantly suppress the enhanced diffusion of B in annealing [2]. On the other hand, the activation ratio of B after FLA is lower than rapid thermal annealing. The simulation on annealing should capture this character. In this paper, the localized MD method [3] is applied to simulate boron cluster implantation. Attractive interaction within B cluster is considered. The simulated concentration profiles of B and H agree with SIMS data. The difference of B monomer and cluster implantation can be reproduced well. It is demonstrated that the KMC method can be extended to the simulation on annealing of B cluster implantation. The activation of B implanted at 0.5keV and annealed at 900 o C~1200 o C as well as FLA can be simulated by KMC model.