Gap modification of atomically thin boron nitride by phonon mediated interactions

A theory is presented for the modification of bandgaps in atomically thin boron nitride (BN) by attractive interactions mediated through phonons in a polarizable substrate, or in the BN plane. Gap equations are solved, and gap enhancements are found to range up to 70% for dimensionless electron-phonon coupling λ =1, indicating that a proportion of the measured BN bandgap may have a phonon origin.

dominant cause of a gap of order 2 . Tight-binding fits to results from ab-initio simulations of monolayer BN have established the hopping t=2.33 eV [10], with an estimate of = 1.96 eV=0.84 t. The experiments indicate larger gaps: bulk h-BN has 5.971 eV [11], and monolayer h-BN has a gap of 5.56 eV [12] corresponding to = 2.78 eV=1.20 t. There is significant variation in phonon energies, q , in h-BN [13]. Longitudinal acoustic (LA) phonon energies range up to around 140 meV at the M point, and transverse acoustic (TA) modes to around 110 meV at the K point. Optical phonon energies range between 160 and 200 meV. Coupling, f n (m), between electrons and phonons in either a polarisable substrate or the BN monolayer, is possible and the corresponding Hamiltonian is as follows: The Hamiltonian terms are shown schematically in Figure 1 (a). a † nσ creates electrons of spin σ on B sites and c † n σ on N sites. Vectors n are to atoms in the monolayer, and m to atoms in the substrate. N m and n nσ are the number operators for phonons and electrons respectively. ξ m is the atomic displacement. The Hamiltonian is also approximately valid for interactions in the plane, and Figure 1  on B sites and vice versa). The diagram indicates that the strongest interaction is between electrons and optical phonon modes.
For simplicity, the Holstein electron-phonon interaction was used, f n (m) ∝ δ n,n , which qualitatively captures the physics. There may be quantitative changes to the results for longer range Fröhlich interactions and from modulation of the electron-phonon interaction due to incommensurability of the substrate, which was estimated at around ±8% of the average value [8].

Results and discussion
The low order perturbation theory is applicable for low phonon frequency and weak coupling. A set of gap equations was derived by symmetrizing the self energy, The local approximation used here is a good starting point because the modulated potential is large, and electrons are well localized. Off-diagonal terms do not feature in the lowest order perturbation theory for the Holstein model since the interaction is site diagonal. Z n is the quasi-particle weight and¯ n is the gap function. For bosonic quantities, ω s = 2π k B Ts and for fermions, ω n = 2π k B T(n + 1/2). T is the temperature and n and s are integers.
The full Green function can be established using Substituting the expression for the Green function into the lowest order contribution to the self energy, Here, the phonon propagator, d 0 (q, ω s ) = δ ij 2 /( 2 + ω 2 s ); hence, there are no off-diagonal elements of the lowest order self energy. The use of a single averaged and λ here is consistent with the mean-field approximation. At half-filling, it is reasonable to assume that λ A = λ B so that A = B and Z A = Z B . This leads to the gap equations, where the full gap is n =¯ n + . The density of states for a tight binding hexagonal lattice in the absence of a gap, D( ), has the form given in reference [14]. The equations http://www.nanoscalereslett.com/content/7/1/303 may be solved self-consistently by performing a truncated sum on Matsubara frequencies. Gap and quasi-particle weight functions only have a weak Matsubara frequency dependence (< 0.3% for λ = 1, k B T = = 0.01t). The local gap enhancement factors / , are shown in Figure 2 (a) for various λ, showing a modest increase of around 70% for λ = 1. The enhancement factor increases slightly with decreasing but is essentially unchanged by modifications to phonon frequency and temperature for the parameter values used here. The temperature dependence of the gap was also calculated in Figure 2   = 0.06t = 140 meV, = 0.09t = 210 meV, covering the full range of phonon frequencies in reference [13]. k B T = 0.01t (T=268K) and λ ≤ 1. (b) Variation of the gap with temperature, = 0.05t = 117 meV and λ = 0.2. There is a weak temperature dependence due to the large , consistent with the measurements in reference [15], with gap starting to close only for extremely high temperatures T > 8, 000 K, presumably above the melting point of the material. Red circles show the size of the gap, red lines are a guide to the eye.

Conclusions
A theory for the modification of BN band-gaps by interaction with phonons was presented here. It is of interest to make a comparison between the bandgaps of bulk h-BN, nanotubes, monolayer h-BN, and the theory presented here. Measured bandgaps of bulk h-BN are of between 5.8 eV [15] and 5.971 eV [11], indicating that interaction between layers increases the bandgap, consistent with the theory here. The bulk gap is also higher than that for nanotubes (5 eV) [16]. On the other hand, Song et al. [12] claim that the gap is reduced as BN thickness increases. The above discussion is presented with a warning that the theory requires that hopping between the substrate and the BN monolayer is small. Interlayer hopping will affect the bandwidth and bandgap, and the direct Coulomb interaction with strongly ionic substrates could also affect the band structure if the charge density at the surface of the substrate varies dramatically.
It is also of interest to estimate the magnitude of the bandgap modification due to electron-phonon interaction in isolated monolayers of BN. Ab initio calculations have attempted to quantify the magnitude of the interaction between electrons and acoustic phonons for small momentum excitations [17]. Extrapolating the interaction and taking a mean-field average (assuming mean momentum magnitude of 4π/9a), the electron-phonon coupling can be estimated as λ = (4π/9) 2 E 2 1 /2a 2 tM 2 , taking E 1 = 3.66 eV from reference [17],M ≈ 12.5 amu, a = 2.5Å. The mean energy of longitudinal acoustic phonons lies in the range of 50 to 75 meV, giving a range of λ = 0.05 to 0.12, so the contribution of phonons to the bandgap is estimated as 3% to 7%. I would expect BN to have stronger interaction with optical phonons, since the pattern of distortions around an electronic defect is consistent with optical modes (see Figure 1).
The BN gap is too wide for digital applications. Recently, it has become possible to manufacture silicene, an atomically thick layer of silicon with similar properties to graphene [18], so it may be possible to make GaAs or AlP analogues to BN. Smaller gaps could be available from those materials, which might be used to create tunable bandgaps for atomically thick transistors. http://www.nanoscalereslett.com/content/7/1/303