基于密度泛函理论的微气泡臭氧化诺氟沙星降解途径与毒性评估

张静, 于倩, 丁亮亮, 刘春, 闫晨旭, 牛建瑞, 张瑞娜. 基于密度泛函理论的微气泡臭氧化诺氟沙星降解途径与毒性评估[J]. 环境化学, 2025, 44(1): 288-298. doi: 10.7524/j.issn.0254-6108.2023122102
引用本文: 张静, 于倩, 丁亮亮, 刘春, 闫晨旭, 牛建瑞, 张瑞娜. 基于密度泛函理论的微气泡臭氧化诺氟沙星降解途径与毒性评估[J]. 环境化学, 2025, 44(1): 288-298. doi: 10.7524/j.issn.0254-6108.2023122102
ZHANG Jing, YU Qian, DING Liangliang, LIU Chun, YAN Chenxu, NIU Jianrui, ZHANG Ruina. Degradation pathway and toxicity assessment of microbubble ozonation of Norfloxacin based on Density Functional Theory[J]. Environmental Chemistry, 2025, 44(1): 288-298. doi: 10.7524/j.issn.0254-6108.2023122102
Citation: ZHANG Jing, YU Qian, DING Liangliang, LIU Chun, YAN Chenxu, NIU Jianrui, ZHANG Ruina. Degradation pathway and toxicity assessment of microbubble ozonation of Norfloxacin based on Density Functional Theory[J]. Environmental Chemistry, 2025, 44(1): 288-298. doi: 10.7524/j.issn.0254-6108.2023122102

基于密度泛函理论的微气泡臭氧化诺氟沙星降解途径与毒性评估

    通讯作者: E-mail:liuchun@hebust.edu.cn
  • 基金项目:
    国家重点研发计划项目(2022YFE0119000)资助.
  • 中图分类号: X-1;O6

  • CSTR: 32061.14.hjhx.2023122102

