三个镧系金属氮氧自由基配合物的合成、结构及磁性

胡鹏 高媛媛 肖凤仪 邓肖娟 黄国洪 张淼 苏芬 王莉娜

引用本文: 胡鹏, 高媛媛, 肖凤仪, 邓肖娟, 黄国洪, 张淼, 苏芬, 王莉娜. 三个镧系金属氮氧自由基配合物的合成、结构及磁性[J]. 无机化学学报, 2017, 33(1): 33-40. doi: 10.11862/CJIC.2017.019 shu
Citation:  HU Peng, GAO Yuan-Yuan, XIAO Feng-Yi, DENG Xiao-Juan, HUANG Guo-Hong, ZHANG Miao, SU Feng, WANG Li-Na. Three Lanthanide Nitronyl Nitroxide Radical Compounds: Synthesises, Structures and Magnetic Properties[J]. Chinese Journal of Inorganic Chemistry, 2017, 33(1): 33-40. doi: 10.11862/CJIC.2017.019 shu

三个镧系金属氮氧自由基配合物的合成、结构及磁性

    通讯作者: 高媛媛,E-mail:gaoyuanyuan@imut.edu.cn;肖凤仪,E-mail:michellexiao2014@163.com
  • 基金项目:

    内蒙古自然科学基金 2016BS0206

    广东省大创项目 201510580046, 201610580046

    广东省教育厅青年创新人才项目 CQ2014064

