The electronic and geometrical structures, in particular the electric field gradients (EFGs), of ͓CdD͔ 2 (D P, As, Sb) acceptor-donor pairs in Si and Ge are studied using the full potential Korringa-KohnRostoker Green's function method. In addition, also neutral complexes (͓CdD͔ 0 ) and trimers (͓CdD 2 ͔ 0 ) in Si are investigated. The EFG depends very sensitively on the large lattice relaxations induced by the defects and can be understood by a simple hybridization model. Our calculations are in good agreement with experimental results and provide a consistent picture of acceptor-donor complexes in Si and Ge. PACS numbers: 71.55.Cn, 61.72.Tt, The donor-acceptor pairs in silicon and germanium have attracted intensive research interest for a long time [1][2][3][4][5]. If impurities with opposite charges are either intentionally or unintentionally introduced in semiconductors, the Coulomb attraction leads to the formation of donoracceptor complexes which can cause serious technological problems such as an uncontrolled annihilation or creation of charge carriers. To study the behavior of these defects, mostly nuclear techniques such as perturbed gg angular correlation (PAC) experiments were successfully applied to obtain microscopic information on an atomic scale. A very popular and well-known probe atom is the radioactive 111 In͞ 111 Cd nucleus. In these experiments initially a dimer complex is formed consisting of the In acceptor and a donor atom such as P, As, or Sb, so that neutral ͓InD͔ 0 (D P, As, Sb) pairs are produced. Since such neutral pairs have no states in the gap and thus are not magnetically or electrically active, very little information can be obtained with other experimental methods. After the decay 111 In ! 111 Cd, the PAC measurement takes place at the daughter complex CdD. The fact that the electric field gradient is independent from the InD mother complex has been proven by PAC experiments starting from the 111m Cd nucleus as mother probe atom [3]. Because of the pair formation the cubic surrounding of the probe atom 111 Cd is disturbed, leading to an electric field gradient (EFG) at the Cd site which can be determined in particular for nonmagnetic systems. This "fingerprint" of the defect system gives some important information about the defect structure. For instance, the main axis of the EFG is found to point along the [111] direction (h 0), which strongly suggests a nearest neighbor (nn) pair configuration. However, for a fundamental understanding of the complexes electronic structure calculations are needed [3], which can be performed with good precision by density-functional methods [6,7]. The use of pseudopotential methods, which usually give excellent results in semiconductor physics, is not suitable in this case since the EFG is determined by the p charge density in the inner core region [6], so that up to now no theoretical treatment of these fundamental acceptor-donor complexes in Si and Ge exists.In this paper we calculate the electronic and geometrical structures and the EFGs of s...