Literature DB >> 35071914

Temperature Dependence of Ultrathin Mixed-Phase Ga2O3 Films Grown on the α-Al2O3 Substrate via Mist-CVD.

Abhay Kumar Mondal1, Loh Kean Ping1, Muhammad Aniq Shazni Mohammad Haniff1, Mohd Arif Mohd Sarjidan2, Boon Tong Goh2, Mohd Ambri Mohamed1.   

Abstract

Alpha (α)- and beta (β)-phase gallium oxide (Ga2O3), emerging as ultrawide-band gap semiconductors, have been paid a great deal of attention in optoelectronics and high-performance power semiconductor devices owing to their ultrawide band gap ranging from 4.4 to 5.3 eV. The hot-wall mist chemical vapor deposition (mist-CVD) method has been shown to be effective for the growth of pure α- and β-phase Ga2O3 thin films on the α-Al2O3 substrate. However, challenges to preserve their intrinsic properties at a critical growth temperature for robust applications still remain a concern. Here, we report a convenient route to grow a mixed α- and β-phase Ga2O3 ultrathin film on the α-Al2O3 substrate via mist-CVD using a mixture of the gallium precursor and oxygen gas at growth temperatures, ranging from 470 to 700 °C. The influence of growth temperature on the film characteristics was systematically investigated. The results revealed that the as-grown Ga2O3 film possesses a mixed α- and β-phase with an average value of dislocation density of 1010 cm-2 for all growth temperatures, indicating a high lattice mismatch between the film and the substrate. At 600 °C, the ultrathin and smooth Ga2O3 film exhibited a good surface roughness of 1.84 nm and an excellent optical band gap of 5.2 eV. The results here suggest that the mixed α- and β-phase Ga2O3 ultrathin film can have great potential in developing future high-power electronic devices.
© 2022 The Authors. Published by American Chemical Society.

Entities:  

Year:  2022        PMID: 35071914      PMCID: PMC8771688          DOI: 10.1021/acsomega.1c05859

Source DB:  PubMed          Journal:  ACS Omega        ISSN: 2470-1343


Introduction

In recent days, ultrawide band gap Ga2O3 has been recognized as a prevalent fourth-generation power device material, owing to its excellent intrinsic physical properties such as high dielectric constant, high breakdown field, and high Baliga’s figure of merit.[1] IIIA–VIA oxide family of Ga2O3 is composed of five phases: α, β, ε, δ, and λ[2,3] The most popular and highly studied phase in this polymorph is monoclinic β-Ga2O3 with an ultrawide band gap of 4.4–4.8 eV.[4,5] Another α-phase with an ultrawide band gap is 5.3 eV larger than the β-phase, and it is the most well-known power semiconductor material due to its superior band gap tuning and multifunctional alloy properties.[6−9] The α and β-Ga2O3 have emerged as promising candidates for novel power and optoelectronic devices.[2,3,10] It should be noted that Ga2O3 shows drastic changes in the electronic and optical properties at lower to higher transition growth temperatures. In general, the α-phase is stable up to a growth temperature of 550 °C, while the β-phase is thermodynamically stable, and the growth temperature lies between 650 and 700 °C.[11] Several methods have been reported to grow α- and β-phase Ga2O3 epilayers on the sapphire substrate, such as pulsed-laser deposition,[12−14] molecular beam epitaxy,[13,15−17] atomic layer deposition,[18] metal–organic chemical vapor deposition (CVD),[13,19] metal–organic vapor phase epitaxy,[20−23] low-pressure CVD,[24] mist-CVD,[25] halide vapor phase epitaxy,[26] and so on. α-Ga2O3 epilayer growth expedition still lags far behind β-Ga2O3 due to its instability in achieving highly crystalline epilayer for device fabrication.[27] However, mist-CVD is a low-cost, highly scalable, non-vacuum, solution-processed method to grow α-Ga2O3 thin films on the sapphire α-Al2O3 substrate.[25,28−31] Shinohara and Fujita first reported highly crystalline mist-CVD grown α-Ga2O3 heteroepitaxy on the α-Al2O3 substrate (3.4% lattice mismatch) at 400–500 °C.[25] Another study by Son et al. showed that the epilayer was stable till 600 °C using the Halide vapor epitaxy method.[32] Recently, Lu et al. demonstrated the high responsivity mixed-phase Ga2O3-based solar-blind photodetector by the pulse laser deposition method.[14] Until now, no studies have reported the ultrathin film mixed oxide phase via the mist CVD method at higher growth temperature, which could be more promising than a pure phase for high-power devices.[14] Although in the perspective of the ultrathin film, breakdown strength Ec of the thicker layer is lower than that of the thin layer due to significant lattice dislocation.[2] This mixed-phase polycrystalline Ga2O3 is the advantage of a longer recovery time, and its ultrathin nature has suitable capacitance and breakdown strength behavior. The optimum growth temperature can mostly control the material crystal quality as it provides sufficient atomic kinetics to maintain phase stability during the growth process. In this work, an ultrathin mixed-phase Ga2O3 film was successfully grown on the α-Al2O3 substrate via mist-CVD. By varying the growth temperature and keeping the other growth parameters constant, the film characteristics such as crystallinity, dislocation density, surface roughness, elemental composition, and optical properties were carefully examined using various characterization techniques, including X-ray diffraction (XRD), atomic force microscopy (AFM), and UV–visible transmission spectroscopy, Raman, X-ray photoelectron spectroscopy (XPS), and high-resolution transmission electron microscopy (HR-TEM).

