Ab initio atomic and electronic structures of neutral ͑V 0 2 ͒ and charged ͑V 1 2 , V 2 2 , V 22 2 ͒ Si divacancies are investigated using bulk-terminated clusters with up to ഠ320 Si atoms. For the first time, the relaxed structures for V 1 2 , V 0 2 , and V 2 2 are found to exhibit large pairing Jahn-Teller distortions consistent with electron paramagnetic resonance experiments. Atomic relaxations, Jahn-Teller and relaxation energies, and hyperfine parameters are calculated for various charge states and compared with available experimental data. The cluster size and the anisotropic nature of the relaxations play key roles in establishing the stability of the structures with large pairing distortions. PACS numbers: 61.72.Ji, 71.55.CnSince the pioneering electron paramagnetic resonance (EPR) studies of Watkins and Corbett [1], there have been controversies about the electronic and atomic structures of the divacancy V 2 in crystalline Si. These controversies have revolved around the exact nature of the symmetrylowering Jahn-Teller (JT) distortions that split the degenerate deep levels. To date, the sense and the magnitude of these distortions, as inferred from the EPR data and theoretical calculations, have been at variance. This disagreement has created an ongoing debate for the ground state structures of V 2 in Si and presented a fundamental challenge to understanding this important semiconductor defect.The ideal simple divacancy, which is created by removing two neighboring Si atoms, has a D 3d symmetry and two doubly degenerate deep levels, labeled e u and e g (Fig. 1). While the e g level is empty, the lower e u level is occupied by one, two, and three electrons for the positively charged ͑V 1 2 ͒, neutral ͑V 0 2 ͒, and negatively charged ͑V 2 2 ͒ divacancies, respectively. As a result of this degeneracy, the lattice distorts to gain electronic energy. Upon symmetry lowering to C 2h , each e level splits into levels labeled a and b. Of these four levels, a g and b u have amplitudes, while a u and b g have nodes on the mirror plane of the C 2h distortion that contains the vacant sites and atoms 3 and 6 in Fig. 1(a). The ordering of these four levels inside the band gap depends on the particular sense and magnitude of the distortion.The EPR experiments of Watkins and Corbett indicated that V 1 2 and V 2 2 have the low symmetry of C 2h , and the highest occupied levels for both charge states have amplitudes on the mirror plane. They also established that the JT distortions occur in a pairing sense, which results in d 12 , d 13 d 23 , where d ij is the distance between atoms i and j [ Fig. 1(b)]. The JT energies E JT were estimated to be 1.3 eV and 2.4 eV for V 1 2 and V 2 2 , respectively. From these experimental findings, Watkins and Corbett concluded that the JT distortion had to be large enough for the a g level split from the upper e g level to be lowered below the a u level, as shown in Fig. 1(b). This would give configurations of V 1 2 :b 1 u (or a 1 g ) and V 2 2 :b 2 u a 1 g (or a 2 g b 1 u ), explaining ...