Literature DB >> 33863720

Nonreciprocal second harmonic generation in a magnetoelectric material.

Shingo Toyoda1, Manfred Fiebig2,3, Taka-Hisa Arima2,4, Yoshinori Tokura2,5,6, Naoki Ogawa2,6,7.   

Abstract

Mirror symmetries are of particular importance because they are connected to fundamental properties and conservation laws. Spatial inversion and time reversal are typically associated to charge and spin phenomena, respectively. When both are broken, magnetoelectric cross-coupling can arise. In the optical regime, a difference between forward and backward propagation of light may result. Usually, this nonreciprocal response is small. We show that a giant nonreciprocal optical response can occur when transferring from linear to nonlinear optics, specifically second harmonic generation (SHG). CuB2O4 exhibits SHG transmission changes by almost 100% upon reversal of a magnetic field of just ±10 mT. The observed nonreciprocity results from an interference between magnetic-dipole and electric-dipole SHG. Although the former is inherently weaker than the latter, a resonantly enhanced magnetic-dipole transition has a comparable amplitude as a nonresonant electric-dipole transition, thus maximizing the nonreciprocity. Multiferroics and magnetoelectrics are an obvious materials platform to exhibit nonreciprocal nonlinear optical functionalities.
Copyright © 2021 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution License 4.0 (CC BY).

Entities:  

Year:  2021        PMID: 33863720      PMCID: PMC8051877          DOI: 10.1126/sciadv.abe2793

Source DB:  PubMed          Journal:  Sci Adv        ISSN: 2375-2548            Impact factor:   14.136


INTRODUCTION

Unidirectional manipulation of photons is a key issue in modern information technology as exemplified by optical isolators used for lasers and in optical networks. While conventional optical isolators are the composite of a magneto-optical medium, magnets, and polarizers, recent studies show that such a one-way flow of photons can be achieved in single bulk magnetoelectric materials, where time-reversal and space-inversion symmetries are simultaneously broken. Such nonreciprocal propagation occurs because the optical properties change with the reversal of the propagation direction of light (–). In the past, a wide variety of nonreciprocal phenomena has been found in linear optics, including the absorption (–), emission (, ), refraction (), and diffraction () of light. In stark contrast, there have been few reports on nonreciprocal phenomena in the nonlinear optical regime (, ), although they are naturally expected (). Second harmonic generation (SHG), which denotes the frequency-doubling of a light wave in a material, is one of the simplest nonlinear optical processes. SHG can be classified into two types depending on its origin. One is electric-dipole (ED)–SHG, in which the frequency-doubled electric polarization in a material is considered, while the other is magnetic-dipole (MD)-SHG, which refers to a frequency-doubled magnetization (Fig. 1A). Here, we neglect electric quadrupole contributions because they are forbidden for the d-d transition of Cu2+ holes in CuB2O4 (). The source term of these SHG processes is given by ()
Fig. 1

Schematic illustration of conventional SHG and nonreciprocal SHG.

(A) SHG intensity generally remains unchanged after reversal of the propagation direction of the light because the SHG signal originates from either ED-SHG or MD-SHG. The electronic transition processes for resonant MD-SHG and nonresonant ED-SHG are sketched on the right-hand side. (B) When the ED-SHG interferes with MD-SHG, the total intensity can depend on the propagating direction of the light. In particular, when the ED-SHG and the MD-SHG yield become similar in both amplitude and phase, their constructive interference substantially enhances the total SHG signal for light propagating in one particular direction. The situation reverses for the opposite propagation direction due to destructive interference, resulting in the suppression or even extinction of the SHG output. (C) Crystal and magnetic structures of CuB2O4 in the canted antiferromagnetic phase, projected onto the (001) plane. The red arrows indicate the magnetic moments of the Cu2+ holes in a magnetic field μ0H ⊥ [001]. (D) Magnetic phase diagram of CuB2O4 within a magnetic field along the ab plane. Canted AFM and PM represent the canted antiferromagnetic and paramagnetic phases, respectively.

Schematic illustration of conventional SHG and nonreciprocal SHG.

