Determination of the thickness and orientation of few-layer tungsten ...

12 downloads 34625 Views 944KB Size Report
Technology, Ulsan 44919, Korea. E-mail: [email protected] ... Polarization dependence of Raman spectrum of bulk WTe2 sample. (a) Polarized Raman ...
Determination of the thickness and orientation of few-layer tungsten ditelluride using polarized Raman spectroscopy Minjung Kim1, Songhee Han1, Jung Hwa Kim2, Jae-Ung Lee1, Zonghoon Lee2 and Hyeonsik Cheong1,* 1

Department of Physics, Sogang University, Seoul 04107, Korea

2

School of Materials Science and Engineering, Ulsan National Institute of Science and Technology, Ulsan 44919, Korea

E-mail: [email protected] Keywords: Tungsten ditelluride (WTe2), Raman spectroscopy, polarization, orientation, thickness

Abstract. Orthorhombic tungsten ditelluride (WTe2), with a distorted 1T structure, exhibits a large magnetoresistance that depends on the orientation, and its electrical characteristics changes from semimetallic to insulating as the thickness decreases. Through polarized Raman spectroscopy in combination with transmission electron diffraction, we establish a reliable method to determine the thickness and crystallographic orientation of few-layer WTe2. The Raman spectrum shows a pronounced dependence on the polarization of the excitation laser. We found that the separation between two Raman peaks at ~90 cm−1 and at 80–86 cm−1, depending on thickness, is a reliable fingerprint for determination of the thickness. For determination of the crystallographic orientation, the polarization dependence of the

1

A1 modes, measured with the 632.8-nm excitation, turns out to be the most reliable. We also discovered that the polarization behaviors of some of the Raman peaks depend on the excitation wavelength as well as thickness, indicating a close interplay between the band structure and anisotropic Raman scattering cross section.

2

1. Introduction Tungsten ditelluride (WTe2), a layered 2-dimensional material in the transition metal dichalcogenide (TMD) family, has attracted much attention since a large, non-saturating magnetoresistance (MR) was reported [1–9]. Unlike other TMD materials such as MoS2 which are semiconductors, WTe2 is a semimetal in bulk and becomes insulating in the few-layer limit [8]. Furthermore, because WTe2 has an orthorhombic, distorted-1T structure, its physical properties show pronounced anisotropy in the plane of the layers. For example, the MR behavior shows a strong anisotropy: when the magnetic field is applied to the c axis, MR is maximum when the current is parallel to the a axis and minimum along the b axis [1,3,5,6]. For these reasons, it is important to determine the thickness and the crystallographic orientation of WTe2 when this material is used in MR devices, etc. Polarized Raman spectroscopy is a powerful method for determining a number of layers and crystallographic orientation of layered 2-dimensional materials [2,8,10–27]. Several Raman studies on few-layer WTe2 have been reported by several groups to correlated the Raman spectrum with the number of layers [8,23,26–28]. However, as will be shown below, some of the peaks in the spectrum are convolutions of two peaks and the line shape changes depending on the polarization. If this polarization dependence is not considered, determination of the layer thickness based on the position of the Raman peaks can be inaccurate. Polarized Raman spectroscopy was carried out on bulk WTe2 to assign Raman modes [2], but no such work is reported yet for few-layer WTe2. In this work, we investigated the polarization dependence of the Raman spectrum as a function of the thickness of WTe 2 and the excitation wavelength. We propose a reliable method to determine the number of layers and the crystallographic orientation of WTe2 using polarized Raman spectroscopy.

3

2. Results and Discussion 2.1. Determination of crystallographic orientation

Figure 1. Crystallographic orientation of bulk WTe2 sample. (a) Crystal structure of WTe2. (b) Low magnification TEM image and (c) TED pattern of a bulk WTe2 sample. Figure 1(a) shows the crystal structure of WTe2 with orthorhombic, distorted-1T structure with space group Pmn21 [2,21,23,26,28,29]. The a-axis is defined as the zigzag direction of the W-W chains. Figure 1(b) is a low magnification transmission electron microscope (TEM) image of a thick (bulk) WTe2 sample on TEM grid, and Figure 1(c) is a corresponding transmission electron diffraction (TED) pattern. The a-axis direction determined from the TED pattern is indicated in Figure 1(b). Unlike the case of black phosphorus, another anisotropic 2-dimensional material, long straight edges do not necessarily correspond to a principal axis (Supplementary Figure 1).

4

