CN105701340A - 预测气态含硫化合物常温下在活性炭上的吸附速率常数的方法 - Google Patents

预测气态含硫化合物常温下在活性炭上的吸附速率常数的方法 Download PDF

Info

Publication number
CN105701340A
CN105701340A CN201610002449.0A CN201610002449A CN105701340A CN 105701340 A CN105701340 A CN 105701340A CN 201610002449 A CN201610002449 A CN 201610002449A CN 105701340 A CN105701340 A CN 105701340A
Authority
CN
China
Prior art keywords
activated carbon
gaseous state
containing compound
molecular
state sulfur
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201610002449.0A
Other languages
English (en)
Other versions
CN105701340B (zh
Inventor
汤立红
张贵剑
李凯
林强
宁平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Kunming University of Science and Technology
Original Assignee
Kunming University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Kunming University of Science and Technology filed Critical Kunming University of Science and Technology
Priority to CN201610002449.0A priority Critical patent/CN105701340B/zh
Publication of CN105701340A publication Critical patent/CN105701340A/zh
Application granted granted Critical
Publication of CN105701340B publication Critical patent/CN105701340B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/10Analysis or design of chemical reactions, syntheses or processes

Landscapes

  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Computing Systems (AREA)
  • Theoretical Computer Science (AREA)
  • Organic Low-Molecular-Weight Compounds And Preparation Thereof (AREA)
  • Treating Waste Gases (AREA)

Abstract

本发明公开了一种通过建立定量构效关系模型(Quantitative Structure-Activity Relationship,QSAR)预测气态含硫化合物常温下在活性炭上的吸附速率常数的方法;在已知气态含硫化合物分子结构的基础上,通过计算具有结构特征的分子描述符,并通过多元线性回归方法,构建了常温下活性炭对气态含硫化合物的吸附速率的QSAR模型,可快速、高效地预测气态含硫化合物在常温下在活性炭上的吸附速率常数。该方法简单、快捷、成本低,且能节省实验测试所需的人力、物力和财力;本预测模型简明、易于程序化;具有明确的应用领域、良好的拟合能力、稳健性和预测能力。

Description

预测气态含硫化合物常温下在活性炭上的吸附速率常数的方法
技术领域
本发明涉及一种通过建立定量构效关系模型(QSAR)预测气态含硫化合物常温下在活性炭上的吸附速率常数的方法,属于生态风险评价的定量结构与活性关系技术领域。
背景技术
有机化合物结构-活性定量相关的研究,最初作为定量药物设计的一个研究分支,是为了适应合理设计生物活性分支的需要而发展起来的。定量构效关系研究是应用化学计量学方法研究有机物的分子结构与理化性质或活性之间的定量相关关系,通过选取分子的理化参数或结构参数,用化学计量学和数理统计法研究有机化合物结构与其理化性质或生物活性之间的定量关系,建立定量构效模型方程来预测有机化合物分子理化性质或生物活性。它对于设计和筛选生物活性显著的药物,以及阐明药物的作用机理等均具有指导作用。特别是近二三十年来,由于计算机技术的发展和应用,使QSAR不仅已成为定量药物设计的一种重要方法,而且在环境化学、环境毒理学等领域中也得到了广泛的应用。许多环境科学研究者通过各种污染物结构-毒性定理关系的研究,建立了多种具有毒性预测能力的环境模型,如大连理工大学发明的专利“通过定量构效关系模型预测大气中有机物与羟基反应速率常数的方法”(中国专利申请号201310307098.0)和“一种通过定量构效关系模型预测有机物液相蒸汽压的方法”(中国专利申请号201110410088.0)及山东大学发明的专利“一种通过定量构效关系模型预测有机磷农药对水生生物急性毒性的方法”(中国专利申请号201410053184.8)。这对已进入环境的污染物及尚未投放市场的新化合物的生物活性、毒性乃至环境行为进行了成功的预测、评价和筛选,这些都说明QSAR在环境领域中已显示出极其广阔的应用前景。
硫的有机化合物广泛存在于大气环境中,来源可分为自然源和人为源。气态含硫有机化合物的人为源主要来自工业废气的排放,气态含硫有机化合物排放到大气环境中,可以进行一些物理和化学过程,从而导致它们在大气中消除或在大气中进一步转化,会对环境和生物造成非常严重的污染和危害,例如COS和CS2扩散到大气圈的平流层时,会通过光解-氧化作用生产SO2气体,这是酸雨的主要来源之一,与此同时有可能转化为硫酸盐的气溶胶,引起大气层中臭氧的损耗,加剧全球气候变化;而且当大气环境中的气态含硫有机化合物含量达到一定浓度时,可以侵袭人类的神经系统,会带来巨大的危害,它通过呼吸道、消化道和皮肤进入人体,作用于人体的各个器官,产生致畸、神经衰弱、神经性麻痹、胚胎发育障碍和子代先天性缺陷等症状,危及人体健康。有效去除气态含硫有机物迫在眉睫,特别是持久性气态有机硫。目前去除方法包括吸附、还原、氧化等方法,利用活性炭吸附去除的较为常见。因此,利用定量构效关系模型预测气态含硫化合物在常温下在活性炭上的吸附速率常数能够有效、及时的帮助人们了解活性炭对某种有机硫的吸附特性,以便提出适当的解决方式。但经检索,利用建立定量构效关系模型预测气态含硫化合物在常温下在活性炭上的吸附速率常数的方法还未见报道。
发明内容
针对现有技术上的不足,本发明的目的在于提供一种预测气态含硫化合物常温下在活性炭上的吸附速率常数的方法。
本发明方法具体步骤如下:
(1)通过查阅相关文献及书籍获得待测气态含硫化合物的分子结构信息,利用量子化学软件对待测气态含硫化合物进行几何结构优化,获得最优构型,从而得到相对分子质量M、轨道能量差、分子平衡电负性三个量子化学参数;其中为最低空轨道能量,为最高占据轨道能量;
(2)运用多元线性回归分析建立QSAR模型,获得如下回归方程,通过回归方程计算气态含硫化合物在常温下在活性炭上的吸附速率常数
,其中,N为分子中总原子数,为原子电负性,为分子中某个原子的原子数;
拟合能力:R2=0.9312。
本发明技术方案的原理是利用已知气态含硫有机化合物分子,运用量子化学软件对其进行几何全优化,得出与吸附相关的一些分子结构参数,并查得其他无法直接计算的参数。然后结合吸附速率实验值,利用多元线性回归分析方法,建立各种分子描述符与吸附速率常数之间的定量关系拟合方程,并对方程的拟合能力,预测能力进行验证。最后对模型的适合应用范围进行规范。由此,可以快捷、有效的气态含硫化合物在常温下在活性炭上的吸附速率常数。
本发明预测气态含硫化合物在常温下在活性炭上的吸附速率常数的方法通过如下步骤构建:
(1)通过进行低温水解实验或通过查阅相关数据库和文献,获得气态含硫有机化合物的吸附速率常数()数据(本工作共收集了14个气态含硫有机化合物的相关数据);
(2)利用量子化学软件对所要研究的气态含硫有机化合物进行几何结构优化,获得相对分子质量(M)、轨道能量差(为最高占据轨道能量)、分子平衡电负性()三个量子化学参数作为分子描述符;
(3)将步骤(1)得到的吸附速率数据中抽取1/4作为验证集数据,其余为训练集数据,训练集用来构建预测模型,验证集用来验证模型的预测能力;
(4)以步骤(2)获得的分子描述符为自变量,气态含硫化合物在常温下在活性炭上的吸附速率常数为因变量,运用多元线性回归分析建立QSAR模型,最后获得如下回归方程:
拟合能力:R2=0.9312;
如图1所示,对于训练集来说,实验值与预测值误差较小,说明此模型具有较好的预测性及准确性,因此,此模型可以用于预测气态含硫化合物在活性炭上的吸附速率常数。
(5)将验证集数据带入获得的回归方程(图2),得到待测气态含硫有机化合物常温下在活性炭上的吸附速率常数预测值,然后根据外部预测能力评价系数的值判定外部预测能力的好坏:当大于0.7时,表示建立的模型具有良好的外部预测能力,越大,外部预测能力越好。综合实验测定,最终得其外部预测能力=0.831,说明模型具有良好的外部预测能力。其中,上述外部预测能力评价系数的计算公式如下:
为验证集实验值,为验证集预测值,为训练集实验值均值,n为验证集个数,i表示第i个验证集)。
6)本预测模型适用于所有已知分子结构的气态含硫有机化合物。
本发明方法的优点和技术效果:
在已知气态含硫化合物分子结构的基础上,通过计算具有结构特征的分子描述符,并通过多元线性回归方法,构建了常温下活性炭对气态含硫化合物的吸附速率的QSAR模型,可快速、高效地预测气态含硫化合物在常温下在活性炭上的吸附速率常数;该方法简单、快捷、成本低,且能节省实验测试所需的人力、物力和财力;本预测模型简明、易于程序化;具有明确的应用领域、良好的拟合能力、稳健性和预测能力。
附图说明
图1为训练集气态含硫化合物在常温下在活性炭上的吸附速率常数的实验值与预测值的拟合图;
图2为验证集气态含硫化合物在常温下在活性炭上的吸附速率常数的实验值与预测值的拟合图。
具体实施方式
下面结合具体实施例进一步详细描述本发明,但本发明保护范围并不限于如下所述内容。
实施例1:预测二甲基三硫在常温下在活性炭上的吸附速率常数方法
首先查得二甲基三硫的分子结构信息(H3C-S-S-S-CH3),然后利用量子化学软件Gaussian09对分子结构进行优化,获得二甲基三硫结构的最优构型(对称构型,C-S的键长为1.89748?,S-S的键长为2.25417?,∠C-S-S=101.4815°,∠S-S-S=108.22126°),然后获得模型所需的描述符:相对分子质量(M为126)和轨道能量差(为4.18636eV,,其中为最低空轨道能量,为最高占据轨道能量)及分子平衡电负性(为2.3533)三个量子化学参数作为分子描述符;最后得到的吸附常数预测值为20.38647gm/g,而查得的实验值为21gm/g,误差仅为0.61353,与实验值非常相符。
实施例2:预测甲硫乙酸脂在常温下在活性炭上的吸附速率常数方法
首先查得甲硫乙酸脂的分子结构信息(S-C(H2)-C(H2)-O-C(O)-CH3),然后利用量子化学软件Gaussian09对分子结构进行优化,获得甲硫乙酸脂结构的最优构型(S-C的键长为1.88277?,∠S-C-C=113.42298°),然后获得模型所需的描述符:相对分子质量(M为120)和轨道能量差(为2.5976eV,,其中为最低空轨道能量,为最高占据轨道能量)及分子平衡电负性(为2.4295)三个量子化学参数作为分子描述符。最后得到的吸附常数预测值为6.265791gm/g,而查得的实验值为5.8gm/g,误差仅为0.465791,与实验值非常相符。
实施例3:预测乙硫醇在常温下在活性炭上的吸附速率常数方法
首先查得乙硫醇的分子结构信息(H3C-C(H2)-SH),然后利用量子化学软件Gaussian09对分子结构进行优化,获得乙硫醇结构的最优构型(C-C的键长为1.52244?,C-S的键长为1.91543?,∠C-C-S=113.40251°),然后获得模型所需的描述符:相对分子质量(M为62)和轨道能量差(为6.82584eV,,其中为最低空轨道能量,为最高占据轨道能量)及分子平衡电负性(为2.3669)三个量子化学参数作为分子描述符。最后得到的吸附常数预测值为30.74464gm/g,而查得的实验值为30gm/g,误差仅为0.744641,与实验值非常相符。
实施例4:预测3-甲基-2-丁烯-1-硫醇在常温下在活性炭上的吸附速率常数方法
首先查得3-甲基-2-丁烯-1-硫醇的分子结构信息(HS-C(H2)-C(H)=C(CH3)-CH3),然后利用量子化学软件Gaussian09对分子结构进行优化,获得3-甲基-2-丁烯-1-硫醇结构的最优构型(S-C的键长为1.94020?,∠S-C-C=112.77367°),然后获得模型所需的描述符:相对分子质量(M为102)和轨道能量差(为6.02181eV,,其中为最低空轨道能量,为最高占据轨道能量)及分子平衡电负性(为2.3209)三个量子化学参数作为分子描述符。最后得到的吸附常数预测值为31.1919gm/g,而查得的实验值为31.2gm/g,误差仅为0.0081,与实验值非常相符。