Results and Discussion

Figure a shows the number of XRD crystalline peaks of the Ga2O3 epilayer grown at 470–700 °C. In these XRD crystallinity patterns, black, red, blue, and green colors indicate the crystalline peaks of different growth temperature thin films. All the samples are seen to have a mixed α- and β-phase, except at a higher growth temperature of 700 °C. At this elevated temperature, a crystalline diffraction peak appeared, which is 2θ = 37.50°, belonging to the monoclinic beta phase (401) plane, which is consistent with other temperatures. It is noticeable that this β-phase Ga2O3 exists at all growth temperatures. At 700 °C growth temperature, the most significant crystalline peak is 38.36°, corresponding to the (−402) plane monoclinic β-phase. However, here we notice that for all other growth temperatures, this diffraction peak is located at 39.41°, which belongs to the beta monoclinic phase (−402) plane lattice parameter of a = 12.23 Å, b = 3.04 Å, c = 5.80 Å, α = 90°, β = 103.71, and γ = 90° (JCPDS 43-1012).[33] This peak shift could be attributed to varying crystalline sizes and defects of the epi layer. Moreover, additional crystalline diffraction peaks located at 34.56 and 59° could belong to β-phase (11–1) and (−603) planes, respectively.[14,30] Another minor diffraction peak located at 2θ = 40.23° corresponds to the α-Ga2O3 (0006) plane with the hexagonal phase lattice constant a = b = 4.98 Å, and c = 13.43 Å, validated with COD no. 96-101-0584.[34] This peak intensity suggests that minor amounts of the α phase are present in these samples. All the samples have the diffraction peak at 41.98° belonging to the corundum structure of the α-Al2O3 (0006) sapphire substrate, followed by COD no. 96-900-9675. Figure b,c shows the unit cell view of β and α of (−402) and (0006) orientation and their atomic arrangement. However, diffraction peak intensity may vary slightly at different growth temperatures. α- and β-phase peak intensity exhibits a comparatively higher value with increasing growth temperature, indicating that more crystalline structures appear at the higher temperature. It is worth noticing that the β-phase exists at a much lower temperature, ranging from 470 to 600 °C, while a very unstable low-crystalline α-phase exists at 600 °C. Here, a mixed-phase can be seen up to 600 °C, and a pure β-phase is only noticeable at 700 °C. A previous study stated that α-Ga2O3 is usually unstable at higher temperatures. This work showed the noticeably existence of α- and β-phase thin film with a large value of dislocation density >1010 cm–2. Table shows the calculated dislocation density of Ga2O3 grown at different temperatures. Table shows the dislocation density and the surface roughness value of the Ga2O3 epilayer grown at different deposition parameters. Dislocation density was calculated from the given below Williamson Smallman relation in eq .[35]Here, D is the crystalline size that can be calculated from Debye–Scherrer’s eq , as given below[36]where K is the crystalline shape factor of 0.9, λ refers to the Cu XRD source wavelength, and β is the full width at half-maximum (fwhm) value of the corresponding crystalline peak. The dislocation density was then calculated using eq . From the calculation, the sample grown at 600 °C has a high dislocation density of 2.63 × 1010 cm–2 due to the formation of high strain over the surface area. It is noted that stress on the particular surface area is much higher than the average stress in the lattice leading to dislocations. Son et al. reported the epitaxial growth of Ga2O3 using the Halide vapor epitaxy method, and the threading dislocation density was estimated to be 1.7 × 1010 cm–2.[32]
Figure 1

(a) XRD spectra of the Ga2O3 film grown on an α-Al2O3 substrate at 470–700 °C. Unit cell of the (b) β-Ga2O3 (−402) orientation and (c) α-Ga2O3 (0006) orientation.

Table 1

Calculation of Dislocation Density of the Ga2O3 Thin Films at Different Growth Temperatures

deposition temp (°C)peak position, 2θ (deg)fwhm, β (deg)cystalline size,D (nm) D = Kλ/β cos θdislocation density 10–4 × (nm–2) δ = 1/D2average dislocation density (1010cm–2)
47037.518220.097789.691508211.243075042.82
 39.858530.1296268.090696042.156872398 
 40.723520.1995344.502013415.049410508 
50037.519520.0983189.13532821.2586363322.70
 39.859630.1302767.7511842.178543427 
 40.727070.1917446.310585044.662721381 
60037.516880.0957291.546443281.193210452.63
 39.855140.1269169.543939472.067670934 
 40.723050.1907346.402646674.644238301 
Table 2

Deposition Conditions of the Ga2O3 Thin Films

phaseprecursor concentration (M)Ar and O2 flow rate (sccm)deposition temp (°C) and time (h)average dislocation density (1010× cm–2)rms (nm)
α and β0.01350 & 50470 & 12.8216.49
α and β0.01350 & 50550 & 12.7014.79
α and β0.01350 & 50600 & 12.632.40
(a) XRD spectra of the Ga2O3 film grown on an α-Al2O3 substrate at 470–700 °C. Unit cell of the (b) β-Ga2O3 (−402) orientation and (c) α-Ga2O3 (0006) orientation. As shown in Figure a–c, Ga2O3 thin-film surface roughness and morphology were characterized by AFM with the contact operation mode. 470, 550, and 600 °C grown thin film root-mean-square (rms) surface roughness values of 16.49, 14.79, and 2.40 nm, respectively, were analyzed within an area of 5 × 5 μm2. The rms value represents the standard deviation of the surface height profile over the area.[37] It can be seen that the thin-film dislocation density reflects on the surface roughness. This dislocation density has a relation with the crystalline/grain size, which can be found out in XRD analysis. Many grains are present in the different orientations of crystallographic axes. It is shown that 600 °C grown thin films exhibit relatively lower (rms) average surface roughness in comparison of lower growth temperatures, and the surface roughness is uniform in the entire area (Figure c). As a result, temperature could play a significant role in the stabilized crystal phase. So far, gaining the higher thermal energy of thin films, β-phase crystallinity reconstruction enables ad-atoms mobility and interface strain slightly less between two crystal materials. The surface nature appears to be consisting of several nuclei growing independently and randomly distributed three-dimensional crystal particles with uniform height. The morphology of the film became more uniform by reducing its surface roughness and increasing its density. The larger lattice mismatch and residual strain caused misfit dislocation to appear on the thin-film surface. It was shown that the sample grown at 600 °C has better surface morphology than the other samples. Therefore, further characterization and analysis were conducted for the sample grown at 600 °C only.
Figure 2

AFM surface topography characteristics of the Ga2O3 film grown at (a) 470, (b) 550, and (c) 600 °C.

AFM surface topography characteristics of the Ga2O3 film grown at (a) 470, (b) 550, and (c) 600 °C. The optical transmittance spectra (250–800 nm) of Ga2O3 samples at 600 °C deposition temperatures were performed using a double-beam spectrometer with a reference α-Al2O3 substrate (Figure ). Due to the light diffraction, wave curves occur in the UV–visible range. This Ga2O3 sample showed high transmission of about 95% in the UV–visible regions, and their transmissibility began to increase significantly from UV to the visible area. The α- and β-phase Ga2O3 behave as direct band gap semiconductors, and the Tauc equation can determine the Ga2O3 band gap. Figure b shows the Tauc equation where the absorption coefficient α is directly proportional to the optical band gap, Eg.where n = 2 (direct allowed transition) and n = 1/2 (indirect allowed transition), A is the substant-dependent constant, and hυ is the photon energy.[33,38] The Eg can be determined by extrapolating the linear region of Tauc’s plot. A is the slope of the linear absorption edge associated with the structural disorder of the Ga2O3. The calculated direct band gap was 5.2 eV for the corresponding grown temperature at 600 °C. This optical band gap lies between 4.9 and 5.1 eV, as described by Kumar et al.[39] Also, these band gap characteristics show a redshift absorption edge in this deposition temperature.
Figure 3

