Excited states of a phosphorus pair in silicon: Combining valley-orbital interaction and electron-electron interactions

Abstract

Excitations of impurity complexes in semiconductors can not only provide a route to fill the terahertz gap in optical technologies, but can also connect local quantum bits to scale up solid-state quantum-computing devices. However, taking into account both the interactions among electrons/holes, and the host band structures, is challenging. Here we combine first-principles band-structure calculations with quantum-chemistry methodology to evaluate the ground and excited states of a pair of phosphorous donors in silicon within s single framework. We use a broken-symmetry Hartree-Fock approach, followed by a time-dependent Hartree-Fock method to compute the excited states. Our Hamiltonian for each valley includes an anisotropic kinetic energy term, which splits the 2p_0 and 2p_+- transitions of isolated donors by ~4 meV, in good agreement with experiments. Our single-valley calculations show the optical response is a strong function of the optical polarisation, and suggest the use of valley polarisation to control optics and reduce oscillations in exchange interactions. When taking into account all valleys, including valley-orbital interactions, we find a gap opens between the 1s to 2p transition and low-energy charge-transfer states within 1s manifolds (which become optically allowed because of inter-donor interactions). In contrast to the single-valley case, we find charge-transfer excited states also in the triplet sector, thanks to the valley degrees of freedom. These states have a qualitatively correct energy as compared with the previous experiments; additionally, we predict new excitations below 20 meV that have not been analysed previously. A statistical average of nearest-neighbour pairs at different separations suggests that THz radiation could be used to excite pairs spin-selectively. Our approach can readily be extended to other donors and to other semiconducting hosts.