(A) SHG intensity generally remains unchanged after reversal of the propagation direction of the light because the SHG signal originates from either ED-SHG or MD-SHG. The electronic transition processes for resonant MD-SHG and nonresonant ED-SHG are sketched on the right-hand side. (B) When the ED-SHG interferes with MD-SHG, the total intensity can depend on the propagating direction of the light. In particular, when the ED-SHG and the MD-SHG yield become similar in both amplitude and phase, their constructive interference substantially enhances the total SHG signal for light propagating in one particular direction. The situation reverses for the opposite propagation direction due to destructive interference, resulting in the suppression or even extinction of the SHG output. (C) Crystal and magnetic structures of CuB2O4 in the canted antiferromagnetic phase, projected onto the (001) plane. The red arrows indicate the magnetic moments of the Cu2+ holes in a magnetic field μ0H ⊥ [001]. (D) Magnetic phase diagram of CuB2O4 within a magnetic field along the ab plane. Canted AFM and PM represent the canted antiferromagnetic and paramagnetic phases, respectively. The SHG intensity on the detector will be , indicating that the intensity changes with the reversal of when both processes coexist and interfere with each other (Fig. 1B). Obviously, for a given , if the first term (ED-SHG) has the same amplitude and phase as the second term (MD-SHG), the total intensity I(2ω) is enhanced by constructive interference, but can be extinguished just by the reversal of . Note that we can also control the SHG interference by reversing or , which, for example, can be realized by the reversal of the static magnetization of the system with an external magnetic field. In the present study, ED-SHG is an i-type tensor and MD-SHG is a c-type tensor (see the Supplementary Materials); the latter changes sign under the reversal of the magnetic field. However, experimental realization of such a huge modulation of the SHG yield (i.e., nonreciprocal SHG of ~100% efficiency) remains a challenge. ED-SHG usually dominates MD-SHG because in the expansion of the electromagnetic field potential, with respect to , the two terms are of zeroth and first order, respectively. Therefore, MD-SHG is less commonly observed because usually it is visible in centrosymmetric materials only, where ED-SHG is forbidden or permitted only through inversion symmetry–breaking magnetically ordered phases (, , ). This unbalance in the ED and MD transition dipole moments hampers the manifestation of large optical nonreciprocity. Here, we report on the observation of nonreciprocal SHG in noncentrosymmetric CuB2O4, where the SHG intensity changes by 97% upon the reversal of an external magnetic field. The effect is most pronounced for light around 1.4 eV, which is resonant with the lowest intra-atomic d-d transition of the Cu2+ ions. We found that the MD-SHG process is resonantly enhanced at this d-d transition, in contrast to the ED-SHG process. The resultant resonant MD-SHG and nonresonant ED-SHG contributions are of the same order of magnitude, yielding a large nonreciprocal signal via their interference. CuB2O4 crystalizes in a noncentrosymmetric tetragonal structure with the space group symmetry I2d (point group 2m) (). The Cu2+ (d, S = 1/2) ions occupy two inequivalent sites denoted as A and B (Fig. 1C), where Cu2+ ions on the A site are responsible for the linear nonreciprocal optical properties (). This material undergoes successive magnetic transitions at TN = 21 K and T* = 9 K (, ) (see Fig. 1D). Below T*, magnetic moments at both A and B sites exhibit incommensurate helical order. Between TN and T*, magnetic moments of the Cu2+ ions show an easy-plane-type canted antiferromagnetic order, as schematically shown in Fig. 1D. By the application of an in-plane external magnetic field of the order of 10 mT, we can align the in-plane magnetization and obtain a single-domain state. In the canted antiferromagnetic phase, where the time-reversal and space-inversion symmetries are simultaneously broken, CuB2O4 shows a magnetoelectric effect, which is explained by the modification of the metal-ligand hybridization with the Cu2+ moments (, ). The magnetic-field-induced electric polarization appears along the c axis (P ∝ sin 2θ) in an external magnetic field in the ab plane, where θ denotes the angle between the crystal’s [010] axis and the direction of the magnetic field ().

RESULTS AND DISCUSSION

Figure 2A shows the experimental setup to detect the nonreciprocal SHG signal. The SHG intensity was measured in a transmission geometry in Voigt configuration (external magnetic field ). The light source was a regenerative amplifier, which produced 190-fs laser pulses at 6 kHz. The energy of the fundamental light was tuned by an optical parametric amplifier (OPA) to ℏω = 0.703 eV, for which the SHG energy 2ℏω = 1.406 eV was resonant with the d-d transition of Cu2+ holes between the d and d orbitals (, , –). Note that the spectral width of the electronic transition of ~1 meV was one order of magnitude smaller than the energy spread of the incident laser pulse. We used a spectrometer to resolve the fine structure of the SHG spectrum. The thickness of the sample was 50 μm, with the widest crystal faces exhibiting (100) orientation. The sample was mounted on a copper holder, which could rotate around the c axis by the angle θ (Fig. 2A). We tilt the sample around the c axis, and therefore, the electric field of the fundamental light Eω is always in the ab plane and that of the SHG field E2ω remains parallel to the c axis. This suggests that the refractive index and, hence, the phase matching condition are not affected by the tilting of the sample because the crystal structure of CuB2O4 is uniaxial. Figure 2 (B to D) shows the tilt-angle dependence of the SHG contribution polarized along the c axis in a magnetic field μ0H of 50 mT. When the sample is not tilted (θ = 0°), the electric polarization is absent for H parallel to the [010] axis (). In this situation, the SHG signal originates solely from MD-SHG, because the complementing ED-SHG contribution with Eω ∥ [010] is zero (see the Supplementary Materials). This results in the absence of ED-MD interference and, hence, of a change in SHG intensity upon the reversal of the magnetic field (see Fig. 2C). When the sample is tilted around the c axis, the [100] component of the electric field of light and the associated ED-SHG from χ and χ become allowed. Figure 2 (B and D) shows the SHG spectra in the magnetic field μ0H = 50 mT for the sample tilted by θ = −15° and θ = +15°, respectively. The associated SHG spectra exhibit drastic changes with the reversal of magnetic field, demonstrating the nonlinear optical nonreciprocity. Notably, the SHG spectra show a Fano resonance–like asymmetric shape. A Fano resonance is a signature of the interference between a resonant process and a nonresonant background (, ), in the present case represented by MD-SHG for the former and ED-SHG for the latter.
Fig. 2

Experimental observation of nonreciprocal SHG.

(A) Optical setup to detect nonreciprocal SHG emission. The light pulses (photon energy ℏω = 0.703 eV) from an optical parametric amplifier were focused onto a CuB2O4 single crystal after setting their polarization to Eω ⊥ [001] and Hω ∥ [001] by using a Glan-laser prism. The sample was tilted around the c axis by an angle θ. SHG spectra for linearly polarized light (E2ω ∥ [001] and H2ω ⊥ [001]) were measured in a transmission geometry with a spectrometer and a charge-coupled device camera as detector. An external magnetic field was applied in the ab plane normal to the propagation direction of the light (Voigt geometry). (B to D) Spectral dependence of the SHG for tilt angles (B) θ = −15°, (C) θ = 0°, and (D) θ = +15°, measured at T = 12 K. Green and yellow lines show the spectra in a magnetic field of μ0H = −50 mT and μ0H = +50 mT, respectively.

Experimental observation of nonreciprocal SHG.

(A) Optical setup to detect nonreciprocal SHG emission. The light pulses (photon energy ℏω = 0.703 eV) from an optical parametric amplifier were focused onto a CuB2O4 single crystal after setting their polarization to Eω ⊥ [001] and Hω ∥ [001] by using a Glan-laser prism. The sample was tilted around the c axis by an angle θ. SHG spectra for linearly polarized light (E2ω ∥ [001] and H2ω ⊥ [001]) were measured in a transmission geometry with a spectrometer and a charge-coupled device camera as detector. An external magnetic field was applied in the ab plane normal to the propagation direction of the light (Voigt geometry). (B to D) Spectral dependence of the SHG for tilt angles (B) θ = −15°, (C) θ = 0°, and (D) θ = +15°, measured at T = 12 K. Green and yellow lines show the spectra in a magnetic field of μ0H = −50 mT and μ0H = +50 mT, respectively. Next, we investigate the temperature dependence of the nonreciprocal SHG signal for the sample tilted by θ = + 15°. Figure 3 (A to C) shows the SHG spectra at different temperatures in a magnetic field μ0H = ± 50 mT. Whereas the nonreciprocal signal is clearly observed in the canted antiferromagnetic phase, it disappears in the helical and the paramagnetic phases, where the time-reversal symmetry is preserved. This result confirms that the space- and time-inversion symmetry breaking is essential to the emergence of the nonreciprocal behavior. The SHG signal shows a broad and featureless spectral distribution in the helical and paramagnetic phases, which, referring to our analysis of the Fano-resonant behavior above, points to its ED-SHG origin. ED-SHG can be expressed by , where and are the SHG susceptibility and the electric field of the incident light, respectively. The 2m point-group symmetry allows the SHG tensor components χ and χ for SHG light polarized along the c axis (, ), which becomes accessible in our experiment, where Eω ⊥ [001] and E2ω ∥ [001]. The ED-SHG polarization associated to χ and χ changes sign with the change of rotation of the sample from +15° to −15° because of the associated reversal in the sign of . This explains the sign reversal in the magnetic-field dependence in Fig. 2 (B and D).
Fig. 3