UV–visible optical transmittance spectrum (250–800 nm) of the Ga2O3 thin films grown at 600 °C. The insets show Tauc plots of (αhν)2 as a function of photon energy (E).

UV–visible optical transmittance spectrum (250–800 nm) of the Ga2O3 thin films grown at 600 °C. The insets show Tauc plots of (αhν)2 as a function of photon energy (E). Figure shows the measured Raman spectrum of Ga2O3 thin films grown at 600 °C, taken in the range of wavenumber 300–800 cm–1. This Raman spectroscopy tool used a powerful 532 nm excitation source and yielded a 1 μm spot size. In the Raman shift, 419 and 750 cm–1 peaks indicate the sapphire peaks.[40] The Raman peak, 451 cm–1, belongs to the C32h space group β-Ga2O3, representing the β-phase bending of GaO4 tetrahedral.[41] The lower Raman peak at 431 cm–1 belongs to α-Ga2O3, which is attributed to the Raman allowed vibrational modes Eg. This mode is attributed to Ga atoms vibrating against each other along the c-axis. Another peak A1g mode at 577 cm–1 is due to the lighter oxygen atom vibrations perpendicular to the c-axis.[42,43] These obtained Raman modes agree with the experimental and theoretical values reported in the literature.[42,43] The low intensity of Raman modes of α-Ga2O3 may be due to little amount of the α phase content in this sample.
Figure 4

Raman spectrum of the mixed-phase Ga2O3 thin films grown at 600 °C.

Raman spectrum of the mixed-phase Ga2O3 thin films grown at 600 °C. XPS measurements were performed to characterize the surface elemental chemical state of the epilayer. The XPS survey scan exhibited Al, Ga, and O element of the Ga2O3 thin film grown on the α-Al2O3 substrate, as shown in Figure a. Figure b presents the Ga 2p spectrum that can be split into two main peaks such as Ga 2p3/2 and Ga 2p1/2, having binding energy appearing at 1115.79 and 1142.62 eV. The spin split between these two peaks is approximately 27 eV, consistent with the literature.[44] Meanwhile, the binding energy of the Ga 3d peak is 21 eV, which is in good agreement with the literature.[45]Figure c shows that the O 1s can be deconvoluted into two peaks using the Gaussian fitting method. The two peaks located at about 528.1 and 529.2 eV correspond to the Ga–O bond (peak I) and the oxygen-related vacancies (peak II), respectively.[46] The hybridization state of surface chemical compositions from the Ga 3d, O 1s, and C 1s peaks was estimated based on the respective sensitivity factors.[47] This is consistent with the results of recent findings on the relationship between oxidation and crystal structures.
Figure 5

(a) XPS survey scan of the as-grown Ga2O3 thin films at 600 °C. (b) Ga 2p and (c) O 1s peaks.

(a) XPS survey scan of the as-grown Ga2O3 thin films at 600 °C. (b) Ga 2p and (c) O 1s peaks. The cross-sectional TEM images revealed the thin-film thickness and interfacial structure of the thin heteroepitaxial mixed-phase Ga2O3. Figure a–c shows the cross-sectional TEM images of the as-grown Ga2O3 thin films on the α-Al2O3 substrate at 600 °C, from which 20.35 nm thick Ga2O3 thin films were obtained with a growth rate of 0.339 nm/min. The Ga2O3 microstructure formed an amorphous and crystalline region, and the Ga2O3/Al2O3 interface is not obvious, resulting in high dislocation in-between layers. These misfit dislocations are caused by the interfacial strain formed in the lattice mismatch, as observed by Ma et al.[31] The FFT measured interplanar d spacing value is 0.28 nm, corresponding to the (−202) plane of single-crystalline β-Ga2O3 (see Figure d), which is identical to the XRD pattern. Figure e shows the fast Fourier transform (FFT) pattern of Ga2O3/Al2O3. This Ga2O3 layer is not so obvious and uniform due to the presence of dislocations. The FFT pattern shows the rectangular diffraction patterns in the dark field virtual image. Nevertheless, these diffraction spots are so clear with periodic order of uniform distance. The α-Al2O3 diffraction spots are beneath Ga2O3.
Figure 6

(a) Cross-sectional TEM of the as-grown Ga2O3 thin films on α-Al2O3 at 600 °C. (b) Cross-sectional high-resolution TEM image showing 20.35 nm thick Ga2O3 films. (c) Close-up view of the Ga2O3 thin films in the yellow dashed line box in (b). (d) Filtered HR-TEM image in the green dashed line box in (c). (e) FFT pattern of Ga2O3/α-Al2O3.

(a) Cross-sectional TEM of the as-grown Ga2O3 thin films on α-Al2O3 at 600 °C. (b) Cross-sectional high-resolution TEM image showing 20.35 nm thick Ga2O3 films. (c) Close-up view of the Ga2O3 thin films in the yellow dashed line box in (b). (d) Filtered HR-TEM image in the green dashed line box in (c). (e) FFT pattern of Ga2O3/α-Al2O3. Figure shows cross-sectional analysis of elemental mapping of the Ga2O3 thin films obtained using the energy-dispersive X-ray spectroscopy (EDX) system from the orange box, as shown in Figure a. The EDX results provide the elemental composition and mapping for each Ga, O, Al, and Pt, as shown in Figure b–f. Here, the results revealed that the Ga has a low concentration at about 3.05 atomic %. It is interesting to note that EDX analysis provides the bulk concentration of elements present in the sample. The high O concentration at 36.60 atomic % may represent both Ga2O3 thin films and α-Al2O3 substrates. Hence, the values of O concentration are inaccurate and not comparable to that of Ga concentration. From the EDX analysis (distribution of O and Ga) here, we can confirm the presence of Ga2O3 thin films grown on the α-Al2O3 substrate.
Figure 7

TEM–EDX measurements of the as-grown Ga2O3 thin films on α-Al2O3 at 600 °C. (a) Cross-sectional TEM image of the Ga2O3 thin films on the α-Al2O3 substrate. (b) Integrated EDX spectrum of the sample scanned in the orange box in (a). The inset shows the elemental composition in atomic and weight percentage. (e–f) EDX mapping for Ga, O, Al, and Pt element, respectively.

TEM–EDX measurements of the as-grown Ga2O3 thin films on α-Al2O3 at 600 °C. (a) Cross-sectional TEM image of the Ga2O3 thin films on the α-Al2O3 substrate. (b) Integrated EDX spectrum of the sample scanned in the orange box in (a). The inset shows the elemental composition in atomic and weight percentage. (e–f) EDX mapping for Ga, O, Al, and Pt element, respectively.

Conclusions

The low-cost mist-CVD method has been employed to grow the mixed-phase Ga2O3 thin films on the α-Al2O3 substrate at deposition temperatures of 470–600 °C. Above 600 °C temperature, only a pure β-phase formed on Ga2O3 thin films on the α-Al2O3 substrate. It was found that the Ga2O3 thin films grown at 600 °C have high dislocation density associated with smooth surface roughness. The UV–vis characterization shows that the excellent transmission nature with an optical band gap of 5.2 eV and an ultrathin film of 20.23 nm at 600 °C may be due to the lower precursor concentration with the lower oxygen flow rate. These mixed-phase Ga2O3 thin films would help lead to new promising properties for transistors and deep UV photodetectors.

Experimental Studies

