Structural and electronic properties of epitaxial multilayer h-BN on Ni(111) for spintronics applications.

Hexagonal boron nitride (h-BN) is a promising material for implementation in spintronics due to a large band gap, low spin-orbit coupling, and a small lattice mismatch to graphene and to close-packed surfaces of fcc-Ni(111) and hcp-Co(0001). Epitaxial deposition of h-BN on ferromagnetic metals is aimed at small interface scattering of charge and spin carriers. We report on the controlled growth of h-BN/Ni(111) by means of molecular beam epitaxy (MBE). Structural and electronic properties of this system are investigated using cross-section transmission electron microscopy (TEM) and electron spectroscopies which confirm good agreement with the properties of bulk h-BN. The latter are also corroborated by density functional theory (DFT) calculations, revealing that the first h-BN layer at the interface to Ni is metallic. Our investigations demonstrate that MBE is a promising, versatile alternative to both the exfoliation approach and chemical vapour deposition of h-BN.

The fast development of the modern information technology raises the quest for advanced materials and new methods of nano-devices fabrication. A prominent topical example is graphene, a two-dimensional (2D) crystalline atomic layer system with charge carrier mobilities exceeding by far those of Si1. This ideal two-dimensional electron gas (2DEG) system exhibits unique properties for charge23 and spin4 transport. The spin property prompts applications in memory, sensor, and active spin-current based devices. For the latter the detrimental spin transport parameters like spin relaxation time and relaxation length have called for materials with negligible intrinsic spin-orbit and hyperfine interactions. Here graphene appears as a prime spintronics candidate being gate tunable at room temperature456. Thereby, graphene has also spurred the conjectured change in paradigm from charge-based electronics, i. e. transfer and storage of electron charges, to spintronics7 based on the electron’s spin.

The progress in graphene-based devices8910 triggered the search for a suitable insulating substrate. The crystalline structure of wide band gap hexagonal boron nitride (h-BN) is isomorphic to the structure of graphite making h-BN an ideal match to graphene. Weak van der Waals bonds between honeycomb h-BN sheets and the absence of surface states allow easy manipulations with exfoliated single or multilayers of h-BN. Owing to these properties new physical phenomena11121314 and new devices15161718 became possible with exfoliated h-BN. In particular, superior charge and spin transport properties in graphene on h-BN1920 over graphene on amorphous SiO221 are reported. In this case, the smoothening of graphene flakes by h-BN flakes has enabled significant improvements of spin lifetimes, spin relaxation lengths, and electron mobilities in Co/MgO/graphene/h-BN spin-valve devices20.

Moreover, for spin injection into graphene it has been realized that h-BN with a large band gap of ≈6 eV2223 has the potential of replacing defective oxide tunnel barriers2425. In view of early predictions by Kirczenow26 for metal/h-BN/semiconductor interfaces the analogous barrier replacement could be assumed for magnetic tunnel junctions (MTJs). Currently, MTJs are made of ferromagnetic electrodes separated by Al2O3 or a hygroscopic MgO barrier27. An attempt of making MTJs with chemical vapour deposition (CVD) grown h-BN tunnel barriers has resulted in a rather low tunneling magnetoresistance (TMR) of less than 1%28. Such TMR is far below predicted values29 due to contaminations during h-BN transfer.

All these examples of excellent physical properties and specific device demonstrations by applying exfoliated h-BN are preceding the next stage in the development of van der Waals materials: large scale and integrated device fabrication, requiring wafer scale fabrication methods of h-BN. Recently, several approaches of epitaxial h-BN growth on close-packed surfaces of polycrystalline transition metals have been proposed including borasine-free chemical vapour deposition3031, sputtering32, and molecular beam epitaxy (MBE)33. In these publications mainly structural properties of epitaxial h-BN films have been addressed, while their electronic structures have not been reported yet. Our work addresses this topic by presenting a comprehensive characterization of large scale epitaxial h-BN films complemented by density-functional theory (DFT) modeling.