Temperature dependence of nonreciprocal SHG.

SHG spectra measured at (A) 6 K in the helical, (B) 12 K in the canted AFM, and (C) 23 K in the paramagnetic phase with the sample tilted by θ = +15°. Green and yellow lines indicate the spectra in a magnetic field of μ0H = −50 mT and μ0H = +50 mT, respectively. The magnetic field is applied in the ab plane in the Voigt geometry (see Fig. 2A). (D to F) Temperature dependence of the SHG intensities (D) below the resonance energy (Δ = −0.7 meV), (E) at the resonance energy (Δ = 0 meV), and (F) above the resonance energy (Δ = +1.1 meV), where Δ denotes the energy deviation from the resonance energy 1.4066 eV. These photon energies are indicated by the dashed lines in (B).

Temperature dependence of nonreciprocal SHG.

SHG spectra measured at (A) 6 K in the helical, (B) 12 K in the canted AFM, and (C) 23 K in the paramagnetic phase with the sample tilted by θ = +15°. Green and yellow lines indicate the spectra in a magnetic field of μ0H = −50 mT and μ0H = +50 mT, respectively. The magnetic field is applied in the ab plane in the Voigt geometry (see Fig. 2A). (D to F) Temperature dependence of the SHG intensities (D) below the resonance energy (Δ = −0.7 meV), (E) at the resonance energy (Δ = 0 meV), and (F) above the resonance energy (Δ = +1.1 meV), where Δ denotes the energy deviation from the resonance energy 1.4066 eV. These photon energies are indicated by the dashed lines in (B). Figure 3 (D to F) shows the temperature dependence of the SHG intensity at photon energies of 1.4059, 1.4066, and 1.4077 eV, which are indicated by the dashed lines in Fig. 3B. At off-resonant photon energies (Fig. 3, D and F), the ED-MD interference leads to a pronounced nonreciprocal signal in the canted antiferromagnetic phase. In contrast, the SHG intensity shows little change with the reversal of magnetic field at the resonant energy (see Fig. 3E). Here, we defined the resonance energy as the peak position of the pure MD-SHG spectrum in Fig. 2C, which is measured at T = 12 K without tilting the sample (θ = 0°). The remaining dependence on the direction of the magnetic field is caused by the small shift of the resonance energy experiences with temperature (). Thus, at the resonance energy, the nonreciprocity of the signal in the canted antiferromagnetic phase formally disappears. To further elucidate the origin of the nonreciprocal signal, we investigate the SHG spectrum across a broader spectral range (see Fig. 4A). The ED- and MD-SHG spectra can be measured separately by properly selecting the temperature and the tilt angle θ. Pure ED-SHG is obtained at 25 K in the paramagnetic phase for the sample tilted by 5°, whereas pure MD-SHG is measured at 12 K in the canted antiferromagnetic phase without tilting the sample. The ED-SHG signal does not show any peak at the d-d transition of 1.4 eV and increases with small variations as the photon energy is raised. We ascribe this to the nonresonant virtual excitation of the charge transfer transition between Cu2+ and O2− ions. On the other hand, the MD-SHG signal is resonantly enhanced at the d-d transitions between the d and d orbitals of Cu2+ holes at the A site. These results explain why MD- and ED-SHG can be of comparable amplitude, which then leads to the strongly nonreciprocal signal.
Fig. 4

Origin of nonreciprocal SHG.

