1
|
Hendrikse RL, Amador C, Wilson MR. A many-body dissipative particle dynamics parametrisation scheme to study behaviour at air-water interfaces. SOFT MATTER 2023; 19:3590-3604. [PMID: 37161599 DOI: 10.1039/d3sm00276d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/11/2023]
Abstract
In this article, we present a general parametrisation scheme for many-body dissipative particle dynamics (MDPD). The scheme is based on matching model components to experimental surface tensions and chemical potentials. This allows us to obtain the correct surface and mixing behaviours of complex, multicomponent systems. The methodology is tested by modelling the behaviour of nonionic polyoxyethylene alkyl ether surfactants at an air/water interface. In particular, the influence of the number of ethylene oxide units in the surfactant head group is investigated. We find good agreement with many experimentally obtained parameters, such as minimum surface area per molecule; and a decrease in the surface tension with increasing surfactant surface density. Moreover, we observe an orientational transition, from surfactants lying directly on the water surface at low surface coverage, to surfactants lying parallel or tilted with respect to the surface normal at high surface coverage. The parametrisation scheme is also extended to cover the zwitterionic surfactant lauryldimethylamine oxide (LDAO), where we provide good predictions for the surface tension at maximum surface coverage. Here, if we exceed this coverage, we are able to demonstrate the spontaneous production of micelles from the surface surfactant layer.
Collapse
Affiliation(s)
| | - Carlos Amador
- Department of Chemistry, Durham University, Durham, DH1 3LE, UK.
| | - Mark R Wilson
- Department of Chemistry, Durham University, Durham, DH1 3LE, UK.
| |
Collapse
|
2
|
Anderson RL, Gunn DSD, Taddese T, Lavagnini E, Warren PB, Bray DJ. Phase Behavior of Alkyl Ethoxylate Surfactants in a Dissipative Particle Dynamics Model. J Phys Chem B 2023; 127:1674-1687. [PMID: 36786752 PMCID: PMC9969514 DOI: 10.1021/acs.jpcb.2c08834] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/15/2023]
Abstract
We present a dissipative particle dynamics (DPD) model capable of capturing the liquid state phase behavior of nonionic surfactants from the alkyl ethoxylate (CnEm) family. The model is based upon our recent work [Anderson et al. J. Chem. Phys. 2017, 147, 094503] but adopts tighter control of the molecular structure by setting the bond angles with guidance from molecular dynamics simulations. Changes to the geometry of the surfactants were shown to have little effect on the predicted micelle properties of sampled surfactants, or the water-octanol partition coefficients of small molecules, when compared to the original work. With these modifications the model is capable of reproducing the binary water-surfactant phase behavior of nine surfactants (C8E4, C8E5, C8E6, C10E4, C10E6, C10E8, C12E6, C12E8, and C12E12) with a good degree of accuracy.
Collapse
|
3
|
Minkara MS, Josephson TR, Venteicher CL, Greenvall BR, Lindsey RK, Koenig PH, Siepmann JI. Nonane and Hexanol Adsorption in the Lamellar Phase of a Nonionic Surfactant: Molecular Simulations and Comparison to Ideal Adsorbed Solution Theory. J Phys Chem B 2022; 126:3940-3949. [PMID: 35594369 DOI: 10.1021/acs.jpcb.2c02871] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Adsorption of n-nonane/1-hexanol (C9/C6OH) mixtures into the lamellar phase formed by a 50/50 w/w triethylene glycol mono-n-decyl ether (C10E3)/water system was studied using configurational-bias Monte Carlo simulations in the osmotic Gibbs ensemble. The interactions were described by the Shinoda-Devane-Klein coarse-grained force field. Prior simulations probing single-component adsorption indicated that C9 molecules preferentially load near the center of the bilayer, increasing the bilayer thickness, whereas C6OH molecules are more likely to be found near the interface of the polar and nonpolar moieties, swelling the bilayer in the lateral dimension. Here, we extend this work to binary C9/C6OH adsorption to probe whether the difference in the spatial preferences may lead to a synergistic effect and enhanced loadings for the mixture. Comparing loading trends and the thermodynamics of binary adsorption to unary adsorption reveals that C9-C9 interactions lead to the largest enhancement, whereas C9-C6OH and C6OH-C6OH interactions are less favorable for this bilayer system. Ideal adsorbed solution theory yields satisfactory predictions of the binary loading.
Collapse
Affiliation(s)
- Mona S Minkara
- Department of Bioengineering, Northeastern University, 360 Huntington Avenue, Boston, Massachusetts 02115, United States.,Department of Chemistry and Chemical Theory Center, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, United States
| | - Tyler R Josephson
- Department of Chemistry and Chemical Theory Center, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, United States.,Department of Chemical, Biochemical, and Environmental Engineering, University of Maryland, Baltimore County, 1000 Hilltop Circle, Baltimore, Maryland 21250, United States
| | - Connor L Venteicher
- Department of Chemistry and Chemical Theory Center, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, United States
| | - Benjamin R Greenvall
- Department of Bioengineering, Northeastern University, 360 Huntington Avenue, Boston, Massachusetts 02115, United States
| | - Rebecca K Lindsey
- Department of Chemistry and Chemical Theory Center, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, United States
| | - Peter H Koenig
- Computational Chemistry, Modeling and Simulation, The Procter and Gamble Company, 8256 Union Centre Blvd, West Chester, Ohio 45069, United States
| | - J Ilja Siepmann
- Department of Chemistry and Chemical Theory Center, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, United States.,Department of Chemical Engineering and Materials Science, University of Minnesota, 421 Washington Avenue SE, Minneapolis, Minnesota 55455-0132, United States
| |
Collapse
|
4
|
Guruge AG, Warren DB, Benameur H, Pouton CW, Chalmers DK. Aqueous phase behavior of the PEO-containing non-ionic surfactant C 12E 6: A molecular dynamics simulation study. J Colloid Interface Sci 2021; 588:257-268. [PMID: 33388586 DOI: 10.1016/j.jcis.2020.12.032] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2020] [Revised: 12/09/2020] [Accepted: 12/10/2020] [Indexed: 10/22/2022]
Abstract
HYPOTHESIS Non-ionic surfactants containing polyethylene oxide (PEO) chains are widely used in drug formulations, cosmetics, paints, textiles and detergents. High quality molecular dynamics models for PEO surfactants can give us detailed, atomic-scale information about the behavior of surfactant/water mixtures. SIMULATIONS We used two molecular dynamics force fields (FFs), 2016H66 and 53A6DBW, to model the simple non-ionic PEO surfactant, hexaoxyethylene dodecyl ether (C12E6). We investigated surfactant/water mixtures that span the phase diagram of starting from randomly distributed arrangements. In some cases, we also started with prebuilt, approximate models. The simulations results were compared with the experimentally observed phase behavior. FINDINGS Overall, this study shows that the spontaneous self-assembly of PEO non-ionic surfactants into different colloidal structures can be accurately modeled with MD simulations using the 2016H66 FF although transitions to well-formed hexagonal phase are slow. Of the two FFs investigated, the 2016H66 FF better reproduces the experimental phase behavior across all regions of the C12E6/water phase diagram.
Collapse
Affiliation(s)
- Amali G Guruge
- Medicinal Chemistry Monash Institute of Pharmaceutical Sciences, Monash University, 381 Royal Parade, Parkville, Victoria 3052, Australia
| | - Dallas B Warren
- Drug Delivery, Disposition and Dynamics, Monash Institute of Pharmaceutical Sciences, Monash University, 381 Royal Parade, Parkville, Victoria 3052, Australia
| | - Hassan Benameur
- Oral Drug Delivery Innovation, Chemical Division, Lonza, Strasbourg, France
| | - Colin W Pouton
- Drug Delivery, Disposition and Dynamics, Monash Institute of Pharmaceutical Sciences, Monash University, 381 Royal Parade, Parkville, Victoria 3052, Australia
| | - David K Chalmers
- Medicinal Chemistry Monash Institute of Pharmaceutical Sciences, Monash University, 381 Royal Parade, Parkville, Victoria 3052, Australia.
| |
Collapse
|
5
|
Talley Edwards A, Javidialesaadi A, Weigandt KM, Stan G, Eads CD. Structure and Dynamics of Spherical and Rodlike Alkyl Ethoxylate Surfactant Micelles Investigated Using NMR Relaxation and Atomistic Molecular Dynamics Simulations. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2019; 35:13880-13892. [PMID: 31573205 PMCID: PMC10552554 DOI: 10.1021/acs.langmuir.9b01345] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/22/2023]
Abstract
Predicting and controlling the properties of amphiphile aggregate mixtures require understanding the arrangements and dynamics of the constituent molecules. To explore these topics, we study molecular arrangements and dynamics in alkyl ethoxylate nonionic surfactant micelles by combining NMR relaxation measurements with large-scale atomistic molecular dynamics simulations. We calculate parameters that determine relaxation rates directly from simulated trajectories, without introducing specific functional forms to describe the dynamics. NMR relaxation rates, which depend on relative motions of interacting atom pairs, are influenced by wide distributions of dynamic time scales. We find that relative motions of neighboring atom pairs are rapid and liquidlike but are subject to structural constraints imposed by micelle morphology. Relative motions of distant atom pairs are slower than nearby atom pairs because changes in distances and angles are smaller when the moving atoms are further apart. Large numbers of atom pairs undergoing these slow relative motions contribute to predominantly negative cross-relaxation rates. For spherical micelles, but not for cylindrical micelles, cross-relaxation rates are positive only for surfactant tail atoms connected to the hydrophilic headgroup. This effect is related to the lower packing density of these atoms at the hydrophilic-hydrophobic boundary in spherical vs cylindrical arrangements, with correspondingly rapid and less constrained motion of atoms at the boundary.
Collapse
Affiliation(s)
- Allison Talley Edwards
- Department of Chemistry, University of Cincinnati, Cincinnati, Ohio 45221, United States
- Corporate Research & Development, The Procter & Gamble Company, Mason, Ohio 45040, United States
| | | | - Katie M. Weigandt
- National Institute of Standards and Technology, 100 Bureau Drive, MS 6102, Gaithersburg, Maryland 20899, United States
| | - George Stan
- Department of Chemistry, University of Cincinnati, Cincinnati, Ohio 45221, United States
| | - Charles D. Eads
- Corporate Research & Development, The Procter & Gamble Company, Mason, Ohio 45040, United States
| |
Collapse
|
6
|
Minkara MS, Lindsey RK, Noether CO, Venteicher CL, Jamadagni SN, Eike DM, Ghobadi AF, Koenig PH, Siepmann JI. Probing Additive Loading in the Lamellar Phase of a Nonionic Surfactant: Gibbs Ensemble Monte Carlo Simulations Using the SDK Force Field. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2018; 34:8245-8254. [PMID: 29902016 DOI: 10.1021/acs.langmuir.8b00687] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Understanding solute uptake into soft microstructured materials, such as bilayers and worm-like and spherical micelles, is of interest in the pharmaceutical, agricultural, and personal care industries. To obtain molecular-level insight on the effects of solutes loading into a lamellar phase, we utilize the Shinoda-Devane-Klein (SDK) coarse-grained force field in conjunction with configurational-bias Monte Carlo simulations in the osmotic Gibbs ensemble. The lamellar phase is comprised of a bilayer formed by triethylene glycol mono- n-decyl ether (C10E3) surfactants surrounded by water with a 50:50 surfactant/water weight ratio. We study both the unary adsorption isotherm and the effects on bilayer structure and stability caused by n-nonane, 1-hexanol, and ethyl butyrate at several different reduced reservoir pressures. The nonpolar n-nonane molecules load near the center of the bilayer. In contrast, the polar 1-hexanol and ethyl butyrate molecules both load with their polar bead close to the surfactant head groups. Near the center of the bilayer, none of the solute molecules exhibits a significant orientational preference. Solute molecules adsorbed near the polar groups of the surfactant chains show a preference for orientations perpendicular to the interface, and this alignment with the long axis of the surfactant molecules is most pronounced for 1-hexanol. Loading of n-nonane leads to an increase of the bilayer thickness, but does not affect the surface area per surfactant. Loading of polar additives leads to both lateral and transverse swelling. The reduced Henry's law constants of adsorption (expressed as a molar ratio of additive to surfactant per reduced pressure) are 0.23, 1.4, and 14 for n-nonane, 1-hexanol, and ethyl butyrate, respectively, and it appears that the SDK force field significantly overestimates the ethyl butyrate-surfactant interactions.
Collapse
Affiliation(s)
- Mona S Minkara
- Department of Chemistry and Chemical Theory Center , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455 , United States
| | - Rebecca K Lindsey
- Department of Chemistry and Chemical Theory Center , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455 , United States
| | - Celeste O. Noether
- Department of Chemistry and Chemical Theory Center , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455 , United States
| | - Connor L Venteicher
- Department of Chemistry and Chemical Theory Center , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455 , United States
| | - Sumanth N Jamadagni
- Computational Chemistry, Modeling and Simulation , The Procter & Gamble Company , 8256 Union Centre Blvd , West Chester , Ohio 45069 , United States
| | - David M Eike
- Computational Chemistry, Modeling and Simulation , The Procter & Gamble Company , 8256 Union Centre Blvd , West Chester , Ohio 45069 , United States
| | - Ahmad F Ghobadi
- Computational Chemistry, Modeling and Simulation , The Procter & Gamble Company , 8256 Union Centre Blvd , West Chester , Ohio 45069 , United States
| | - Peter H Koenig
- Computational Chemistry, Modeling and Simulation , The Procter & Gamble Company , 8256 Union Centre Blvd , West Chester , Ohio 45069 , United States
| | - J Ilja Siepmann
- Department of Chemistry and Chemical Theory Center , University of Minnesota , 207 Pleasant Street SE , Minneapolis , Minnesota 55455 , United States
- Department of Chemical Engineering and Materials Science , University of Minnesota , 421 Washington Avenue SE , Minneapolis , Minnesota 55455 , United States
| |
Collapse
|
7
|
Vu TV, Papavassiliou DV. Oil-water interfaces with surfactants: A systematic approach to determine coarse-grained model parameters. J Chem Phys 2018; 148:204704. [DOI: 10.1063/1.5022798] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Tuan V. Vu
- School of Chemical, Biological, and Material Engineering, University of Oklahoma, Norman, Oklahoma 73019-1004, USA
| | - Dimitrios V. Papavassiliou
- School of Chemical, Biological, and Material Engineering, University of Oklahoma, Norman, Oklahoma 73019-1004, USA
| |
Collapse
|
8
|
Vo MD, Shiau B, Harwell JH, Papavassiliou DV. Adsorption of anionic and non-ionic surfactants on carbon nanotubes in water with dissipative particle dynamics simulation. J Chem Phys 2016; 144:204701. [DOI: 10.1063/1.4949364] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
9
|
Rezaei H, Modarress H. Dissipative Particle Dynamics Study of Interfacial Properties and the Effects of Nonionic Surfactants on Hydrocarbon/Water Microemulsions. J DISPER SCI TECHNOL 2015. [DOI: 10.1080/01932691.2015.1077453] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
10
|
Suttipong M, Grady BP, Striolo A. Self-assembled surfactants on patterned surfaces: confinement and cooperative effects on aggregate morphology. Phys Chem Chem Phys 2014; 16:16388-98. [DOI: 10.1039/c4cp00882k] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
|
11
|
Ferreira TM, Topgaard D, Ollila OHS. Molecular conformation and bilayer pores in a nonionic surfactant lamellar phase studied with 1H-13C solid-state NMR and molecular dynamics simulations. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2014; 30:461-469. [PMID: 24372090 DOI: 10.1021/la404684r] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/03/2023]
Abstract
The structure of the lamellar phase of aqueous pentaethylene glycol mono-n-dodecyl ether (C12E5) surfactant at various temperatures and molar fractions is studied by using united atom molecular dynamics simulations and nuclear magnetic resonance measurements. Namely, the simulation model is used to interpret the magnitude and temperature dependence of experimental C-H order parameter profiles in terms of the molecular conformation and orientation. Our simulations suggest that the low order parameters that are generally measured in poly(ethylene oxide) surfactant bilayers are due to the presence of bilayer pores throughout the entire lamellar phase region.
Collapse
Affiliation(s)
- Tiago M Ferreira
- Division of Physical Chemistry, Centre for Chemistry and Chemical Engineering, Lund University , Lund, Sweden
| | | | | |
Collapse
|
12
|
Levine BG, LeBard DN, DeVane R, Shinoda W, Kohlmeyer A, Klein ML. Micellization Studied by GPU-Accelerated Coarse-Grained Molecular Dynamics. J Chem Theory Comput 2011; 7:4135-45. [DOI: 10.1021/ct2005193] [Citation(s) in RCA: 55] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Affiliation(s)
- Benjamin G. Levine
- Institute for Computational Molecular Science and Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - David N. LeBard
- Institute for Computational Molecular Science and Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Russell DeVane
- The Procter & Gamble Company, Cincinnati, Ohio 45253, United States
| | - Wataru Shinoda
- National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba, Ibaraki, Japan
| | - Axel Kohlmeyer
- Institute for Computational Molecular Science and Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Michael L. Klein
- Institute for Computational Molecular Science and Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| |
Collapse
|
13
|
Rodgers TL, Mihailova O, Siperstein FR. Dissolution of Lamellar Phases. J Phys Chem B 2011; 115:10218-27. [DOI: 10.1021/jp111464b] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Affiliation(s)
- Thomas L. Rodgers
- SCEAS, The University of Manchester, Manchester, United Kingdom, M13 9PL
| | - Olga Mihailova
- SCEAS, The University of Manchester, Manchester, United Kingdom, M13 9PL
| | - Flor R. Siperstein
- SCEAS, The University of Manchester, Manchester, United Kingdom, M13 9PL
| |
Collapse
|