1
|
Wu Z, Müller-Plathe F. Slip-Spring Hybrid Particle-Field Molecular Dynamics for Coarse-Graining Branched Polymer Melts: Polystyrene Melts as an Example. J Chem Theory Comput 2022; 18:3814-3828. [PMID: 35617016 DOI: 10.1021/acs.jctc.2c00107] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The topology of chains significantly modifies the dynamical properties of polymer melts. Here, we extend a recently developed efficient simulation method, namely the slip-spring hybrid particle-field (SS-hPF) model, to study the structural and dynamical properties of branched polymer melts over large spatial-temporal scales. In the coarse-grained SS-hPF simulation of polymers, the bonded potentials are derived by iterative Boltzmann inversion from the underlying fine-grained model. The nonbonded potentials are computed from a density functional field instead of pairwise interactions used in standard molecular dynamics simulations, which increases the computational efficiency by a factor of 10-20. The entangled dynamics is lost due to the soft-core nature of density functional field interactions. It is recovered by a multichain slip-spring model that is rigorously parametrized from existing experimental or simulation data. To quantitatively predict the relaxation and diffusion of branched polymers, which are dominated by arm retraction rather than chain reptation, the slip-spring algorithm is augmented to improve the polymer dynamics near the branch point. Multiple dynamical observables, e.g., diffusion coefficients, arm relaxations, and tube survival probabilities, are characterized in an example coarse-grained model of symmetric and asymmetric star-shaped polystyrene melts. Consistent dynamical behaviors are identified and compared with theoretical predictions. With a single rescaling factor, the prediction of diffusion coefficients agrees well with the available experimental measurements. In this work, an efficient approach is provided to build chemistry-specific coarse-grained models for predicting the dynamics of branched polymers.
Collapse
Affiliation(s)
- Zhenghao Wu
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, Alarich-Weiss-Strasse 8, 64287 Darmstadt, Germany
| | - Florian Müller-Plathe
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, Alarich-Weiss-Strasse 8, 64287 Darmstadt, Germany
| |
Collapse
|
2
|
Masubuchi Y, Yang L, Uneyama T, Doi Y. Analysis of Elongational Viscosity of Entangled Poly (Propylene Carbonate) Melts by Primitive Chain Network Simulations. Polymers (Basel) 2022; 14:741. [PMID: 35215654 PMCID: PMC8874545 DOI: 10.3390/polym14040741] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2021] [Revised: 02/07/2022] [Accepted: 02/11/2022] [Indexed: 11/28/2022] Open
Abstract
It has been established that the elongational rheology of polymers depends on their chemistry. However, the analysis of experimental data has been reported for only a few polymers. In this study, we analyzed the elongational viscosity of poly (propylene carbonate) (PPC) melts in terms of monomeric friction via primitive chain network simulations. By incorporating a small polydispersity of materials, the linear viscoelastic response was semi-quantitatively reproduced. Owing to this agreement, we determined units of time and modulus to carry out elongational simulations. The simulation with constant monomeric friction overestimated elongational viscosity, whereas it nicely captured the experimental data if friction decreased with increasing segment orientation. To see the effect of chemistry, we also conducted the simulation for a polystyrene (PS) melt, which has a similar entanglement number per chain and a polydispersity index. The results imply that PPC and PS behave similarly in terms of the reduction of friction under fast deformations.
Collapse
Affiliation(s)
- Yuichi Masubuchi
- Department of Materials Physics, Nagoya University, Nagoya 4648603, Japan; (L.Y.); (T.U.); (Y.D.)
| | | | | | | |
Collapse
|
3
|
Huang LH, Wu CH, Hua CC, Huang TJ. Multiscale simulations of coupled composition-stress-morphology of binary polymer blend. POLYMER 2020. [DOI: 10.1016/j.polymer.2020.122366] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
4
|
Masubuchi Y, Ianniruberto G, Marrucci G. Primitive chain network simulations for H-polymers under fast shear. SOFT MATTER 2020; 16:1056-1065. [PMID: 31859310 DOI: 10.1039/c9sm01971e] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Branchpoint Withdrawal (BPW) has been recognized as one of the important molecular mechanisms for the description of the dynamics of entangled branched polymers under fast flows. However, the relation to the other known molecular mechanisms has not been fully elucidated yet. In this study we performed primitive chain network (i.e., multi-chain slip-link) Brownian simulations for a melt of a well-characterized monodisperse polystyrene H-polymer, for which the linear viscoelasticity and shear viscosity growth curves at several shear rates are available in the literature. After confirming the consistency of the simulations with the rheological data, we used the simulations to analyze the molecular motion in detail. The results reveal that molecular tumbling occurs in branched polymers just as in linear ones, and that it is accelerated by BPW. Furthermore, BPW not only mitigates backbone stretch, as expected, but also arm stretch. However, because the transient startup viscosity is anyhow dominated by chain stretch dynamics rather than by molecular tumbling, our results rationalize the fact that pom-pom theories successfully ignore tumbling in shear flows.
Collapse
Affiliation(s)
- Yuichi Masubuchi
- Department of Materials Physics, Nagoya University, Nagoya 4648603, Japan.
| | - Giovanni Ianniruberto
- Dipartimento di Ingegneria Chimica, dei Materiali e della Produzione Industriale, Università degli Studi di Napoli "Federico II", Piazzale Tecchio 80-80125, Napoli, Italy
| | - Giuseppe Marrucci
- Dipartimento di Ingegneria Chimica, dei Materiali e della Produzione Industriale, Università degli Studi di Napoli "Federico II", Piazzale Tecchio 80-80125, Napoli, Italy
| |
Collapse
|
5
|
Cao J, Wang Z, Likhtman AE. Determining Tube Theory Parameters by Slip-Spring Model Simulations of Entangled Star Polymers in Fixed Networks. Polymers (Basel) 2019; 11:E496. [PMID: 30960480 PMCID: PMC6473678 DOI: 10.3390/polym11030496] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2019] [Revised: 03/08/2019] [Accepted: 03/11/2019] [Indexed: 11/26/2022] Open
Abstract
Dynamical properties of branched polymer melts are determined by the polymer molecular weights and architectures containing junction points. Relaxation of entangled symmetric star polymers proceeds via arm-retraction and constraint release (CR). In this work, we investigate arm-retraction dynamics in the framework of a single-chain slip-spring model without CR effect where entanglements are treated as binary contacts, conveniently modeled as virtual "slip-links", each involving two neighboring strands. The model systems are analogous to isolated star polymers confined in a permanent network or a melt of very long linear polymers. We find that the distributions of the effective primitive path lengths are Gaussian, from which the entanglement molecular weight N e , a key tube theory parameter, can be extracted. The procured N e value is in good agreement with that obtained from mapping the middle monomer mean-square displacements of entangled linear chains in slip-spring model to the tube model prediction. Furthermore, the mean first-passage (FP) times of destruction of original tube segments by the retracting arm end are collected in simulations and examined quantitatively using a theory recently developed in our group for describing FP problems of one-dimensional Rouse chains with improbable extensions. The asymptotic values of N e as obtained from the static (primitive path length) and dynamical (FP time) analysis are consistent with each other. Additionally, we manage to determine the tube survival function of star arms μ ( t ) , or equivalently arm end-to-end vector relaxation function ϕ ( t ) , through the mean FP time spectrum τ ( s ) of the tube segments after careful consideration of the inner-most entanglements, which shows reasonably good agreement with experimental data on dielectric relaxation.
Collapse
Affiliation(s)
- Jing Cao
- School of Mathematical, Physical and Computational Sciences, University of Reading, Reading RG6 6AX, UK.
| | - Zuowei Wang
- School of Mathematical, Physical and Computational Sciences, University of Reading, Reading RG6 6AX, UK.
| | - Alexei E Likhtman
- School of Mathematical, Physical and Computational Sciences, University of Reading, Reading RG6 6AX, UK
| |
Collapse
|
6
|
Affiliation(s)
- Yuichi Masubuchi
- Department of Materials Physics, Nagoya University, Nagoya, Japan
| |
Collapse
|
7
|
Masubuchi Y, Pandey A, Amamoto Y, Uneyama T. Orientational cross correlations between entangled branch polymers in primitive chain network simulations. J Chem Phys 2017; 147:184903. [DOI: 10.1063/1.5001960] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Affiliation(s)
- Yuichi Masubuchi
- Department of Materials Physics, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8603, Japan
| | - Ankita Pandey
- Department of Materials Physics, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8603, Japan
| | - Yoshifumi Amamoto
- Department of Materials Physics, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8603, Japan
| | - Takashi Uneyama
- Institute of Science and Engineering, Kanazawa University, Kakuma, Kanazawa 920-1192, Japan
| |
Collapse
|
8
|
Cao J, Wang Z. Microscopic Picture of Constraint Release Effects in Entangled Star Polymer Melts. Macromolecules 2016. [DOI: 10.1021/acs.macromol.6b00554] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Affiliation(s)
- Jing Cao
- Department of Mathematics
and Statistics, University of Reading, Whiteknights, PO Box 220, Reading RG6 6AX, U.K
| | - Zuowei Wang
- Department of Mathematics
and Statistics, University of Reading, Whiteknights, PO Box 220, Reading RG6 6AX, U.K
| |
Collapse
|
9
|
Abstract
To optimize automation for polymer processing, attempts have been made to simulate the flow of entangled polymers. In industry, fluid dynamics simulations with phenomenological constitutive equations have been practically established. However, to account for molecular characteristics, a method to obtain the constitutive relationship from the molecular structure is required. Molecular dynamics simulations with atomic description are not practical for this purpose; accordingly, coarse-grained models with reduced degrees of freedom have been developed. Although the modeling of entanglement is still a challenge, mesoscopic models with a priori settings to reproduce entangled polymer dynamics, such as tube models, have achieved remarkable success. To use the mesoscopic models as staging posts between atomistic and fluid dynamics simulations, studies have been undertaken to establish links from the coarse-grained model to the atomistic and macroscopic simulations. Consequently, integrated simulations from materials chemistry to predict the macroscopic flow in polymer processing are forthcoming.
Collapse
Affiliation(s)
- Yuichi Masubuchi
- Institute for Chemical Research, Kyoto University, Gokasho Uji-City, Japan 611-0011
| |
Collapse
|
10
|
Masubuchi Y, Matsumiya Y, Watanabe H, Marrucci G, Ianniruberto G. Primitive Chain Network Simulations for Pom-Pom Polymers in Uniaxial Elongational Flows. Macromolecules 2014. [DOI: 10.1021/ma500357g] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Yuichi Masubuchi
- Institute
for Chemical Research, Kyoto University, Gokasyo, Uji, Kyoto 611-0011, Japan
| | - Yumi Matsumiya
- Institute
for Chemical Research, Kyoto University, Gokasyo, Uji, Kyoto 611-0011, Japan
| | - Hiroshi Watanabe
- Institute
for Chemical Research, Kyoto University, Gokasyo, Uji, Kyoto 611-0011, Japan
| | - Giuseppe Marrucci
- Dipartimento
di Ingegneria Chimica, dei Materiali e della Produzione Industriale, Università degli Studi di Napoli “Federico II”, Piazzale
Tecchio 80, 80125 Napoli, Italy
| | - Giovanni Ianniruberto
- Dipartimento
di Ingegneria Chimica, dei Materiali e della Produzione Industriale, Università degli Studi di Napoli “Federico II”, Piazzale
Tecchio 80, 80125 Napoli, Italy
| |
Collapse
|
11
|
Bačová P, Lentzakis H, Read DJ, Moreno AJ, Vlassopoulos D, Das C. Branch-Point Motion in Architecturally Complex Polymers: Estimation of Hopping Parameters from Computer Simulations and Experiments. Macromolecules 2014. [DOI: 10.1021/ma5003936] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Petra Bačová
- Centro de Física
de Materiales (CSIC, UPV/EHU) and Materials Physics Center MPC, Paseo Manuel de Lardizabal 5, E-20018 San Sebastián, Spain
| | - Helen Lentzakis
- Institute of Electronic Structure & Laser, Foundation for Research and Technology Hellas (FORTH), Heraklion, Crete 71110, Greece
- Department
of Materials Science and Technology, Univerisity of Crete, Heraklion, Crete 71300, Greece
| | - Daniel J. Read
- Department
of Applied Mathematics, University of Leeds, Leeds LS2 9JT, United Kingdom
| | - Angel J. Moreno
- Centro de Física
de Materiales (CSIC, UPV/EHU) and Materials Physics Center MPC, Paseo Manuel de Lardizabal 5, E-20018 San Sebastián, Spain
- Donostia International
Physics Center, Paseo Manuel de Lardizabal 4, E-20018 San Sebastián, Spain
| | - Dimitris Vlassopoulos
- Institute of Electronic Structure & Laser, Foundation for Research and Technology Hellas (FORTH), Heraklion, Crete 71110, Greece
- Department
of Materials Science and Technology, Univerisity of Crete, Heraklion, Crete 71300, Greece
| | - Chinmay Das
- School of
Physics and Astronomy, University of Leeds, Leeds LS2 9JT, United Kingdom
| |
Collapse
|
12
|
Milchev A, Müller M, Klushin L. Arm Retraction Dynamics and Bistability of a Three-Arm Star Polymer in a Nanopore. Macromolecules 2014. [DOI: 10.1021/ma500053n] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- A. Milchev
- Institute
for Physical Chemistry, Bulgarian Academy of Sciences, 1113 Sofia, Bulgaria
- Institut für Theoretische Physik, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| | - M. Müller
- Institut für Theoretische Physik, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
| | - L. Klushin
- Department
of Physics, American University of Beirut, Beirut, Lebanon
| |
Collapse
|
13
|
Bačová P, Hawke LGD, Read DJ, Moreno AJ. Dynamics of Branched Polymers: A Combined Study by Molecular Dynamics Simulations and Tube Theory. Macromolecules 2013. [DOI: 10.1021/ma4005988] [Citation(s) in RCA: 38] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Petra Bačová
- Departamento de Física de Materiales, Universidad del País Vasco (UPV/EHU), Apartado
1072, 20080 San Sebastián, Spain
| | | | - Daniel J. Read
- Department of Applied Mathematics, University of Leeds, LS2 9JT Leeds, U.K
| | - Angel J. Moreno
- Centro de Física de Materiales (CSIC, UPV/EHU) and Materials Physics Center MPC, Paseo
Manuel de Lardizabal 5, 20018 San Sebastián, Spain
- Donostia International Physics Center, Paseo Manuel de Lardizabal 4, 20018
San Sebastián, Spain
| |
Collapse
|
14
|
Uneyama T, Masubuchi Y. Multi-chain slip-spring model for entangled polymer dynamics. J Chem Phys 2012; 137:154902. [DOI: 10.1063/1.4758320] [Citation(s) in RCA: 90] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
15
|
Uneyama T, Masubuchi Y. Detailed balance condition and effective free energy in the primitive chain network model. J Chem Phys 2011; 135:184904. [DOI: 10.1063/1.3658775] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|