Literature DB >> 30081483

Feasible Route for a Large Area Few-Layer MoS₂ with Magnetron Sputtering.

Wei Zhong1, Sunbin Deng2, Kai Wang3, Guijun Li4, Guoyuan Li5, Rongsheng Chen6,7, Hoi-Sing Kwok7.   

Abstract

In this article, we report continuous and large-area molybdenum disulfide (MoS₂) growth on a SiO₂/Si substrate by radio frequency magnetron sputtering (RFMS) combined with sulfurization. The MoS₂ film was synthesized using a two-step method. In the first step, a thin MoS₂ film was deposited by radio frequency (RF) magnetron sputtering at 400 °C with different sputtering powers. Following, the as-sputtered MoS₂ film was further subjected to the sulfurization process at 600 °C for 60 min. Sputtering combined with sulfurization is a viable route for large-area few-layer MoS₂ by controlling the radio-frequency magnetron sputtering power. A relatively simple growth strategy is demonstrated here that simultaneously enhances thin film quality physically and chemically. Few-layers of MoS₂ are established using Raman spectroscopy, X-ray diffractometer, high-resolution field emission transmission electron microscope, and X-ray photoelectron spectroscopy measurements. Spectroscopic and microscopic results reveal that these MoS₂ layers are of low disorder and well crystallized. Moreover, high quality few-layered MoS₂ on a large-area can be achieved by controlling the radio-frequency magnetron sputtering power.

Entities:  

Keywords:  disorder; few-layer MoS2; magnetron sputtering; magnetron sputtering power; raman spectroscopy

Year:  2018        PMID: 30081483      PMCID: PMC6116247          DOI: 10.3390/nano8080590

Source DB:  PubMed          Journal:  Nanomaterials (Basel)        ISSN: 2079-4991            Impact factor:   5.076


1. Introduction

The emergence of monolayer graphene [1,2] and transition metal dichalcogenides (TMDs) [3,4] has inspired a series of high-profile discoveries in the electronic and optoelectronic fields [5,6,7], and has initiated potentially new areas [8,9]. Thus, two-dimensional (2D) materials have recently been intensively studied. In the context of 2D TMDs, molybdenum disulfide (MoS2) is one of the attractive embodiments due to its stable form in few- and single-layers [10,11] as well as its desirable electrical and optical properties [12]. Few- and single-layer MoS2 is firstly obtained via top-down mechanical stripping techniques [13], which are commonly used for graphene exfoliation. Although this method of exfoliation and dip coating [14,15,16,17] has its own the advantages to achieve high-quality 2D materials, none of these are proper solutions for radical large-scale commercial manufacturing, where the mass-producible growth of large-area, continuous, and high-quality 2D MoS2 thin films on dielectrics is a pre-requisite. From this point of view, the bottom-up strategies for thin film growth, including chemical vapor deposition (CVD) and physical vapor deposition (PVD), are better choices when compared with the top-down methods mentioned above. CVD methods have already been successfully demonstrated [18,19,20], but the control of thin film thickness, purity and uniformity on a large scale needs to be further enhanced [21]. On the other hand, PVD methods, especially magnetron sputtering, have been broadly employed in large-scale commercial manufacturing at low cost and with easy control. However, the exploration of 2D MoS2 thin film growth using magnetron sputtering technology is quite insufficient [21]. In recent years, there have been few attempts at sputtering techniques for the growth of MOS2 thin films. Muratore et al. [22], Kaindl et al. [23], and Samassekou et al. [24] reported the synthesis of continuous few-layer MoS2 by sputtering method using a MoS2 target, and the sputtered MoS2 films was annealed in an argon atmosphere. Tao et al. [21] and Santoni et al. [25] reported MoS2 film using Mo target sputtered in vaporized sulfur ambient. However, the reported films either are relatively thick or have poor crystal quality and optical properties [21,22,23,24,25]. The main obstacles for large-area sputtered high quality MoS2 films are possibly located in the difficulty of disorder control during thin film deposition, and the lack of metrics to evaluate the deposited thin films. Studies over the past decade have shown that Raman spectroscopy has historically played an important role in the structural characterization of graphitic materials [26,27,28] and has also become a powerful tool to understand the effects of process parameters on the production of high quality graphene by monitoring changes in disordered peaks [28,29,30,31,32]. Recently, Raman spectroscopy has also been used to study the effects of disorder on the MoS2 [33]. On the other hand, since many scholars have previously investigated the crystallinity of MoS2 thin film by high-temperature vulcanization or deposition on different substrates [21,34,35], few studies have investigated the effect of sputtering power on the MoS2 thin film. Therefore, in this work, we use a Raman spectroscopy approach to describe the quality and to study the effect of RF power on the deposition of the large-scale few-layer MoS2 films. We also investigated the crystalline structure of the films through an X-ray diffractometer (XRD) and a high-resolution field emission transmission electron microscope (HRTEM). The binding energies of Mo in the MoS2 film grown by radio frequency magnetron sputtering (RFMS) were further analyzed by X-ray photoelectron spectroscopy (XPS). For precise detection, the surface 5-nm-thick thin films were etched using Ar ions before XPS characterization.