(A) Spectra of the MD-SHG and ED-SHG contributions. The MD-SHG spectrum is obtained without tilting the sample (θ = 0°) in a magnetic field of μ0H = +50 mT along the [010] axis at T = 12 K in the canted AFM phase. The ED-SHG spectrum is measured at T = 25 K in the paramagnetic phase for a tilt angle θ = 5°. Both spectra are normalized to the maximum values. (B) Schematic of the origin of maximum nonreciprocity. Blue and red lines (dashed lines) show the amplitude (phase) of ED-SHG and MD-SHG contributions, respectively. The resonantly enhanced MD-SHG contribution shows a comparable amplitude as the nonresonant ED-SHG contribution and has the same amplitude and the same (the opposite) phase below (above) the resonant photon energy as indicated by dotted circles. (C and D) Magnetic-field dependence of the SHG intensity (T = 12 K) for the sample tilted by +15° (C) below (1.4059 eV) and (D) above (1.4077 eV) the resonance energies. These photon energies are indicated by dashed circles in (B).

Origin of nonreciprocal SHG.

(A) Spectra of the MD-SHG and ED-SHG contributions. The MD-SHG spectrum is obtained without tilting the sample (θ = 0°) in a magnetic field of μ0H = +50 mT along the [010] axis at T = 12 K in the canted AFM phase. The ED-SHG spectrum is measured at T = 25 K in the paramagnetic phase for a tilt angle θ = 5°. Both spectra are normalized to the maximum values. (B) Schematic of the origin of maximum nonreciprocity. Blue and red lines (dashed lines) show the amplitude (phase) of ED-SHG and MD-SHG contributions, respectively. The resonantly enhanced MD-SHG contribution shows a comparable amplitude as the nonresonant ED-SHG contribution and has the same amplitude and the same (the opposite) phase below (above) the resonant photon energy as indicated by dotted circles. (C and D) Magnetic-field dependence of the SHG intensity (T = 12 K) for the sample tilted by +15° (C) below (1.4059 eV) and (D) above (1.4077 eV) the resonance energies. These photon energies are indicated by dashed circles in (B). Next, we discuss the observed spectral shape. Figure 4B illustrates a schematic of the MD- and ED-SHG susceptibilities around the resonance energy. Both the amplitude and phase of the ED-SHG contribution are almost unchanged around the resonance, because its origin is unrelated to the d-d transitions. On the other hand, phase and amplitude of the MD-SHG contribution show pronounced changes in this region. The MD-SHG wave experiences a 180° phase shift across the resonance energy and the amplitude sharply increases to show the resonance peak. Therefore, the ED- and MD-SHG light fields interfere constructively (destructively) at photon energies below (above) the MD resonance, which explains the sign change of the nonreciprocal effect across the resonance. Notably, the amplitudes of the MD- and ED-SHG contributions become the same at the positions indicated by dotted circles in Fig. 4B. Here, we expect a maximum of the nonreciprocal behavior in the SHG response by the ED-MD interference. Figure 4C shows the magnetic-field dependence of the SHG intensity at 1.4059 eV (that is, slightly below the resonance). The SHG intensity almost disappears for the negative magnetic field, whereas a strong SHG signal shows up for the same yet positive magnetic-field value, indicating that MD-SHG and ED-SHG light waves of the same amplitude interfere. Slightly above the resonance (1.4077 eV), the situation is reversed, as shown in Fig. 4D. This is in excellent agreement with our model of the interference between resonant MD and nonresonant ED components. The asymmetry in the extinction of the SHG signal above (Fig. 4D) and below (Fig. 4C) the resonance photon energy can be explained by the phase between the ED-SHG and MD-SHG contributions. Below the resonance, ED- and MD-SHG are in phase, which allows for near-perfect nonreciprocity. However, since the phase change of the MD-SHG contribution across the resonance is less than 180°, the phase between ED- and MD-SHG slightly differs above the resonance, resulting in the imperfect interference and, hence, extinction of the SHG in Fig. 4D. To conclude, we experimentally demonstrate the nearly perfectly unidirectional propagation of SHG light waves in CuB2O4. The SHG intensity changes by almost 100% upon the reversal of a small magnetic field of 10 mT. The observed nonreciprocity originates from the interference between ED-SHG and MD-SHG. The MD-SHG process is resonantly enhanced at the d-d transition of Cu2+ holes and becomes comparable in amplitude to the nonresonant broadband ED-SHG process. This work establishes that nonreciprocal wave-propagation effects are not limited to linear optics but are also present as nonlinear optical effects and exhibit a magnitude up to the optical-diode-like unidirectional propagation of light.