Figure 2. Polarization dependence of Raman spectrum of bulk WTe2 sample. (a) Polarized Raman spectra of bulk WTe2 in parallel polarization. Polarization angles with respect to the a-axis are indicated. The excitation wavelength is 632.8 nm. (b) Polarization dependence of intensities of peaks P2, P10, and P11. The radial scale of all polar plots in this paper starts at zero. Figure 2(a) shows polarized Raman spectra of this sample in parallel polarization, using the excitation wavelength of 632.8 nm. The polarization angles are measured with respect to the a-axis direction determined from the TED measurement. The spectrum varies significantly depending on the polarization, and so the spectra taken in three representative polarizations are shown. For convenience, those peaks that are clearly resolved are numbered from the lowest frequency peak (P1). The optical phonon modes of space group Pmn21 are A1, A2, B1, and B2 [2,21,23,26,28–30], but B1 and B2 modes are forbidden in the back-scattering geometry. The complex Raman tensors, which take absorption into account [31–35], of A1 and A2 modes can be written as

5

 a eia  A1   0  0 

0   0 , c eic 

0 be

ib

0

 0  A2   d eid  0 

d eid 0 0

0  0 . 0 

(1)

The Raman intensity in parallel polarization is given by

I ( A1 )  a cos 4   b sin 4   2 a b cos 2  sin 2  cos ab ,

(2)

I ( A2 )  d cos 2  sin 2  ,

(3)

2

2

and

2

where  ab  a  b , and the angle  is measured with respect to the a-axis. By comparing with Eqs. (2) and (3), we can assign P1, P2, P5, P6, P7, P8, P10, P11, and P12 as A1 modes, and P3, P4, and P9 as A2 modes (Supplementary Figure 2), which is consistent with theoretical predictions [2,21,23,26,28]. The atomic displacements of some of these modes can be found in Refs. [21,23,26]. Figure 2(b) shows that the polarization dependence of the A1 modes can be used to determine the crystallographic orientation of the sample: the polarization direction in which the intensity of P2 is maximum corresponds to the b-axis.

6

2.2. Determination of thickness

Figure 3. Thickness dependence of Raman spectrum of WTe2. (a) Raman spectrum as a function of thickness of WTe2, measured with excitation wavelength of 632.8-nm. For the low-frequency region, the spectra measured with 514.5-nm excitation are also shown. Shear (○) and breathing (□) modes are indicated. The Raman spectrum of bulk WTe2 is magnified 10 times. Thickness dependence of (b) P2 and P3 and (c) P10, P11, and P12 peak positions. (d) Thickness dependence of peak position differences between P2 and P3; P10 and P11; and P10 and P12. Figure 3(a) shows the Raman spectra as a function of thickness. For the low-frequency range (< 35 cm‒1), the Raman spectra for the polarization angle of 0° are shown because the low-frequency peaks have the maximum intensity at this polarization angle (Supplementary Figure 3). For the higher frequency range (>70 cm‒1), the spectra for 45° polarization are shown because some closely spaced peaks (P9 and P10 or P11 and P12) are best resolved at this polarization angle. (See Supplementary Figures 4‒6 for spectra in different polarizations.) The low-frequency region of the spectrum shows the most pronounced

7

thickness dependence. Both shear and breathing modes are observed. The breathing mode is stronger for the 632.8-nm excitation, but the shear mode is only slightly stronger for the 514.5-nm excitation. Both modes have the same polarization dependence (Supplementary Figure 3). The sharp peak at ~9 cm‒1 for bulk WTe2 seems to be an overlap of the breathing mode and a shear mode, but the sharpness of the peak seems to indicate that the contribution from the shear mode is stronger. In the higher frequency range, several peaks show clear shift with thickness. There have been several attempts to correlate the Raman spectrum with the thickness of WTe2 [8,23,26–28]. The difference between the two Raman peaks positioned at ~160 cm‒1 (P10) and ~210 cm‒1 (P11 and P12) has been proposed as an indicator of the thickness of WTe2. However, because the Raman peak near ~210 cm‒1 is in fact a superposition of two separate Raman peaks, P11 and P12, with strong polarization dependence (Supplementary Figure 4), its position is not well defined unless the polarization is specified. The significant scatter of the reported separation between the peak at ~210 cm‒1 and the peak at ~160 cm‒1 (Supplementary Figure 7) can be thus explained. We find that the separation between P2 and P3 also depends on the thickness. Figures 3(b) and (c) summarize the positions of P2, P3, P10, P11, and P12 as a function of thickness. Figure 3(d) compares the difference between these peaks also as a function of thickness. It is clear that all three differences, P3‒P2, P11‒P10, and P12‒P10, correlate well with thickness. However, since one has to deconvolute the peak at ~210 cm‒1 to determine the positions of P11 and P12, their positions are less reliable. A small difference between the data from 632.8-nm and 514.5-nm excitations is a result of this limitation. For this reason, we propose that P3‒P2 is a more reliable measure of the thickness, because these peaks are not convolution of multiple peaks. Also, because P3 is very weak for 514.5 nm but strong and clearly resolved for 632.8 nm, we suggest that the excitation wavelength of 632.8 nm is more preferable (Supplementary Figure 5).