2. Materials and Methods

The large-scale few layer MoS2 thin films were deposited using the RFMS technique. Firstly, both silicon substrates coated with thermally grown SiO2 and glass substrates (Eagle 2000, Corning) were ultrasonically cleaned in acetone and then isopropanol (IPA). After being rinsed in deionized (DI) water and dried, the substrates were loaded into the chamber of an RFMS system (AJA International Inc., Scituate, MA, USA) and heated to 400 °C. The distance between the substrate and the MoS2 target (99.99%, 2 Inc., Plasmaterials Inc., Livermore, CA, USA) was 10 cm. In order to suppress the influence of oxygen and moisture on the deposited thin films, the base pressure of the chamber should be pumped down as low as possible. In this work, the value was 2.67 × 10−5 Pa. During deposition, Ar gas (99.999%) was allowed to flow into the chamber with a stable flow rate of 20 sccm, and the working pressure of the chamber was maintained at 3 × 10−3 Torr. The RF sputtering power applied to the MoS2 target was varied from 10 W to 150 W in order to investigate the relationship between the large-scale few layer MoS2 growth and RF power. The deposition time was dependent on the required thickness of the thin films and the thickness of the film is maintained at 15 nm, the thickness of the thin films is confirmed by HRTEM (High-Resolution Transmission Electron Microscopy). When the RFMS process was completed, the heater and the Ar gas pipeline were switched off, and the substrates were cooled down to room temperature naturally. After loading out of the chamber, the SiO2/Si substrates and glass substrates with as-deposited MoS2 thin film were immediately subjected to post-annealing in a sulfurization atmosphere for 60 min at 600 °C to enhance thin film quality physically and chemically (Figure 1).
Figure 1

This post-deposition annealing treatment was performed to further enhance crystalline quality in as-sputtered MoS2 on glass substrates under sulfur environment.

In order to identify the layered structure and obtain the disorder information, the deposited MoS2 thin films were analyzed using Raman spectroscopy (Renishaw invia RE04, 514 nm Ar laser with a 1 μm spot size, Renishaw plc, Gloucestershire, United Kingdom). For the crystalline phase characterizations of the thin films, an X-ray diffractometer (XRD, Empyrean, PANalytical, Almelo, The Netherlands) with Cu Kα radiation was used. It was operated in thin-film mode, and the angle between the X-ray and the thin film surface was fixed at 0.5 degrees. Furthermore, a high-resolution transmission electron microscopy (HRTEM, JEM-2010HR, JEOL, Tokyo, Japan) was applied to identify the number of layers and atomic structure of the MoS2 film. To analyze the material composition of the deposited thin films as well as the chemical environment of the atoms in the thin films, X-ray photoelectron spectroscopy (XPS) measurement was conducted on the Physical Electronics 5600 multi-technique system (Physical Electronics Inc., Chanhassen, MN, USA).

3. Results and Discussion