Claims (3)

1.一种预测气态含硫化合物常温下在活性炭上的吸附速率常数的方法,其特征在于步骤如下:
(1)通过查阅书籍及文献获得待测气态含硫化合物的分子结构信息,利用量子化学软件对待测气态含硫化合物进行几何结构优化,获得最优构型,从而得到相对分子质量M、轨道能量差、分子平衡电负性三个量子化学参数;其中为最低空轨道能量,为最高占据轨道能量;
(2)运用多元线性回归分析建立QSAR模型,获得如下回归方程,通过回归方程计算气态含硫化合物在常温下在活性炭上的吸附速率常数
,其中,N为分子中总原子数,为原子电负性,为分子中某个原子的原子数;
拟合能力:R2=0.9312。
2.根据权利要求1所述预测气态含硫化合物在常温下在活性炭上的吸附速率常数的方法,其特征在于:气态含硫化合物为有机气态含硫化合物。
3.根据权利要求2所述预测气态含硫化合物在常温下在活性炭上的吸附速率常数的方法,其特征在于:气态含硫化合物包括硫醚、硫醇、二甲基二硫、甲硫乙酸脂、磺酸。
CN201610002449.0A 2016-01-06 2016-01-06 预测气态含硫化合物常温下在活性炭上的吸附速率常数的方法 Active CN105701340B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610002449.0A CN105701340B (zh) 2016-01-06 2016-01-06 预测气态含硫化合物常温下在活性炭上的吸附速率常数的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610002449.0A CN105701340B (zh) 2016-01-06 2016-01-06 预测气态含硫化合物常温下在活性炭上的吸附速率常数的方法

