Literature DB >> 31888184

Functionalized rGO Interlayers Improve the Fill Factor and Current Density in PbS QDs-Based Solar Cells.

Anton A Babaev1, Peter S Parfenov1, Dmitry A Onishchuk1, Aliaksei Dubavik1, Sergei A Cherevkov1, Andrei V Rybin1, Mikhail A Baranov1, Alexander V Baranov1, Aleksandr P Litvin1, Anatoly V Fedorov1.   

Abstract

Graphene-quantum dot nanocomposites attract significant attention for novel optoelectronic devices, such as ultrafast photodetectors and third-generation solar cells. Combining the remarkable optical properties of quantum dots (QDs) with the exceptional electrical properties of graphene derivatives opens a vast perspective for further growth in solar cell efficiency. Here, we applied (3-mercaptopropyl) trimethoxysilane functionalized reduced graphene oxide (f-rGO) to improve the QDs-based solar cell active layer. The different strategies of f-rGO embedding are explored. When f-rGO interlayers are inserted between PbS QD layers, the solar cells demonstrate a higher current density and a better fill factor. A combined study of the morphological and electrical parameters of the solar cells shows that the improved efficiency is associated with better layer homogeneity, lower trap-state densities, higher charge carrier concentrations, and the blocking of the minor charge carriers.

Entities:  

Keywords:  impedance spectroscopy; quantum dots; reduced graphene oxide; solar cells

Year:  2019        PMID: 31888184      PMCID: PMC6947317          DOI: 10.3390/ma12244221

Source DB:  PubMed          Journal:  Materials (Basel)        ISSN: 1996-1944            Impact factor:   3.623


1. Introduction

Many research groups are attempting to improve the efficiency and stability of lead sulfide colloidal quantum dots (QDs) based solar cells [1,2,3,4]. In addition to being tunable throughout the NIR region bandgap, low-cost solution processing and the capacity for multiple exciton generation make lead sulfide QDs a highly suitable material for solar cells. Currently, 12% power conversion efficiency (PCE) is the record for heterojunction PbS QDs-based solar cells [5]. However, it is still beyond the theoretical limits of this material [6]. Defects in the active layer are among the main sources of charge carrier losses in solar cells. In solar cells based on QDs, the process of ligand replacement causes a rise in the defects in the active layer. Long ligands are replaced by short ones, which leads to the rearrangement of QDs and, as a result, to the formation of cracks and pinholes [7]. Nonetheless, the interdot spacing and layer homogeneity are crucial for layer conductivity [8]. In a number of works, graphene and its derivatives were applied to improve the parameters of solar cells’ active layers [9,10,11,12,13], including morphology. Graphene and its derivatives possess an excellent charge mobility and adjustable properties, which make them a promising material for photovoltaic applications [14]. Recently, Beatriz Martín-García et al., presented a noticeable technique [15] of reduced graphene oxide (rGO) functionalization by (3-Mercaptopropyl) trimethoxysilane (MPTS) for rGO-PbS QD linking that has vast perspectives for device fabrication. Later, we showed that a preliminary colloidal-phase ligand exchange allows a higher efficiency of charge transfer between MPTS-bonded QDs and rGO to be achieved [16]. Additionally, it was shown that a thin layer of MPTS-bonded rGO-PbS deposited on top of a hole-transporting layer improves the stability of a solar cell [17]. Here, we show that MPTS-functionalized rGO (f-rGO) can be inserted into the active layer of a PbS QDs-based solar cell. The interlayers of f-rGO drastically boost the performance of a heterojunction PbS QDs-based solar cell due to increased current density and improved fill factor (FF).

2. Materials and Methods

2.1. Materials