摘要: 以氮氧自由基为配体, 合成了3例氮氧自由基-稀土三自旋单核配合物[Ln(hfac)3(NIT-Ph-4-Br)2](Ln=Gd(1), Tb(2), Dy(3), hfac=六氟乙酰丙酮, NIT-Ph-4-Br=4, 4, 5, 5-四甲基-2-(4'-溴)-咪唑啉-3-氧化-1-氧基自由基。单晶结构分析表明3个配合物均属单斜晶系P21/c空间群, 配合物中的Ln (III) 离子为八配位模式, 并且拥有相似的自由基-稀土-自由基单核结构。对配合物的磁性测试结果表明, 配合物1中自由基与Gd (III) 离子之间存在着铁磁相互作用, 自由基与自由基之间存在着反铁磁相互作用;配合物2, 3中, 稀土离子与自由基之间存在弱的反铁磁相互作用

English

  • 

    0    Introduction

    Molecular magnetic materials have attracted scientists attention in the past two decades due to their potential applications in high-density magnetic memories, quantum computing devices and molecular spintronics[1-6].Among the different approaches to prepare molecular magnetic materials the metal-radical strategy that consists of matching paramagnetic organic molecules with transition metal complex gives rise to a variety of compounds with different structural and magnetic dimensionalities.Up to now various organic paramagnetic molecules such as verdazyl, semiquinone and nitronyl nitroxide (NIT) radicals have been widely studied in the field of molecular magnetism[7-14].Among them, nitronyl nitroxide radicals have received noteworthy attention because this type of radicals can act as bidentate ligands with identical N-O coordination groups.Besides, nitronyl nitroxide family of radicals are relatively stable and easy to obtain derivatives with substituents containing donor atoms.However, NIT radicals are poorly donating ligands, thus utilization of strong electron withdrawing coligands such as hexafluoroacetylacetonate (hfac) [15-17] and trifluoroacetylacetonate (tfac) [18-20] in the metal to improve weak coordination ability is necessary.The steric demand of hfac or tfac restrict the dimensionality of the resulting metal-radical complex.Thus numbers of zero- and one-dimensional complexes were prepared by this strategy.Recently lanthanide ions based molecular magnetic materials have been extensively studied because lanthanide ions such as terbium (III) and dysprosium (III) are good candidates for the construction of SMMs due to their significant magnetic anisotropy arising from the large, unquenched orbital angular momentum[21-24].

    To further study the magnetic properties of NIT radical-lanthanide compounds, in this paper we report a nitronyl nitroxide radical (Scheme 1) and its corres-ponding Ln-nitronyl nitroxide compounds [Ln (hfac) 3 (NIT -Ph-4-Br) 2] (Ln=Gd (1) , Tb (2) , Dy (3) , hfac=hexafluo-roacetylacetonate, NIT-Ph-4-Br=2- (4′-bromine phenyl) -4, 4, 5, 5-tetramethyl-imidazoline-1-oxyl-3-oxide) , their crystal structures and magnetic properties were described in detail.

    图Scheme 1 Molecular structure of NIT-Ph-4-Br

    1    Experimental

    1.1    Materials and measurements

    All reagents and solvents were purchased from Aladdin and used without further purification.The radical ligand NIT-Ph-4-Br was synthesized according to literature [25].Elemental analyses (C, H and N) were performed on a Perkin-Elmer 240C elemental analyzer.IR spectra were recorded on a Nicolet IS10 IR spectrometer using KBr pellets in the range of 4 000~500 cm-1.The magnetic measurements were carried out with MPMSXL-7 SQUID magnetometer.Diamagnetic corrections were made with Pascals constants for all the constituent atoms.

    1.2    Synthesis of complex 1

    A suspension of Gd (hfac) 3·2H2O (0.05 mmol) in 15 mL dry boiling heptane was heated to reflux for about 1 h.Then the solution was cooled to 65 ℃, a solution of NIT-Ph-4-Br 0.1 mmol) in 2 mL of CHCl3 was added.The resulting solution was stirred for about 2 min and then cooled to room temperature.The filtrate was allowed to stand at room temperature for slow evaporation.Slow evaporation of the final solution for about four days yielded dark-blue block crystals suitable for single-crystal X-ray analysis.Yield: 43.5% based on rare-earth.Elemental analysis calculated for C41H35Br2GdF18N4O10 (%) : C: 35.10; H: 2.51; N: 3.99.Found (%) : C: 33.88; H: 2.79; N: 4.12.FTIR (KBr, cm-1) : 1 656 (s) , 1527 (w) , 1385 (w) , 1350 (w) , 1255 (s) , 1198 (s) , 1096 (s) , 796 (w) , 660 (w) .

    1.3    Synthesis of Complex 2

    Complex 2 was synthesized using the same procedure for complex 1 with Tb (hfac) 3·2H2O instead of Gd (hfac) 3·2H2O.Yield: 47.9%.Elemental analysis calculated for C41H35Br2GdF18N4O10 (%) : C: 35.06; H: 2.51; N: 3.99.Found (%) : C: 35.69; H: 2.57; N: 4.11.FTIR (KBr, cm-1) : 1 655 (s) , 1599 (w) , 1528 (w) , 1399 (w) , 1352 (w) , 1255 (s) , 1199 (s) , 1096 (w) , 797 (w) , 660 (w) .

    1.4    Synthesis of Complex 3

    Complex 3 was synthesized using the same procedure for complex 1 with Dy (hfac) 3·2H2O instead of Gd (hfac) 3·2H2O.Yield: 42.9%.Elemental analysis calculated for C41H35Br2GdF18N4O10 (%) : C: 34.97; H: 2.51; N: 3.98.Found (%) : C: 34.91; H: 2.77; N: 3.87 (%) .FTIR (KBr, cm-1) : 1 655 (s) , 1600 (w) , 1527 (w) , 1399 (w) , 1352 (w) , 1255 (s) , 1199 (s) , 1096 (w) , 798 (w) , 662 (w) .

    1.5    X-ray Crystallographic Study

    表1 Crystal data and structure refinement for 1, 2 and 3
    表1 Crystal data and structure refinement for 1, 2 and 3

    The crystal structure data of complexes 1, 2 and 3 were collected using a Rigaku Saturn CCD diffractometer equipped with graphite-monochromated Mo Kα radiation (λ=0.071 073 nm) .Crystal size of complexes 1, 2 and 3 are 0.15 mm×0.12 mm× 0.11 mm, 0.17 mm×0.11 mm×0.10 mm, 0.12 mm×0.11 mm×0.10 mm respectively.The structures were solved by the direct methods with SHELXS-97[26] and refined by full-matrix least-squares methods on F 2 with SHELXL-97 program package[27].Anisotropic thermal parameters were assigned to all non-hydrogen atoms.The hydrogen atoms were set in calculated positions and refined as riding atoms with a common fixed isotropic thermal parameter.The details of the crystal parameters, data collection, and refinements for these complexes were listed in Table 1 and 2.

    CCDC: 1496095, 1; 1496096, 2; 1496097, 3.

    表2 Selected bond lengths (nm) and angles (°) for 1 , 2 and 3
    表2 Selected bond lengths (nm) and angles (°) for 1 , 2 and 3

    2    Results and discussion

    2.1    Crystal structures of complexes 1~3

    As is shown in Fig.1, Compound 1 is a mononuclear coordination compound crystallizing in the monoclinic space group P21/c with Z=4.The central Gd (III) ion is eight coordinated in slightly distorted triangular dodecahedral GdO8 geometry completed by two non-bridged NO groups from two separate organic radicals and three bischelate hfac- anions.The distances of Gd-O bonds range from 0.235 6 (4) to 0.240 7 (4) nm.The coordinated N (3) -O (2) and N (2) -O (3) bond lengths of nitronyl nitroxide radicals are 0.130 7 (5) nm and 0.129 9 (6) nm respec-tively and the uncoordinated N (4) -O (1) and N (1) -O (4) bond lengths are 0.126 4 (6) nm and 0.126 0 (7) nm respectively, which are comparable to those of reported tri-spin radical-Ln (III) -radical complexes[28-33].The nearest Gd…Gd distance between adjacent molecules is 1.071 6 (5) nm (Fig.1) .

    图1 Molecular structure (left) with thermal ellipsoids drawn at 30% probability and crystal packing diagram (right) of complex 1

    Compound 3 is also isostructural to compound 1 and the bond lengths of Dy-O are in the range of 0.233 0 (2) ~0.238 6 (2) nm, which are a little bit shorter than the bond lengths of Gd-O and Tb-O.The nearest Dy…Dy distance between adjacent molecules is 1.060 8 (3) nm.

    Compound 2 is isostructural to compound 1 and the bond lengths of Tb-O are in the range of 0.233 8 (5) ~0.239 4 (5) nm, which are a little bit shorter than the bond lengths of Gd-O.The nearest Tb…Tb distance between adjacent molecules is 1.081 2 (4) nm.

    2.2    Magnetic property of complex 1

    图2 Temperature dependence of χMT (left) and field dependence of magnetization at 2.0 K (right) for complex 1

    The temperature dependence of the magnetic susceptibilities 1, 2, and 3 were measured from 300 to 2.0 K in an applied field of 1 kOe, and the magnetic behaviors of complex 1 are shown in Fig.2.At 300 K, the χMT value is 8.61 cm3·K·mol-1.The values are in good agreement with the theoretical value of 8.63 cm3·K·mol-1 (uncoupled one Gd (III) ion ( f 7 electron confi-guration, χMT=7.88 cm3·K·mol-1) plus two organic radicals (S=1/2, χMT=0.375 cm3·K·mol-1) ) .Upon cooling, the χMT value of complex 1 increase steadily to a maximum of 10.03 cm3·K·mol-1 at 13.9 K, afterward decreases to 9.13 cm3·K·mol-1 at 2.0 K.

    Furthermore, the field dependence of magnetization of complex 1 has been determined at 2 K in the range of 0~70 kOe (Fig.2) .Upon increasing in the applied field, M increases up to 8.87Nβ at 70 kOe, which corresponds well to the value expected for a ground state with a spin multiplicity of S=9/2 in the case of one Gd (III) and two radical ferromagnetically coupled.At the lower fields the value is smaller than the magnetization calculated with the Brillouin function for noncoupled S=7/2 and two S=1/2 spin centers (g=2.0, T=2 K) , which suggests the dominant antiferromagnetic interaction between the Gd (III) ion and the coordinated NIT radical.

    The magnetic interactions between Gd (III) and the radicals can be described by isotropic exchange interaction.Therefore the experimental data for complex 1 can be analyzed with an expression derived from a spin Hamiltonian.Considering the g value range of the radical and Gd (III) ion, we assume that the radical and Gd (III) ion have the same g value.Thus, the variable-temperature magnetic susceptibility data for complex 1 can be analyzed by a theoretical expression (Eq. (2) ) deduced from a spin Hamiltonian (Eq. (1) ) .The JRad-Gd represent the magnetic coupling for the Gd-radical and JRad-Rad for radical-radical, and the zJ′ in Eq. (2) representing the intermolecular interactions.

    The best fitting results give coupling parameters g=1.99, JRad-Gd=2.81 cm-1, JRad-Rad=-10.85 cm-1, zJ′=-0.02 cm-1, R=1.64×10-5, where R is defined as R= ( χM, obs- χM, calc) 2/ ( χM, obs) 2 for complex 1.The positive value of JRad-Gd indicates that there is a weak ferromagnetic interaction between the Gd (III) and the radicals in the molecule.The negative JRad-Rad indicates the antiferromagnetic interaction between the two intramolecular radicals.It is worth noting that the value of zJ′is much smaller than that of JRad-Rad.The obtained J value is comparable with the previously reported Gd (III) -radicals compounds[28-34].

    There are two kinds of magnetic interactions in this radical-Gd (III) -radical complex at the same time.The first one is Gd (III) -radical interaction and the second one is radical-radical interaction.

    2.3    Magnetic properties of complex 2 andcomplex 3

    While for complex 2 (Fig.3) , at 300 K, the χMT value is 12.46 cm3·K·mol-1, and the values are in good agreement with to the theoretical value 12.57 cm3·K·mol-1 in uncoupled system of one Tb (III) ion ( f 9 electron configuration, χMT=11.82 cm3·K· mol-1) and two organic radical (S=1/2, χMT=0.375 cm3·K·mol-1) .Upon cooling, the χMT values of complex 2 maintain a constant behavior down to about 50 K then the value decrease gradually and reach a minimum of 10.37 cm3·K·mol-1.

    图4 Field dependence of magnetization of 2 and 3 at 2.0 K

    In the expression, Δ is the zero-field-splitting parameter, g is the Lande factor, k is the Boltzmann constant, β is the Bohr magneton constant and N is Avogadros constant. The zJ′ parameter based on the molecular field approximation in Eq. (7) is introduced to simulate the magnetic interactions between all the paramagnetic species in the system.Thus the magnetic data can be analyzed by the following approximate treatment of Eq. (4) ~ (7) [29, 31].Giving the best fitting parameters for complex 2 are g=1.49, =3.89×10-2 cm-1, zJ′=1.63×10-2 cm-1, R=1.12×10-5 and for complex 3 are g=1.36, Δ=3.17×10-2 cm-1 , zJ′=2.14×10-2 cm-1, R=1.92×10-5, where R is defined as R= ( χM, obs-χM, calc) 2/ ( χM, obs) 2.The very small positive zJ′ values are indicative of very weak ferromagnetic interaction between Ln (III) ions and the coordinated nitronyl nitroxide, which is consistent with the reported heavy lanthanide-nitronyl nitroxide complexes.

    图3 Temperature dependence of χMT for complex 2 (left) and complex 3 (right)

    Complex 3 shows similar magnetic properties with complex 2 (Fig.3) .At 300 K, the χMT value is 14.31 cm3·K·mol-1, close to the theoretical value of 14.92 cm3·K·mol-1 (one Dy (III) ion ( f 9 electron configuration, χMT=14.17 cm3·K·mol-1) plus two organic radicals (S=1/2, χMT=0.375 cm3·K·mol-1) ) .Upon cooling, the χMT values of complex 3 maintain a constant behavior down to about 50 K then the value decrease gradually and reach a minimum of 8.21 cm3·K·mol-1.

    There is no available expression to determine the magnetic susceptibilities of Ln (III) systems with large anisotropy.To obtain a rough quantitative estimation of the magnetic interaction between Ln (III) and radicals, the magnetic susceptibility χtotal of the complex can be assumed as the sum of χLn of the isolated Dy (III) or Tb (III) ion and χRad of the radical (Eq. (4) ) .The χTb and χDy can be described as Eq. (5) and (6) , respectively.

    The field dependences of magnetizations for complexes 2 and 3 have been determined at 2 K in the range of 0~70 kOe (Fig.4) .Upon increasing in the applied field, M increases up to 5.87Nβ and 7.73Nβ at 70 kOe for 2 and 3, respectively, which is much lower than the saturation value of 11.0Nβ (one Tb (III) ion with g=3/2 and J=6 (9.0Nβ) plus two radicals with g=2.0 and S=1/2) and 12.0Nβ (one Dy (III) ion with g=4/3 and J=15/2 (10.0Nβ) plus two radicals with g=2.0 and S=1/2) .Considering the strong spin-orbit coupling in Ln (III) ions, the large gaps between experimental data and theoretical saturation values for compounds 2 and 3 can be ascribed to the presence of a magnetic anisotropy and/or low-lying excited states in the system.

    2.4    Dynamic magnetic properties for 2 and 3

    Alternating current (ac) susceptibility measurements for 2 and 3 were carried out in the low temperature region under a zero dc field with frequency of 111 and 511Hz.The result (Fig.5) shows that there are no obvious frequency dependent in-phase ( χ′) and out-of-phase ( χ″) signals for both complex 2 and 3, they do not express SMMs behavior at low temperature.

    图5 Temperature dependence of the in-phase and out-of-phase components of ac susceptibility for 2 (left) and 3 (right) in zero dc field with an oscillation of 3.5 Oe

    3    Conclusions

    In conclusion, we report three new complexes based on nitronyl nitroxide radicals and lanthanide ions.These three compounds have similar structures, in which two radical ligands are coordinated to the Ln (III) ions via the oxygen atoms of the nitroxide to form the three spin system.The magnetic studies reveal that ferromagnetic interactions (between the intramolecular Ln and radical) and antiferromagnetic interactions (between the intramolecular radicals) coexist in complex 1.Complexes 2 and 3 show very weak ferromagnetic interaction between Ln (III) ions and the coordinated nitronyl nitroxide.Both complex 2 and 3 do not have SMMs behavior at low temperature, this may due to the small energy barrier which could not prevent the inversion of spin.

    1. [1]

      Kahn O. Molecular Magnetism. New York:VCH Publishers Inc., 1993:1-23

    2. [2]

      Moller S, Perlov C, Jackson W, et al. Nature, 2003,426:166-169 doi: 10.1038/nature02070

    3. [3]

      Kahn M L, Sutter J P, Golhen S, et al. J. Am. Chem. Soc., 2000,122:3413-3421 doi: 10.1021/ja994175o

    4. [4]

      Zhang P, Guo Y N, Tang J K. Coord. Chem. Rev., 2013,257: 1728-1737 doi: 10.1016/j.ccr.2013.01.012

    5. [5]

      Graham M J, Zadrozny J M, Shiddiq M, et al. J. Am. Chem. Soc., 2014,136:7623-7626 doi: 10.1021/ja5037397

    6. [6]

      Bar A K, Pichon C, Gogoi N. Chem. Commun., 2015,51: 3616-3619 doi: 10.1039/C4CC10182K

    7. [7]

      Rinehart J D, Fang M, Evans W J, et al. J. Am. Chem. Soc., 2011,133:14236-14239 doi: 10.1021/ja206286h

    8. [8]

      Rinehart J D, Fang M, Evans W J, et al. Nat. Chem., 2011, 3:538-542 doi: 10.1038/nchem.1063

    9. [9]

      Wang X F, Hu P, Li Y G, et al. Chem.-Asian J., 2015,10: 325-330 doi: 10.1002/asia.v10.2

    10. [10]

      Adugna S, Revunova K, Djukic B, et al. Inorg. Chem. 2010, 49:10183-10190 doi: 10.1021/ic101471u

    11. [11]

      Chernick E T, Casillas R, Zirzlmeier J, et al. J. Am. Chem. Soc., 2015,137:857-863 doi: 10.1021/ja510958k

    12. [12]

      Mailman A, Winter S M, Wong J W L, et al. J. Am. Chem. Soc., 2015,137:1044-1049 doi: 10.1021/ja512235h

    13. [13]

      Brook D J R, Richardson C J, Haller B C, et al. Chem. Commun., 2010,46:6590-6592 doi: 10.1039/c0cc02233k

    14. [14]

      Norel L, Chamoreau L M, Journaux Y, et al. Chem. Commun., 2010,45:2381-2383

    15. [15]

      Murakami R, Ishida T, Yoshii S, et al. Dalton Trans. 2013, 42:13968-13973 doi: 10.1039/c3dt51784e

    16. [16]

      Bernot K, Bogani L, Caneschi A, et al. J. Am. Chem. Soc., 2006,128:7947-7956 doi: 10.1021/ja061101l

    17. [17]

      Zhu M, Hu P, Li Y, et al. Chem. Eur. J., 2014,20:13356-13364 doi: 10.1002/chem.201403581

    18. [18]

      Hu P, Zhang C, Gao Y, et al. Inorg. Chim. Acta, 2013,398: 136-140 doi: 10.1016/j.ica.2012.12.025

    19. [19]

      Mei X L, Ma Y, Li L C. Dalton Trans., 2012,41:505-510 doi: 10.1039/C1DT11795E

    20. [20]

      Mei X L, Liu R N, Wang C. Dalton Trans., 2012,41:2904-2909 doi: 10.1039/c2dt11671e

    21. [21]

      Bernot K, Luzon J, Bogani L, et al. J. Am. Chem. Soc., 2009, 131:5573-5579 doi: 10.1021/ja8100038

    22. [22]

      Liu J L, Wu J Y, Chen Y C, et al. Angew. Chem. Int. Ed., 2014,53:12966-12970 doi: 10.1002/anie.201407799

    23. [23]

      Chatelain L, Walsh J P S, Pecaut J, et al. Angew. Chem. Int. Ed., 2014,53:13434-13439 doi: 10.1002/anie.201407334

    24. [24]

      Zhang P, Zhang L, Wang C. J. Am. Chem. Soc., 2014,136: 4484-4489 doi: 10.1021/ja500793x

    25. [25]

      Ullman E F, Osiecki J H, Boocock D G B, et al. J. Am. Chem. Soc., 1972,94:7049-7059 doi: 10.1021/ja00775a031

    26. [26]

      Sheldrick G M. SHELXS-97:Program for the Solution of Crystal Structures, University of Göttingen, Germany, 1997.

    27. [27]

      Sheldrick G M. SHELXL-97:Program for the Refinement of Crystal Structures, University of Göttingen, Germany, 1997.

    28. [28]

      Wang C, Wang Y L, Qin Z X, et al. Inorg. Chem. Commun., 2012,20:112-117 doi: 10.1016/j.inoche.2012.02.030

    29. [29]

      Zhang C X, Qiao X M, Kong Y K, et al. J. Mol. Struct., 2015,108:1348-1354

    30. [30]

      Du F X, Hu P, Gao Y Y, et al. Inorg. Chem. Commun., 2014, 48:166-170 doi: 10.1016/j.inoche.2014.08.035

    31. [31]

      Zhang C X, Chen H W, Wang W M, et al. Inorg. Chem. Commun., 2012,24:177-181 doi: 10.1016/j.inoche.2012.07.003

    32. [32]

      Zhou N, Ma Y, Wang C, et al. Dalton Trans., 2009:8489-8492

    33. [33]

      Li L L, Liu S, Zhang Y, et al. Dalton Trans., 2015,44:6118-6125 doi: 10.1039/C4DT03941F

    34. [34]

      Sutter J-P, Golhen S, Kahn O, et al. Chem. Eur. J., 1998,4: 571-576 doi: 10.1002/ (ISSN) 1521-3765

  • Scheme 1  Molecular structure of NIT-Ph-4-Br

    图 1  Molecular structure (left) with thermal ellipsoids drawn at 30% probability and crystal packing diagram (right) of complex 1

    图 2  Temperature dependence of χMT (left) and field dependence of magnetization at 2.0 K (right) for complex 1

    Left: The solid lines represent the theoretical values based on the corresponding equations; Right: the solid line represents the theoretical curve for the sum of the Brillouin function for one Gd (III) (S=7/2) and two radical (S=1/2) with g=2

    图 3  Temperature dependence of χMT for complex 2 (left) and complex 3 (right)

    The solid lines represent the theoretical values based on the corresponding equations

    图 4  Field dependence of magnetization of 2 and 3 at 2.0 K

    图 5  Temperature dependence of the in-phase and out-of-phase components of ac susceptibility for 2 (left) and 3 (right) in zero dc field with an oscillation of 3.5 Oe

    表 1  Crystal data and structure refinement for 1, 2 and 3

    下载: 导出CSV

    表 2  Selected bond lengths (nm) and angles (°) for 1 , 2 and 3

    下载: 导出CSV
  • 加载中
计量
  • PDF下载量:  0
  • 文章访问数:  729
  • HTML全文浏览量:  179
文章相关
  • 收稿日期:  2016-08-13
  • 修回日期:  2016-11-09
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

/

返回文章