Σχέσεις δομής και ιξωδοελαστικών, μηχανικών και συγκολλητικών ιδιοτήτων πολυακρυλικών σε στερεά υποστρώματα μέσω ατομιστικών προσομοιώσεων

Abstract

The present Doctoral Thesis focuses on the investigation, characterization and influence of polyacrylic materials in different scientific and technological disciplines via a detailed computer simulation using the Molecular Dynamics (MD) technique, in conjunction with the very accurate, all-atom Dreiding force-field. The main research concepts and objectives are discussed and analyzed in three separate parts. In the first part, atomistic configurations of two model pressure-sensitive acrylic adhesives (PSAs), the atactic homopolymer poly(n-BA) [poly(n-butyl acrylate)] and the atactic copolymer poly(n-BA-co-AA) [poly(n-butyl acrylate-co-acrylic acid)] in the bulk phase or confined between two selected substrates, glassy silica (SiO2) and metallic α-ferrite (α-Fe), were built and simulated by MD in the NPT statistical ensemble. First, an equilibration cycle consisting of temperature annealings and coolings was followed, in order to generate well-equilibrated configurations of the PSA systems. Detailed results from the atomistic simulations are presented concerning their volumetric behavior, glass transition temperature, conformational, structural, viscoelastic and dynamic properties. Particular emphasis was given to the analysis and characterization of the hydrogen bonds that form in the poly(n-BA-co-AA) system. By analyzing the MD trajectories, poly(n-BA-co-AA) was found to exhibit a higher density than poly(n-BA) by about 7% at all temperatures, to be characterized by smaller-size chains for a given molecular weight (MW), to exhibit significantly slower terminal and segmental dynamics properties, and to be characterized by a glass transition temperature that was approximately 40% higher than that of poly(n-BA). We also examined the type and degree of adsorption of the two acrylic systems on the selected substrates by analyzing the MD results for the local mass density as a function of distance from the solid plane and the distribution of adsorbed chain segments in train, loop, and tail conformations, and by computing the work of adhesion at the two substrates. The results revealed a stronger adsorption for both acrylics on the SiO2 surface due to highly attractive interactions between polymer molecules and substrate atoms, and as a consequence a higher value for the work of adhesion compared to that on the α-Fe surface. Furthermore, we have developed a generalized non-equilibrium molecular dynamics (NEMD) algorithm to simulate the mechanical response of the two adhesives under a uniaxial stretching deformation. In the second part of the Thesis, results have been obtained from a hierarchical simulation methodology that led to the prediction of the thermodynamic, conformational, structural, dynamic and mechanical properties of two polymer nanocomposites based on syndiotactic poly(methyl methacrylate) or sPMMA. The first was reinforced with uniformly dispersed graphene sheets and the second with fullerene particles. How graphene functionalization affects the elastic constants of the resulting nanocomposite has also been examined. The phase behavior of the nanocomposite (in particular as we varied the relative size between the sPMMA chains and the diameter of fullerene molecules) has also been studied as a function of fullerene volume fraction. The simulation strategy entailed three steps: 1) Generation of an initial structure, which was then subjected to potential energy minimization and detailed molecular dynamics (MD) simulations at T = 500K and P = 1atm to obtain well relaxed melt configurations of the nanocomposite. 2) Gradual cooling of selected configurations down to room temperature to obtain a good number of structures representative of the glassy phase of the polymer nanocomposite. 3) Molecular mechanics (MM) calculations of its mechanical properties following the method originally proposed by Theodorou and Suter. By analyzing the results under constant temperature and pressure, all nanocomposite systems were found to exhibit slower terminal and segmental relaxation dynamics than the pure polymer matrices. The addition of a small fraction of graphene sheets led in all cases to the enhancement of the elastic constants; this was significantly more pronounced in the case of functionalized graphene sheets. We further mention that, for all polymer/fullerene nanocomposites addressed here, no phase separation or variation of polymer chain dimensions was observed as a function of fullerene size and/or fullerene volume fraction.In the third part of the Thesis, and motivated by the use of acrylic polymers for the design of membranes with aligned carbon nanotubes (CNTs) for several separation technologies (such as water desalination and wastewater treatment), we report results from a detailed computer simulation study for the nano-sorption and mobility of four different small molecules (water, tyrosol, vanillic acid, and p-coumaric acid) inside smooth single-wall CNTs (SWCNTs). Most of the results have been obtained with the molecular dynamics (MD) method, but especially for the most narrow of the CNTs considered, the results for water molecule were further confirmed through an additional Grand Canonical (μVT) Monte Carlo (GCMC) simulation using a value for the water chemical potential μ pre-computed with the particle deletion method. Issues addressed in the Thesis include molecular packing and ordering inside the nanotube for the four molecules, average number of sorbed molecules per unit length of the tube, and mean residence time and effective axial diffusivities, all as a function of tube diameter and tube length. In all cases, a strong dependence of the results on carbon nanotube diameter was observed, especially in the way the different molecules are packed and organized inside the CNT. For water for which predictions of properties such as local structure and packing were computed with both methods (MD and GCMC), the two sets of results were found to be fully self-consistent for all types of SWCNTs considered. Water diffusivity inside the CNT (although, strongly dependent on the CNT diameter) was computed with two different methods, both of which gave identical results. For large enough CNT diameters (larger than about 13 Å), this was found to be higher than the corresponding experimental value in the bulk by about 55%. Surprisingly enough, for the rest of the (phenolic) molecules simulated in this Thesis, the simulations revealed no signs of mobility inside nanotubes with a diameter smaller than the (20, 20) tube. This has been attributed to strong phenyl-phenyl attractive interactions, also to favorable interactions of these molecules with the CNT walls, which cause them to form highly ordered, very stable structures inside the nanotube, especially under strong confinement. The interaction, in particular, of the methyl group (present in tyrosol, vanillic acid, and p-coumaric acid) with the CNT walls seems to play a key role in all these compounds causing them to remain practically immobile inside nanotubes characterized by diameters smaller than about 26 Å. It was only for larger-diameter CNTs that tyrosol, vanillic acid, and p-coumaric acid were observed to demonstrate appreciable mobility.