十二烷基苯磺酸钠在SiO2表面聚集的分子动力学模拟

采用分子动力学方法研究了阴离子表面活性剂十二烷基苯磺酸钠(SDBS)在无定形SiO2固体表面的吸附.设置不同的水层厚度,观察同液界面和气液界面吸附的差异.模拟发现表面活性剂分子能够在短时间内吸附到SiO2表面,受碳链和固体表面之间相互作用的影响形成表面活性剂分子层,并依据吸附量的大小形成不同的聚集结构;在水层足够厚的情况下,由于有较多的表面活性剂分子吸附在固体表面,从而形成带有疏水核心的半胶束结构;计算得到的成对势表明极性头与钠离子或水分子之间的结合或解离与二者之间的能垒有关,解离能垒远大于结合能垒,引起更多Na+聚集在极性头周围而只有少数Na+存在于溶液中;无论气液还是固液界面,极性头均伸向水相,与水分子形成不同类型的氢键.模拟表明,分子动力学方法可以作为实验的一种补充,为实验提供必要的微观结构信息.     

十二烷基苯磺酸钠在SiO2表面聚集的分子动力学模拟

采用分子动力学方法研究了阴离子表面活性剂十二烷基苯磺酸钠(SDBS)在无定形SiO2固体表面的吸附. 设置不同的水层厚度, 观察固液界面和气液界面吸附的差异. 模拟发现表面活性剂分子能够在短时间内吸附到SiO2表面, 受碳链和固体表面之间相互作用的影响形成表面活性剂分子层, 并依据吸附量的大小形成不同的聚集结构; 在水层足够厚的情况下, 由于有较多的表面活性剂分子吸附在固体表面,从而形成带有疏水核心的半胶束结构; 计算得到的成对势表明极性头与钠离子或水分子之间的结合或解离与二者之间的能垒有关, 解离能垒远大于结合能垒, 引起更多Na+聚集在极性头周围而只有少数Na+存在于溶液中; 无论气液还是固液界面, 极性头均伸向水相, 与水分子形成不同类型的氢键. 模拟表明, 分子动力学方法可以作为实验的一种补充, 为实验提供必要的微观结构信息.     

抗衡离子对硫酸盐表面活性剂气/液界面性质影响的分子动力学模拟

采用全原子分子动力学方法研究了抗衡离子为第一主族离子(Li⁺、Na⁺、K⁺、Rb⁺和Cs⁺)的十二烷基硫酸盐表面活性剂的气/液界面性质. 通过分析体系中各组分的密度分布曲线, 考察表面活性剂单分子层在界面的聚集形态, 并利用径向分布函数分析了表面活性剂极性头基与抗衡离子间的相互作用. 研究结果表明: 随着抗衡离子半径的增大, 不同体系的界面水层厚度依次增加, 表面活性剂极性头基与抗衡离子形成的Stern和扩散层厚度也相应增加. 但表面活性剂吸附层的抗衡离子缔合度以及体系表面张力却随抗衡离子半径的增大而减小. 研究表明抗衡离子的差异对十二烷基硫酸盐表面活性剂气/液界面性质有很大影响.     

阴离子表面活性剂与阳离子的相互作用

用密度泛函理论, 在B3LYP/6-31G水平上, 对十二烷基磺酸盐和羧酸盐阴离子表面活性剂与阳离子(Na+, Ca2+, Mg2+)形成的离子对进行结构优化, 从分子水平上研究表面活性剂与阳离子之间的相互作用. 计算结果表明: 磺酸盐和羧酸盐表面活性剂均采用2:1型, 即极性头中两个氧原子与阳离子发生稳定结合; 在与阳离子结合之前, 表面活性剂分子上的α-亚甲基带有明显的负电荷, 因此将其归为极性头; 但在阳离子电荷诱导下, α-亚甲基转而带有部分弱正电荷, 使极性头范围缩小. 计算也发现, 表面活性剂尾链带有弱正电荷, 使胶束内核带有了部分极性, 利于表面活性剂在溶液中的聚集, 此种极性介于烷烃油相和水相的极性之间.     

分子动力学模拟在三次采油中的应用