Precursors PbO (99.999%), hexamethyldisilathiane (TMS), zinc acetate dehydrate (99.99% Aldrich, Saint Louis, MO, USA), MAI (methylammonium iodide), EDT 1,2-ethanedithiol, TBAI (tetrabutylammonium iodide), and DMF (dimethylformamide anhydrous) were purchased from Sigma-Aldrich (Saint Louis, MO, USA); acetonitrile from AppliChem (Darmstadt, Germany); octadecene (ODE) was obtained from Acros (Geel, Belgium); oleic acid (OlAc) from Fisher Chemicals (Waltham, MA, USA); EtOH (ethanol), MeOH (methanol), and chloroform were purchased from Vekton (St. Petersburg, Russia). Reduced graphene oxide (rGO) stabilized with poly (sodium 4-styrenesulfonate) was purchased from Sigma-Aldrich as an aqueous solution and used as received. The 95% (3-Mercaptopropyl) trimethoxysilane (MPTS) was purchased from Sigma-Aldrich. PbS QDs were synthesized using the hot-injection method described elsewhere [18,19]. Briefly, 0.088 g of PbO was dissolved with 0.26 mL OlAc in 3.8 mL ODE at 90 °C under vacuum as a lead precursor. The sulfur precursor was prepared by dissolving 70 μL TMS in 2 mL ODE. The sulfur precursor was swiftly injected in the lead precursor at 85 °C, heated for 2 min, and cooled using a water-bath. The resulting QDs with a diameter of 3.5 nm were obtained by washing them with acetone and centrifugation at 6000 RPM. The PL and absorption spectra of QDs are shown in Supplementary Materials Figure S1. For the MAI phase ligand-exchange process, 1.5 mL of a 0.04 M MAI in 2:1 toluene-DMF solution was added at the rate of (3 µL·s−1) to 4 mL of the PbS QD solution with a concentration of 50 mg/mL, under gentle stirring. The resulting solution was kept under an inert atmosphere for 12 h to complete the exchange. Following this, MeOH was used to precipitate the solution, and finally, the solution was redispersed in octane [20]. The solution in octane was ultrasonicated for 15 min and filtered through a 0.22 µm syringe hydrophobic filter prior to use. The method applied by Beatriz Martin-Garcia et al. [15] was used for rGO functionalization. Briefly, 10 mL of 0.5 mg/mL rGO dispersion in EtOH was sonicated for 30 min and refluxed at 60 °C for 15 h with 1.25 mL of MPTS. The resulting mixture was washed with EtOH by centrifugation and finally dispersed in EtOH by sonication for 30 min. The solution was additionally sonicated for at least 15 min just before use. ZnO nanoparticles were synthesized according to the literature with some modifications [21]. Thirteen millimoles of zinc acetate dehydrate was dissolved in 125 mL of methanol at 60 °C. Twenty-six millimoles of potassium hydroxide was dissolved in 65 mL of methanol. The potassium hydroxide solution was slowly added to the zinc acetate solution and the solution was left stirring at 65 °C for 3 h. ZnO nanocrystals were extracted by centrifugation and then washed twice with methanol followed by centrifugation. Finally, 10 mL of chloroform and 10 mL of methanol were added to the precipitates and the solution was filtered with a 0.45 mm filter.

2.2. Device Fabrication

The devices were fabricated following the layer-by-layer spin-coating method according to [22]. Fifty nanometers-thick ZnO layers were deposited on the clean pre-patterned ITO-coated glasses at 3000 rpm and annealed at 150 °C for 30 min. After cooling to room temperature, 6 layers of MAI-treated QDs from 50 mg/mL octane solution were deposited at 2500 rpm on the ZnO film with a 2-step ligand-exchange procedure for each layer, using 15 mg/mL TBAI in methanol (MeOH) solution for exchange and acetonitrile (ACN) and hexane (HEX) for washing. Furthermore, 2 layers of OA-capped PbS QDs were deposited from a HEX solution at 2500 rpm on the TBAI layers, with ligand exchange using 0.1% EDT in ACN as an exchange solution and ACN and HEX for washing. Finally, a 100 nm Au electrode was deposited by magnetron plasma sputtering.

2.3. Measurements

The atomic force microscopy (AFM) measurements were provided using a Solver PRO-M microscope (NT-MST, Moscow, Russia) in the semi-contact mode. The scanning electron microscopy (SEM) measurements were performed using a Merlin Zeiss electron microscope (Carl Zeiss, Oberkochen, Germany) in high vacuum mode at 10 kV accelerating voltages. The light and dark J–V characteristics were measured by an Ossila Solar Cell I-V Test System. A solar simulator was used based on an OSRAM XBO-150W/1 xenon lamp (OSRAM, Munich, Germany). The light intensity adjusted was 100 mW/cm2. Capacitive and frequency characteristics were measured by a Keysight E4980A LCR Meter (Keysight Technologies, Santa Rosa, CA, USA).

3. Results and Discussion