Figure 1 shows an image of MoS2 thin layer grown on glass substrates with different deposition time (15 s, 30 s, and 45 s). The deposited MoS2 layer was light gray, and after annealing under a sulfur atmosphere, the MoS2 layer was pale yellow and was found to have specular reflection of ambient light. This result is consistent with previous reports [21]. The size of the synthesized films is limited by the dimensions of our sample heating holder. Figure 2 shows the Raman spectra of the MoS2 thin films on SiO2/Si substrates deposited under different RF powers of 10 W, 80 W, 120 W and 150 W, respectively. From Figure 2, it can be seen that all of the thin films exhibit two specific Raman characteristic peaks of MoS2, namely the in-plane () peak at ~381 cm-1 and the out-of-plane (A1g) peak at ~405 cm−1. Their peak positions herein are quite consistent with the results in other reports [24,36,37] , which the and A1g mode peaks appear at ~ 381.5 cm−1 and ~404.8 cm−1. Besides this, a relatively small LA(M) peak (a defect peak generated by LA phonons at the M point in the Brillouin zone) is also observed at ~227 cm−1. Theoretically, the mode corresponds to the S and Mo atoms oscillating in antiphase parallel to the crystal plane, while the A1g mode corresponds to the S atoms oscillating in antiphase out-of-plane, as shown in the insets of Figure 2 [29,38]. In addition to the absolute positions of these two peaks, the frequency difference (Δk) between the A1g peak and the peak is a good indicator of the layer number in MoS2 thin films. In this work, Δk is around 24 cm-1, which is larger than that in monolayer MoS2 thin films (~18–19 cm−1) [19,33,39], but smaller than the typical value in bulk MoS2 (~26 cm−1) [37,38]. This indicates the existence of few-layer MoS2 [36]. Table 1 lists the peak positions corresponding to the and A1g mode as well as the Δk values for all samples under various RF sputtering powers in this work. With the continuous increase of RF power from 10 W to 150 W, the Δk values remain in the vicinity of 24 cm−1. This means that all of the thin films under different RF powers are MoS2 with a few layers [36,37]. Another indicator of film quality is the full width at half maxima (FWHM) of the observed vibration modes. FMHM values for the sputtered FL-MoS2 film with different RF powers is compared in Table 1. In general, higher FMHM values mean more disorder [24,25]. From the Table 1, it can be seen that the MoS2 films with RF powers of 120 W has the lowest FMHM values, meaning it has the least disorder. It also can be seen that the A1g peak and the peak all have a blue shift as the RF powers increased from 10 W to 120 W and a red shift as the RF powers increased from 120 W to 150 W. The red shift is attributed to the high RF power, resulting in an increase in the residual-stress of the film [40].While the blue shift is attributed to O2-doping of MoS2, which will be shown on the XPS result [41].
Figure 2

Raman spectra of the MoS2 thin films deposited on SiO2/Si substrates under different radio frequency (RF) powers. The insets illustrate the oscillating mode of the and A1g peak.

Table 1

The - and A1g-related Raman peak information of MoS2 thin films deposited using radio frequency magnetron sputtering (RFMS) under various RF powers.

RF Power (W)A1g (cm−1)E2g1 (cm−1)Δk (A1g-E2g1) (cm−1)Full Width at Half-Maximum (cm−1)LA(M) to A1g Peak Intensity Ration
A1g E2g1
10405.2 ± 0.1381.2 ± 0.124.0 ± 0.0210.64 ± 0.4610.87 ± 0.560.219 ± 0.010
80405.4 ± 0.4381.4 ± 0.424.0 ± 0.0110.17 ± 0.3410.31 ± 0.210.195 ± 0.004
120407.0 ± 0.1383.0 ± 0.124.0 ± 0.029.55 ± 0.029.57 ± 0.200.180 ± 0.003
150403.7 ± 0.4379.7 ± 0.424.0 ± 0.0210.49 ± 0.1310.56 ± 0.020.204 ± 0.006
According to Raman fundamental selection rules, only phonons with wave vector q ≌ 0 are Raman active around the center of the Brillouin zone. However, this rule will be broken by defects, which cause the appearance of peaks away from the zone center [33,42]. Monitoring the evolution of disorder-related sub-peaks in the Raman spectrum enables us to understand the effects of process parameters and has allowed great strides to be made in the CVD preparation of high quality graphene [29,43,44]. Similarly, in Figure 2, apart from two major peaks with regard to the basic and A1g vibration modes, several sub-peaks related to the defects can also be observed. Among them, the sub-peak at 227 cm−1 is the most intense, which is attributed to the longitudinal phonon branch at the point M (LA (M)) of the Brillouin region [33,45]. Therefore, this sub-peak is able to form a very clear marker for disorder in the system, especially for the few- and single-layer MoS2. The intensity ratio of the LA(M) sub-peak to A1g peak is plotted as a function of RF power in Figure 3. It can be clearly observed that the intensity ratio reaches a minimum when the RF power climbs to 120 W, indicating the improvement of thin film quality [46]. In general, higher RF power could assist with the formation of crystalline films with lower disorder (namely, lower LA(M)/ A1g intensity ratio), but excessively high RF power is not welcomed. For instance, the LA(M)/A1g intensity ratio rises when the RF power increases from 120 W to 150 W. It reveals the increase of disorder in the deposited MoS2 thin films. A possible explanation for this phenomenon could be the increased defect generation caused by ion bombardment under over-high RF power.
Figure 3