8

2.3. Anomalous behaviors of some Raman peaks

Figure 4. Anomalous polarization behaviors of some Raman peaks. Polarization dependence of intensities of peaks P2, P10, and P11 for different thicknesses measured with (a) 632.8-nm and (b) 514.5nm excitation wavelengths. (c) Polarization dependence of intensities of peaks P2, P10, and P11 for mono-layer WTe2 measured with 632.8-nm and 514.5-nm excitation wavelengths. The points in all polar plots are the normalized Raman intensity data, and the curves are results of fitting to Eq. (2). Figure 4(a) shows the polarization behaviors of some representative peaks measured with the 632.8nm excitation. We take the a-axis direction (0°) of few-layer WTe2 so that the polarization dependences of the Raman peaks coincide with those of bulk WTe2. The polarization behaviors are essentially identical for all thicknesses, which give credence to the above assignment of the a-axis direction. On the other

9

hand, Figure 4(b) shows that for 514.5-nm excitation, there is a strong thickness dependence of the polarization behavior. For example, the polarization behavior of P2 is similar to that for 632.8-nm excitation for all thicknesses except for mono-layer WTe2, which shows an opposite polarization dependence. In the case of P11, the polarization behavior varies gradually from mono-layer to bulk WTe2. (See Supplementary Figures. 6, 8, and 9 for complete data.) Figure 4(c) compares the polarization behaviors of some of the Raman peaks from mono-layer WTe2, measured with excitation wavelengths of 632.8 nm and 514.5 nm. Similar anomalous polarization behavior that depends on the excitation wavelength or the sample thickness was observed in black phosphorus. The anisotropic electron-phonon interaction as well as the interference effect combined with birefringence and dichroism due to anisotropy could partly explain this anomalous behavior [32–34,36]. Unlike black phosphorus, WTe2 has only weak birefringence (Supplementary Figure 10), and so its contribution to the anomalous polarization behavior, if any, would not be large. We propose that dichroism is the major factor that contributes to this anomalous behavior. If we compare the polarization dependence of the Raman intensities with Eq. (2), we find that the variation of the polarization behaviors with the thickness or the excitation wavelength can be accounted for by varying the phase factor, ab  a  b . The phase factors in the Raman tensor elements are related to the imaginary part of the dielectric function, which is directly correlated with optical absorption. Because the band structure determines the optical absorption through the joint density of states, the phase factors reflect the differences in the band structure. The band structure of WTe2 depends sensitively on the number of layers [4,22], the observed thickness dependence can be explained in terms of the thicknessdependent dielectric function. The excitation wavelength dependence is understood in the same way because the anisotropy in the dielectric function would depend on the wavelength for a given thickness. Supplementary Figures 8 and 9 summarize the polarization dependence of all the peaks. The A2 modes do not show any dependence on the excitation wavelength or the thickness because the phase factor does not affect the intensity in Eq. (3). The polarization behavior of each A1 mode except P5 does not vary much

10

with thickness for the 632.8-nm excitation but shows significant variation for the 514.5-nm excitation. This is somewhat different from the case of black phosphorus, in which a shorter-wavelength excitation showed less variation with thickness [33]. An exception is P9, which appears to be an A2 mode for 632.8 nm but an A1 mode for 514.5 nm. Because P9 is relatively weak for the 514.5-nm excitation, its deconvolution from P10 is less reliable than for the 632.8-nm excitation. We tentatively assign P9 as an A2 mode which is also consistent with theoretical calculations [2,21,26,28].