采用分子动力学模拟在微观尺度上解释三次采油领域的相关实验现象和过程机制,对后续提高原油采收率的现场应用有着重要的理论指导意义.本文综述了本课题组近年来利用分子动力学模拟针对三次采油中若干基本科学问题的研究,模拟表明表面活性剂的耐盐能力主要由极性头周围极化出的水化层的性质决定;泡沫的稳定性直接取决于假乳液膜和油桥结构的稳定性,两者分别受表面活性剂亲水基与水分子之间的结合能力以及疏水基与油相间的相互作用强弱影响;水通道的形成是影响固体表面的润湿翻转的重要因素;凝胶颗粒在纳米孔内的运移受孔径、表面结构等多种因素影响;稠油的高黏度主要来源于沥青质和胶质的堆积形成的网状结构,不同类型表面活性剂通过影响内部和界面结构起到降黏效果.同时应当指出,模拟结果受制于分子动力学方法中模型大小和模拟时间、力场的精确性等因素的影响,仍然需要结合实验方法进一步证实.     

分子动力学模拟电解质对阴非离子表面活性剂界面行为的影响

通过分子动力学(MD)方法研究了不同类型电解质对阴非离子表面活性剂C₁₂EO₃C油水界面性能的影响。运用z轴质量密度分布、径向分布函数、分子间相互作用配位数、空间分布函数及均方根位移五种模拟参数来分析电解质与阴非离子表面活性剂的相互作用情况。研究表明,三种离子的加入均对水分子与表面活性剂亲水基形成的水化层结构产生影响,且从微观层面验证三种离子对表面活性剂亲水基相互作用强度大小顺序为Na⁺ < Ca²⁺ < Mg²⁺。通过扩散模拟结果可以较好地解释离子加入对界面张力平衡时间的影响情况。这对指导实验方向、制订最佳复配方案具有重要意义。     

阴离子Gemini表面活性剂在油/水界面行为的分子动力学模拟

采用分子动力学方法研究了磺酸盐型阴离子Gemini表面活性剂在油/水界面的吸附行为, 考察了不同长度的连接基(Spacer)对表面活性剂在界面的聚集形态及界面性质的影响. 密度分布和微观结构信息显示, Gemini表面活性剂能在油/水界面形成单层膜结构. Gemini表面活性剂能使油/水界面的厚度显著增大, 并使界面形成能降低. 当连接基为6个碳时, 此类磺酸盐型Gemini表面活性剂的界面厚度最大, 形成的界面最稳定. 连接基长度对Gemini表面活性剂单层膜周围的水分子和Na⁺的吸附结构影响不大, 但是能影响水分子的扩散行为.     

气液界面上阴离子表面活性剂单层膜的分子动力学模拟

用分子动力学方法研究了阴离子表面活性剂十二烷基硫酸钠(SDS)在气液界面上的结构和动力学性质. 选择单分子占有面积分别为0.45和0.68 nm²的两个模拟体系, 通过径向分布函数表征了单层膜的厚度, 并根据疏水链中碳原子与极性头中硫原子之间组成的矢量分布和取向函数, 对比了不同界面单层膜的有序排列情况. 结果表明在分子占有面积较小达到饱和吸附的情况下, 界面上的SDS具有较好的有序性. 通过计算气液界面附近水分子的扩散系数发现: 由于氢键和静电作用的影响, 界面区域内的水分子较本体溶液中的水分子有较弱的迁移能力.     

表面活性剂在非极性有机溶剂中的自组装

本文总结了表面活性剂在非极性有机溶剂中复杂反相聚集的研究进展。首先突破了表面活性剂在非极性溶剂(油)中溶解的难题,通过设计合成大头基的新表面活性剂,或引入合适添加剂使之与表面活性剂头基相互作用,由此增大头基有效尺寸,这些均能有效促进表面活性剂形成核-壳完整的聚集体,进而带动表面活性剂分散(溶解)在非极性溶剂中。基于聚集体带动溶解的思路,建立了制备表面活性剂/油均相溶液的直接溶解方法,讨论了制备方法的关键要素,它比文献常用的甲醇预溶解法方便且有效。列举了若干典型的表面活性剂/环己烷均相体系,以此评述了聚集体带动溶解的方法,也展现了丰富多样的反相聚集形貌,讨论了表面活性剂头基尺寸对聚集结构的影响。     

阴离子表面活性剂在水/正烷烃界面的分子动力学模拟(英文)