LA(M) to A1g peak intensity ratio of the MoS2 films with different RF powers deposited on SiO2/Si substrates.

The crystalline information of the samples was characterized using XRD. The X-ray diffraction patterns of MoS2 thin films on SiO2/Si substrates under different RF powers are shown in Figure 4. The MoS2 thin films deposited under different RF powers all exhibit three obvious diffraction peaks, which are located at 14.0°, 21.6°, and 51.5° respectively. The two weak diffraction peaks at 21.6° and 51.5° correspond to the SiO2 (JCPDS: 27-0605) (111) and (400) planes, respectively. The strong diffraction peak at 14.0° is an indicator of the MoS2 (JCPDS: 37-1492) (002) plane. For the MoS2 thin films deposited under an RF power of 10 W, the broad and weak diffraction peak at 14.0° indicates the amorphous structure of the film. However, when the RF power increases to 80 W, 120 W and then 150 W, all of the MoS2 thin films exhibit a strong and narrow diffraction peak at 14.0°, which indicate the well crystallization of MoS2 thin films. Since the (002) plane is parallel to the surface of the substrates, the deposited MoS2 thin films using the RFMS technique under an RF power of 80 W, 120 W and 150 W could grow along the c-axis. Meanwhile, the exclusive diffraction peak also suggests the thin films are highly oriented. These properties are quite helpful for the formation of stacked microstructures in MoS2 thin films. Besides this, it should be noted that the intensity of the (002) peak for MoS2 thin films under an RF power of 150 W is lower compared to the other two samples. This is possibly related to the degradation of crystalline quality under such high RF powers [47]. At the same time, we also noticed that under the RF power of 150 W, the sample showed a weak diffraction peak at 28.4°, which corresponds to the Si (JCPDS: 27-1402) (100) plane.
Figure 4

X-ray diffraction patterns of the MoS2 thin films on SiO2/Si substrates under different RF sputtering powers.

To further elucidate the crystalline structure, the MoS2 thin film deposited by 120 W RF power was transferred onto a lacey copper grid for HRTEM characterizations. Figure 5a presents its cross-sectional HRTEM image, it can be seen that there are about 20 layers of the MoS2 thin films on the SiO2/Si substrate, and the interlayer spacing (0.68 nm) of the MoS2 thin films is consistent with the previous results [21]. Moreover, it is also verified that the deposited MoS2 thin films are stacked a few layers in parallel, which is in agreement with the results extracted from the Raman spectra above. As far as the typical high resolution TEM image in Figure 5b is concerned, the first-order diffraction spots of the FFT image (inset of Figure 5b) on a selected area are shown according to a regular hexagonal symmetry, thus indicating the presence of the MoS2 thin film made of single crystal domains (without Moireé patterns). Moreover, HRTEM images of the selected area in Figure 5b after FFT filtering are shown in Figure 5c. The IFFT-filtered image (Figure 5c) shows a regular honeycomb pattern due to the atomic arrangements of the Mo and S atoms. In Figure 5d, the calculated profile along the selected direction in Figure 5c is drawn, which reveals the (004) plane of MoS2 with a lattice spacing of 0.31 nm that is observed in Figure 5c. Figure 5e is the magnified image of the square-surrounded region in Figure 5c. The periodic atom arrangement for Mo confirms that the MoS2 thin films deposited under an RF power of 120 W own a crystalline structure. This is also consistent with the XRD results in Figure 4.
Figure 5

(a) Cross-sectional high-resolution field emission transmission electron microscope (HRTEM) image of samples deposited under an RF power of 120 W. (b) High resolution TEM image of the MoS2 film deposited by 120 W RF power transferred onto a lacey carbon grid. The inset shows the fast Fourier transformation (FFT) image corresponding to the TEM image selected area of a portion of (b), showing the hexagonal symmetry of the MoS2 structure. (c) Inverse FFT images corresponding to the TEM image selected area of a portion of (b). (d) Atomic spacing along the selected direction of the basal plane. (e) Zoom-in image of the area highlighted in (c). The hexagonal structure formed by Mo atoms is indicated.

