留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

B/Al/Ga-MOR分子筛催化甲醇/二甲醚羰基化反应机理的理论计算研究

任鹏宇 刘卓 权燕红 郭军军 马宏 武建兵 王永钊

任鹏宇, 刘卓, 权燕红, 郭军军, 马宏, 武建兵, 王永钊. B/Al/Ga-MOR分子筛催化甲醇/二甲醚羰基化反应机理的理论计算研究[J]. 燃料化学学报(中英文). doi: 10.1016/S1872-5813(23)60395-0
引用本文: 任鹏宇, 刘卓, 权燕红, 郭军军, 马宏, 武建兵, 王永钊. B/Al/Ga-MOR分子筛催化甲醇/二甲醚羰基化反应机理的理论计算研究[J]. 燃料化学学报(中英文). doi: 10.1016/S1872-5813(23)60395-0
REN Pengyu, LIU Zhuo, QUAN Yanhong, GUO Jun-jun, MA Hong, WU Jianbing, WANG Yongzhao. Theoretical calculation study on the reaction mechanism of methanol/dimethyl ether carbonylation catalyzed by B/Al/Ga-MOR Zeolites[J]. Journal of Fuel Chemistry and Technology. doi: 10.1016/S1872-5813(23)60395-0
Citation: REN Pengyu, LIU Zhuo, QUAN Yanhong, GUO Jun-jun, MA Hong, WU Jianbing, WANG Yongzhao. Theoretical calculation study on the reaction mechanism of methanol/dimethyl ether carbonylation catalyzed by B/Al/Ga-MOR Zeolites[J]. Journal of Fuel Chemistry and Technology. doi: 10.1016/S1872-5813(23)60395-0

B/Al/Ga-MOR分子筛催化甲醇/二甲醚羰基化反应机理的理论计算研究

doi: 10.1016/S1872-5813(23)60395-0
基金项目: 国家自然科学基金(22072079, 22302115)太原理工大学煤科学与技术教育部重点实验室开放基金 (MKX202103)
详细信息
    通讯作者:

    Tel:15234060735;E-mail:wujianbing@sxu.edu.cn

    catalyst@sxu.edu.cn

  • 中图分类号: O64

