色谱 ›› 2022, Vol. 40 ›› Issue (3): 281-288.DOI: 10.3724/SP.J.1123.2021.10011

• 研究论文 • 上一篇    下一篇

基于随机扩散理论的气相色谱分离模拟

孙寅璐1, 王琳1, 银芷玉1, 赵健伟2,*()   

  1. 1.辽宁大学化学院, 辽宁 沈阳 110036
    2.嘉兴学院材料与纺织工程学院, 浙江 嘉兴 314001
  • 收稿日期:2021-10-14 出版日期:2022-03-08 发布日期:2022-03-04
  • 通讯作者: 赵健伟
  • 基金资助:
    辽宁省自然科学基金项目(LJKZ0091);浙江省自然科学基金项目(LY19B030006);嘉兴市科技计划项目(2019AD32017);嘉兴市科技计划项目(2020AD10017)

Simulation of gas chromatographic separation based on random diffusion

SUN Yinlu1, WANG Lin1, YIN Zhiyu1, ZHAO Jianwei2,*()   

  1. 1. College of Chemistry, Liaoning University, Shenyang 110036, China
    2. College of Materials and Textile Engineering, Jiaxing University, Jiaxing 314001, China
  • Received:2021-10-14 Online:2022-03-08 Published:2022-03-04
  • Contact: ZHAO Jianwei
  • Supported by:
    Liaoning Province Natural Science Foundation, China(LJKZ0091);Zhejiang Province Natural Science Foundation, China(LY19B030006);Science and Technology Program of Jiaxing, China(2019AD32017);Science and Technology Program of Jiaxing, China(2020AD10017)

摘要:

色谱分离过程中的粒子扩散问题是色谱动力学研究的基础,深入理解粒子的扩散行为对优化分离操作条件、提升色谱性能和开发新型色谱柱尤为关键。现有的模拟方法多集中于局部过程的热力学研究,而整体的扩散分离过程报道并不多见。为此,该文基于微尺度受限空间内随机扩散的方法,通过动态追踪粒子的运动轨迹,实现粒子在气相色谱开管柱内的扩散全过程模拟。基于前期烷烃同系物的分离模拟研究,结合Kovats保留指数,分别建立了吸附步数与温度、吸附步数与碳数的函数关系,由此获得不同类型的同系物在不同温度条件下的分离参数系统。以醇类同系物的分离验证模拟的可靠性,结果表明保留时间的相对误差基本控制在5%以内,而峰宽相对误差在0.75%~60%之间。峰宽误差较大的原因在于:(1)参数化计算过程中未能充分迭代以及使用外推法;(2)模型中忽略了醇分子之间的氢键作用。该文提出的模拟方法虽然可以准确地预测色谱保留时间以及合理描述色谱峰的基本形貌特征,但尚有进一步发展空间,特别是增加对分子间相互作用的细节处理,例如可参考分子力学的方法建立分子间势函数和吸附步数的关系,利用分子力学计算的能量来取代参数化的吸附步数,从而实现更为精确的分离过程模拟。总体而言,该文所提出的模拟方法为优化色谱分离操作条件和开发新型色谱分离技术提供了有价值的参考。

关键词: 气相色谱, 分离, 扩散, 随机模拟

Abstract:

Understanding the diffusion behavior of particles during chromatographic analysis is critical for optimizing the operation conditions, improving the chromatographic performance, and designing a new separation device. Most of the existing simulations focus on chromatographic thermodynamics, while very few consider the overall diffusion and separation process. Herein, a new simulation method for gas chromatography separation was developed based on the random diffusion theory in microscale restricted space. This method retained the key information for controlling the diffusion behavior, simplified the interaction between the particles to be separated, and enlarged the time scale of each step one molecule walks. Thus, the operational efficiency could be significantly improved due to reduced calculation, and the entire diffusion process in the gas chromatography capillary column could be simulated. In this model, the capillary column was represented by a two-dimensional long and narrow rectangle where the particles to be separated randomly diffused. Besides, a directional velocity along the axis of the chromatographic column exerted on the particle represented the driving force of the mobile phase. If a particle collided with the inner wall of the column, it would remain at the collision position even after some time lapsed. When desorption occured, the particle would flow along with the mobile phase until its next adsorption on the surface. The interaction between the particle and the inner wall was expressed by the equivalent adsorption time. By dynamically tracking the trajectories of particles, the statistical distribution of time for the residence of the particles in the chromatographic column could be obtained, that is, the detection signal (retention time). Based on the previous simulation studies on the separation of n-alkane homologues, combined with the Kovats Retention Index, the functional relationships between adsorption steps and temperature as well as carbon number were established. As a result, the separation parameter systems for various homologues at different temperatures were set up. The separation of alcohol homologues at different temperatures was considered as an example to verify the reliability of the simulation method. The relative errors between the measured and simulated values were within 5% for the retention time and 0.75%-60% for the peak width. The reasons for the large relative errors in the peak width are summarized as follows. On the one hand, parameterization of alcohol homologues was performed on the basis of a previous study on the separation law of n-alkane. Given the limitations of the current computing capability, the insufficient iteration in the parameterized process led to large errors. In addition, the errors at different temperatures further accumulated in extrapolated approximations. On the other hand, the strong hydrogen bond force between the alcohol molecules was simplified with the elastic collision, which increased the magnitude of the errors. Although the simulation method proposed in this paper can accurately predict the retention time and reasonably describe the morphological characteristics of chromatographic peaks, there is still room for improvement. In particular, the description of the detailed interactions between molecules must be improved. For example, the method of molecular mechanics may be assistant with the investigation of the functional relationship between interaction potential and adsorption steps. The interaction potential calculated on the basis of molecular mechanics replaces the parameterized adsorption steps, yielding more accurate simulation results. In general, the simulation method proposed in this study is a valuable reference for the optimization of chromatographic operating conditions and for the development of new chromatographic techniques.

Key words: gas chromatography (GC), separation, diffusion, random simulation

中图分类号: