Simulation of the plastic behavior of amorphous glassy bis-phenol-A- polycairbonate

Journal of Chemical Physics

Volumn 121, Issue 10, 2004, Pages 4941-4950

  Fortunelli, Alessandro ^a   Geloni, Claudio ^b   Lazzeri, Andrea ^b  

^a CNR   (Italy)

^b UNIVERSITY OF PISA   (Italy)

Author keywords

[No Author keywords available]

Indexed keywords

ACTIVATION ENERGY; CHEMICAL BONDS; COMPUTATIONAL METHODS; COMPUTER SIMULATION; ELASTIC MODULI; GLASS; KINETIC THEORY; MOLECULAR DYNAMICS; OPTIMIZATION; PERCOLATION (FLUIDS); PHENOLS; RELAXATION PROCESSES; SHEAR FLOW; STRAIN RATE;

AMORPHOUS GLASSY POLYMERS; ENERGY MINIMIZATION; MOLECULAR MECHANICS; PLASTIC TRANSFORMATION;

POLYCARBONATES;

EID: 4544313845     PISSN: 00219606     EISSN: None     Source Type: Journal    
DOI: 10.1063/1.1778157     Document Type: Article

References (86)

1. P. Ballone, B. Montanari, R. O. Jones, and O. Hahn, J. Phys. Chem. A103, 5387 (1999).
2. S. L. Mayo, B. D. Olafson, and W. A. Goddard, J. Phys. Chem. 94, 8897(1990).
3. S. W. Bunte and H. Sun, J. Phys. Chem. B 104, 2477 (2000).
4. H. Sun, J. Phys. Chem. B 102, 7338 (1998).
5. H. Sun, P. Ren, and J. R. Fried, Comput. Theor. Polym. Sci. 8, 229 (1998).
6. H. Sun, Macromolecules 28, 701 (1995).
7. H. Sun, Macromolecules 26, 5924 (1993).
8. H. Sun, J. Comput. Chem. 15, 752 (1994).
9. D. N. Theodorou and U. W. Suter, Macromolecules 18, 1467 (1985).
10. D. N. Theodorou and U. W. Suter, Macromolecules 19, 139 (1986).
11. D. N. Theodorou and U. W. Suter, Macromolecules 19, 379 (1986).
12. P. H. Mott, A. S. Argon, and U. W. Suter, Philos. Mag. A 67, 931 (1993).
13. M. Hutnik, A. S. Argon, and U. W. Suter, Macromolecules 24, 5956 (1991).
14. M. Hutnik, F. T. Gentile, P. J. Ludovice, U. W. Suter, and A. S. Argon, Macromolecules 24, 5962 (1991).
15. M. Hutnik, A. S. Argon, and U. W. Suter, Macromolecules 24, 5970 (1991).
16. M. Hutnik, A. S. Argon, and U. W. Suter, Macromolecules 26, 1097 (1993).
17. M. Utz, P. Robyr, and U. W. Suter, Macromolecules 33, 6808 (2000).
18. C. F. Fan, T. Caǧin, Z. M. Chen, and K. A. Smith, Macromolecules 27, 2383 (1994).
19. C. F. Fan, Macromolecules 28, 5215 (1995).
20. A. F. Yee, L. Yang, and D. J. Srolovitz, J. Chem. Phys. 107, 4396 (1997).
21. S. S. Jang and W. H. Jo, J. Chem. Phys. 110, 7524 (1999).
22. F. M. Capaldi, M. C. Boyce, and G. C. Rutledge, Phys. Rev. Lett. 89, 175505 (2002).
23. F. M. Capaldi, M. C. Boyce, and G. C. Rutledge, Polymer 45, 1391 (2004).
24. M. Goldstein, J. Chem. Phys. 51, 3728 (1969).
25. E. Orowan, in Proceedings of the First National Congress of Applied Mechanics, (American Society of Mechanical Engineers, New York, 1952), pp. 453-472.
26. F. H. Stillinger and T. A. Weber, Science 225, 983 (1984).
27. F. H. Stillinger, Science 267, 1935 (1995).
28. N. P. Kopsias and D. N. Theodorou, J. Chem. Phys. 109, 8573 (1998).
29. M. L. Greenfield and D. N. Theodorou, Mol. Simul. 19, 329 (1997).
30. T. F. Middleton and D. J. Wales, Phys. Rev. B 64, 024205 (2001).
31. F. H. Stillinger and P. G. DeBenedetti, J. Chem. Phys. 116, 3353 (2002).
32. D. L. Malandro and D. J. Lacks, J, Chem. Phys. 110, 4593 (1999).
33. M. Parrinello and A. Rahman, Phys. Rev. Lett. 45, 1196 (1980).
34. M. Parrinello and A. Rahman, J. Appl. Phys. 52, 7182 (1981).
35. M. Parrinello and A. Rahman, J. Chem. Phys. 76, 2662 (1982).
36. A. S. Argon and M. I. Bessanov, Philos. Mag. A 35, 917 (1977)
37. C. F. Abrams and K. Kremer, Macromolecules 36, 260 (2003).
38. Aecelrys Inc., CERIUS2, Release 4.0, San Diego: Accelrys Inc., 1999.
39. P. Procacci, T. Darden, and M. Marchi, J. Comp. Chem. 18, 1848 (1997).
40. M. Clark, R. D. Cramer III, and N. Van Opdenbosch, J. Comput. Chem. 10, 982 (1989).
41. The multistep algorithm of the ORAC program proved very useful in speeding up the MD simulations, allowing one to use a time step of 10 fsec at 300 K and 4.5 fsec at 1500 K.
42. Courtesy of Dr. G. Marchetti and Dr. M. Marchi (Orsay, France).
43. Accelrys Inc., MATERIALS STUDIO, Release 2.1 (Accelrys Inc., San Diego, 2002).
44. Results analogous to those described below were found when using different starting configurations, obtained by annealing the system at 1500 K for longer times.
45. The Discover® program does not allow one to perform "strained" energy minimizations - i.e., minimizations in which one value of the cell parameter is held fixed, while the others are optimized, or "strained" MD simulations - i.e., in which one value of the cell parameters is held fixed while the others are allowed to vary, following the Parrinello-Rahman dynamics (see Ref, 20).
46. We chose a small value of the fictitious cell mass to increase the speed of approaching equilibrium in the MD runs. A discussion of the optimal value of the fictitious cell mass is deferred to future work.
47. This may be thought of as corresponding to a sequence of creep steps at increasing stress levels (Ref. 48) conducted at a very high strain speed.
48. N. McCrum, C. P. Buckley, and C. B. Bucknall, Principles of Polymer Engineering (Oxford University Press, Oxford, 1988), Chap. 4.
49. I. M. Ward, Mechanical Properties of Solid Polymers (Wiley, New York,1985), Chaps. 5 and 6.
50. Handbook of Chemistry and Physics, edited by D. R. Lide, 72th ed. (CRC, Boca Raton, 1991).
51. I. V. Yannas and R. R. Luise, J. Macromol. Sci., Phys. B21, 443 (1982).
52. R. J. Yauag and P. A. Lovell, Introduction to Polymers (Chapman & Hall, New York, 1991), Chap. 5.
53. R. N. Haward and R. J. Young, The Physics of Glassy Polymers (Chapman & Hall, New York, 1997).
54. C. G'Sell and J. J. Jonas, J. Mater. Sci. 14, 583 (1979).
55. J. L. Bucaille, E. Felder, and G. Hochstetter, J. Mater, Sci. 37, 3999 (2002).
56. D. D. Dlott, Acc. Chem. Res. 33, 37 (2000).
57. Other structures exibit strain hardening at larger values of strain.
58. J. H. Weiner, Statistical Mechanics of Elasticity (Dover, New York, 2002), Chaps. 1 and 4.
59. The new structure obtained by bringing a given ideal structure outside its range of elastic behavior can be used as the configuration of subsequent MD simulations plus "sweeping procedures" to generate other ideal structures, somewhat related to (but not identical with) the original one.
60. One can easily devise more rigorous procedures to attain this result, such as local variants of global optimization techniques, etc. (Ref. 30). The "ideal structures" are in fact representative of "catchment basins" within which ergodicity is recovered in the glassy state, with a stability (entropic component of the free energy) proportional to the volume of the phase space associated with them. The sweeping procedure can be seen as a rough but very economical and (in the present case) rather efficient approximate version of these general approaches.
61. Those rearrangements which are only weakly coupled to the overall strain cannot be discriminated by the sweeping procedure, but they are also less interesting in this context, as they represent "floppy" local modes of "soft" regions of the system.
62. R. Zallen, The Physics of Amorphous Solids (Wiley, New York, 1983).
63. The Polak-Ribière (Réf. 64) algorithm (a conjugate gradient method) was used for the energy minimizations, with a stopping criterion for maximum derivative of 0.01 kcal/mol/Å.
64. J. C. Gilbert and J. Nocedal, SIAM J. Optim. 1, 21 (1992).
65. The Discover® program does not allow one to calculate the Hessian matrix for a system of 10 000 atoms, so that we were unable to implement the normal-mode analysis of Ref. 32.
66. H. Eyring, J. Chem. Phys. 4, 283 (1936).
67. Actually, this reasonable assumption brings us outside the original formulation of the Eyring model (Ref. 66).
68. L. P. Chen, A. F. Yee, and E. J. Moskala, Macromolecules 32, 5944 (1999).
69. S. S. Lee and A. F. Yee, Macromolecules 36, 6791 (2003).
70. E. Oleynik, Prog. Colloid Polym. Sci. 80, 140 (1989).
71. M. Nandagopal and M. Utz, J. Chem. Phys. 118, 8373 (2003).
72. J. Lubliner, Plasticity Theory (Macmillan, New York, 1990), p. 76.
73. C. G'Sell, H. El Bari, J. Perez, J. Y. Cavaille and G. P. Johari, Mater. Sci. Eng., A 110, 223 (1989).
74. In the case that the cell parameters are different, the choice of a threshold value allows us to get rid of spurious effects due to the difference between affine transformation and elastic response. Better results could be obtained by keeping stiff intramolecular degrees of freedom rigid.
75. A rough "coarse-grained" check was also performed, by monitoring which phenyl-carbonate-phenyl7 segments underwent a syn-anti conformational transition.
76. It can be noted that the Discover® program does not calculate atomic stresses (as defined in Ref. 11). They can be obtained via numerical differentiation of the atomic energies, without, however, conveying substantially new information in the present case.
77. The average length of these segments was about two to three monomeric units, with a maximum length of ten units.
78. From previous work (Ref. 16), it is known that in BPA-PC each barrier typically has an activation energy of several kcal/mol, and there are 80 of such motions.
79. A. Fortunelli, C. Geloni, A. Lazzeri, and G. Marchetti (unpublished).
80. G. Adam and J. H. Gibbs, J. Chem. Phys. 43, 139 (1965).
81. S.Müller, in Calculus of Variations and Geometric Evolution Problems, edited by F. Bethuel, G. Huisken, S. Müller, K. Steffen, K. Hildebrandt, and M. Struwe (Springer, Berlin, 1999).
82. M. Utz, P. G. DeBenedetti, and F. H. Stillinger, Phys. Rev. Lett. 84, 1471 (2000).
83. M. Ortiz, J. Appl. Mech. 66, 289 (1999).
84. G. B. McKenna, J. Phys.: Condens. Matter 15, S737 (2003).
85. M. Doi and S. F. Edwards, The Theory of Polymer Dynamics (Oxford Science Publications, Oxford, 2001), Chap. 6.
86. D. L. Malandro and D. J. Lacks, Comput. Theor. Polym. Sci. 110, 4593 (1999).