利用分子动力学模拟方法研究了阴离子表面活性剂在水/正烷烃(壬烷,癸烷和十一碳烷)界面的结构和动力学特点.十六烷基苯磺酸钠作为研究对象,其中苯磺酸基团在十六碳烷的第4号碳原子上,记作4-C16.分析了不同油相和特定盐度条件下正烷烃-表面活性剂-水体系的界面特点(如密度剖面图、界面张力和径向分布函数).模拟结果表明平衡模型体系展现了一个很好的水/正烷烃界面.当加氯化钠到水溶液中,正烷烃-表面活性剂-水体系的界面张力有微小的变化,有趣的是表面活性剂二面角的反式结构分数的变化联系着界面张力的微小变化.可见,表面活性剂在界面处的结构对降低界面张力起到重要的作用.此外,还发现表面活性剂的极性头与钠离子和水分子存在较强的相互作用.     

Molecular dynamics study of the effect of calcium ions on the monolayer of SDC and SDSn surfactants at the vapor/liquid interface

The effect of Ca(2+) ions on the hydration shell of sodium dodecyl carboxylate (SDC) and sodium dodecyl sulfonate (SDSn) monolayer at vapor/liquid interfaces was studied using molecular dynamics simulations. For each surfactant, two different surface concentrations were used to perform the simulations, and the aggregation morphologies and structural details have been reported. The results showed that the aggregation structures relate to both the surface coverage and the calcium ions. The divalent ions can screen the interaction between the polar head and Na(+) ions. Thus, Ca(2+) ions locate near the vapor/liquid interface to bind to the headgroup, making the aggregations much more compact via the salt bridge. The potential of mean force (PMF) between Ca(2+) and the headgroups shows that the interaction is decided by a stabilizing solvent-separated minimum in the PMF. To bind to the headgroup, Ca(2+) should overcome the energy barrier. Among contributions to the PMF, the major repulsive interaction was due to the rearrangement of the hydration shell after the calcium ions entered into the hydration shell of the headgroup. The PMFs between the headgroup and Ca(2+) in the SDSn systems showed higher energy barriers than those in the SDC systems. This result indicated that SDSn binds the divalent ions with more difficulty compared with SDC, so the ions have a strong effect on the hydration shell of SDC. That is why sulfonate surfactants have better efficiency in salt solutions with Ca(2+) ions for enhanced oil recovery.     

髓鞘碱性蛋白对细胞膜脂质分子DPPC、DPPE单层膜影响机理研究

本文通过Langmuir单层膜的表面压力-平均分子面积(π-A)曲线的测定与分析,分别对髓鞘碱性蛋白(MBP)与细胞膜中不同头部基团脂质分子二棕榈酰基磷脂胆碱(DPPC)和二棕榈酰基磷脂酰乙醇胺(DPPE)在空气/液体界面上的相互作用过程进行了系统研究.实验结果表明:(1)当界面上脂质含量一定时,亚相中随着MBP浓度的增大,DPPC、DPPE单层膜的等温线向平均分子面积较大的方向移动;(2)在单层膜表面压力为10 mN/m时,一个MBP分子分别结合140±3个DPPC分子和100±3个DPPE分子,随着表面压力增大,当MBP分子分别与两种磷脂分子相互作用时,MBP插入到磷脂单层界面的个数逐渐减少;(3)随着蛋白质浓度的增加,脂分子形成的单层膜变得较为疏松,且MBP分子易于插入到分子头部较小的DPPE单层膜中;(4)蛋白质的存在使DPPC单层膜的表面压力逐渐减小,且蛋白质浓度越大表面压力降低越多,DPPC被MBP带入到亚相中越多;(5)对于DPPE单层膜,蛋白质通过与DPPE相互作用插入到界面膜中,引起表面压力增大,且蛋白质浓度越高,压力变化量越大.     

Formation and interaction of hydrated alkali metal ions at the graphite-water interface

Ion hydration at a solid surface ubiquitously exists in nature and plays important roles in many natural processes and technological applications. Aiming at obtaining a microscopic insight into the formation of such systems and interactions therein, we have investigated the hydration of alkali metal ions at a prototype surface-graphite (0001), using first-principles molecular dynamics simulations. At low water coverage, the alkali metal ions form two-dimensional hydration shells accommodating at most four (Li, Na) and three (K, Rb, Cs) waters in the first shell. These two-dimensional shells generally evolve into three-dimensional structures at higher water coverage, due to the competition between hydration and ion-surface interactions. Exceptionally K was found to reside at the graphite-water interface for water coverages up to bulk water limit, where it forms an "umbrellalike" surface hydration shell with an average water-ion-surface angle of 115 degrees . Interactions between the hydrated K and Na ions at the interface have also been studied. Water molecules seem to mediate an effective ion-ion interaction, which favors the aggregation of Na ions but prevents nucleation of K. These results agree with experimental observations in electron energy loss spectroscopy, desorption spectroscopy, and work function measurement. In addition, the sensitive dependence of charge transfer on dynamical structure evolution during the hydration process, implies the necessity to describe surface ion hydration from electronic structure calculations.     