We report on the epitaxial growth of multilayered h-BN (ml-h-BN) on the close-packed Ni(111) surface of single crystal Ni by means of MBE. The obtained films of ml-h-BN demonstrate very high structural perfection with A-B stacking of BN layers as observed by means of cross-sectional transmission electron microscopy (TEM). Photoelectron spectroscopy experiments demonstrate clearly the energy dispersion of the valence band states below the Fermi level (E) of the system. These data in combination with near-edge x-ray absorption spectroscopy (NEXAFS) results and DFT calculations allow for the systematic mapping of the electronic structure of the system. Our findings demonstrate that MBE growth of ml-h-BN is a perspective method that can be easily adapted to electronic and spintronic applications of h-BN such as, e. g., advanced devices for spin injection into graphene or for MTJs. These results demonstrate the advanced approach compared to previous attempts to grow epitaxial h-BN on ferromagnetic substrates by means of CVD from borasine as a precursor34. In this work authors performed studies of the electronic structure of 1 ML-thick h-BN on Ni(111) by means of photoelectron spectroscopy of core levels and valence band states. It was also found that thicker (>1 ML) BN films do not show epitaxial growth as was indicated in low-energy electron diffraction (LEED) experiments.

Results

h-BN/Ni(111) growth and characterization

The epitaxial ml-h-BN films were grown using the MBE technique [Fig. 1(a)], by which boron atoms were deposited in the presence of nitrogen atoms, generated in an RF-plasma source, on Ni(111) substrates kept at 700 °C (for details, see section Methods). The results of the structural characterization of the ml-h-BN/Ni(111) stack are compiled in Fig. 1 for a 10 ML-thick h-BN film. A reflection high-energy electron diffraction (RHEED) pattern taken in situ immediately after BN deposition reveals well-defined stripes [see inset of Fig. 1(a) and Fig. S1 of the Supplementary material]. No signatures of amorphous (diffuse-type pattern) or polycrystalline (ring-shaped pattern) phases are observed. The elongated shape of the stripes could be attributed to a rough surface, which is typical for commercially available single crystals of metals. This pattern allows a preliminary conclusion that the epitaxial growth of BN takes place and that the surface morphology resembles that of single crystal Ni(111). High-resolution cross-sectional TEM images are presented in panels (b) and (c). These images were collected in the <110> direction of bulk Ni. One can clearly see that the epitaxial growth of h-BN on Ni(111) is obtained with the A-B stacking of single BN layers [see inset of Fig. 1(c)]. The TEM-measured interlayer distance in h-BN is 3.86 ± 0.32 Å while the h-BN-Ni separation equals 1.87 ± 0.12 Å. The corresponding extracted in-plane lattice constant of BN is 2.44 ± 0.2 Å. These values are close to the respective parameters of bulk h-BN: 2.504 Å for in-plane lattice constant and 3.33 Å for the distance between layers35. After synthesis, the ml-h-BN/Ni(111) sample was transferred under atmospheric conditions to other experimental stations for surface-sensitive experiments. The high inertness of BN to the ambient environment allowed to avoid irreversible oxidation and contamination of the substrate36.

Figure 1

Structural characterization of h-BN/Ni(111).

(a) Step-by-step preparation of the h-BN/Ni(111) structures by MBE. RHEED image of this system is shown as an inset. (b,c) High-resolution TEM images of the 10 ML-thick BN on Ni(111). Inset in (c) shows a zoomed image of the ml-h-BN demonstrating the A-B stacking of BN layers. (d) LEED image of the same sample obtained after atmospheric transfer of ml-h-BN/Ni(111) into a UHV chamber and annealing at 400 °C. The energy of the primary electron beam is 74 eV.

Figure 1(d) shows a LEED image of the ml-h-BN/Ni(111) system obtained after atmospheric transfer of the sample into a UHV chamber and the respective annealing step at 400 °C in order to remove possible air-adsorbates (water and CO). The hexagonal LEED pattern measured at the primary electron beam energy of 74 eV indicates the in-plane long-range order of the studied sample.

Electronic structure of ml-h-BN/Ni(111)

The chemical composition of the sample was verified by x-ray photoelectron spectroscopy (XPS). The results are shown in Fig. 2(a). This figure shows XPS spectra (survey spectra and insets with detailed N 1s and B 1s spectra) recorded (i) after the sample was transferred to the UHV chamber (brown lines) and (ii) after it was annealed in UHV at 400 °C for 30 min (blue lines). This simple temperature treatment leads to the complete disappearance of the C 1s peak and thus to desorption of the carbon-containing molecules (CO, CO2) from the surface of h-BN, and the XPS signal originating from oxygen is largely decreased.

Figure 2

XPS and NEXAFS of h-BN/Ni(111).

(a) XPS survey spectra obtained with Al Kα radiation for as-transferred into UHV h-BN/Ni(111) sample (i) and after annealing the same sample at 400 °C (ii). The main core-level and Auger emission lines are marked in the Figure. Insets show the corresponding highly resolved XPS spectra of the N 1s and B 1s core levels. (b) TEY NEXAFS spectra of h-BN/Ni(111) collected as functions of photon energy at the K absorption edges of N and B for different incident angles α between the sample plane and the incoming light (angles are marked in the figure).

