留言板

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

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

基于EBT3胶片开展电子FLASH放疗剂量测量研究

王诗岚 羊奕伟 程德琪 唐镭迅 吴岱

孙列鹏, 袁震宇, 张诚, 等. 强流RFQ的固态功率源模块故障分析[J]. 强激光与粒子束, 2019, 31: 065103. doi: 10.11884/HPLPB201931.180245
引用本文: 王诗岚, 羊奕伟, 程德琪, 等. 基于EBT3胶片开展电子FLASH放疗剂量测量研究[J]. 强激光与粒子束, 2024, 36: 126005. doi: 10.11884/HPLPB202436.240095
Sun Liepeng, Yuan Zhenyu, Zhang Cheng, et al. Analysis of modules failure in solid-state amplifier for high current RFQ[J]. High Power Laser and Particle Beams, 2019, 31: 065103. doi: 10.11884/HPLPB201931.180245
Citation: Wang Shilan, Yang Yiwei, Cheng Deqi, et al. Study of electron FLASH radiotherapy dose measurement based on EBT3 film[J]. High Power Laser and Particle Beams, 2024, 36: 126005. doi: 10.11884/HPLPB202436.240095

基于EBT3胶片开展电子FLASH放疗剂量测量研究

doi: 10.11884/HPLPB202436.240095
基金项目: 国家自然科学基金项目(12375318)
详细信息
    作者简介:

    王诗岚,wangshilan2024@163.com

    通讯作者:

    羊奕伟,winfield1920@126.com

    吴 岱,wudai04@163.com

  • 中图分类号: R144.1