Solvation of calcium ions in methanol-water mixtures: molecular dynamics simulation

Molecular dynamics simulations of CaCl2 solutions in water and methanol-water mixtures, with methanol concentrations of 5, 10, 50, and 90 mol %, at room temperature, have been performed. The methanol and water molecules have been modeled as flexible three-site bodies. Solvation of the calcium ions has been discussed on the basis of the radial and angular distribution functions, the orientation of the solvent molecules, and their geometrical arrangement in the coordination shells. Analysis of the H-bonds of the solvent molecules coordinated by Ca2+ has been done. Residence time of the solvent molecules in the coordination shell has been calculated. The preferential hydration of the calcium ions has been found over the whole range of the mixture composition. The water concentration in the first and second coordination shells of Ca2+ significantly exceeds the water content in the solution, despite the very similar interaction energy of the calcium ion with water and methanol. In aqueous solution and methanol-water mixtures, the first coordination shell of Ca2+ is irregular and long-living. The solvent molecules prefer the anti-dipole arrangement, but, in aqueous solutions and water-rich mixtures, the water molecules in the primary shell have only one H-bonded neighbor.     

Revisiting the structure and dynamics of hydrated Cd2+ in aqueous solutions: Insights from the RI-SCS-MP2/MM molecular dynamics simulation

The spin component scale MP2/molecular mechanics molecular dynamics simulation investigated the hydration shell formation and hydrated Cd²⁺ dynamics in the water environment. At the first hydration shell, six water molecules with 2.27 Å for the average distance between water and Cd²⁺. Dynamical properties were analyzed by computing the water molecule's mean residence time (MRT) in its first and second hydration shells. The MRT of each shell was determined to be 31.8 and 1.92 ps, suggesting the strong influence of Cd²⁺ in the first hydration shell. The second shell was labile, with an average number of water molecules being 18. Despite the strong interaction between Cd²⁺ and water molecules in the first shell, the influence of ions in the second hydration shell remained weak.     

Effect of Na+ and Ca2+ Ions on a Lipid Langmuir Monolayer: An Atomistic Description by Molecular Dynamics Simulations

Studying the effect of alkali and alkaline‐earth metal cations on Langmuir monolayers is relevant from biophysical and nanotechnological points of view. In this work, the effect of Na⁺ and Ca²⁺ on a model of an anionic Langmuir lipid monolayer of dimyristoylphosphatidate (DMPA^?) is studied by molecular dynamics simulations. The influence of the type of cation on lipid structure, lipid–lipid interactions, and lipid ordering is analyzed in terms of electrostatic interactions. It is found that for a lipid monolayer in its solid phase, the effect of the cations on the properties of the lipid monolayer can be neglected. The influence of the cations is enhanced for the lipid monolayer in its gas phase, where sodium ions show a high degree of dehydration compared with calcium ions. This loss of hydration shell is partly compensated by the formation of lipid–ion–lipid bridges. This difference is ascribed to the higher charge‐to‐radius ratio q/r for Ca²⁺, which makes ion dehydration less favorable compared to Na⁺. Owing to the different dehydration behavior of sodium and calcium ions, diminished lipid–lipid coordination, lipid–ion coordination, and lipid ordering are observed for Ca²⁺ compared to Na⁺. Furthermore, for both gas and solid phases of the lipid Langmuir monolayers, lipid conformation and ion dehydration across the lipid/water interface are studied.     

Dielectric saturation of the ion hydration shell and interaction between two double helices of DNA in mono- and multivalent electrolyte solutions: foundations of the epsilon-modified Poisson-Boltzmann theory