The positions of the N 1s and B 1s XPS emission lines of the clean surface area are 399.25 eV and 191.35 eV, respectively, and these values are in a very good agreement with previously published data for the thick h-BN layer formed on polycrystalline surfaces of Cu and Ni3637. The analysis of the h-BN film composition on the basis of the corresponding XPS line intensities gives a stoichiometric compound within the error of the XPS measurements (±1%). The subsequent angular-resolved XPS (AR-XPS) measurements of the intensity of Ni 2p3/2 and N 1s emission lines has been carried out (see Fig. S2 of Supplementary material). In this measurement the probed thickness of the BN film is varied as a function of the angle θ between the sample plane and the direction to the energy analyzer. As a result the AR-XPS method gives a thickness of the studied h-BN film of 40 ± 5 Å, which is in a good agreement with TEM data (34.7 ± 3 Å) in Fig. 1.

Further confirmation of the high crystalline quality of the obtained h-BN is revealed by NEXAFS experiments. The NEXAFS spectra were acquired in the total electron yield mode (TEY, sample drain current) for the annealed samples at the K absorption edges of N and B as a function of angle α between sample plane and the direction of the incoming light (Fig. 2(b), bottom). These spectra demonstrate an example of the so-called search-light-like effect38, when the absorption intensity depends on the relative orientation of the electric field vector of the incoming light and the alignment of the orbital, onto which the electron is transferred from the respective core level. Therefore, the NEXAFS peaks at 398–405 eV for N K and 190–194 eV for B K are assigned to the transitions of 1s core electrons of corresponding atoms into the antibonding π* band of h-BN39. On the other hand, the high energy tails of both NEXAFS spectra are assigned to the transitions of 1s electrons into the σ* bands of h-BN39. The presented spectra are very similar to the ones of bulk h-BN40. These signatures underscore the ordered crystalline structure of our epitaxial h-BN on Ni(111) having the electronic structure similar to that of bulk h-BN. The pre-edge structure in the NEXAFS spectra for both absorption edges originates from the interface states between h-BN and Ni(111) as a result of the orbital mixing between Ni 3d and BN π states40.

An additional investigation of the electronic structure below E of ml-h-BN/Ni(111) was carried out by means of angle-resolved PES (ARPES). Figure 3(a) shows constant energy cuts obtained via combination of two data sets collected as a series of angle-resolved photoemission maps along Γ − M and Γ − K directions of the Brillouine zone (BZ) of the studied system. The emission from the h-BN π states is found in the valence band. This emission band demonstrates the characteristic trigonal warping around the K point [Fig. 3(a)]. The corresponding emission intensity map taken as a wave vector cut along the black line passing through the K point, marked in the panel (a), is shown in (b) as intensity distribution. This map and the corresponding one around the M point allow to extract the respective intensity profiles shown in Fig. 3(c). These data allow to find positions of the valence band maximum for different high symmetry points of the BZ: E = 3.70 ± 0.02 eV for the K point and E = 4.73 ± 0.02 eV for the M point. The weak band dispersion observed at higher binding energies might be assigned to the faint photoemission intensity from the deep interface BN layer (see section Discussion).

Figure 3

ARPES and DFT for h-BN/Ni(111).

(a) Series of the ARPES constant binding energy cuts for ml-h-BN/Ni(111). Data are extracted from the two combined data sets consecutively acquired via polar angle scanning of the sample along the Γ − K and Γ − M directions of the hexagonal BZ. The corresponding binding energies (E) for the energy cuts and the high symmetry points of the BZ are marked in the figure. (b) ARPES intensity map extracted from the 3D data set along the black line in the k-space as marked in (a). (c) PES intensity profiles corresponding to the K and M points extracted from the 3D data sets. (d) Calculated binding energies of the π band of ml-h-BN on Ni(111) as a function of BN-thickness at the K (squares) and M (circles) points obtained in DFT calculations. Dashed horizontal lines are experimental values from (c).

Discussion

The experimental results presented above were compared with the DFT calculations of electronic structure modifications of ml-h-BN/Ni(111) depending on the thickness of h-BN (for technical details, see section Methods). Figure 3(d) shows the evolution of the binding energy of the BN π band at the K and M points with the h-BN thickness. Despite the difficulties in the modeling of insulating compounds by means of DFT (underestimation of the band gap41), a careful selection of hybrid density functionals4243 resulted in reasonable theoretical values of K and M point binding energies for thick h-BN layer (>2 ML), which are very close to the experimental ones (shown by the dashed lines of the corresponding colour in Fig. 3(d)).

A single h-BN layer is adsorbed on Ni(111) in the top-fcc arrangement as deduced from our DFT calculations (see section Methods), where the N-atom of h-BN is placed above the Ni-atom of the interface layer, while the B-atom occupies the fcc hollow adsorption site of Ni(111) (see Supplementary material, Fig. S3). The calculated mean distance between h-BN and Ni(111) is 2.04 Å. Note: N- and B- atoms are separated from the Ni interface differently:  Å, d(B − Ni − top) = 1.98 Å, the corrugation is 0.12 Å. This result is consistent with the ones of the previously published works3944. The resulting band structure of 1 ML-h-BN/Ni(111) along the high symmetry directions of the hexagonal BZ is shown in Fig. 4(a), where it is overlaid with the weight of the BN π band shown by blue lines. The binding energy of BN π states is 5.04 eV and 6.33 eV for the K and M points, respectively. The interaction between the monolayer-thick h-BN layer and Ni(111) leads to the formation of several interface states as a result of energy-, space-, and k-overlapping of the BN π orbitals and Ni 3d states having z-components in the angular part of the wave function. This effect leads to the metallic character of 1 ML-thick h-BN (sizeable weight of BN π states at E) and pushing down the π bands to larger binding energies around K and M points.

Figure 4

DFT results for h-BN on Ni(111).

Calculated spin-up band structures of (a) 1 ML-thick h-BN and (b) 2 ML-thick h-BN on Ni(111) along high-symmetry directions of the Brillouine zone (spin-up channel). Blue lines mark the weight of the π band of the top BN layer.

In case of the A-B stacked double layer of h-BN on Ni(111), the second layer is completely decoupled from the underlying substrate, that leads to the restoring of the insulating nature of h-BN (Fig. 4(b)). Increasing the number of h-BN layers does not lead to qualitative changes in the DFT-obtained band structure. The A-B stacking of h-BN is found in our high-resolution TEM measurements [see Fig. 1(b)] and is confirmed in the respective DFT calculations as the most stable configuration with the distance between BN layers of 3.07 Å. The binding energy of BN π states is 3.1 eV and 4.1 eV for the K and M points, respectively. For these points the calculated band gap for 2 ML-h-BN/Ni(111) is 6.26 eV and 6.79 eV. These values are very close to the experimentally measured band gap for h-BN222345. It is interesting to note that while the structure and binding properties (like e. g. interlayer interaction energy) of h-BN as well as hybrid materials on the basis of h-BN can get reliable description nowadays when employing standard DFT functionals augmented by the dispersion correction (see e.g. refs 46,47), the reported DFT band gaps are often underestimated with respect to the experiment (e. g. ref. 48 predict E = 4.56 eV for the h-BN sheet49). The Heyd-Scuseria-Ernzerhof (HSE) functional used in the present study predicts the h-BN gap very close to experiment4950.

Conclusion

We have shown that the epitaxial growth of an h-BN multilayer (ml-h-BN) is feasible on the close-packed Ni(111) surface of a Ni single crystal by means of MBE. The obtained epitaxial films of ml-h-BN demonstrate very high structural perfection with A-B stacking of BN layers as observed using cross-sectional TEM. Photoelectron spectroscopy experiments demonstrate clearly the energy dispersion of the valence band states below the Fermi level. These data in combination with results of NEXAFS and DFT calculations allow for the systematic mapping of the electronic structure of the system. The DFT calculations show good agreement with the experimental data and with bulk electronic properties of h-BN. DFT also reveals that the first h-BN layer at the interface of Ni/ml-h-BN is metallic, thus placing the metal/insulator interface inside of ml-h-BN. Our findings demonstrate that MBE growth of ml-h-BN is a perspective method that can be easily adapted to electronic and spintronic applications of h-BN. We propose that this epitaxial Ni(111)/ml-h-BN system, can be used as an effective substrate for the fabrication of FM/BN/graphene spin injection devices or FM/BN/graphene/BN/FM spin switches (FM = ferromagnet), where graphene is effectively decoupled from the metallic substrate. This allows to preserve the unique band structure of graphene in the vicinity of E.

Methods

Samples preparation