MATERIALS AND METHODS

Sample preparation

Single crystals of CuB2O4 were grown by a flux method. Powders of CuO (7.250 g), B2O3 (15.228 g), and LiCO3 (4.041 g) were mixed. They were heated in air at 1020°C, subsequently cooled down to 800°C at a rate of 1.1°C/hour, and then cooled to room temperature at a rate of 16°C/hour. A single crystal was oriented by using a Laue x-ray diffractometer and cut into thin plates of 50 μm thickness with the widest (100) faces. The sample surfaces were specularly polished by alumina lapping films. Note added in proof: We have recently become aware of the related work by J. Mund et al. (arXiv: 2005.05393 2020).
  10 in total

1.  Second harmonic generation in the centrosymmetric antiferromagnet NiO.

Authors:  M Fiebig; D Fröhlich; T Lottermoser; V V Pavlov; R V Pisarev; H J Weber
Journal:  Phys Rev Lett       Date:  2001-09-04       Impact factor: 9.161

2.  Second harmonic generation and magnetic-dipole-electric-dipole interference in antiferromagnetic Cr2O3.

Authors: 
Journal:  Phys Rev Lett       Date:  1994-10-10       Impact factor: 9.161

3.  One-Way Transparency of Light in Multiferroic CuB(2)O(4).

Authors:  S Toyoda; N Abe; S Kimura; Y H Matsuda; T Nomura; A Ikeda; S Takeyama; T Arima
Journal:  Phys Rev Lett       Date:  2015-12-29       Impact factor: 9.161

4.  Enhanced optical magnetoelectric effect in a patterned polar ferrimagnet.

Authors:  N Kida; Y Kaneko; J P He; M Matsubara; H Sato; T Arima; H Akoh; Y Tokura
Journal:  Phys Rev Lett       Date:  2006-04-26       Impact factor: 9.161

5.  Giant nonreciprocal second-harmonic generation from antiferromagnetic bilayer CrI3.

Authors:  Zeyuan Sun; Yangfan Yi; Tiancheng Song; Genevieve Clark; Bevin Huang; Yuwei Shan; Shuang Wu; Di Huang; Chunlei Gao; Zhanghai Chen; Michael McGuire; Ting Cao; Di Xiao; Wei-Tao Liu; Wang Yao; Xiaodong Xu; Shiwei Wu
Journal:  Nature       Date:  2019-07-31       Impact factor: 49.962

6.  Antiferromagnetic Dichroism in a Complex Multisublattice Magnetoelectric CuB(2)O(4).

Authors:  K N Boldyrev; R V Pisarev; L N Bezmaternykh; M N Popova
Journal:  Phys Rev Lett       Date:  2015-06-19       Impact factor: 9.161

7.  Nonreciprocal Refraction of Light in a Magnetoelectric Material.

Authors:  S Toyoda; N Abe; T Arima
Journal:  Phys Rev Lett       Date:  2019-08-16       Impact factor: 9.161

8.  Enhanced directional dichroism of terahertz light in resonance with magnetic excitations of the multiferroic Ba2CoGe2O7 oxide compound.

Authors:  I Kézsmárki; N Kida; H Murakawa; S Bordács; Y Onose; Y Tokura
Journal:  Phys Rev Lett       Date:  2011-02-04       Impact factor: 9.161

9.  Magnetic-field induced second harmonic generation in CuB2O4.

Authors:  R V Pisarev; I Sänger; G A Petrakovskii; M Fiebig
Journal:  Phys Rev Lett       Date:  2004-07-15       Impact factor: 9.161

Review 10.  Nonreciprocal responses from non-centrosymmetric quantum materials.

Authors:  Yoshinori Tokura; Naoto Nagaosa
Journal:  Nat Commun       Date:  2018-09-14       Impact factor: 14.919

  10 in total
  1 in total

1.  Magnetoelectric dissociation of Alzheimer's β-amyloid aggregates.

Authors:  Jinhyeong Jang; Chan Beum Park
Journal:  Sci Adv       Date:  2022-05-11       Impact factor: 14.957

  1 in total

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