Publications (2)

Publication Number Publication Date
CN105701340A true CN105701340A (zh) 2016-06-22
CN105701340B CN105701340B (zh) 2018-10-23

Family

ID=56226979

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610002449.0A Active CN105701340B (zh) 2016-01-06 2016-01-06 预测气态含硫化合物常温下在活性炭上的吸附速率常数的方法

Country Status (1)

Country Link
CN (1) CN105701340B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114496112A (zh) * 2022-01-21 2022-05-13 内蒙古工业大学 一种基于多目标优化的抗乳腺癌药物成分智能量化方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020052694A1 (en) * 1998-10-28 2002-05-02 Mcgregor Malcolm J. Pharmacophore fingerprinting in primary library design
CN103425872A (zh) * 2013-07-20 2013-12-04 大连理工大学 通过定量构效关系模型预测大气中有机物与羟基反应速率常数的方法
CN105069315A (zh) * 2015-08-26 2015-11-18 中国环境科学研究院 基于金属形态和有效性的水生生物毒性预测方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020052694A1 (en) * 1998-10-28 2002-05-02 Mcgregor Malcolm J. Pharmacophore fingerprinting in primary library design
CN103425872A (zh) * 2013-07-20 2013-12-04 大连理工大学 通过定量构效关系模型预测大气中有机物与羟基反应速率常数的方法
CN105069315A (zh) * 2015-08-26 2015-11-18 中国环境科学研究院 基于金属形态和有效性的水生生物毒性预测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
段强林,等: "含硫苯衍生物对发光细菌生物毒性的QSAR研究", 《南华大学学报(自然科学版)》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114496112A (zh) * 2022-01-21 2022-05-13 内蒙古工业大学 一种基于多目标优化的抗乳腺癌药物成分智能量化方法
CN114496112B (zh) * 2022-01-21 2023-10-31 内蒙古工业大学 一种基于多目标优化的抗乳腺癌药物成分智能量化方法