In our work, Ga2O3 thin films were grown on the α-Al2O3 substrates via the low-cost mist-CVD method. As shown in Figure a, the mist-CVD method utilizes three main components: mist formation, gas supply and reaction furnace. Table shows that the Ga2O3 thin films were grown on the α-Al2O3 substrate by varying the growth temperature with the minimal oxygen flow rate. A gallium precursor was prepared by dissolving 0.01 mol of gallium acetylacetonate (0.367 gm) into 100 mL of distilled water (DI). Then, a small amount of 1 ml HCl was dropped into the solution to have better dissolution of gallium acetylacetonate during the stirring process. The gallium precursor was atomized by a 1.7 MHz piezoelectric ultrasonication transducer to produce mist in the bubbler unit. The schematic diagram of the overall growth process is shown in Figure b. Before growing the Ga2O3 thin films, the α-Al2O3 substrates were cleaned by rinsing with ethanol, acetone, and DI water for 10 min each at room temperature. Then, α-Al2O3 substrates were placed manually in the quartz tube and heated up to the desired growth temperature (470, 550, 600, and 700 °C) in an Argon (Ar) atmosphere for 40 min (t1). After the furnace reached the desired growth temperature, the Ga-containing mist that was kept in a bubbler was transported into the reaction furnace with a 350 sccm Ar flow as carrier gas while another 50 sccm oxygen (O2) flow was transported separately. The O2 gas was used to compensate oxygen vacancies for enhancing the crystallinity of Ga2O3. The typical growth time was about 60 min (t2) at ambient pressure. After completion of the growth, the as-grown Ga2O3 samples were cooled to room temperature in an Ar atmosphere for the first 20 min and in air afterward.
Figure 8

(a) Schematic representations of the mist-CVD system to grow mixed-phase Ga2O3 thin films on the α-Al2O3 substrate. (b) Schematic diagram of the growth temperature as a function of time for the growth of Ga2O3.

(a) Schematic representations of the mist-CVD system to grow mixed-phase Ga2O3 thin films on the α-Al2O3 substrate. (b) Schematic diagram of the growth temperature as a function of time for the growth of Ga2O3. Structural crystallinity information of Ga2O3 thin films was analyzed by thin-film high-resolution XRD (SmartLab, Rigaku, Japan) with an X-ray source Cu Kα1 (λ = 1.54056 Å). The thin-film surface roughness profile measurement was performed using Park systems’ AFM. The optical transmission spectra of the samples were characterized using a dual-beam LAMBDA 950 PerkinElmer UV Vis-spectrometer. Raman spectroscopy was carried out using the Thermo Scientific Model: DXR2xi system with 532 nm solid-state laser diode excitation. XPS (Axis Ultra DLD) with monochromatic Al Kα (1486.6 eV) characterized binding energy, and surface elemental analysis was used. Thin-film cross-sectional thickness and elemental percentage were studied by transmission electron microscopy (HR-TEM) of Talos L120C.
  3 in total

1.  Influence of Carrier Gases on the Quality of Epitaxial Corundum-Structured α-Ga2O3 Films Grown by Mist Chemical Vapor Deposition Method.

Authors:  Yu Xu; Chunfu Zhang; Yaolin Cheng; Zhe Li; Ya'nan Cheng; Qian Feng; Dazheng Chen; Jincheng Zhang; Yue Hao
Journal:  Materials (Basel)       Date:  2019-11-07       Impact factor: 3.623

2.  First-Principles Studies for Electronic Structure and Optical Properties of p-Type Calcium Doped α-Ga2O3.

Authors:  Abhay Kumar Mondal; Mohd Ambri Mohamed; Loh Kean Ping; Mohamad Fariz Mohamad Taib; Mohd Hazrie Samat; Muhammad Aniq Shazni Mohammad Haniff; Raihana Bahru
Journal:  Materials (Basel)       Date:  2021-01-28       Impact factor: 3.623

3.  First-Principles Studies for Electronic Structure and Optical Properties of Strontium Doped β-Ga2O3.

Authors:  Loh Kean Ping; Mohd Ambri Mohamed; Abhay Kumar Mondal; Mohamad Fariz Mohamad Taib; Mohd Hazrie Samat; Dilla Duryha Berhanuddin; P Susthitha Menon; Raihana Bahru
Journal:  Micromachines (Basel)       Date:  2021-03-24       Impact factor: 2.891

  3 in total
  1 in total

1.  Characteristics of Vertical Ga2O3 Schottky Junctions with the Interfacial Hexagonal Boron Nitride Film.

Authors:  Venkata Krishna Rao Rama; Ajinkya K Ranade; Pradeep Desai; Bhagyashri Todankar; Golap Kalita; Hiroo Suzuki; Masaki Tanemura; Yasuhiko Hayashi
Journal:  ACS Omega       Date:  2022-07-22
  1 in total

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