1
|
Yungerman I, Starodumov I, Fulati A, Uto K, Ebara M, Moskovitz Y. Full-Atomistic Optimized Potentials for Liquid Simulations and Polymer Consistent Force Field Models for Biocompatible Shape-Memory Poly(ε-caprolactone). J Phys Chem B 2022; 126:3961-3972. [PMID: 35605974 DOI: 10.1021/acs.jpcb.2c01973] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Thermally induced shape memory poly(ε-caprolactone) (PCL)-based polymers are one of the most extensively researched families of biocompatible materials. They are degradable under physiological conditions and have high applicability in general biomedical engineering, with cross-linked PCL networks being particularly useful for tissue engineering. In this study, we used the optimized potentials for liquid simulations (OPLS) force field, which is well suited for describing intermolecular interactions in biomolecules, and the class II polymer consistent force field (PCFF) to investigate the properties of telechelic PCL with diacrylates as reactive functionalities on its end groups. PCFF has been specifically parameterized for simulating synthetic polymeric materials. We compare the findings of all-atom molecular dynamics simulations with known experimental data and theoretical assumptions to verify the applicability of both these force fields. We estimated the melt density, volume, transition temperatures, and mechanical characteristics of two-branched PCL diacrylates with a molecular weight of 2481 Da. Our findings point to the utility of the aforementioned force fields in predicting the properties of PCL-based polymers. It also opens avenues for developing PCL cross-linked polymer models and employing OPLS to investigate the interactions of synthetic polymers with biomolecules.
Collapse
Affiliation(s)
- Irena Yungerman
- Department of Chemical Engineering, Technion-Israel Institute of Technology, Haifa 32000, Israel
| | - Ilya Starodumov
- Department of Theoretical and Mathematical Physics, Laboratory of Multi-Scale Mathematical Modeling, Ural Federal University, Ekaterinburg 620000, Russian Federation.,Ural State Medical University, Ekaterinburg 620000, Russian Federation
| | - Ailifeire Fulati
- Research Center for Functional Materials, National Institute for Materials Science (NIMS), 1-1Namiki, Tsukuba, Ibaraki 305-0044, Japan
| | - Koichiro Uto
- Research Center for Functional Materials, National Institute for Materials Science (NIMS), 1-1Namiki, Tsukuba, Ibaraki 305-0044, Japan
| | - Mitsuhiro Ebara
- Research Center for Functional Materials, National Institute for Materials Science (NIMS), 1-1Namiki, Tsukuba, Ibaraki 305-0044, Japan
| | - Yevgeny Moskovitz
- Department of Chemical Engineering, Technion-Israel Institute of Technology, Haifa 32000, Israel.,Department of Theoretical and Mathematical Physics, Laboratory of Multi-Scale Mathematical Modeling, Ural Federal University, Ekaterinburg 620000, Russian Federation
| |
Collapse
|
2
|
Bi FL, Xi ZH, Zhao L. Reaction Mechanisms and Kinetics of the Melt Transesterification of Bisphenol-A and Diphenyl Carbonate. INT J CHEM KINET 2018. [DOI: 10.1002/kin.21151] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- Feng-Lei Bi
- Shanghai Key Laboratory of Multiphase Materials Chemical Engineering; East China University of Science and Technology; Shanghai 200237 People's Republic of China
| | - Zhen-Hao Xi
- Shanghai Key Laboratory of Multiphase Materials Chemical Engineering; East China University of Science and Technology; Shanghai 200237 People's Republic of China
| | - Ling Zhao
- Shanghai Key Laboratory of Multiphase Materials Chemical Engineering; East China University of Science and Technology; Shanghai 200237 People's Republic of China
| |
Collapse
|
3
|
Bendler JT, Boyles DA, Edmondson CA, Filipova T, Fontanella JJ, Westgate MA, Wintersgill MC. Dielectric Properties of Bisphenol A Polycarbonate and Its Tethered Nitrile Analogue. Macromolecules 2013. [DOI: 10.1021/ma4002269] [Citation(s) in RCA: 55] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- John T. Bendler
- BSC, Inc., 3046 Player Drive, Rapid City,
South Dakota 57702, United States
| | - David A. Boyles
- Department of Chemistry, South Dakota School of Mines and Technology, Rapid City, South Dakota 57701, United
States
| | - Charles A. Edmondson
- Physics Department, U.S. Naval Academy, Annapolis, Maryland 21402, United
States
| | - Tsvetanka Filipova
- Department of Chemistry, South Dakota School of Mines and Technology, Rapid City, South Dakota 57701, United
States
| | - John J. Fontanella
- Physics Department, U.S. Naval Academy, Annapolis, Maryland 21402, United
States
| | - Mark A. Westgate
- Physics Department, U.S. Naval Academy, Annapolis, Maryland 21402, United
States
| | - M. C. Wintersgill
- Physics Department, U.S. Naval Academy, Annapolis, Maryland 21402, United
States
| |
Collapse
|
4
|
Sulatha MS, Natarajan U. Effect of chemical substituents on the structure of glassy diphenyl polycarbonates. J Phys Chem B 2011; 115:1579-89. [PMID: 21275412 DOI: 10.1021/jp105954z] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Polycarbonates offer a wide variety of physical property behavior that is difficult to predict due to complexities at the molecular scale. Here, the physical structure of amorphous glassy polycarbonates having aliphatic and cycloaliphatic chemical groups is explored through atomistic simulations. The influence of chemical structure on solubility parameter, torsion distributions, radial distribution function, scattering structure factor, orientation distributions of phenylene rings and carbonate groups, and free volume distributions, leading to interchain packing effects, are shown. The effect of the cyclohexyl ring at the isopropylidene carbon as compared to the effect of the methyl groups positioned on the phenylene rings results in a larger reduction in the solubility parameter (δ). The interchain distance estimated for polycarbonates in this work is in the range of 5-5.8 Å. The o-methyl groups on the phenylene rings, as compared to a cyclohexyl ring, lead to higher interchain distances. The highest interchain distance is observed with a trimethylcyclohexylidene group at the isopropylidene carbon. Atomistic simulations reveal two different types of packing arrangement of nearest-neighbor chains in the glassy state, one type of which agrees with the NMR experimental data. The fundamental insights provided here can be utilized for design of chemical structures for tailored macroscopic properties.
Collapse
Affiliation(s)
- M S Sulatha
- Macromolecular Modeling and Simulation Laboratory, Polymer Chemistry Division, National Chemical Laboratory (NCL), Pune 411 008, India
| | | |
Collapse
|
5
|
Fortunelli A, Geloni C, Lazzeri A. Simulation of the plastic behavior of amorphous glassy bis-phenol-A-polycarbonate. J Chem Phys 2004; 121:4941-50. [PMID: 15332930 DOI: 10.1063/1.1778157] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A protocol for studying the plastic deformation of amorphous glassy polymers is presented. The protocol is based on a viable computational procedure which combines constant-stress molecular dynamics simulations and fixed-cell energy minimizations, followed by kinetic, configurational, and energy analyses. It is shown that the computational results can be accounted for within a "potential energy landscape" theoretical framework, in which the plastic transition is interpreted as a crossing between and a collapse onto each other of "ideal (thermodynamic) structures." The procedure is applied to bis-phenol-A-polycarbonate (BPA-PC), but is equally valid for a wide variety of polymeric species. Allowing for the limited size of the simulation cell, the high strain rate, and the fact that the simulation are conducted at low temperature, the values of the density, Young's modulus, yield strain, yield stress, activation energy, and activation volume are in fair agreement with the experimental data on BPA-PC. The analysis of the results shows that the plastic relaxation for this polymer has both a collective and cooperative character (as in classical percolation theories), involves a significant fraction of the simulation cell, and can be viewed as a "nanoscopic shear band."
Collapse
Affiliation(s)
- Alessandro Fortunelli
- Molecular Modeling Laboratory, Istituto per i Processi Chimico-Fisici del C.N.R. (IPCF-CNR), via G. Moruzzi 1, 56124 Pisa, Italy.
| | | | | |
Collapse
|
6
|
Delle Site L, Leon S, Kremer K. BPA-PC on a Ni111 surface: the interplay between adsorption energy and conformational entropy for different chain-end modifications. J Am Chem Soc 2004; 126:2944-55. [PMID: 14995212 DOI: 10.1021/ja0387406] [Citation(s) in RCA: 53] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We extend a previous dual scale modeling approach for the behavior of polymers near a metal surface to a variety of end groups. Our approach combines a coarse-grained polymer model with ab initio DFT calculations. Such a procedure was applied to a melt of phenolic-like terminated Bisphenol A-polycarbonate (BPA-PC) interacting with a (111) nickel surface (Delle Site, L.; Abrams, C. F.; Alavi, A.; Kremer, K. Phys. Rev. Lett. 2002, 89, 156103. Abrams, C. F.; Delle Site, L.; Kremer, K. Phys. Rev. E 2003, 67, 021807). This work extends this study to different chain-end modifications of BPA-PC, p-tert-butylphenolic, p-tetramethylpropylphenolic, and p-cumylphenolic. We show how the interplay between adsorption energies and conformational entropy selects different morphologies for the various melts at the interface. Implications of these results for realistic technical materials are finally discussed.
Collapse
Affiliation(s)
- Luigi Delle Site
- Max-Planck-Institute for Polymer Research, Ackermannweg 10, D 55021 Mainz, Germany.
| | | | | |
Collapse
|
7
|
|
8
|
Molecular dynamics and thermal analysis study of anomalous thermodynamic behavior of poly (ether imide)/polycarbonate blends. POLYMER 2003. [DOI: 10.1016/s0032-3861(03)00054-5] [Citation(s) in RCA: 40] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
9
|
Abrams CF, Delle Site L, Kremer K. Dual-resolution coarse-grained simulation of the bisphenol-A-polycarbonate/nickel interface. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2003; 67:021807. [PMID: 12636708 DOI: 10.1103/physreve.67.021807] [Citation(s) in RCA: 33] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/19/2002] [Indexed: 05/24/2023]
Abstract
We present a dual-resolution coarse-graining scheme for efficient molecular dynamics simulations of bisphenol-A-polycarbonate (BP-A-PC) liquids in contact with a (111) nickel surface. The essential feature of this model is the strong adsorption of phenoxy chain ends, and the absence of adsorption by other parts of the chains. Details of how phenoxy chain ends interact with the nickel surface were extracted from Car-Parrinello molecular dynamics calculations of adsorption of phenol on nickel. These calculations show that phenol adsorption on nickel is short ranged (<3 A) and strongly dependent on the C1-C4 orientation of the ring. The structure of BP-A-PC prevents internal phenylene groups from interacting with the surface, due to steric hindrances from the noninteracting isopropylidenes. These dependencies are incorporated in the coarse-grained model of the BP-A-PC chain by resolving chain-terminating carbonate groups with atomistic detail, while the rest of the chain is represented by coarsened "beads." This allows specification of the C1-C4 orientation of the terminal phenoxy groups, while overall allowing for system equilibration with reasonable computer time. We simulate liquids of up to 240 chains of ten chemical repeat units, confined in a slit pore formed by two frozen (111) planes of atoms with the lattice spacing of nickel. We find that the strong adsorption of chain ends has a large effect on the liquid structure through a distance of more than two bulk radii of gyration from the surface. These effects are explained by a competition among single- and double-end adsorption, and dense packing. The structure of the interface less than 10 A from the wall is greatly sensitive to the orientational dependence of the phenoxy adsorption.
Collapse
Affiliation(s)
- Cameron F Abrams
- Max-Planck-Institut für Polymerforschung, Ackermannweg 10, 55128 Mainz, Germany.
| | | | | |
Collapse
|
10
|
Meyer H, Hahn O, Müller-Plathe F. Structure and Dynamics of Liquid Diphenyl Carbonate Investigated by Molecular Dynamics Simulations. J Phys Chem B 1999. [DOI: 10.1021/jp991733c] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Hendrik Meyer
- Max-Planck-Institut für Polymerforschung, 55021 Mainz, Germany
| | - Oliver Hahn
- Max-Planck-Institut für Polymerforschung, 55021 Mainz, Germany
| | | |
Collapse
|