No Arabic abstract
Formation energies of charged point defects in semiconductors are calculated using periodic supercells, which entail a divergence arising from long-range Coulombic interactions. The divergence is typically removed by the so-called jellium approach. Recently, Wu, Zhang and Pantelides [WZP, Phys. Rev. Lett. 119, 105501 (2017)] traced the origin of the divergence to the assumption that charged defects are formed by physically removing electrons from or adding electrons to the crystal, violating charge neutrality, a key principle of statistical mechanics that determines the Fermi level. An alternative theory was constructed by recognizing that charged defects form by trading carriers with the energy bands, whereby supercells are always charge-neutral so that no divergence is present and no ad-hoc procedures need to be adopted for calculations. Here we give a more detailed exposition of the foundations of both methods and show that the jellium approach can be derived from the statistical-mechanics-backed WZP definition by steps whose validity cannot be assessed a priori. In particular, the divergence appears when the charge density of band carriers is dropped, leaving a supercharged crystal. In the case of charged defects in two-dimensional (2D) materials, unphysical fields appear in vacuum regions. None of these pathological features are present in the reformulated theory. Finally, we report new calculations in both bulk and 2D materials. The WZP approach yields formation energies that differ from jellium values by up to ~1 eV. By analyzing the spatial distribution of wave functions and defect potentials, we provide insights into the inner workings of both methods and demonstrate that the failure of the jellium approach to include the neutralizing electron density of band carriers, as is the case in the physical system, is responsible for the numerical differences between the two methods.
Two-dimensional (2D) materials are strongly affected by the dielectric environment including substrates, making it an important factor in designing materials for quantum and electronic technologies. Yet, first-principles evaluation of charged defect energetics in 2D materials typically do not include substrates due to the high computational cost. We present a general continuum model approach to incorporate substrate effects directly in density-functional theory calculations of charged defects in the 2D material alone. We show that this technique accurately predicts charge defect energies compared to much more expensive explicit substrate calculations, but with the computational expediency of calculating defects in free-standing 2D materials. Using this technique, we rapidly predict the substantial modification of charge transition levels of two defects in MoS$_2$ and ten defects promising for quantum technologies in hBN, due to SiO$_2$ and diamond substrates. This establishes a foundation for high-throughput computational screening of new quantum defects in 2D materials that critically accounts for substrate effects.
Defects in 2D materials are becoming prominent candidates for quantum emitters and scalable optoelectronic applications. However, several physical properties that characterize their behavior, such as charged defect ionization energies, are difficult to simulate with conventional first-principles methods, mainly because of the weak and anisotropic dielectric screening caused by the reduced dimensionality. We establish fundamental principles for accurate and efficient calculations of charged defect ionization energies and electronic structure in ultrathin 2D materials. We propose to use the vacuum level as the reference for defect charge transition levels (CTLs) because it gives robust results insensitive to the level of theory, unlike commonly used band edge positions. Furthermore, we determine the fraction of Fock exchange in hybrid functionals for accurate band gaps and band edge positions of 2D materials by enforcing the generalized Koopmans condition of localized defect states. We found the obtained fractions of Fock exchange vary significantly from 0.2 for bulk $h$-BN to 0.4 for monolayer $h$-BN, whose band gaps are also in good agreement with experimental results and calculated GW results. The combination of these methods allows for reliable and efficient prediction of defect ionization energies (difference between CTLs and band edge positions). We motivate and generalize these findings with several examples including different defects in monolayer to few-layer hexagonal boron nitride ($h$-BN), monolayer MoS$_2$ and graphane. Finally, we show that increasing the number of layers of $h$-BN systematically lowers defect ionization energies, mainly through CTLs shifting towards vacuum, with conduction band minima kept almost unchanged.
Charged defects in 2D materials have emerging applications in quantum technologies such as quantum emitters and quantum computation. Advancement of these technologies requires rational design of ideal defect centers, demanding reliable computation methods for quantitatively accurate prediction of defect properties. We present an accurate, parameter-free and efficient procedure to evaluate quasiparticle defect states and thermodynamic charge transition levels of defects in 2D materials. Importantly, we solve critical issues that stem from the strongly anisotropic screening in 2D materials, that have so far precluded accurate prediction of charge transition levels in these materials. Using this procedure, we investigate various defects in monolayer hexagonal boron nitride (h-BN) for their charge transition levels, stable spin states and optical excitations. We identify $C_BN_V$ (nitrogen vacancy adjacent to carbon substitution of boron) to be the most promising defect candidate for scalable quantum bit and emitter applications.
A detailed understanding of charged defects in two-dimensional semiconductors is needed for the development of ultrathin electronic devices. Here, we study negatively charged acceptor impurities in monolayer WS$_2$ using a combination of scanning tunnelling spectroscopy and large-scale atomistic electronic structure calculations. We observe several localized defect states of hydrogenic wave function character in the vicinity of the valence band edge. Some of these defect states are bound, while others are resonant. The resonant states result from the multi-valley valence band structure of WS$_2$, whereby localized states originating from the secondary valence band maximum at $Gamma$ hybridize with continuum states from the primary valence band maximum at K/K$^{prime}$. Resonant states have important consequences for electron transport as they can trap mobile carriers for several tens of picoseconds.
We provide a quantitative analysis of all kinds of topological defects present in 2D passive and active repulsive disk systems. We show that the passage from the solid to the hexatic is driven by the unbinding of dislocations. Instead, although we see dissociation of disclinations as soon as the liquid phase appears, extended clusters of defects largely dominate below the solid-hexatic critical line. The latter percolate at the hexatic-liquid transition in continuous cases or within the coexistence region in discontinuous ones, and their form gets more ramified for increasing activity.