Combining classical molecular dynamics and first-principles DFT calculations, we perfom an extensive investigation of low energy configurations for He n V m complexes in silicon. The optimal helium fillings are hence determined for V 1 , V 2 , and V 6 (figure on the right), and the structures formed by helium atoms arrangements in the vacancy defect are analyzed. For V 1 and V 2 , the He atoms structure is mainly controled by the host silicon matrix, whereas a high density helium packing is obtained for V 6 . For the latter, we estimate a helium density of about 170 He nm −3 in the center of the hexa-vacancy at the optimal helium filling.Relaxed structures obtained from DFT calculations for configurations with the lowest formation energies: (a) He 14 V 1 , (b) He 20 V 2 , and (c) He 40 V 6 .