Degradation pathway and toxicity assessment of microbubble ozonation of Norfloxacin based on Density Functional Theory

    Corresponding author: LIU Chun, liuchun@hebust.edu.cn
  • Fund Project: the National Key Research and Development Program (2022YFE0119000).
  • 摘要: 诺氟沙星(Norfloxacin,NOF)为第三代喹诺酮类抗菌药被广泛应用,易在环境中积累,威胁人类健康和环境安全. 本研究采用Gaussian与Toxicity Estimation Software Tool(T.E.S.T.)软件预测在微气泡臭氧化(MB/O3)体系中NOF的降解途径与降解产物的生态毒性. 结果表明,MB/O3体系中的羟基自由基等活性氧物质易对NOF分子结构上键级较小、平均局部离子化能较低以及简缩Fukui函数较高的反应位点进行攻击;液相色谱-串联质谱法(LC-MS/MS)分析结果结合模型计算推测,MB/O3处理NOF降解途径包括:脱氟反应、哌嗪环转化、羟基取代反应和喹诺酮环氧化;利用T.E.S.T.软件评估NOF及降解产物的毒性,发现NOF对口服大鼠、胖头鱼、水蚤均存在负面影响,具有发育毒性,无致突变性,且降解过程中可能会产生比NOF毒性大的产物. 本研究采用模型预测NOF降解途径与生态毒性的方法便捷有效,可为MB/O3体系降解复杂有机分子的降解途径和生态毒性研究提供理论依据.
  • 2030年前确保全球范围内实现人人都能获得安全和负担得起的饮用水,是联合国可持续发展目标(Sustainable Development Goals, SDGs)中第六项“清洁饮用水和卫生设施”(SDG6)的主要内容之一. 现今水污染形势日益严峻,可供人们直接利用的液态水源均存在不同程度的污染,全球约有22亿人的日常饮用水无法得到安全保障[1],预计到21世纪中期甚至将有近5亿人常年面临严重的淡水资源短缺问题[2]. 基于KIM模型[3],谷歌公司与世界卫生组织、联合国儿童基金委员会的共同研究表明,不同于海水淡化技术的高费用低回报和技术的区域性限制,空气取水技术没有地域限制,能经济有效地在全球范围内提供液态水源,同时不会对全球水生态水循环造成较大影响[4]. 以我国西北部干旱地区为例,1960—2015年间,区域内平均相对湿度达50.34%[5],当夜晚温度为15 ℃时,空气中水分子含量约为6.4 g·m−3. 空气中丰富的水资源(总量约12.9万亿吨)保证了空气取水技术能为干旱少雨地区、灾害重建地区和军队野外行动有效提供液态水源.

    为有效实现空气取水,基于吸附材料的吸附式空气取水技术是当前领域的研究热点,金属有机框架(Metal-Organic Frameworks, MOFs)材料因其优异的水吸附性能和循环特性受到了广泛关注. 21世纪初,Yaghi课题组首次了报道MOF-5材料的合成[6],并用二级结构单元(Second Building Unit, SBU)多金属簇理论成功解释MOFs材料的合成,为可控合成MOFs材料指明了方向;之后软硬酸碱理论(Hard-Soft-Acid-Base, HASB)的成功应用则为水稳定性MOFs材料的制备提供了理论基础[7],,极大拓宽了MOFs材料的应用性能. MOFs材料因具有高比表面积、高孔隙率和功能高度可控等优异特性,逐渐从水吸附性能表征拓展到空气取水实际应用[8]. 2012年,Seo等[9]率先以铁基MIL-100(Fe)和铬基MIL-101(Cr)两类MOFs材料为基础进行水吸附实验,结果表明MOFs材料在空气除湿和淡水获取领域具有极大的应用前景. 近年来以Yaghi课题组为代表的研究人员以MOF-303(Al) [10]、MOF-801(Zr) [11]、MIL-160(Al)[12]和PC-MOF(Cr)[13]等材料为研究对象,对MOFs吸附材料进行了一系列空气取水技术和装置的试验研究,包括在沙漠等干旱低湿地区的现场试验,充分证实了基于MOFs材料空气取水潜力. 与传统吸附材料相比,基于MOFs材料的吸附式空气取水装置不仅可在环境湿度较低时仍能稳定持续地从空气中获取水分,而且MOFs材料因其稳定的框架结构可作为传统吸湿盐的良好载体,抑制吸湿盐的潮解问题,增强复合材料的水吸附性能以满足实际环境的工作需求.

    本文通过梳理近年来面向空气取水的MOFs材料研究现状,简述空气取水的基本原理和方法,总结MOFs材料的主要发展历程、研究和应用进展,并结合实际需求对面向空气取水的MOFs材料未来发展进行展望.

    大气中水分子主要以液态(如雾滴和露水)和气态两种方式存在,当环境中水蒸气分压达到饱和或超出饱和时,空气中水分子会在凝结核上凝结成云滴或冰晶[14]. 空气取水技术是通过外部作用,如降低环境温度或增加水分子总量,在一定区域环境内使空气相对湿度(relative humidity, RH)趋近或大于100%,液化空气中水分子以达到空气取水目的[15].

    基于空气取水原理,目前空气取水技术的基本方法主要分为3种,如表1所示. 1)自然聚雾法直接收集空气中已存在的雾滴或露水[16];2)冷凝集露法基于表面冷却技术在装置表面收集空气中的水分子[17];3) 吸附-解吸法基于吸附材料富集空气中的水分子,在冷凝管中实现水分子收集[17].

    表 1  空气取水技术主要方法
    Table 1.  The main method of Atmosphere water harvester
    主要方法Main method主要类别Main category工作特性Operating characteristics
    自然聚雾法网状装置[16]气候条件决定取水性能
    冷凝集露法主动冷凝[24]外部供能决定装置运行
    辐射冷却[25]无需外部供能可自动工作
    吸附解吸法吸湿盐[32]吸附性能好但存在潮解泄漏危险
    水凝胶[35]吸附性能好但相对湿度要求较高
    MOFs[7]吸附性能好且相对湿度要求较低
     | Show Table
    DownLoad: CSV

    自然聚雾法主要是利用网状装置收集空气中的雾滴和露珠 [16]. 在秘鲁,Lummerich等比较了5种形状各异且材料不同的雾收集装置,基于拉舍尔面料(Raschel)的“埃菲尔”状雾收集器(Eiffel Fog Collectors)的实际运行效果最好[18]. 在雾收集装置的基础上,Park等的研究表明,控制装置的网格大小、增加装置的阴影面积,可有效提高装置的雾收集效率[19];将材料由单根铜丝替换成具有纳米疏水和超亲水特性的三维铜丝,能有效减少网状装置在雾收集过程中的堵塞和液滴脱落问题,雾滴传输效率提高了20%以上 [20]. Feng等基于仿生技术,制备了光滑的超亲水表面的雾自动收集装置(nephrolepis auriculata-inspired patterned surface, NAIPS),该装置取水效率可达(519.84±17.04) g·m−2·d−1,较常规均质光滑表面雾收集装置提高了139%[21]. 由于应用场景和工作条件的限制,雾收集装置只能在常年潮湿多雾地区实现高效率取水,并不能满足大部分地区的空气取水要求.

    冷凝集露法是通过将装置表面温度降低至环境露点以下,收集冷凝液化的水分子[22]. 以一个多山岛屿(Punaauia, Tahiti Island)和一个环礁(Tikehau, Tuamotu Archipelago)为例,Clus等[23]评估了在潮湿热带地区的旱季露水作为饮用水源的潜力,发现旱季岛屿的露水量相当于降水量的两倍,具有极大的开发潜力. 根据冷凝驱动力的不同来源,冷凝方法分为主动冷凝和被动冷凝. 主动冷凝是依靠外部能源输入,通过蒸汽压缩或者基于珀耳帖效应(Peltier effect)[24]的热电转化设备进行可控的空气制冷;被动冷凝(又称辐射冷却)无需外部能源输入,通过材料选择使装置的整体辐射能量大于吸收热量,进而自动降低其表面温度[25]. 研究表明,多层辐射冷却装置较单层装置能有效减少其工作过程中的热损失[26-27],通过将辐射冷却与超疏水冷凝采集器协同运行,露水采集量高达1200 g·m−2·d−1[28]. 然而,目前该方法的实际应用大多限于露点温度与环境温度相差较小的温暖潮湿地区[29],需进一步开发适于露点温度与环境温度差值较大情形的取水装置,以满足当地群众或沙漠绿化的用水需求将是未来的研究重点.

    吸附解吸法是通过吸附材料充分吸附空气中的水分子,待吸附饱和后依靠外部能源输入(如太阳能)实现吸附材料中水分子的解吸,释放的水分子通过冷凝管完成液化收集[30-31]图1).

    图 1  基于太阳能的吸附式空气取水装置工作示意图[31]
    Figure 1.  Schematic diagram of the adsorption-atmosphere water harvester based on solar energy[31]

    吸附材料决定了吸附式空气取水装置的理论效率[30]. 赵亚等的研究表明,以无水氯化钙为代表的吸湿盐材料的水吸附效率与容量均较为突出[32],将吸湿盐与硅胶[33]、活性炭纤维[34]等进行复合可明显提升材料的水吸附性能,但吸湿盐遇水易潮解泄漏的问题并未得到解决,这极大限制了其实际应用潜力. 水凝胶材料的水吸附容量高且结构稳定性好[35],Shi等将一种新型分层三维微结构水凝胶材料应用于汽态雾滴收集,户外测试结果表明产水量达到了34 L·m−2·d−1,但水凝胶对其水吸附过程中环境相对湿度的要求较高[36],大都要求70%以上的相对湿度. MOFs吸附材料的结构可控且水吸附性能较好[7],Hyunho等[11]进行MOF-801(Zr)材料的模拟计算,结果显示仅依靠太阳能供能,其产水量可达2.8 L·kg−1·d−1;Yaghi课题组[10]报道的MOF-303(Al)材料,不仅能在白天持续工作实现多次集水循环,而且在莫哈韦沙漠(The Mojave Desert, Northern America)相对湿度低至10%的极端条件下其产水量仍能达到0.7 L·kg−1·d−1. 上述研究结果充分说明基于MOFs吸附材料的空气取水技术在低湿度环境下具有巨大的应用潜力. 考虑到我国先进的光伏发电技术及其发展态势,通过技术集成,基于MOFs材料的空气取水技术不仅完全有可能实现能量自给[37]实现自动化运行,而且能赋予光伏发电技术新的生态环境效益.

    MOFs材料是指金属离子和有机配体通过化学键配位作用实现自组装连接的周期性三维网状结构[38],自Yaghi基于多元金属簇理论的二级结构单元(SBU)概念成功解释MOF-5[39]的合成后,水稳定性MOFs材料的选择合成愈益受到关注[40]. 软硬酸碱理论的成功应用则为金属源和有机配体的选择提供了理论支撑[7], 以溶剂热法[41-42]、微波辅助加热[43]、电化学法[44]和机械化学法[45]等为主的合成方法为MOFs材料研发提供了简便制备途径,并形成多种MOFs系列材料,表2列举了部分MOFs系列材料名称及相关信息. 以2022年5月在Web of Science 核心数据库的数据为例,对MOF和 application所有字段进行搜索,VOSviewer软件[46]的关键词共性分析结果表明,吸附性能是当前MOFs材料研究的主要方向,而其吸附性能也是吸附式空气取水装置的关键.

    表 2  MOFs系列材料
    Table 2.  The series of MOFs materials
    系列名称Name代表性材料Representative materials金属离子Metal ions有机配体Organic ligands首发单位Inventive institution主要应用领域Main applications
    IRMOFIRMOF-1[39]Zn2+对苯二甲酸等美国Yaghi课题组氢气吸附
    ZIFZIF-100[47]Zn2+、Co2+咪唑等美国Yaghi课题组CO2储存
    MILMIL-53[48]Al3+、Cr3+、Fe3+对苯二甲酸等法国Ferey课题组甲醇吸附
    CPLCPL-1[49]Cu2+吡嗪等日本Kitagawa研究组甲烷吸附
    UiOUiO-66[50]Zr4+对苯二甲酸等挪威奥斯陆大学反应催化
    PCNPCN-9[51]Cu2+三嗪等美国迈阿密大学吸附分离
    HKUSTHKUST-1[52]Cu2+间苯三甲酸等香港科技大学CO2吸附
     | Show Table
    DownLoad: CSV

    同样以2022年05月份数据为例,在Web of Science 核心数据库以water harvesting、water harvester、AWH和AWG所有字段进行搜索并采用VOSviewer软件[41]进行关键词共性分析,如图2所示,结果主要分为4部分,water harvesting 与metal-organic frameworks二者占据其中两部分的中心位置,且二者与adsorption和atmosphere water harvesting两节点均存在较大的共性关系,这表明基于MOFs吸附材料的空气取水技术已受到广泛关注. 对于吸附式空气取水装置,吸附材料需满足水热稳定性好、水吸附容量高、解吸温度低、吸附-解吸过程迅速和循环特性好等特性[53]. 研究表明,MOFs材料因其优异的水热稳定性、超高的比表面积孔隙率和快速的吸附动力学不仅可满足空气取水的应用需求[54],而且其功能可控性能够实现在不同环境条件下定向合成适宜的吸附材料[55].

    图 2  空气取水装置相关文献关键词的共性关系
    Figure 2.  Common relationships of keywords in the literature related to atmosphere water harvester

    如上所述,二级结构单元(SBU)的引入和软硬酸碱理论(HSAB)的应用为水稳定性MOFs材料合成指明了方向. 现今常用的羧酸盐有机配体(对苯二甲酸、反丁烯二酸、吡唑二羧酸等)去质子化后可被视为硬碱,而高价金属离子(Cr3+、Fe3+、Al3+、Zr4+、Ti4+等)是硬酸,自组装过程中,金属离子会率先形成SBU结构(如Zn4O、Zr6O8等),然后与有机配体通过强作用键连接形成水稳定性MOF材料;同时,亲水性好的羧酸盐有机配体利于实现水分子的选择性吸附,金属位点的不饱和性与材料的多孔特性能给予水分子充分的吸附位点,促进水分子的吸附[7]. 此外,Yaghi提出的网状化学(Reticular Chemistry)理论[56-57],能在确保MOFs材料整体结构稳定性的同时实现多元官能团的掺杂和金属置换的设计合成,从而将不同类型官能团和金属离子稳定结合到MOFs结构中,实现MOFs材料性能的精准调控.

    在水分子吸附过程中,MOFs材料通过SBU结构中多元金属簇和有机配体之间的强作用力键维持材料整体结构的稳定[58],水分子则通过与有机配体和不饱和金属位点的结合在材料孔隙表面实现物理吸附和化学吸附[53,59]图3),MOFs多孔材料同时也会发生毛细管冷凝以促进水分子的进一步吸附[60]. 基于化学吸附和毛细管冷凝储存的水分子在解吸过程中往往需要更大的能量输入[61-62],因此,为减少解吸过程的能量输入,增强水分子的物理吸附(如亲水性官能团的掺杂)、控制开放金属位点比例和孔隙结构是平衡材料稳定结构和优异性能的重要方法. 同时,MOFs材料巨大的比表面积和超高的孔隙率[63]为水分子提供了足够的吸附位点和存储空间. 已有研究表明,环境相对湿度较低时,材料水吸附容量主要由MOFs材料的比表面积决定;而环境相对湿度较高时,其则由材料的孔隙率和总容积控制[64].

    图 3  水分子吸附机理:物理吸附(a)、化学吸附(b)[53]
    Figure 3.  Adsorption mechanism of water molecules: Physical adsorption(a), chemical adsorption(b) [53]

    2021年10月,Yaghi课题组通过单晶衍射和密度函数理论 (density functional theory, DFT)模拟成功证明了MOF-303(Al)在水分子吸附过程中,首个水分子吸附位点是在材料骨架内的吡唑之间,水分子通过与两个吡唑基团和一个u2-OH基团形成3个氢键实现吸附 [65],这充分说明有机配体的选择对水分子吸附过程具有重要作用. 亲水性有机配体能有效增强材料对水分子的选择性吸附[66],以MIL-101(Cr)为例,材料内部的水传输动力学强烈依赖于自身独特的孔结构和亲水性空间分布特征[67],选择适宜的亲水性官能团(如氨基:—NH2[66])对材料结构进行修饰优化,可调控MOFs材料的水吸附性能以适应不同的应用环境[68-69].

    为有效实现空气取水, MOFs材料的水吸附等温线呈现明显的S型和较小的迟滞回线,这保证了材料在较窄湿度范围内具有水分子的快速吸附-解吸并保持良好的循环特性[62]. 在低湿度地区,为保障有效吸附空气中水分子,MOFs材料的水吸附动力学曲线拐点值必须低于环境湿度值[70]. Rieth等研究证明,将Ni2Cl2BTDD材料桥联结构中的氯离子替换为半径更大的溴离子,离子替换后合成的Ni2Br2BTDD材料能在25%低相对湿度下达到0.64 g·g−1的高吸附量,并能在400次水吸附-解吸循环中表现出持久的稳定性[71]. 进一步的研究证明,当MOFs材料存在强成核位点时,增加材料的亲水性能只能增强吸附的初始速率,并不能改变水吸附曲线拐点,而适当缩小材料的结构孔径或改变孔隙中心水分子的极性,可有效提升材料在低相对湿度下的吸水性能[72]. 将氯离子替换为溴离子后,MOFs材料的水吸附动力曲线基本保持不变(符合国际定义的Ⅳ类吸附曲线[73]),但其吸附拐点会向低相对湿度值移动(图4).

    图 4  MOFs水吸附动力学曲线[71]
    Figure 4.  The water adsorption kinetic curve of MOFs[71]
    (a), 在25 ℃下测量的Ni2Br2BTDD、Ni2Cl2BTDD、Ni2F0.83Cl0.17BTDD和Ni2(OH)2BTDD的水蒸气吸附(实心点)和解吸(空心点)等温线;(b) 使用材料密度将吸附等温线转换为体积单位
    (a) Water vapor adsorption (closed symbols) and desorption (open symbols) isotherms of Ni2Br2BTDD, Ni2Cl2BTDD, Ni2F0.83Cl0.17BTDD, and Ni2(OH)2BTDD measured at 25 ℃;(b) Adsorption isotherms converted to volumetric units using the material density

    将MOFs材料与传统高效吸湿盐(如LiCl、CaCl2等)复合,利用MOF材料稳定的骨架结构负载吸湿盐,可有效提升材料的整体吸附性能. Hu等[74]成功将CaCl2纳米晶体负载在水吸附性能较弱的铁基MOFs材料内部, MOFs不仅能有效防止CaCl2在水吸附过程中的潮解,而且提高了材料的水吸附性能,其水吸附容量高达2.685 g·g−1(以CaCl2为主),是原始CaCl2粉末水吸附容量的2.23倍. 将吸湿性MOFs材料富马酸铝与CaCl2进行复合,虽然富马酸铝基体的比表面积在加入CaCl2后急剧下降,但材料的水吸附容量从纯富马酸铝的0.4 g·g−1增加到0.68 g·g-1[75]. 上海交通大学王如竹课题组[76]成功将LiCl负载在MIL-101(Cr)基质中,该复合吸附剂整合了吸湿盐典型的三步水吸附性能:化学吸附、潮解和溶液吸收. 在典型低相对湿度工作条件下(30 ℃,30%RH),复合吸附剂的吸水容量高达0.77 g·g−1. 尽管MOFs材料能较好抑制吸湿盐的潮解泄漏,但为保障用水安全,基于复合吸附剂获得的液态水源的水质检测仍必不可少,例如,以满足生活饮用水的总硬度为例,水中Ca2+浓度需低于180 mg·L−1.

    在水分子吸附过程中,尽管MOFs材料能保持较好的孔隙率和结晶度,但由于客体-主体相互作用(如水分子诱导的键重排),材料内部仍会发生重要的分子水平结构变化[77]. Choi等基于MOF-801(Zr)开展了结构缺陷对比实验,发现MOFs材料结构中的缺陷密度是决定水吸附性能优劣的因素之一,水分子在纳米多孔MOF-801(Zr)材料上的吸附过程会优先沿着(110)方向进行[78],但通过结构缺陷控制材料性能往往会严重破坏材料的整体稳定性[79],进而影响材料的循环特性. Krajnc等发现,LTA(Linde Type A)型拓扑结构的沸石状磷酸铝不仅能够在10%—15%低相对湿度下表现出优异的水吸附性能,吸附容量高(0.42 g·g−1),而且材料的解吸温度与MOF-801(Zr)相比存在10–15 ℃的降幅,解吸能量需求较低. 40次水吸附-解吸循环的测试结果表明,材料总吸附容量下降程度不到2%,这说明LTA型拓扑结构的沸石状磷酸铝的循环稳定特性较好[80],具有长期工作的可能. 基于多金属氧酸盐(Polyoxometalates,POM)结构,Zhu等首次将3种不同类型的有机配体同时引入到金属骨架结构,发现材料拥有非常优异的吸水性能,同时多种有机配体赋予的功能性也增强了MOFs材料的环境适应性[81]. 在网状化学基础上,Towsif等构建了一种全新的稳定性MOF材料“Cr-soc-MOF-1”,材料具有极高的孔隙率(比表面积达4549 m2·g−1),在70%相对湿度下水吸附容量高达1.95 g·g−1,且在100多次吸附-解吸循环后,材料水吸附容量未出现较大降幅[82].

    在MOFs材料合成过程中,样品最终大都呈现粉末状态. 为应用于实际,粉末状样品需均匀负载在吸附床等基底上才能保证空气取水装置的正常运行[83]. 为增强MOFs材料的适用性能和均匀负载,样品制备过程可选择合适的基底(如铝箔[84]),将MOFs材料直接在基体进行沉淀以实现均匀负载的整体式吸附材料,并将其直接应用于实际环境中. Tan等以富马酸铝为例,以玻璃纤维纸作为前驱体溶液的沉积基底,在富马酸/水物质的量比为0.02、50 ℃反应温度、60 min反应时间下,成功获得了BET比表面积为740.15 m2·g−1和吸附容量0.3906 g·g−1的高效吸附材料,检测结果表明样品中含有76%的富马酸铝功能性粉末[85]. 除去玻璃纤维纸,水凝胶也可作为MOFs材料的合成基质形成复合材料以便于直接应用,同时由于水凝胶自身良好的吸水性能,这也有利于材料在高湿度范围下的应用[13].

    总之,基于MOFs材料的吸附式空气取水装置能为部分地区液态水源缺乏问题提供一种极具潜力的解决方案[86],部分已被证明可用于空气取水的MOFs材料如表3所示. 但目前MOFs材料在制备过程中经常使用二甲基酰胺(N,N-Dimethylformamide,DMF)等有毒有机溶剂和有毒重金属元素,因此,筛选、开发绿色低成本的金属源和生物安全性高的有机配体是实现MOFs材料绿色合成、经济化转型和商业化发展的研究重点,无有机溶剂参与的水热合成法、无溶剂参与的机械合成法和整体式吸附材料的制备是未来MOFs材料合成的主要趋势.

    表 3  可用于空气取水的MOFs材料
    Table 3.  MOFs materials that can be used for Atmosphere water harvester
    名称Name吸附转折点Adsorption turning point水吸附容量/(g.g−1)Water adsorption capacity孔体积/(cm3.g−1)Pore volume比表面积/(m2.g−1)Specific surface area稳定性Stability
    Cr-soc-MOF-1 [82]0.691.952.14549无损失(>100次循环)
    MIL-101(Cr)[9]0.41.554150无损失(100 ℃水热处理7 d)
    MIL-100(Fe)[9]0.350.790.822300无损失(100 ℃水热处理7 d)
    NiBr2(BTDD)[71]0.250.641770无损失(>400次循环)
    CoCl2(BTDD)[72]0.280.82191263%损失(>30次循环)
    MOF-801(Zr)[87]0.090.360.45无损失(>5次循环)
    MOF-808(Zr)[88]0.30.590.84不稳定
    MOF-841(Zr)[11]0.220.510.537%损失(>5次循环)
    Al-flumarate[89]0.270.450.481080无损失(>4500次循环)
    MOF-303(Al)[10]0.130.450.541119无损失(>5次循环)
    MOF-333(Al)[65]0.220.440.4812803%损失(>2000次循环)
    CAU-23(Al)[90]0.300.370.481250无损失(>5000次循环)
    MIL-160(Al)[91]0.090.380.401150无损失(>10次循环)
     | Show Table
    DownLoad: CSV

    吸附性空气取水装置(adsorption-based atmosphere water harvester, ABAWH)的工作步骤主要分为3部分[92-93]:1)水分子吸附存储:水分子在吸附材料表面和孔隙内的吸附存储;2)水分子解吸:外部能源输入促使吸附材料内存储的水分子解吸;3)水分子液化收集:水分子在冷凝装置中的液化收集. 研究表明,在温暖潮湿条件下,空气取水装置工作的平均能耗为1.02 kWh·L−1,在寒冷潮湿气候则高达6.23 kWh·L-1[94]. 为降低能耗,依托光伏发电技术的空气取水装置(AWH)是当前研究的主要方向[95-96],同时以太阳能技术等为基础确保取水装置的能量自给是实现吸附式空气取水装置工作自动化的尝试方向.

    太阳能吸附式取水装置在日照时间充足的情况下每天至少能实现一次取水循环(图5a):夜晚温度低(10 ℃—15 ℃)、相对湿度高(40%), MOFs吸附材料能够快速吸附空气中的水分子达到饱和状态;白天温度较高(35 ℃—40 ℃)、相对湿度低(5%),利用太阳能直接加热吸附剂,可促进储存水分子的快速解吸脱附. Yaghi课题组基于锆离子(Zr4+)和富马酸有机配体,制备得到了水吸附性能优异的MOF-801(Zr)[11,87],基于循环工作状态的吸附取水装置模拟计算结果表明,装置日产水量约为2.8 L·kg−1·d−1. 当装置安装在沙漠地区时,为最大限度吸收太阳能和提高装置的取水效率,需将其安装在坡度37°的支架上,并以高热容量和低热导率的土壤作为冷凝器外部隔热材料以降低冷凝器表面温度. 基于昼夜单循环的吸附数据,进一步模拟结果表明若吸附初始相对湿度(RH)为20%,装置可实现0.2 L·kg−1·d−1的日产水量. 为进一步扩展装置的适用性能,进行储能装置(如光伏储能)和表面冷却技术的集成,可实现装置取水效率和循环工作能力的提升,在实际光照较弱时持续工作.

    图 5  (a)基于太阳能的单循环系统(MOF-801)[70];(b)基于太阳能的多循环系统(MOF-303(Al))[10];(c)基于太阳能的翻转式系统(UiO-66NH2(TUN))[98];(d, e)基于太阳能的连续取水系统(PC-MOF)[13]
    Figure 5.  (a) Solar-based single-cycle system(MOF-801)[70]; (b) Solar-based multi-cycle system(MOF-303(Al))[10] ; (c) Solar-based flip system(UiO-66NH2(TUN))[98]; (d, e) Continuous water harvester based on solar energy(PC-MOF)[13]

    MOFs材料的等温吸附曲线呈现明显的S型,能在较窄相对湿度范围内实现水蒸汽快速吸附解吸,若装置仅支持每日单循环工作,则浪费了材料优异的水蒸气吸附解吸性能. 进一步研究表明,基于水吸附动力学曲线控制好材料的吸附-解吸时间,最大化利用材料水蒸汽吸附的快速动力学特征,吸附取水装置完全可以在一天内实现多次集水循环[70]. Yaghi课题组以三价铝离子(Al3+)和3,5-吡唑二羧酸有机配体为原料制备了新型MOFs材料:MOF-303(Al)[10, 97]. 为充分利用材料水吸附性能,吸附装置在强对流风扇作用下快速实现水分子的吸附解吸,结果表明,在室内稳定环境为32%相对湿度和27 ℃的条件下连续工作24 h,装置完成了9次取水循环,日产水量为1.3 L·kg−1·d−1. 将装置安装在北美的莫哈韦沙漠中进行为期3 d 72h的取水工作(图5b),结果表明,即使在85%的时间内环境露点温度都在5 ℃以下,装置仍能持续保持每天7次取水循环和0.7 L·kg−1·d−1的日产水量. 但装置由于部件较多,占用空间较大,若依托光伏板进行集成,开发光伏板下的多级连续空气取水装置是一种潜在的应用方向,既能有效促进光伏厂区生产生活用水的自给,也能提升厂区的经济效益,甚至有利于厂区周围的土地灌溉或沙漠绿化.

    为充分利用装置的吸附-解吸循环的双步骤工作状态,使之连续产水,Wu等基于Ti3C2掺杂的UiO-66NH2(TUN)材料设计了一种新型可翻转吸附层的吸附式空气取水装置(图5c). 模拟室内低湿环境(20%相对湿度,298 K),装置能在阳光照射下保持0.6 L·kg−1·d−1的日产水量持续产水 [98]. 一般而言,材料饱和吸附与完全解吸所用时间并不一致,基于水吸附动力学曲线,计算出吸附时间/解吸时间的比值,以此确定吸附层滚筒的数量,可最大化利用材料吸附解吸性能,实现高效持续取水. 进一步研究表明,复合材料具有显著的水吸附能力和可通过坚固的多孔水通道诱导的直接释水性能,其中MOFs材料作为催化剂能加快装置的吸附-解吸速率,同时作为吸附剂可有效增加装置的水吸附容量[76]. 为实现真正意义上的连续供水,Yilmaz等将连有水活性位点的聚合物(基于吸湿盐接枝改性)与水稳定MOF材料进行复合得到了PC-MOF[13],结果表明在相对湿度低于60%时,装置以水分子的吸附储存为主;当相对湿度大于60%时,装置可连续获取液态水(图5d);在相对湿度(RH)达到 90%时,吸附装置能以95%的总液体输送效率和71%的自主液体输送效率实现6 g·g−1·d−1的日产水量,并连续供水1440 h(图5e). 为保证液态水的连续获取,装置存在相对湿度的限制(>60%),选择研发气体分离膜等方法提升环境湿度以满足装置需求,将出口高湿度空气作为PC-MOF的工作氛围,可实现低湿度条件下水分子的液化收集.

    尽管基于MOFs材料的吸附式空气取水装置已取得了较多成功的实验室和实际现场数据(如表4),以美国Yaghi教授和国内上海交通大学王如竹教授为代表的研究团队在空气取水领域做了诸多开拓性的研究,但迄今为止,基于MOFs材料的吸附式空气取水装置并未规模化应用[99]. 吸附装置与表面冷却技术的集成是一种值得深入尝试的设计方案,将辐射冷却[100-101]和热电冷却[102-103]方法与装置进行结合,能够在不损耗材料吸附性能的前提下形成相对较高的湿度环境,进而促进水分子在材料表面的吸附. 同时,空气取水技术实际持续工作的前提是满足环境相对湿度的要求,基于MOFs材料的膜分离技术是一种理想方法:通过膜分离法富集环境空气中的水分子,环境湿度提升后,可直接采用冷凝装置收集液化的蒸汽,也可基于高湿度吸附材料(如水凝胶)完成连续取水,以实现低湿度环境下的空气取水目标. 从能耗角度出发,膜分离法取水装置不仅比一般吸附装置可降低50%甚至更大比例的能耗[104],而且可充分利用低压再循环吹扫膜表面,从而以最节能的方式获得较大的产水率[105].

    表 4  空气取水装置的实际应用
    Table 4.  Practical devices for Atmosphere Water Harvesting
    名称Name吸附材料Adsorption material工作特点Characteristics实际取水效率Actual water harvesting efficiency不足Inadequacies潜在改进方法Potential improvements
    单循环装置MOF-801(Zr) [11,87]夜晚湿度较高充分吸附白天依靠光照实现解吸0.2 L·kg−1·d−1(20%RH)吸附时间较长长,装置只能在白天进行解吸集水集成储能装置、表面冷却技术的结合
    LiCl@MIL-101(Cr) [76]0.45 L·kg−1·d−1(30%RH)
    多循环装置MOF-303(Al) [10, 97]循环时间较短,24 h内可进行多次取水循环1.3 L·kg−1·d−1(27%RH)装置较大,各部件较为分散基于光伏产区实现装置集成
    连续取水装置Ti3C2@ UiO-66NH2(TUN) [98]吸附床轮流吸附解吸,实现连续取水0.6 L·kg−1·d−1(20%RH)吸附解吸速率间并不一致,简单分为双侧易存在空窗期计算吸附解吸时间比值,增加吸附床
    PC-MOF[13]直接水释放,实现连续取水6 L·kg−1·d−1(90%RH)水释放湿度要求高,低湿度地区不适宜增加空气富集装置,提升空气相对湿度以适用低湿度区域
     | Show Table
    DownLoad: CSV

    MOF材料的高比表面积和大孔隙率为水分子充分吸附提供了充足的负载位点,材料的多金属簇结构确保在吸附过程中整体的稳定性. 不同于吸湿盐材料的近线性吸附,MOFs材料的水吸附等温线呈现明显的“S”形且吸附拐点所需的环境湿度相对较低,这有利于材料在低湿度地区充分吸附水分子获取液态水源,但同时也限制了单一MOFs材料在高湿度区域的应用前景. 为增强MOFs材料的适用性能和推广应用,通过前文分析,本文从材料构建和修饰、复合材料制备和装置优化设计等几个角度出发进行分析.

    (1)材料构建和修饰. 在MOFs材料的水吸附过程,对称性较低的杂环有机配体在材料合成过程中更能凸显有机配体的极性并产生较多的结构缺陷,在保持结构稳定基础上有效增强材料的水吸附性能,因此可在制备过程中多进行杂环有机配体的选择尝试. 同时,绿色金属源、廉价有机配体的选择、有机溶剂的减少使用可大幅削减材料制备的成本;基于金属离子的替换和亲水性官能团的引入修饰,可增强材料水吸附性能并实现精准调控,以适应多种不同环境下的性能要求. 在材料结构的选择上,可结合同类型多孔材料的研究,如共价有机框架(covalent organic framework,COF)[106]和多孔聚丙烯酸钠-石墨烯框架(porous sodium polyacrylate graphene framework,PGF)[107],补充完善MOFs材料的合成理论,为制备高性能吸水材料奠定基础.

    (2)复合材料制备. 单一MOFs因其典型的S型吸附曲线,导致材料吸附性能只能在较窄湿度范围内实现高效取水(如Yaghi教授研制的MOF-801(Zr)、MOF-303(Al)). 因此,将MOFs材料与吸湿盐、水凝胶等材料的进行复合(如王如竹教授研发的LiCl@MIL-101(Cr)、Yilmaz教授研发的PC-MOF),可有效解决单一材料的应用局限. 在水热法制备吸湿盐复合材料过程中,往往需要加入分散剂(如LiOH)以促进有机配体的质子化,可考虑减少常规的固液分离过程直接获取沉淀物,将MOFs材料与吸湿盐(如LiCl)进行复合. 同时,将粉末状MOFs材料与水凝胶材料等基质材料进行复合,实现整体式材料的合成以有效避免粉末材料负载的问题,便于材料的推广应用.

    (3)装置运行过程的优化设计. 在装置设计过程,新能源技术(如太阳能)和表面冷却技术的充分结合是一种实现空气取水装置连续自动化高效供水的有效方法. 进一步研究表明,冷凝表面的粗糙化效应能够快速收集直径小于20 mm水滴 [108-109],在高温干燥条件下表面温度的降低可有效提高MOFs吸附材料周边的相对湿度,相应的降低材料的水吸附拐点,促进水分子吸附. 同时,在装置运行过程中,实现空气过滤、强化传质传热过程,这可有效减少能量损耗,提升取水装置的工作效率[110]. 现今空气取水技术的核心问题仍在于解决环境湿度与材料性能之间的缺口,大部分研究集中于调控材料的性能以匹配环境湿度已达到空气取水目的;可基于MOFs材料制备水分子气体选择透过膜,据此搭建多级膜组件装置以富集环境水分子,完成低湿度气流向高湿度气流的转变,将膜组件装置的出口高湿气流作为工作气体以匹配合适的取水方法和材料选择,如:直接通过冷凝装置收集液态水、基于高吸湿性材料(如水凝胶)自动获取水资源.

    (4)水质检测. 吸附式空气取水装置在取水过程中,水分子伴随着解吸和冷凝过程会不可避免粘附部分杂质并影响到收集水的水质. 因此,空气取水的水质检测尤其是毒理性问题不容忽视. 可根据用水的目的不同,例如,生活饮用水、工业生产、土地灌溉、沙漠绿化等,因地制宜地选择不同类型的空气取水装置.

    以吸附式空气取水装置为代表的空气取水技术可在全球范围内无区域限制地补给水资源,可为内陆少雨地区、灾后重建地区和部队野外训练等提供液态水资源,保障生活饮用水、土地灌溉及沙漠绿化用水的需求. MOFs材料是吸附式空气取水领域的主要研究方向之一,通过金属离子、有机配体的选择,基于不同的制备方法可得到性能各异的MOFs材料. 基于材料的框架结构,MOFs与吸湿盐等材料的高效复合可得到高性能的复合材料;可据此将其设计成单循环、多批次或连续循环工作等多种取水装置. 为推进吸附式空气取水技术的产业化应用,绿色无害化的金属源和有机配体的选择、简易可控制备方法的尝试、MOFs材料和传统吸湿盐(如氯化钙CaCl2)及新型聚合物(如水凝胶材料)的复合、空气取水装置的集成化设计等是未来的主要发展趋势. 综上所述,面向空气取水的MOFs材料研究目前仍处于探索阶段,仍需要持续的深入研究,以推进MOF吸附取水装置的实际应用.

  • 图 1  MB/O3体系的实验装置图

    Figure 1.  Experimental setup diagram of MB/O3 system

    图 2  NOF的优化结构

    Figure 2.  Optimized structure of NOF

    图 3  NOF的ALIE极小值点分布

    Figure 3.  Distribution of ALIE minima of NOF

    图 4  NOF分子的CFF着色图

    Figure 4.  CFF coloring of NOF molecule

    图 5  MB/O3体系中NOF的降解途径

    Figure 5.  Degradation pathway of NOF in MB/O3 system

    图 6  在MB/O3体系下NOF及降解产物的毒性

    Figure 6.  Toxicity of NOF and degradation products in MB/O3 system

    表 1  梯度洗脱程序

    Table 1.  Gradient elution program

    时间/minTimeA相比例/%Proportion of mobile phase AB相比例/%Proportion of mobile phase B
    0955
    1955
    8595
    13595
    14955
    15955
    时间/minTimeA相比例/%Proportion of mobile phase AB相比例/%Proportion of mobile phase B
    0955
    1955
    8595
    13595
    14955
    15955
    下载: 导出CSV

    表 2  NOF不同位点键级分布

    Table 2.  Bond level distribution of different sites of NOF

    键BondLBO键级Laplacian bond order键BondLBO键级Laplacian bond order键BondLBO键级Laplacian bond order
    1C-19F0.1206C-15N0.9572C-3C1.317
    21C-23O0.54716C-17C0.9953C-4C1.332
    24N-25C0.54813C-14C0.9981C-6C1.344
    14C-15N0.57025C-26C1.0045C-6C1.350
    15N-16C0.61621C-22O1.0294C-5C1.377
    13C-18N0.66310C-24N1.11910C-11C1.460
    17C-18N0.66511C-21C1.1431C-2C1.601
    4C-24N0.8403C-8C1.151
    8C-20O0.8938C-11C1.193
    键BondLBO键级Laplacian bond order键BondLBO键级Laplacian bond order键BondLBO键级Laplacian bond order
    1C-19F0.1206C-15N0.9572C-3C1.317
    21C-23O0.54716C-17C0.9953C-4C1.332
    24N-25C0.54813C-14C0.9981C-6C1.344
    14C-15N0.57025C-26C1.0045C-6C1.350
    15N-16C0.61621C-22O1.0294C-5C1.377
    13C-18N0.66310C-24N1.11910C-11C1.460
    17C-18N0.66511C-21C1.1431C-2C1.601
    4C-24N0.8403C-8C1.151
    8C-20O0.8938C-11C1.193
    下载: 导出CSV

    表 3  NOF不同位点的CFF

    Table 3.  CFF at different sites of NOF

    原子Atom 序号Serial number qN qN+1) qN−1) f f+ f0
    C 1 0.074 0.021 0.131 0.056 0.053 0.055
    C 2 −0.054 −0.110 −0.018 0.036 0.056 0.046
    C 3 −0.048 −0.071 0.026 0.073 0.024 0.048
    C 4 0.046 0.024 0.075 0.029 0.023 0.026
    C 5 −0.074 −0.125 0.001 0.075 0.051 0.063
    C 6 0.047 0.015 0.087 0.040 0.032 0.036
    C 8 0.091 0.029 0.115 0.024 0.061 0.043
    C 10 0.061 −0.084 0.079 0.019 0.144 0.081
    C 11 −0.075 −0.126 −0.059 0.016 0.052 0.034
    C 13 −0.009 −0.011 0.006 0.015 0.003 0.009
    C 14 −0.001 −0.006 0.026 0.028 0.004 0.016
    N 15 −0.045 −0.060 0.103 0.148 0.015 0.081
    C 16 0.002 −0.003 0.029 0.027 0.005 0.016
    C 17 −0.009 −0.011 0.002 0.010 0.002 0.006
    N 18 −0.183 −0.185 −0.175 0.008 0.002 0.005
    F 19 −0.105 −0.130 −0.071 0.033 0.025 0.029
    O 20 −0.396 −0.468 −0.349 0.047 0.072 0.060
    C 21 0.195 0.157 0.203 0.008 0.039 0.023
    O 22 −0.352 −0.392 −0.341 0.012 0.040 0.026
    O 23 −0.186 −0.211 −0.180 0.006 0.025 0.015
    N 24 0.028 −0.021 0.035 0.007 0.049 0.028
    C 25 0.028 0.016 0.033 0.004 0.013 0.009
    C 26 −0.070 −0.078 −0.067 0.003 0.008 0.005
    原子Atom 序号Serial number qN qN+1) qN−1) f f+ f0
    C 1 0.074 0.021 0.131 0.056 0.053 0.055
    C 2 −0.054 −0.110 −0.018 0.036 0.056 0.046
    C 3 −0.048 −0.071 0.026 0.073 0.024 0.048
    C 4 0.046 0.024 0.075 0.029 0.023 0.026
    C 5 −0.074 −0.125 0.001 0.075 0.051 0.063
    C 6 0.047 0.015 0.087 0.040 0.032 0.036
    C 8 0.091 0.029 0.115 0.024 0.061 0.043
    C 10 0.061 −0.084 0.079 0.019 0.144 0.081
    C 11 −0.075 −0.126 −0.059 0.016 0.052 0.034
    C 13 −0.009 −0.011 0.006 0.015 0.003 0.009
    C 14 −0.001 −0.006 0.026 0.028 0.004 0.016
    N 15 −0.045 −0.060 0.103 0.148 0.015 0.081
    C 16 0.002 −0.003 0.029 0.027 0.005 0.016
    C 17 −0.009 −0.011 0.002 0.010 0.002 0.006
    N 18 −0.183 −0.185 −0.175 0.008 0.002 0.005
    F 19 −0.105 −0.130 −0.071 0.033 0.025 0.029
    O 20 −0.396 −0.468 −0.349 0.047 0.072 0.060
    C 21 0.195 0.157 0.203 0.008 0.039 0.023
    O 22 −0.352 −0.392 −0.341 0.012 0.040 0.026
    O 23 −0.186 −0.211 −0.180 0.006 0.025 0.015
    N 24 0.028 −0.021 0.035 0.007 0.049 0.028
    C 25 0.028 0.016 0.033 0.004 0.013 0.009
    C 26 −0.070 −0.078 −0.067 0.003 0.008 0.005
    下载: 导出CSV

    表 4  MB/O3体系降解NOF过程中主要中间产物

    Table 4.  The main intermediate products in the degradation of NOF in MB/O3 system

    产物Product 保留时间/mintR 质核比Mass (m/z) 峰强度Peak intensity 分子式Molecular formula 结构式Structural formula
    NOF 4.33 320.14 6.85×106 C16H18FN3O3
    D1 4.33 302.13 5.75×105 C16H19N3O3
    D2 7.55 318.30 3.11×106 C16H19N3O4
    D3 7.55 274.28 6.13×106 C15H19N3O2
    D4 4.80 350.12 1.89×105 C16H16FN3O5
    D5 4.33 322.25 1.98×105 C15H16FN3O4
    D6 4.06 294.13 2.04×105 C14H16FN3O3
    D7 5.34 279.08 6.21×105 C13H11FN2O4
    D8 7.06 251.18 1.46×105 C12H11FN2O3
    D9 0.62 190.12 8.47×104 C11H12N2O
    D10 4.06 336.14 1.19×106 C16H18FN3O4
    D11 10.64 338.34 5.33×106 C16H20FN3O4
    D12 3.25 296.14 3.73×106 C14H18FN3O3
    D13 3.25 268.15 3.59×105 C13H18FN3O2
    D14 2.98 240.05 1.25×105 C12H18FN3O
    D15 0.62 172.11 6.58×104 C8H10FNO2
    D16 3.72 305.16 1.55×105 C16H20FN3O2
    D17 3.25 296.14 3.73×106 C14H18FN3O3
    D18 2.98 269.15 6.64×104 C12H16FN3O3
    D19 4.93 227.18 1.90×105 C10H11FN2O3
    产物Product 保留时间/mintR 质核比Mass (m/z) 峰强度Peak intensity 分子式Molecular formula 结构式Structural formula
    NOF 4.33 320.14 6.85×106 C16H18FN3O3
    D1 4.33 302.13 5.75×105 C16H19N3O3
    D2 7.55 318.30 3.11×106 C16H19N3O4
    D3 7.55 274.28 6.13×106 C15H19N3O2
    D4 4.80 350.12 1.89×105 C16H16FN3O5
    D5 4.33 322.25 1.98×105 C15H16FN3O4
    D6 4.06 294.13 2.04×105 C14H16FN3O3
    D7 5.34 279.08 6.21×105 C13H11FN2O4
    D8 7.06 251.18 1.46×105 C12H11FN2O3
    D9 0.62 190.12 8.47×104 C11H12N2O
    D10 4.06 336.14 1.19×106 C16H18FN3O4
    D11 10.64 338.34 5.33×106 C16H20FN3O4
    D12 3.25 296.14 3.73×106 C14H18FN3O3
    D13 3.25 268.15 3.59×105 C13H18FN3O2
    D14 2.98 240.05 1.25×105 C12H18FN3O
    D15 0.62 172.11 6.58×104 C8H10FNO2
    D16 3.72 305.16 1.55×105 C16H20FN3O2
    D17 3.25 296.14 3.73×106 C14H18FN3O3
    D18 2.98 269.15 6.64×104 C12H16FN3O3
    D19 4.93 227.18 1.90×105 C10H11FN2O3
    下载: 导出CSV
  • [1] KOVALAKOVA P, CIZMAS L, MCDONALD T J, et al. Occurrence and toxicity of antibiotics in the aquatic environment: A review[J]. Chemosphere, 2020, 251: 126351. doi: 10.1016/j.chemosphere.2020.126351
    [2] LI S M, LI J, LI Z, et al. Toxic effects of norfloxacin in soil on fed and unfed Folsomia candida (Isotomidae: Collembola) and on gut and soil microbiota[J]. The Science of the Total Environment, 2021, 788: 147793. doi: 10.1016/j.scitotenv.2021.147793
    [3] ZHANG Q F, ZHANG H, ZHANG Q X, et al. Degradation of norfloxacin in aqueous solution by atmospheric-pressure non-thermal plasma: Mechanism and degradation pathways[J]. Chemosphere, 2018, 210: 433-439. doi: 10.1016/j.chemosphere.2018.07.035
    [4] ZHU F, WU Y Y, LIANG Y K, et al. Degradation mechanism of norfloxacin in water using persulfate activated by BC@nZVI/Ni[J]. Chemical Engineering Journal, 2020, 389: 124276. doi: 10.1016/j.cej.2020.124276
    [5] GHOSH S, KAR S, PAL T, et al. Sunlight-driven photocatalytic degradation of Norfloxacin antibiotic in wastewater by ZnSe microsphere functionalized RGO composite[J]. Sustainable Chemistry for the Environment, 2023, 4: 100038. doi: 10.1016/j.scenv.2023.100038
    [6] WANG X J, HU S S, MAO H T, et al. Facile fabrication of AgVO3/rGO/BiVO4 hetero junction for efficient degradation and detoxification of norfloxacin[J]. Environmental Research, 2023, 227: 115623. doi: 10.1016/j.envres.2023.115623
    [7] FANG C, XU H B, WANG S H, et al. Simultaneous removal of norfloxacin and chloramphenicol using cold atmospheric plasma jet (CAPJ): Enhanced performance, synergistic effect, plasma-activated water (PAW) contribution, mechanism and toxicity evaluation[J]. Journal of Hazardous Materials, 2023, 452: 131306. doi: 10.1016/j.jhazmat.2023.131306
    [8] KIM S H, PARK S Y, KIM G E, et al. Effect of pH and temperature on the biodegradation of oxytetracycline, streptomycin, and validamycin A in soil[J]. Applied Biological Chemistry, 2023, 66(1): 63. doi: 10.1186/s13765-023-00822-1
    [9] SUMESH C K, PAREKH K. Nanocatalytic physicochemical adsorption and degradation of organic dyes[J]. Pramana, 2019, 92(6): 87. doi: 10.1007/s12043-019-1760-0
    [10] WANG C K, LIN C Y, LIAO G Y. Degradation of antibiotic tetracycline by ultrafine-bubble ozonation process[J]. Journal of Water Process Engineering, 2020, 37: 101463. doi: 10.1016/j.jwpe.2020.101463
    [11] WANG P, ZHANG H X, WU Z H, et al. A data-based review on norfloxacin degradation by persulfate-based advanced oxidation processes: Systematic evaluation and mechanisms[J]. Chinese Chemical Letters, 2023, 34(12): 108722. doi: 10.1016/j.cclet.2023.108722
    [12] YUAN Q B, QU S Y, LI R, et al. Degradation of antibiotics by electrochemical advanced oxidation processes (EAOPs): Performance, mechanisms, and perspectives[J]. The Science of the Total Environment, 2023, 856(Pt 2): 159092.
    [13] ANJALI R, SHANTHAKUMAR S. Insights on the current status of occurrence and removal of antibiotics in wastewater by advanced oxidation processes[J]. Journal of Environmental Management, 2019, 246: 51-62.
    [14] LI S, WU Y N, ZHENG H S, et al. Antibiotics degradation by advanced oxidation process (AOPs): Recent advances in ecotoxicity and antibiotic-resistance genes induction of degradation products[J]. Chemosphere, 2023, 311(Pt 2): 136977.
    [15] ZHANG J, LV S Y, YU Q, et al. Degradation of sulfamethoxazole in microbubble ozonation process: Performance, reaction mechanism and toxicity assessment[J]. Separation and Purification Technology, 2023, 311: 123262 doi: 10.1016/j.seppur.2023.123262
    [16] LIU T, ZHANG B, LI W Q, et al. The catalytic oxidation process of atrazine by ozone microbubbles: Bubble formation, ozone mass transfer and hydroxyl radical generation[J]. Chemosphere, 2023, 325: 138361. doi: 10.1016/j.chemosphere.2023.138361
    [17] VERINDA S B, MUNIROH M, YULIANTO E, et al. Degradation of ciprofloxacin in aqueous solution using ozone microbubbles: Spectroscopic, kinetics, and antibacterial analysis[J]. Heliyon, 2022, 8(8): e10137. doi: 10.1016/j.heliyon.2022.e10137
    [18] 杨旭. 微气泡催化臭氧化降解制药废水过程及机制研究[D]. 石家庄: 河北科技大学, 2022.

    YANG X. Study on the process and mechanism of catalytic ozonation of microbubbles in degradation of pharmaceutical wastewater[D]. Shijiazhuang: Hebei University of Science and Technology, 2022 (in Chinese).

    [19] ZHANG J, LIU M, PANG B W, et al. Ciprofloxacin degradation in microbubble ozonation combined with electro-generated H2O2 process: Operational parameters and oxidation mechanism[J]. Separation and Purification Technology, 2023, 325: 124676 doi: 10.1016/j.seppur.2023.124676
    [20] 胡英, 刘洪来. 密度泛函理论[M]. 北京: 科学出版社, 2016.

    HU Y, LIU H L. Density functional theory[M]. Beijing: Science Press, 2016(in Chinese).

    [21] 卢志磊, 范勇杰, 陈洁洁, 等. 基于Gaussian、ECOSAR模型的紫外/次氯酸体系降解含卤阻燃剂的产物预测与毒性评估[J]. 环境化学, 2024, 43(1): 82-91. doi: 10.7524/j.issn.0254-6108.2022063002

    LU Z L, FAN Y J, CHEN J J, et al. Product prediction and toxicity assessment of halogen-containing flame retardants degraded by UV/hypochlorous acid systems based on Gaussian and ECOSAR models[J]. Environmental Chemistry, 2024, 43(1): 82-91 (in Chinese). doi: 10.7524/j.issn.0254-6108.2022063002

    [22] GUTIÉRREZ-SÁNCHEZ P, ÁLVAREZ-TORRELLAS S, LARRIBA M, et al. Influence of transition metal-based activating agent on the properties and catalytic activity of sewage sludge-derived catalysts. Insights on mechanism, DFT calculation and degradation pathways[J]. Journal of Molecular Liquids, 2023, 381: 121840. doi: 10.1016/j.molliq.2023.121840
    [23] YEOM Y, HAN J R, ZHANG X R, et al. A review on the degradation efficiency, DBP formation, and toxicity variation in the UV/chlorine treatment of micropollutants[J]. Chemical Engineering Journal, 2021, 424: 130053. doi: 10.1016/j.cej.2021.130053
    [24] 陈洁洁, 范勇杰, 杨婧, 等. 基于Guassian、Multiwfn及ECOSAR的氮磷系阻燃剂在UV/PDS体系中的降解路径及产物毒性预测研究[J]. 环境科学学报, 2022, 42(12): 39-48.

    CHEN J J, FAN Y J, YANG J, et al. Degradation pathway and product toxicity prediction of nitrogen phosphorus flame retardants in UV/PDS system based on Guassian, Multiwfn and ECOSAR[J]. Acta Scientiae Circumstantiae, 2022, 42(12): 39-48 (in Chinese).

    [25] MARENICH A V, CRAMER C J, TRUHLAR D G. Universal solvation model based on solute electron density and on a continuum model of the solvent defined by the bulk dielectric constant and atomic surface tensions[J]. The Journal of Physical Chemistry. B, 2009, 113(18): 6378-6396. doi: 10.1021/jp810292n
    [26] FRISCH M J, TRUCKS G W, SCHLEGEL H B, et al. Gaussian 16 Rev. C. 01[CP/OL].
    [27] LU T, CHEN F W. Bond order analysis based on the Laplacian of electron density in fuzzy overlap space[J]. The Journal of Physical Chemistry. A, 2013, 117(14): 3100-3108. doi: 10.1021/jp4010345
    [28] LU T, CHEN F W. Multiwfn: A multifunctional wavefunction analyzer[J]. Journal of Computational Chemistry, 2012, 33(5): 580-592. doi: 10.1002/jcc.22885
    [29] ALAM M S, LEE D U. Molecular structure, spectral (FT-IR, FT-Raman, Uv-Vis, and fluorescent) properties and quantum chemical analyses of azomethine derivative of 4-aminoantipyrine[J]. Journal of Molecular Structure, 2021, 1227: 129512. doi: 10.1016/j.molstruc.2020.129512
    [30] HUMPHREY W, DALKE A, SCHULTEN K. VMD: Visual molecular dynamics[J]. Journal of Molecular Graphics, 1996, 14(1): 33-38,27-28. doi: 10.1016/0263-7855(96)00018-5
    [31] GEERLINGS P, de PROFT F, LANGENAEKER W. Conceptual density functional theory[J]. Chemical Reviews, 2003, 103(5): 1793-1873. doi: 10.1021/cr990029p
    [32] FENG X F, LONG R X, LIU C C, et al. Flexible Z-scheme heterojunction membrane reactors for visible-light-driven antibiotic degradation and oil-water separation[J]. Chemical Engineering Journal, 2023, 471: 144447. doi: 10.1016/j.cej.2023.144447
    [33] LYKOS C, KOURKOUTA T, KONSTANTINOU I. Study on the photocatalytic degradation of metronidazole antibiotic in aqueous media with TiO2 under lab and pilot scale[J]. The Science of the Total Environment, 2023, 870: 161877. doi: 10.1016/j.scitotenv.2023.161877
    [34] 王福生, 孙玮, 张渝, 等. 过渡金属离子促进煤自燃机理的量子化学计算[J]. 煤炭学报, 2024, 49(5) : 2347-2359.

    WANG F S, SUN W, ZHANG Y, et al. Quantum chemical calculation of transition metal ions promoting coal spontaneous combustion mechanism[J]. Journal of China Coal Society, 2024, 49(5) : 2347-2359 (in Chinese).

    [35] PARLAK C, ALVER Ö. A density functional theory investigation on amantadine drug interaction with pristine and B, Al, Si, Ga, Ge doped C60 fullerenes[J]. Chemical Physics Letters, 2017, 678: 85-90. doi: 10.1016/j.cplett.2017.04.025
    [36] KURBAN M, MUZ İ. Theoretical investigation of the adsorption behaviors of fluorouracil as an anticancer drug on pristine and B-, Al-, Ga-doped C36 nanotube[J]. Journal of Molecular Liquids, 2020, 309: 113209. doi: 10.1016/j.molliq.2020.113209
    [37] 付蓉, 卢天, 陈飞武. 亲电取代反应中活性位点预测方法的比较[J]. 物理化学学报, 2014, 30(4): 628-639. doi: 10.3866/PKU.WHXB201401211

    FU R, LU T, CHEN F W. Comparing methods for predicting the reactive site of electrophilic substitution[J]. Acta Physico-Chimica Sinica, 2014, 30(4): 628-639 (in Chinese). doi: 10.3866/PKU.WHXB201401211

    [38] WU Y H, LI Y L, HE J Y, et al. Nano-hybrids of needle-like MnO2 on graphene oxide coupled with peroxymonosulfate for enhanced degradation of norfloxacin: A comparative study and probable degradation pathway[J]. Journal of Colloid and Interface Science, 2020, 562: 1-11. doi: 10.1016/j.jcis.2019.11.121
    [39] DENG J, XU M Y, CHEN Y J, et al. Highly-efficient removal of norfloxacin with nanoscale zero-valent copper activated persulfate at mild temperature[J]. Chemical Engineering Journal, 2019, 366: 491-503. doi: 10.1016/j.cej.2019.02.073
    [40] XUE L, ZHAO C X, LIU J M, et al. Microwave electrodeless UV light source combine ozone generation with photocatalytic simultaneous degradation of norfloxacin[J]. Chemical Engineering and Processing - Process Intensification, 2023, 186: 109325. doi: 10.1016/j.cep.2023.109325
    [41] XU Z M, XUE X J, HU S H, et al. Degradation effect and mechanism of gas-liquid phase dielectric barrier discharge on norfloxacin combined with H2O2 or Fe2+[J]. Separation and Purification Technology, 2020, 230: 115862. doi: 10.1016/j.seppur.2019.115862
    [42] ZHAO Y P, SUN Q, ZHANG J, et al. Construction of Fe/N/C nano-clusters anchored on porous diatomite for efficient removal of norfloxacin via the adsorption-PMS activation[J]. Separation and Purification Technology, 2023, 310: 123127. doi: 10.1016/j.seppur.2023.123127
    [43] ZHOU Y, HE J, LU J, et al. Enhanced removal of bisphenol A by cyclodextrin in photocatalytic systems: Degradation intermediates and toxicity evaluation[J]. Chinese Chemical Letters, 2020, 31(10): 2623-2626. doi: 10.1016/j.cclet.2020.02.008
    [44] LIU H Y, CHEN H J, ADDISON F, et al. Insights into electrocatalytic oxidation of aqueous ampicillin: Degradation mechanism and potential toxicity from intermediates[J]. Journal of Environmental Chemical Engineering, 2022, 10(6): 108673. doi: 10.1016/j.jece.2022.108673
    [45] YANG J X, LV G C, LI T T, et al. Theoretical insight into the degradation of diclofenac by hydroxyl and sulfate radicals in aqueous-phase: Mechanisms, kinetics and eco-toxicity[J]. Journal of Environmental Chemical Engineering, 2022, 10(5): 108311. doi: 10.1016/j.jece.2022.108311
    [46] HE B, SONG L X, ZHAO Z X, et al. CuFe2O4/CuO magnetic nano-composite activates PMS to remove ciprofloxacin: Ecotoxicity and DFT calculation[J]. Chemical Engineering Journal, 2022, 446: 137183. doi: 10.1016/j.cej.2022.137183
  • 加载中
