1
|
Siachouli P, Karadima KS, Mavrantzas VG, Pandis SN. The effect of functional groups on the glass transition temperature of atmospheric organic compounds: a molecular dynamics study. SOFT MATTER 2024; 20:4783-4794. [PMID: 38847330 DOI: 10.1039/d4sm00405a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/21/2024]
Abstract
Organic compounds constitute a substantial part of atmospheric particulate matter not only in terms of mass concentration but also in terms of distinct functional groups. The glass transition temperature provides an indirect way to investigate the phase state of the organic compounds, playing a crucial role in understanding their behavior and influence on aerosol processes. Molecular dynamics (MD) simulations were implemented here to predict the glass transition temperature (Tg) of atmospherically relevant organic compounds as well as the influence of their functional groups and length of their carbon chain. The cooling step used in the simulations was chosen to be neither too low (to supress crystallization) nor too high (to avoid Tg overprediction). According to the MD simulations, the predicted Tg is sensitive to the functional groups as follows: carboxylic acid (-COOH) > hydroxyl (-OH) and (-COOH) > carbonyls (-CO). Increasing the number of carbon atoms leads to higher Tg for the linearly structured compounds. Linear compounds with lower molecular weight were found to exhibit a lower Tg. No clear correlation between O : C and Tg was observed. The architecture of the carbon chain (linear, or branched, or ring) was also found to impact the glass transition temperature. Compounds containing a non-aromatic carbon ring are characterized by a higher Tg compared to linear and branched ones with the same number of carbon atoms.
Collapse
Affiliation(s)
- Panagiota Siachouli
- Department of Chemical Engineering, University of Patras, Patras, GR 26504, Greece.
- Institute of Chemical Engineering Sciences (ICE-HT/FORTH), Patras, GR 26504, Greece
| | - Katerina S Karadima
- Department of Chemical Engineering, University of Patras, Patras, GR 26504, Greece.
- Institute of Chemical Engineering Sciences (ICE-HT/FORTH), Patras, GR 26504, Greece
| | - Vlasis G Mavrantzas
- Department of Chemical Engineering, University of Patras, Patras, GR 26504, Greece.
- Institute of Chemical Engineering Sciences (ICE-HT/FORTH), Patras, GR 26504, Greece
- Particle Technology Laboratory, Department of Mechanical and Process Engineering, ETH Zürich, CH-8092 Zürich, Switzerland
| | - Spyros N Pandis
- Department of Chemical Engineering, University of Patras, Patras, GR 26504, Greece.
- Institute of Chemical Engineering Sciences (ICE-HT/FORTH), Patras, GR 26504, Greece
| |
Collapse
|
2
|
Shamsieva A, Piyanzina I, Minisini B. Amorphous cis-1,4-polybutadiene P-V-T properties from atomistic simulations. J Mol Model 2023; 29:249. [PMID: 37452231 DOI: 10.1007/s00894-023-05658-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2023] [Accepted: 07/08/2023] [Indexed: 07/18/2023]
Abstract
CONTEXT As a result of the diversity of microstructures encountered in cis-1,4-polybutadiene and the variety of measurement methods used, experimental values of variation of glass transition temperature (Tg) with pressure are relatively dispersed. However, atomistic simulations enable access to valuable information for very well-controlled chemistry and structures with a well-defined and systematic acquisition protocol. By varying the temperature and pressure, the specific volume of the melt was computed, yielding results that deviated by only 2% from experimental data. A linear relationship between Tg and pressure was observed, with Tg predicted to be 162 K at zero pressure and a rate of change of Tg with respect to pressure (dTg/dP) of 0.24 K/MPa. METHOD The atomistic dilatometry experiments were conducted on a model of amorphous cis-1,4 polybutadiene with an approximate molecular weight of 5400 g/mol using the LAMMPS code and the all-atom forcefield pcff + . The dilatometry process involved cooling and heating at a rate of 9 × 1012 K/min. The specific volume was calculated by averaging over seven independent configurations for each temperature. The Tait equation was employed to fit the specific volume evolution within the temperature range of 10 to 700 K under different pressures of 0, 60, and 100 MPa.
Collapse
Affiliation(s)
- Aigul Shamsieva
- Institute of Physics, Kazan Federal University, Kremlyovskaya St. 18, Kazan, 420008, Republic of Tatarstan, Russia
| | - Irina Piyanzina
- Institute of Physics, Kazan Federal University, Kremlyovskaya St. 18, Kazan, 420008, Republic of Tatarstan, Russia
| | - Benoit Minisini
- Materials Design SARL, 42 Avenue Verdier, 92120, Paris, Montrouge, France.
| |
Collapse
|
3
|
Phillies GDJ. Simulational Tests of the Rouse Model. Polymers (Basel) 2023; 15:2615. [PMID: 37376261 DOI: 10.3390/polym15122615] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2023] [Revised: 05/24/2023] [Accepted: 06/03/2023] [Indexed: 06/29/2023] Open
Abstract
An extensive review of literature simulations of quiescent polymer melts is given, considering results that test aspects of the Rouse model in the melt. We focus on Rouse model predictions for the mean-square amplitudes ⟨(Xp(0))2⟩ and time correlation functions ⟨Xp(0)Xp(t)⟩ of the Rouse mode Xp(t). The simulations conclusively demonstrate that the Rouse model is invalid in polymer melts. In particular, and contrary to the Rouse model, (i) mean-square Rouse mode amplitudes ⟨(Xp(0))2⟩ do not scale as sin-2(pπ/2N), N being the number of beads in the polymer. For small p (say, p≤3) ⟨(Xp(0))2⟩ scales with p as p-2; for larger p, it scales as p-3. (ii) Rouse mode time correlation functions ⟨Xp(t)Xp(0)⟩ do not decay with time as exponentials; they instead decay as stretched exponentials exp(-αtβ). β depends on p, typically with a minimum near N/2 or N/4. (iii) Polymer bead displacements are not described by independent Gaussian random processes. (iv) For p≠q, ⟨Xp(t)Xq(0)⟩ is sometimes non-zero. (v) The response of a polymer coil to a shear flow is a rotation, not the affine deformation predicted by Rouse. We also briefly consider the Kirkwood-Riseman polymer model.
Collapse
|
4
|
Rissanou A, Chazirakis A, Polinska P, Burkhart C, Doxastakis M, Harmandaris V. Polybutadiene Copolymers via Atomistic and Systematic Coarse-Grained Simulations. Macromolecules 2022. [DOI: 10.1021/acs.macromol.1c01939] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Anastassia Rissanou
- Institute of Applied and Computational Mathematics (IACM), Foundation for Research and Technology Hellas, (FORTH), IACM/FORTH, GR-71110 Heraklion, Greece
- Department of Mathematics and Applied Mathematics, University of Crete, GR-71409 Heraklion, Crete, Greece
| | - Antonis Chazirakis
- Institute of Applied and Computational Mathematics (IACM), Foundation for Research and Technology Hellas, (FORTH), IACM/FORTH, GR-71110 Heraklion, Greece
- Department of Mathematics and Applied Mathematics, University of Crete, GR-71409 Heraklion, Crete, Greece
| | | | - Craig Burkhart
- The Goodyear Tire and Rubber Company, 142 Goodyear Blvd., 44305 Akron, Ohio, United States
| | - Manolis Doxastakis
- Department of Chemical and Biomolecular Engineering, University of Tennessee, 37996 Knoxville, Tennessee, United States
| | - Vagelis Harmandaris
- Institute of Applied and Computational Mathematics (IACM), Foundation for Research and Technology Hellas, (FORTH), IACM/FORTH, GR-71110 Heraklion, Greece
- Department of Mathematics and Applied Mathematics, University of Crete, GR-71409 Heraklion, Crete, Greece
- Computation-Based Science and Technology Research Center, The Cyprus Institute, 2121 Nicosia, Cyprus
| |
Collapse
|
5
|
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
|
6
|
Tow GM, Maginn EJ. Fully Atomistic Molecular Dynamics Simulations of Hydroxyl-Terminated Polybutadiene with Insights into Hydroxyl Aggregation. Macromolecules 2020. [DOI: 10.1021/acs.macromol.9b02632] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Garrett M. Tow
- Department of Chemical and Biomolecular Engineering, University of Notre Dame, Notre Dame, Indiana 46556, United States
| | - Edward J. Maginn
- Department of Chemical and Biomolecular Engineering, University of Notre Dame, Notre Dame, Indiana 46556, United States
| |
Collapse
|
7
|
Kempfer K, Devémy J, Dequidt A, Couty M, Malfreyt P. Realistic Coarse-Grain Model of cis-1,4-Polybutadiene: From Chemistry to Rheology. Macromolecules 2019. [DOI: 10.1021/acs.macromol.8b02750] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- K. Kempfer
- CNRS, SIGMA Clermont, Institut de Chimie de Clermont-Ferrand, Université Clermont Auvergne, F-63000 Clermont-Ferrand, France
- Manufacture Française des Pneumatiques Michelin, 23, Place des Carmes, 63040 Clermont-Ferrand, France
| | - J. Devémy
- CNRS, SIGMA Clermont, Institut de Chimie de Clermont-Ferrand, Université Clermont Auvergne, F-63000 Clermont-Ferrand, France
| | - A. Dequidt
- CNRS, SIGMA Clermont, Institut de Chimie de Clermont-Ferrand, Université Clermont Auvergne, F-63000 Clermont-Ferrand, France
| | - M. Couty
- Manufacture Française des Pneumatiques Michelin, 23, Place des Carmes, 63040 Clermont-Ferrand, France
| | - P. Malfreyt
- CNRS, SIGMA Clermont, Institut de Chimie de Clermont-Ferrand, Université Clermont Auvergne, F-63000 Clermont-Ferrand, France
| |
Collapse
|
8
|
F Behbahani A, Vaez Allaei SM, H Motlagh G, Eslami H, Harmandaris VA. Structure and dynamics of stereo-regular poly(methyl-methacrylate) melts through atomistic molecular dynamics simulations. SOFT MATTER 2018; 14:1449-1464. [PMID: 29393331 DOI: 10.1039/c7sm02008b] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Poly(methyl-methacrylate), PMMA, is a disubstituted vinyl polymer whose properties depend significantly on its tacticity. Here we present a detailed study of the structure, conformation, and dynamics of syndiotactic, atactic, and isotactic PMMA melts at various temperatures (580, 550, 520, and 490 K) via all-atom molecular dynamics simulations. The calculated volumetric properties are close to experimental data. The Tg and chain dimensions of PMMA model systems are found to depend strongly on tacticity in agreement with experimental findings. The backbone bonds in trans (t), diads in tt, and inter-diads in t|t torsional states are the most populated for all PMMA stereo-chemistries and their fractions increase with the number of syndiotactic sequences. Also, the effective torsional barrier heights for the backbone, ester side group, and α-methyl group are larger for syndiotactic PMMA compared to the isotactic one. The structure of the PMMA chains is studied by computing the intra- and inter-chain static structure factors, S(q), and the radial pair distribution functions. In the first peak of S(q), both intra- and inter-chain components contribute, whereas the second and third peaks mainly come from inter- and intra-chain parts, respectively. For all PMMA stereo-isomers a clear tendency of ester-methyl groups to aggregate is observed. The local dynamics are studied by analyzing torsional autocorrelation functions for various dihedral angles. A wide spectrum of correlation times and different activation energies are observed for the motions of different parts of PMMA chains. The stereo-chemistry affects the backbone, ester side group, and α-methyl motions, whereas the ester-methyl rotation remains unaffected. The dynamic heterogeneity of the PMMA chains is also studied in detail for the different stereo-chemistries via the temperature dependence of the stretching exponent. Furthermore, the reorientational dynamics at the chain level and translational dynamics for monomer and chain centers-of-mass are analyzed.
Collapse
Affiliation(s)
- Alireza F Behbahani
- Advanced Polymer Materials and Processing Lab, School of Chemical Engineering, College of Engineering, University of Tehran, Tehran 11155-4563, Iran
| | | | | | | | | |
Collapse
|
9
|
Sharma P, Roy S, Karimi-Varzaneh HA. Validation of Force Fields of Rubber through Glass-Transition Temperature Calculation by Microsecond Atomic-Scale Molecular Dynamics Simulation. J Phys Chem B 2016; 120:1367-79. [DOI: 10.1021/acs.jpcb.5b10789] [Citation(s) in RCA: 44] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Pragati Sharma
- Physical
Chemistry Division, National Chemical Laboratory, Pune, India
| | - Sudip Roy
- Physical
Chemistry Division, National Chemical Laboratory, Pune, India
| | | |
Collapse
|
10
|
Lyubimov IY, Guenza MG. Theoretical reconstruction of realistic dynamics of highly coarse-grained cis-1,4-polybutadiene melts. J Chem Phys 2013; 138:12A546. [DOI: 10.1063/1.4792367] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
11
|
Gkourmpis T, Mitchell GR. Three Dimensional Picture of the Local Structure of 1,4-Polybutadiene from a Complete Atomistic Model and Neutron Scattering Data. Macromolecules 2011. [DOI: 10.1021/ma102840p] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Affiliation(s)
- Thomas Gkourmpis
- Polymer Science Centre, J. J. Thomson Physical Laboratory, Department of Physics, University of Reading, Reading, RG6 6AF, U.K
| | - Geoffrey R. Mitchell
- Polymer Science Centre, J. J. Thomson Physical Laboratory, Department of Physics, University of Reading, Reading, RG6 6AF, U.K
| |
Collapse
|
12
|
Harmandaris VA, Floudas G, Kremer K. Temperature and Pressure Dependence of Polystyrene Dynamics through Molecular Dynamics Simulations and Experiments. Macromolecules 2010. [DOI: 10.1021/ma102179b] [Citation(s) in RCA: 55] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Vagelis A. Harmandaris
- Department of Applied Mathematics, University of Crete, GR-71110 Heraklion, Greece
- IACM FORTH, GR-71110 Heraklion, Greece
- Max Planck Institute for Polymer Research, D-55128 Mainz, Germany
| | - George Floudas
- Department of Physics, University of Ioannina, GR-45110, Ioannina, Greece
- Foundation for Research and Technology-Biomedical Research Institute, Ioannina, Greece
| | - Kurt Kremer
- Max Planck Institute for Polymer Research, D-55128 Mainz, Germany
| |
Collapse
|
13
|
Turgut D, Atilgan AR, Atilgan C. Assortative mixing in close-packed spatial networks. PLoS One 2010; 5:e15551. [PMID: 21179578 PMCID: PMC3002975 DOI: 10.1371/journal.pone.0015551] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2010] [Accepted: 11/11/2010] [Indexed: 12/02/2022] Open
Abstract
Background In recent years, there is aroused interest in expressing complex systems as networks of interacting nodes. Using descriptors from graph theory, it has been possible to classify many diverse systems derived from social and physical sciences alike. In particular, folded proteins as examples of self-assembled complex molecules have also been investigated intensely using these tools. However, we need to develop additional measures to classify different systems, in order to dissect the underlying hierarchy. Methodology and Principal Findings In this study, a general analytical relation for the dependence of nearest neighbor degree correlations on degree is derived. Dependence of local clustering on degree is shown to be the sole determining factor of assortative versus disassortative mixing in networks. The characteristics of networks constructed from spatial atomic/molecular systems exemplified by self-organized residue networks built from folded protein structures and block copolymers, atomic clusters and well-compressed polymeric melts are studied. Distributions of statistical properties of the networks are presented. For these densely-packed systems, assortative mixing in the network construction is found to apply, and conditions are derived for a simple linear dependence. Conclusions Our analyses (i) reveal patterns that are common to close-packed clusters of atoms/molecules, (ii) identify the type of surface effects prominent in different close-packed systems, and (iii) associate fingerprints that may be used to classify networks with varying types of correlations.
Collapse
Affiliation(s)
- Deniz Turgut
- Faculty of Engineering and Natural Sciences, Sabanci University, Istanbul, Turkey
| | - Ali Rana Atilgan
- Faculty of Engineering and Natural Sciences, Sabanci University, Istanbul, Turkey
| | - Canan Atilgan
- Faculty of Engineering and Natural Sciences, Sabanci University, Istanbul, Turkey
- * E-mail:
| |
Collapse
|
14
|
Tsolou G, Stratikis N, Baig C, Stephanou PS, Mavrantzas VG. Melt Structure and Dynamics of Unentangled Polyethylene Rings: Rouse Theory, Atomistic Molecular Dynamics Simulation, and Comparison with the Linear Analogues. Macromolecules 2010. [DOI: 10.1021/ma1017555] [Citation(s) in RCA: 101] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Georgia Tsolou
- Department of Chemical Engineering, University of Patras & FORTH-ICE/HT, Patras, GR 26504, Greece
| | - Nikos Stratikis
- Department of Chemical Engineering, University of Patras & FORTH-ICE/HT, Patras, GR 26504, Greece
| | - Chunggi Baig
- Department of Chemical Engineering, University of Patras & FORTH-ICE/HT, Patras, GR 26504, Greece
| | - Pavlos S. Stephanou
- Department of Chemical Engineering, University of Patras & FORTH-ICE/HT, Patras, GR 26504, Greece
| | - Vlasis G. Mavrantzas
- Department of Chemical Engineering, University of Patras & FORTH-ICE/HT, Patras, GR 26504, Greece
| |
Collapse
|
15
|
Pawlus S, Sokolov AP, Paluch M, Mierzwa M. Influence of Pressure on Chain and Segmental Dynamics in Polyisoprene. Macromolecules 2010. [DOI: 10.1021/ma100383r] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Sebastian Pawlus
- Institute of Physics, University of Silesia, Uniwersytecka 4, 40-007 Katowice, Poland
| | - Alexei P. Sokolov
- Chemical Sciences Division, ORNL, and Department of Chemistry, University of Tennessee, Knoxville, Tennessee 37996-1600
| | - Marian Paluch
- Institute of Physics, University of Silesia, Uniwersytecka 4, 40-007 Katowice, Poland
| | - Michal Mierzwa
- Institute of Physics, University of Silesia, Uniwersytecka 4, 40-007 Katowice, Poland
| |
Collapse
|
16
|
Liu J, Wu S, Cao D, Zhang L. Effects of pressure on structure and dynamics of model elastomers: A molecular dynamics study. J Chem Phys 2008; 129:154905. [DOI: 10.1063/1.2996009] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
17
|
Tsolou G, Mavrantzas VG, Makrodimitri ZA, Economou IG, Gani R. Atomistic Simulation of the Sorption of Small Gas Molecules in Polyisobutylene. Macromolecules 2008. [DOI: 10.1021/ma8007652] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Affiliation(s)
| | | | | | | | - Rafiqul Gani
- CAPEC, Department of Chemical and Biochemical Engineering, Technical University of Denmark, DK-2800 Lyngby, Denmark
| |
Collapse
|
18
|
Mpoukouvalas K, Floudas G. Effect of Pressure on the Dynamic Heterogeneity in Miscible Blends of Poly(methyl methacrylate) with Poly(ethylene oxide). Macromolecules 2008. [DOI: 10.1021/ma702299n] [Citation(s) in RCA: 29] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- K. Mpoukouvalas
- Department of Physics, University of Ioannina, P.O. Box 1186, 451 10 Ioannina, Greece, and Foundation for Research and Technology-Hellas (FORTH), Biomedical Research Institute (BRI)
| | - G. Floudas
- Department of Physics, University of Ioannina, P.O. Box 1186, 451 10 Ioannina, Greece, and Foundation for Research and Technology-Hellas (FORTH), Biomedical Research Institute (BRI)
| |
Collapse
|