Badhe Yogesh, Kumar Dharmendr, Gupta Rakesh, Jain Vinay, Rai Beena
Physical Science Research Area, TCS Research, Sahyadri Park 2, Phase 3, Hinjewadi Rajiv Gandhi Infotech Park, Hinjewadi, Pimpri-Chinchwad, Pune, Maharashtra, 411057, India.
J Mol Model. 2024 May 9;30(6):162. doi: 10.1007/s00894-024-05952-x.
This study involves simulating the process of inhibiting corrosion through the formation of micelles by surfactants and their deposition on iron (Fe) surfaces. The primary focus is on examining CTAB/SDS mixtures in aqueous solutions with different concentrations. Micelle properties, including size, shape, aggregation number, cluster size, and surfactant diffusion, were calculated and validated with experimental data. The coarse-grained Fe surface was modeled and validated against experimental water contact-angle data. Subsequently, the deposition of CTAB/SDS mixtures on the Fe surface and air-water interface was studied systematically. We found that the relative ratio of CTAB/SDS in the solution directly influences surfactant deposition behavior, which might impact the corrosion inhibition efficiency.
All the MD simulations were performed using the GROMACS software with MARTINI2 force field and Martini polar water. The molecules are packed using PACKMOL software. Both NVT and NPT simulations are caried out at temperature and pressure of 303 K and 1 bar respectively, with a nonbonded interaction cut-off (r) of 1.1 nm. The LJ potential was shifted from 0.9 nm to r, while the electrostatic potential was shifted from 0.0 nm to r. For electrostatics, reaction-field coulomb type is used, relative dielectric constant (epsilon-r) and the reaction field dielectric constant (epsilon-rf) are equal to 2.5 and infinity respectively. The dielectric constant below r is epsilon-r, and beyond the cut-off is epsilon-rf. Coulomb-modifier used as potential-shift which leads to shift in the coulomb potential by a constant such that it is zero at the r. This makes the potential of the integral of the force . The neighbor list was updated every 10 steps, employing a neighbor list cut-off equal to r. Using a polar water model, we used a constant time step of 0.02 ps throughout the simulation. The used epsilon-r = 2.5, is recommended for polar water.
本研究涉及模拟表面活性剂通过形成胶束来抑制腐蚀的过程及其在铁(Fe)表面的沉积。主要重点是研究不同浓度水溶液中的十六烷基三甲基溴化铵(CTAB)/十二烷基硫酸钠(SDS)混合物。计算了胶束的性质,包括大小、形状、聚集数、簇尺寸和表面活性剂扩散,并与实验数据进行了验证。对粗粒度的铁表面进行了建模,并根据实验水接触角数据进行了验证。随后,系统地研究了CTAB/SDS混合物在铁表面和空气-水界面的沉积。我们发现溶液中CTAB/SDS的相对比例直接影响表面活性剂的沉积行为,这可能会影响缓蚀效率。
所有分子动力学(MD)模拟均使用带有MARTINI2力场和Martini极性水的GROMACS软件进行。使用PACKMOL软件对分子进行填充。分别在303 K和1巴的温度和压力下进行等温等容(NVT)和等温等压(NPT)模拟,非键相互作用截止值(r)为1.1 nm。Lennard-Jones(LJ)势从0.9 nm移至r,而静电势从0.0 nm移至r。对于静电作用,使用反应场库仑类型,相对介电常数(ε-r)和反应场介电常数(ε-rf)分别等于2.5和无穷大。低于r的介电常数为ε-r,超过截止值的为ε-rf。使用库仑修饰符作为势移,导致库仑势通过一个常数发生位移,使其在r处为零。这使得力的积分势为零。每10步更新一次邻居列表,采用等于r的邻居列表截止值。使用极性水模型,在整个模拟过程中我们使用0.02 ps的恒定时间步长。所使用的ε-r = 2.5,是极性水推荐使用的值。