1
|
Wang Y, Hernandez R. Construction of Multiscale Dissipative Particle Dynamics (DPD) Models from Other Coarse-Grained Models. ACS OMEGA 2024; 9:17667-17680. [PMID: 38645334 PMCID: PMC11025104 DOI: 10.1021/acsomega.4c01868] [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: 02/26/2024] [Revised: 03/06/2024] [Accepted: 03/15/2024] [Indexed: 04/23/2024]
Abstract
We present a general scheme for converting coarse-grained models into Dissipative Particle Dynamics (DPD) models. We build the corresponding DPD models by analogy with the de novo DPD coarse-graining scheme suggested by Groot and Warren (J. Chem. Phys., 1997). Electrostatic interactions between charged DPD particles are represented though the addition of a long-range Slater Coulomb potential as suggested by González-Melchor et al. (J. Chem. Phys., 2006). The construction is illustrated by converting MARTINI models for various proteins into a DPD representation, but it not restricted to the usual potential form in the MARTINI model-viz., Lennard-Jones potentials. We further extended the DPD scheme away from the typical use of homogeneous particle sizes, therefore faithfully representing the variations in the particle sizes seen in the underlying MARTINI model. The accuracy of the resulting construction of our generalized DPD models with respect to several structural observables has been benchmarked favorably against all-atom and MARTINI models for a selected set of peptides and proteins, and variations in the scales of the coarse-graining of the water solvent.
Collapse
Affiliation(s)
- Yinhan Wang
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Rigoberto Hernandez
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, United States
| |
Collapse
|
2
|
Leister N, Götz V, Jan Bachmann S, Nachtigall S, Hosseinpour S, Peukert W, Karbstein H. A comprehensive methodology to study double emulsion stability. J Colloid Interface Sci 2023; 630:534-548. [DOI: 10.1016/j.jcis.2022.10.119] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2022] [Revised: 10/15/2022] [Accepted: 10/22/2022] [Indexed: 11/05/2022]
|
3
|
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
|
4
|
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
|
5
|
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: 35] [Impact Index Per Article: 11.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
|
6
|
Fayaz‐Torshizi M, Müller EA. Coarse‐Grained Molecular Simulation of Polymers Supported by the Use of the SAFT‐γ$\gamma$ Mie Equation of State. MACROMOL THEOR SIMUL 2021. [DOI: 10.1002/mats.202100031] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
| | - Erich A. Müller
- Department of Chemical Engineering Imperial College London London SW7 2AZ UK
| |
Collapse
|
7
|
van der Haven DLH, Köhler S, Schreiner E, In 't Veld PJ. Closed-Form Coexistence Equation for Phase Separation of Polymeric Mixtures in Dissipative Particle Dynamics. J Phys Chem B 2021; 125:7485-7498. [PMID: 34196184 DOI: 10.1021/acs.jpcb.0c11274] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
To date, no extensive study of the phase diagram for binary fluid mixtures in dissipative particle dynamics (DPD) has been published. This is especially pertinent for newer parameterization schemes where the self-self interaction, or the effective volume, of different particle types is varied. This work presents an exhaustive study of the parameter space concerning DPD particles with soft interaction potentials. Moreover, we propose a closed-form coexistence equation or binodal curve that is inspired by the Flory-Huggins model. This equation describes the phase diagram of all binary mixtures made up out of monomers, homopolymers, and the mixtures thereof when self-self interactions are varied. The mean absolute percentage error (MAPE) of the equation on simulated data, including validation simulations, is 1.02%. The equation can a priori predict the phase separation of mixtures using only DPD interaction parameters. The proposed coexistence equation can therefore be used to directly validate interaction parameters resulting from novel parameterization schemes, including coarse graining and equations of state, without the need for additional simulations. Finally, it is shown that the choice of bond potential markedly influences phase behavior.
Collapse
Affiliation(s)
- Dingeman L H van der Haven
- Department of Chemical Engineering and Chemistry, Eindhoven University of Technology, Eindhoven 5600 MB, The Netherlands
| | - Stephan Köhler
- Polymer Physics, BASF SE, Ludwigshafen am Rhein 67056, Germany
| | | | | |
Collapse
|
8
|
Sherck N, Shen K, Nguyen M, Yoo B, Köhler S, Speros JC, Delaney KT, Shell MS, Fredrickson GH. Molecularly Informed Field Theories from Bottom-up Coarse-Graining. ACS Macro Lett 2021; 10:576-583. [PMID: 35570772 DOI: 10.1021/acsmacrolett.1c00013] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Polymer formulations possessing mesostructures or phase coexistence are challenging to simulate using atomistic particle-explicit approaches due to the disparate time and length scales, while the predictive capability of field-based simulations is hampered by the need to specify interactions at a coarser scale (e.g., χ-parameters). To overcome the weaknesses of both, we introduce a bottom-up coarse-graining methodology that leverages all-atom molecular dynamics to molecularly inform coarser field-theoretic models. Specifically, we use relative-entropy coarse-graining to parametrize particle models that are directly and analytically transformable into statistical field theories. We demonstrate the predictive capability of this approach by reproducing experimental aqueous poly(ethylene oxide) (PEO) cloud-point curves with no parameters fit to experimental data. This synergistic approach to multiscale polymer simulations opens the door to de novo exploration of phase behavior across a wide variety of polymer solutions and melt formulations.
Collapse
Affiliation(s)
- Nicholas Sherck
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, United States
| | - Kevin Shen
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, United States
- Materials Research Laboratory, University of California, Santa Barbara, California 93106, United States
| | - My Nguyen
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, United States
| | - Brian Yoo
- BASF Corporation, Tarrytown, New York 10591, United States
| | | | - Joshua C. Speros
- California Research Alliance (CARA) by BASF, Berkeley, California 94720, United States
| | - Kris T. Delaney
- Materials Research Laboratory, University of California, Santa Barbara, California 93106, United States
| | - M. Scott Shell
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, United States
| | - Glenn H. Fredrickson
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, United States
- Materials Research Laboratory, University of California, Santa Barbara, California 93106, United States
- Department of Materials, University of California, Santa Barbara, California 93106, United States
| |
Collapse
|
9
|
Lavagnini E, Cook JL, Warren PB, Hunter CA. Translation of Chemical Structure into Dissipative Particle Dynamics Parameters for Simulation of Surfactant Self-Assembly. J Phys Chem B 2021; 125:3942-3952. [PMID: 33848165 PMCID: PMC8154614 DOI: 10.1021/acs.jpcb.1c00480] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2021] [Revised: 03/08/2021] [Indexed: 01/28/2023]
Abstract
Dissipative particle dynamics (DPD) can be used to simulate the self-assembly properties of surfactants in aqueous solutions, but in order to simulate a new compound, a large number of new parameters are required. New methods for the calculation of reliable DPD parameters directly from chemical structure are described, allowing the DPD approach to be applied to a much wider range of organic compounds. The parameters required to describe the bonded interactions between DPD beads were calculated from molecular mechanics structures. The parameters required to describe the nonbonded interactions were calculated from surface site interaction point (SSIP) descriptions of molecular fragments that represent individual beads. The SSIPs were obtained from molecular electrostatic potential surfaces calculated using density functional theory and used in the SSIMPLE algorithm to calculate transfer free energies between different bead liquids. This approach was used to calculate DPD parameters for a range of different types of surfactants, which include ester, amide, and sugar moieties. The parameters were used to simulate the self-assembly properties in aqueous solutions, and comparison of the results for 27 surfactants with the available experimental data shows that these DPD simulations accurately predict critical micelle concentrations, aggregation numbers, and the shapes of the supramolecular assemblies formed. The methods described here provide a general approach to determining DPD parameters for neutral organic compounds of arbitrary structure.
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.
| | - Christopher A. Hunter
- Department
of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, U.
K.
| |
Collapse
|
10
|
Deaton TA, Aydin F, Li NK, Chu X, Dutt M, Yingling YG. Dissipative Particle Dynamics Approaches to Modeling the Self-Assembly and Morphology of Neutral and Ionic Block Copolymers in Solution. FOUNDATIONS OF MOLECULAR MODELING AND SIMULATION 2021. [DOI: 10.1007/978-981-33-6639-8_4] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/09/2022]
|
11
|
Tzorín A, Zamarripa AL, Goicochea AG, Vallejo-Montesinos J. Effect of increasing the number of amino groups in the solubility of Copolysiloxanes using dissipative particle dynamics. MOLECULAR SIMULATION 2020. [DOI: 10.1080/08927022.2020.1800006] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Affiliation(s)
- Amanda Tzorín
- Edificio T-12, Facultad de Ciencias Químicas y Farmacia, Ciudad Universitaria, Guatemala, Guatemala
| | - Ana L. Zamarripa
- División de Ciencias Naturales y Exactas, Universidad de Guanajuato, Guanajuato, Mexico
| | - Armando Gama Goicochea
- División de Ingeniería Química y Bioquímica, Tecnológico de Estudios Superiores de Ecatepec, Ecatepec, Estado de México, Mexico
| | | |
Collapse
|
12
|
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
|
13
|
Xu Z, Gao L, Chen P, Yan LT. Diffusive transport of nanoscale objects through cell membranes: a computational perspective. SOFT MATTER 2020; 16:3869-3881. [PMID: 32236197 DOI: 10.1039/c9sm02338k] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Diffusion is an essential and fundamental means of transport of substances on cell membranes, and the dynamics of biomembranes plays a crucial role in the regulation of numerous cellular processes. The understanding of the complex mechanisms and the nature of particle diffusion have a bearing on establishing guidelines for the design of efficient transport materials and unique therapeutic approaches. Herein, this review article highlights the most recent advances in investigating diffusion dynamics of nanoscale objects on biological membranes, focusing on the approaches of tailored computer simulations and theoretical analysis. Due to the presence of the complicated and heterogeneous environment on native cell membranes, the diffusive transport behaviors of nanoparticles exhibit unique and variable characteristics. The general aspects and basic theories of normal diffusion and anomalous diffusion have been introduced. In addition, the influence of a series of external and internal factors on the diffusion behaviors is discussed, including particle size, membrane curvature, particle-membrane interactions or particle-inclusion, and the crowding degree of membranes. Finally, we seek to identify open problems in the existing experimental, simulation, and theoretical research studies, and to propose challenges for future development.
Collapse
Affiliation(s)
- Ziyang Xu
- Key Laboratory of Advanced Materials (MOE), Department of Chemical Engineering, Tsinghua University, Beijing 100084, P. R. China.
| | - Lijuan Gao
- Key Laboratory of Advanced Materials (MOE), Department of Chemical Engineering, Tsinghua University, Beijing 100084, P. R. China.
| | - Pengyu Chen
- Key Laboratory of Advanced Materials (MOE), Department of Chemical Engineering, Tsinghua University, Beijing 100084, P. R. China.
| | - Li-Tang Yan
- Key Laboratory of Advanced Materials (MOE), Department of Chemical Engineering, Tsinghua University, Beijing 100084, P. R. China.
| |
Collapse
|
14
|
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
|
15
|
Vanya P, Sharman J, Elliott JA. Invariance of experimental observables with respect to coarse-graining in standard and many-body dissipative particle dynamics. J Chem Phys 2019; 150:064101. [PMID: 30770006 DOI: 10.1063/1.5046851] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Dissipative particle dynamics (DPD) is a well-established mesoscale simulation method. However, there have been long-standing ambiguities regarding the dependence of its (purely repulsive) force field parameter on temperature as well as the variation of the resulting experimental observables, such as diffusivity or surface tension, with coarse-graining (CG) degree. Here, we rederive the temperature dependence of DPD interaction parameter and revisit the role of the CG degree in standard DPD simulations. Consequently, we derive a scaling of the input variables that renders the system properties invariant with respect to CG degree and illustrate the versatility of the method by computing the surface tensions of binary solvent mixtures. We then extend this procedure to many-body dissipative particle dynamics and, by computing surface tensions of the same mixtures at a range of CG degrees, demonstrate that this newer method, which has not been widely applied so far, is also capable of simulating complex fluids of practical interest.
Collapse
Affiliation(s)
- Peter Vanya
- Department of Materials Science and Metallurgy, University of Cambridge, 27 Charles Babbage Road, Cambridge CB3 0FS, United Kingdom
| | - Jonathan Sharman
- Johnson Matthey Technology Centre, Blounts Court Road, Sonning Common, Reading RG4 9NH, United Kingdom
| | - James A Elliott
- Department of Materials Science and Metallurgy, University of Cambridge, 27 Charles Babbage Road, Cambridge CB3 0FS, United Kingdom
| |
Collapse
|
16
|
Xu Z, Yang Y, Zhu G, Chen P, Huang Z, Dai X, Hou C, Yan L. Simulating Transport of Soft Matter in Micro/Nano Channel Flows with Dissipative Particle Dynamics. ADVANCED THEORY AND SIMULATIONS 2018. [DOI: 10.1002/adts.201800160] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Ziyang Xu
- State Key Laboratory of Chemical EngineeringDepartment of Chemical EngineeringTsinghua University Beijing 100084 China
| | - Ye Yang
- State Key Laboratory of Chemical EngineeringDepartment of Chemical EngineeringTsinghua University Beijing 100084 China
| | - Guolong Zhu
- State Key Laboratory of Chemical EngineeringDepartment of Chemical EngineeringTsinghua University Beijing 100084 China
| | - Pengyu Chen
- State Key Laboratory of Chemical EngineeringDepartment of Chemical EngineeringTsinghua University Beijing 100084 China
| | - Zihan Huang
- State Key Laboratory of Chemical EngineeringDepartment of Chemical EngineeringTsinghua University Beijing 100084 China
| | - Xiaobin Dai
- State Key Laboratory of Chemical EngineeringDepartment of Chemical EngineeringTsinghua University Beijing 100084 China
| | - Cuiling Hou
- State Key Laboratory of Chemical EngineeringDepartment of Chemical EngineeringTsinghua University Beijing 100084 China
| | - Li‐Tang Yan
- State Key Laboratory of Chemical EngineeringDepartment of Chemical EngineeringTsinghua University Beijing 100084 China
| |
Collapse
|
17
|
Vallejo-Montesinos J, Villegas A, Cervantes J, Pérez E, Goicochea AG. Study of Polymer-Solvent Interactions of Complex Polysiloxanes Using Dissipative Particle Dynamics. J MACROMOL SCI B 2018. [DOI: 10.1080/00222348.2018.1503336] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Affiliation(s)
| | - Antonio Villegas
- División de Ciencias Naturales y Exactas, Universidad de Guanajuato, Guanajuato, Mexico
| | - Jorge Cervantes
- División de Ciencias Naturales y Exactas, Universidad de Guanajuato, Guanajuato, Mexico
| | - Elías Pérez
- Instituto de Física, Universidad Autónoma de San Luis Potosí, San Luis Potosí, Mexico
| | - Armando Gama Goicochea
- División de Ingeniería Química y Bioquímica, Tecnológico de Estudios Superiores de Ecatepec, Ecatepec, Mexico
| |
Collapse
|
18
|
Soroush Barhaghi M, Torabi K, Nejahi Y, Schwiebert L, Potoff JJ. Molecular exchange Monte Carlo: A generalized method for identity exchanges in grand canonical Monte Carlo simulations. J Chem Phys 2018; 149:072318. [DOI: 10.1063/1.5025184] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023] Open
Affiliation(s)
- Mohammad Soroush Barhaghi
- Department of Chemical Engineering and Materials Science, Wayne State University, Detroit, Michigan 48202, USA
| | - Korosh Torabi
- Department of Chemical Engineering and Materials Science, Wayne State University, Detroit, Michigan 48202, USA
| | - Younes Nejahi
- Department of Computer Science, Wayne State University, Detroit, Michigan 48202, USA
| | - Loren Schwiebert
- Department of Computer Science, Wayne State University, Detroit, Michigan 48202, USA
| | - Jeffrey J. Potoff
- Department of Chemical Engineering and Materials Science, Wayne State University, Detroit, Michigan 48202, USA
| |
Collapse
|
19
|
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
|
20
|
Hideo D, Koji O, Takamitsu N, Sona S, Yuji M. A portable code for dissipative particle dynamics (DPD) simulations with additional specific interactions. CHEM-BIO INFORMATICS JOURNAL 2018. [DOI: 10.1273/cbij.18.70] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Affiliation(s)
- Doi Hideo
- Department of Chemistry and Research Center for Smart Molecules, Faculty of Science, Rikkyo University
| | - Okuwaki Koji
- Department of Chemistry and Research Center for Smart Molecules, Faculty of Science, Rikkyo University
| | - Naito Takamitsu
- Department of Chemistry and Research Center for Smart Molecules, Faculty of Science, Rikkyo University
| | - Saitou Sona
- Department of Chemistry and Research Center for Smart Molecules, Faculty of Science, Rikkyo University
| | - Mochizuki Yuji
- Department of Chemistry and Research Center for Smart Molecules, Faculty of Science, Rikkyo University
- Institute of Industrial Science, The University of Tokyo
| |
Collapse
|
21
|
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: 37] [Impact Index Per Article: 6.2] [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
|
22
|
Investigation of oil-in-water emulsion stability with relevant interfacial characteristics simulated by dissipative particle dynamics. Colloids Surf A Physicochem Eng Asp 2018. [DOI: 10.1016/j.colsurfa.2018.02.063] [Citation(s) in RCA: 38] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
23
|
Ma D, Keten S. Stable micelles based on a mixture of coiled-coils: the role of different oligomeric states. NANOSCALE 2018; 10:7589-7596. [PMID: 29637959 DOI: 10.1039/c7nr09695j] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Homomeric micelles with tunable size, shape and stability have been extensively studied for biomedical applications such as drug carriers. However, designing the local valency and self-assembled morphology of nanophase-separated multicomponent micelles with varied ligand binding possibilities remains challenging. Here, we present micelles self-assembled from amphiphilic peptide-PEG-lipid hybrid conjugates, where the peptides can be either a 3-helix or 4-helix coiled-coil. We demonstrate that the micelle size and sphericity can be controlled based on the coiled-coil oligomeric state. Using theory and coarse-grained dissipative particle dynamics (DPD) simulations in an explicit solvent simulation, we studied the distribution of 3-helix and 4-helix conjugates within the mixed micelles and observed self-organization into nanodomains within the mixed micelle. We discovered that the phase separation behavior is dictated by the geometry mismatch in the alkyl chain length from different coiled-coil oligomeric states. Our analyses of the self-assembly tendency and drug delivery potency of mixed micelles with controlled multivalency provide important insights into the assembly and formation of nanophase-separated micelles.
Collapse
Affiliation(s)
- Dan Ma
- Department of Mechanical Engineering, Northwestern University, Evanston, Illinois 60208, USA.
| | | |
Collapse
|
24
|
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
|
25
|
Hou M, Liu T, Zhang B, Chen B, Li Q, Liu X, Lu C, Wang Z, Dang L. Alpha-linolenic acid-loaded oil/water microemulsion: Effects of phase behaviour simulation and environmental stress on phase stabilizing and anti-oxidation capacity. Food Chem 2018; 256:311-318. [PMID: 29606454 DOI: 10.1016/j.foodchem.2018.02.100] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2017] [Revised: 01/13/2018] [Accepted: 02/19/2018] [Indexed: 12/29/2022]
Abstract
α-Linolenic acid (ALA)-loaded microemulsion (ME) was prepared from isoamyl acetate, polyoxyethylene ether 35 (EL-35), ethanol and water. The dynamic phase behaviour was simulated using dissipative particle dynamics (DPD), which showed that spherical ME was formed at water/oil ratios of 1:9 and 9:1, while a lamellar structure with distinctive water-course and oil layer appeared at ratios of 3:7, 5:5, and 7:3. Phase stabilizing and anti-oxidation effect of environmental stresses on ALA-loaded microemulsion were investigated. Results showed that the ME region was large and had good environmental tolerance. Subsequently, the investigation of anti-oxidation stability revealed that more than 60% ALA of ALA-loaded ME could be protected from oxidation under environmental stresses. Furthermore, ALA-loaded ME was applied in aqueous-based foods. The transparency, precipitate, stratification and phase separation were used to evaluate influence of ME on product properties, confirming great feasibility and stability of ALA-loaded ME for practical applications.
Collapse
Affiliation(s)
- Mengna Hou
- School of Chemical Engineering and Technology, Tianjin University, Tianjin 300072, People's Republic of China
| | - Tiankuo Liu
- School of Chemical Engineering and Technology, Tianjin University, Tianjin 300072, People's Republic of China
| | - Bo Zhang
- School of Chemical Engineering and Technology, Tianjin University, Tianjin 300072, People's Republic of China
| | - Boru Chen
- School of Chemical Engineering and Technology, Tianjin University, Tianjin 300072, People's Republic of China
| | - Qing Li
- School of Chemical Engineering and Technology, Tianjin University, Tianjin 300072, People's Republic of China
| | - Xiaoxue Liu
- School of Chemical Engineering and Technology, Tianjin University, Tianjin 300072, People's Republic of China
| | - Chao Lu
- School of Chemical Engineering and Technology, Tianjin University, Tianjin 300072, People's Republic of China
| | - Zhanzhong Wang
- School of Chemical Engineering and Technology, Tianjin University, Tianjin 300072, People's Republic of China.
| | - Leping Dang
- School of Chemical Engineering and Technology, Tianjin University, Tianjin 300072, People's Republic of China.
| |
Collapse
|
26
|
Okuwaki K, Mochizuki Y, Doi H, Ozawa T. Fragment Molecular Orbital Based Parametrization Procedure for Mesoscopic Structure Prediction of Polymeric Materials. J Phys Chem B 2018; 122:338-347. [PMID: 29285920 DOI: 10.1021/acs.jpcb.7b08461] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In the analyses of miscibility behaviors of macromolecules and polymers, dissipative particle dynamics (DPD) simulations are generally performed. In these simulations, the so-called χ parameters describing the effective interactions among particles are crucial. It has been known that such parameters can be obtained within the classical or empirical force field frameworks. However, there is a potential problem that charge transfer and polarization occasionally occur. Additionally, satisfactory reference parameters are not available for some cases. Therefore, we developed a new procedure to evaluate the set of parameters by using the ab initio fragment molecular orbital (FMO) method which can provide the set of interaction energies among segments as polymer units. Moreover, we evaluated the anisotropy of molecules by using the FMO-based effective interaction parameters for three standard binary mixture systems (hexane-nitrobenzene, polyisobutylene-diisobutyl ketone, and polyisoprene-polystyrene). The calculated values showed good agreement with the experimental values with about 10% errors.
Collapse
Affiliation(s)
- Koji Okuwaki
- Department of Chemistry and Research Center for Smart Molecules, Faculty of Science, Rikkyo University , 3-34-1 Nishi-ikebukuro, Toshima-ku, Tokyo 171-8501, Japan
| | - Yuji Mochizuki
- Department of Chemistry and Research Center for Smart Molecules, Faculty of Science, Rikkyo University , 3-34-1 Nishi-ikebukuro, Toshima-ku, Tokyo 171-8501, Japan.,Institute of Industrial Science, The University of Tokyo , 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| | - Hideo Doi
- Department of Chemistry and Research Center for Smart Molecules, Faculty of Science, Rikkyo University , 3-34-1 Nishi-ikebukuro, Toshima-ku, Tokyo 171-8501, Japan
| | - Taku Ozawa
- JSOL Corporation , 2-5-24 Harumi, Chuo-ku, Tokyo 104-0053, Japan
| |
Collapse
|
27
|
OKUWAKI K, DOI H, MOCHIZUKI Y. An automated framework to evaluate effective interactionparameters for dissipative particle dynamics simulations basedon the fragment molecular orbital (FMO) method. JOURNAL OF COMPUTER CHEMISTRY-JAPAN 2018. [DOI: 10.2477/jccj.2017-0048] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Affiliation(s)
- Koji OKUWAKI
- Department of Chemistry and Research Center for Smart Molecules, College of Science, Rikkyo University, 3-34-1 Nishi-ikebukuro, Toshima-ku, Tokyo 171-8501, Japan
| | - Hideo DOI
- Department of Chemistry and Research Center for Smart Molecules, College of Science, Rikkyo University, 3-34-1 Nishi-ikebukuro, Toshima-ku, Tokyo 171-8501, Japan
| | - Yuji MOCHIZUKI
- Department of Chemistry and Research Center for Smart Molecules, College of Science, Rikkyo University, 3-34-1 Nishi-ikebukuro, Toshima-ku, Tokyo 171-8501, Japan
- Institute of Industrial Science, The University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8503, Japan
| |
Collapse
|
28
|
Okuwaki K, Mochizuki Y, Doi H, Kawada S, Ozawa T, Yasuoka K. Theoretical analyses on water cluster structures in polymer electrolyte membrane by using dissipative particle dynamics simulations with fragment molecular orbital based effective parameters. RSC Adv 2018; 8:34582-34595. [PMID: 35548624 PMCID: PMC9086946 DOI: 10.1039/c8ra07428c] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2018] [Accepted: 10/01/2018] [Indexed: 12/04/2022] Open
Abstract
The mesoscopic structures of polymer electrolyte membrane (PEM) affect the performances of fuel cells. Nafion® with the Teflon® backbone has been the most widely used of all PEMs, but sulfonated poly-ether ether-ketone (SPEEK) having an aromatic backbone has drawn interest as an alternative to Nafion. In the present study, a series of dissipative particle dynamics (DPD) simulations were performed to compare Nafion and SPEEK. These PEM polymers were modeled by connected particles corresponding to the hydrophobic backbone and the hydrophilic moiety of sulfonic acid group. The water particle interacting with Nafion particles was prepared as well. The crucial interaction parameters among DPD particles were evaluated by a series of calculations based on the fragment molecular orbital (FMO) method in a non-empirical way (Okuwaki et al., J. Phys. Chem. B, 2018, 122, 338–347). Through the DPD simulations, the water and hydrophilic particles aggregated, forming cluster networks surrounded by the hydrophobic phase. The structural features of formed water clusters were investigated in detail. Furthermore, the differences in percolation behaviors between Nafion and SPEEK revealed much better connectivity among water clusters by Nafion. The present FMO-DPD simulation results were in good agreement with available experimental data. The mesoscopic structures of polymer electrolyte membrane (PEM) affect the performances of fuel cells.![]()
Collapse
Affiliation(s)
- Koji Okuwaki
- Department of Chemistry and Research Center for Smart Molecules
- Faculty of Science
- Rikkyo University
- Toshima-ku
- Japan
| | - Yuji Mochizuki
- Department of Chemistry and Research Center for Smart Molecules
- Faculty of Science
- Rikkyo University
- Toshima-ku
- Japan
| | - Hideo Doi
- Department of Chemistry and Research Center for Smart Molecules
- Faculty of Science
- Rikkyo University
- Toshima-ku
- Japan
| | - Shutaro Kawada
- Department of Chemistry and Research Center for Smart Molecules
- Faculty of Science
- Rikkyo University
- Toshima-ku
- Japan
| | | | - Kenji Yasuoka
- Department of Mechanical Engineering
- Keio University
- Yokohama 223-8522
- Japan
| |
Collapse
|
29
|
Zhang Q, Lin J, Wang L, Xu Z. Theoretical modeling and simulations of self-assembly of copolymers in solution. Prog Polym Sci 2017. [DOI: 10.1016/j.progpolymsci.2017.04.003] [Citation(s) in RCA: 77] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
|
30
|
Guo Z, Kindt JT. Gibbs ensemble Monte Carlo with solvent repacking: phase coexistence of size–asymmetrical binary Lennard-Jones mixtures. MOLECULAR SIMULATION 2017. [DOI: 10.1080/08927022.2017.1373192] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Affiliation(s)
- Ziwei Guo
- Department of Chemistry, Emory University, Atlanta, GA, USA
| | - James T. Kindt
- Department of Chemistry, Emory University, Atlanta, GA, USA
| |
Collapse
|
31
|
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
|
32
|
Doi H, Okuwaki K, Mochizuki Y, Ozawa T, Yasuoka K. Dissipative particle dynamics (DPD) simulations with fragment molecular orbital (FMO) based effective parameters for 1-Palmitoyl-2-oleoyl phosphatidyl choline (POPC) membrane. Chem Phys Lett 2017. [DOI: 10.1016/j.cplett.2017.07.032] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
|
33
|
Fraaije JGEM, van Male J, Becherer P, Serral Gracià R. Coarse-Grained Models for Automated Fragmentation and Parametrization of Molecular Databases. J Chem Inf Model 2016; 56:2361-2377. [DOI: 10.1021/acs.jcim.6b00003] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Johannes G. E. M. Fraaije
- Leiden
Institute of Chemistry, Leiden University, Einsteinweg 55, 2300 RA Leiden, The Netherlands
- Culgi BV, Galileiweg 8, 2333 BD Leiden, The Netherlands
| | - Jan van Male
- Culgi BV, Galileiweg 8, 2333 BD Leiden, The Netherlands
| | - Paul Becherer
- Culgi BV, Galileiweg 8, 2333 BD Leiden, The Netherlands
| | | |
Collapse
|
34
|
Lee MT, Mao R, Vishnyakov A, Neimark AV. Parametrization of Chain Molecules in Dissipative Particle Dynamics. J Phys Chem B 2016; 120:4980-91. [DOI: 10.1021/acs.jpcb.6b00031] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Ming-Tsung Lee
- Department of Chemical and
Biochemical Engineering, Rutgers, The State University of New Jersey, 98 Brett Road, Piscataway, New Jersey 08854-8058, United States
| | - Runfang Mao
- Department of Chemical and
Biochemical Engineering, Rutgers, The State University of New Jersey, 98 Brett Road, Piscataway, New Jersey 08854-8058, United States
| | - Aleksey Vishnyakov
- Department of Chemical and
Biochemical Engineering, Rutgers, The State University of New Jersey, 98 Brett Road, Piscataway, New Jersey 08854-8058, United States
| | - Alexander V. Neimark
- Department of Chemical and
Biochemical Engineering, Rutgers, The State University of New Jersey, 98 Brett Road, Piscataway, New Jersey 08854-8058, United States
| |
Collapse
|
35
|
Li X, Gao L, Fang W. Dissipative Particle Dynamics Simulations for Phospholipid Membranes Based on a Four-To-One Coarse-Grained Mapping Scheme. PLoS One 2016; 11:e0154568. [PMID: 27137463 PMCID: PMC4854440 DOI: 10.1371/journal.pone.0154568] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2016] [Accepted: 04/15/2016] [Indexed: 11/19/2022] Open
Abstract
In this article, a new set of parameters compatible with the dissipative particle dynamics (DPD) force field is developed for phospholipids. The coarse-grained (CG) models of these molecules are constructed by mapping four heavy atoms and their attached hydrogen atoms to one bead. The beads are divided into types distinguished by charge type, polarizability, and hydrogen-bonding capacity. First, we derive the relationship between the DPD repulsive force and Flory-Huggins χ-parameters based on this four-to-one CG mapping scheme. Then, we optimize the DPD force parameters for phospholipids. The feasibility of this model is demonstrated by simulating the structural and thermodynamic properties of lipid bilayer membranes, including the membrane thickness, the area per lipid, the lipid tail orientation, the bending rigidity, the rupture behavior, and the potential of mean force for lipid flip-flop.
Collapse
Affiliation(s)
- Xiaoxu Li
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Lianghui Gao
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Weihai Fang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| |
Collapse
|
36
|
|
37
|
Oyarzún Rivera B, van Westen T, Vlugt TJH. Liquid-crystal phase equilibria of Lennard-Jones chains. Mol Phys 2016. [DOI: 10.1080/00268976.2015.1134824] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
Affiliation(s)
| | - Thijs van Westen
- Process and Energy Laboratory, Delft University of Technology, Delft, The Netherlands
| | - Thijs J. H. Vlugt
- Process and Energy Laboratory, Delft University of Technology, Delft, The Netherlands
| |
Collapse
|
38
|
Li NK, Fuss WH, Tang L, Gu R, Chilkoti A, Zauscher S, Yingling YG. Prediction of solvent-induced morphological changes of polyelectrolyte diblock copolymer micelles. SOFT MATTER 2015; 11:8236-8245. [PMID: 26315065 DOI: 10.1039/c5sm01742d] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/04/2023]
Abstract
Self-assembly processes of polyelectrolyte block copolymers are ubiquitous in industrial and biological processes; understanding their physical properties can also provide insights into the design of polyelectrolyte materials with novel and tailored properties. Here, we report systematic analysis on how the ionic strength of the solvent and the length of the polyelectrolyte block affect the self-assembly and morphology of the polyelectrolyte block copolymer materials by constructing a salt-dependent morphological phase diagram using an implicit solvent ionic strength (ISIS) method for dissipative particle dynamics (DPD) simulations. This diagram permits the determination of the conditions for the morphological transition into a specific shape, namely vesicles or lamellar aggregates, wormlike/cylindrical micelles, and spherical micelles. The scaling behavior for the size of spherical micelles is predicted, in terms of radius of gyration (R(g,m)) and thickness of corona (Hcorona), as a function of solvent ionic strength (c(s)) and polyelectrolyte length (NA), which are R(g,m) ∼ c(s)(-0.06)N(A)(0.54) and Hcorona ∼ c(s)(-0.11)N(A)(0.75). The simulation results were corroborated through AFM and static light scattering measurements on the example of the self-assembly of monodisperse, single-stranded DNA block-copolynucleotides (polyT50-b-F-dUTP). Overall, we were able to predict the salt-responsive morphology of polyelectrolyte materials in aqueous solution and show that a spherical-cylindrical-lamellar change in morphology can be obtained through an increase in solvent ionic strength or a decrease of polyelectrolyte length.
Collapse
Affiliation(s)
- Nan K Li
- Department of Materials Science and Engineering, North Carolina State University, 911 Partners Way, Raleigh, NC 27695, USA.
| | - William H Fuss
- Department of Materials Science and Engineering, North Carolina State University, 911 Partners Way, Raleigh, NC 27695, USA.
| | - Lei Tang
- Department of Mechanical Engineering and Materials Science, Duke University, 144 Hudson Hall, Durham, NC 27708, USA
| | - Renpeng Gu
- Department of Mechanical Engineering and Materials Science, Duke University, 144 Hudson Hall, Durham, NC 27708, USA
| | - Ashutosh Chilkoti
- Department of Biomedical Engineering, Duke University, 136 Hudson Hall, Durham, NC 27708, USA
| | - Stefan Zauscher
- Department of Mechanical Engineering and Materials Science, Duke University, 144 Hudson Hall, Durham, NC 27708, USA
| | - Yaroslava G Yingling
- Department of Materials Science and Engineering, North Carolina State University, 911 Partners Way, Raleigh, NC 27695, USA.
| |
Collapse
|
39
|
Theoretical description of 2D-cluster formation of nonionic surfactants at the air/water interface. Colloid Polym Sci 2015. [DOI: 10.1007/s00396-015-3630-8] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
40
|
Oyarzún B, van Westen T, Vlugt TJH. Isotropic-nematic phase equilibria of hard-sphere chain fluids—Pure components and binary mixtures. J Chem Phys 2015; 142:064903. [DOI: 10.1063/1.4907639] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023] Open
|
41
|
Ren Y, Liu B, Kiryutina T, Xi H, Qian Y. Investigation of structure formation mechanism of a mesoporous ZSM-5 zeolite by mesoscopic simulation. Chem Phys 2015. [DOI: 10.1016/j.chemphys.2014.12.011] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
42
|
Liyana-Arachchi TP, Jamadagni SN, Eike D, Koenig PH, Siepmann JI. Liquid–liquid equilibria for soft-repulsive particles: Improved equation of state and methodology for representing molecules of different sizes and chemistry in dissipative particle dynamics. J Chem Phys 2015; 142:044902. [DOI: 10.1063/1.4905918] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Thilanga P. Liyana-Arachchi
- 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
| | - David Eike
- Computational Chemistry, Modeling and Simulation, The Procter & Gamble Company, 8611 Beckett Road, West Chester, Ohio 45069, USA
| | - Peter H. Koenig
- 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
- Department of Chemical Engineering and Materials Science, University of Minnesota, 421 Washington Avenue SE, Minneapolis, Minnesota 55455, USA
| |
Collapse
|
43
|
Sethuraman V, Nguyen BH, Ganesan V. Coarse-graining in simulations of multicomponent polymer systems. J Chem Phys 2014; 141:244904. [DOI: 10.1063/1.4904390] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
| | - Bryan H. Nguyen
- 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
|
44
|
Cheng J, Vishnyakov A, Neimark AV. Morphological transformations in polymer brushes in binary mixtures: DPD study. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2014; 30:12932-40. [PMID: 25295697 DOI: 10.1021/la503520e] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/26/2023]
Abstract
Morphological transformations in polymer brushes in a binary mixture of good and bad solvents are studied using dissipative particle dynamics simulations drawing on a characteristic example of polyisoprene natural rubber in an acetone-benzene mixture. A coarse-grained DPD model of this system is built based on the experimental data in the literature. We focus on the transformation of dense, collapsed brush in bad solvent (acetone) to expanded brush solvated in good solvent (benzene) as the concentration of benzene increases. Compared to a sharp globule-to-coil transition observed in individual tethered chains, the collapsed-to-expanded transformation in brushes is found to be gradual without a prominent transition point. The transformation becomes more leveled as the brush density increases. At low densities, the collapsed brush is highly inhomogeneous and patterned into bunches composed of neighboring chains due to favorable polymer-polymer interaction. At high densities, the brush is expanded even in bad solvent due to steric restrictions. In addition, we considered a model system similar to the PINR-acetone-benzene system, but with the interactions between the solvent components worsened to the limit of miscibility. Enhanced contrast between good and bad solvents facilitates absorption of the good solvent by the brush, shifting the collapsed-to-expanded transformation to lower concentrations of good solvent. This effect is especially pronounced for higher brush densities.
Collapse
Affiliation(s)
- Jianli Cheng
- Department of Chemical and Biochemical Engineering, Rutgers, The State University of New Jersey , 98 Brett Road, Piscataway New Jersey 08854, United States
| | | | | |
Collapse
|
45
|
Li NK, Fuss WH, Yingling YG. An Implicit Solvent Ionic Strength (ISIS) Method to Model Polyelectrolyte Systems with Dissipative Particle Dynamics. MACROMOL THEOR SIMUL 2014. [DOI: 10.1002/mats.201400043] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- Nan K. Li
- Department of Materials Science and Engineering; North Carolina State University; 911 Partners Way Raleigh NC 27695 USA
| | - William H. Fuss
- Department of Materials Science and Engineering; North Carolina State University; 911 Partners Way Raleigh NC 27695 USA
| | - Yaroslava G. Yingling
- Department of Materials Science and Engineering; North Carolina State University; 911 Partners Way Raleigh NC 27695 USA
| |
Collapse
|
46
|
Lee MT, Vishnyakov A, Neimark AV. Calculations of Critical Micelle Concentration by Dissipative Particle Dynamics Simulations: The Role of Chain Rigidity. J Phys Chem B 2013; 117:10304-10. [DOI: 10.1021/jp4042028] [Citation(s) in RCA: 81] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Affiliation(s)
- Ming-Tsung Lee
- Department of Chemical
Engineering, Rutgers, The State University of New Jersey, 98 Brett Road, Piscataway,
New Jersey 08854, United States
| | - Aleksey Vishnyakov
- Department of Chemical
Engineering, Rutgers, The State University of New Jersey, 98 Brett Road, Piscataway,
New Jersey 08854, United States
| | - Alexander V. Neimark
- Department of Chemical
Engineering, Rutgers, The State University of New Jersey, 98 Brett Road, Piscataway,
New Jersey 08854, United States
| |
Collapse
|
47
|
Sandhu P, Zong J, Yang D, Wang Q. On the comparisons between dissipative particle dynamics simulations and self-consistent field calculations of diblock copolymer microphase separation. J Chem Phys 2013; 138:194904. [DOI: 10.1063/1.4804608] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
48
|
Schulz SG, Frieske U, Kuhn H, Schmid G, Müller F, Mund C, Venzmer J. The Self-Assembly of an Amphiphilic Block Copolymer: A Dissipative Particle Dynamics Study. TENSIDE SURFACT DET 2013. [DOI: 10.3139/113.100254] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
Abstract
Abstract
The self-assembly of a poly(ethylene oxide)-block-poly(propylene oxide)-block-poly(ethylene oxide) blockcopolymer (EO13PO30EO13) in the presence of water was investigated. The phases have been studied with the mesoscopic simulation technique Dissipative Particle Dynamics (DPD). The micellar, hexagonal, lamellar, bicontinuous and inverse micellar phase have been identified with this method and showed remarkable agreement with the experimental phase behavior. This work showed that the Dissipative Particle Dynamics simulations are also suitable for complicated polymer systems such as the EO13PO30EO13 triblock copolymer.
Collapse
Affiliation(s)
- S. G. Schulz
- CAM-D Technologies GmbH, Essen, Germany
- University Duisburg-Essen, Essen, Germany
| | - U. Frieske
- CAM-D Technologies GmbH, Essen, Germany
- University Duisburg-Essen, Essen, Germany
| | - H. Kuhn
- CAM-D Technologies GmbH, Essen, Germany
- University Duisburg-Essen, Essen, Germany
| | - G. Schmid
- University Duisburg-Essen, Essen, Germany
| | - F. Müller
- Degussa, Goldschmidt GmbH, Home Care, Essen, Germany
| | - C. Mund
- Degussa, Goldschmidt GmbH, Oligomeres & Silicones, Essen, Germany
| | - J. Venzmer
- Degussa, Goldschmidt GmbH, Oligomeres & Silicones, Essen, Germany
| |
Collapse
|
49
|
Vishnyakov A, Lee MT, Neimark AV. Prediction of the Critical Micelle Concentration of Nonionic Surfactants by Dissipative Particle Dynamics Simulations. J Phys Chem Lett 2013; 4:797-802. [PMID: 26281935 DOI: 10.1021/jz400066k] [Citation(s) in RCA: 118] [Impact Index Per Article: 10.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/08/2023]
Abstract
Micellization of surfactant solutions is a ubiquitous phenomenon in natural systems and technological processes, and its theoretical description represents one of the cornerstone problems in the physical chemistry of colloidal systems. However, successful attempts of quantitative modeling confirmed by experimental data remains limited. We show, for the first time, that the dissipative particle dynamics with rigorously defined soft repulsion interaction and rigidity parameters is capable of predicting micellar self-assembly of nonionic surfactants. This is achieved due to a novel approach suggested for defining the interaction parameters by fitting to the infinite dilution activity coefficients of binary solutions formed by reference compounds that represent coarse-grained fragments of surfactant molecules. Using this new parametrization scheme, we obtained quantitative agreement with the experimental critical micelle concentration and aggregation number for several typical surfactants of different chemical structures. The proposed approach can be extended to various colloidal and polymeric systems beyond nonionic surfactant solutions.
Collapse
Affiliation(s)
- Aleksey Vishnyakov
- Department of Chemical Engineering, Rutgers, the State University of New Jersey, 98 Brett Road, Piscataway New Jersey 08854, United States
| | - Ming-Tsung Lee
- Department of Chemical Engineering, Rutgers, the State University of New Jersey, 98 Brett Road, Piscataway New Jersey 08854, United States
| | - Alexander V Neimark
- Department of Chemical Engineering, Rutgers, the State University of New Jersey, 98 Brett Road, Piscataway New Jersey 08854, United States
| |
Collapse
|
50
|
Zhang C, Hu J, Chen D, Zhu Y, Fan Y, Liu Y. Does competitive hydrogen bonding induce self-assembly of A-b-B/C blend? POLYMER 2012. [DOI: 10.1016/j.polymer.2012.08.027] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|