Study of electron FLASH radiotherapy dose measurement based on EBT3 film

  • 摘要: 闪光放疗使用超高剂量率在毫秒时间内将剂量全部注入靶区,其超高剂量率使现有的在线剂量计基本失效,目前通常使用辐射显色胶片来测量剂量。基于中国工程物理研究院应用电子学研究所研制的电子加速器搭建了电子FLASH放疗平台,基于EBT3胶片的快速读出方法,研究了此平台的剂量率范围及剂量分布。实验结果表明,EBT3胶片的快速读出方法可用于电子FLASH放疗的剂量测量,在源皮距100 cm及深度1 cm处剂量率在240~290 Gy/s之间;电子束到达模体表面的平均能量的波动会导致靶区约±5%的剂量波动;面剂量分布满足平坦性在±5%以内和对称性在±3%以内的要求。
  • 固体材料按照其导电特性可以分为导体、半导体和绝缘体三大类,既可以通过电子的填充能带理解,也可以通过材料的电阻率特性来区分。从第一代的硅、锗材料发展到第二代的砷化镓和磷化铟等材料,以及目前市场需求巨大的第三代宽禁带半导体碳化硅、氮化镓、金刚石等材料,说明半导体材料向宽禁带的发展越来越快,应用场合也从晶体管、集成电路等,拓展到通讯、激光、电力电子以及微波等领域。在高温环境下由于禁带宽度大,其本征激发所需要的温度越高,这就大大提高了材料的工作温度,以硅和碳化硅作为比较,硅材料的极限温度为300℃,而碳化硅可以达到600 ℃以上[1, 2]。在其他特性方面比如其热学特性,碳化硅的热导率为4.9 W·cm-1·K-1,而硅材料仅有1.5 W·cm-1·K-1,这就使得碳化硅等高热导率材料在一定工作温度范围内,无需增加散热装置,减小器件体积[3]。目前对于宽禁带半导体材料主要为碳化硅、氮化镓、金刚石、氧化镓、氮化铝等材料,其禁带宽度都大于3 eV。由于其不同的材料特性,因此其应用环境有所侧重。碳化硅作为第三代宽禁带半导体材料中具有代表性的一员,正是在禁带宽度、载流子饱和迁移率、热导率、临界击穿场强、抗辐射能力等性质的优异表现,越来越受到产业界的广泛关注。目前利用碳化硅材料制备的SBD和MOSFET已经实现商品化,并广泛应用于航空航天、军事装备、新能源汽车、超高压输变电、智能电网、高速轨道交通等领域。因此近年来全球在碳化硅材料生长和器件制备投入的研发费用和生产规模迅速增长,对碳化硅衬底提出了迫切需要,大尺寸高质量碳化硅单晶生长和碳化硅器件制备一直是国际上的研究热点。本文介绍碳化硅单晶衬底的发展以及利用碳化硅材料制备的低导通光导开关,分析单晶生长技术以及器件性能等,并对碳化硅材料的发展趋势进行展望。

    碳化硅材料具有200多种同质多形体[4-5],其化学计量成分相同,但是结构不同,这是由碳化硅中Si-C双原子层的堆叠方式不同所决定的。碳化硅晶型结构可以分成六方、立方、菱方结构,但是所有的结构都可以用六方系统来描述[6]。碳化硅材料的晶格常数a可以看作常数,而晶格常数c可以有不同的数值,通过晶型结构和晶格常数的影响使得碳化硅材料可以拥有多种同质多形体。这些同质多形体形成的自由能非常相近,但是它们之间具有很高的能量势垒,因此在基本成分相同的情况下,多种同质多形体的物理性质是独立存在的。

    多种的碳化硅晶型来说,3C-SiC,6H-SiC,4H-SiC以及15R-SiC最为常见,图 1给出了不同构型的碳化硅单晶的硅碳双原子层占据三种位置的示意图[7]。密堆积有三种不同的位置,记为A,B,C,如果第一层为A位置,根据密排结构原则,第二层将占据B或者C位置,如果第二层占据B位置,那么第三层将占据A或者C位置;如果第二层占据C位置,那么第三层占据A或者B位置,以此类推,因此由于不同的堆叠方式,使得碳化硅材料产生了多种多样的晶型。

    图  1  不同构型的SiC单晶的硅碳双原子层排列顺序
    Figure  1.  Stacking sequence of double atomic layers of different SiC polytypes

    SiC晶体结构决定了其独特的物理、化学性质(表 1)[3]。相对于Si和GaAs,SiC拥有更为优越的物理性质;禁带宽度大,接近于Si的3倍,保证器件在高温工作下的长期可靠;击穿场强高,是Si的10倍,保证器件耐压容量,提高器件的耐压值;饱和电子速率大,是Si的2倍,增加器件的工作频率和功率密度;热导率高,超过Si的3倍,增加了器件的散热能力,实现设备的小型化等。

    表  1  SiC与Si和GaAs的物理特性参数比较
    Table  1.  Properties comparison between Si, GaAs and SiC
    material band gap/eV dielectric constant breakdown field/(MV·cm-1) saturated electron drift velocity/(cm·s-1) intrinsic carrier concentration/cm-3 electron mobility/(cm2·V·s-1) thermal conductivity/(W·cm-1·K-1)
    Si 1.12 11.8 0.3 1.0×107 1.5×1010 1400 1.50
    GaAs 1.43 12.8 0.6 1.0×107 1.8×106 8500 0.46
    6H-SiC 3.03 9.6 3.2 2.0×107 2.3×10-6 400 4.90
    4H-SiC 3.26 9.7 3.0 2.0×107 8.2×10-9 1140 4.90
    下载: 导出CSV 
    | 显示表格

    目前碳化硅单晶的生长主要集中到以下两种方法:液相生长法和物理气相传输法(PVT)[8]。液相生长法主要集中在日本的高校和科研院所,其生长如图 2(a)所示。采用中频加热,高纯石墨坩埚作为容器,同时提供碳源。溶液加热到1500~1900 ℃保温数小时,黏在石墨棒上的籽晶跟随着石墨棒一同浸入溶液中,由于石墨坩埚中的温差,提供了晶体生长的过冷度,进而在籽晶上生长晶体[8-9]。物理气相传输法(如图 2(b)所示)生长碳化硅单晶是目前生长大直径、高质量碳化硅单晶最常用的方法[6]。该方法使用感应线圈进行加热。由于法拉第电磁感应定律,坩埚外壁会产生感应电动势,感应电动势产生后由于石墨的导电性会在石墨坩埚中形成涡流,涡流产生焦耳热使石墨坩埚温度升高。通过热辐射、热传导和对流等方式,将热量传递给坩埚内的碳化硅粉料和籽晶,最终建立起单晶生长温度场。

    图  2  碳化硅单晶生长示意图
    Figure  2.  Schematic of SiC growth geometry models

    随着半导体材料应用到商用的功率器件,衬底直径的大小成为影响器件成本的重要因素,随着碳化硅直径的增加,其功率器件的成本将大大降低。由于碳化硅材料的结构以及生长环境的限制,早期在大尺寸碳化硅单晶生长方面相对困难。1978年,改进Lely法的提出,许多国家开始研究碳化硅材料和器件。1991年美国Cree公司采用升华法生长出商品化2英寸6H-SiC衬底。2013年,Cree公司推出商品化6英寸的碳化硅衬底,逐渐代替4英寸衬底成为市场上的主流产品。2015年,Cree公司在意大利召开的SiC国际会议上“ICSCRM2015”,展示了其研究级8英寸口径衬底试样(如图 3所示)。目前国外从事SiC单晶生长研究的企业主要包括美国的Cree公司,Dow Corning公司、II-VI公司;德国的SiCrystal AG公司;瑞典的Okmetic公司;日本的Nippon Steel公司。

    图  3  Cree公司SiC产品:8 inch试样
    Figure  3.  SiC substrates made by Cree: 8inch wafer

    国内也相继开展了碳化硅单晶生长研究,主要包括山东大学、中国科学院物理研究所、中国科学院硅酸盐研究所、中国电子科技集团公司第46研究所等单位。以相关的技术为基础,能批量生产SiC单晶衬底的公司包括:山东天岳先进材料科技有限公司、北京天科合达半导体股份有限公司等。到目前为止,国内高质量的4英寸碳化硅单晶衬底已经实现商品化,6英寸仍处于研发和试生产阶段。山东大学于2013年10月成功生长出6英寸SiC单晶(图 6),于2015年6月生长的掺V半绝缘6英寸SiC单晶[10],整片电阻率均大于5×109 Ω·cm,X射线摇摆曲线半峰宽小于30″(弧秒),晶体结晶质量好(如图 4所示)。

    图  4  山东大学6 inch SiC单晶及衬底
    Figure  4.  6 inch SiC single crystal and substrate grown by Shandong University

    光导开关是一种通过光控制载流子的产生与复合来实现开关的导通和关闭状态的新型超快半导体器件,其工作时域可达到ps量级,是产生高功率超短脉冲的关键器件。作为脉冲功率系统中的新型开关器件,相对于传统的脉冲功率器件,光导开关具有重复频率性能好、闭合时间短(ps量级)、时间抖动小(ps量级)、开关电感低(亚nH)、同步精度高(ps量级)、电磁兼容性强等优点,使其在固态紧凑型脉冲功率源上有着较为广阔的应用前景[11]。1995年,P. S. Cho等人研制了横向结构的6H-SiC光导开关[12]。1997年,Senpeng Sheng等人采用CVD方法在石墨烯衬底上沉积3C-SiC薄膜,并在薄膜上制作了光导开关,获得250 kV的击穿电压和10 kA/cm2的峰值电流,电流脉宽为15 ns[13]。2003年,S.Dogan等人制作了4H-SiC光导开关[14]。2005年,K. Zhu在4H-SiC上外延高掺杂的n-GaN,使器件的欧姆接触大大提高,增加了光导开关的寿命[15]。2015年,Mauch等人分析了在瞬时高电场强度下,光导开关的失效机理[16]。同年,Tiskumara等人分析了半绝缘4H-SiC光导开关的稳态下的电流-电压响应,通过仿真模拟研究了基于电场驱动效应的内部机制[17]。2016年,Chowdhury等人通过仿真模拟增加高介质层材料氧化铪,降低电极边缘场强高的问题[18]。国内2012年,中国工程物理研究院流体物理研究所袁建强等人,制作了导通电阻为22 Ω的体结构PCSS[19]。2014年,上海硅酸盐研究所刘学超等人报道的10 mm间隙的6H-SiC光导开关,工作电压14 kV时,电流为185 A,导通电阻为22 Ω[20]。2018年,西安电子科技大学利用Al掺杂的ZnO作为窗口,并利用银镜反射,制备了异面电极4H-SiC光导开关,耐压值为12 kV,导通电阻为7.5 Ω[21]。2018年,山东大学通过制作同面电极高纯半绝缘4H-SiC光导开关,在触发波长355 nm,触发能量10 mJ,偏置电压6 kV的情况下,得到导通电阻小于1 Ω[22]。同年,中国工程物理研究院流体物理研究所与山东大学合作,制备了4H-SiC光导开关,并发现一种微波振荡现象,振荡频率为1 GHz,这是目前基于碳化硅光导开关首次发现由于材料本身特性引起的微波振荡[23]

    根据器件特性以及应用环境的不同,需要选择不同的补偿机制制备半绝缘碳化硅衬底材料。碳化硅光导开关一般工作在线性模式,光生载流子浓度对于器件性能起到至关重要的作用。对于高纯半绝缘4H-SiC光导开关主要利用本征吸收方式,利用波长为355 nm,脉冲宽度为17 ns的脉冲激光器进行触发,该方式通过带带跃迁可以产生大量的光生载流子。

    图  5  电流值随不同偏置电压的变化
    Figure  5.  Current value with different bias voltage

    通过对上述制备的1 mm间距的高纯半绝缘4H碳化硅光导开关进行不同偏置电压下电流波形测试,并且得到触发能量10 mJ偏置电压与电流的关系,以及偏置电压为6 kV,触发激光能量与电流的关系(如图 5所示)。可以发现随着偏置电压的增加,产生的电流值逐渐增加,偏压值与电流值呈线性增加的关系。对于工作在线性区的碳化硅光导开关来说,当光能量从4.5 mJ开始时,随着光能量的变换,电流值变化不大,分析认为此时光生载流子浓度趋于稳定(如图 6所示)。因此结合偏置电压和激光能量与电流值的关系,可以得出当光生载流子浓度一定时,偏置电压决定了PCSS的导通电流值。

    图  6  电流值随不同激光能量的变化
    Figure  6.  Current value with different laser energy

    对上述PCSS工作在偏压6 kV,触发能量10 mJ,导通电阻小于1 Ω,工作频率10 Hz的条件下,测试其电流变化情况,器件一直工作到失效状态。其电流值的变化如图 7所示,当开关触发次数为651次时,器件彻底失效,利用激光共聚焦显微镜以及扫描电子显微镜对PCSS进行表面失效性分析,具体结果如图 8所示,可以看出,在电极之间形成了如树枝状的贯穿性损伤通道,认为这是影响器件失效的根本原因。

    图  7  器件寿命图
    Figure  7.  Device lifetime diagram
    图  8  PCSS失效形貌图
    Figure  8.  Images of cracks in the failed PCSS

    大直径高质量的碳化硅衬底依然是目前国内外研究的热点,在当前节能减排的国际形势下,快速稳定地突破衬底尺寸和质量等关键问题,才能够更好地占据碳化硅市场,并且降低碳化硅器件的制作成本,真正将碳化硅材料商用到各个领域;此外在高功率微波领域,利用高纯半绝缘4H-SiC材料制备的低导通电阻光导开关也是当前研究的热点,结合碳化硅材料生长调控手段,制备适用于光导开关的特定杂质浓度的材料也应该作为目前碳化硅光导开关商用化的研究重点。随着碳化硅单晶生长技术及器件制造技术的进一步完善,碳化硅基功率器件会在成品率、可靠性和价格方面获得大大改善,将对电力电子带来重要的技术革新。

  • 图  1  电子FLASH辐照平台示意图(俯视图)

    Figure  1.  Schematic diagram of electron FLASH radiotherapy platform (top view)

    图  2  胶片标定结果

    Figure  2.  Film calibration result

    图  3  不同标定曲线之间测量剂量的相对偏差

    Figure  3.  Relative deviation in measured dose between different calibration curves

    图  4  不同剂量曲线的剂量分布测量结果

    Figure  4.  Measurement results of dose distribution for different dose curves

    图  5  电子束剂量分布特征

    Figure  5.  Electron beam dose distribution characteristics

  • [1] 高峰, 曹璐璐, 羊奕伟, 等. 基于PARTER开展肿瘤Flash-RT研究——设计及计算[J]. 中国医学物理学杂志, 2020, 37(9):1081-1087 doi: 10.3969/j.issn.1005-202X.2020.09.001

    Gao Feng, Cao Lulu, Yang Yiwei, et al. Design and calculation of Flash-RT based on PARTER[J]. Chinese Journal of Medical Physics, 2020, 37(9): 1081-1087 doi: 10.3969/j.issn.1005-202X.2020.09.001
    [2] Ramish Ashraf M, Rahman M, Zhang Rongxiao, et al. Dosimetry for FLASH radiotherapy: a review of tools and the role of radioluminescence and cherenkov emission[J]. Frontiers in Physics, 2020, 8: 328. doi: 10.3389/fphy.2020.00328
    [3] McManus M, Romano F, Lee N D, et al. The challenge of ionisation chamber dosimetry in ultra-short pulsed high dose-rate Very High Energy Electron beams[J]. Scientific Reports, 2020, 10: 9089. doi: 10.1038/s41598-020-65819-y
    [4] Favaudon V, Lentz J M, Heinrich S, et al. Time-resolved dosimetry of pulsed electron beams in very high dose-rate, FLASH irradiation for radiotherapy preclinical studies[J]. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 2019, 944: 162537. doi: 10.1016/j.nima.2019.162537
    [5] Chabi S, Van To T H, Leavitt R, et al. Ultra-high-dose-rate FLASH and conventional-dose-rate irradiation differentially affect human acute lymphoblastic leukemia and normal hematopoiesis[J]. International Journal of Radiation Oncology·Biology·Physics, 2021, 109(3): 819-829.
    [6] Zhu Hongyu, Xie Dehuan, Wang Ying, et al. Comparison of intratumor and local immune response between MV X-ray FLASH and conventional radiotherapies[J]. Clinical and Translational Radiation Oncology, 2023, 38: 138-146. doi: 10.1016/j.ctro.2022.11.005
    [7] Vozenin M C, De Fornel P, Petersson K, et al. The advantage of FLASH radiotherapy confirmed in mini-pig and cat-cancer patients[J]. Clinical Cancer Research, 2019, 25(1): 35-42. doi: 10.1158/1078-0432.CCR-17-3375
    [8] Schüler E, Trovati S, King G, et al. Experimental platform for ultra-high dose rate flash irradiation of small animals using a clinical linear accelerator[J]. International Journal of Radiation Oncology·Biology·Physics, 2017, 97(1): 195-203.
    [9] Montay-Gruel P, Petersson K, Jaccard M, et al. Irradiation in a flash: unique sparing of memory in mice after whole brain irradiation with dose rates above 100 Gy/s[J]. Radiotherapy and Oncology, 2017, 124(3): 365-369. doi: 10.1016/j.radonc.2017.05.003
    [10] Petersson K, Jaccard M, Germond J F, et al. High dose-per-pulse electron beam dosimetry - A model to correct for the ion recombination in the Advanced Markus ionization chamber[J]. Medical Physics, 2017, 44(3): 11571167.
    [11] Gotz M, Karsch L, Pawelke J. A new model for volume recombination in plane-parallel chambers in pulsed fields of high dose-per-pulse[J]. Physics in Medicine & Biology, 2017, 62(22): 8634-8654.
    [12] Marinelli M, Felici G, Galante F, et al. Design, realization, and characterization of a novel diamond detector prototype for FLASH radiotherapy dosimetry[J]. Medical Physics, 2022, 49(3): 1902-1910. doi: 10.1002/mp.15473
    [13] Velthuis J J, Page R F, Purves T M, et al. Toward pulse by pulse dosimetry using an SC CVD diamond detector[J]. IEEE Transactions on Radiation and Plasma Medical Sciences, 2017, 1(6): 527-533. doi: 10.1109/TRPMS.2017.2750799
    [14] Khanna R, James R, Hugtenburg R. Diamond dosimeter development for real-time microdosimetry and its use in FLASH RT[J]. Physica Medica, 2022, 94 Suppl 1: S102.
    [15] Bourgouin A, Schüller A, Hackel T, et al. Calorimeter for real-time dosimetry of pulsed ultra-high dose rate electron beams[J]. Frontiers in Physics, 2020, 8: 567340. doi: 10.3389/fphy.2020.567340
    [16] Bisogni M G. Inoriganic scintillators development and test of LYSO detector prototype[R]. FRPT 2021.
    [17] Favaudon V, Caplier L, Monceau V, et al. Ultrahigh dose-rate FLASH irradiation increases the differential response between normal and tumor tissue in mice[J]. Science Translational Medicine, 2014, 6: 245ra93.
    [18] Niroomand-Rad A, Chiu-Tsao S T, Grams M P, et al. Report of AAPM Task Group 235 radiochromic film dosimetry: an update to TG-55[J]. Medical Physics, 2020, 47(12): 5986-6025. doi: 10.1002/mp.14497
    [19] Liu K, Jorge P G, Tailor R, et al. Comprehensive evaluation and new recommendations in the use of Gafchromic EBT3 film[J]. Medical Physics, 2023, 50(11): 7252-7262. doi: 10.1002/mp.16593
    [20] Lewis D F, Chan M F. Technical Note: on GAFChromic EBT-XD film and the lateral response artifact[J]. Medical Physics, 2016, 43(2): 643-649. doi: 10.1118/1.4939226
    [21] Palmer A L, Bradley D A, Nisbet A. Evaluation and mitigation of potential errors in radiochromic film dosimetry due to film curvature at scanning[J]. Journal of Applied Clinical Medical Physics, 2015, 16(2): 425-431. doi: 10.1120/jacmp.v16i2.5141
    [22] Borca V C, Pasquino M, Russo G, et al. Dosimetric characterization and use of GAFCHROMIC EBT3 film for IMRT dose verification[J]. Journal of Applied Clinical Medical Physics, 2013, 14(2): 158-171. doi: 10.1120/jacmp.v14i2.4111
    [23] van Battum L J, Hoffmans D, Piersma H, et al. Accurate dosimetry with GafChromic™ EBT film of a 6 MV photon beam in water: what level is achievable?[J]. Medical Physics, 2008, 35(2): 704-716. doi: 10.1118/1.2828196
    [24] Devic S, Seuntjens J, Hegyi G, et al. Dosimetric properties of improved GafChromic films for seven different digitizers[J]. Medical Physics, 2004, 31(9): 2392-2401. doi: 10.1118/1.1776691
    [25] Sukhikh E, Sukhikh L, Malikov E, et al. Uncertainty of measurement absorbed dose by GAFCHROMIC EBT3 dosimeter for clinical electron and photon beams of medical accelerators[J]. Medical Radiology and Radiation Safety, 2019, 64(4): 56-63.
    [26] Sharma M, Singh R, Robert N, et al. Beam quality and dose rate dependency of Gafchromic EBT3 film irradiated with therapeutic megavolt photon beams[J]. Radiation Measurements, 2021, 146: 106632. doi: 10.1016/j.radmeas.2021.106632
    [27] Ferreira B C, Lopes M C, Capela M. Evaluation of an Epson flatbed scanner to read Gafchromic EBT films for radiation dosimetry[J]. Physics in Medicine & Biology, 2009, 54(4): 1073-1085.
    [28] 周婉仪, 宫辉, 邱睿, 等. 金刚石探测器在Flash照射实时剂量测量中的应用[J]. 中华放射医学与防护杂志, 2023, 43(9):729-735 doi: 10.3760/cma.j.cn112271-20230216-00036

    Zhou Wanyi, Gong Hui, Qiu Rui, et al. Feasibility of diamond detector on Flash radiation dosimetry online[J]. Chinese Journal of Radiological Medicine and Protection, 2023, 43(9): 729-735 doi: 10.3760/cma.j.cn112271-20230216-00036
    [29] 翁邓胡, 徐海荣. 高能电子束放射治疗的研究进展[J]. 中国辐射卫生, 2011, 20(3):375-378

    Weng Denghu, Xu Hairong. Research progress in high-energy electron beam radiation therapy[J]. Chinese Journal of Radiological Health, 2011, 20(3): 375-378
    [30] Musolino S V. Absorbed dose determination in external beam radiotherapy: an international code of practice for dosimetry based on standards of absorbed dose to water; Technical Reports Series No. 398[J]. Health Physics, 2001, 81(5): 592-593.
  • 加载中
图(5)
计量
  • 文章访问数:  270
  • HTML全文浏览量:  140
  • PDF下载量:  33
  • 被引次数: 0
出版历程
  • 收稿日期:  2024-03-14
  • 修回日期:  2024-09-14
  • 录用日期:  2024-09-02
  • 网络出版日期:  2024-09-23
  • 刊出日期:  2024-11-08

目录

/

返回文章
返回