Also Published As

Publication number Publication date
CN105701340B (zh) 2018-10-23

Similar Documents

Publication Publication Date Title
Kinney Interactions of climate change, air pollution, and human health
Harshman et al. Storage stability of exhaled breath on Tenax TA
Cozma et al. Modeling and simulation of high pressure water scrubbing technology applied for biogas upgrading
Awang et al. Multivariate methods to predict ground level ozone during daytime, nighttime, and critical conversion time in urban areas
Carrico et al. Hygroscopic growth behavior of a carbon-dominated aerosol in Yosemite National Park
Pathak et al. Nighttime enhancement of PM2. 5 nitrate in ammonia-poor atmospheric conditions in Beijing and Shanghai: Plausible contributions of heterogeneous hydrolysis of N2O5 and HNO3 partitioning
Gu et al. Multi-scale volatile organic compound (VOC) source apportionment in Tianjin, China, using a receptor model coupled with 1-hr resolution data
Schneider et al. Combined computational and experimental investigation of SO x adsorption on MgO
Papurello et al. Biogas purification: a comparison of adsorption performance in D4 siloxane removal between commercial activated carbons and waste wood-derived char using isotherm equations
CN105678069A (zh) 预测气态含硫化合物在低温水解条件下消除速率常数的方法
Chung et al. Input–output surrogate models for efficient economic evaluation of amine scrubbing CO2 capture processes
Rosati et al. New particle formation and growth from dimethyl sulfide oxidation by hydroxyl radicals
Sciare et al. Diurnal and seasonal variation of atmospheric dimethylsulfoxide at Amsterdam Island in the southern Indian Ocean
Rodríguez-Sevilla et al. Absorption equilibria of dilute SO2 in seawater
CN105701340A (zh) 预测气态含硫化合物常温下在活性炭上的吸附速率常数的方法
Metz et al. Structure evolution of chemically degraded ZIF-8
Molino et al. Design and performance of an adsorption bed with activated carbons for biogas purification
Orhan et al. Accelerating the prediction of CO2 capture at low partial pressures in metal-organic frameworks using new machine learning descriptors
CN110038389A (zh) 一种有害废气的处理设备及处理方法
Tursumbayeva et al. Development of time-weighted average sampling of odorous volatile organic compounds in air with solid-phase microextraction fiber housed inside a GC glass liner: proof of concept
Widziewicz-Rzońca et al. Strongly and Loosely Bound Water in Ambient Particulate Matter—Qualitative and Quantitative Determination by Karl Fischer Coulometric Method
Qi et al. Multi-model Evaluation and Bayesian Model Averaging in Quantitative Air Quality Forecasting in Central China
Li et al. Reaction kinetics of hydrogen chloride with calcium oxide by fourier transform infrared spectroscopy
Sivtsova et al. Kinetics of carbon dioxide absorption from air in a flow reactor with a fixed bed of K 2 CO 3-based sorbent
Gorecki et al. The methods and stands for testing fixed sorbent and sorbent polymer composite materials for the removal of mercury from flue gases

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant