SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

Open Access Nano Express

Gap modification of atomically thin boron nitride by phonon mediated interactions

James P Hague

Author affiliations

Department of Physical Sciences, The Open University, Walton Hall, Milton Keynes, MK7 6AA, UK

Citation and License

Nanoscale Research Letters 2012, 7:303  doi:10.1186/1556-276X-7-303


The electronic version of this article is the complete one and can be found online at: http://www.nanoscalereslett.com/content/7/1/303


Received:17 January 2012
Accepted:23 April 2012
Published:14 June 2012

© 2012 Hague; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License(http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

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.

Keywords:
Boron nitride; Electron-phonon interactions; Semiconductors; Two-dimensional materials; Graphene

Background

The need for bandgaps in graphene on electronvolt scales has led to a number of proposals, such as the use of bilayer graphene [1], creation of nanoribbons [2], and manipulation through substrates [3,4]. Recently, it has become possible to manipulate atomically thin layers of boron nitride (BN) and other materials with structure similar to graphene [5]. This may lead to a complimentary method of manipulating bandgaps to make digital transistors.

In low dimensional materials, strong effective electron-electron interactions can be induced via an interaction between electrons confined to a plane and phonons in a polarizable neighboring layer [6]. The theory has shown that similar interactions account for the transport properties of graphene on polarizable substrates [7] and that sandwiching graphene between polarisable superstrates and gap opening substrates can cause gap enhancement [8]. This paper examines similar gap changes in atomically thin BN due to interactions mediated through substrates.

Methods

Atomically thick hexagonal BN (h-BN) has similar chemistry to graphene: bonding occurs through sp2 hybridization, and electrons with energies close to the chemical potential are in unhybridized π orbitals [9]. A key difference is that the electronic charge is not completely screened by the sp2 hybridization, shifting π orbitals by Δ n = + Δ on N sites, and −Δ on B sites. This shift is the 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, <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/303/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/303/mathml/M1">View MathML</a>, 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:

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/303/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/303/mathml/M2">View MathML</a>

(1)

The Hamiltonian terms are shown schematically in Figure 1 (a). <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/303/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/303/mathml/M3">View MathML</a> creates electrons of spin σ on B sites and <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/303/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/303/mathml/M4">View MathML</a> 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 (b) shows the forces on ions from an increase in electron density at a B site. The largest forces are on the near-neighbor sites, so that the effective interaction is mainly site diagonal (electrons on A sites self-interact through phonons on B sites and vice versa). The diagram indicates that the strongest interaction is between electrons and optical phonon modes.

thumbnailFigure 1 . Bn substrate system and interactions in a monolayer of BN.(a) BN substrate system annotated with interactions. Electron-phonon interactions between the BN layer and substrate are poorly screened, and large interactions of strength f n (m ) are possible. Ions in the substrate oscillate with frequency Ω. N sites have energy + Δ and B sites −Δ , opening a gap. The attractive phonon-mediated electronic interaction f binds electrons onto the same site, effectively enhancing the gap. (b) Interactions in a monolayer of BN. Red circles represent N atoms and black circles B atoms. Light blue arrows represent distortions expected from an excess of charge on the site labelled e .

For simplicity, the Holstein electron-phonon interaction was used, <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/303/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/303/mathml/M5">View MathML</a>, 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,

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/303/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/303/mathml/M6">View MathML</a>

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 <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/303/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/303/mathml/M7">View MathML</a> is the gap function. For bosonic quantities, <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/303/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/303/mathml/M8">View MathML</a> and for fermions, <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/303/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/303/mathml/M9">View MathML</a>. T is the temperature and n and s are integers.

The full Green function can be established using Dyson’s equation <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/303/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/303/mathml/M10">View MathML</a>, leading to,

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/303/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/303/mathml/M11">View MathML</a>

(2)

Substituting the expression for the Green function into the lowest order contribution to the self energy,

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/303/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/303/mathml/M12">View MathML</a>

(3)

Here, the phonon propagator, <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/303/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/303/mathml/M13">View MathML</a>; 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,

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/303/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/303/mathml/M14">View MathML</a>

(4)

<a onClick="popup('http://www.nanoscalereslett.com/content/7/1/303/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/303/mathml/M15">View MathML</a>

(5)

where the full gap is <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/303/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/303/mathml/M16">View MathML</a>. 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 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 <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/303/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/303/mathml/M17">View MathML</a>). 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 (b). For very large temperatures, where kBT approaches Δ , there is a drop in the gap size. Below approximately 8,000 K (0.3t ), this levels off, and the gap becomes relatively constant.

thumbnailFigure 2 . Modification of the BN bandgap.(a) The gap enhancement depends mainly on λ , is weakly dependent on Δ and shows almost no change with Ω. Calculations are made for Δ =t corresponding to a BN gap of 2Δ =4.66 eV, Δ =1.20t (2Δ =5.6 eV), and Δ =0.84t (2Δ =3.92 eV, the tight binding fit from reference [10]). t =2.33 eV, <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/303/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/303/mathml/M18">View MathML</a> meV, <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/303/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/303/mathml/M19">View MathML</a> meV, <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/303/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/303/mathml/M20">View MathML</a> meV, <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/303/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/303/mathml/M21">View MathML</a> 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, <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/303/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/303/mathml/M22">View MathML</a> 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 <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/303/mathml/M23','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/303/mathml/M23">View MathML</a>, taking E 1=3.66 eV from reference [17], <a onClick="popup('http://www.nanoscalereslett.com/content/7/1/303/mathml/M24','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/7/1/303/mathml/M24">View MathML</a> 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.

