1
|
Williams AP, Faber JM, Recsei C, de Campo L, Darwish TA, Tuck KL, Dagastine RR, Tabor RF. Salt-Induced Linker Dehydration Modulates Micellar Structure in Ether-Linked Sulfate Surfactants. J Phys Chem B 2024; 128:6648-6653. [PMID: 38935971 DOI: 10.1021/acs.jpcb.4c03429] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/29/2024]
Abstract
Ether-linked surfactants are widely used in formulations such as liquid soaps, but despite their ubiquity, it is unclear how n-ethylene glycol linkers in surfactants, such as sodium lauryl n-(ethylene glycol) sulfate (SLEnS), influence micellar packing in the presence of NaCl. In the present work, we probe the structure and hydration of ether linkers in micelles comprising monodisperse SLEnS surfactants using contrast-variation small-angle neutron scattering (CV-SANS) and small-angle X-ray scattering (SAXS). Using SAXS, changes in micellar structure were observed for SLEnS (n = 1, 2, or 3) arising from the extent of ethoxylation. Scattering profiles indicated a clear transition from elongated cylindrical micelles to shorter ellipsoidal micelles with increasing ethoxylation. With CV-SANS, micellar structure and linker geometries of SLE3S were able to be resolved, indicating that a change in micellar architecture is modulated by dehydration of the tri(ethylene glycol) linker, offering new insights into the role of water and ions in the self-assembly of this key class of surfactants.
Collapse
Affiliation(s)
- Ashley P Williams
- School of Chemistry, Monash University, Clayton, Victoria 3800, Australia
| | - Jonathan M Faber
- School of Chemistry, Monash University, Clayton, Victoria 3800, Australia
| | - Carl Recsei
- National Deuteration Facility, ANSTO, Lucas Heights, New South Wales 2234, Australia
| | - Liliana de Campo
- Australian Centre for Neutron Scattering, ANSTO, Lucas Heights, New South Wales 2234, Australia
| | - Tamim A Darwish
- National Deuteration Facility, ANSTO, Lucas Heights, New South Wales 2234, Australia
- Faculty of Science and Technology, University of Canberra, Bruce, Australian Capital Territory 2617, Australia
| | - Kellie L Tuck
- School of Chemistry, Monash University, Clayton, Victoria 3800, Australia
| | - Raymond R Dagastine
- Department of Chemical and Biological Engineering, University of Melbourne, Parkville 3052, Australia
| | - Rico F Tabor
- School of Chemistry, Monash University, Clayton, Victoria 3800, Australia
| |
Collapse
|
2
|
Li J, Amador C, Wilson MR. Computational predictions of interfacial tension, surface tension, and surfactant adsorption isotherms. Phys Chem Chem Phys 2024; 26:12107-12120. [PMID: 38587476 DOI: 10.1039/d3cp06170a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/09/2024]
Abstract
All-atom (AA) molecular dynamics (MD) simulations are employed to predict interfacial tensions (IFT) and surface tensions (ST) of both ionic and non-ionic surfactants. The general AMBER force field (GAFF) and variants are examined in terms of their performance in predicting accurate IFT/ST, γ, values for chosen water models, together with the hydration free energy, ΔGhyd, and density, ρ, predictions for organic bulk phases. A strong correlation is observed between the quality of ρ and γ predictions. Based on the results, the GAFF-LIPID force field, which provides improved ρ predictions is selected for simulating surfactant tail groups. Good γ predictions are obtained with GAFF/GAFF-LIPID parameters and the TIP3P water model for IFT simulations at a water-triolein interface, and for GAFF/GAFF-LIPID parameters together with the OPC4 water model for ST simulations at a water-vacuum interface. Using a combined molecular dynamics-molecular thermodynamics theory (MD-MTT) framework, a mole fraction of C12E6 molecule of 1.477 × 10-6 (from the experimental critical micelle concentration, CMC) gives a simulated surface excess concentration, ΓMAX, of 76 C12E6 molecules at a 36 nm2 water-vacuum surface (3.5 × 10-10 mol cm-2), which corresponds to a simulated ST of 35 mN m-1. The results compare favourably with an experimental ΓMAX of C12E6 of 3.7 × 10-10 mol cm-2 (80 surfactants for a 36 nm2 surface) and experimental ST of C12E6 of 32 mN m-1 at the CMC.
Collapse
Affiliation(s)
- Jing Li
- Department of Chemistry, Durham University, Stockton Road, Durham, DH1 3LE, UK.
| | - Carlos Amador
- Newcastle Innovation Centre, Procter & Gamble Ltd, Newcastle Upon Tyne, NE12 9BZ, UK
| | - Mark R Wilson
- Department of Chemistry, Durham University, Stockton Road, Durham, DH1 3LE, UK.
| |
Collapse
|
3
|
Shen K, Nguyen M, Sherck N, Yoo B, Köhler S, Speros J, Delaney KT, Shell MS, Fredrickson GH. Predicting surfactant phase behavior with a molecularly informed field theory. J Colloid Interface Sci 2023; 638:84-98. [PMID: 36736121 DOI: 10.1016/j.jcis.2023.01.015] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/24/2022] [Revised: 12/24/2022] [Accepted: 01/03/2023] [Indexed: 01/09/2023]
Abstract
HYPOTHESIS The computational study of surfactants and self-assembly is challenging because 1) models need to reflect chemistry-specific interactions, and 2) self-assembled structures are difficult to equilibrate with conventional molecular dynamics. We propose to overcome these challenges with a multiscale simulation approach where relative entropy minimization transfers chemically-detailed information from all-atom (AA) simulations to coarse-grained (CG) models that can be simulated using field-theoretic methods. Field-theoretic simulations are not limited by intrinsic physical time scales like diffusion and allow for rigorous equilibration via free energy minimization. This approach should enable the study of properties that are difficult to obtain by particle-based simulations. SIMULATION WORK We apply this workflow to sodium dodecylsulfate. To ensure chemical fidelity we present an AA force field calibrated against interfacial tension experiments. We generate CG models from AA simulation trajectories and show that particle-based and field-theoretic simulations of the CG model reproduce AA simulations and experimental measurements. FINDINGS The workflow captures the complex balance of interactions in a multicomponent system ultimately described by an atomistic model. The resulting CG models can study complex 3D phases like double or alternating gyroids, and reproduce salt effects on properties like aggregation number and shape transitions.
Collapse
Affiliation(s)
- Kevin Shen
- Department of Chemical Engineering, University of California, Santa Barbara, Santa Barbara 93106, CA, United States; Materials Research Laboratory, University of California, Santa Barbara, Santa Barbara 93106, CA, United States.
| | - My Nguyen
- Department of Chemical Engineering, University of California, Santa Barbara, Santa Barbara 93106, CA, United States
| | - Nicholas Sherck
- Department of Chemical Engineering, University of California, Santa Barbara, Santa Barbara 93106, CA, United States
| | - Brian Yoo
- BASF Corporation, Tarrytown 10591, NY, United States
| | | | - Joshua Speros
- California Research Alliance (CARA) by BASF, Berkeley 94720, CA, United States
| | - Kris T Delaney
- Materials Research Laboratory, University of California, Santa Barbara, Santa Barbara 93106, CA, United States
| | - M Scott Shell
- Department of Chemical Engineering, University of California, Santa Barbara, Santa Barbara 93106, CA, United States.
| | - Glenn H Fredrickson
- Department of Chemical Engineering, University of California, Santa Barbara, Santa Barbara 93106, CA, United States; Materials Research Laboratory, University of California, Santa Barbara, Santa Barbara 93106, CA, United States; Department of Materials Engineering, University of California, Santa Barbara, Santa Barbara 93106, CA, United States.
| |
Collapse
|
4
|
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
|
5
|
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
|
6
|
Navarro-Aquino D, Medeiros M. Molecular simulation of pentaethylene glycol monooctyl ether micelles in water. Prediction of the micellization Gibbs energy. Colloids Surf A Physicochem Eng Asp 2022. [DOI: 10.1016/j.colsurfa.2022.128406] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/03/2022]
|
7
|
Inverse ISAsomes in Bio-Compatible Oils—Exploring Formulations in Squalane, Triolein and Olive Oil. NANOMATERIALS 2022; 12:nano12071133. [PMID: 35407249 PMCID: PMC9000821 DOI: 10.3390/nano12071133] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/12/2022] [Revised: 03/26/2022] [Accepted: 03/28/2022] [Indexed: 11/17/2022]
Abstract
In contrast to their more common counterparts in aqueous solutions, inverse ISAsomes (internally self-assembled somes/particles) are formulated as kinetically stabilised dispersions of hydrophilic, lyotropic liquid-crystalline (LC) phases in non-polar oils. This contribution reports on their formation in bio-compatible oils. We found that it is possible to create inverse hexosomes, inverse micellar cubosomes (Fd3m) and an inverse emulsified microemulsion (EME) in excess squalane with a polyethylene glycol alkyl ether as the primary surfactant forming the LC phase and to stabilise them with hydrophobised silica nanoparticles. Furthermore, an emulsified -phase and inverse hexosomes were formed in excess triolein with the triblock-copolymer Pluronic® P94 as the primary surfactant. Stabilisation was achieved with a molecular stabiliser of type polyethylene glycol (PEG)-dipolyhydroxystearate. For the inverse hexosomes in triolein, the possibility of a formulation without any additional stabiliser was explored. It was found that a sufficiently strong stabilisation effect was created by the primary surfactant alone. Finally, triolein was replaced with olive oil which also led to the successful formation of inverse hexosomes. As far as we know, there exists no previous contribution about inverse ISAsomes in complex oils such as triolein or plant oils, and the existence of stabiliser-free (i.e., self-stabilising) inverse hexosomes has also not been reported until now.
Collapse
|
8
|
Kroll P, Benke J, Enders S, Brandenbusch C, Sadowski G. Influence of Temperature and Concentration on the Self-Assembly of Nonionic C iE j Surfactants: A Light Scattering Study. ACS OMEGA 2022; 7:7057-7065. [PMID: 35252696 PMCID: PMC8892478 DOI: 10.1021/acsomega.1c06766] [Citation(s) in RCA: 16] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/30/2021] [Accepted: 02/02/2022] [Indexed: 06/07/2023]
Abstract
Nonionic poly(ethylene oxide) alkyl ether (CiEj) surfactants self-assemble into aggregates of various sizes and shapes above their critical micelle concentration (CMC). Knowledge on solution attributes such as CMC as well as aggregate characteristics is crucial to choose the appropriate surfactant for a given application, e.g., as a micellar solvent system. In this work, we used static and dynamic light scattering to measure the CMC, aggregation number (N agg), and hydrodynamic radius (R h) of four different CiEj surfactants (C8E5, C8E6, C10E6, and C10E8). We examined the influence of temperature, concentration, and molecular structure on the self-assembly in the vicinity of the CMC. A minimum in the CMC vs temperature curve was identified for all surfactants investigated. Further, extending the hydrophilic and hydrophobic chain lengths leads to an increase and decrease of the CMC, respectively. The size of the aggregates strongly depends on temperature. N agg and R h increase with increasing temperature for all surfactants investigated. Additionally, N agg and R h both increase with increasing surfactant concentration. The data obtained in this work further improve the understanding of the influence of temperature and molecular structure on the self-assembly of CiEj surfactants and will further foster their use in micellar solvent systems.
Collapse
Affiliation(s)
- Peter Kroll
- Laboratory
of Thermodynamics, Department of Biochemical and Chemical Engineering, TU Dortmund University, Emil-Figge-Str. 70, 44227 Dortmund, Germany
| | - Julius Benke
- Laboratory
of Thermodynamics, Department of Biochemical and Chemical Engineering, TU Dortmund University, Emil-Figge-Str. 70, 44227 Dortmund, Germany
| | - Sabine Enders
- Institute
for Technical Thermodynamics and Refrigeration, Karlsruhe Institute of Technology, Engler-Bunte-Ring 21, 76131 Karlsruhe, Germany
| | - Christoph Brandenbusch
- Laboratory
of Thermodynamics, Department of Biochemical and Chemical Engineering, TU Dortmund University, Emil-Figge-Str. 70, 44227 Dortmund, Germany
| | - Gabriele Sadowski
- Laboratory
of Thermodynamics, Department of Biochemical and Chemical Engineering, TU Dortmund University, Emil-Figge-Str. 70, 44227 Dortmund, Germany
| |
Collapse
|
9
|
Affiliation(s)
- Nilmoni Sarkar
- Department of Chemistry, Indian Institute of Technology, Kharagpur 721302, WB India
| |
Collapse
|
10
|
Lee MT. Micellization of Rhamnolipid Biosurfactants and Their Applications in Oil Recovery: Insights from Mesoscale Simulations. J Phys Chem B 2021; 125:9895-9909. [PMID: 34423979 DOI: 10.1021/acs.jpcb.1c05802] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
The dissipative particle dynamics (DPD) mesoscopic method is used to investigate the self-assembly of rhamnolipid congeners and their aggregation behaviors with paraffins including nonane and pentadecane. The coarse-grained force field is parameterized by combining molecular dynamics (MD) simulations, COSMOtherm calculations, and available experimental data. This model reproduces the vesicular formation of α-l-rhamnopyranosyl-β-hydroxydecanoyl-β-hydroxydecanoate (Rha-C10-C10) reported by all-atom MD simulations. The vesicle composed of Rha-C10-C10 is found to be most stable at a surfactant concentration of 100-146 mM based on asphericity analysis. The architecture of rhamnolipid congeners affects the morphology of their aggregates. Di-rhamno-di-lipidic dRha-C16-C16 forms vesicles with a thicker unilamellar layer of 3.2 nm. Rha-C16-C16 forms vesicles at a lower concentration of 70 mM, but the enclosed water space collapses when the surfactant concentration increases. dRha-C10-C10 forms wormlike micelles, which agglomerate into a torus and interconnected network at higher concentrations. In the presence of alkane molecules, dRha-C10-C10 maintains its wormlike micellar morphology with alkane molecules wrapped inside the aggregates. For Rha-C10-C10, Rha-C16-C16, and dRha-C16-C16, nonane molecules are distributed in the hydrophobic subdomain formed by rhamnolipid molecules. Spherical vesicles are formed at a surfactant concentration of 50 mM and then develop into ellipsoidal vesicles when the concentration increases to 125 mM. When mixed with pentadecane, the alkane molecules are aggregated and surrounded by surfactants forming a core-shell structure at a low surfactant concentration of 20 mM. At higher alkane and surfactant concentrations, the morphologies develop into disk micelles, wormlike micelles, and vesicles, with pentadecane molecules being distributed and packed with rhamnolipids. The obtained simulation results suggest that these biosurfactants have potential as environmental remediation agents.
Collapse
Affiliation(s)
- Ming-Tsung Lee
- Department of Chemical Engineering and Biotechnology, National Taipei University of Technology, Taipei 10608, Taiwan
| |
Collapse
|
11
|
Moulik SP, Rakshit AK, Naskar B. Evaluation of Non‐Ambiguous Critical Micelle Concentration of Surfactants in Relation to Solution Behaviors of Pure and Mixed Surfactant Systems: A Physicochemical Documentary and Analysis. J SURFACTANTS DETERG 2021. [DOI: 10.1002/jsde.12503] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- Satya P. Moulik
- Centre for Surface Science, Department of Chemistry Jadavpur University Kolkata 700032 India
| | - Animesh K. Rakshit
- Indian Society for Surface Science & Technology, Department of Chemistry Jadavpur University Kolkata 700032 India
| | - Bappaditya Naskar
- Department of Chemistry, Sundarban Hazi Desarat College University of Calcutta Pathankhali 743611 India
| |
Collapse
|
12
|
Synthesis and application of geminal dicationic ionic liquids and poly (ionic liquids) combined imidazolium and pyridinium cations as demulsifiers for petroleum crude oil saline water emulsions. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2020.115264] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
13
|
Šindelka K, Lísal M. Interplay between surfactant self-assembly and adsorption at hydrophobic surfaces: insights from dissipative particle dynamics. Mol Phys 2020. [DOI: 10.1080/00268976.2020.1857863] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
Affiliation(s)
- Karel Šindelka
- Department of Molecular and Mesoscopic Modelling, The Czech Academy of Sciences, Institute of Chemical Process Fundamentals, Prague, Czech Republic
| | - Martin Lísal
- Department of Molecular and Mesoscopic Modelling, The Czech Academy of Sciences, Institute of Chemical Process Fundamentals, Prague, Czech Republic
- Department of Physics, Faculty of Science, Jan Evangelista Purkyně University in Ústí nad Labem, Ústí nad Labem, Czech Republic
| |
Collapse
|
14
|
Johnston MA, Duff AI, Anderson RL, Swope WC. Model for the Simulation of the C nE m Nonionic Surfactant Family Derived from Recent Experimental Results. J Phys Chem B 2020; 124:9701-9721. [PMID: 32986421 DOI: 10.1021/acs.jpcb.0c06132] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Using a comprehensive set of recently published experimental results for training and validation, we have developed computational models appropriate for simulations of aqueous solutions of poly(ethylene oxide) alkyl ethers, an important class of micelle-forming nonionic surfactants, usually denoted CnEm. These models are suitable for use in simulations that employ a moderate amount of coarse graining and especially for dissipative particle dynamics (DPD), which we adopt in this work. The experimental data used for training and validation were reported earlier and produced in our laboratory using dynamic light scattering (DLS) measurements performed on 12 members of the CnEm compound family yielding micelle size distribution functions and mass-weighted mean aggregation numbers at each of several surfactant concentrations. The range of compounds and quality of the experimental results were designed to support the development of computational models. An essential feature of this work is that all simulation results were analyzed in a way that is consistent with the experimental data. Proper account is taken of the fact that a broad distribution of micelle sizes exists, so mass-weighted averages (rather than number-weighted averages) over this distribution are required for the proper comparison of simulation and experimental results. The resulting DPD force field reproduces several important trends seen in the experimental critical micelle concentrations and mass-averaged mean aggregation numbers with respect to surfactant characteristics and concentration. We feel it can be used to investigate a number of open questions regarding micelle sizes and shapes and their dependence on surfactant concentration for this important class of nonionic surfactants.
Collapse
Affiliation(s)
| | - Andrew Ian Duff
- STFC Hartree Centre, SciTech Daresbury, Warrington, Cheshire WA4 4AD, U.K
| | - Richard L Anderson
- STFC Hartree Centre, SciTech Daresbury, Warrington, Cheshire WA4 4AD, U.K
| | - William C Swope
- IBM Almaden Research Center, 650 Harry Road, San Jose, California 95120, United States
| |
Collapse
|
15
|
|
16
|
McDonagh JL, Swope WC, Anderson RL, Johnston MA, Bray DJ. What can digitisation do for formulated product innovation and development? POLYM INT 2020. [DOI: 10.1002/pi.6056] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
| | | | | | | | - David J Bray
- The Hartree Centre STFC Daresbury Laboratory Warrington WA4 4AD UK
| |
Collapse
|
17
|
Lavagnini E, Cook JL, Warren PB, Williamson MJ, Hunter CA. A Surface Site Interaction Point Method for Dissipative Particle Dynamics Parametrization: Application to Alkyl Ethoxylate Surfactant Self-Assembly. J Phys Chem B 2020; 124:5047-5055. [PMID: 32510951 PMCID: PMC7309324 DOI: 10.1021/acs.jpcb.0c01895] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
![]()
Dissipative
particle dynamics (DPD) is a coarse-grained approach
to the simulation of large supramolecular systems, but one limitation
has been that the parameters required to describe the noncovalent
interactions between beads are not readily accessible. A first-principles
computational method has been developed so that bead interaction parameters
can be calculated directly from ab initio gas-phase
molecular electrostatic potential surfaces of the molecular fragments
that represent the beads. A footprinting algorithm converts the molecular
electrostatic potential surfaces into a discrete set of surface site
interaction points (SSIPs), and these SSIPs are used in the SSIMPLE
(surface site interaction model for the properties of liquids at equilibrium)
algorithm to calculate the free energies of transfer of one bead into
a solution of any other bead. The bead transfer free energies are
then converted into the required DPD interaction parameters for all
pairwise combinations of different beads. The reliability of the parameters
was demonstrated using DPD simulations of a range of alkyl ethoxylate
surfactants. The simulations reproduce the experimentally determined
values of the critical micelle concentration and mean aggregation
number well for all 22 surfactants studied.
Collapse
Affiliation(s)
- Ennio Lavagnini
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, U.K
| | - Joanne L Cook
- Unilever R&D Port Sunlight, Quarry Road East, Bebington CH63 3JW, U.K
| | - Patrick B Warren
- Unilever R&D Port Sunlight, Quarry Road East, Bebington CH63 3JW, U.K.,The Hartree Centre, STFC Daresbury Laboratory, Warrington WA4 4AD, U.K
| | - Mark J Williamson
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, U.K
| | - Christopher A Hunter
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, U.K
| |
Collapse
|
18
|
Constructing the phase diagram of sodium laurylethoxysulfate using dissipative particle dynamics. J Colloid Interface Sci 2019; 557:34-44. [DOI: 10.1016/j.jcis.2019.08.091] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2019] [Revised: 08/23/2019] [Accepted: 08/24/2019] [Indexed: 01/08/2023]
|
19
|
McDonagh JL, Shkurti A, Bray DJ, Anderson RL, Pyzer-Knapp EO. Utilizing Machine Learning for Efficient Parameterization of Coarse Grained Molecular Force Fields. J Chem Inf Model 2019; 59:4278-4288. [PMID: 31549507 DOI: 10.1021/acs.jcim.9b00646] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
We present a machine learning approach to automated force field development in dissipative particle dynamics (DPD). The approach employs Bayesian optimization to parametrize a DPD force field against experimentally determined partition coefficients. The optimization process covers a discrete space of over 40 000 000 points, where each point represents the set of potentials that jointly forms a force field. We find that Bayesian optimization is capable of reaching a force field of comparable performance to the current state-of-the-art within 40 iterations. The best iteration during the optimization achieves an R2 of 0.78 and an RMSE of 0.63 log units on the training set of data, these metrics are maintained when a validation set is included, giving R2 of 0.8 and an RMSE of 0.65 log units. This work hence provides a proof-of-concept, expounding the utility of coupling automated and efficient global optimization with a top down data driven approach to force field parametrization. Compared to commonly employed alternative methods, Bayesian optimization offers global parameter searching and a low time to solution.
Collapse
Affiliation(s)
- James L McDonagh
- IBM Research U.K. , Hartree Centre, Daresbury WA4 4AD , United Kingdom
| | - Ardita Shkurti
- STFC Daresbury Laboratories , Daresbury WA4 4AD , United Kingdom
| | - David J Bray
- STFC Daresbury Laboratories , Daresbury WA4 4AD , United Kingdom
| | | | | |
Collapse
|
20
|
Wu X, Chen L, Fan Y, Fu F, Li J, Zhang J. Water Solubility and Surface Property of Alkyl Di-/Tri-/Tetraoxyethyl β-d-Xylopyranosides. JOURNAL OF AGRICULTURAL AND FOOD CHEMISTRY 2019; 67:10361-10372. [PMID: 31487173 DOI: 10.1021/acs.jafc.9b03435] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Alkyl di-/tri-/tetraoxyethyl β-d-xylopyranosides as derivatives of alkyl xylosides are a class of non-ionic sugar-based surfactants. They were stereoselectively synthesized by the Helferich method. Their properties including hydrophilic-lipophilic balance number, water solubility, surface property, foam property, emulsifying property, and thermotropic liquid crystal property were mainly investigated. The results showed that their water solubility decreased with increasing the alkyl chain length and increasing the number of the oligooxyethyl fragment. The critical micelle concentration had a monotonous decreasing trend with increasing the alkyl chain length. Nonyl di-/tri-/tetraoxyethyl β-d-xylopyranosides [-(OCH2CH2)m-, where m = 2, 3, and 4] exhibited the most excellent foaming ability and foam stability. In the n-octane/water system, dodecyl tetraoxyethyl β-d-xylopyranosides and tetradecyl tetraoxyethyl β-d-xylopyranosides had the strongest emulsion ability. In addition, some alkyl di-/tri-/tetraoxyethyl β-d-xylopyranosides had thermotropic liquid crystal properties. Such sugar-based surfactants, alkyl di-/tri-/tetraoxyethyl β-d-xylopyranosides, will be expected to develop for a variety of practical application.
Collapse
Affiliation(s)
- Xiubing Wu
- Key Laboratory of Environmentally Friendly Chemistry and Applications of Ministry of Education, College of Chemistry , Xiangtan University , Xiangtan , Hunan 411105 , People's Republic of China
| | - Langqiu Chen
- Key Laboratory of Environmentally Friendly Chemistry and Applications of Ministry of Education, College of Chemistry , Xiangtan University , Xiangtan , Hunan 411105 , People's Republic of China
| | - Yulin Fan
- Key Laboratory of Environmentally Friendly Chemistry and Applications of Ministry of Education, College of Chemistry , Xiangtan University , Xiangtan , Hunan 411105 , People's Republic of China
| | - Fang Fu
- Key Laboratory of Environmentally Friendly Chemistry and Applications of Ministry of Education, College of Chemistry , Xiangtan University , Xiangtan , Hunan 411105 , People's Republic of China
| | - Jiping Li
- Key Laboratory of Environmentally Friendly Chemistry and Applications of Ministry of Education, College of Chemistry , Xiangtan University , Xiangtan , Hunan 411105 , People's Republic of China
| | - Jing Zhang
- Key Laboratory of Environmentally Friendly Chemistry and Applications of Ministry of Education, College of Chemistry , Xiangtan University , Xiangtan , Hunan 411105 , People's Republic of China
| |
Collapse
|