XPS was used to analyze the chemical environment of Mo in the MoS2 thin films, and to detect any impurity (particularly oxygen) involvement during preparation. The high-resolution XPS spectra of the Mo 3d and S 2s region are shown in Figure 6. Since the peaks of Mo 3d and S 2s are too close to distinguish, in order to obtain the chemical environment information of the Mo species, both the S 2s and Mo 3d spectra are taken into account. Moreover, since the main Mo doublet peak signals of the thin films under an RF power of 10 W are so weak that they reach noise level, the analysis in terms of such samples is not reliable. Hence, only the data with an RF power of 80 W, 120 W and 150 W have been fitted by a 20% Lorentzian–Gaussian ration fit and the according results together with the Shirley background are presented in Figure 6b–d. According to [35], the two peaks at 229.1 eV and 232.2 eV for the MoS2 thin film are attributed to the doublet Mo 3d5/2 and Mo 3d3/2 orbitals, respectively. Meanwhile, the fitting shows that there is a second characteristic at the lower binding energy. For species originating from lower binding energy (BE) Mo3d peaks, the lower BE Mo species is Mo still associated with the Mo-S lattice, or it reflects a single amorphous MoSx phase, where Mo has a different number of nearest neighbor S atoms [48,49]. According to [25], the lower BE Mo species can be assigned to zero-valent Mo (Mo(0)) occurring in small aggregates dispersed in a MoS2 matrix. From the principle of sputter coating, we can know that the presence of Mo(0) cannot be avoided, and the purpose of annealing under a sulfur atmosphere is to convert Mo(0) to Mo(IV), which can reduce disorder and defect formation. However, when the film is etched by an ion beam, there is also a chemical shift of its binding energy toward smaller values [25,35]. Therefore, compared to XPS characterization, the advantages of strong non-destructive characterization of Raman spectra are even more pronounced. For accurate detection, 5-nano-thin-thick-surface films should be etched with Ar ions prior to XPS characterization. However, this will lead to sample destruction and will have a certain impact on the test results. Although this changed the chemical state of the Mo atom, all the samples were processed in the same way, and the overall trend of its variation with sputtering power did not change. The relative ratio of Mo species calculated from the deconvolved Mo 3d spectrum results is shown in Figure 7. The Mo(IV) fraction is defined as the area under the Mo(IV) peak divided by the total Mo peak area. It can clearly be seen that the Mo (IV) fraction reaches a maximum when the RF power is increased to 120 W, indicating the improvement of thin film quality. In addition, smaller peaks at 230.2 eV and 233.3 eV are assigned to Mo5+ [50], indicating the presence of chemisorbed oxygen at sulfur vacancies [51,52] or sub-stoichiometric oxides MoOx [53]. The origin of MoOx is probably due to the sulfurization process. We believe that MoOx has probably formed during the deposition and/or during the sulfurization process at 600 °C by reaction with the oxygen diffused from the SiO2/Si substrate. As shown in Figure 7, the Mo (V) fraction is around 0.1, 0.094 and 0.102 for an RF power of 80 W, 120 W, and 150 W, respectively. For different RF powers, the Mo(V) fraction is basically the same, indicating that the power change has no effect on the Mo(V) fraction.
Figure 6

High-resolution X-ray photoelectron spectroscopy (XPS) spectra of (a) Mo 3d/S 2s, and (b), (c), and (d) Mo3d core-level spectra for an RF power of 80 W, 120 W, and 150 W, respectively. The background is shown with a black line at the bottom. The black dots represent the raw data. The red line is the total least-squares fit. The orange lines indicate the Mo(V) 3d components. The blue lines are the Mo3d components linked to Mo (0) and Mo (IV), the Mo (0) are the lower BE doublet. The olive lines are the S2s components. (See text for more explanation).

Figure 7

Relative ratio of Mo species with various chemical states at different RF powers.

4. Conclusions