The conventional ITO-ZnO-PbS(TBAI)-PbS(EDT)-Au solar cell architecture, with MAI passivation QDs for the TBAI layer, was employed for our study. The preliminary colloidal-phase MAI of PbS QDs surface has been recently used to enhance the performance of PbS-based solar cells [20]. Additionally, we have recently shown that such a ligand exchange provides better charge transfer in rGO-PbS QD composites [16]. Devices with the three types of the active layer were explored, as shown in Figure 1a. For the reference device (Device 1), only MAI-treated QDs were employed for the active layer fabrication. For the f-rGO interlayers device (Device 2), the same conditions for Device 1 were used, but two drops of the 0.05 mg/mL f-rGO in EtOH solution were additionally spin-coated on the interfaces between the TBAI-treated PbS layers at 3000 rpm. For the device with an f-rGO-PbS hybrid in the active layer (Device 3), 0.5 mL of 50 mg/mL MAI-treated PbS colloidal solution in octane was mixed with 0.05 mL of 0.5 mg/mL f-rGO in EtOH solution by stirring for 30 min at 2000 rpm. These hybrid inks were employed instead of QD inks for the fabrication of six TBAI layers in Device 3 with the same deposition and ligand-exchange conditions as for Device 1. The PbS(EDT) hole transport layer remained the same for all the devices. The light parameters of the devices are represented in Table 1.
Figure 1

(a) Schematic illustration of the reference (Device 1), interlayers (Device 2), and hybrid structures (Device 3); (b) J–V curves of the devices.

Table 1

The fill factor (FF), JSC, and VOC values of the reference, functionalized reduced graphene oxide (f-rGO) interlayers, and f-rGO mix devices.

DeviceFF (%)JSC (mA⋅cm−2)VOC (V)PCE (%)
Reference (1)48.4 ± 3.513.5 ± 1.20.57 ± 0.0053.5 ± 0.4
Interlayers (2)54.2 ± 0.314.2 ± 1.30.55 ± 0.0054.2 ± 0.35
Hybrid (3)42.6 ± 610.8 ± 0.20.54 ± 0.0182.5 ± 0.4
Figure 1b presents the J–V curves obtained for the three types of devices. As seen, the devices with f-rGO interlayers possess better current, while the f-rGO mixture device suffers current quenching, and both devices exhibit slight VOC reduction. Moreover, Device 2 has a significant FF improvement. The FF, short circuit current density (JSC) and open circuit voltage (VOC) values are given in Table S1. The AFM, capacitance–voltage (C–V), and capacitance–frequency (C–F) measurements were performed to reveal the nature of the changes. The AFM measurements were performed to examine the influence of the rGO on the film quality. The following samples were compared: 2-layer sample of TBAI-treated QDs; 2-layer sample of TBAI-treated QDs with rGO interlayer; 2-layer sample of TBAI-treated QDs with f-rGO interlayer; 2-layer sample of TBAI-treated layers from f-rGO-PbS hybrid inks. The samples obtained were measured on 20 × 20 µm of 10 randomly selected areas on each sample, and average roughness (Ra) data were calculated by NT-MDT NOVA image analysis software and averaged on the measured points. The examples of AFM images are shown in Figure 2, and the average roughness data are given in Table S2. Additionally, the listed structures were studied using SEM (Figure S2). The lower average roughness indicates a lower amount of cracks and pinholes. From the data, it follows that the f-rGO interlayers help minimize the quantity of the defects in the planar PbS layer-by-layer film, while using f-rGO–PbS hybrid inks results in huge aggregates in the film. In contrast, the non-functionalized rGO layer corrupts the film quality, which confirms that the MPTS linker attached to rGO plays a key role in the improvement of the films’ morphology.
Figure 2

Typical atomic force microscopy (AFM) images (a) tetrabutylammonium iodide (TBAI) treated quantum dots (QDs); (b) TBAI-treated QDs with rGO interlayer; (c) TBAI-treated QDs with f-rGO interlayer, and (d) TBAI-treated layers from f-rGO–PbS hybrid inks.

