The modeling of the atom-by-atom growth of extended defects is coupled to the diffusion equations of boron by transferring the free interstitial supersaturation calculated with a defect model into a process simulator. Two methods to achieve this coupling (equilibrium method and fully coupled method, respectively) are presented and tested against a variety of experimental conditions. They are first applied to a transient enhanced diffusion experiment carried out on a structure containing several B delta-doped layers, in which the amount of diffusion of the different layers is accurately predicted. The fully coupled method is then used to simulate the diffusion of ultrashallow B-implanted profiles. This work definitely demonstrates the relevance of accurate physical defect models for the successful design of ultrashallow junctions in future generations of integrated circuits.