图( 6) 表( 4)
计量
  • 文章访问数:  874
  • HTML全文浏览数:  874
  • PDF下载数:  20
  • 施引文献:  0
出版历程
  • 收稿日期:  2023-12-21
  • 录用日期:  2024-03-26
  • 刊出日期:  2025-01-27
张静, 于倩, 丁亮亮, 刘春, 闫晨旭, 牛建瑞, 张瑞娜. 基于密度泛函理论的微气泡臭氧化诺氟沙星降解途径与毒性评估[J]. 环境化学, 2025, 44(1): 288-298. doi: 10.7524/j.issn.0254-6108.2023122102
引用本文: 张静, 于倩, 丁亮亮, 刘春, 闫晨旭, 牛建瑞, 张瑞娜. 基于密度泛函理论的微气泡臭氧化诺氟沙星降解途径与毒性评估[J]. 环境化学, 2025, 44(1): 288-298. doi: 10.7524/j.issn.0254-6108.2023122102
ZHANG Jing, YU Qian, DING Liangliang, LIU Chun, YAN Chenxu, NIU Jianrui, ZHANG Ruina. Degradation pathway and toxicity assessment of microbubble ozonation of Norfloxacin based on Density Functional Theory[J]. Environmental Chemistry, 2025, 44(1): 288-298. doi: 10.7524/j.issn.0254-6108.2023122102
Citation: ZHANG Jing, YU Qian, DING Liangliang, LIU Chun, YAN Chenxu, NIU Jianrui, ZHANG Ruina. Degradation pathway and toxicity assessment of microbubble ozonation of Norfloxacin based on Density Functional Theory[J]. Environmental Chemistry, 2025, 44(1): 288-298. doi: 10.7524/j.issn.0254-6108.2023122102