3. Conclusions By comparing TED measurements with polarization of the Raman spectrum, we found a reliable method to determine the thickness and crystallographic orientation of few-layer WTe2. The polarization dependence of A1 modes, P2, P10, and P11, measured with the 632.8-nm excitation can be used to determine the crystallographic orientation of 1 to 4-layer and bulk WTe2. On the other hand, the polarization dependence changes with thickness when the 514.5 nm excitation is used. We also propose to use the position difference between P3 and P2, measured with the 632.8-nm excitation, for determination of the thickness, because P3 is enhanced for this excitation wavelength. The difference in the polarization dependence of the Raman peak intensities for different excitation wavelengths is explained in terms of energy-dependent optical dichroism due to the anisotropic structure of WTe2.

4. Methods 4.1. Sample preparation and Raman measurements WTe2 samples were mechanically exfoliated on p-type silicon substrates covered with 90-nm thick silicon dioxide. Because WTe2 can get degraded in air [8,21,23] (Supplementary Figure 11), we placed the samples in an optical chamber with ultrapure N2 environment immediately after exfoliation. After suitable

11

thin WTe2 samples were identified using an optical microscope, the substrate with WTe 2 samples was transferred to a micro vacuum chamber for Raman measurements. All Raman measurements were carried out in a vacuum of 10–5 Torr in order to minimize degradation. A long-working-distance objective (40×, N.A. 0.6) with a correction ring was used to focus the laser beam of 514.5-nm and 632.8-nm wavelengths to a spot size of about 1 m. The Raman signal was dispersed by a spectrometer with a focal length of 550 mm with a grating with 2400 grooves/mm and detected using a charge-coupled device (CCD). For low frequency Raman measurements in the range below 100 cm‒1, volume holographic filters were used to clean the laser lines and to eliminate the Rayleigh-scattered light (Supplementary Figure 12). A halfwave plate and a polarizer were used to control the polarization of the excitation laser and another halfwave plate was used to compensate for the polarization dependent throughput of the spectrometer (Supplementary Figure 12). Raman scattered light with the same polarization as the excitation laser (parallel polarization) was measured. The laser power was kept at 300 W to avoid laser-induced damages to the sample. The slight day-to-day fluctuations in the optical alignment, etc., are corrected for by normalizing against the intensity of the Raman peak from a piece of Si, by measuring the same piece in the same polarization. The number of layers of WTe2 was determined by comparing the optical contrast and atomic force microscopy (AFM) measurements after all the Raman measurements were completed (Supplementary Figure 13). The inter-layer spacing of bulk WTe2 is about 0.7 nm, but the height of monolayer WTe2 measured by AFM is about 1 nm, similar to other TMD materials.

4.2. Transmission electron diffraction (TED) measurements For TED measurements, exfoliated WTe2 flakes on the SiO2/Si substrate were transferred onto a TEM grid by a wet direct transfer method. In detail, for comparison between Raman and TED data, TED analysis should be performed at the particular flakes characterized by Raman spectroscopy. The wet direct transfer method without polymer coating is adopted for deterministic transfer because this method

12

is capable of fixing the grid on a particular flake through optical microscopy. To bond WTe 2 flake onto the amorphous carbon film on TEM grids, isopropyl alcohol (IPA) was dropped on top of the grid. As IPA evaporates, surface tension induces tight bonding between amorphous carbon and the WTe 2 flake [37]. Then, for detaching the grid bonded with the flake from the SiO2/Si substrate, a KOH solution was used for etching away SiO2. In order to minimize damages on WTe2, we soaked the grid/flake/SiO2/Si stack in the KOH solution to etch the surface of SiO2 rapidly. After a few seconds, the grid with the WTe2 flake was detached from the SiO2/Si substrate. Rinsing the KOH residue was the final stage for the preparation of TEM specimens. Electron diffraction pattern analysis was carried out using a transmission electron microscope, FEI Titan3 G2 60-300, operated at 80 kV acceleration voltage.

Acknowledgements This work was supported by a grant from the Center for Advanced Soft Electronics under the Global Frontier Research Program of the MSIP (No. 2011-0031630) and the National Research Foundation of Korea grants funded by the Korea government (MSIP) (No. 2015R1A2A2A01006992).

13

References [1]

Ali M N, Xiong J, Flynn S, Tao J, Gibson Q D, Schoop L M, Liang T, Haldolaarachchige N, Hirschberger M, Ong N P and Cava R J 2014 Large, non-saturating magnetoresistance in WTe2 Nature 514 205–8

[2]

Kong W-D, Wu S-F, Richard P, Lian C-S, Wang J-T, Yang C-L, Shi Y-G and Ding H 2015 Raman scattering investigation of large positive magnetoresistance material WTe2 Appl. Phys. Lett. 106 081906

[3]

