1
|
Hendrikse RL, Amador C, Wilson MR. Many-body dissipative particle dynamics simulations of micellization of sodium alkyl sulfates. SOFT MATTER 2024. [PMID: 39034768 DOI: 10.1039/d4sm00533c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/23/2024]
Abstract
We present a study of micelle formation in alkyl sulfate surfactants using the simulation method of many-body dissipative particle dynamics (MDPD). We parametrise our model by tuning the intermolecular interactions in order to reproduce experimental values for the chemical potential and density at room temperature. Using this approach, we find that our model shows good agreement with experimental values for the critical micelle concentration (CMC). Furthermore, we show that our model can accurately predict CMC trends, which result from varying properties such as surfactant tail length and the salt concentration. We apply our model to investigate the effect of aggregation number on various micellar properties, such as the shape of individual micelles and the fraction of bound counterions. We show that micelles become aspherical at large aggregation numbers, in line with experimental predictions, and that longer tail surfactants are generally more spherical at all aggregation numbers compared to those which are shorter. We find excellent agreement between our simulations and experimental values for the degree of counterion binding, a factor that is crucial to accurately studying micellar shape, but one that is typically overlooked in the existing literature.
Collapse
Affiliation(s)
| | - Carlos Amador
- Procter and Gamble, Newcastle Innovation Centre, Whitley Road, Newcastle upon Tyne, NE12 9BZ, UK
| | - Mark R Wilson
- Department of Chemistry, Durham University, Durham, DH1 3LE, UK.
| |
Collapse
|
2
|
Zhang B, Guan B, Tao Y, Liu W, Peng B, Lv K. Surfactant Synergistic Effect and Interfacial Properties of Microemulsions Compounded with Anionic and Nonionic Surfactants Using Dissipative Particle Dynamics. ACS OMEGA 2024; 9:23903-23916. [PMID: 38854575 PMCID: PMC11154924 DOI: 10.1021/acsomega.4c01933] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/01/2024] [Revised: 05/08/2024] [Accepted: 05/13/2024] [Indexed: 06/11/2024]
Abstract
Microemulsions are one of the most promising directions in enhanced oil recovery, but conventional screening methods are time-consuming and labor-intensive and lack the means to analyze them at the microscopic level. In this paper, we used the Clint model to predict the changes in the synergistic effect of the mixed system of anionic surfactant sodium dodecyl benzenesulfonate and nonionic surfactant polyethoxylated fatty alcohols (C12E6), generated microemulsions using surfactant systems with different mole fractions, and used particle size to analyze the performance and stability of microemulsions, analyze the properties and stability of microemulsions using particle size, and analyze the interfacial behaviors and changes of microemulsions when different systems constitute microemulsions from the point of view of mesoscopic microemulsion self-assembly behaviors by combining with dissipative particle dynamics. It has been shown that microemulsion systems generated from anionic and nonanionic surfactants with a synergistic effect, based on the Clint model, exhibit excellent performance and stability at the microscopic level. The method proposed in this paper can dramatically improve the screening efficiency of microemulsions of anionic and nonanionic surfactants and accurately analyze the properties of microemulsions, so as to provide a theoretical basis for the subsequent research on microemulsions.
Collapse
Affiliation(s)
- Biao Zhang
- Research
Institute of Percolation Fluids Mechanics, Chinese Academy of Sciences, Beijing 100010, China
| | - Baoshan Guan
- Research
Institute of Percolation Fluids Mechanics, Chinese Academy of Sciences, Beijing 100010, China
- Research
Institute of Petroleum Exploration and Development, Beijing 100010, China
| | - Yufan Tao
- Research
Institute of Percolation Fluids Mechanics, Chinese Academy of Sciences, Beijing 100010, China
| | - Weidong Liu
- Research
Institute of Percolation Fluids Mechanics, Chinese Academy of Sciences, Beijing 100010, China
- Research
Institute of Petroleum Exploration and Development, Beijing 100010, China
| | - Baoliang Peng
- Research
Institute of Petroleum Exploration and Development, Beijing 100010, China
| | - Kai Lv
- Research
Institute of Percolation Fluids Mechanics, Chinese Academy of Sciences, Beijing 100010, China
| |
Collapse
|
3
|
Dorenbos G. How fork-length asymmetry affects solvent connectivity and diffusion in grafted polymeric model membranes. J Chem Phys 2024; 160:064901. [PMID: 38341779 DOI: 10.1063/5.0193120] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2023] [Accepted: 01/15/2024] [Indexed: 02/13/2024] Open
Abstract
The hydrophilic pore morphology and solvent diffusion within model (amphiphilic) polymer membranes are simulated by dissipative particle dynamics (DPD). The polymers are composed of a backbone of 18 covalently bonded A beads to which at regular intervals side chains are attached. The side chains are composed of linear Ap chains (i.e., -A1-A2…Ap) from which two branches, [AsC] and [ArC], split off (s ≤ r). C beads serve as functionalized hydrophilic pendent sites. The branch lengths (s + 1 and r + 1) are varied. Five repeat unit designs (with general formula A3[Ap[AsC][ArC]]) are considered: A2[A3C][A3C] (symmetric branching), A2[A2C][A4C], A2[AC][A5C], A2[C][A6C] (highly asymmetric branching), and A4[AC][A3C]. The distribution of water (W) and W diffusion through nanophase segregated hydrophilic pores is studied. For similar primary length p, an increase in side chain symmetry favors hydrophilic pore connectivity and long-range water transport. C beads located on the longer [ArC] branches reveal the highest C bead mobility and are more strongly associated with water than the C beads on the shorter [AsC] branches. The connectivity of hydrophilic (W and W + C) phases through mapped replica of selected snapshots obtained from Monte Carlo tracer diffusion simulations is in line with trends found from the W bead diffusivities during DPD simulations. The diffusive pathways for protons (H+) in proton exchange membranes and for hydronium (OH-) in anion exchange membranes are the same as for solvents. Therefore, control of the side chain architecture is an interesting design parameter for optimizing membrane conductivities.
Collapse
Affiliation(s)
- G Dorenbos
- Private research, Sano 1107-2, Belle Crea 502, 410-1118 Susono, Japan
| |
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
|
Palkar V, Thakar D, Kuksenok O. Nanogel Degradation at Soft Interfaces and in Bulk: Tracking Shape Changes and Interfacial Spreading. Macromolecules 2023. [DOI: 10.1021/acs.macromol.2c02470] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/18/2023]
Affiliation(s)
- Vaibhav Palkar
- Department of Materials Science and Engineering, Clemson University, Clemson, South Carolina 29634, United States
| | - Devanshu Thakar
- Department of Materials Science and Engineering, Clemson University, Clemson, South Carolina 29634, United States
- Department of Chemical Engineering, Indian Institute of Technology, Gandhinagar 382055, India
| | - Olga Kuksenok
- Department of Materials Science and Engineering, Clemson University, Clemson, South Carolina 29634, United States
| |
Collapse
|
6
|
Kopanichuk I, Scerbacova A, Ivanova A, Cheremisin A, Vishnyakov A. The effect of the molecular structure of alkyl ether carboxylate surfactants on the oil–water interfacial tension. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.119525] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
7
|
Koide Y, Goto S. Flow-induced scission of wormlike micelles in nonionic surfactant solutions under shear flow. J Chem Phys 2022; 157:084903. [DOI: 10.1063/5.0096830] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We investigate flow-induced scission of wormlike micelles with dissipative particle dynamics simulations of nonionic surfactant solutions under shear flow. To understand flow-induced scission in terms of micellar timescales, we propose a method to evaluate the longest relaxation time of unentangled surfactant micelles from the rotational relaxation time and the average lifetime at equilibrium. The mean squared displacement of surfactant molecules provides evidence that the longest relaxation time estimated by the proposed method serves as the characteristic timescale at equilibrium. We also demonstrate that the longest relaxation time plays an essential role in flow-induced scission. Using conditional statistics based on the aggregation number of micelles, we examine the statistical properties of the lifetime of wormlike micelles. We then conclude that flow-induced scission occurs when the Weissenberg number defined as the product of the longest relaxation time and the shear rate is larger than a threshold value.
Collapse
Affiliation(s)
- Yusuke Koide
- Osaka University Graduate School of Engineering Science Department of Mechanical Science and Bioengineering, Japan
| | - Susumu Goto
- Graduate School of Engineering Science, Osaka University, Japan
| |
Collapse
|
8
|
Chen QG, Lee MT. Anion Exchange Membranes for Fuel Cells Based on Quaternized Polystyrene-b-poly(ethylene-co-butylene)-b-polystyrene Triblock Copolymers with Spacer-Sidechain Design. Polymers (Basel) 2022; 14:polym14142860. [PMID: 35890636 PMCID: PMC9317406 DOI: 10.3390/polym14142860] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2022] [Revised: 07/08/2022] [Accepted: 07/12/2022] [Indexed: 02/04/2023] Open
Abstract
This work studied the polystyrene-b-poly(ethylene-co-butylene)-b-polystyrene (SEBS) triblock copolymers functionalized by butyl quaternary ammonium (C4Q) groups and alkyl side chains of different chain lengths (Cn, n = 0 to 24). The hydrated membrane morphology was modeled by dissipative particle dynamics simulation at hydration levels from 10 to 30. A hydroxide model was devised to characterize the diffusivity of anions under the coarse-grained framework. In general, the ionomers with alkyl side chains provided ion conductivity of a similar level at a lower ion exchange capacity. All hydrated SEBS–C4Q–Cn ionomers showed clear phase separation of the hydrophobic and hydrophilic domains, featuring 18.6 mS/cm to 36.8 mS/cm ion conductivity. The hydrophilic channels expanded as the water content increased, forming more effective ion conductive pathways. Introducing excess alkyl side chains enhanced the nano-segregation, leading to more ordered structures and longer correlation lengths of the aqueous phase. The membrane morphology was controlled by the length of alkyl side-chains as well as their tethering positions. Ionomers with functionalized side chains tethered on the same block resulted in well-connective water networks and higher conductivities. The detailed structural analysis provides synthesis guidelines to fabricate anion exchange membranes with improved performances.
Collapse
|
9
|
Faria BF, Vishnyakov AM. Simulation of surfactant adsorption at liquid-liquid interface: what we may expect from soft-core models?. J Chem Phys 2022; 157:094706. [DOI: 10.1063/5.0087363] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The present work attempt to explore systematically the surfactant sorption at liquid-liquid interfaces with coarse-grained models targeting thermodynamic properties of reference liquid solutions. We employ dissipative particle dynamics with soft-core forcefield tested against experimental data on micellization of surfactants in water, and the previous results are reproduced in this work. We consider three different nonionic surfactants: hexaethylene glycol monododecyl ether (C12E6), 2-[4-(2,4,4-trimethylpentan-2-yl)phenoxy]ethanol) knows as Triton X-100 (TХ-100), and two alkyl glucoside surfactants (CnG1) with n-alkane tail fragments and a saccharide hydrophilic head at decane-water and toluene-water interfaces. For TX-100, we composed a model based on the literature forcefield and found a good agreement with the experimental CMC. The head-head interactions are of different origins for different surfactant groups: entropic repulsion between ethylene oxide chains of C12E6 and TX-100, and more chemically specific and complex interactions between the maltose heads of alkyl glucosides. We interpret our results with the Redlich-Peterson equation of monolayer adsorption in order to relate the adsorption to the bulk concentration of the surfactant and the interfacial tension. The densities of the adsorbed monolayer at CMC mostly agree with the experimental data, and a reasonable agreement was obtained for the interfacial tension at CMC. At the same time, we found significant discrepancies between the simulated and experimental adsorption isotherms. We explain them by the oversimplified forcefield: when the parameters are fitted to the free energies of bulk solutions, they may not correctly reproduce the interfacial free energies.
Collapse
Affiliation(s)
| | - Aleksey M Vishnyakov
- Computational and Data-Intensive Science and Engineering, Skolkovo Institute of Science and Technology, Russia
| |
Collapse
|
10
|
Chen IC, Lee MT. Rhamnolipid Biosurfactants for Oil Recovery: Salt Effects on the Structural Properties Investigated by Mesoscale Simulations. ACS OMEGA 2022; 7:6223-6237. [PMID: 35224385 PMCID: PMC8867548 DOI: 10.1021/acsomega.1c06741] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/29/2021] [Accepted: 01/28/2022] [Indexed: 06/14/2023]
Abstract
Rhamnolipids (RLs) are biosurfactants produced by Pseudomonas. The biodegradability and the variety of their functionality make them suitable for environmental remediation and oil recovery. We use dissipative particle dynamics simulations to investigate the aggregation behaviors of ionic RL congeners with nonane in various operating conditions. Under zero-salinity conditions, all RL congeners studied here form small ellipsoidal clusters with detectable free surfactants. When salt ions are present, the electrostatic repulsion between the ionized heads is overcome, resulting in the formation of larger aggregates of unique structures. RLs with C10-alkyl tails tend to form elongated wormlike micelles, while RLs with C16-alkyl tails tend to form clusters in spherical symmetry, including vesicles. Di-rhamnolipids (dRLs) require stronger solvation than monorhamnolipids (mRLs) to form clusters, and the resulting size of micelles is decreased. The morphology of the mixed dRL/mRL/oil systems is controlled based on the type of the congeners and the oil contents. In addition, the divalent calcium ions are found to be influential to the structure of the micelles through different mechanisms. For 5 wt % salinity, the ionic RLs can form oil-swollen micelles up to a 1:1 surfactant-to-oil ratio, suggesting that ionic RLs are superb to act as cleaning agents for petroleum hydrocarbons in the marine area. These key findings may guide the design for RL-based washing techniques in enhanced oil recovery.
Collapse
|
11
|
Vishnyakov A, Mao R, Kam K, Potanin A, Neimark AV. Interactions of Crosslinked Polyacrylic Acid Polyelectrolyte Gels with Nonionic and Ionic Surfactants. J Phys Chem B 2021; 125:13817-13828. [PMID: 34905689 DOI: 10.1021/acs.jpcb.1c08638] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The morphology and stability of surfactant-loaded polyelectrolyte gels are of great interest for a variety of personal care, cosmetic, and pharmaceutical products. However, the mechanisms of surfactant interactions with gel-forming polymers are poorly understood and experimentally challenging. The aim of this work is to explore in silico the specifics of surfactant absorption within polyelectrolyte gels drawing on the examples of typical non-ionic octaethylene glycol monooctyl ether (C8E8) and anionic sodium dodecyl sulfate (SDS) surfactants and polyacrylic acid modified with hydrophobic sidechains mimicking the practically important Carbopol polymer. Using the systematically parameterized coarse-grained dissipative particle dynamics models, we generate and characterize the equilibrium conformations and swelling of the polymer films in aqueous solutions with the surfactant concentrations varied up to the critical micelle concentration (cmc). We discover the striking difference in interactions of Carbopol-like polymers with nonionic and ionic surfactants under mildly acidic conditions. The sorption of C8E8 within the polymer film is found substantial. As the surfactant concentration increases, the polymer film swells and, close to cmc, becomes unstable due to the formation and growth of water pockets filled with surfactant micelles. Sorption of SDS at the same bulk concentrations is found much lower, with only about 1% of surfactant mass fraction achieved at cmc. As the SDS concentration increases further, a lamellae structure is formed within the film, which remains stable. Reduced swelling and higher stability indicate better prospects of using SDS-type surfactants with Carbopol-based gels in formulations for detergents and personal care products.
Collapse
Affiliation(s)
- Aleksey Vishnyakov
- Department of Chemical and Biochemical Engineering, Rutgers University, New Brunswick, New Jersey 08854, United States.,Skolkovo Institute of Technology, Moscow 143005, Russia
| | - Runfang Mao
- Department of Chemical and Biochemical Engineering, Rutgers University, New Brunswick, New Jersey 08854, United States
| | - Kimberly Kam
- Department of Chemical and Biochemical Engineering, Rutgers University, New Brunswick, New Jersey 08854, United States
| | - Andrei Potanin
- Colgate-Palmolive, Piscataway, New Jersey 08855, United States
| | - Alexander V Neimark
- Department of Chemical and Biochemical Engineering, Rutgers University, New Brunswick, New Jersey 08854, United States
| |
Collapse
|
12
|
Santo KP, Neimark AV. Dissipative particle dynamics simulations in colloid and Interface science: a review. Adv Colloid Interface Sci 2021; 298:102545. [PMID: 34757286 DOI: 10.1016/j.cis.2021.102545] [Citation(s) in RCA: 32] [Impact Index Per Article: 10.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2021] [Revised: 10/07/2021] [Accepted: 10/08/2021] [Indexed: 12/31/2022]
Abstract
Dissipative particle dynamics (DPD) is one of the most efficient mesoscale coarse-grained methodologies for modeling soft matter systems. Here, we comprehensively review the progress in theoretical formulations, parametrization strategies, and applications of DPD over the last two decades. DPD bridges the gap between the microscopic atomistic and macroscopic continuum length and time scales. Numerous efforts have been performed to improve the computational efficiency and to develop advanced versions and modifications of the original DPD framework. The progress in the parametrization techniques that can reproduce the engineering properties of experimental systems attracted a lot of interest from the industrial community longing to use DPD to characterize, help design and optimize the practical products. While there are still areas for improvements, DPD has been efficiently applied to numerous colloidal and interfacial phenomena involving phase separations, self-assembly, and transport in polymeric, surfactant, nanoparticle, and biomolecules systems.
Collapse
Affiliation(s)
- Kolattukudy P Santo
- Department of Chemical and Biochemical Engineering, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, United States
| | - Alexander V Neimark
- Department of Chemical and Biochemical Engineering, Rutgers, The State University of New Jersey, Piscataway, NJ 08854, United States.
| |
Collapse
|
13
|
Santo KP, Neimark AV. Effects of metal-polymer complexation on structure and transport properties of metal-substituted polyelectrolyte membranes. J Colloid Interface Sci 2021; 602:654-668. [PMID: 34147755 DOI: 10.1016/j.jcis.2021.06.018] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2021] [Revised: 06/01/2021] [Accepted: 06/02/2021] [Indexed: 11/30/2022]
Abstract
Morphological and transport properties of hydrated metal-substituted Nafion membranes doped with metal ions of different valency and coordination strength are explored using coarse-grained dissipative particle dynamics simulations. To incorporate the effects of metal-polymer complexation, we introduce a novel metal ion complexation model, in which the charged central metal ion is surrounded by dummy sites that coordinate with ligands. The model parameters are determined by matching the metal-ligand running coordination numbers and the diffusion coefficients obtained from atomistic simulations and/or experiments. The increase of valency and coordination strength is found to strongly influence both the morphology and transport characteristics of the membrane at all hydration levels. The membrane segregation into hydrophobic and hydrophilic sub-phases is affected by metal-sulphonate coordination induced crosslinking at the hydrophilic/hydrophobic interface. The simulation results indicate that the interfacial crosslinking influences the interfacial tension and thereby affect the growth and coalescence of water clusters upon the increase of hydration. Multivalent complexation hinders water and ion mobility and causes anomalous sub-diffusion and dramatic decrease of the water permeability and ionic conductivity. Our DPD model is found efficient in elucidating the mechanisms of coordination-induced cross-linking and complexation and predicting on a semi-quantitative level the morphological and transport properties of metal-substituted Nafion membranes depending on the ion valency and coordination strength. The proposed model can be further advanced and adopted for other polyelectrolyte systems, such as sulfonated block-copolymers, polysaccharide solutions and composites, and biopolymer assemblies.
Collapse
Affiliation(s)
- Kolattukudy P Santo
- Department of Chemical and Biochemical Engineering, Rutgers, The State University of New Jersey, 98 Brett Rd, Piscataway, NJ 08854, USA
| | - Alexander V Neimark
- Department of Chemical and Biochemical Engineering, Rutgers, The State University of New Jersey, 98 Brett Rd, Piscataway, NJ 08854, USA.
| |
Collapse
|
14
|
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: 3] [Impact Index Per Article: 1.0] [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
|
15
|
Lee MT. Designing Highly Conductive Block Copolymer-Based Anion Exchange Membranes by Mesoscale Simulations. J Phys Chem B 2021; 125:2729-2740. [PMID: 33719456 DOI: 10.1021/acs.jpcb.0c10909] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
Hydroxide ion conductivity is a key aspect of anion exchange membranes and is mainly determined by the nanoscale membrane morphologies. Fundamental understanding of the structural and transport properties of membranes in terms of polymer architectures is crucial for future development of membrane-based applications. Using mesoscale simulations, this work predicts the mesostructure of the hydrated triblock copolymers; the designed polymers are composed of aromatic (polyphenylene oxide, PPO) or aliphatic (polystyrene-ethylene-butylene-styrene, SEBS) backbones, with cationic side chains being modified by hydrophobic or hydrophilic spacers. For PPO-based polymers, using octyl spacers creates a meshlike water network, yielding ion conductivity equal to 30.6 mS/cm at room temperature. For SEBS-based polymers, the nonmodified form is sufficient to produce ion-conducting pathways. Adding hydrophobic spacers further enhances the nanosegregation, and the membranes provide similar conductivity at a lower ion exchange capacity and water content. Adding hydrophilic spacers, however, has negative impacts on the ion transport. The side chains are in the stretched configurations, which sterically hinder the mobility of water and hydroxide ions. Such a resistance can be overcome by adapting multication side-chain designs, where large water channels are formed, yielding ion conductivity as high as 32.8 mS/cm.
Collapse
Affiliation(s)
- Ming-Tsung Lee
- Department of Chemical Engineering and Biotechnology, National Taipei University of Technology, Taipei 10608, Taiwan
| |
Collapse
|
16
|
Wang X, Santo KP, Neimark AV. Modeling Gas-Liquid Interfaces by Dissipative Particle Dynamics: Adsorption and Surface Tension of Cetyl Trimethyl Ammonium Bromide at the Air-Water Interface. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2020; 36:14686-14698. [PMID: 33216560 DOI: 10.1021/acs.langmuir.0c02572] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/27/2023]
Abstract
Adsorption of surfactants at gas-liquid interfaces that causes reduction in the surface tension is a classical problem in colloid and interface science with multiple practical applications in oil and gas recovery, separations, cosmetics, personal care, and biomedicine. Here, we develop an original coarse-grained model of the liquid-gas interface within the conventional dissipative particle dynamics (DPD) framework with the goal of quantitatively predicting the surface tension in the presence of surfactants. As a practical case-study example, we explore the adsorption of the cationic surfactant cetyl trimethyl ammonium bromide (CTAB) on the air-water interface. The gas phase is modeled as a DPD fluid composed of fictitious hard-core "gas" beads with exponentially decaying repulsive potentials to prevent penetration of the liquid phase components. A rigorous parametrization scheme is proposed based on matching the bulk and interfacial properties of water and octane taken as the reference compounds. Quantitative agreement between the simulated and experimental surface tension of CTAB solutions is found for a wide range of bulk surfactant concentrations (∼10-3 to ∼1 mmol/L) with the reduction of the surface tension from ∼72 mN/m (pure water) to the limiting value of ∼37.5 mN/m at the critical micelle concentration. The gas phase DPD model with the proposed parametrization scheme can be extended and applied to modeling various gas-liquid interfaces with surfactant and lipid monolayers, such as bubble suspensions, foams, froths, etc.
Collapse
Affiliation(s)
- Xinyang Wang
- Chemical and Biochemical Engineering Department, Rutgers, The State University of New Jersey, Piscataway, New Jersey 08854, United States
- School of Resources and Civil Engineering, Northeastern University, Shenyang, Liaoning 110819, China
| | - Kolattukudy P Santo
- Chemical and Biochemical Engineering Department, Rutgers, The State University of New Jersey, Piscataway, New Jersey 08854, United States
| | - Alexander V Neimark
- Chemical and Biochemical Engineering Department, Rutgers, The State University of New Jersey, Piscataway, New Jersey 08854, United States
| |
Collapse
|
17
|
Bray DJ, Anderson RL, Warren PB, Lewtas K. Wax Formation in Linear and Branched Alkanes with Dissipative Particle Dynamics. J Chem Theory Comput 2020; 16:7109-7122. [PMID: 32857939 DOI: 10.1021/acs.jctc.0c00605] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
We present a dissipative particle dynamics (DPD) model for wax formation (i.e., the freezing transition) in linear and branched alkanes at room temperature (298 K) and atmospheric pressure. We parametrize the model using pure liquid phase densities and the onset of wax formation as a function of alkyl chain length. Significant emphasis is placed on building an accurate representation of the underlying molecular architecture by careful consideration of bond lengths and angles, aided by distributions obtained from molecular dynamics simulation. Using the derived model, we observe wax formation in n-alkanes when the alkyl chain length is greater than 18 (n-octadecane), in excellent agreement with experimental observations. Further, we reproduce the behavior of branched alkanes and mixtures including solubilities of heavy alkanes in light alkane solvents.
Collapse
Affiliation(s)
- David J Bray
- The Hartree Centre, STFC Daresbury Laboratory, Warrington WA4 4AD, United Kingdom
| | - Richard L Anderson
- The Hartree Centre, STFC Daresbury Laboratory, Warrington WA4 4AD, United Kingdom
| | - Patrick B Warren
- The Hartree Centre, STFC Daresbury Laboratory, Warrington WA4 4AD, United Kingdom.,Unilever R&D Port Sunlight, Quarry Road East, Bebington, Wirral CH63 3JW, United Kingdom
| | - Kenneth Lewtas
- Lewtas Science & Technologies Ltd., 246 Banbury Road, Oxford OX2 7DY, United Kingdom.,School of Chemistry, The University of Edinburgh, Joseph Black Building, David Brewster Road, Edinburgh EH9 3FJ, United Kingdom
| |
Collapse
|
18
|
Pieczywek PM, Płaziński W, Zdunek A. Dissipative particle dynamics model of homogalacturonan based on molecular dynamics simulations. Sci Rep 2020; 10:14691. [PMID: 32895471 PMCID: PMC7477560 DOI: 10.1038/s41598-020-71820-2] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2020] [Accepted: 08/19/2020] [Indexed: 12/23/2022] Open
Abstract
In this study we present an alternative dissipative particle dynamics (DPD) parametrization strategy based on data extracted from the united-atom molecular simulations. The model of the homogalacturonan was designed to test the ability of the formation of large-scale structures via hydrogen bonding in water. The extraction of coarse-grained parameters from atomistic molecular dynamics was achieved by means of the proposed molecule aggregation algorithm based on an iterative nearest neighbour search. A novel approach to a time-scale calibration scheme based on matching the average velocities of coarse-grained particles enabled the DPD forcefield to reproduce essential structural features of homogalacturonan molecular chains. The successful application of the proposed parametrization method allowed for the reproduction of the shapes of radial distribution functions, particle velocities and diffusivity of the atomistic molecular dynamics model using DPD force field. The structure of polygalacturonic acid molecules was mapped into the DPD force field by means of the distance and angular bond characteristics, which closely matched the MD results. The resulting DPD trajectories showed that randomly dispersed homogalacturonan chains had a tendency to aggregate into highly organized 3D structures. The final structure resembled a three-dimensional network created by tightly associated homogalacturonan chains organized into thick fibres.
Collapse
Affiliation(s)
- P M Pieczywek
- Institute of Agrophysics, Polish Academy of Sciences, Doświadczalna 4, 20-270, Lublin, Poland.
| | - W Płaziński
- Jerzy Haber Institute of Catalysis and Surface Chemistry, Polish Academy of Sciences, Niezapominajek 8, 30-239, Cracow, Poland
| | - A Zdunek
- Institute of Agrophysics, Polish Academy of Sciences, Doświadczalna 4, 20-270, Lublin, Poland
| |
Collapse
|
19
|
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
|
20
|
Aryal D, Howard MP, Samanta R, Antoine S, Segalman R, Truskett TM, Ganesan V. Influence of pore morphology on the diffusion of water in triblock copolymer membranes. J Chem Phys 2020; 152:014904. [PMID: 31914764 DOI: 10.1063/1.5128119] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Understanding the transport properties of water in self-assembled block copolymer morphologies is important for furthering the use of such materials as water-purifying membranes. In this study, we used coarse-grained dissipative particle dynamics simulations to clarify the influence of pore morphology on the self-diffusion of water in linear-triblock-copolymer membranes. We considered representative lamellar, cylindrical, and gyroid morphologies and present results for both the global and local diffusivities of water in the pores. Our results suggest that the diffusivity of water in the confined, polymer-coated pores differs from that in the unconfined bulk. Explicitly, in confinement, the mobility of water is reduced by the hydrodynamic friction arising from the hydrophilic blocks coating the pore walls. We demonstrate that in lamella and cylindrical morphologies, the latter effects can be rendered as a universal function of the pore size relative to the brush height of the hydrophilic blocks.
Collapse
Affiliation(s)
- Dipak Aryal
- Department of Chemical Engineering, University of Texas at Austin, Austin, Texas 78712, USA
| | - Michael P Howard
- Department of Chemical Engineering, University of Texas at Austin, Austin, Texas 78712, USA
| | - Rituparna Samanta
- Department of Chemical Engineering, University of Texas at Austin, Austin, Texas 78712, USA
| | - Segolene Antoine
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, USA
| | - Rachel Segalman
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, USA
| | - Thomas M Truskett
- Department of Chemical Engineering, University of Texas at Austin, Austin, Texas 78712, USA
| | - Venkat Ganesan
- Department of Chemical Engineering, University of Texas at Austin, Austin, Texas 78712, USA
| |
Collapse
|
21
|
Minkara MS, Hembree RH, Jamadagni SN, Ghobadi AF, Eike DM, Siepmann JI. A new equation of state for homo-polymers in dissipative particle dynamics. J Chem Phys 2019; 150:124104. [PMID: 30927875 DOI: 10.1063/1.5058280] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A chain-revised Groot-Warren equation of state (crGW-EOS) was developed and tested to describe systems of homo-oligomeric chains in the framework of dissipative particle dynamics (DPD). First, thermodynamic perturbation theory is applied to introduce correction terms that account for the reduction in pressure with an increasing number of bonds at constant bead number density. Then, this EOS is modified by introducing a set of switching functions that yields an accurate second virial coefficient in the low-density limit. The crGW-EOS offers several improvements over the revised Groot-Warren equation of state and Groot-Warren equation of state for chain molecules. We tested the crGW-EOS by using it to predict the pressure of oligomeric systems and the B2 virial coefficient of chain DPD particles for a range of bond lengths. Additionally, a method is developed for determining the strength of cross-interaction parameters between chains of different compositions and sizes and for thermal and athermal mixtures. We explored how different levels of coarse-graining affect the upper-critical solution temperature.
Collapse
Affiliation(s)
- Mona S Minkara
- Department of Chemistry and Chemical Theory Center, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455, USA
| | - Robert H Hembree
- Department of Chemistry and Chemical Theory Center, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455, USA
| | - Sumanth N Jamadagni
- Computational Chemistry, Modeling and Simulation, The Procter & Gamble Company, 8611 Beckett Road, West Chester, Ohio 45069, USA
| | - Ahmad F Ghobadi
- Computational Chemistry, Modeling and Simulation, The Procter & Gamble Company, 8611 Beckett Road, West Chester, Ohio 45069, USA
| | - David M Eike
- Computational Chemistry, Modeling and Simulation, The Procter & Gamble Company, 8611 Beckett Road, West Chester, Ohio 45069, USA
| | - J Ilja Siepmann
- Department of Chemistry and Chemical Theory Center, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455, USA
| |
Collapse
|
22
|
Liu L, den Otter WK, Briels WJ. Coarse-Grained Simulations of Three-Armed Star Polymer Melts and Comparison with Linear Chains. J Phys Chem B 2018; 122:10210-10218. [DOI: 10.1021/acs.jpcb.8b03104] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Li Liu
- Department of Information Science and Engineering, Dalian Polytechnic University, Dalian 116034, China
| | | | - Wim J. Briels
- Forschungszentrum Jülich, ICS 3, D-52425 Jülich, Germany
| |
Collapse
|
23
|
Steinmetz D, Creton B, Lachet V, Rousseau B, Nieto-Draghi C. Simulations of Interfacial Tension of Liquid-Liquid Ternary Mixtures Using Optimized Parametrization for Coarse-Grained Models. J Chem Theory Comput 2018; 14:4438-4454. [PMID: 29906108 DOI: 10.1021/acs.jctc.8b00357] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
In this work, liquid-liquid systems are studied by means of coarse-grained Monte Carlo simulations (CG-MC) and Dissipative Particle Dynamics (DPD). A methodology is proposed to reproduce liquid-liquid equilibrium (LLE) and to provide variation of interfacial tension (IFT), as a function of the solute concentration. A key step is the parametrization method based on the use of the Flory-Huggins parameter between DPD beads to calculate solute/solvent interactions. Parameters are determined using a set of experimental compositional data of LLE, following four different approaches. These approaches are evaluated, and the results obtained are compared to analyze advantages/disadvantages of each one. These methodologies have been compared through their application on six systems: water/benzene/1,4-dioxane,water/chloroform/acetone, water/benzene/acetic acid, water/benzene/2-propanol, water/hexane/acetone, and water/hexane/2-propanol. CG-MC simulations in the Gibbs (NVT) ensemble have been used to check the validity of parametrization approaches for LLE reproduction. Then, CG-MC simulations in the osmotic (μsoluteNsolventP zzT) ensemble were carried out considering the two liquid phases with an explicit interface. This step allows one to work at the same bulk concentrations as the experimental data by imposing the precise bulk phase compositions and predicting the interface composition. Finally, DPD simulations were used to predict IFT values for different solute concentrations. Our results on variation of IFT with solute concentration in bulk phases are in good agreement with experimental data, but some deviations can be observed for systems containing hexane molecules.
Collapse
Affiliation(s)
- David Steinmetz
- IFP Energies nouvelles , 1 et 4 avenue de Bois-Préau , 92852 Rueil-Malmaison , France
| | - Benoit Creton
- IFP Energies nouvelles , 1 et 4 avenue de Bois-Préau , 92852 Rueil-Malmaison , France
| | - Véronique Lachet
- IFP Energies nouvelles , 1 et 4 avenue de Bois-Préau , 92852 Rueil-Malmaison , France.,Laboratoire de Chimie Physique , Université Paris-Sud , UMR 8000 CNRS , 91405 Orsay , France
| | - Bernard Rousseau
- Laboratoire de Chimie Physique , Université Paris-Sud , UMR 8000 CNRS , 91405 Orsay , France
| | - Carlos Nieto-Draghi
- IFP Energies nouvelles , 1 et 4 avenue de Bois-Préau , 92852 Rueil-Malmaison , France
| |
Collapse
|
24
|
Vishnyakov A, Mao R, Lee MT, Neimark AV. Coarse-grained model of nanoscale segregation, water diffusion, and proton transport in Nafion membranes. J Chem Phys 2018; 148:024108. [PMID: 29331134 DOI: 10.1063/1.4997401] [Citation(s) in RCA: 39] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a coarse-grained model of the acid form of Nafion membrane that explicitly includes proton transport. This model is based on a soft-core bead representation of the polymer implemented into the dissipative particle dynamics (DPD) simulation framework. The proton is introduced as a separate charged bead that forms dissociable Morse bonds with water beads. Morse bond formation and breakup artificially mimics the Grotthuss hopping mechanism of proton transport. The proposed DPD model is parameterized to account for the specifics of the conformations and flexibility of the Nafion backbone and sidechains; it treats electrostatic interactions in the smeared charge approximation. The simulation results qualitatively, and in many respects quantitatively, predict the specifics of nanoscale segregation in the hydrated Nafion membrane into hydrophobic and hydrophilic subphases, water diffusion, and proton mobility. As the hydration level increases, the hydrophilic subphase exhibits a percolation transition from a collection of isolated water clusters to a 3D network of pores filled with water embedded in the hydrophobic matrix. The segregated morphology is characterized in terms of the pore size distribution with the average size growing with hydration from ∼1 to ∼4 nm. Comparison of the predicted water diffusivity with the experimental data taken from different sources shows good agreement at high and moderate hydration and substantial deviation at low hydration, around and below the percolation threshold. This discrepancy is attributed to the dynamic percolation effects of formation and rupture of merging bridges between the water clusters, which become progressively important at low hydration, when the coarse-grained model is unable to mimic the fine structure of water network that includes singe molecule bridges. Selected simulations of water diffusion are performed for the alkali metal substituted membrane which demonstrate the effects of the counter-ions on membrane self-assembly and transport. The hydration dependence of the proton diffusivity reproduces semi-qualitatively the trend of the diverse experimental data, showing a sharp decrease around the percolation threshold. Overall, the proposed model opens up an opportunity to study self-assembly and water and proton transport in polyelectrolytes using computationally efficient DPD simulations, and, with further refinement, it may become a practical tool for theory informed design and optimization of perm-selective and ion-conducting membranes with improved properties.
Collapse
Affiliation(s)
- Aleksey Vishnyakov
- Department of Chemical and Biochemical Engineering, Rutgers, The State University of New Jersey, 98 Brett Road, Piscataway, New Jersey 08854-8058, USA
| | - Runfang Mao
- Department of Chemical and Biochemical Engineering, Rutgers, The State University of New Jersey, 98 Brett Road, Piscataway, New Jersey 08854-8058, USA
| | - Ming-Tsung Lee
- Department of Chemical and Biochemical Engineering, Rutgers, The State University of New Jersey, 98 Brett Road, Piscataway, New Jersey 08854-8058, USA
| | - Alexander V Neimark
- Department of Chemical and Biochemical Engineering, Rutgers, The State University of New Jersey, 98 Brett Road, Piscataway, New Jersey 08854-8058, USA
| |
Collapse
|
25
|
Saathoff J. Effectively parameterizing dissipative particle dynamics using COSMO-SAC: A partition coefficient study. J Chem Phys 2018; 148:154102. [DOI: 10.1063/1.5019952] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jonathan Saathoff
- ExxonMobil Research and Engineering Company, Annandale, New Jersey 08801-3059, USA
| |
Collapse
|
26
|
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
|
27
|
Anderson RL, Bray DJ, Ferrante AS, Noro MG, Stott IP, Warren PB. Dissipative particle dynamics: Systematic parametrization using water-octanol partition coefficients. J Chem Phys 2017; 147:094503. [DOI: 10.1063/1.4992111] [Citation(s) in RCA: 48] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
| | - David J. Bray
- STFC Hartree Centre, Scitech Daresbury, Warrington WA4 4AD, United Kingdom
| | - Andrea S. Ferrante
- Novidec Ltd., 3 Brook Hey, Parkgate, Neston CH64 6TH, United Kingdom
- Ferrante Scientific Ltd., 5 Croft Lane, Bromborough CH62 2BX, United Kingdom
| | - Massimo G. Noro
- Unilever R&D Port Sunlight, Quarry Road East, Bebington CH63 3JW, United Kingdom
| | - Ian P. Stott
- Unilever R&D Port Sunlight, Quarry Road East, Bebington CH63 3JW, United Kingdom
| | - Patrick B. Warren
- Unilever R&D Port Sunlight, Quarry Road East, Bebington CH63 3JW, United Kingdom
| |
Collapse
|
28
|
Gooneie A, Schuschnigg S, Holzer C. A Review of Multiscale Computational Methods in Polymeric Materials. Polymers (Basel) 2017; 9:E16. [PMID: 30970697 PMCID: PMC6432151 DOI: 10.3390/polym9010016] [Citation(s) in RCA: 113] [Impact Index Per Article: 16.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2016] [Revised: 12/07/2016] [Accepted: 12/22/2016] [Indexed: 11/17/2022] Open
Abstract
Polymeric materials display distinguished characteristics which stem from the interplay of phenomena at various length and time scales. Further development of polymer systems critically relies on a comprehensive understanding of the fundamentals of their hierarchical structure and behaviors. As such, the inherent multiscale nature of polymer systems is only reflected by a multiscale analysis which accounts for all important mechanisms. Since multiscale modelling is a rapidly growing multidisciplinary field, the emerging possibilities and challenges can be of a truly diverse nature. The present review attempts to provide a rather comprehensive overview of the recent developments in the field of multiscale modelling and simulation of polymeric materials. In order to understand the characteristics of the building blocks of multiscale methods, first a brief review of some significant computational methods at individual length and time scales is provided. These methods cover quantum mechanical scale, atomistic domain (Monte Carlo and molecular dynamics), mesoscopic scale (Brownian dynamics, dissipative particle dynamics, and lattice Boltzmann method), and finally macroscopic realm (finite element and volume methods). Afterwards, different prescriptions to envelope these methods in a multiscale strategy are discussed in details. Sequential, concurrent, and adaptive resolution schemes are presented along with the latest updates and ongoing challenges in research. In sequential methods, various systematic coarse-graining and backmapping approaches are addressed. For the concurrent strategy, we aimed to introduce the fundamentals and significant methods including the handshaking concept, energy-based, and force-based coupling approaches. Although such methods are very popular in metals and carbon nanomaterials, their use in polymeric materials is still limited. We have illustrated their applications in polymer science by several examples hoping for raising attention towards the existing possibilities. The relatively new adaptive resolution schemes are then covered including their advantages and shortcomings. Finally, some novel ideas in order to extend the reaches of atomistic techniques are reviewed. We conclude the review by outlining the existing challenges and possibilities for future research.
Collapse
Affiliation(s)
- Ali Gooneie
- Chair of Polymer Processing, Montanuniversitaet Leoben, Otto Gloeckel-Strasse 2, 8700 Leoben, Austria.
| | - Stephan Schuschnigg
- Chair of Polymer Processing, Montanuniversitaet Leoben, Otto Gloeckel-Strasse 2, 8700 Leoben, Austria.
| | - Clemens Holzer
- Chair of Polymer Processing, Montanuniversitaet Leoben, Otto Gloeckel-Strasse 2, 8700 Leoben, Austria.
| |
Collapse
|