The growth of h-BN films was carried out by applying molecular beam epitaxy (MBE, Riber SIVA 45 setup) on a Ni(111) surface. The growth recipe was combined out of two growth approaches described elsewhere3351. A single crystalline Ni(111) was used as a substrate. The cleaning of the substrates consisted of a heating step to the temperature of approximately 800 °C followed by the deposition of a 20 nm thick Ni buffer layer using e-beam evaporation at the substrate temperature of 700 °C. (No sputtering was used during cleaning of Ni(111) substrates.) Then, an equivalent of one monolayer (ML) amount of B was deposited on the substrate followed by the deposition of BN at 700 °C. The flux of boron atoms was generated by using a high-temperature effusion cell operated at 1900 °C. A preliminary calibration of the B-flux was carried out by applying secondary ion mass spectroscopy on a uniformly doped Si:B layer grown on a Si(001) substrate. The flux of nitrogen was obtained by an RF-plasma source (Riber RF-N 50–63) at the RF-power of 250 W and the MBE background pressure of 1.3 × 10−5 Torr. The growth rate of BN was 0.14 ML/min.

Sample characterization

The as-grown h-BN films were characterized in situ by using reflection high-energy electron diffraction (RHEED) and Auger electron spectroscopy (AES). The ex situ characterization was carried out outside the MBE system requiring sample transfer through the atmosphere and, in some cases, additional preparations. The morphology and crystal structure of the films were analyzed by transmission electron microscopy (TEM) and scanning TEM (STEM). The corresponding specimens of specific sample orientations have been generated by Focused Ion Beam (FIB) preparation technique. TEM investigations were carried out in a TITAN 80/300 electron microscope. A low-energy electron diffraction (LEED) was carried out in the SPECS UHV probe station. In some cases of ex situ characterization (LEED, XPS, ARPES) the samples were annealed at 400 °C. A set of surface-sensitive characterization methods of the h-BN electronic structure is carried out and described below.

NEXAFS experiments

These experiments were performed at the Russian-German beamline of the BESSY II storage ring of HZB, Berlin. X-ray absorption spectroscopy experiments were performed in the total electron yield (TEY) mode via measurement of the drain current. Spectra were collected around K absorption edges of N and B with the total resolution of the beamline of 50 meV. Absorption experiments performed at different angles between sample surface and the direction of the incoming light allow to verify spacial orientation of the valence band states via the so-called search-light-like effect.

ARPES experiments

The ARPES measurements with He II radiation (photon energy of hv = 40.81 eV) were performed in the SPECS demo lab using a FlexPS system with PHOIBOS 150 2D-CCD analyzer. In this case a 5-axis motorized manipulator was used, allowing for a precise alignment of the sample in k space. The sample was azimuthally pre-aligned in such a way that the polar scans were performed along the Γ − K or Γ − M directions of the hexagonal Brillouine zone of h-BN/Ni(111) with the photoemission intensity on the channelplate images acquired along the direction perpendicular to Γ − K or Γ − M directions, respectively. The final 3D data sets of the photoemission intensity as a function of kinetic energy and two emission angles, I(E, angle1, angle2), were then carefully analyzed.

DFT calculations

Calculations were performed using the plane-wave projector augmented-wave (PAW)5253 method applying the semi-local Perdew-Burke-Ernzerhof (PBE)48 exchange-correlation functional augmented by the DFT-D2 method of Grimme54 for describing the dispersion interactions and the Heyd-Scuseria-Ernzerhof (HSE)55 hybrid functional as implemented in the Vienna ab initio simulation package (VASP)5657 code. The HSE screening parameter was set to a value of 0.2 Å−1 58. The plane-wave cutoff energy was set to 500 eV. Brillouin-zone integration was performed on Γ-centered symmetry reduced Monkhorst-Pack59 meshes using a Methfessel-Paxton smearing method of first order with σ = 0.2 eV, except for the calculation of total energies and densities of states (DOSs). For those calculations, the tetrahedron method with Blöchl corrections60 was used. A 12 × 12 × 1 k-mesh was used. Although for both computational schemes (HSE and PBE) the qualitative results are found to be similar, the HSE functional predicts the band gaps reasonably well, while the PBE functional yields a significant underestimation of the band gaps. ml-h-BN/Ni(111) was modelled, in case of relaxation (band structure) calculations, as a stack of 13 ML (7 ML) separated by a spacer vacuum of 45 Å with the corresponding number of h-BN layers in A-B stacking on both (one) sides of the Ni-slab.

Additional Information

How to cite this article: Tonkikh, A. A. et al. Structural and electronic properties of epitaxial multilayer h-BN on Ni(111) for spintronics applications. Sci. Rep. 6, 23547; doi: 10.1038/srep23547 (2016).