The Ra correlates with the FF and JSC of the devices and is totally independent from VOC change. VOC is sensitive to the charge carriers’ recombination [23,24], which could be enhanced by the rGO incorporation, while the better the film quality the better the QD layer resistance, thus increasing the current and FF. The dark J–V characteristics were measured for Devices 1–3, as shown in Figure S3. The resistances and ideality factors of the devices were obtained from these curves. The ideality factors are in the range of 3.5 ± 0.3, which is typical for the devices of this type [25,26,27]. The ideality factors and other device parameters are given in Table 2.
Table 2

The calculated characteristics of the devices.

SampleRelative Permittivity εND (cm−3)nVbi (V)RS Dark|Light (Ohm × cm2)RSH Dark|Light (kOhm × cm2)
Referance (1)19.85.3 × 10153.30.572.8|11.00.86|0.12
Interlayers (2)16.410.4 × 10153.80.574.3|17.50.82|0.17
Hybrid (3)18.96.7 × 10153.50.572.4|13.00.75|0.28
Next, the dark C–F characteristics were obtained, as shown in Figure S4. The contribution of deep traps to the device capacitance is minimal in flat areas [28]; therefore, we chose a frequency of 10 kHz—a relatively flat frequency—for all the devices to measure C–V. The C–V characteristics obtained are shown in Figure 3.
Figure 3

Capacitance of the devices at different bias conditions.

From the minimum capacitance in Figure 3, which is observed at a negative bias, the dielectric constant can be determined from the equation for the capacitance of a flat capacitor, [29], where Cdep is the measured capacitance, ε is the relative dielectric constant, S is the pixel area, Wdep is the depletion region width of the device that equals the geometric thickness when the device is fully depleted, and ε0 is the vacuum permittivity. The reverse bias at which the capacitance reaches the minimum was −0.5 V, −0.9 V, and −0.4 V for 1, 2, and 3 devices, respectively. The calculated dielectric constants of f-rGO devices are slightly lower than the value obtained for Device 1, whereas the addition of material with a higher dielectric constant should increase the total dielectric constant. With a large forward bias, the capacitance of Device 2 stays positive, while the capacitance of the remaining devices becomes negative. This indicates the transition of the capacitance into the inductance. The inductance effect is present in diodes, due to the modulation of the layers’ conductivity, when the density of minority charge carriers exceeds a certain part of the main density of charge carriers [30]. Thus, we believe that the difference at a high forward bias may occur due to the blocking of minor charge carriers by f-rGO layers. Figure 4 shows the Mott-Schottky plot, which is described by Equation (1): where ε is the relative dielectric constant; Vbi is the built-in potential; Vr is the applied voltage; ND is assumed to be the concentration of charge carriers. Generally, the ND value in Equation (1) is the density of uncompensated donors or acceptors [31], but in our case, charge carrier transport in PbS-TBAI is dominated by electron transport since PbS-TBAI shows n-type behavior [32,33,34]; therefore, we assume ND as the concentration of free charge carriers (electrons).
Figure 4

Mott-Schottky plot.

