Projection of atomistic simulation data for the dynamics of entangled polymers onto the tube theory: Calculation of the segment survival probability function and comparison with modern tube models Pavlos S. Stephanou, Chunggi Baig, * and Vlasis G. Mavrantzas * Department of Chemical Engineering, University of Patras & FORTH-ICE/HT, Patras, GR 654, Greece Supporting Information * Authors to whom correspondence should be addressed. Email: cbaig@iceht.forth.gr (Chunggi Baig), Tel.: +3-6-96955, FAX: +3-6-9653 Email: vlasis@chemeng.upatras.gr (Vlasis Mavrantzas), Tel.: +3-6-997398, FAX: +3-6-9653, Web site: http://lstm.chemeng.upatras.gr.
In this supplementary material we present the details of the potential models employed in our molecular dynamics (MD) simulations for polyethylene (PE) and cis-,4 and trans-,4 polybutadiene (PB) melts. The reader is also directed to Refs S-S3 for additional information.. Potential model for polyethylene melts All MD simulations of PE melts have been conducted in isobaric-isothermal (NPT) statistical ensemble at T=45 K and P= atm using the LAMMPS (large-scale atomic/molecular massively parallel simulator) software. S4 For the numerical integration of the equations of motion, the reversible Reference System Propagator Algorithm (r-respa) S5 was employed with the small integration time step set as equal to fs and the large one to 5 fs. The Nosé-Hoover thermostat S6 and the Andersen barostat S7 were used to control the temperature and the pressure, respectively. A standard united-atom description is used, according to which all hydrogen atoms attached to a carbon atom (e.g., CH and CH 3 ) are lumped into a single spherically symmetric interaction site. All intramolecular interactions between sites separated by more than three bonds along the chain as well as all intermolecular ones between sites belonging to different chains are described via a -6 Lennard- Jones (LJ) potential: VLJ ( r) 6 σ σ 4ε =, (S) r r with ε and σ for the CH and CH 3 being taken from the well-known TraPPE model: S8 ε=.93 kcal.mol - and σ=3.95 Å. Following Nath et al., S9 bond lengths are allowed to fluctuate via a harmonic potential: V () ( ) stretching l = kl l l, (S) where l =.54 Å and k l =9.55 kcal mol - Berendsen S bending potential: Å -. Bond bending interactions obey the Van der Ploeg- Vbending = k θ, (S3) ( θ ) ( θ θ )
where θ =4 ο and k θ =4.6 kcal mol - rad -. Finally, the torsional potential employed follows Toxvaerd: S V torsional 8 ( φ ) c cos n ( φ ) =, (S4) n= n with c =.9869, c =4.8, c =-.65, c 3 =-7.698, c 4 =4.46, c 5 =3.95, c 6 =-8.97, c 7 =-3.446 and c 8 =5.597 (all in units of kcal mol - ).. Potential model for cis-,4 and trans-,4 polybutadiene melts All MD simulations of cis-,4 and trans-,4 PB melts have also been conducted in the NPT statistical ensemble at T=43 K and P= atm using the LAMMPS (large-scale atomic/molecular massively parallel simulator) software. S4 The r-respa S5 was employed with the small integration time step set as equal to fs and the long to 5 fs. The Nosé-Hoover thermostat and barostat S6 were used to control the temperature and the pressure, respectively. A united-atom description (e.g., CH, CH, and CH 3 ) was employed. All intramolecular interactions between sites separated by more than three bonds along the chain as well as all intermolecular ones are described via a -6 Lennard-Jones potential (eq S) with ε and σ being taken from the work of Smith and Paul S [see (a) in Table S]. Bond lengths fluctuate via the harmonic potential (eq S) where the values of l and k l are obtained from Gee and Boyd S3 and the bond angles again follow the harmonic potential (eq S3) with the values of θ and k θ obtained from Gee and Boyd S3, which are the same as in Smith and Paul S [see (b) in Table S]. Finally, the torsional potential is obtained from Smith and Paul S 6 Vtorsional ( φ ) = cn ( cos( nφ )), (S5) n= where the parameters c n are given in Table S. Table S: Values of the Force Field parameters for the simulation of PB (a) bond-stretching types k l (kcal.mol -. Å - ) l (Å) CH -CH 58.5.54 CH -CH 83.8.5 3
CH=CH 46.9.34 (b) bond-bending types k θ (kcal.mol -. Å - ) θ ( ο ) CH -CH -CH 5.65 CH -CH 89.4 5.89 (c) bond-torsional types k k k 3 k 4 k 5 k 6 (kcal.mol - ) (kcal.mol - ) (kcal.mol - ) (kcal.mol - ) (kcal.mol - ) (kcal.mol - ) CH -CH=CH-CH 4. CH -CH -CH=CH.33 -.47.554.63.346.64 CH- CH -CH -CH -.888 -.69-3.639 -.66 -.47 -.9 (d) non-bonded Lennard-Jones types ε (kcal.mol - ) σ (Å) CH -CH.963 4.5 CH -CH. 3.8 CH=CH.5 4.57 The initial configuration in all the simulations of PE and PB systems was obtained from the Cerius software S4 and was then subjected to a minimization process to minimize the total energy of the initial configuration. Constant-volume NVT and constant-pressure NPT Monte Carlo simulations (e.g. ref S5) were subsequently employed to further equilibrate the initial system configuration; the resulting configuration was then used in the MD simulations in this study. This initial system configuration was carefully chosen to have a uniform chain-length distribution, a density close to the average melt density of the simulated system at the temperature and pressure of interest, and a mean-square chain end-to-end distance R close to the real average value of the simulated system. 5 Conventional periodic boundary conditions were applied in all three directions (x, y, z). 6 References and Notes (S) V. A. Harmandaris, V. G. Mavrantzas and D. N. Theodorou, Macromolecules,, 33, 86. (S) N. Ch. Karayiannis and V. G. Mavrantzas, Macromolecules, 5, 38, 8593. 4
(S3) G. Tsolou, V. G. Mavrantzas and D. N. Theodorou, Macromolecules, 5, 38, 478; G. Tsolou, V. A. Harmandaris and V. G. Mavrantzas, Macromol. Theory Simul., 6, 5, 38 (note that the values of k str reported in Table of the Macromolecules paper are in error and the corrected values are presented in Table of the Macromol. Theory Simul. paper) (S4) S. J. Plimton, J. Comput. Phys., 995, 7,. (S5) M. Tuckerman, B. J. Berne and G. J. Martyna, J. Chem. Phys., 99, 97, 99; G. J. Martyna, M. Tuckerman, D. J. Tobias and M. L. Klein, Mol. Phys., 996, 87, 7. (S6) S. Nosé, Prog. Theor. Phys. Suppl., 99, 3, ; W. G. Hoover, Phys. Rev. A, 986, 3, 695. (S7) H. C. Andersen, J. Chem. Phys., 98, 7, 384. (S8) M. G. Martin and J. I. Siepmann, J. Phys. Chem. B., 998,, 569. (S9) S. K. Nath, F. A. Escobedo and J. J. de Pablo, J. Chem. Phys., 998, 8, 995; S. K. Nath and R. Khare, J. Chem. Phys.,, 5, 837. (S) P. Van der Ploeg and H. J. C. Berendsen, J. Chem. Phys., 98, 76, 37. (S) S. Toxvaerd, J. Chem. Phys., 997, 7, 597. (S) G. D. Smith and W. Paul, J. Phys. Chem. A, 998,,. (S3) R. H. Gee and R. H. Boyd, J. Chem. Phys., 994,, 88. (S4) www.accelrys.com/cerius/ (S5) P. Gestoso, E. Nicol, M. Doxastakis and D. N. Theodorou, Macromolecules., 3, 36, 695. (S6) M. P. Allen and D. J. Tildesley, Computer Simulation of Liquids, Clarendon Press, Oxford, 987. 5