基于密度泛函理论的微气泡臭氧化诺氟沙星降解途径与毒性评估

    通讯作者: E-mail:liuchun@hebust.edu.cn
  • 河北科技大学环境科学与工程学院,河北省污染防治生物技术重点实验室,石家庄,050018
基金项目:
国家重点研发计划项目(2022YFE0119000)资助.

摘要: 诺氟沙星(Norfloxacin,NOF)为第三代喹诺酮类抗菌药被广泛应用,易在环境中积累,威胁人类健康和环境安全. 本研究采用Gaussian与Toxicity Estimation Software Tool(T.E.S.T.)软件预测在微气泡臭氧化(MB/O3)体系中NOF的降解途径与降解产物的生态毒性. 结果表明,MB/O3体系中的羟基自由基等活性氧物质易对NOF分子结构上键级较小、平均局部离子化能较低以及简缩Fukui函数较高的反应位点进行攻击;液相色谱-串联质谱法(LC-MS/MS)分析结果结合模型计算推测,MB/O3处理NOF降解途径包括:脱氟反应、哌嗪环转化、羟基取代反应和喹诺酮环氧化;利用T.E.S.T.软件评估NOF及降解产物的毒性,发现NOF对口服大鼠、胖头鱼、水蚤均存在负面影响,具有发育毒性,无致突变性,且降解过程中可能会产生比NOF毒性大的产物. 本研究采用模型预测NOF降解途径与生态毒性的方法便捷有效,可为MB/O3体系降解复杂有机分子的降解途径和生态毒性研究提供理论依据.

English Abstract

  • 诺氟沙星(NOF)又名氟哌酸,是第三代喹诺酮类抗菌药,是一种难降解性有机物[1]. NOF具有抗菌谱广和作用性强的特点,已广泛用于人类和兽医学,以预防和治疗细菌引起的感染[2]. 而NOF并不被人类和牲畜完全代谢,因此,NOF可能会通过尿液和粪便的方式被排放到城市污水中[3],容易在环境中累积[45],导致人类生活环境风险的增加. 有研究表明,NOF对植物和水生生物具有毒性[67],会对水生环境造成不利的生态影响. 显然,为了保证人类及生态环境健康发展,对NOF进行降解减毒研究具有重要意义. 废水中喹诺酮类抗生素的传统去除方法可分为生物降解[8]、物理吸附[9]、和化学氧化,但存在降解时间长、矿化度不高、易造成二次污染等缺点. 近年来,臭氧[10]、过硫酸盐[11]、电化学[12]等高级氧化过程(AOPs)氧化技术引起了广泛的关注. AOPs已被证明可有效降解抗生素药物[1314]. 其中,臭氧氧化具有无二次污染、反应效率高和成本低的优势. 臭氧在水溶液中可由链式反应[15]分解产生活性氧物质(ROS),如·OH、1O2、·O2-等,·OH具有强氧化性(E0 = 2.80 eV)和非选择性的特点,可以与溶液中大多数抗生素药物发生快速反应. 但单独臭氧氧化在利用率和矿化能力上存在缺点,限制了其应用. 微气泡(MB,平均直径<50 μm)技术的研究为单独臭氧化过程的不足提供了解决办法. MB不同于传统气泡(CB),MB表面张力会引起较高的内部压力,使得臭氧溶解在水中,并发生MB的收缩和塌陷[16]. MB的界面面积更大,导致其在水中上升速度变慢. 因此,这增加了臭氧从气相转移到液相的效率,进而增加抗生素与臭氧微气泡发生氧化反应的可能性[17]. 微气泡臭氧化(MB/O3)体系利用MB与臭氧耦合产生ROS,依据已有报道,该系统主要ROS为·OH[15,1819]. 密度泛函理论(DFT)是一种研究物质结构统计性的理论,既可以用于研究微观结构,也可以用于研究宏观结构,是提供原子、分子以及晶体和表面的结构,并得到相应结构特性的基本工具[20]. 近年来,基于DFT应用Gaussian模型预测优先反应位点并推测降解产物的研究者逐渐增多,卢志磊等[21]和Pablo等[22]使用Gaussian模型预测有机物活性位点推测降解产物路径. 但降解产物毒性鉴定需要昂贵或复杂的仪器进行分析[23].

    因此,本研究基于密度泛函理论(DFT)方法,采用Gaussian模型和Multiwfn程序计算NOF的各拉普拉斯键级(Laplacian bond order, LBO)、平均局部离子化能(average local ionization energy, ALIE)和简缩Fukui函数确定MB/O3体系中NOF优先发生反应的活性位点,与MB/O3工艺降解NOF实验所得LC-MS/MS的结果相结合,从而预测可能的反应路径. 最后,通过定量构效关系(QSAR)考察降解中间产物的毒性[24],利用T.E.S.T.模型评估NOF及其中间体的生态毒性,为NOF在MB/O3体系下的降解途径与中间产物毒性预测提供新思路.

    • 诺氟沙星(NOF,≥98%)购自麦克林生化科技股份有限公司(上海,中国),靛蓝二磺酸钠(90%)从易恩化学技术有限公司购买(上海,中国). 碘化钾和硫代硫酸钠购于永大试剂有限公司(天津,中国). 去离子水由Milli-Q超纯水机制备. 所有化学试剂均为分析纯.

    • 实验装置如图1所示,在臭氧浓度12 mg·min−1、气流量0.3 L·min−1、初始NOF浓度为50 mg·L−1、有效体积为10 L的实验条件下,进行MB/O3降解NOF废水实验. 实验时间90 min,前30 分钟每隔5 min取适量NOF水样,后30 min每隔10 min取样1次,水样经有机滤膜(0.45 μm,有机尼龙66)过滤后,进行LC-MS/MS测定.

    • 本研究采用液相色谱(1290 UPLC,安捷伦)与质谱(Q-TOF 6550,安捷伦)联用测定中间产物的质核比进而与模拟结果相对比,确定降解中间产物.

      色谱条件:Waters BEH C18(2.1 mm×100 mm 1.7 μm);流速0.3 mL·min−1;进样量5μL;流动相A:0.1%甲酸-水溶液;流动相B:甲醇溶液. 梯度洗脱程序见表1.

      质谱条件:ESI+模式电压4000V;一级扫描50—600 m/z;二级扫描起始离子50 m/z;鞘气温度:350 ℃;鞘气流量:12 L·min−1.

    • 借助Gauss View 6.0对NOF分子进行几何结构构建,以DFT为基础,溶剂化效应选取Solvation Model Based on Density (SMD)溶剂模型[25] (solvent=water),使用B3LYP-D3泛函,选择6-311 G (d,p)基组,采用Gaussian 16W[26]进行NOF分子振动频率和几何结构优化计算.

    • 在对NOF分子的化学键键级计算时,选用LBO键级. LBO比Mayer bond order和natural localized molecular orbital bond order更具有区分键强度的能力,这使得LBO可以更准确地描述和预测化学键的性质和行为. LBO与键极性和键解离能之间存在直接的相关性,它已被证明具有合理性和实用性[27]. 因此,使用Multiwfn 3.8[28]程序进行LBO键级计算.

    • ALIE被定义为由轨道密度加权的轨道能量的总和,用于理解分子的化学反应性[29]. ALIE方法可用于分子反应位点的定量分析,其中可能更倾向于亲电攻击. ALIE极值点分布图可以使用Multiwfn程序可视化,利用Visual Molecular Dynamics (VMD 1.9.3)[30]程序对ALIE极值点分布图进行渲染. ALIE的极小值表示电子结合更松散且容易与亲电试剂反应的区域.

    • Fukui函数是指分子中电子密度函数的变化[31]. 为了便于在定量层面上比较各个位点的福井函数,定义了简缩Fukui函数(CFF),CFF是基于Hirshfeld电荷(q)来计算的,是将Fukui函数收缩到各个原子上,使每个原子具有一个确定的数值. 使用Multiwfn采用了Hirshfeld电荷计算NOF分子的CFF. 结合VMD软件绘制了NOF分子的CFF原子着色图来预测NOF分子亲电或自由基反应活性位点.

    • 采用QSAR计算NOF及其降解产物的生态毒性. 通过美国环境保护署开发的T.E.S.T.软件(版本5.1.2)评估氧化产物的急性毒性、生物累积因子、发育毒性和致突变性. 将CAS登录号或物质名称或分子结构输入T.E.S.T.程序中,随后得出LC50(96 h/48 h内对半数暴露于水中的胖头鱼/水蚤致死的浓度)、LD50(大鼠口服该化学品后,半数致死的数量)、生物累积因子、发育毒性和致突变性[32]. 所采用的方法是共识方法,因为此方法在外部验证过程中获得了最佳的预测结果[33].

    • 借助Gauss View 6.0对NOF分子进行几何结构构建,NOF分子几何结构优化如图2所示. 键级又称为键序,属于对原子之间化学键的定量分析,适用于有机体系,可以被用来预测分子的断键强度[34]. LBO与键解离能有直接关系,可以用来评估两个原子之间的化学键强弱[27,35],且LBO与原子之间的键长呈正相关[36]. LBO计算所得数值如表2所示,1C-19F的键级值最小,为0.120. 其次为21C-23O和24N-25C,分别为0.547和0.548. 苯环上的碳碳双键的键级在1.31—1.38的范围内,键级值较高. 因此,C(1)和F(19)之间的键最容易断裂,发生脱氟反应. 其次为O(21)-O(23)和N(24)-C(25)之间的共价键容易受到影响被破坏. 苯环的共价键由于其较高的键级值不容易发生断裂.

    • ALIE是指分子失去电子变成离子的能量,所需能量越小,表示电子解离能越低,更易发生亲电和自由基反应. ALIE可弥补分子表面静电势在分析亲电反应活性位点方面的不足. 如图3所示,极小值点已经标出,以便直观地表达所有位置的情况. 蓝色代表较小的ALIE,绿色代表较大的ALIE,白色充当从蓝色到绿色的过渡颜色. 青色代表等表面上的ALIE极小值,可能为亲电攻击的有利位置. NOF分子的哌嗪环上的N(18)原子、苯环和喹诺酮环上的C(11)原子附近的蓝色区域尤为明显,分子哌嗪环上N(15)原子和喹诺酮环上3个O(20)、O(22)和O(23)原子处也存在蓝色区域. N原子附近存在整体ALIE的最小值是7.39 eV,表明N原子处最容易发生亲电反应或者自由基反应. C(10)和C(11)之间的双键、N(15)和喹诺酮环上3个O(20)、O(22)和O(23)原子附近ALIE数值较小,亲电反应也可能在这些地方发生. 这是由于这些原子存在孤对电子从而容易发生亲电反应. 而苯环为π-π共轭体系,比其他键反应活性弱,不易于发生亲电反应.

    • CFF中f-表示为NOF分子体系中各原子可以作为亲电反应位点的程度,f0体现的是自由基攻击NOF分子上的各原子的程度,数值越正越容易发生相应类型的反应[37]. 如图4所示为根据CFF的数值大小对NOF分子上的原子进行着色,红色为正值,蓝色为负值. 各原子的颜色越红表明数值越大,数值越小颜色越蓝. NOF不同位点的CFF计算结果如表2所示,根据表3描述亲电反应的数值由大到小排序为:15(N) > 5(C) > 3(C) >1(C) > 20(O) >6(C) >2(C) > 19(F) > 4(C) > 14(C)>16(C) > 8(C) > 10(C ) > 11(C)> 13(C) > 22(O) > 17(C) > 18(N) > 21(C) > 24(N) > 23(O) > 27(H) > 25(C) > 26(C). 描述自由基反应的数值由大到小排序为:10(C) > 15(N) > 5(C) > 20(O) > 1(C) > 3(C) > 2(C) > 8(C) > 6(C) > 11(C) > 19(F) > 24(N) > 4(C) > 22(O) > 21(C) > 16(C) > 14(C) > 23(O) > 13(C) > 25(C) > 17(C) > 26(C) > 18(N).

      综上所述,N15(f-=0.148)、C5(f-=0.075)、O20(f-=0.047)、F19(f-=0.033)容易受到亲电试剂的攻击,C10(f0 = 0.081)、N15(f0 = 0.081)、C5(f0 = 0.063)和O20(f0 = 0.060)容易受到自由基攻击.

    • 在MB/O3体系中,·OH起主要降解作用. 采用LC-MS/MS检测MB/O3过程中的NOF降解中间产物和粗略的定量分析,在MB/O3降解实验中检测到19种中间体如表4所示. DFT预测结果结合LC-MS/MS检测的结果,可获得MB/O3降解NOF主要通过3种反应类型:(1)脱氟反应;(2)哌嗪环转化;(3)喹诺酮环氧化. MB/O3降解NOF反应途径如图5所示.

      途径Ⅰ:由于NOF分子的1(C)-19(F)键的不稳定性,该键发生断裂,导致脱氟反应的发生,并生成D1,与上述计算所得出1(C)-19(F)的LBO键级是最小的一致,CFF预测1(C)-19(F)之间的共价键容易断裂. 随后在1(C)上发生羟基化反应生成羟基化产物D2,进一步通过脱羧反应生成D3. Wu等[38]报道了使用MnO2/GO纳米材料活化过一硫酸盐技术降解诺氟沙星过程中有类似的途径. 在本研究中降解30 min和80 min时,通过总离子色谱仪(TIC)模式观察,发现D2的峰有明显的增加趋势. 在80 min时,观察到D3的峰最为显著,这可能进一步证明了途径Ⅰ为主要途径.

      途径Ⅱ:哌嗪环部分被·OH氧化生成中间产物D4. 根据计算CFF的亲电反应,13(C)和14(C)原子处具有较高的f-,易亲电试剂所攻击. 然后哌嗪环被氧化发生裂解和结构重排形成含有醛基产物D5,醛基被氧化失去后形成D6. 消除反应发生在N原子上,16(C)-17(C)被氧化形成带醛基的产物D7. 哌嗪环被活性氧氧化消失,剩下一个氨基生成D8. 然后发生脱羧和脱氟反应生成D9. 在Deng等[39]和Li等[40]的研究结果中,发现了与此途径类似的中间产物.

      途径Ⅲ:10(C)处CFF的f0= 0.0814,数值最大,也最有可能在此处发生自由基反应. ·OH通过抽氢反应和自由基结合形成产物D10,Xu等[41]报道·OH容易攻击NOF喹诺酮环中羧酸基的C=C键. D11是C=C被ROS氧化喹诺酮环开环的产物,喹诺酮环再进行氧化失去羰基生成D12和D13,然后脱去羧基形成D14和CO2,哌嗪环被氧化完全后再发生羟基加成形成D15. 通过30 min时的TIC模式观察,D11的丰度仅次于D3,途径Ⅲ中的各中间产物(除了D15)均在30 min时被观察到. 在30 min降解过程中,途径Ⅲ中的喹诺酮环氧化是MB/O3过程中仅次于途径Ⅰ的降解途径.

      途径Ⅳ:由ALIE方法得出22号O原子可以被ROS所攻击发生亲电反应. NOF发生脱羧反应生成D16[42]. D17是在D16的喹诺酮环上面的C=C处发生氧化反应形成的. D17的哌嗪环上发生氧化开环反应形成D18. D19是哌嗪环全部被氧化完全得到的结果. 通过观察降解实验发生30 min和80 min的TIC,可以发现在所有途径中,途径Ⅳ的产物离子丰度都是最低的,甚至远远低于途径Ⅰ. 相比之下,途径Ⅰ和Ⅲ的产物离子丰度较高,说明在MB/O3过程中它们是主导的反应途径. 因此,根据这些观察结果,可以得出途径Ⅳ是MB/O3过程中最次要的降解途径.

    • 基于QSAR,使用T.E.S.T毒性评估工具预测MB/O3体系中的NOF和其中间产物的生态毒性. QSAR模型根据分子理化性质参数或化学结构参数与其毒理活性之间的定量关系预测相应的毒性数据[43]. 通过T.E.S.T.软件可计算出导致50%的口服大鼠(OR)死亡的化学物质剂量LD50、导致50%的胖头鱼(FM)在96 h后死亡的化学物质浓度LC50、导致50%的水蚤(DM)在48 h后死亡的化学物质浓度LC50、生物富集因子(BCF)、致突变性和发育毒性. 根据全球化学品统一分类和标签制度(GHS)和上述计算所得值,进而评估NOF及降解产物的毒性.

      在MB/O3体系下NOF及降解产物的毒性如图6(a)所示. 可以看到,一些OR的LD50没有预测数值,这可能是缺乏相应的实验数据[44]. 对于OR,NOF被认定具有轻度危害. 而在MB/O3过程中产生的降解中间体的LD50均小于母体化合物(NOF),表明在降解过程中产生了毒性略强的中间体,尤其是D14对于口服大鼠被认定具有中度危害. 导致FM和DM在96 h和48 h后死亡的化学物质浓度小于1 mg·L−1、1—10 mg·L−1、10—100 mg·L−1和大于100 mg·L−1,分别被定义为“剧毒”、“有毒”、“有害”和“无害”. NOF对FM表现为有毒作用,而对DM表现为有害作用. 对于FM,大多数臭氧氧化中间体被定义为“有毒”化合物,对水生生物构成了威胁,这可能是因为这些产物中含有脂肪族胺结构或者酰胺结构. 但这些产物的毒性小于母体化合物,D19对于FM是没有毒害作用的. 对于DM,D6、D8、D9、D14和D15被定义为“有毒”化合物,其毒性大于母体化合物. D11、D17和D18被定义为“无害”化合物,因为它们的急性毒性均高于100 mg·L−1. 此外,BCF可以用来衡量生物从环境之中吸收并积累某种物质的程度,除D3、D9、D13、D14、D15和D16外,其余产物的生物累积因子值均低于母体化合物,表明NOF的降解过程伴随着累积毒性.

      发育毒性是受到环境因素(如药物、环境中的有毒化学物质)的影响导致生物体不能够正常生长、分化和发育[45]. 如图6(b)所示,NOF的发育毒性为1.09,高于毒性阈值0.5. 表明NOF是一种发育毒性化合物. 这与NOF本身药物使用事的注意事项一致,18岁以下儿童患者不能服用,否则会延缓未成年人的骨骼形成,会影响到发育. 尽管NOF的降解中间产物的发育毒性有所降低,但仍处于发育毒性水平(D15除外). 此外,根据致突变性结果,除D14、D15、D18和D19外,NOF及其他降解中间产物致突变性均高于0.5[46],并不具有致突变性. 随着NOF的降解,降解产物中仅少部分为诱变毒性化合物. 急性毒性,生物富集因子,致突变性和发育毒性是综合毒性的决定因素. 一些降解中间产物,如D6、D8、D10、D14和D15,由于其较高的综合毒性,具有较高的环境风险,我们应更加重视降解产物的毒性.

    • (1) 采用Gaussian模型计算典型抗生素NOF的LBO键级,利用Multiwfn程序计算平均局部离子化能和简缩Fukui函数,其中C(1)和F(19)之间的键最容易断裂,发生脱氟反应;哌嗪环上的N(18)原子、苯环和喹诺酮环上的C(11)原子最易发生亲电反应和自由基反应;C5 (f 0 = 0.063)和O20 (f 0 = 0.0598)既容易受到亲电攻击又容易受到自由基攻击. 基于模型预测,结合LC-MS/MS分析推测得到MB/O3体系下NOF的降解途径主要包括:脱氟、哌嗪环转化和喹诺酮环氧化反应.

      (2) 利用T.E.S.T.软件对上述降解产物的毒性进行评估,得出NOF对大鼠、胖头鱼和水蚤有害,NOF具有发育毒性但无致突变性;NOF降解过程中会产生毒性比母体强的产物,并且降解过程中可能会伴随毒性累积.

    参考文献 (46)

返回顶部

目录

/

返回文章
返回