Theoretical calculation study on the reaction mechanism of methanol/dimethyl ether carbonylation catalyzed by B/Al/Ga-MOR Zeolites

  • 摘要: 采用DFT理论计算详细比较了B、Al、Ga分别同晶取代MOR分子筛八元环侧袋T3位点及十二元环孔道T4位点时甲醇及二甲醚羰基化反应机制的共性及差异。研究发现,CO插入甲氧基生成乙酰基的反应遵循SN2机制,且为羰基化反应过程中的决速步;473K条件下,无论甲醇或二甲醚为原料,生成的乙酰基更倾向于与甲醇中的CH3O作用生成乙酸甲酯;对于羰基化反应和由三甲基氧鎓离子生成芳烃致催化剂失活的副反应,T3位点表现出更好的羰基化择形性而T4位点上更倾向于发生副反应。与Al-MOR相比,在T3位点引入B、Ga会导致羰基化反应能垒的升高,降低其催化性能;在T4位点引入B、Ga尤其是B则可大幅提升其生成三甲基氧鎓离子的能垒,进而抑制该过程进行,提升催化剂稳定性。本工作有助于认识MOR分子筛不同孔道内酸性位点发生同晶取代时催化羰基化反应机制的差异,为调控设计高效MOR沸石催化剂提供一定的理论支撑。
  • 图  1  MOR结构示意图(红球代表O原子,黄球代表Si原子)

    Figure  1  MOR-zeolite structure diagram (Red spheres represent O atoms, yellow spheres represent Si atoms)

    图  2  甲醇、二甲醚在MOR分子筛上的羰基化反应路径

    Figure  2  Reaction network for the carbonylation of MeOH and DME over MOR zeolite

    图  3  甲氧基与CO反应的反应物、过渡态和产物的结构:八元环侧袋 T3位点上insertion机理(a)和SN2机理(b)及十二元环孔道T4位点上insertion机理(c)和SN2机理(d)。其中红球代表O原子,黄球代表Si原子,紫球代表Al原子,白球代表H原子,灰球代表C原子;球棍模型代表反应中心原子,棍棒模型代表Brönsted酸位点,线条代表外层框架原子。

    Figure  3  Structures of ISs, TSs and FSs for the reaction of methoxy group with CO: the insertion mechanism (a) and SN2 mechanism (b) on the T3 site of the 8-MR side pocket, and the insertion mechanism (c) and SN2 mechanism (d) on the T4 site of the 12-MR channel. Red spheres represent O atoms, yellow spheres represent Si atoms, purple spheres represent Al atoms, white spheres represent O atoms, gray spheres represent C atoms; ball-and-stick parts represent the reaction center atoms, stick parts represent Brönsted acid sites, and the lines represent the outer frame atoms.

    图  4  473 K下B/Al/Ga-MOR分子筛八元环侧袋内甲醇至乙酸甲酯过程自由能曲线

    Figure  4  Free energy profiles of MeOH to methyl acetate in the 8-MR side pockets of B/Al/Ga-MOR at 473 K

    图  5  473 K下B/Al/Ga-MOR分子筛八元环侧袋内二甲醚至乙酸甲酯过程自由能曲线

    Figure  5  Free energy profiles of DME to methyl acetate in the 8-MR side pockets of B/Al/Ga-MOR at 473 K

    图  6  乙酰基物种与MeOH反应(TS3(M2))的反应物、过渡态和产物的结构:(a)、(b)、(c)依次为B/Al/Ga-MOR分子筛八元环侧袋 T3位点。其中红球代表O原子,黄球代表Si原子,粉球代表B原子,紫球代表Al原子,橙球代表Ga原子,白球代表H原子,灰球代表C原子;球棍模型代表反应中心原子,棍棒模型代表Brönsted酸位点,线条代表外层框架原子

    Figure  6  Structures of ISs, TSs and FSs for the reaction of acetyl species with MeOH (TS3(M2)): (a), (b) and (c) are located at the T3 sites of the 8-MR side pocket of B/Al/Ga-MOR, respectively. Red spheres represent O atoms, yellow spheres represent Si atoms, pink spheres represent B atoms, purple spheres represent Al atoms, orange spheres represent Ga atoms, white spheres represent O atoms, gray spheres represent C atoms; ball-and-stick parts represent the reaction center atoms, stick parts represent Brönsted acid sites, and the lines represent the outer frame atoms

    图  7  473 K下B/Al/Ga-MOR分子筛十二元环孔道内甲醇至乙酸甲酯过程自由能曲线

    Figure  7  Free energy profiles of MeOH to methyl acetate in the 12-MR channels of B/Al/Ga-MOR at 473 K

    图  8  473 K下B/Al/Ga-MOR分子筛十二元环孔道内二甲醚至乙酸甲酯过程自由能曲线

    Figure  8  Free energy profiles of DME to methyl acetate in the 12-MR channels of B/Al/Ga-MOR at 473 K

    表  1  473 K下甲氧基与 CO 在反应形成乙酰基的本征吉布斯自由能垒($ {\Delta \mathrm{G}}_{\mathrm{i}\mathrm{n}\mathrm{t}}^{\ne } $)

    Table  1  Computed intrinsic Gibbs free energy barriers ($ {\Delta \mathrm{G}}_{\mathrm{i}\mathrm{n}\mathrm{t}}^{\ne } $) for the reaction of methoxy group with CO to form an acetyl group at 473 K

    mechanism$ {\Delta \mathrm{G}}_{\mathrm{i}\mathrm{n}\mathrm{t}}^{\ne } $ (kJ mol−1)
    T3T4
    DFTDFT-DDFTDFT-D
    insertion220174220214
    SN2132130128139
    下载: 导出CSV

    表  2  473 K下B/Al/Ga-MOR分子筛八元环侧袋中甲醇及二甲醚羰基化反应的自由能垒($ {\Delta \mathrm{G}}_{\mathrm{i}\mathrm{n}\mathrm{t}}^{\ne } $)、速率常数($ k $)、焓势垒($ {\Delta H}_{{\rm{int}}}^{\ne } $)、熵损失($ {-T\Delta S}_{{\rm{int}}}^{\ne } $)和反应自由能($ {\Delta G}_{R} $)

    Table  2  Calculated kinetic results of free energy barriers ($ {\Delta \mathrm{G}}_{\mathrm{i}\mathrm{n}\mathrm{t}}^{\ne } $), relative rate constants ($ k $), enthalpy barriers ($ {\Delta H}_{{\rm{int}}}^{\ne } $) and entropy losses ($ {-T\Delta S}_{{\rm{int}}}^{\ne } $), and thermodynamic results of reaction free energies ($ {\Delta G}_{R} $) of each reaction step for MeOH and DME carbonylation in the 8-MR side pockets of B/Al/Ga-MOR at 473 K

    Reaction step
    8MR
    $ {\Delta \mathrm{G}}_{\mathrm{i}\mathrm{n}\mathrm{t}}^{\ne } $
    (kJ mol−1)
    $ k $
    (s−1)
    $ {\Delta H}_{{\rm{int}}}^{\ne } $
    (kJ mol−1)
    $ {-T\Delta S}_{{\rm{int}}}^{\ne } $
    (kJ mol−1)
    $ {\Delta G}_{R} $
    (kJ mol−1)
    B-MORTS1(M)1034.16 × 101976-50
    TS1(D)1151.97 × 100116−1−33
    TS21467.41 × 10−41361021
    TS3(M1)1711.29 × 10−615714−35
    TS3(M2)1332.02 × 10−211419−6
    TS3(D)3061.59 × 10−2127135−31
    TS42421.85 × 10−142202296
    Al-MORTS1(M)1133.27 × 100114−137
    TS1(D)872.43 × 103870−7
    TS21304.33 × 10−211317−1
    TS3(M1)971.91 × 1028611−80
    TS3(M2)656.54 × 1054124−14
    TS3(D)1952.87 × 10−91869−55
    TS42123.81 × 10−111872535
    Ga-MORTS1(M)1071.50 × 101108−133
    TS1(D)801.44 × 10481−1−16
    TS21377.31 × 10−3117205
    TS3(M1)1232.57 × 10−11230−69
    TS3(M2)836.72 × 1035033−13
    TS3(D)1872.20 × 10−815730−54
    TS42071.36 × 10−101832445
    下载: 导出CSV

    表  3  B/Al/Ga-MOR分子筛八元环侧袋和十二元环孔道内的质子亲和势和氨气吸附能

    Table  3  Proton affinities (PA, kJ mol−1) and NH3 adsorption energies ($ {\Delta E}_{ads-NH3} $, kJ mol−1) in the 8-MR side pockets and 12-MR channels of B/Al/Ga-MOR

    Acid sitesPA
    (kJ mol−1)
    $ {\Delta E}_{ads-NH3} $
    (kJ mol−1)
    T3B-MOR6.17−1.09
    Al-MOR5.89−1.74
    Ga-MOR6.12−1.59
    T4B-MOR6.28−1.09
    Al-MOR5.90−1.81
    Ga-MOR6.13−1.75
    下载: 导出CSV

    表  4  473 K下B/Al/Ga-MOR分子筛十二元环孔道内甲醇及二甲醚羰基化反应的自由能垒($ {\Delta \mathrm{G}}_{\mathrm{i}\mathrm{n}\mathrm{t}}^{\ne } $ )、速率常数($ k $)、焓势垒($ {\Delta H}_{{\rm{int}}}^{\ne } $)、熵损失($ {-T\Delta S}_{{\rm{int}}}^{\ne } $ )和反应自由能($ {\Delta G}_{R} $)

    Table  4  Calculated kinetic results of free energy barriers ($ {\Delta \mathrm{G}}_{\mathrm{i}\mathrm{n}\mathrm{t}}^{\ne } $), relative rate constants ($ k $), enthalpy barriers ($ {\Delta H}_{{\rm{int}}}^{\ne } $) and entropy losses ($ {-T\Delta S}_{{\rm{int}}}^{\ne } $), and thermodynamic results of reaction free energies ($ {\Delta G}_{R} $) of each reaction step for MeOH and DME carbonylation in the 12-MR channels of B/Al/Ga-MOR at 473 K

    Reaction step
    12MR
    $ {\Delta \mathrm{G}}_{\mathrm{i}\mathrm{n}\mathrm{t}}^{\ne } $
    (kJ mol−1)
    $ k $
    (s−1)
    $ {\Delta H}_{int}^{\ne } $
    (kJ mol−1)
    $ {-T\Delta S}_{int}^{\ne } $
    (kJ mol−1)
    $ {\Delta G}_{R} $
    (kJ mol−1)
    B-MORTS1(M)773.09 × 104752−53
    TS1(D)944.10 × 102868−35
    TS21881.70 × 10−814741−14
    TS3(M1)1287.21 × 10−211216−37
    TS3(M2)962.46 × 1028511−67
    TS3(D)1369.43 × 10−312412−45
    TS41341.57 × 10−21151986
    Al-MORTS1(M)935.29 × 10289422
    TS1(D)918.79 × 10295−410
    TS21394.40 × 10−310633−38
    TS3(M1)394.86 × 108309−53
    TS3(M2)351.34 × 1091025−93
    TS3(D)683.05 × 1054523−51
    TS4738.55 × 104694−34
    Ga-MORTS1(M)1261.20 × 10−1121551
    TS1(D)926.82 × 10293−16
    TS21394.40 × 10−310732−33
    TS3(M1)575.00 × 1062730−40
    TS3(M2)171.31 × 101189−94
    TS3(D)801.44 × 104755−46
    TS4845.21 × 103691537
    下载: 导出CSV
  • [1] LI Y, HE D W, NIU D J, et al. Acetic acid production from food wastes using yeast and acetic acid bacteria micro-aerobic fermentation[J]. Bioprocess Biosyst Eng,2015,38(5):863−869. doi: 10.1007/s00449-014-1329-8
    [2] SANTIAGO M A N, SÁNCHEZ-CASTILLO M A, CORTRIGHT R D, et al. Catalytic reduction of acetic acid, methyl acetate, and ethyl acetate over silica-supported copper[J]. J Catal,2000,193(1):16−28. doi: 10.1006/jcat.2000.2883
    [3] KALCK P, LE BERRE C and SERP P. Recent advances in the methanol carbonylation reaction into acetic acid[J]. Coord Chem Rev,2020,402(1):213078.
    [4] CHENG Z Z, HUANG S Y, LI Y, et al. Deactivation kinetics for the carbonylation of dimethyl ether to methyl acetate on H-MOR[J]. Ind Eng Chem Res,2017,56(46):13618−13627. doi: 10.1021/acs.iecr.7b03500
    [5] FAN B H, ZHANG W N, GAO P, et al. Quantitatively mapping the distribution of intrinsic acid sites in mordenite zeolite by high-field 23Na solid-state nuclear magnetic resonance[J]. J Phys Chem Lett,2022,13(23):5186−5194. doi: 10.1021/acs.jpclett.2c00932
    [6] LIU S P, FANG X D, LIU Y, et al. Dimethyl ether Carbonylation over mordenite zeolite modified by Alkyimidazolium ions[J]. Catal Commun,2020,147:106161. doi: 10.1016/j.catcom.2020.106161
    [7] HAM H, JUNG H S, KIM H S, et al. Gas-phase carbonylation of dimethyl ether on the stable seed-derived ferrierite[J]. ACS Catal,2020,10(9):5135−5146. doi: 10.1021/acscatal.9b05144
    [8] FENG X B, YAO J, LI H J, et al. A brand new zeolite catalyst for carbonylation reaction[J]. Chem Commun,2019,55(8):1048−1051. doi: 10.1039/C8CC08411D
    [9] FUJIMOTO K, SHIKADA T, OMATA K, et al. Vapor phase carbonylation of methanol with solid acid catalysts[J]. Chem Lett,1984,13(12):2047−2050. doi: 10.1246/cl.1984.2047
    [10] CHEUNG P, BHAN A, SUNLEY G J, et al. Selective carbonylation of dimethyl ether to methyl acetate catalyzed by acidic zeolites[J]. Angew Chem,2006,118(10):1647−1650. doi: 10.1002/ange.200503898
    [11] CHEUNG P, BHAN A, SUNLEY G J, et al. Site requirements and elementary steps in dimethyl ether carbonylation catalyzed by acidic zeolites[J]. J Catal,2007,245(1):110−123. doi: 10.1016/j.jcat.2006.09.020
    [12] NARSIMHAN K, MICHAELIS V K, MATHIES G, et al. Methane to acetic acid over Cu-exchanged zeolites: Mechanistic insights from a site-specific carbonylation reaction[J]. J Am Chem Soc,2015,137(5):1825−1832. doi: 10.1021/ja5106927
    [13] BLASCO T, BORONAT M, CONCEPCIÓN P, et al. Carbonylation of methanol on metal-acid zeolites: Evidence for a mechanism involving a multisite active center[J]. Angew Chem,2007,119(21):4012−4015. doi: 10.1002/ange.200700029
    [14] BORONAT M, MARTINEZ-SANCHEZ C, LAW D, et al. Enzyme-like specificity in zeolites: A unique site position in mordenite for selective carbonylation of methanol and dimethyl ether with CO[J]. J Am Chem Soc,2008,130(48):16316−16323. doi: 10.1021/ja805607m
    [15] LI B J, XU J, HAN B, et al. Insight into dimethyl ether carbonylation reaction over mordenite zeolite from in-situ solid-state NMR spectroscopy[J]. J Phys Chem C,2013,117(11):5840−5847. doi: 10.1021/jp400331m
    [16] CHU Y Y, LO A Y, WANG C, et al. Origin of high selectivity of dimethyl ether carbonylation in the 8-membered ring channel of mordenite zeolite[J]. J Phys Chem C,2019,123(25):15503−15512. doi: 10.1021/acs.jpcc.9b01874
    [17] WANG S R, GUO W W, ZHU L J, et al. Methyl acetate synthesis from dimethyl ether carbonylation over mordenite modified by cation exchange[J]. J Phys Chem C,2014,119(1):524−533.
    [18] REULE ALLEN A C, SEMAGINA N. Zinc hinders deactivation of copper-mordenite: Dimethyl ether carbonylation[J]. ACS Catal,2016,6(8):4972−4975. doi: 10.1021/acscatal.6b01464
    [19] ALBERTI A. Location of Brønsted sites in mordenite[J]. Zeolites,1997,19(5-6):411−415. doi: 10.1016/S0144-2449(97)00114-0
    [20] LI Y, SUN Q, HUANG S Y, et al. Dimethyl ether carbonylation over pyridine-modified MOR: Enhanced stability influenced by acidity[J]. Cataly Today,2018,311(1):81−88.
    [21] NI Y M, SHI L, LIU H C, et al. A green route for methanol carbonylation[J]. Catal Sci Technol,2017,7(20):4818−4822. doi: 10.1039/C7CY01621B
    [22] HE T, LIU X C, XU S T, et al. Role of 12-ring channels of mordenite in DME carbonylation investigated by solid-state NMR[J]. J Phys Chem C,2016,120(39):22526−22531. doi: 10.1021/acs.jpcc.6b07958
    [23] LU P, YANG G H, TANAKA Y, et al. Ethanol direct synthesis from dimethyl ether and syngas on the combination of noble metal impregnated zeolite with Cu/ZnO catalyst[J]. Cataly Today,2014,232(1):22−26.
    [24] MA M, ZHAN E S, HUANG X M, et al. Carbonylation of dimethyl ether over Co-HMOR[J]. Catal Sci Technol,2018,8(8):2124−2130. doi: 10.1039/C8CY00407B
    [25] LI Y, HUANG S Y, CHENG Z Z, et al. Promoting the activity of Ce-incorporated MOR in dimethyl ether carbonylation through tailoring the distribution of Brønsted acids[J]. Appl Catal B,2019,256(5):117777.
    [26] KRESSE G, HAFNER J. Ab initio molecular-dynamics simulation of the liquid-metal-amorphous-semiconductor transition in germanium[J]. Phys Rev B:Condens Matter Mater Phys,1994,49(20):14251−14269. doi: 10.1103/PhysRevB.49.14251
    [27] GRIMME S. Semiempirical GGA-type density functional constructed with a long-range dispersion correction[J]. J Comput Chem,2006,27(15):1787−1799. doi: 10.1002/jcc.20495
    [28] WANG V, XU N, LIU J C, et al. VASPKIT: A user-friendly interface facilitating high-throughput computing and analysis using VASP code[J]. Comput Phys Commun,2021,267:108033. doi: 10.1016/j.cpc.2021.108033
    [29] EYRING H. The Activated complex in chemical reactions[J]. J Chem Phys,1935,3:107−115. doi: 10.1063/1.1749604
    [30] YUAN S P, WANG J G, LI Y M, et al. Theoretical studies on the properties of acid site in isomorphously substituted ZSM-5[J]. Mole. Catal. A:Chemical,2002,178(1-2):267. doi: 10.1016/S1381-1169(01)00335-1
    [31] JONES A J, IGLESIA E. The strength of Brønsted acid sites in microporous Aluminosilicates[J]. ACS Catal,2015,5:5741−5755. doi: 10.1021/acscatal.5b01133
    [32] WANG S, LI S Y, ZHANG L, et al. Mechanistic insights into the catalytic role of various acid sites on ZSM-5 zeolite in the carbonylation of methanol and dimethyl ether[J]. Catal Sci Technol,2018,8(12):3193−3204. doi: 10.1039/C8CY00296G
    [33] BORONAT M, MARTINEZ C, CORMA A. Mechanistic differences between methanol and dimethyl ether carbonylation in side pockets and large channels of mordenite[J]. Phys Chem Chem Phys,2011,13(7):2603−2612. doi: 10.1039/c0cp01996h
    [34] RASMUSSEN D B, CHRISTENSEN J M, TEMEL B, et al. Ketene as a reaction intermediate in the carbonylation of dimethyl ether to methyl acetate over mordenite[J]. Angew Chem Int Ed Engl,2015,54(25):7261−7264. doi: 10.1002/anie.201410974
    [35] RASMUSSEN D B, CHRISTENSEN J M, TEMEL B, et al. Reaction mechanism of dimethyl ether carbonylation to methyl acetate over mordenite-a combined DFT/experimental study[J]. Catal Sci Technol,2017,7(5):1141−1152. doi: 10.1039/C6CY01904H
    [36] LIU Z Q, YI X F, WANG G R, et al. Roles of 8-ring and 12-ring channels in mordenite for carbonylation reaction: From the perspective of molecular adsorption and diffusion[J]. J Catal,2019,369:335−344. doi: 10.1016/j.jcat.2018.11.024
    [37] KING S T. Reaction mechanism of oxidative carbonylation of methanol to dimethyl carbonate in Cu–Y zeolite[J]. J Catal,1996,161(2):530−538. doi: 10.1006/jcat.1996.0215
    [38] CALDERAZZO F, COTTON F A. Carbon monoxide insertion reactions. I. the carbonylation of methyl manganese pentacarbonyl and decarbonylation of acetyl manganese pentacarbonyl[J]. Inorg Chem,1962,1(1):30−36. doi: 10.1021/ic50001a008
    [39] MA H, LIAO J, WEI Z H, et al. Trimethyloxonium ion - a zeolite confined mobile and efficient methyl carrier at low temperatures- a DFT study coupled with microkinetic analysis[J]. Catal Sci Technol,2022,12:3328. doi: 10.1039/D2CY00207H
    [40] ZHOU Z Q, LIU H C, NI Y M, et al. Direct conversion of dimethyl ether and CO to acetone via coupling carbonylation and ketonization[J]. J Catal,2021,396:360−373. doi: 10.1016/j.jcat.2021.03.006
    [41] Bhan A, Allian A D, Sunley G J, et al. Specificity of sites within eight-membered ring zeolite channels for carbonylation of methyls to acetyls[J]. J Am Chem Soc,2007,129(16):4919−4924. doi: 10.1021/ja070094d
    [42] YUAN S P, WANG J G, LI Y W, et al. Density functional investigations into the siting of Fe and the acidic properties of isomorphously substituted mordenite by B, Al, Ga and Fe[J]. J Mol Struct:THEOCHEM,2004,674(1-3):267−274. doi: 10.1016/S0166-1280(03)00463-9
    [43] CAI K, HUANG S Y, LI Y, et al. Influence of acid strength on the reactivity of dimethyl ether carbonylation over H-MOR[J]. ACS Sustainable Chem Eng,2018,7(2):2027−2034.
  • 加载中
图(8) / 表(4)
计量
  • 文章访问数:  21
  • HTML全文浏览量:  5
  • PDF下载量:  3
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-10-08
  • 修回日期:  2023-10-24
  • 录用日期:  2023-10-24
  • 网络出版日期:  2023-11-10

目录

    /

    返回文章
    返回