期刊
ACTA PHYSICO-CHIMICA SINICA
卷 34, 期 10, 页码 1163-1170出版社
PEKING UNIV PRESS
DOI: 10.3866/PKU.WHXB201802271
关键词
Free energy change; Aggregation pathway; SDS; Micelle; Molecular dynamics calculation
资金
- FLAGSHIP, MEXT within Priority Study 5 (Development of New Fundamental Technologies for High-Efficiency Energy Creation, Conversion/Storage and Use) Using Computational Resources of the K Computer by the RIKEN Advanced Institute for Computational Science t [hp170241]
- MEXT KAKENHI [17K04758]
- Grants-in-Aid for Scientific Research [17K05748, 17K04758] Funding Source: KAKEN
Surfactant molecules, when dispersed in solution, have been shown to spontaneously form aggregates. Our previous studies on molecular dynamics (MD) calculations have shown that ionic sodium dodecyl sulfate molecules quickly aggregated even when the aggregation number is small. The aggregation rate, however, decreased for larger aggregation numbers. In addition, studies have shown that micelle formation was not completed even after a 100 ns-long MD run (Chem. Phys. Lett. 2016, 646, 36). Herein, we analyze the free energy change of micelle formation based on chemical species model combined with molecular dynamics calculations. First, the free energy landscape of the aggregation, Delta G(i+j)(dagger), where two aggregates with sizes i and j associate to form the (i + j)-mer, was investigated using the free energy of micelle formation of the i-mer, G(i)(dagger), which was obtained through MD calculations. The calculated Delta G(i+j)(dagger) was negative for all the aggregations where the sum of DS ions in the two aggregates was 60 or less. From the viewpoint of chemical equilibrium, aggregation to the stable micelle is desired. Further, the free energy profile along possible aggregation pathways was investigated, starting from small aggregates and ending with the complete thermodynamically stable micelles in solution. The free energy profiles, G(I, k), of the aggregates at /-th aggregation path and k-th state were evaluated by the formation free energy Sigma(i)n(i)(I,K)G(i)(dagger) and the free energy of mixing Sigma(i)n(i)(I,K)K(B)Tin(n(i)(I,K)In(I,K)), where n(l, k) is the number of i-mer in the system at the /-th aggregation path and k-th state, with n(I,k)=Sigma n(i)(I,k). All the aggregation pathways were obtained from the initial state of 12 pentamers to the stable micelle with i = 60. All the calculated G(I, k) values monotonically decreased with increasing k. This indicates that there are no free energy barriers along the pathways. Hence, the slowdown is not due to the thermodynamic stability of the aggregates, but rather the kinetics that inhibit the association of the fragments. The time required for a collision between aggregates, one of the kinetic factors, was evaluated using the fast passage time, t(FPT). The calculated t(FPT )was about 20 ns for the aggregates with N = 31. Therefore, if aggregation is a diffusion-controlled process, it should be completed within the 100 ns-simulation. However, aggregation does not occur due to the free energy barrier between the aggregates, that is, the repulsive force acting on them. This may be caused by electrostatic repulsions produced by the overlap of the electric double layers, which are formed by the negative charge of the hydrophilic groups and counter sodium ions on the surface of the aggregates.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据