Caveats in obtaining high-quality 2D materials and property characterization

<jats:p><jats:fig position="anchor"><jats:graphic xmlns:xlink="http://www.w3.org/1999/xlink" mime-subtype="png" mimetype="image" xlink:href="S0884291419003947_figAb.png" /></jats:fig></jats:p>


Introduction
In 2004, the state-of-the-art technique for realizing atomically thin materials consisted of a roll of tape and a layered crystal held together by van der Waals forces [1,2]. The so-called scotch-tape exfoliation method enabled isolation of molecularly thin layers from dozens of van der Waals [also known as twodimensional (2D)] materials such as molybdenum disulfide (MoS 2 ), tungsten diselenide (WSe 2 ), and hexagonal boron nitride (hBN) [2,3]. On realizing the promise of 2D materials (2DMs) for use in next-generation technologies, the community where he focused on carbon nanotube synthesis and devices. During this time, the fever that is 2DMs was reborn, and Dr. Robinson began working on epitaxial graphene growth of silicon carbide. Subsequently, in 2007, Dr. Robinson moved back to Penn State to work for the Penn State Electro-Optics Center where he established projects in graphene and "beyond graphene" materials. The success of these programs enabled him to transition to the Penn State Materials Science and Engineering Department as an Assistant Professor and was promoted to an Associate Professor in 2015. He co-founded the Center for Two-Dimensional and Layered Materials and the NSF I/UCRC Center for Atomically Thin Multifunctional Coatings (ATOMIC) and currently serves as Co-Director of both. He is a 2015 NSF-CAREER recipient and has authored or coauthored over 200 peer reviewed journal publications with a significant focus on low dimensional electronic materials.
ª Materials Research Society 2020. This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
quickly advanced beyond micrometer-sized flakes of 2D layers produced by the scotch-tape method and toward the realization of continuous films over large-area substrates (up to 300 mm diameter SiO 2 /silicon) [4,5,6,7,8,9], all enabled by simple vapor deposition techniques [10]. The most rudimentary of these techniques utilizes a tube furnace-based approach, including solid source precursors and, more recently, nucleation or seeding agents that can yield 2D films over large areas [4,11,12,13]. The beauty of this approach lies in its logistical simplicityone can merely heat source powders and substrates together in a tube furnace to achieve rapid synthesis of large-area 2D films in nearly any laboratory. This enables researchers to investigate fundamental material properties, basic physics, and proof-ofconcept devices with only a small barrier to creating the films themselves [10]. Following deposition, 2D films are commonly characterized via Raman, photoluminescence (PL), and X-ray photoelectron spectroscopies, which can be used to determine film composition, layer number, and band gap. Despite the ease of these deposition and characterization techniques, some limitations exist in their application to 2DMs [10]. Here, we highlight commonly overlooked scientific details of 2DMs synthesis and characterization that can ultimately impact scientific claims, with focus on solid-source and seed-promoted deposition, as well as Raman, PL, and X-ray photoelectron spectroscopic characterization.

