A metal-catalyzed triplet relay strategy for aza-Pauson-Khand reactions

Quoc Hoang Pham Ivan Sliusarevskyi Argha Saha Lennard Kloene Nico J. Linnartz Claire Empel Iris M. Oppel Debabrata Maiti Rene M. Koenigs

Citation:  Quoc Hoang Pham, Ivan Sliusarevskyi, Argha Saha, Lennard Kloene, Nico J. Linnartz, Claire Empel, Iris M. Oppel, Debabrata Maiti, Rene M. Koenigs. A metal-catalyzed triplet relay strategy for aza-Pauson-Khand reactions[J]. Chinese Chemical Letters, 2026, 37(10): 112306. doi: 10.1016/j.cclet.2025.112306 shu

A metal-catalyzed triplet relay strategy for aza-Pauson-Khand reactions

English

  • Catalytic amination reactions belong to the most important transformations in synthesis with widespread application in pharmaceutical, agrochemical and material research [19]. Not surprisingly, a vast range of methods has been developed in the past years, among which the transfer of a nitrene fragment, a highly reactive low-valent nitrogen species with one substituent, has received significant attention [911]. An important feature of low-valent intermediates lies within their spin state, which recently attracted the interest of chemists in the realm of photochemical carbene and nitrene transfer reactions [3,12,13]. These reactions proceed in the absence of a catalyst and via the respective free unbound species. The spin state of the reagent is controlled by steric, electronic or geometric properties of the reactive intermediate [13]. Recent reports describe significant advances, where conventional reactivity paradigms of such low-valent species can be overcome through the access of a triplet spin state. For instance, depending on the reaction conditions, the reactivity of triplet reagents with unsaturated systems can result in pathways, where conventional (2 + 1) cycloaddition is not taking place. Instead, coupling reactions or complex cycloaddition reactions can be accessed [14].

    The most common method for stabilization of low-valent nitrene intermediates, however, lies within the use of transition metal catalysts for stabilization of the pivotal nitrene intermediate and for control of nitrene transfer reactions. In the reaction with unsaturated compounds, commonly (2 + 1) cycloaddition reactions, or aziridination reactions, occur and catalytic higher-order cycloaddition reactions remain significantly underdeveloped (Fig. 1a).

    Figure 1

    Figure 1.  Cycloaddition reactions with nitrenes and their limitations. (A) (2 + 1) cycloaddition reactions. (B) Photochemical nitrene transfer reactions with olefines. (C) This work: (2 + 2 + 1) cycloaddition reactions via metal-catalyzed relay of triplet nitrene reactivity.

    Recent advances by Koenigs and co-workers demonstrated that free triplet nitrene intermediates can be generated under photochemical conditions enabling (2 + 1) cycloaddition reactions [11,15]. Moreover, such free triplet nitrene intermediates can be involved in more advanced cycloaddition reactions such as photochemical (2 + 2 + 1) cycloadditions (Fig. 1b) [14].

    While photochemical (2 + 2 + 1) cycloadditions rely on triplet nitrene generation, many substrates are unstable under irradiation or undergo competing photodecomposition, limiting their applicability. We hypothesized that identifying strategies to relay triplet nitrene reactivity under metal-catalyzed, light-free conditions would significantly expand the portfolio of advanced nitrene cycloaddition chemistry. A previous report describes two single compounds that arise from such advanced cycloaddition process albeit with severe synthetic limitations, such as stoichiometric amounts of metal salts and low yields [16]. We considered that a key challenge for this approach thus lies in the identification of suitable catalysts enabling a triplet relay strategy and avoiding the outcome of conventional aziridination reactions. To address this challenge, we considered two distinct approaches - either through the formation of an intermediate metal nitrene complex that adopts a doublet spin state possessing radical reactivity, or through formation of suitable metal triplet nitrene complex (Fig. 1c).

    To validate the viability of catalytic relay of the triplet nitrene reactivity, we examined the reaction of diallyl malonate 3a and iminoiodinane 4a. In this case, the traditional reactivity of metal-nitrene complexes results in the formation of an aziridine reaction product. Contrarily, the relay of the triplet reactivity would result in an aza-Pauson-Khand reaction, where the triplet reactivity is relayed by a 5-exo cyclization reaction, ultimately leading to the formation of a newly formed azabicyclic core system. For evaluation of this hypothesis, we therefore examined this transformation in the presence of a selection of archetypal nitrene transfer catalysts. Rhodium paddlewheel complexes, however, failed to give the desired (2 + 2 + 1) cycloaddition product, resulting mainly in the formation of tosyl amide. On the contrary, mixed-oxidation state Ru(Ⅱ,Ⅲ) paddlewheel complexes have been described to adopt a doublet spin state with nitrene reagents [1720]. In this case, we observed a good yield of the bicyclic pyrrolidine 5a, with little by-product formation. Notably, the reaction proceeded with high cis-diastereoselectivity and no formation of the trans-cyclization product was observed. In a similar fashion, a Ag(Ⅰ)Tp complex gave the desired bicyclic pyrrolidine 5a in good yield, presumably via triplet nitrene reactivity [2127]. Early transition metal complexes proved compatible, yet lower yields of the desired (2 + 2 + 1) cycloaddition reaction product were observed. Further reaction optimization consisted of the evaluation of further mixed-valency Ru(Ⅱ,Ⅲ) and AgTp catalysts, solvents, stoichiometries, etc. (for details, please see Tables S2-S5 in Supporting information). The most notable influence on the reaction yield was observed, when 3 equiv. of HFIP were added to the reaction mixture (Scheme 1). HFIP can effectively suppresses catalyst deactivation by engaging in stabilizing interactions with nucleophilic species that would otherwise poison the metal center [28]. It is further important to note that counterion exchange from chloride to BArF resulted in a significant improvement of catalyst reactivity, which may be rationalized by the tight binding of the chloride counterion, as also observed by single crystal structure analysis (Scheme 1).

    Scheme 1

    Scheme 1.  Realization of metal-catalyzed (2 + 2 + 1) cycloaddition reactions. Catalyst screening and reaction optimization. Reaction conditions: 1,6-diene (0.1 mmol), iminoiodinane (0.3 mmol), [M] ([M] = Ru or Rh 2.5 mol%/all other metal catalysts 5 mol%), CH2Cl2 (0.1 mol/L), r.t., 3 h; yields were determined by 1H NMR using mesitylene as an internal. a Reaction was carried out in CHCl3 (0.2 mol/L), r.t., 16 h. b Reaction was carried out in CH2Cl2 (0.025 mol/L), r.t., 16 h.

    The different reactivities of the BArF and chloride catalysts may further be explained by the differences in their respective crystal structures. While bond lengths and angles remain almost unaffected, exchanging the chloride ion for the weakly coordinating BArF leads to a significant increase in catalyst activity (Scheme 1). This can be explained by the tight binding of the chloride counterion to the ruthenium centers compared to BArF, which does not coordinate to the metal centers at all. It may also be rationalized from scXRD data that the free ligand sites around the ruthenium atoms can be occupied by solvent molecules, as exemplified by the coordination of two κO-DMF molecules and one methanol for the BArF and chloride structure, respectively (for details please see Supporting information). This can potentially elucidate as to why some solvents such as THF or trace of water appear to inhibit the reaction entirely, whereas others like MeCN do not (for details please see Table S2 in Supporting information). Solvents which can act as κO-donor ligands may block the active sites in a similar vein to chloride ions. In contrast, the catalyst appears less sensitive towards potential κN-coordinating solvents, which still afforded the product in a 50% yield when using MeCN as solvent (for details please see Table S2).

    However, it should be noted, that this observation relates to the bulk reaction since small amounts of HFIP actually enhance the reactivity and lead to more product formation (Scheme 1).

    A robustness screen further indicates that the reaction outcome is significantly altered, if the reaction is performed in air or in the presence of water. Less pronounced effects were observed with regards reaction temperature or concentration (for details, please see Figs. S1 and S2 in Supporting information).

    With the optimized conditions in hand for both catalysts, we then explored the generality of these catalytic (2 + 2 + 1)-cycloaddition protocols (Scheme 2). For the Ru(Ⅱ,Ⅲ)-catalyzed process, a wide array of iminoiodinanes bearing diversely substituted sulfonyl moieties were all found to be well-tolerated, all of which reacted smoothly to furnish the desired bicyclic pyrrolidine products in moderate to good yields (5a-l). The scope of the diallyl reaction partner was examined next. Various 1,6-heptadienes were first subjected to the established reaction conditions, which afforded the corresponding product in good yields (5m-x) except hepta-1,6-diene, which did not show any product formation in the ruthenium-catalyzed reaction. In comparison to the silver-catalyzed reaction, a significant decline in yield was observed when with different unsymmetrically substituted counterparts, although good to high diastereo-selectivity were observed in all cases (5m and 5r). Furthermore, the optimized protocol was found to tolerate various N-protected diallylamines as well as diallyl ether as substrates, all of which provided the corresponding products in moderate to good yields, even on a large scale (5a and 5u). Introduction of a methyl substituent to one of the olefin moieties—either at the terminal or internal position—resulted in only a slight decrease in yields, with modest diastereoselectivity observed in both cases (8 and 10).

    Scheme 2

    Scheme 2.  Investigations on the substrate scope and limitations of the present methodology. Reaction conditions: A: 1,6-diene (0.2 mmol), iminoiodinane (0.6 mmol), Ru2esp2BArF (2.5 mol%), HFIP (0.6 mmol), CHCl3 (0.2 mol/L), r.t., 16 h. Reaction conditions B: 1,6-diene (0.2 mmol), iminoiodinane (0.6 mmol), AgTpBr3 (2.5 mol%), HFIP (0.6 mmol), CH2Cl2 (0.025 mol/L), r.t., 16 h. Isolated yields are reported.

    Regarding the scope of the Ag(Ⅰ)Tp-catalyzed process, while the majority of iminoiodinanes employed were well tolerated, substrates bearing strong electron-withdrawing moieties proved to be completely inactive (for example compound 5e). No side products were detected in these cases, suggesting that strongly electron-withdrawing substituents hinder formation of the key Ag–nitrene complex required to initiate the C–N bond-forming step. Compared to the aforementioned Ru(Ⅱ,Ⅲ) protocol, overall, the optimized silver-catalyzed conditions proved to better tolerate different 1,6-diene derivatives, all of which afforded the desired product in moderate to good yields.

    However, a significant decline in yields was observed with various N-protected diallyl amines and no identifiable side products were detected (5t-w). Employment of core-substituted derivatives, similar to that previously observed with the Ru2-catalyzed procedure, had little effect on the efficacy of the transformation.

    The wide array of bicyclic bioisosteres prepared herein can serve as useful synthetic building blocks, particularly in a medicinal chemistry setting. Piperazine bioisoster 5u, bearing two orthogonal N-protecting groups, readily undergoes Boc-deprotection under mild acidic conditions to afford the free NH product 12 in 90% yield, which can be further derivatized to access valuable medicinally relevant scaffolds (Scheme 3A). Derivative 12 undergoes SNAr with 2,4-difluoronitrobenzene to afford 13, an analogue of several prospective mosquitocides, in 59% yield [29]. Antimalarial analogue 14 can also be accessed in 40% yield from 12 via simple tethering with 2-chloroethanol [30]. Furthermore, Clopidogrel analogue 15 can also be conveniently prepared via a N—H insertion reaction with the corresponding diazo precursor under photochemical conditions (39% yield). In addition, removal of tosyl group of 5a was also successfully realized with phenol in HBr/AcOH to afford the free NH analogue 16 in 39% yield (Scheme 3B).

    Scheme 3

    Scheme 3.  (A) Synthetic applications. Reaction conditions: (i) TFA (15 equiv.), CH2Cl2 (0.1 mol/L), r.t., overnight. (ii) Methyl 2-(2-chlorophenyl)-2-diazoacetate (1.1 equiv.), CH2Cl2 (0.1 mol/L), r.t., blue LED, 12 h. (iii) 2-Chloroethanol (1.5 equiv.), KI (1.5 equiv.), K2CO3 (1.2 equiv.), acetone (0.1 mol/L), reflux, 72 h. (iv) 2,4-Difluoro-1-nitrobenzene (1.8 equiv.), K2CO3 (2.8 equiv.), MeCN (0.2 mol/L), reflux, 18 h. (B) Removal of Ts-protecting group.

    Several control experiments were performed to probe the mechanisms of these catalytic processes (Scheme 4). For both protocols, the presence of a radical trapping agent such as TEMPO or DMPO, fully suppressed product formation, confirming the radical nature of both transformations. When investigating vinyl cyclopropanes as reaction parter, we only observed the decomposition of 4a and no product formation arising from a ring-opening reaction of the cyclopropane unit was detected (for details please see Supporting information). The involvement of an aziridine intermediate was then probed, in which 17 was subjected to the reaction conditions. In both cases, 1H NMR analysis of the crude reaction mixtures showed no trace of the desired bicyclic product 5a, and decomposition of the aziridine species was observed.

    Scheme 4

    Scheme 4.  Control experiments. (A) In the presence of radical trapping agents. (B) Investigation of 17. Detailed conditions information of A or B, see Scheme 2.

    To gain insight into the reaction mechanism of this transition metal-catalyzed (2 + 2 + 1)-cycloaddition, we conducted DFT studies. From a mechanism perspective, we consider that the relay of the triplet nitrene reactivity is achieved by through two related reaction mechanisms.

    For the silver-catalyzed variant, we consider the formation of a silver nitrene complex, as previously described by Pérez and co-workers [25]. This complex preferentially adopts a triplet spin state. Following initial C–N bond formation via the low-energy transition state (∆G = 2.9 kcal/mol), the reaction proceeds through cis-cyclization (∆Gcis = 3.9 kcal/mol vs. ∆Gtrans = 8.6 kcal/mol). After an intersystem crossing (ISC) via a low-lying minimum energy crossing point (MECP), the desired product is formed. Notably, the silver center remains coordinated to the nitrogen atom, while release of the free product is thermodynamically uphill by 8.2 kcal/mol. This observation provides a rationale for the importance of HFIP as an additive: by coordinating to the basic nitrogen, it facilitates a more facile dissociation of the silver catalyst and therefore enables a more efficient catalysis (Scheme 5A and Fig. S3 in Supporting information) [28].

    Scheme 5

    Scheme 5.  Computational studies. (A) Calculated reaction mechanism in the presence of AgTpBr3. (B) Formation of doublet spin state ruthenium nitrene complex.

    In the case of diruthenium(Ⅱ,Ⅲ) paddlewheel complexes, the respective metal nitrene complex takes up a doublet spin state as described by Du Bois and co-workers [17]. This stabilization is favored due to the mixed valency and results in a radical mechanism for C–N bond formation (Scheme 5B). Based on the radical nature of this nitrene complex, we consider that the observed (2 + 2 + 1) cycloaddition reaction therefore occurs through a radical mechanism.

    In summary, we herein report on transition metal-catalyzed (2 + 2 + 1) cycloaddition reactions employing 1,6-dienes and iminoiodinanes as reagents and HFIP as a key additive to improve reactivity of the metal catalyst. We identified diruthenium(Ⅱ,Ⅲ) paddlewheel complexes and silver scorpionate complexes as privileged catalysts to access this advanced cycloaddition reaction, while more commonly used dirhodium paddlewheel complexes proved incompatible. This concept was exemplified in the synthesis of a range of different bicyclic pyrrolidines and their application in drug analogue synthesis. This work describes a conceptually distinct strategy in accessing cycloaddition reactions with nitrenes beyond aziridination chemistry and it can be expected that future research will significantly expand the scope of complex cycloaddition reactions with nitrenes from both fundamental and applied perspective.

    Quoc Hoang Pham: Visualization, Validation, Methodology, Investigation. Ivan Sliusarevskyi: Validation, Methodology, Investigation. Argha Saha: Validation, Methodology, Investigation. Lennard Kloene: Validation, Methodology, Investigation. Nico J. Linnartz: Validation, Methodology, Investigation. Claire Empel: Writing – review & editing, Visualization, Validation, Resources, Formal analysis, Conceptualization. Iris M. Oppel: Validation, Supervision, Methodology, Investigation, Formal analysis. Debabrata Maiti: Writing – review & editing, Writing – original draft, Visualization, Validation, Supervision, Resources, Project administration, Formal analysis, Conceptualization. Rene M. Koenigs: Writing – review & editing, Writing – original draft, Visualization, Supervision, Resources, Project administration, Methodology, Funding acquisition, Formal analysis, Conceptualization.

    The authors declare no competing interests.

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


    1. [1]

      G. Dequirez, V. Pons, P. Dauban, Angew. Chem. Int. Ed. 51 (2012) 7384–7395. doi: 10.1002/anie.201201945

    2. [2]

      H. Hayashi, T. Uchida, Eur. J. Org. Chem. 2020 (2020) 909–916. doi: 10.1002/ejoc.201901562

    3. [3]

      C. Empel, R.M. Koenigs, Chem. Catal. 2 (2022) 2506–2514.

    4. [4]

      H.J. Dequina, C.L. Jones, J.M. Schomaker, Chem 9 (2023) 1658–1701.

    5. [5]

      D.N. Zalatan, J.Du Bois, Top. Curr. Chem. 292 (2009) 347–378. doi: 10.1007/128_2009_19

    6. [6]

      J.W.W. Chang, T.M.U. Ton, P.W.H. Chan, Chem. Rec. 11 (2011) 311–357. doi: 10.3938/NPSM.61.311

    7. [7]

      Y. Park, Y. Kim, S. Chang, Chem. Rev. 117 (2017) 9247–9301. doi: 10.1021/acs.chemrev.6b00644

    8. [8]

      K.M. Van Vliet, B. De Bruin, ACS Catal. 10 (2020) 4751–4769. doi: 10.1021/acscatal.0c00961

    9. [9]

      R. Singh, A. Mukherjee, ACS Catal. 9 (2019) 3604–3617. doi: 10.1021/acscatal.9b00009

    10. [10]

      B. Plietker, A. Röske, Catal. Sci. Technol. 9 (2019) 4188–4197. doi: 10.1039/c9cy00675c

    11. [11]

      Y. Guo, C. Pei, R.M. Koenigs, Nature Commun. 13 (2022) 86. doi: 10.58168/amsp2022_86-91

    12. [12]

      M. Bao, S. Dong, K. Chen, X. Xu, Chem. Asian J. 20 (2025) e202401521. doi: 10.1002/asia.202401521

    13. [13]

      C. Empel, Q.H. Pham, R.M. Koenigs, Acc. Chem. Res. 57 (2024) 2717–2727. doi: 10.1021/acs.accounts.4c00405

    14. [14]

      F. Li, W.F. Zhu, C. Empel, et al., Science 383 (2024) 498–503. doi: 10.1126/science.adm8095

    15. [15]

      Y. Guo, C. Pei, S. Jana, R.M. Koenigs, ACS Catal. 11 (2021) 337–342. doi: 10.1021/acscatal.0c04564

    16. [16]

      S. Minakata, D. Kano, R. Fukuoka, Y. Oderaotoshi, M. Komatsu, Heterocycles 60 (2003) 289–298. doi: 10.3987/COM-02-9648

    17. [17]

      M.E. Harvey, D.G. Musaev, J.Du Bois, J. Am. Chem. Soc. 133 (2011) 17207–17216. doi: 10.1021/ja203576p

    18. [18]

      T. Miyazawa, T. Suzuki, Y. Kumagai, et al., Nat. Catal. 3 (2020) 851–858. doi: 10.1038/s41929-020-00513-w

    19. [19]

      J.K. Sailer, J.C. Sharland, J. Bacsa, et al., Organometallics 42 (2023) 2122–2133. doi: 10.1021/acs.organomet.3c00268

    20. [20]

      K. Makino, Y. Kumagai, T. Yoshino, M. Kojima, S. Matsunaga, Org. Lett. 25 (2023) 3234–3238. doi: 10.1021/acs.orglett.3c00940

    21. [21]

      M.R. Rodríguez, M. Besora, F. Molina, et al., J. Am. Chem. Soc. 142 (2020) 13062–13071. doi: 10.1021/jacs.0c04395

    22. [22]

      S. Harada, M. Kobayashi, M. Kono, T. Nemoto, ACS Catal. 10 (2020) 13296–13304. doi: 10.1021/acscatal.0c04057

    23. [23]

      N.S. Dolan, R.J. Scamp, T. Yang, J.F. Berry, J.M. Schomaker, J. Am. Chem. Soc. 138, (2016) 14658–14667. doi: 10.1021/jacs.6b07981

    24. [24]

      E.G. Burke, J.M. Schomaker, Angew. Chem. Int. Ed. 54 (2015) 12097–12101. doi: 10.1002/anie.201504723

    25. [25]

      J. Llaveria, A. Beltrán, W.M.C. Sameera, et al., J. Am. Chem. Soc. 136 (2014) 5342–5350. doi: 10.1021/ja412547r

    26. [26]

      L. Maestre, W.M.C. Sameera, M.M. Díaz-Requejo, F. Maseras, P.J. Pérez, J. Am. Chem. Soc. 135 (2013) 1338–1348. doi: 10.1021/ja307229e

    27. [27]

      J.W. Rigoli, C.D. Weatherly, J.M. Alderson, B.T. Vo, J.M. Schomaker, J. Am. Chem. Soc. 135 (2013) 17238–17241. doi: 10.1021/ja406654y

    28. [28]

      J.C. Sharland, D. Dunstan, D. Majumdar, et al., ACS Catal. 12 (2022) 12530–12542. doi: 10.1021/acscatal.2c03909

    29. [29]

      C.D. Aretz, S.V. Kharade, K. Chronister, et al., ChemMedChem 16 (2021) 319–327. doi: 10.1002/cmdc.202000598

    30. [30]

      J.T. Kilbile, Y. Tamboli, S.S. Gadekar, et al., J. Mol. Struct. 1278 (2023) 134971.

  • Figure 1  Cycloaddition reactions with nitrenes and their limitations. (A) (2 + 1) cycloaddition reactions. (B) Photochemical nitrene transfer reactions with olefines. (C) This work: (2 + 2 + 1) cycloaddition reactions via metal-catalyzed relay of triplet nitrene reactivity.

    Scheme 1  Realization of metal-catalyzed (2 + 2 + 1) cycloaddition reactions. Catalyst screening and reaction optimization. Reaction conditions: 1,6-diene (0.1 mmol), iminoiodinane (0.3 mmol), [M] ([M] = Ru or Rh 2.5 mol%/all other metal catalysts 5 mol%), CH2Cl2 (0.1 mol/L), r.t., 3 h; yields were determined by 1H NMR using mesitylene as an internal. a Reaction was carried out in CHCl3 (0.2 mol/L), r.t., 16 h. b Reaction was carried out in CH2Cl2 (0.025 mol/L), r.t., 16 h.

    Scheme 2  Investigations on the substrate scope and limitations of the present methodology. Reaction conditions: A: 1,6-diene (0.2 mmol), iminoiodinane (0.6 mmol), Ru2esp2BArF (2.5 mol%), HFIP (0.6 mmol), CHCl3 (0.2 mol/L), r.t., 16 h. Reaction conditions B: 1,6-diene (0.2 mmol), iminoiodinane (0.6 mmol), AgTpBr3 (2.5 mol%), HFIP (0.6 mmol), CH2Cl2 (0.025 mol/L), r.t., 16 h. Isolated yields are reported.

    Scheme 3  (A) Synthetic applications. Reaction conditions: (i) TFA (15 equiv.), CH2Cl2 (0.1 mol/L), r.t., overnight. (ii) Methyl 2-(2-chlorophenyl)-2-diazoacetate (1.1 equiv.), CH2Cl2 (0.1 mol/L), r.t., blue LED, 12 h. (iii) 2-Chloroethanol (1.5 equiv.), KI (1.5 equiv.), K2CO3 (1.2 equiv.), acetone (0.1 mol/L), reflux, 72 h. (iv) 2,4-Difluoro-1-nitrobenzene (1.8 equiv.), K2CO3 (2.8 equiv.), MeCN (0.2 mol/L), reflux, 18 h. (B) Removal of Ts-protecting group.

    Scheme 4  Control experiments. (A) In the presence of radical trapping agents. (B) Investigation of 17. Detailed conditions information of A or B, see Scheme 2.

    Scheme 5  Computational studies. (A) Calculated reaction mechanism in the presence of AgTpBr3. (B) Formation of doublet spin state ruthenium nitrene complex.

  • 加载中
计量
  • PDF下载量:  0
  • 文章访问数:  8
  • HTML全文浏览量:  0
文章相关
  • 发布日期:  2026-10-15
  • 收稿日期:  2025-09-04
  • 接受日期:  2025-12-21
  • 修回日期:  2025-12-05
  • 网络出版日期:  2025-12-22
通讯作者: 陈斌, bchen63@163.com
  • 1. 

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

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

/

返回文章