When C−2 tends to 0, Vbi = Vr; thus, extrapolating the curve C−2 to zero, we can obtain the values of Vbi. As one can see from Figure 4, the Vbi obtained this way for Device 2 is much larger than for the other devices. However, this is unlikely because the Vbi value is usually determined by the material of the electrodes and the interfaces on the electrodes, which remain the same for all the devices. Additionally, a real change in Vbi should induce a change in the position of the capacitance peak [34]. Apparently, as seen in Figure 4, the position of the peak is approximately the same for all the devices. For this reason, Vbi was taken as the constant and was determined directly from the Mott-Schottky plot of Device 1. The error in the determination of Vbi is associated with the rise in the curve of Device 2 on the Mott-Schottky graph, which can be explained by the appearance of additional capacitance, which is connected to the main junction in series. If the additional capacitance is in the active layer, the effect of series capacitance should maximally emerge at a forward bias and be minimal at a reverse bias, since the energy levels at which the charge accumulates are also depleted in a reverse bias and filled in a forward bias. This is confirmed by Figure 4, where the minimum difference is observed in the fully depleted mode, and the maximum difference is observed when the depletion zone has a minimum width. Additional capacitance can come from different sources, for example, an oxide layer [35] or deep traps [36]. We believe that it arises due to the introduction of f-rGO layers. Thus, the appearance of additional capacitance leads to the underestimation of the total capacitance measured for Device 2, which explains the lower calculated value of the dielectric constants. The C−2 curve of Device 2 in Figure 4 demonstrates a flatter slope at the reverse bias. Such a slope is usually associated with the existence of deep traps, which inject charge carriers even with the depletion region reaching the electrodes [35]. Thus, the sequential f-rGO deposition leads to the appearance of deep traps, which act as an additional capacitance. Another effect of the f-rGO introduction into the active layer of devices is an increase in the concentration of charge carriers, which can be determined from the angle of inclination of the curve C−2 from the graph in Figure 4. The values calculated are shown in Table 2. The increase in the concentration of charge carriers correlates with an increase in the efficiency of the devices observed from the light J–V characteristics. The density of states inside the band gap can be calculated from the С–F dependence shown in Figure S4 as follows: where C is the capacitance; ω is the angular frequency; w is the width of the depleted zone; kB is the Boltzmann constant; T is the temperature in Kelvins; E is the energy associated with the applied frequency by the following Equation (3): where ω0 = 2πv0, and v0 is a constant, which for the devices can be considered equal to ν0 = 2.88 × 1012 s−1 [37]; Eω is the energy level depending on the type of material (EV for p-type and EC for n-type). The plot of g(E) is shown in Figure 5a. The pronounced peak in the plot of g(E) in the region of 0.30–0.35 eV corresponds with the deep trap states of PbS QDs [37]. These states corrupt the efficiency of the device and their contribution increases with the degradation of QDs caused by atmospheric air [34]. As seen from Figure 5a, the density of states in the region of 0.30–0.35 eV is much lower for both devices with f-rGO. It is known that rGO [17] and other allotropic forms of carbon [38] increase the stability of photovoltaic devices with employed architecture due to the protection from atmospheric moisture. This correlates with our results because the devices were fabricated and measured in ambient conditions.
Figure 5

(a) Density of the states and profile of the devices (g(E)), and (b) energy levels and structure of the devices.

For Device 2, a significant increase in the density of trap states in the energy region above 0.45 eV is observed. This is in agreement with the results obtained from C–V measurements, which indicate the appearance of additional deep traps in Sample 2. It should be noted that the influence of the middle bandgap states on the device efficiency is lower than the influence of trap states between 0.30 and 0.35 eV, since the rate of trap state emission exponentially depends on their depth [28]. Nevertheless, such states affect the VOC [23,39] via trap-assisted recombination, which correlates with the results obtained for our devices from the light J–V characteristics and slight n growth.

4. Conclusions

We have demonstrated different approaches to embed the functionalized rGO in the active layer of PbS QDs-based solar cells. Employed as interlayers between TBAI-treated layers, functionalized rGO improves the fill factor and current density in PbS QDs-based solar cells. By using AFM, C–V, and C–F analysis, we have shown that such an increase in device performance is induced by the optimization of the active layer quality via arranging the QDs’ layer growth, having a higher concentration of charge carriers, and by blocking the minor charge carriers. Additionally, analyzing the density of states revealed the reduction of trap states between 0.3 and 0.35 eV associated with QD degradation. The method demonstrated is universal and can be applied elsewhere to fabricate QD multilayer solids.
  19 in total

1.  Colloidal Quantum Dot Solar Cells.

Authors:  Graham H Carey; Ahmed L Abdelhady; Zhijun Ning; Susanna M Thon; Osman M Bakr; Edward H Sargent
Journal:  Chem Rev       Date:  2015-06-24       Impact factor: 60.622

2.  Enhanced mobility and controlled transparency in multilayered reduced graphene oxide quantum dots: a charge transport study.

Authors:  Mawin J M Jimenez; Rafael F de Oliveira; Carlos César B Bufon; Marcelo A Pereira-da-Silva; Varlei Rodrigues; Ângelo L Gobbi; Maria H O Piazzetta; Fernando Alvarez; Carlos Lenz Cesar; Antonio Riul
Journal:  Nanotechnology       Date:  2019-03-20       Impact factor: 3.874

3.  Optoelectronic Properties of Semiconductor Quantum Dot Solids for Photovoltaic Applications.

Authors:  A A Chistyakov; M A Zvaigzne; V R Nikitenko; A R Tameev; I L Martynov; O V Prezhdo
Journal:  J Phys Chem Lett       Date:  2017-08-21       Impact factor: 6.475

4.  Hybrid organic-inorganic inks flatten the energy landscape in colloidal quantum dot solids.

