Growth of ‘W’ doped molybdenum disulfide on graphene transferred molybdenum substrate

In the present study, a novel method has been carried out to grow tungsten (W) doped molybdenum disulfide (MoS2) on the graphene transferred TEM grid in a chemical vapor deposition (CVD) setup. Tungsten trioxide (WO3) has been used as a source for ‘W’ while ‘Mo’ has been derived from Mo based substrate. Different experimental parameters were used in this experiment. Higher gas flow rate decreases the size of the sample flake and on other side increases the dopant concentrations. The interaction mechanism between Mo, S, W and oxygen (O) have been explored. The influence of oxygen seems to be not avoidable completely which also imposes effective growth condition for the reaction of Mo with incoming sulfur atoms. The difference in the migration energies of Mo, WO3, S clusters on the graphene and the higher reactivity of Mo clusters over other possibly formed atomic clusters on the graphene leads to the growth of W doped MoS2 monolayers. Formation of MoS2 monolayer and the nature of edge doping of ‘W’ is explained well with the crystal model using underlying nucleation principles. We believe our result provide a special route to prepare W doped MoS2 on graphene substrate in the future.

Previous researches have been conducted using Mo solid thin films prepared by thermal evaporation 29,30 , electron beam evaporation 31,32 , sputtering 33 and then sulfurization continued either in the second step of experiment or in a specially designed system 33 , in all the cases the resulting MoS 2 layers were dependent on the thickness of deposited thin films and so far, produced only MoS 2 monolayers without any possible doping. In this work, the sulfurization of Mo and then doping of W into MoS 2 on the graphene have been carried out in a single and simple CVD set-up (refer Fig. 1a). We have used Mo based substrate as a source for 'Mo' and WO 3 powder as a source for 'W' . MoS 2 monolayers are expected to grow first in this process, due to the higher chance for Mo to react with incoming 'S' atoms compared to 'W' atoms, and hence, 'W' atoms reach after the formation of MoS 2 monolayer will find its place as dopant atoms at Mo lattice site. Direct observation of such alloy structure at the atomic scale is essential to understand and to further control the wide range of its behavior, and hence, graphene has been chosen as a substrate and placed on TEM grid to prevent the loss of small atomic clusters during the traditional transfer process when it grows on SiO 2 , or sapphire or mica 34 . In addition, the formation of graphene and transition metal doped heterostructures is very interesting, owing to the combination of high mobility and transparency of graphene and the exceptional optical properties of MoS 2 18,35 .