Few-layer MoS2 thin film deposition on large-area thermally oxidized silicon substrates was demonstrated using the RFMS technique. Raman analysis verified the achievement of few-layer MoS2 thin films. Meanwhile, it was proposed that the disorder inside the thin films could be monitored using Raman spectra, and controlled by adjusting the RF sputtering power. Furthermore, the XRD spectra and cross-sectional TEM images confirmed the high quality of few-layer MoS2 thin films, implying that RFMS was suitable for layered MoS2 growth. Additionally, the XPS characterizations on RFMS grown few-layer MoS2 thin films revealed the RF power has a great effect on the binding energies of Mo atoms. Our work illustrates that sputtering combined with sulfurization is a viable route for the high quality of large-area few-layer MoS2 by controlling the radio-frequency magnetron sputtering power.
  30 in total

1.  Control of valley polarization in monolayer MoS2 by optical helicity.

Authors:  Kin Fai Mak; Keliang He; Jie Shan; Tony F Heinz
Journal:  Nat Nanotechnol       Date:  2012-06-17       Impact factor: 39.213

2.  Valley polarization in MoS2 monolayers by optical pumping.

Authors:  Hualing Zeng; Junfeng Dai; Wang Yao; Di Xiao; Xiaodong Cui
Journal:  Nat Nanotechnol       Date:  2012-06-17       Impact factor: 39.213

3.  Spatially resolved Raman spectroscopy of single- and few-layer graphene.

Authors:  D Graf; F Molitor; K Ensslin; C Stampfer; A Jungen; C Hierold; L Wirtz
Journal:  Nano Lett       Date:  2007-02       Impact factor: 11.189

4.  High performance multilayer MoS2 transistors with scandium contacts.

Authors:  Saptarshi Das; Hong-Yan Chen; Ashish Verma Penumatcha; Joerg Appenzeller
Journal:  Nano Lett       Date:  2012-12-19       Impact factor: 11.189

5.  High-mobility three-atom-thick semiconducting films with wafer-scale homogeneity.

Authors:  Kibum Kang; Saien Xie; Lujie Huang; Yimo Han; Pinshane Y Huang; Kin Fai Mak; Cheol-Joo Kim; David Muller; Jiwoong Park
Journal:  Nature       Date:  2015-04-30       Impact factor: 49.962

6.  Review of chemical vapor deposition of graphene and related applications.

Authors:  Yi Zhang; Luyao Zhang; Chongwu Zhou
Journal:  Acc Chem Res       Date:  2013-10-15       Impact factor: 22.384

7.  An atlas of two-dimensional materials.

Authors:  Pere Miró; Martha Audiffred; Thomas Heine
Journal:  Chem Soc Rev       Date:  2014-09-21       Impact factor: 54.564

8.  Few-layer MoS2: a promising layered semiconductor.

Authors:  Rudren Ganatra; Qing Zhang
Journal:  ACS Nano       Date:  2014-04-14       Impact factor: 15.881

Review 9.  Studying disorder in graphite-based systems by Raman spectroscopy.

Authors:  M A Pimenta; G Dresselhaus; M S Dresselhaus; L G Cançado; A Jorio; R Saito
Journal:  Phys Chem Chem Phys       Date:  2007-01-11       Impact factor: 3.676

10.  Growth, structure and stability of sputter-deposited MoS2 thin films.

Authors:  Reinhard Kaindl; Bernhard C Bayer; Roland Resel; Thomas Müller; Viera Skakalova; Gerlinde Habler; Rainer Abart; Alexey S Cherevan; Dominik Eder; Maxime Blatter; Fabian Fischer; Jannik C Meyer; Dmitry K Polyushkin; Wolfgang Waldhauser
Journal:  Beilstein J Nanotechnol       Date:  2017-05-22       Impact factor: 3.649

View more
  2 in total

1.  A study on the tribological property of MoS2/Ti-MoS2/Si multilayer nanocomposite coating deposited by magnetron sputtering.

Authors:  Ning Kong; Boyu Wei; Dongshan Li; Yuan Zhuang; Guopeng Sun; Bo Wang
Journal:  RSC Adv       Date:  2020-03-05       Impact factor: 4.036

2.  Preparation and Properties of 2D Materials.

Authors:  Byungjin Cho; Yonghun Kim
Journal:  Nanomaterials (Basel)       Date:  2020-04-16       Impact factor: 5.076

  2 in total

北京卡尤迪生物科技股份有限公司 © 2022-2023.