Authors:  Mengxia Liu; Oleksandr Voznyy; Randy Sabatini; F Pelayo García de Arquer; Rahim Munir; Ahmed Hesham Balawi; Xinzheng Lan; Fengjia Fan; Grant Walters; Ahmad R Kirmani; Sjoerd Hoogland; Frédéric Laquai; Aram Amassian; Edward H Sargent
Journal:  Nat Mater       Date:  2016-11-14       Impact factor: 43.841

5.  In Situ Passivation for Efficient PbS Quantum Dot Solar Cells by Precursor Engineering.

Authors:  Yongjie Wang; Kunyuan Lu; Lu Han; Zeke Liu; Guozheng Shi; Honghua Fang; Si Chen; Tian Wu; Fan Yang; Mengfan Gu; Sijie Zhou; Xufeng Ling; Xun Tang; Jiawei Zheng; Maria Antonietta Loi; Wanli Ma
Journal:  Adv Mater       Date:  2018-03-15       Impact factor: 30.849

6.  Hole transport layer selection toward efficient colloidal PbS quantum dot solar cells.

Authors:  Gang Yang; Yongsheng Zhu; Jinshu Huang; Xiumei Xu; Shaobo Cui; Zhiwen Lu
Journal:  Opt Express       Date:  2019-09-30       Impact factor: 3.894

7.  Understanding the High Performance of over 15% Efficiency in Single-Junction Bulk Heterojunction Organic Solar Cells.

Authors:  Akchheta Karki; Joachim Vollbrecht; Alana L Dixon; Nora Schopp; Max Schrock; G N Manjunatha Reddy; Thuc-Quyen Nguyen
Journal:  Adv Mater       Date:  2019-10-09       Impact factor: 30.849

8.  High-Efficiency PbS Quantum-Dot Solar Cells with Greatly Simplified Fabrication Processing via "Solvent-Curing".

Authors:  Kunyuan Lu; Yongjie Wang; Zeke Liu; Lu Han; Guozheng Shi; Honghua Fang; Jun Chen; Xingchen Ye; Si Chen; Fan Yang; Artem G Shulga; Tian Wu; Mengfan Gu; Sijie Zhou; Jian Fan; Maria Antonietta Loi; Wanli Ma
Journal:  Adv Mater       Date:  2018-05-02       Impact factor: 30.849

9.  Enhanced open-circuit voltage of PbS nanocrystal quantum dot solar cells.

Authors:  Woojun Yoon; Janice E Boercker; Matthew P Lumb; Diogenes Placencia; Edward E Foos; Joseph G Tischler
Journal:  Sci Rep       Date:  2013       Impact factor: 4.379

10.  High Performance PbS Quantum Dot/Graphene Hybrid Solar Cell with Efficient Charge Extraction.

Authors:  Byung-Sung Kim; Darren C J Neo; Bo Hou; Jong Bae Park; Yuljae Cho; Nanlin Zhang; John Hong; Sangyeon Pak; Sanghyo Lee; Jung Inn Sohn; Hazel E Assender; Andrew A R Watt; SeungNam Cha; Jong Min Kim
Journal:  ACS Appl Mater Interfaces       Date:  2016-05-23       Impact factor: 9.229

View more
  2 in total

1.  Ligand-Assisted Formation of Graphene/Quantum Dot Monolayers with Improved Morphological and Electrical Properties.

Authors:  Aleksandr P Litvin; Anton A Babaev; Peter S Parfenov; Aliaksei Dubavik; Sergei A Cherevkov; Mikhail A Baranov; Kirill V Bogdanov; Ivan A Reznik; Pavel O Ilin; Xiaoyu Zhang; Finn Purcell-Milton; Yurii K Gun'ko; Anatoly V Fedorov; Alexander V Baranov
Journal:  Nanomaterials (Basel)       Date:  2020-04-11       Impact factor: 5.076

2.  Beyond Charge Transfer: The Impact of Auger Recombination and FRET on PL Quenching in an rGO-QDs System.

Authors:  Anton A Babaev; Anastasiia V Sokolova; Sergei A Cherevkov; Kevin Berwick; Alexander V Baranov; Anatoly V Fedorov; Aleksandr P Litvin
Journal:  Nanomaterials (Basel)       Date:  2021-06-21       Impact factor: 5.076

  2 in total

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