Results and Discussions
To verify the possibilities of 'Mo' to be derived from Mo based substrate and to evaluate the interaction mechanism of Mo with S and Mo with W, the experiments with different parameters (refer Table 1) were conducted using typical CVD set-up shown in Fig. 1a. A typical low magnification annular dark-field (ADF) scanning transmission electron microscope (STEM) images of MWS-1a, MWS-1b and MWS-1c are displayed in Fig. 1(b-d), respectively. The flakes in MWS-1a (Fig. 1b) are identified to grown at a scale of 70-200 nm, and most of them are in uniform orientation.
MWS1-b (Fig. 1c) contains the flakes around 75-250 nm. MWS-1b contains nanostrips and other polygons, while MWS-1a contains mostly with triangles while other polygons are rarely found. Transformation of triangular, hexagonal flakes to a layer with different morphologies through the collision of neighboring structures usually observed in 2D materials growth 25,36 .
MWS-1c (Fig. 1d) shows it contains triangles in a large distribution, while almost formed in multi-layer or layer over layer growth. The flakes are grown at a scale of 50-80 nm, much smaller compared to MWS1-a. Flow rate of carrier gas, Ar 400 sccm, reduces the growth propagation while increases the nucleation density (Fig. 1d), compared to the flow rate, Ar 300 sccm (Fig. 1b). The growth temperature 875 °C seems to impose the best condition for the evolution of triangles, however, the argon flow rate decides the high-quality growth.
Growth of monolayers on graphene also seems to be influenced by the underlying substrate and on the graphene synthesis techniques 37 . In all our present experiments, the underlying substrate is TEM grids and the graphene for all these experiments were transferred to TEM grids grown on the same copper foil, and hence there were no variation of influence based on the substrate principles.
At relatively lower temperature, that is at 825 °C, for the sample MWS-1b (Fig. 1c), the sulfurization in the initial period is difficult, leading to the variation of Mo:S ratio and the MoS 2 becomes unsaturated 38 . Mo-terminated zig-zag (Mo-zz) will always grow faster than S-terminated zig-zag (S-zz) [39][40][41] , and the relatively lower temperature would increases the differences in the rate of edge growth, and leads to the other morphologies 38 than the triangles and hexagons, in this case its nanostrips (Fig. 1c). At the same time, with the increase in flow rate of carrier gas, S-concentration is relatively higher 40 , and with the variation of Mo:S ratio, with S-rich environment, the rate of growth between two edges, Mo-zz and S-zz would be small, results in a small triangle with truncated shape 39,41 .
On an average, 20-25 flakes have been counted to every 5 μm circumference region of TEM grid in MWS-1a samples, while MWS-1b sample contains less than 10 flakes, rest is filled with polygons and MWS1-c sample contains more than 75-100 flakes, however identifying individual flakes are hard as it contains mostly a layer over layer growth.
An atomic-scale resolution ADF-STEM images of MWS-1a, MWS-1b and MWS-1c are displayed in Fig. 2(a-c). All three images show the distribution of W doping within a monolayer, effective doping around the corners of a flake, while the interior part without any effective doping. The experimental image (Fig. 2d) of a selected region of W doped MoS 2 structure within MWS-1a (Fig. 2a). And the corresponding line profile together with simulation data (Fig. 2(e-g)). The intensity of line profile in Fig. 2d matches with the line profile in Fig. 2f confirms that the W atoms are doped into the MoS 2 monolayer. During the STEM characterization, it is already evaluated that there is no displacement of atoms of higher intensity confirms those atoms are alloy atoms and not probably 'Mo' or 'W' add atoms. QSTEM simulation 42 have also been employed to validate it. Intensity profile of marked region (green arrow) in Fig. 2d illustrating MWS1-a sample matches well with Fig. 2g intensity profile, confirms that those bright atoms are W alloy into MoS 2 monolayers.
An atomic-scale resolution ADF-STEM images of MWS-1b is displayed in Fig. 3a and an approximate area within green rectangular and blue square regions marked in Fig. 3a are illustrated in Fig. 3b and c respectively. Inset in Fig. 3b   Atomic-scale ADF-STEM images and crystal model of some small atomic clusters of MoS 2 (MWS-1c) are displayed in Fig. 4(a,b). Of more than 30 samples we checked, most samples in small size are MoS 2 monolayer clusters without W atom doping. Only one W atom was found which was marked in Fig. 4a (Arrowed). This result confirms that in this growth process, MoS 2 was prepared first. The formation of MoS 2 nanoflakes and monolayers are also possible in the growth temperature range of 700-750 °C, as evident from the reported elsewhere 13,41 , we suspect these small clusters might be formed during the cooling stage, at which point, Mo and S are still energetically active species to react and grow further, while W atoms does not have a chance to find its energetically active position in the already formed MoS 2 structures.
To further elucidate the contribution of Mo atoms (clusters) derived from Mo-based substrate into the reaction and to discuss the role of oxygen in the monolayer formation mechanism, low magnification ADF-STEM images and Energy dispersive X-ray spectroscopy (EDS) mapping of as-synthesized MoS (Fig. S1) and MoWO (Fig. S2). MoS sample possess irregular and small sized particles (Fig. S1a), while MoWO sample consists mostly of regular shaped bright particles (Fig. S2a). Sulfurization of Mo/TEM grid without WO 3 precursor does not produce any MoS 2 monolayer flakes, suggests partial pressure environment formed by the addition of W-O (x) /WS x O (3−x) is necessary for optimum experimental condition to grow MoS 2 . EDS mappings (refer Fig. S2(b-e) for MoS and Fig.  S2(b-e) for MoWO) shows the common trend in both the samples, Mo signals are stronger in the core. S and W signals for respective MoS and MoWO samples are detected in the larger area, extended to outer shell.
According to these results, we could discuss the microscopic growth process of W doped MoS 2 (Fig. 5). The availability of 'O' may play an important role during this special CVD growth process. Oxygen (O) in these reactions are probably derived from the residue happened to form during chemical-wet transfer of graphene into TEM grid from copper substrate, residual atmosphere and the oxidation of molybdenum. In the traditional CVD process, the excess oxygen may cause the etching of MoS 2 43 and play a major role in the particles or fullerene structure formation as reported earlier in few works 41,44 . In this process, more oxygen is required to form MoO x clusters first. After sulfurization (in the case of MoS sample), Mo clusters possibly forms an amorphous alloy which possess octahedral or tetrahedral co-ordination 45,46 . Due to this complexity in its structure, it requires large energy to change its geometry to evolve as two-dimensional hexagonal network 45,47 . Speaking in terms of partial pressure environments, higher S partial pressure environment would lead to the fullerene structure, lower S partial pressure leads to oxy-sulfide particles 48 . In this case, promoter such as WO 3 in its reduced form W-O (x) /WS x O (3−x) probably needed to reduce the complex structure as in the case of MoS and provides suitable S partial pressure environment to form two-dimensional hexagonal network of MoS 2 .
In step 2, at the high temperatures between 825-875 °C, Mo get oxidized and forms thermally active MoO 3 structures together with Mo clusters which is not oxidized. Due to the differences in the surface free energy of graphene (around 46.7 mJ/m 2 ) 49 and Mo (around 2.92-3.34 J/m 2 ) 50 , metal atoms migrating on graphene is expected to form molecular clusters on graphene 51,52 , which is more reactive with sulfur and oxygen than the Mo atoms and W-O (x) clusters, and hence readily reacts with the incoming sulfur atoms (Step 3) 31   In summary, W-doped MoS 2 monolayer structures on graphene have been synthesized using Mo atoms (clusters) derived from Mo based substrate. Comparing to the growth temperatures, 875 °C and 825 °C, the former case provides high-quality and uniform triangular MoS 2 monolayers with less concentrated W doping, while in the latter case, even though the concentration of W doping among the individual triangular monolayer structures is similar as in the case of former one, merged monolayer structures with irregular shape have been found in majority. Nucleation density is large at lower temperature, however, the probability of finding individual triangles is lesser. The size of flake grown in lower temperature experiment is also larger. Variation in carrier gas flow conditions (400 sccm) have been done to verify the tuning possibilities of dopant metal atoms into the MoS 2 monolayers and found that the dopant atoms are distributed up to 60 atomic rows, while in other experiments with carrier gas flow condition of 300 sccm, the dopant atoms distribution is up to maximum of 40 atomic rows (in both temperatures, 875 °C and 825 °C). Higher gas flow rate also decreases the size of the flakes and on other side increased the dopant concentrations. The interaction mechanism between Mo, S, W and O on the graphene and the formation of MoS 2 monolayer and the nature of edge doping of 'W' have been briefly discussed. The influence of oxygen seems to be not avoidable completely and also impose effective growth condition for the reaction of Mo with incoming sulfur atoms. Doping is effective in all the growth conditions and in all structures irrespective of the morphologies, from small clusters to nanostrips indicate the successful nature of doping mechanism. This method of using graphene transferred TEM grid substrate could give us a new idea to prepare TMDs alloys (Re, Ta) in the future.

Methods
Sample preparation. Graphene films grown on polycrystalline copper foils were transferred onto molybdenum (Mo) based TEM grids using a standard PMMA-assisted wet-chemistry transfer process. The TEM grids with graphene facing downwards were mounted onto a home-built ceramic carrier and 1 mg of Tungsten trioxide (WO 3 ) precursors (Sigma-Aldrich, 99%) placed on the silicon wafer. Quartz boat carrying silicon wafer and TEM grid at its center were loaded into the CVD chamber.
Within a typical CVD process, the furnace was firstly heated to 300 °C in 10 mins and hold for additional 10 mins, and then increased to its growth temperature, 825 °C/875 °C in 45 mins. Just 2 mins before the main CVD furnace reaches its growth temperature, 300 mg of sulfur (S) source (Aladdin, 99.999%) were heated (reaches 250 °C in 2 mins) from the separate heating belt in upstream. The whole system was maintained at their respective temperatures for 15-20 mins (growth time) and then cooled down to room temperature naturally. Argon (Ar) flow was maintained throughout an experiment at the rate of 300-400 sccm. The experimental conditions to prepare set of samples is shown in Table 1. Here, typically followed a process as did for MWS-1, and carried out a sulfurization for 15 mins without using WO 3 source material and the sample is named MoS and another experiment conducted using WO 3 alone without sulfurization, that is, after the main temperature reaching 875 °C and the system was kept for 15 mins and naturally cooled down and the sample is named MoWO. STEM characterization. Low magnification annular dark field -scanning transmission electron microscopy (ADF-STEM) images were first obtained with FEI Technai F20 to have a preliminary check over the samples and high resolution atomic-scale ADF-STEM images were captured using FEI Chemi-STEM Titan G 2 which was equipped with a probe-side spherical aberration-corrector and operated at an acceleration voltage of 200 kV. The convergent semi-angle for illumination was set to 24 mrads with a probe current of 50-70 pA, and the collection angle was 50-100 mrads. Energy dispersive X-ray spectroscopy (EDS) was carried out on a Bruker super-X detection system. Data availability. The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request. Table 1. Parameters used for the preparation of MWS-1 samples are given here. *No WO 3 is used, only sulfurization of graphene/Mo-TEM grid was carried out. $ No sulfur is used. Only WO 3 was heated over which graphene/Mo-TEM grid placed face down.