1
|
Peroukidis SD, Stott IP, Mavrantzas VG. Coarse-Grained Model Incorporating Short- and Long-Range Effective Potentials for the Fast Simulation of Micelle Formation in Solutions of Ionic Surfactants. J Phys Chem B 2022; 126:5555-5569. [PMID: 35838193 DOI: 10.1021/acs.jpcb.2c02751] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A coarse-grained model comprising short- and long-range effective potentials, parametrized with the iterative Boltzmann inversion (IBI) method, is presented for capturing micelle formation in aqueous solutions of ionic surfactants using as a model system sodium dodecyl sulfate (SDS). In the coarse-grained (CG) model, each SDS molecule is represented as a sequence of four beads while each water molecule is modeled as a single bead. The proposed CG scheme involves ten potential energy functions: four of them describe bonded interactions and control the distribution functions of intramolecular degrees of freedom (bond lengths, valence angles, and dihedrals) along an SDS molecule while the other six account for intermolecular interactions between pairs of SDS and water beads and control the radial distribution functions. The nonbonded effective potentials between coarse-grained SDS molecules extend up to about 12 nm and capture structural and morphological features of the micellar solution both at short and long distances. The long-range component of these potentials, in particular, captures correlations between surfactant molecules belonging to different micelles and is essential to describe ordering associated with micelle formation. A new strategy is introduced for determining the effective potentials through IBI by using information (target distribution functions) extracted from independent atomistic simulations of a micellar reference system (a salt-free SDS solution at total surfactant concentration cT equal to 103 mM, temperature T equal to 300 K, and pressure P equal to 1 atm) obtained through a multiscale approach described in an earlier study. It employs several optimization steps for bonded and nonbonded interactions and a gradual parametrization of the short- and long-range components of the latter, followed by reparametrization of the bonded ones. The proposed CG model can reproduce remarkably accurately the microstructure and morphology of the reference system within only a few hours of computational time. It is therefore very promising for future studies of structural and morphological behavior of various liquid surfactant formulations.
Collapse
Affiliation(s)
- Stavros D Peroukidis
- Department of Chemical Engineering, University of Patras and FORTH-ICE/HT, GR 26504, Patras, Greece
| | - Ian P Stott
- Unilever Research and Development Port Sunlight, Bebington CH63 3JW, United Kingdom
| | - Vlasis G Mavrantzas
- Department of Chemical Engineering, University of Patras and FORTH-ICE/HT, GR 26504, Patras, Greece.,Particle Technology Laboratory, Department of Mechanical and Process Engineering, ETH Zürich, CH-8092 Zürich, Switzerland
| |
Collapse
|
2
|
Zhou J, Ranjith PG. Self-assembly and viscosity changes of binary surfactant solutions: A molecular dynamics study. J Colloid Interface Sci 2020; 585:250-257. [PMID: 33285463 DOI: 10.1016/j.jcis.2020.11.022] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2020] [Revised: 11/04/2020] [Accepted: 11/05/2020] [Indexed: 11/15/2022]
Abstract
HYPOTHESIS Structure and self-assembly of surfactants in the solution shows a fundamental influence on its viscosity. Through molecular simulations using Martini force field, synergistic effects in aggregation as well as the viscosity changes of a binary ionic surfactant system can be modelled. Simulations: Coarse-grained molecular dynamics simulations are performed to model the SDS/CAPB binary surfactant solution, and both equilibrium and non-equilibrium methods are used to calculate the viscosity of the equilibrated micellar systems. FINDINGS Our simulation results indicate that the new version of the Martini force field can provide more reasonable self-assembly of surfactant, both single and binary system. Synergistic effects in micelle formation for SDS/CAPB have been successfully reproduced, that is, the formation of cylindrical micelles or even wormlike micelles at a lower concentration when compared with the pure system. Meanwhile, both equilibrium and non-equilibrium methods provide quantitatively comparable viscosity for each system. For pure micellar system, the viscosity linearly increases with the total concentration. Nevertheless, our simulation fails to capture the viscosity enhancement of the solution in corresponding with the formation of rodlike or wormlike micelles, and a full parameter optimization of force field is still necessary.
Collapse
Affiliation(s)
- Jun Zhou
- Department of Civil Engineering, Monash University, Building 60, Melbourne, Victoria 3800, Australia
| | - P G Ranjith
- Department of Civil Engineering, Monash University, Building 60, Melbourne, Victoria 3800, Australia.
| |
Collapse
|
3
|
Peroukidis SD, Tsalikis DG, Noro MG, Stott IP, Mavrantzas VG. Quantitative Prediction of the Structure and Viscosity of Aqueous Micellar Solutions of Ionic Surfactants: A Combined Approach Based on Coarse-Grained MARTINI Simulations Followed by Reverse-Mapped All-Atom Molecular Dynamics Simulations. J Chem Theory Comput 2020; 16:3363-3372. [DOI: 10.1021/acs.jctc.0c00229] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Stavros D. Peroukidis
- Department of Chemical Engineering, University of Patras and FORTH-ICE/HT, Patras, GR 26504, Greece
- Hellenic Open University, Patras, GR 26222, Greece
| | - Dimitrios G. Tsalikis
- Department of Chemical Engineering, University of Patras and FORTH-ICE/HT, Patras, GR 26504, Greece
| | - Massimo G. Noro
- UKRI Science and Technology Facilities Council, Daresbury WA4 4AD, U.K
| | - Ian P. Stott
- Unilever Research & Development Port Sunlight, Bebington CH63 3JW, U.K
| | - Vlasis G. Mavrantzas
- Department of Chemical Engineering, University of Patras and FORTH-ICE/HT, Patras, GR 26504, Greece
- Department of Mechanical and Process Engineering, Particle Technology Laboratory, ETH Zürich, CH-8092 Zürich, Switzerland
| |
Collapse
|
4
|
Eslami H, Khani M, Müller-Plathe F. Gaussian Charge Distributions for Incorporation of Electrostatic Interactions in Dissipative Particle Dynamics: Application to Self-Assembly of Surfactants. J Chem Theory Comput 2019; 15:4197-4207. [DOI: 10.1021/acs.jctc.9b00174] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Hossein Eslami
- Department of Chemistry, College of Sciences, Persian Gulf University, Boushehr 75168, Iran
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, Alarich-Weiss-Str. 8, 64287 Darmstadt, Germany
| | - Marzieh Khani
- Department of Chemistry, College of Sciences, Persian Gulf University, Boushehr 75168, Iran
| | - Florian Müller-Plathe
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, Alarich-Weiss-Str. 8, 64287 Darmstadt, Germany
| |
Collapse
|
5
|
Anderson RL, Bray DJ, Del Regno A, Seaton MA, Ferrante AS, Warren PB. Micelle Formation in Alkyl Sulfate Surfactants Using Dissipative Particle Dynamics. J Chem Theory Comput 2018; 14:2633-2643. [PMID: 29570296 DOI: 10.1021/acs.jctc.8b00075] [Citation(s) in RCA: 62] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We use dissipative particle dynamics (DPD) to study micelle formation in alkyl sulfate surfactants, with alkyl chain lengths ranging from 6 to 12 carbon atoms. We extend our recent DPD force field [ J. Chem. Phys. 2017 , 147 , 094503 ] to include a charged sulfate chemical group and aqueous sodium ions. With this model, we achieve good agreement with the experimentally reported critical micelle concentrations (CMCs) and can match the trend in mean aggregation numbers versus alkyl chain length. We determine the CMC by fitting a charged pseudophase model to the dependence of the free surfactant on the total surfactant concentration above the CMC and compare it with a direct operational definition of the CMC as the point at which half of the surfactant is classed as micellar and half as monomers and submicellar aggregates. We find that the latter provides the best agreement with experimental results. Finally, with the same model, we are able to observe the sphere-to-rod morphological transition for sodium dodecyl sulfate (SDS) micelles and determine that it corresponds to SDS concentrations in the region of 300-500 mM.
Collapse
Affiliation(s)
- Richard L Anderson
- STFC Hartree Centre, Scitech Daresbury , Warrington WA4 4AD , United Kingdom
| | - David J Bray
- STFC Hartree Centre, Scitech Daresbury , Warrington WA4 4AD , United Kingdom
| | - Annalaura Del Regno
- STFC Hartree Centre, Scitech Daresbury , Warrington WA4 4AD , United Kingdom
| | - Michael A Seaton
- STFC Hartree Centre, Scitech Daresbury , Warrington WA4 4AD , United Kingdom
| | - Andrea S Ferrante
- Ferrante Scientific Ltd. , 5 Croft Lane , Bromborough CH62 2BX , United Kingdom
| | - Patrick B Warren
- Unilever R&D Port Sunlight , Quarry Road East , Bebington CH63 3JW , United Kingdom
| |
Collapse
|
6
|
Cosby T, Vicars Z, Wang Y, Sangoro J. Dynamic-Mechanical and Dielectric Evidence of Long-Lived Mesoscale Organization in Ionic Liquids. J Phys Chem Lett 2017; 8:3544-3548. [PMID: 28715220 DOI: 10.1021/acs.jpclett.7b01516] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Experimental evidence of the dynamics of mesoscopic structure in room-temperature ionic liquids-a feature expected to correlate with many physicochemical properties of these materials-remains limited. Here, we report the observation of slow, sub-α relaxations corresponding to dynamics of nanoscale hydrophobic aggregates in a systematic series of 1-alkyl-3-methylimidazolium-based ionic liquids from detailed analysis of dynamic-mechanical and broad-band dielectric spectra. The emergence of the sub-α relaxations correlates with increases in the zero-shear viscosity and static dielectric permittivity, constituting direct evidence of the influence of mesoscale aggregation on the physicochemical properties of ionic liquids.
Collapse
Affiliation(s)
- Tyler Cosby
- Department of Chemical and Biomolecular Engineering, University of Tennessee , Knoxville, Tennessee 37996, United States
| | - Zachariah Vicars
- Department of Chemical and Biomolecular Engineering, University of Tennessee , Knoxville, Tennessee 37996, United States
| | - Yangyang Wang
- Center for Nanophase Materials Science, Oak Ridge National Laboratory , Oak Ridge, Tennessee 37831, United States
| | - Joshua Sangoro
- Department of Chemical and Biomolecular Engineering, University of Tennessee , Knoxville, Tennessee 37996, United States
| |
Collapse
|
7
|
Mao R, Lee MT, Vishnyakov A, Neimark AV. Modeling Aggregation of Ionic Surfactants Using a Smeared Charge Approximation in Dissipative Particle Dynamics Simulations. J Phys Chem B 2015; 119:11673-83. [PMID: 26241704 DOI: 10.1021/acs.jpcb.5b05630] [Citation(s) in RCA: 68] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Using dissipative particle dynamics (DPD) simulations, we explore the specifics of micellization in the solutions of anionic and cationic surfactants and their mixtures. Anionic surfactant sodium dodecyl sulfate (SDS) and cationic surfactant cetyltrimethylammonium bromide (CTAB) are chosen as characteristic examples. Coarse-grained models of the surfactants are constructed and parameterized using a combination of atomistic molecular simulation and infinite dilution activity coefficient calibration. Electrostatic interactions of charged beads are treated using a smeared charge approximation: the surfactant heads and dissociated counterions are modeled as beads with charges distributed around the bead center in an implicit dielectric medium. The proposed models semiquantitatively describe self-assembly in solutions of SDS and CTAB at various surfactant concentrations and molarities of added electrolyte. In particular, the model predicts a decline in the free surfactant concentration with the increase of the total surfactant loading, as well as characteristic aggregation transitions in single-component surfactant solutions caused by the addition of salt. The calculated values of the critical micelle concentration reasonably agree with experimental observations. Modeling of catanionic SDS-CTAB mixtures show consecutive transitions to worm-like micelles and then to vesicles caused by the addition of CTAB to micellar solution of SDS.
Collapse
Affiliation(s)
- Runfang Mao
- Department of Chemical Engineering, Rutgers, the State University of New Jersey , 98 Brett Road, Piscataway, New Jersey 08854, United States
| | - Ming-Tsung Lee
- Department of Chemical Engineering, Rutgers, the State University of New Jersey , 98 Brett Road, Piscataway, New Jersey 08854, United States
| | - Aleksey Vishnyakov
- Department of Chemical Engineering, Rutgers, the State University of New Jersey , 98 Brett Road, Piscataway, New Jersey 08854, United States
| | - Alexander V Neimark
- Department of Chemical Engineering, Rutgers, the State University of New Jersey , 98 Brett Road, Piscataway, New Jersey 08854, United States
| |
Collapse
|
8
|
Ravichandran G, Gopinath D. Ultrasonic relaxation studies on micelle formation in aqueous solutions of some bile salts. J Mol Liq 2014. [DOI: 10.1016/j.molliq.2014.06.014] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
9
|
Bhattacharjee JK, Kaatze U. Fluctuations near the critical micelle concentration. II. Ultrasonic attenuation spectra and scaling. J Phys Chem B 2013; 117:3798-805. [PMID: 23506227 DOI: 10.1021/jp401120x] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Based on a thermodynamic model of amphiphile solutions derived in the first part of the paper, the ultrasonic attenuation of such systems has been considered theoretically, including fluctuations of local concentrations and micelle sizes. At amphiphile concentrations smaller than the critical micelle concentration (cmc), scaling behavior in terms of the concentration distance to the cmc is predicted by theory, in fair agreement with experimental evidence. The scaling function in the sound attenuation below the cmc reveals the unsymmetric broadening in the spectra that clearly emerges from measurements when approaching the cmc. The shape of the scaling function corresponds to the experimental spectra of solutions with comparatively large cmc as well as with the relaxation spectral function of the unifying model of non-critical concentration fluctuations. Above the cmc, an additional relaxation term is predicted in correspondence with the Landau-Khalatnikov term in the sound attenuation of superfluid helium. This term is difficult to verify by measurements because, in the relevant frequency range, other processes may also contribute the ultrasonic attenuation spectra.
Collapse
|
10
|
Bhattacharjee JK, Kaatze U. Fluctuations Near the Critical Micelle Concentration. I. Premicellar Aggregation, Relaxation Rate, and Isentropic Compressibility. J Phys Chem B 2013; 117:3790-7. [DOI: 10.1021/jp4011185] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
| | - Udo Kaatze
- Drittes Physikalisches
Institut, Georg-August-Universität Göttingen, Friedrich-Hund-Platz
1, 37077 Göttingen, Germany
| |
Collapse
|
11
|
Tojo T, Fujisawa K, Muramatsu H, Hayashi T, Kim YA, Endo M, Terrones M, Dresselhaus MS. Controlled interlayer spacing of scrolled reduced graphene nanotubes by thermal annealing. RSC Adv 2013. [DOI: 10.1039/c3ra22976a] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
|
12
|
Sanders SA, Sammalkorpi M, Panagiotopoulos AZ. Atomistic Simulations of Micellization of Sodium Hexyl, Heptyl, Octyl, and Nonyl Sulfates. J Phys Chem B 2012; 116:2430-7. [DOI: 10.1021/jp209207p] [Citation(s) in RCA: 71] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Samantha A. Sanders
- Department of Chemical and Biological
Engineering and Institute for the Science and Technology of Materials, Princeton University, Princeton, New Jersey 08544,
United States
| | | | - Athanassios Z. Panagiotopoulos
- Department of Chemical and Biological
Engineering and Institute for the Science and Technology of Materials, Princeton University, Princeton, New Jersey 08544,
United States
| |
Collapse
|
13
|
Kaatze U. Kinetics of Micelle Formation and Concentration Fluctuations in Solutions of Short-Chain Surfactants. J Phys Chem B 2011; 115:10470-7. [DOI: 10.1021/jp205372g] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- U. Kaatze
- Drittes Physikalisches Institut, Georg-August-Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| |
Collapse
|
14
|
Bhattacharjee JK, Kaatze U. Anomalous Dynamics of Solutions of Nonionic Micelles in Water. J Phys Chem B 2011; 115:6069-75. [DOI: 10.1021/jp201015s] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
| | - Udo Kaatze
- Drittes Physikalisches Institut, Georg-August-Universität, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| |
Collapse
|
15
|
Polacek R, Kaatze U. Micelle and oligomer kinetics in aqueous solutions of n-octylammonium chloride: Monomer exchange, protrusion and chain isomerization. J Mol Liq 2010. [DOI: 10.1016/j.molliq.2010.09.001] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
16
|
Haller J, Kaatze U. Monomer Exchange and Rotational Isomerization of Alkyl Monoglycosides in Water. J Phys Chem B 2009; 113:12283-92. [DOI: 10.1021/jp905523p] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Julian Haller
- Drittes Physikalisches Institut, Georg-August-Universität, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| | - Udo Kaatze
- Drittes Physikalisches Institut, Georg-August-Universität, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| |
Collapse
|
17
|
Haller J, Kaatze U. Complexation versus micelle formation: α-Cyclodextrin+n-decyltrimethylammonium bromide aqueous solutions. Chem Phys Lett 2008. [DOI: 10.1016/j.cplett.2008.08.035] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|