Potentials of mean force between single Na+, Ca2+, and Mg2+ cations and a highly charged spherical macroion in SPC/E water have been determined using molecular dynamics simulations. Results are compared to the electrostatic energy calculations for the primitive polarization model (PPM) of hydrated cations describing the ion hydration shell as a dielectric sphere of low permittivity (Gavryushov, S.; Linse, P. J. Phys. Chem. B 2003, 107, 7135). Parameters of the ion dielectric sphere and radius of the macroion/water dielectric boundary were extracted by means of this comparison to approximate the short-range repulsion of ions near the interface. To explore the counterion distributions around a simplified model of DNA, the obtained PPM parameters for Na+ and Ca2+ have been substituted into the modified Poisson-Boltzmann (MPB) equations derived for the PPM and named the epsilon-MPB (epsilon-MPB) theory. epsilon-MPB results for DNA suggest that such polarization effects are important in the case of 2:1 electrolyte and highly charged macromolecules. The three-dimensional implementation of the epsilon-MPB theory was also applied to calculation of the energies of interaction between two parallel macromolecules of DNA in solutions of NaCl and CaCl2. Being compared to results of MPB calculations without the ion polarization effects, it suggests that the ion hydration shell polarization and inhomogeneous solvent permittivity might be essential factors in the experimentally known hydration forces acting between charged macromolecules and bilayers at separations of less than 20 A between their surfaces.     

Electrostatics of B-DNA in NaCl and CaCl2 solutions: ion size, interionic correlation, and solvent dielectric saturation effects

The epsilon-modified Poisson-Boltzmann (-MPB) equations ( J. Phys. Chem. B, 2007, 111, 5264) have been solved on a three-dimensional grid for an all-atom geometry model of B-DNA. The approach is based on the implicit solvent model including finite sizes of hydrated ions and a dielectric approximation of the ion hydration shell. Results were obtained for the detailed geometry model of B-DNA in dilute and moderately concentrated solutions of NaCl and CaCl(2). All -MPB parameters of ions and dielectric medium were extracted from published results of all-atom molecular dynamics simulations. The study allows evaluations of the ion size, interionic correlation, and the solvent dielectric saturation effects on the ion distributions around DNA. It unambiguously suggests that the difference between the -MPB and Poisson-Boltzmann distributions of ions is low for Na(+) counterions. Such a difference in the case of divalent counterions Ca(2+) is dramatic: the dielectric saturation of the ion hydration shell leads to point-like adsorption of Ca(2+) on the phosphate groups of DNA. The -MPB equations were also applied to calculate the energy of interaction between two B-DNA molecules. Results agree with previously published simulations and experimental data. Some aspects of ion specificity of polyelectrolyte properties are discussed.     

Construction and molecular modeling of phospholipid surfaces

A protocol is given for the construction of phospholipid surfaces that possess variable head groups and thus variable net charge. Ab initio quantum mechanical calculations are performed to establish the necessary force field (AMBER) parameters. The charge distribution is defined by an electrostatic potential method consistent with the ab initio wave function. As a model calculation, a monolayer surface with head groups of phosphatidylserine and phosphatidylcholine derived from the crystal structure of 1,2-dilauroyl-DL-phosphatidylethanolamine (DLPE) is placed in a water bath with two Ca(II) ions present. The resultant surface is energy-optimized followed by 64 ps of molecular dynamics integration. Evaluation of calcium ion coordination environments, characterization of the P-N dipole inclination with respect ot the plane of the monolayer, and calculation of molecular surface area is performed and compared with experimental data.     

Exploration of the Ca2+ interaction modes of the nifedipine calcium channel antagonist.

A comprehensive study is carried out using quantum chemical computation and molecular dynamics (MD) simulations to gain insight into the interaction between Ca(2+) ions and the most important class of calcium channel antagonists--nifedipine. First, the chelating structures and energetic characters of nifedipine-Ca(2+) in the gas phase are explored, and 25 isomers are found. The most favorable chelating mode is a tridentate one, that is, Ca(2+) binds to two carbonyl O atoms and one nitryl O atom, where Ca(2+) is above the plane of the three O atoms to form a pyramidal structure. Accurate geometric structures, relative stabilities, vertical and adiabatic binding energies, and charge distributions are discussed. The differences in the geometries and energies among these isomers are analyzed from the contributions of chelating sites, electrostatics and polarizations, steric repulsions, and charge distributions. The interconversions among isomers with similar geometries and energies are also investigated because of the importance of the geometric transformation in the biological system. Furthermore, certain numbers of water molecules are added to the nifedipine-Ca(2+) system to probe the effect of water. A detailed study is performed on the hydrated geometries on the basis of the most stable isomer 1. Stepwise hydration can weaken the nifedipine-Ca(2+) interaction, and the chelating sites of nifedipine are gradually replaced by the added water molecules. Hexacoordination is found to be the most favorable geometry no matter how many water molecules were added, which can be verified by the MD simulations. The transfer of water molecules from the inner shell to the outer shell is also supported by MD simulations of the hexahydrated complexes.