Solid-source chemical vapor deposition
The rapid development of the field of 2DMs synthesis research has been possible due to its relatively low cost of entry. For only a few thousand dollars, solid source precursors such as sulfur, tungsten-oxide, and molybdenum-oxide powders, and commonly available inventory, one can set up a simple tube furnace and begin growing 2DMs. While this tube furnace infrastructure dominates 2DM synthesis research today [14,15,16], realization of truly repeatable, high-quality synthetic 2DMs requires significant modifications to solid-source (S) chemical vapor deposition (CVD) processes or completely different syntheses approaches. This is due to the poor reproducibility of S-CVD synthesis recipes, which often vary widely between furnaces and even between individual researchers in a single laboratory. Several variables of paramount importance can easily be overlooked during S-CVD processes, such as reactor and crucible geometry, furnace heating and cooling rates [17,18], and substrate cleaning or surface preparation [5]. These variables are often unreported or overlooked in the literature, with many researchers claiming S-CVD to be a viable technique for uniform, large-scale 2DM synthesis. While some studies demonstrate quasi-uniform deposition on centimeter-scale samples, scaling to industrially relevant substrate sizes (.200 mm) using S-CVD is unlikely. Similarly, 2DM literature contains little mention of reactor "seasoning" or the effects of precursor buildup on crucibles and the interior of quartz tubes when S-CVD is utilized. This buildup significantly impacts material growth in hot wall S-CVD studies because it leads to a continuously changing partial pressure of precursors with each additional run, making every run unique [19]. Beyond incomplete reporting of reactor and process details in 2DM S-CVD literature, S-CVD offers little control over precursor materials or concentration at the substrate surface. One commonly observed consequence of the horizontal tube furnace-based S-CVD approach is a deposition plume [Figs. 1(a) and 1(b)], containing a gradient of 2D layer domain morphologies across the substrate surface. The quintessential triangular monolayer is often only observed only in a specific region of this deposition plume [ Fig. 1(c)], due to variation in the precursor concentration across the substrate surface, and as a result, varying 2D crystal edge terminations and growth rates [20]. Control over precursor concentration in S-CVD studies is primarily achieved by varying source temperature and distance from substrate [20].
While S-CVD is a common technique for rapid realization of millimeter-scale 2DMs, the complexity of tube furnacebased S-CVD can lead to poor reproducibility even between personnel within a single research group. Therefore, process parameters reported for S-CVD techniques should be viewed merely as qualitative guidelines, where reported recipes are not likely to directly translate between laboratories but can serve as a valuable starting point for 2DM synthesis studies. One effort to reduce variability in S-CVD of 2DMs lies in computational investigations to understand flow patterns and precursor concentration variation at substrates surfaces [17,18,21]. While these efforts are relatively few, they can greatly aid understanding of the requirements for 2DM growth. Beyond S-CVD growth of 2DMs, it is important to consider the use of alternative gas source deposition methods, such as atomic layer deposition [22], traditional CVD [4], and metal organic (MO) CVD [10], which offer controllable introduction of precursors into the reactors. Such techniques may require greater initial capital but can yield highly reproducible results and uniform, large-area 2D films. Efforts to achieve large-area films of 2DMs have matured in recent years, with demonstrations of coalesced films of transition metal dichalcogenides (TMDs) like tungsten diselenide (WSe 2 ) and MoS 2 grown over centimeter scales [4,5,6,7]. These large-area films have been possible through the use of MOCVD, where the concentration and ratios of metal and chalcogen species can be precisely controlled via bubblers, and sequential nucleation and ripening steps may be performed to control TMD domain sizes and layer number [5,6].
Beyond the use of gas source deposition methods, substrate pretreatment can be used to control 2DM growth. This is shown in both MOCVD and S-CVD [Figs. 1(d)-1(g)] studies of TMD growth on c-plane sapphire [20]. While many 2DM synthesis studies utilize SiO 2 /Si substrates, c-plane sapphire is an increasingly popular substrate, due to potential for aligned 2DM domains resulting from long-range commensurability with sapphire [5,6,7,23]. Direct growth of 2DMs on asreceived c-plane sapphire, however, may not be sufficient in achieving oriented 2DM domains from S-CVD growth. Annealing sapphire to reduce the contact angle of the substrates is shown to enable dramatic improvement in domain orientation due to the formation of terraces in the sapphire surface, which promote certain domain orientations [5,23,24]. The use of c-plane sapphire as a substrate for epitaxial growth has enabled understanding of the requirements for lateral growth of domains and control of nominal layer number [6]. Observation of this domain orientation in both S-CVD and MOCVD studies reveals that substrate engineering can be an effective strategy to control the orientation and growth of 2D layers.

Seed-promoted CVD
In addition to the use of MOCVD for large-area 2DM growth, nucleation agents and seeding techniques have gained traction as tools to achieve increased 2DM domain sizes [4,11,12,13,25]. Kang et al. [4] popularized the use of alkali-based compounds (e.g., NaCl) to achieve large-domain 2D TMD layers via MOCVD, and the approach has since been expanded to a large library of 2D layers grown via S-CVD [26]. While salt-assisted growth does result in large-area TMD films with fewer grain boundaries, recent studies show Na atoms present at the TMD/substrate interface, which can lead to variation in coupling between the growth substrate and TMD layers. This causes spatial variation in electronic properties of the grown TMD [27]. Such variations in coupling can also affect the quality of the grown material significantly, resulting in a tensile strain in the layers and a reduction of the PL intensity of MoS 2 layers by roughly 100Â. Furthermore, these studies indicate that while NaCl may serve to control nucleation and promote large area growth of MoS 2 , NaCl can also degrade the electrical performance of these grown semiconducting layers even when the films are transferred from the growth substrate [27]. More importantly, however, is the 501 years of research in the semiconductor industry that provides overwhelming evidence that the presence of alkali ions leads to higher oxide leakage currents and reliability problems in semiconductor devices when operated at high temperatures and voltages [28,29,30,31,32,33].
This was demonstrated by Snow et al. [32] in 1965 (Fig. 2) through nonuniform shifts in capacitance-voltage (C-V) curves in a Na-contaminated metal-oxide-silicon (MOS) capacitor. The oxide leakage currents observed in the MOS structures are associated with the drift in C-V curves. The ionic  drift within the oxide leads to significant shifts in the silicon charge, where upon application of a positive bias at 150°C for 5 min, a .100 V bias shift in the C-V curve occurs [ Fig. 2(a)], corresponding to a buildup of negative charge. This shift is independent of substrate doping. When Na and other alkali ions are not present in the MOS capacitors, the C-V curves do not drift. As shown in Fig. 2(b), alkali-contaminated devices can also exhibit more complicated C-V curves after drift. These distorted curves are caused by nonuniform alkali contamination over the surface of the device, causing some areas to drift more than others, leading to adjacent devices displaying different characteristics. Considering many integrated circuits operate between 55 and 85°C, one can expect significant diffusion of alkali ions and reduced transistor stability if such ions are present, even at low concentrations. Any commercial semiconductor fabrication facility will strongly object to a process that utilizes any alkali ion-based process, indicating that alkali ions must be avoided at all cost when considering 2DMs for augmenting semiconductor technology for complimentary metal-oxide-semiconductor (CMOS) applications. Other applications including low power devices for Internet of Things technology and 2DM-based chemical and biological sensors may have a higher tolerance for these mobile ions.

Raman and PL spectroscopies
Detailed investigation of the quality of synthesized 2DMs requires a technique that examines the sample on an appropriate length scale and provides relevant structural or chemical information. Raman and PL spectroscopies are excellent characterization tools for 2DMs because they require little to no sample preparation and provide a wealth of information within short acquisition times [34,35,36,37]. However, this extensive information can be easily misinterpreted if one is inexperienced in spectral analysis of inelastic scattering in 2DMs. One commonly misinterpreted result of spectral (Raman and PL) data relates to the identification of the number of layers in 2DMs. Generally, this misinterpretation is a result of limited spectral analysis and an unawareness that certain spectral features can be affected by factors unrelated to layer number (Fig. 3).
The most common method of reporting 2DM layer number is via Raman spectroscopy, based on the spectral separation of the E 1 2g and A 1g vibrational modes in MoS 2 [ Fig. 3(a)], which increases with increasing layer number [38]. This approach is challenging due to the factors such as film strain [ Fig. 3(b)] [39], concentration of sulfur vacancies [ Fig. 3(c)] [40], and electronic carrier concentrations [ Fig. 3(d)] [41]-all of which impact the absolute spacing between E 2g and A 1g in MoS 2 . Some of these characteristics like sulfur vacancy concentration can change with time or environment, whereas others like strain and carrier concentration are predominantly determined by the growth conditions. There are a few additional ways to utilize Raman for identifying 2D layer number, as noted in Fig. 3(e). Given access to a Raman spectrometer with ultra-low frequency Rayleigh filters and thus the ability to analyze THz frequencies (approximately spectral regions below 50 cm À1 ), one may use THz modes to confirm 2DM layer number. Raman peaks in this spectral range for van der Waals materials are associated with shear [C in Fig. 3 Fig. 3(e)] modes and are therefore the most direct way of counting the number of 2D layers [42]. In addition, the THz bands are less susceptible to peak shift as a function of strain [43].

(e)] and layer breathing [LB in
Another approach to layer number measurement-often more accessible than THz Raman-is the measurement of 2DM PL. Semiconducting monolayer TMDs typically exhibit an optical band gap in the visible spectrum. Measuring the PL of these materials requires the use of an excitation laser energy that is above the semiconductor optical band gap and a spectrometer that enables access to the PL energy [ Fig. 3(f)]. Most research-grade Raman spectrometers match these conditions for MoS 2 , which has an indirect band gap of ;1.3 eV in bulk and a direct band gap of ;1.8 eV at monolayer thickness. Therefore, the PL peak position of MoS 2 may be used as an additional indication of the layer number [44]. However, the PL response of the monolayer material is so sensitive that the peak rapidly and dynamically evolves in response to the environment-even when one merely breathes onto the sample (Fig. 4). Molecules adsorbed onto MoS 2 layers from breathing onto the surface change the electronic band structure of the MoS 2 . One breath onto the sample changes the PL peak position by 20 meV. This peak position change is greater than that between MoS 2 flakes with thicknesses of 1.3 and 7.6 nm, as shown in Fig. 3(f) [44]. This "breath" experiment highlights the importance of understanding sample history and the environmental considerations when utilizing PL for quantification of layer thickness. Conversely, many groups have investigated the PL response of 2DMs as chemical sensors, leveraging this surface sensitivity while considering irradiation effects from the excitation laser [45,46].

X-ray photoelectron spectroscopy
In addition to Raman and PL, X-ray photoelectron spectroscopy (XPS) can serve as a valuable 2DM characterization technique. XPS provides direct evidence of elemental bonding  at the surfaces of materials, making it attractive for understanding the composition of 2DMs. This technique has a unique combination of attributes that align particularly well with 2DMs, including surface sensitivity, accurate quantification with little or no matrix effects, sensitivity to every element in the periodic table except H and He, and the ability to determine the nearest neighbor chemical environment in many instances. This technique is commonly used in the 2DM field to investigate chemistry and stoichiometry in 2DMs, as well as dopant concentrations and impacts on valence band maxima (VBMs). To effectively utilize XPS for investigation of 2DMs, several considerations regarding stoichiometry and valence band measurements must be made. Perhaps the attribute that XPS is best known for is its ability to determine the nearest neighbor chemical environment. The binding energy of core electrons emitted and detected via XPS is affected by formal valence state and other differences in chemical environment. Thus, XPS can be used to explore surface and interfacial reactions between 2DMs and to investigate relative elemental composition. Figure 5 shows a typical spectrum of WSe 2 grown on sapphire (Al 2 O 3 ), where W, Se, Al, and O peaks are clearly evident. The inset shows high-resolution spectra of the W 4f and Se 3d regions, and the measured Se:W ratio is found to be 1.96 (close to the expected 2.0). As a check of repeatability and accuracy of these measurements, a series of bulk 2DMs were repeatedly exfoliated and measured for composition (for flakes with acceptably low levels of C and O contamination). The standard deviation of the chalcogen-to-metal ratio varied by ,2% for the TMDs MoS 2 , WS 2 , WSe 2 , and the layered material Bi 2 Te 3 . The actual standard deviation of thin 2DMs on substrates is expected to be worse due to lower signal strength but will remain ,5%.
The accuracy of these measured stoichiometries, however, is another matter entirely. In all cases, the bulk material composition and crystal chemistry were confirmed using Xray fluorescence and X-ray diffraction, respectively. The XPSmeasured stoichiometries of MoS 2 and WS 2 were within 3% of the true values. However, when this same methodology was applied to WSe 2 , the Se:W ratio was 1.69 instead of 2.00, a 16% error. The Te:Bi ratio on Bi 2 Te 3 was 1.19 instead of 1.50. These results, obtained under ideal conditions, should provide a cautionary tale that determining the absolute stoichiometry of 2DMs may require derivation of relative sensitivity factors (RSFs) from reference samples rather than relying on RSFs supplied by the XPS instrument vendor. Under the best conditions of a continuous thin film of MoS 2 on a substrate with no interfering peaks, the 2r standard deviation (95% confidence interval) of XPS-measured stoichiometry suggests that the best that can be achieved is MoS 2.060.1 (5%). Therefore,  one should avoid claims of a stoichiometric 2DM based on XPS, as the error is likely 3-5%-a variation large enough to dramatically affect the optical and electronic properties.
Beyond measuring the stoichiometry of 2DMs, there is intense interest in doping 2DMs to alter their electronic properties [10]. XPS is ideal for understanding the bonding environment of the dopants (e.g., dopant bonding to chalcogen or transition metal) but can also shed light on how the dopant impacts the electronic structure via measurement of VBMs. While the data used to identify bonding in XPS are based on core electrons, weakly bound valence electrons are also accessible. Measurements of the VBM can be combined with careful core level measurements to yield valence band offsets, which are important for understanding and predicting potential heterojunctions in 2D devices [47,48]. However, the X-ray cross section of these electrons is typically quite low, resulting in weak XPS signals [ Fig. 6(a)]. This is evident from the spectra acquired from V-doped WS 2 on SiO 2 [ Fig. 6(a)]. The sum of W, S, and V signals from this sample is ,7% of the total signal, and nearly 90% of the signal originates from the SiO 2 substrate (with the balance coming from adsorbed organics). Furthermore, the V 2p 3/2 spectra of undoped, lightly doped, and highly doped WS 2 [ Fig. 6(a), inset] show evidence of multiple chemical states, requiring curve fitting into three overlapping components despite the fact that there is ,0.3% total vanadium present. The peak at 513.5 eV is due to VS x species with the higher binding energy peaks due to various oxides of vanadium. The lightly doped sample exhibits a W:V of ;41, while the more heavily doped sample has a ratio of 7.7. In addition to these weak XPS signals, when films are placed on dielectric substrates, there is inevitable substrate charging that must be accounted for.
When determining VBMs in XPS, care must be taken to ensure that the acquired data are statistically significant. This is exemplified in Fig. 6(b), which shows measured VBMs for Redoped (W:Re 5 4) and undoped WSe 2 grown on sapphire. Freshly exfoliated bulk WSe 2 is also shown for comparison. This provides a prime example of how Se vacancies can also lead to shifting in the VBM, even in undoped samples, where there is a 0.3 eV shift between exfoliated WSe 2 grown via chemical vapor transport and MOCVD grown WSe 2 (undoped). The VBM is typically determined by finding the intersection of two linear fits of background near 0 eV and the near-linear portion of the leading edge of the valence band [ Fig. 6(b) inset]. Through this approach, a very small (0.05 eV) difference is found between the doped and undoped WSe 2 [ Fig. 6(b)]. However, repeated measurements (n 5 4) reveal that these values are not statistically different. By employing high-energy resolution scans with signal-to-noise ratio .30:1, careful charge correction and careful fitting, a 2r standard deviation 60.10 eV, can be achieved. Thus, researchers must take care to acquire sufficiently precise and statistically significant measurements when investigating doping and VBM in 2DMs as the noise levels in these measurements can be much greater than the shifts themselves.

Conclusion
This article accumulates commonly overlooked scientific details in the synthesis and characterization 2DMs, with a focus on TMDs. It is intended to be used as a cautionary guide for researchers working in vapor-phase synthesis and spectroscopic characterization of 2DMs. Specifically, researchers should be aware of variability associated with S-CVD-grown materials and potential influences of nucleation agents on the properties of the synthesized 2DMs. When characterizing grown 2DMs, an understanding of the impacts of growth conditions and the local environment is crucial. Furthermore, reporting of 2DM layer number and stoichiometry must be done with caution if Raman and XPS techniques are used, and statistical significance of VBM measurements should be reported. As the 2DM field continues to progress, synthesis and characterization details must be scrutinized and communicated to researchers in the field at large to ensure consensus in scientific reporting and drawn conclusions.