Reaction crystallization method based on deep eutectic solvents: A novel, green and efficient cocrystal synthesis approach

Xia-Lin Dai Yu-Hang Yao Jian-Feng Zhen Wei Gao Jia-Mei Chen Tong-Bu Lu

Citation:  Xia-Lin Dai, Yu-Hang Yao, Jian-Feng Zhen, Wei Gao, Jia-Mei Chen, Tong-Bu Lu. Reaction crystallization method based on deep eutectic solvents: A novel, green and efficient cocrystal synthesis approach[J]. Chinese Chemical Letters, 2025, 36(11): 110413. doi: 10.1016/j.cclet.2024.110413 shu

Reaction crystallization method based on deep eutectic solvents: A novel, green and efficient cocrystal synthesis approach

English

  • Cocrystal strategy is to construct multi-component crystals by assembling stoichiometric parent molecule and cocrystal co-former (CCF) in the same lattice through non-covalent interactions [1,2]. By manipulating the intermolecular interactions and packing patterns in the crystal, cocrystallization offers an opportunity to optimize the property and performance of functional molecules in an effective, designable and economical manner, which has been widely applied in various fields such as pharmaceutical, agrochemical, explosive and optoelectronic materials [3-8].

    Slurry method is a commonly used approach of cocrystal synthesis, involving suspending the parent molecule and CCF in a solvent and facilitating cocrystal formation via solvent-mediated phase transitions [9]. However, the cocrystallization efficiency is restricted by the isolation of less soluble component due to solubility difference of the two components of cocrystal, and the formation of unwanted solvates. To avoid individual component crystallization, a reaction crystallization method (RCM) was developed in the pharmaceutical field by using saturated solution of CCF to which an amount of active pharmaceutical ingredient (API) exceeding its solubility is added to eliminate the solubility difference between the API and CCF (Fig. 1a) [10,11]. Nevertheless, it also faces the challenge that even slight volatilization of organic solvent during the reaction process can induce substantial precipitation of soluble components, which is especially pronounced in large-scale cocrystal production [12]. Moreover, the issue of solvate formation remains unsolved with RCM [13]. Improper handling of organic solvents also poses risks to public health and environment [14]. Therefore, there is an urgent need to develop an innovate efficient and green cocrystallization approach to ensure the sustainable application of cocrystals.

    Figure 1

    Figure 1.  The cocrystal formation by (a) RCM, (b) DSM and (c) DRCM.

    Deep eutectic solvents (DESs) are room temperature liquids composed of hydrogen bond donors (HBD) and hydrogen bond acceptors (HBA) primarily connected by hydrogen bonds [15]. DESs are new generation of green and sustainable solvents due to unique advantages of biodegradability, reproducible, low cost, non or low toxicity, and low vapor pressure, making them a promising alternative to organic solvents and highly versatile for various industrial applications like organic synthesis, drug delivery, extraction and separation [16-21]. Recently, our group proposed a novel DESs based slurry method (DSM) for green synthesis of cocrystals (Fig. 1b), using DESs as reaction media and the HBD of DESs (such as carboxylic acids, amides and amino acids) as CCFs. The feasibility of DSM has been confirmed in a series of choline chloride DESs, which effectively avoid individual component crystallization and undesirable solvate formation, greatly improving cocrystal synthesis efficiency [22].

    To further verify the universal applicability of DSM, we attempted to use choline (ChOH) DESs as reaction media to synthesize some cocrystals of nicotinamide (NIC), carbamazepine (CBZ) and theophylline (THE). Fifteen DESs were prepared using ChOH as HBA, and a series of CCFs, which have been reported cocrystals with NIC, CBZ and THE, as HBD (Fig. 2) [23-40]. The DESs were characterized by fourier transform infrared spectroscopy (Fig. S1 in Supporting information) and differential scanning calorimetry (Fig. S2 in Supporting information) analyses, respectively, demonstrating the hydrogen bonds formation between HBA and HBD in DESs, and their liquid state at room temperature. Then, DSM experiments were carried out for a library of 33 cocrystal systems. However, NIC, CBZ and THE did not form any cocrystals but instead yielded pure APIs in all experiments (Table 1).

    Figure 2

    Figure 2.  The chemical structures of DESs materials and the model drugs.

    Table 1

    Table 1.  The cocrystal synthesis results of NIC, CBZ and THE using DSM and DRCM.
    DownLoad: CSV
    No. HBD (CCF) NIC: CCF DSM DRCM CBZ: CCF DSM DRCM THE: CCF DSM DRCM
    1 OA 1:1 × 1:1 × 2:1 ×
    2 MA 2:1 × 2:1 × 1:1 ×
    3 SA 2:1 × 2:1 × / / /
    4 GA 1:1 × 1:1 × 1:1 ×
    5 AA 1:1 × 2:1 × 2:1 ×
    6 PA 1:1 × / / / / / /
    7 SuA 1:1 × / / / / / /
    8 BA 1:1 × 1:1 × 1:1 ×
    9 4HBA 1:1 × 1:1 × 1:1 ×
    10 SaA 1:1 × 1:1 × 1:1 ×
    11 4AA / / / / / / 1:1 ×
    12 CA 1:1 × 1:1 × 1:1 ×
    13 MaA 1:1 × / / / / / /
    14 24DHBA 1:1 × / / / 1:1 ×
    15 26DHBA 1:1 × / / / / / /
    √ cocrystal formation, × no cocrystal formation, / no reported cocrystal.

    Previous study has shown that the cocrystal formation in DESs depends on the competitive hydrogen bonding interactions between the API, HBA and HBD molecules, and the dominant hydrogen bonds between the API and HBD could promote their cocrystal formation [22]. To explain why cocrystals of NIC, CBZ and THE cannot be obtained by DSM in ChOH DESs, molecular electrostatic potential surface (MEPS) analysis was conducted to analyze nucleophilic and electrophilic points [41,42], and predict non-covalent interactions between ChOH/API and HBD. ChOH exhibits much lower Epmin value (−63.42 kcal/mol) than CBZ (−53.22 kcal/mol), NIC (−42.52 kcal/mol), and THE (−39.83 kcal/mol), showing stronger competitiveness in forming non-covalent interactions with HBD (Fig. S3 in Supporting information). Artificial Bee Colony (ABC) algorithm was then employed to search for the low-energy stable molecular clusters in three representative THE-DESs systems, which can be approximated as the precursors of cocrystals [43]. As shown in Fig. S4 (Supporting information), no hydrogen bonded supra-molecular synthons between HBD and THE were observed in the first five stable clusters of each system. Combined the MEP analysis and ABC algorithm, it can be inferred that the presence of ChOH prevents mutual recognition of APIs and HBDs, making it difficult for the systems to generate ordered hetero supramolecular structures. This may be one reason for the experimental results of DSM.

    Inspired by conventional RCM in organic solvents, we envision that for compounds that cannot produce cocrystals by DSM, whether we can obtain cocrystals by stirring stoichiometric API and CCF using DESs which were saturated with the reactants in advance (Fig. 1c). It is essentially a DESs based reaction crystallization method (DRCM). The library of 33 cocrystal systems were also treated by DRCM. For each cocrystal system, a stoichiometric ratio mixture of API and HBD (as CCF) was added to 1 mL of DES. The obtained slurries were then stirred at 600 rpm under ambient conditions for 24 h. The solid phases were isolated and examined by powder X-ray diffraction (PXRD) and compared to the reported or simulated PXRD patterns. For those solids that the PXRD patterns are not consistent with the reported or simulated ones, the formation of the cocrystals was confirmed by single crystal structure determination. The outcome of all DRCM experiments is presented in Table 1. As we can see, all reported cocrystals that failed by DSM were successfully synthesized by DRCM, with the PXRD patterns well matched the simulated patterns (Fig. S5 in Supporting information). Additionally, a new cocrystal of THE with adipic acid (THE/AA) was discovered. The single crystal was obtained by volatilization method in acetonitrile and determined by single crystal X-ray diffraction analysis (Tables S1 and S2 in Supporting information), demonstrating it is a 2:1 cocrystal containing THE and adipic acid. Such results confirmed the feasibility of preparing cocrystal of the API and HBD by DRCM. The cocrystal formation process of DRCM is similar to traditional RCM, where supersaturated reactants in DESs increase the reaction activity and reduce the solubility of cocrystal, ultimately promoting the formation of cocrystals.

    To further examine the effect of slurry duration time and slurry speed on outcome and yield of cocrystals, THE cocrystal systems were selected as the research model and DRCM experiments in ten ChOH DESs were further performed at 600 rpm for 6, 12 h and at 300 and 900 rpm for 6 h. The results are compared with those of experiments performed at 600 rpm for 24 h and shown in Table 2 and Fig. S6 (Supporting information). When slurry at 600 rpm for 6 h, only OA, CA and 24DHBA completely converted to the respective cocrystals, while other systems exhibited no or incomplete conversion. By 12 and 24 h, all systems had achieved cocrystal formation. To verify whether different slurry speed could affect our DRCM experiments, we have also conducted 6 h slurry experiments of THE cocrystal systems at 300 and 900 rpm. The results reveal that the conversion to the cocrystals is complete at 900 rpm for all systems except for THE/GA system, while all cocrystal systems at 300 rpm exhibited no conversion except for THE/OA and THE/CA systems. The yields of all cocrystals are between 70%−85%, and extending the slurry duration time did not result in an increase in yield. In a word, slurry duration time affected the outcome of DRCM experiments and slurry speed affected the reaction rate.

    Table 2

    Table 2.  The cocrystal synthesis results of THE using DRCM with different slurry duration time and slurry speed.
    DownLoad: CSV
    No. Cocrystal Slurry duration time (h) and slurry speed (rpm)*
    6 h, 300 rpm 6 h, 600 rpm 6 h, 900 rpm 12 h, 600 rpm 24 h, 600 rpm
    1 THE/OA cocrystal cocrystal cocrystal cocrystal cocrystal
    78.6% 77.9% 80.8% 81.2% 79.3%
    2 THE/MA THE and CCF THE and CCF cocrystal cocrystal cocrystal
    82.3% 81.8% 83.1%
    3 THE/GA THE and CCF THE and CCF THE and CCF cocrystal cocrystal
    80.1% 77.8%
    4 THE/AA THE and CCF THE and CCF cocrystal cocrystal cocrystal
    78.3% 76.3% 73.2%
    5 THE/BA THE and CCF THE and CCF cocrystal cocrystal cocrystal
    79.4% 79.4% 80.5%
    6 THE/4HBA THE and CCF THE, CCF and cocrystal cocrystal cocrystal cocrystal
    79.4% 77.8% 74.2%
    7 THE/SaA THE and CCF THE and CCF cocrystal cocrystal cocrystal
    74.5% 70.5% 72.8%
    8 THE/4AA THE and CCF THE, CCF and cocrystal cocrystal cocrystal cocrystal
    85.2% 84.2% 83.0%
    9 THE/CA cocrystal cocrystal cocrystal cocrystal cocrystal
    83.7% 81.5% 85.5% 85.4% 84.1%
    10 THE/24DHBA THE and CCF cocrystal cocrystal cocrystal cocrystal
    81.4% 79.0% 78.9% 80.8%
    * Percentage values indicate the yields of cocrystals.

    Conventional RCM experiments for ten THE cocrystal systems using organic solvents, such as methanol, acetonitrile, ethyl acetate and acetone, were also conducted for 24 h to compare with DRCM. As shown in Table S3 (Supporting information), for the 40 RCM experiments in organic solvents saturated with THE and CCF, only 40% experiments generate high-purity cocrystals of THE, while 50% experiments produce cocrystal mixed with CCF due to organic solvent volatilization. Non-reaction and incomplete reaction were observed for THE/GA and THE/4AA cocrystal systems, which may be due to inappropriate solvent polarity or insufficient reaction time (Fig. S7 in Supporting information). By comparison, DRCM process produced pure cocrystal in all systems without soluble component precipitation attributed to the low volatility of DESs, indicating it is a more effective cocrystal synthesis method than traditional RCM in these systems.

    In this work, fifteen ChOH DESs were synthesized to further verify the universal applicability of DSM for cocrystal synthesis of NIC, CBZ and THE. However, due to the strong competitiveness of ChOH in forming non-covalent interactions with HBD, the desired API-HBD cocrystals cannot be obtained. In this situation, DRCM should be a viable alternative manner, which produce cocrystals of APIs and HBDs by increasing the reaction activity and reducing the solubility of cocrystals in DESs that were saturated with APIs and HBDs. DRCM not only inherits the advantages of DSM such as being green and safe, avoiding unwanted solvate formation and individual component crystallization events, but also avoids the defect of precipitation of soluble components during traditional RCM process caused by organic solvent volatilization. Undeniably, the DRCM also faces challenges such as high viscosity and solvent residue that need more in-depth and extensive research. But overall, the DRCM we proposed combines the advantages of DSM and RCM, providing a novel, green and economical approach for efficient and high-quality synthesis of cocrystals.

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    Xia-Lin Dai: Writing – original draft. Yu-Hang Yao: Data curation. Jian-Feng Zhen: Investigation. Wei Gao: Methodology. Jia-Mei Chen: Funding acquisition. Tong-Bu Lu: Visualization, Supervision.

    This work was financially supported by National Natural Science Foundation of China (Nos. 22101204 and 22271220). We also gratefully acknowledge the project grant from the Innovation and Strengthening Project of Guangdong Pharmaceutical University-Special Project of the Guangdong Education Commission (No. 2020KZDZX1128), the Research Projects of the Chinese Medicine Council of Guangdong Province (No. 20231209), and the Key Laboratory of Tropical Medicinal Resource Chemistry of the Ministry of Education at Hainan Normal University (No. RDZH2023001).

    Supplementary material associated with this article can be found, in the online version, at doi:10.1016/j.cclet.2024.110413.


    1. [1]

      J. Aaltonen, M. Allesø, S. Mirza, et al., Eur. J. Pharm. Biopharm. 71 (2009) 23–37.

    2. [2]

      M. Malamatari, S.A. Ross, D. Douroumis, et al., Adv. Drug Deliv. Rev. 117 (2017) 162–177.

    3. [3]

      M. Boterashvili, M. Lahav, S. Shankar, et al., J. Am. Chem. Soc. 136 (2014) 11926–11929. doi: 10.1021/ja5066587

    4. [4]

      H.M. Yu, B.X. Zhang, M.J. Liu, et al., Chin. Chem. Lett. 35 (2024) 109032.

    5. [5]

      Z.H. Li, J.M. Zhou, X. Zhang, et al., Chin. Chem. Lett. 34 (2023) 107983.

    6. [6]

      D. Yang, H. Wang, Q. Liu, et al., Chin. Chem. Lett. 33 (2022) 3207–3211.

    7. [7]

      X. Chen, D. Li, Z. Deng, et al., Cryst. Growth Des. 20 (2020) 6973–6982. doi: 10.1021/acs.cgd.0c01014

    8. [8]

      G. Bolla, B. Sarma, A.K. Nangia, Chem. Rev. 122 (2022) 11514–11603. doi: 10.1021/acs.chemrev.1c00987

    9. [9]

      A. Saha, A.A. Ahangar, A.A. Dar, et al., Cryst. Growth Des. 23 (2023) 7558–7581. doi: 10.1021/acs.cgd.2c01553

    10. [10]

      D. Herrmannsdörfer, T.M. Klapötke, Cryst. Growth Des. 21 (2021) 1708–1717. doi: 10.1021/acs.cgd.0c01611

    11. [11]

      N. Rodríguez-Hornedo, S.J. Nehm, K.F. Seefeldt, et al., Mol. Pharmaceutics 3 (2006) 362–367. doi: 10.1021/mp050099m

    12. [12]

      M.M. Haskins, M.J. Zaworotko, Cryst. Growth Des. 21 (2021) 4141–4150. doi: 10.1021/acs.cgd.1c00418

    13. [13]

      I.F.B. Biscaia, S.N. Gomes, L.S. Bernardi, et al., Pharmaceutics 13 (2021) 898. doi: 10.3390/pharmaceutics13060898

    14. [14]

      S. Wang, Y. Liu, S. Guo, et al., ACS Sustain. Chem. Eng. 10 (2022) 12411–12422. doi: 10.1021/acssuschemeng.2c03982

    15. [15]

      A.P. Abbott, G. Capper, D.L. Davies, et al., Chem. Commun. (2003) 70–71.

    16. [16]

      A.P. Abbott, D. Boothby, G. Capper, et al., J. Am. Chem. Soc. 126 (2004) 9142–9147.

    17. [17]

      Q.H. Zhang, K.D. Vigier, S. Royer, et al., Chem. Soc. Rev. 41 (2012) 7108–7146. doi: 10.1039/c2cs35178a

    18. [18]

      M. Ito, I. Yuriko, I. Hiroyasu, et al., Cryst. Growth Des. 23 (2023) 6298–6307. doi: 10.1021/acs.cgd.3c00155

    19. [19]

      A. Paiva, R. Craveiro, I. Aroso, et al., ACS Sustain. Chem. Eng. 2 (2014) 1063–1071. doi: 10.1021/sc500096j

    20. [20]

      Q. Zaib, M.J. Eckelman, Y. Yang, et al., Green Chem. 24 (2022) 7924–7930. doi: 10.1039/d2gc01752k

    21. [21]

      E.L. Smith, A.P. Abbott, K.S. Ryder, Chem. Rev. 114 (2014) 11060–11082. doi: 10.1021/cr300162p

    22. [22]

      J.F. Zhen, Y.H. Yao, W. Gao, et al., Cryst. Growth Des. 24 (2024) 5593–5602. doi: 10.1021/acs.cgd.4c00374

    23. [23]

      S.L. Childs, P.A. Wood, N. Rodriguez-Hornedo, et al., Cryst. Growth Des. 9 (2009) 1869–1888. doi: 10.1021/cg801056c

    24. [24]

      E. Lu, N. Rodríguez-Hornedo, R. Suryanarayanan, CrystEngComm 10 (2008) 665–668. doi: 10.1039/b801713c

    25. [25]

      J.A. McMahon, J.A. Bis, P. Vishweshwar, et al., Z. Kristallogr. 220 (2005) 340–350. doi: 10.1524/zkri.220.4.340.61624

    26. [26]

      V.R. Hathwar, R. Pal, T. Guru Row, Cryst. Growth Des. 10 (2010) 3306–3310. doi: 10.1021/cg100457r

    27. [27]

      A. Lemmerer, D.A. Adsmond, C. Esterhuysen, et al., Cryst. Growth Des. 13 (2013) 3935–3952. doi: 10.1021/cg4006357

    28. [28]

      L.J. Thompson, R.S. Voguri, A. Cowell, et al., Acta Crystallogr. Sect. C: Cryst. Struct. Commun. 66 (2010) o421.

    29. [29]

      S. Karki, T. Friščić, W. Jones, CrystEngComm 11 (2009) 470–481.

    30. [30]

      S. Lukin, T. Stolar, M. Tireli, et al., Chem. Eur. J. 23 (2017) 13941–13949. doi: 10.1002/chem.201702489

    31. [31]

      S.W. Zhang, I.A. Guzei, M.M. de Villiers, et al., Cryst. Growth Des. 12 (2012) 4090–4097. doi: 10.1021/cg3005757

    32. [32]

      J.I. Arenas-Garcia, D. Herrera-Ruiz, K. Mondragón-Vásquez, et al., Cryst. Growth Des. 10 (2010) 3732–3742. doi: 10.1021/cg1005693

    33. [33]

      D.A. Lorenzo, S.J.K. Forrest, H.A. Sparkes, Acta Crystallogr. Sect. B: Struct. Sci., Cryst. Eng. Mater. 72 (2016) 87–95.

    34. [34]

      D.J. Berry, C.C. Seaton, W. Clegg, et al., Cryst. Growth Des. 8 (2008) 1697–1712. doi: 10.1021/cg800035w

    35. [35]

      D.K. Bucar, R.F. Henry, G.G.Z. Zhang et al., Cryst. Growth Des. 14 (2014) 5318–5328. doi: 10.1021/cg501204k

    36. [36]

      S. Heiden, L. Trobs, K.J. Wenzel, et al., CrystEngComm 14 (2012) 5128–5129. doi: 10.1039/c2ce25236h

    37. [37]

      A.V. Trask, W.D. Motherwell, W. Jones, Int. J. Pharm. 320 (2006) 114–123.

    38. [38]

      B. Saikia, P. Bora, R. Khatioda, et al., Cryst. Growth Des. 15 (2015) 5593–5603. doi: 10.1021/acs.cgd.5b01293

    39. [39]

      S.L. Childs, G.P. Stahly, A. Park, Mol. Pharmaceutics 4 (2007) 323–338. doi: 10.1021/mp0601345

    40. [40]

      B. Sarma, B. Saikia, CrystEngComm 16 (2014) 4753–4765.

    41. [41]

      D. Musumeci, C.A. Hunter, R. Prohens, et al., Chem. Sci. 2 (2011) 883–890. doi: 10.1039/c0sc00555j

    42. [42]

      T. Grecu, C.A. Hunter, E.J. Gardiner, et al., Cryst. Growth Des. 14 (2014) 165–171. doi: 10.1021/cg401339v

    43. [43]

      D. Karaboga, B. Basturk, Appl. Soft Comput. 8 (2008) 687–697.

  • Figure 1  The cocrystal formation by (a) RCM, (b) DSM and (c) DRCM.

    Figure 2  The chemical structures of DESs materials and the model drugs.

    Table 1.  The cocrystal synthesis results of NIC, CBZ and THE using DSM and DRCM.

    No. HBD (CCF) NIC: CCF DSM DRCM CBZ: CCF DSM DRCM THE: CCF DSM DRCM
    1 OA 1:1 × 1:1 × 2:1 ×
    2 MA 2:1 × 2:1 × 1:1 ×
    3 SA 2:1 × 2:1 × / / /
    4 GA 1:1 × 1:1 × 1:1 ×
    5 AA 1:1 × 2:1 × 2:1 ×
    6 PA 1:1 × / / / / / /
    7 SuA 1:1 × / / / / / /
    8 BA 1:1 × 1:1 × 1:1 ×
    9 4HBA 1:1 × 1:1 × 1:1 ×
    10 SaA 1:1 × 1:1 × 1:1 ×
    11 4AA / / / / / / 1:1 ×
    12 CA 1:1 × 1:1 × 1:1 ×
    13 MaA 1:1 × / / / / / /
    14 24DHBA 1:1 × / / / 1:1 ×
    15 26DHBA 1:1 × / / / / / /
    √ cocrystal formation, × no cocrystal formation, / no reported cocrystal.
    下载: 导出CSV

    Table 2.  The cocrystal synthesis results of THE using DRCM with different slurry duration time and slurry speed.

    No. Cocrystal Slurry duration time (h) and slurry speed (rpm)*
    6 h, 300 rpm 6 h, 600 rpm 6 h, 900 rpm 12 h, 600 rpm 24 h, 600 rpm
    1 THE/OA cocrystal cocrystal cocrystal cocrystal cocrystal
    78.6% 77.9% 80.8% 81.2% 79.3%
    2 THE/MA THE and CCF THE and CCF cocrystal cocrystal cocrystal
    82.3% 81.8% 83.1%
    3 THE/GA THE and CCF THE and CCF THE and CCF cocrystal cocrystal
    80.1% 77.8%
    4 THE/AA THE and CCF THE and CCF cocrystal cocrystal cocrystal
    78.3% 76.3% 73.2%
    5 THE/BA THE and CCF THE and CCF cocrystal cocrystal cocrystal
    79.4% 79.4% 80.5%
    6 THE/4HBA THE and CCF THE, CCF and cocrystal cocrystal cocrystal cocrystal
    79.4% 77.8% 74.2%
    7 THE/SaA THE and CCF THE and CCF cocrystal cocrystal cocrystal
    74.5% 70.5% 72.8%
    8 THE/4AA THE and CCF THE, CCF and cocrystal cocrystal cocrystal cocrystal
    85.2% 84.2% 83.0%
    9 THE/CA cocrystal cocrystal cocrystal cocrystal cocrystal
    83.7% 81.5% 85.5% 85.4% 84.1%
    10 THE/24DHBA THE and CCF cocrystal cocrystal cocrystal cocrystal
    81.4% 79.0% 78.9% 80.8%
    * Percentage values indicate the yields of cocrystals.
    下载: 导出CSV
  • 加载中
计量
  • PDF下载量:  4
  • 文章访问数:  143
  • HTML全文浏览量:  9
文章相关
  • 发布日期:  2025-11-15
  • 收稿日期:  2024-07-24
  • 接受日期:  2024-09-05
  • 修回日期:  2024-08-23
  • 网络出版日期:  2024-09-06
通讯作者: 陈斌, bchen63@163.com
  • 1. 

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

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

/

返回文章