Competing interest

The author declares that he has no competing interest.

Author’s information

JPH is a lecturer from the Faculty of Science, Department of Physical Sciences, The Open University, Walton Hall, Milton Keynes, UK.

Acknowledgements

The author acknowledges EPSRC grant EP/H015655/1 for funding and useful discussions with A Ilie and A Davenport.

References

  1. McCann E, Fal’ko VI: Landau-level degeneracy and quantum hall effect in a graphite bilayer.

    Phys Rev Lett 2006, 96:086805. PubMed Abstract | Publisher Full Text OpenURL

  2. Brey L, Fertig HA: Electronic states of graphene nanoribbons studied with the Dirac equation.

    Phys Rev B 2006, 73:235411. OpenURL

  3. Zhou SY, Gweon GH, Fedorov AV, First PN, Heer WAD, Lee DH, Guinea F, Neto AHC, Lanzara A: Substrate-induced bandgap opening in epitaxial graphene.

    Nature Materials 2007, 6:770. PubMed Abstract | Publisher Full Text OpenURL

  4. Enderlein C, Kim YS, Bostwick A, Rotenberg E, Horn K: The formation of an energy gap in graphene on ruthenium by controlling the interface.

    New J Phys 2010, 12:033014. Publisher Full Text OpenURL

  5. Novoselov KS, Jiang D, Schedin F, Booth TJ, Khotkevich VV, Morozov SV, Geim AK: Two-dimensional atomic crystals.

    Proc Natl Acad Sci USA 2005, 102:10451. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  6. Alexandrov AS, Kornilovitch PE: The Fröhlich-Coulomb model of high-temperature superconductivity and charge segregation in the cuprates.

    J Phys: Condens Matter 2002, 14:5337. Publisher Full Text OpenURL

  7. Fratini S, Guinea F: Substrate-limited electron dynamics in graphene.

    Phys Rev B 2008, 77:195415. OpenURL

  8. Hague J: Tunable graphene band gaps from superstrate-mediated interactions.

    Phys Rev B 2011, 84:155438. OpenURL

  9. Alem N, Erni R, Kisielowski C, Rossell M, Gannett W, Zettl A: Band-gap engineering and ballistic transport in edge-corrugated graphene nanoribbons.

    Phys Rev B 2009, 80:155425. OpenURL

  10. Ribeiro RM, Peres NMR: Stability of boron nitride bilayers: ground-state energies, interlayer distances, and tight-binding description.

    Phys Rev B 2011, 83:235312. OpenURL

  11. Watanabe K, Taniguchi T, Kanda H: Direct-bandgap properties and evidence for ultraviolet lasing of hexagonal boron nitride single crystal.

    Nat Mater 2004, 3:404. PubMed Abstract | Publisher Full Text OpenURL

  12. Song L, Ci L, Lu H, Sorokin P, Jin C, Ni J, Kvashnin A, Kvashnin D, Lou J, Yakobson B, Ajayan P: Large scale growth and characterization of atomic hexagonal boron nitride layers.

    Nano Lett 2010, 10:3209. PubMed Abstract | Publisher Full Text OpenURL

  13. Serrano J, Bosak A, Arenal R, Krisch M, Watanabe K, Taniguchi T, Kanda H, Rubio A, Wirtz L: Vibrational properties of hexagonal boron nitride: inelastic X-Ray scattering and ab initio calculations.

    Phys Rev Lett 2007, 98:095503. PubMed Abstract | Publisher Full Text OpenURL

  14. Neto AHC, Guinea F, Peres NMR, Novoselov KS, Geim AK: The electronic properties of graphene.

    Rev Mod Phys 2009, 81:109. Publisher Full Text OpenURL

  15. Zunger A, Katzir A, Halperin A: Optical properties of hexagonal boron nitride.

    Phys Rev B 1976, 13:5560. Publisher Full Text OpenURL

  16. Terauchi M, Tanaka M, Matsumoto T, Saito Y: Electron energy-loss spectroscopy study of the electronic structure of boron nitride nanotubes.

    J Electron Microscopy 1998, 47:319. Publisher Full Text OpenURL

  17. Bruzzone S, Fiori G: Ab-initio simulations of deformation potentials and electron mobility in chemically modified graphene and two-dimensional hexagonal boron-nitride.

    Appl Phys Lett 2011, 99:222108. Publisher Full Text OpenURL

  18. Padova PD, Quaresima C, Ottaviani C, Sheverdyaeva P, Moras P, Carbone C, Topwal D, Olivieri B, Kara A, Oughaddou H, Aufray B, Lay GL: Evidence of graphene-like electronic signature in silicene nanoribbons.

    Appl Phys Lett 2010, 96:261905. Publisher Full Text OpenURL