Zhu Z, Lin X, Liu J, Fauqué B, Tao Q, Yang C, Shi Y and Behnia K 2015 Quantum oscillations, thermoelectric coefficients, and the Fermi surface of semimetallic WTe2 Phys. Rev. Lett. 114 176601

[4]

Lv H Y, Lu W J, Shao D F, Liu Y, Tan S G and Sun Y P 2015 Perfect charge compensation in WTe 2 for the extraordinary magnetoresistance: From bulk to monolayer Europhys. Lett. 110 37004

[5]

Zhao Y, Liu H, Yan J, An W, Liu J, Zhang X, Wang H, Liu Y, Jiang H, Li Q, Wang Y, Li X-Z, Mandrus D, Xie X C, Pan M and Wang J 2015 Anisotropic magnetotransport and exotic longitudinal linear magnetoresistance in WTe2 crystals Phys. Rev. B 92 041104(R)

[6]

Thoutam L R, Wang Y L, Xiao Z L, Das S, Luican-Mayer A, Divan R, Crabtree G W and Kwok W K 2015 Temperature-dependent three-dimensional anisotropy of the magnetoresistance in WTe2 Phys. Rev. Lett. 115 046602

[7]

Luo Y, Li H, Dai Y M, Miao H, Shi Y G, Ding H, Taylor A J, Yarotski D A, Prasankumar R P and Thompson J D 2015 Hall effect in the extremely large magnetoresistance semimetal WTe2 Appl. Phys. Lett. 107 182411

[8]

Wang L, Gutiérrez-Lezama I, Barreteau C, Ubrig N, Giannini E and Morpurgo A F 2015 Tuning magnetotransport in a compensated semimetal at the atomic scale Nat. Commun. 6 8892

[9]

Wang Y, Wang K, Reutt-Robey J, Paglione J and Fuhrer M S 2016 Breakdown of compensation and persistence of nonsaturating magnetoresistance in gated WTe2 thin flakes Phys. Rev. B 93 121108(R)

[10]

Yoon D, Moon H, Cheong H, Choi J S, Choi J A and Park B H 2009 Variations in the Raman spectrum as a function of the number of graphene layers J. Korean Phys. Soc. 55 1299–303

[11]

Lee C, Yan H, Brus L E, Heinz T F, Hone J and Ryu S 2010 Anomalous lattice vibrations of single- and few-layer MoS2 ACS Nano 4 2695–700

[12]

Yoon D, Son Y-W and Cheong H 2011 Strain-dependent splitting of the double-resonance Raman scattering band in graphene Phys. Rev. Lett. 106 155502

[13]

Li H, Zhang Q, Yap C C R, Tay B K, Edwin T H T, Olivier A and Baillargeat D 2012 From bulk to monolayer MoS2: Evolution of Raman scattering Adv. Funct. Mater. 22 1385–90

[14]

Zhao Y, Luo X, Li H, Zhang J, Araujo P T, Gan C K, Wu J, Zhang H, Quek S Y, Dresselhaus M S and Xiong Q 2013 Interlayer breathing and shear modes in few-trilayer MoS2 and WSe2 Nano Lett. 13 1007–15

[15]

Corro E d., Terrones H, Elias A, Fantini C, Feng S, Nguyen M A, Mallouk T E, Terrones M and Pimenta M A 2014 Excited excitonic states in 1L, 2L, 3L, and bulk WSe2 observed by resonant Raman spectroscopy ACS Nano 8 9629–35

[16]

Nguyen T A, Lee J-U, Yoon D and Cheong H 2014 Excitation energy dependent Raman signatures of ABA- and ABC-stacked few-layer graphene. Sci. Rep. 4 4630

14

[17]

Wolverson D, Crampin S, Kazemi A S, Ilie A and Bending S J 2014 Raman spectra of monolayer, few-layer, and bulk ReSe2 : An anisotropic layered semiconductor ACS Nano 8 11154–64

[18]

Lee J-U, Park J, Son Y-W and Cheong H 2015 Anomalous excitonic resonance Raman effects in few-layered MoS2 Nanoscale 7 3229–36

[19]

Zhang X, Qiao X-F, Shi W, Wu J-B, Jiang D-S and Tan P-H 2015 Phonon and Raman scattering of two-dimensional transition metal dichalcogenides from monolayer, multilayer to bulk material Chem. Soc. Rev. 44 2757–85

[20]

Guo H, Yang T, Yamamoto M, Zhou L, Ishikawa R, Ueno K, Tsukagoshi K, Zhang Z, Dresselhaus M S and Saito R 2015 Double resonance Raman modes in monolayer and few-layer MoTe2 Phys. Rev. B 91 205415

[21]

Lee C-H, Silva E C, Calderin L, Nguyen M A T, Hollander M J, Bersch B, Mallouk T E and Robinson J A 2015 Tungsten Ditelluride: a layered semimetal Sci. Rep. 5 10013

[22]

Jana M K, Singh A, Late D J, Rajamathi C R, Biswas K, Felser C, Waghmare U V. and Rao C N R 2015 A combined experimental and theoretical study of the structural, electronic and vibrational properties of bulk and few-layer Td-WTe2 J. Phys. Condens. Matter 27 285401

[23]

Kim Y, Jhon Y I, Park J, Kim J H, Lee S and Jhon Y M 2015 Anomalous Raman scattering and lattice dynamics in mono- and few-layer WTe2 Nanoscale 8 2309–16

[24]

Lee J-U, Kim K, Han S, Ryu G H, Lee Z and Cheong H 2016 Raman signatures of polytypism in molybdenum disulfide ACS Nano 10 1948–53

[25]

O’Brien M, McEvoy N, Hanlon D, Hallam T, Coleman J N and Duesberg G S 2016 Mapping of low-frequency Raman modes in CVD-grown transition metal dichalcogenides: Layer number, stacking orientation and resonant effects Sci. Rep. 6 19476

[26]

Jiang Y C, Gao J and Wang L 2016 Raman fingerprint for semi-metal WTe2 evolving from bulk to monolayer Sci. Rep. 6 19624

[27]

Lee J, Ye F, Wang Z, Yang R, Hu J, Mao Z, Wei J and Feng P X-L 2016 Single- and few-layer WTe2 and their suspended nanostructures: Raman signatures and nanomechanical resonances Nanoscale 8 7854–60

[28]

Jana M K, Singh A, Late D J, Rajamathi C R, Biswas K, Felser C, Waghmare U V. and Rao C N R 2015 A combined experimental and theoretical study of the structural, electronic and vibrational properties of bulk and few-layer Td-WTe2 J. Phys. Condens. Matter 27 285401

[29]

Dawson W G and Bullett D W 1987 Electronic structure and crystallography of MoTe2 and WTe2 J. Phys. C Solid State Phys. 20 6159–74

[30]

Aroyo M I, Perez-Mato J M, Capillas C, Kroumova E, Ivantchev S, Madariaga G, Kirov A and Wondratschek H 2006 Bilbao Crystallographic Server: I. Databases and crystallographic computing programs Zeitschrift fur Krist. 221 15–27

[31]

Strach T, Brunen J, Lederle B, Zegenhagen J and Cardona M 1998 Determination of the phase difference between the Raman tensor elements of the A1g-like phonons in SmBa2Cu3O7-δ Phys. Rev. B 57 1292–7

[32]

Ribeiro H B, Pimenta M A, De Matos C J S, Moreira R L, Rodin A S, Zapata J D, De Souza E A T and Castro Neto A H 2015 Unusual angular dependence of the Raman response in black phosphorus ACS Nano 9 4270–6

[33]

Kim J, Lee J-U, Lee J, Park H J, Lee Z, Lee C and Cheong H 2015 Anomalous polarization dependence of Raman scattering and crystallographic orientation of black phosphorus Nanoscale 7 18708–15

15

[34]

Mao N, Wu J, Han B, Lin J, Tong L and Zhang J 2016 Birefringence-directed Raman selection rules in 2D black phosphorus crystals Small 12 2627–33

[35]

Kranert C, Sturm C, Schmidt-Grund R and Grundmann M 2016 Raman tensor formalism for optically anisotropic crystals Phys. Rev. Lett. 116 127401

[36]

Ling X, Huang S, Hasdeo E H, Liang L, Parkin W M, Tatsumi Y, Nugraha A R T, Puretzky A A, Das P M, Sumpter B G, Geohegan D B, Kong J, Saito R, Drndic M, Meunier V and Dresselhaus M S 2016 Anisotropic electron-photon and electron-phonon interactions in black phosphorus Nano Lett. 16 2260–7

[37]

Lin W-H, Chen T-H, Chang J-K, Taur J-I, Lo Y-Y, Lee W-L, Chang C-S, Su W-B and Wu C-I 2014 A direct and polymer-free method for transferring graphene grown by chemical vapor deposition to any substrate ACS Nano 8 1784–91

16