Mechanical Properties of a Polymer at the Interface Structurally

Hess , C.; Kutzner; van der Spoel , D.; Lindahl , E. J. GROMACS 4: Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation J...
1 downloads 0 Views 1MB Size
Subscriber access provided by GAZI UNIV

Article

Mechanical Properties of a Polymer at the Interface Structurally Ordered by Graphene Stanislav G. Falkovich, Victor M. Nazarychev, Sergey V. Larin, Josè M. Kenny, and Sergey V. Lyulin J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.5b11028 • Publication Date (Web): 04 Mar 2016 Downloaded from http://pubs.acs.org on March 7, 2016

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry C is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Mechanical properties of a polymer at the interface structurally ordered by graphene Stanislav G. Falkovich,a Victor M. Nazarychev,a Sergey V. Larin,a José M. Kenny,a,b Sergey V. Lyulin a,c,* a

Institute of Macromolecular Compounds, Russian Academy of Sciences, Bol’shoi pr. 31 (V.O.), St. Petersburg, 199004 Russia.

b

Materials Science and Technology Centre, University of Perugia, Loc. Pentima, 4, Terni, 05100 Italy.

c

Department of Physics, St. Petersburg State University, Ul’yanovskaya str. 1, Petrodvorets, St. Petersburg, 198504 Russia.

*Corresponding author. Sergey Lyulin Fax: +7 (812) 328686; Tel: +7 (812) 3285601; E-mail: [email protected]

ACS Paragon Plus Environment

1

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 32

KEYWORDS Nanocomposite, Graphene, Molecular dynamics, Computer Simulation, Mechanical properties, Polyimides

ABSTRACT

In this study, for the first time, molecular-dynamics simulations helped identifying the substantial mechanical properties changes of a crystallizable polyimide related to the ordering of the polymer chains along the graphene surface. It is demonstrated that graphene-induced polymer ordering at the nanofiller-polymer interface leads to mechanical properties exceeding those of the unfilled amorphous polymer. Moreover, simulation results confirm that the mechanical properties anisotropy is related to the polymer ordering anisotropy induced by graphene-polymer interactions.

INTRODUCTION Polymer composites are widely used in many industrial applications. Recently, researchers have increased the focus on nanocomposites since the addition of low contents of nano-scale filler particles may significantly improve the thermal and mechanical properties of the polymer matrix.1-5 Due to the high specific area of dispersed nanoparticles, nanocomposites feature a higher volume fraction of the polymer at the interface than conventional composites. The properties of the nanoparticle-polymer interface, including the mechanical ones, depend both on

ACS Paragon Plus Environment

2

Page 3 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

the properties of the nanoparticles themselves and of the polymeric matrix in the near-surface area (so-called interface). It is almost impossible to experimentally distinguish contributions between nanoparticles themselves and that of the polymer at interface to the change of the nanocomposite mechanical properties. The question that arises is how important is the contribution of the structurally ordered polymer at the interface? This question can be worded differently: do the properties of the polymer at interface differ significantly from the properties of the polymer without filler? If yes, how big is this difference and how thick is the interface? Answers to these questions are relevant for the development of new polymer nanocomposites with controlled properties. They are also important for solving engineering tasks during the simulation using the finite elements approach, which may account for specific properties of the polymer at interface.6 The investigation reported here deals with structural and mechanical properties of the polymer at the interface using the atomistic molecular-dynamics simulation which considers the polymer chemical structure in detail. In order to answer these questions the heat-resistant crystallizable polyimide R-BAPB (based on 1,3-bis-(3’,4-dicarboxyphenoxy)benzene (dianhydride R) and diamine 4,4’-bis-(4’’-aminophenoxy)-diphenyl (diamine BAPB)) has been taken in consideration (Fig. 1). The objective of this study is to clear up the influence of polymer precrystallization ordering induced by the interaction with the graphene nanoparticle on the mechanical properties of nanocomposites. There is a number of studies devoted to semi-crystalline polymer nanocomposites with carbon nanofillers and their mechanical properties, both experimental and simulation. Most of the atomistic modeling of such systems is devoted to nanocomposites based on polymers with relatively simple chemical structure, such as alkanes. It was shown in the study by Wang et al.7, by atomistic molecular dynamics simulations, that the crystallization of a polyethylene single

ACS Paragon Plus Environment

3

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 32

chain on the surface of graphene occurs into two steps: the adsorption of the chains and further crystallization. A similar result was obtained in the Monte Carlo simulation in Nie et al. 8 and in Anwar et al.9,10 The ordering of one or several alkane chains on the surface of nanotubes or graphene was studied by Yang et al.11and Anwar et al.9,10by atomistic molecular dynamics method. It was shown that the chain of alkanes immediately pulled along a surface of the CNT (which is a kind of director for their ordering), while in the interaction with graphene alkane chains first adsorb on graphene flat surface and only after that ordering starts. Therefore, crystallization induced by nanotubes is faster than the crystallization induced by graphene. A similar result was obtained by Xu et al.12 in the experimental study of the crystallization of polylactic acid, which occurs in the presence of CNTs faster than in the presence of graphene. It should also be noted that there are experimental data for semi-crystalline polyimide R-BAPB,13 showing that this polymer crystallizes faster in the presence of graphitic carbon than in the presence of CNTs possibly due to an orientating effect of the graphene like surface. A potential orienting effect of graphene has been indirectly confirmed by the results of Wei, who showed by atomistic molecular dynamics that, depending on the chirality of CNT, alkanes crystallization occurs with varying degrees of effectiveness. Furthermore, it was shown in works of Rissanou et al.14and Bacova et al.15 by atomistic molecular dynamics simulations that a mobility of polyethylene near the surface of the graphene is slowed down. The authors also showed anisotropy of molecular diffusion, which was lower along the direction perpendicular to the nanotube axis than along a direction parallel to it. This observation contradicts the results obtained by atomistic molecular dynamics simulations for polyimides16,17 where a stretching of the polymer chains along the axis of the nanotube or graphene surface when ordering was shown.

ACS Paragon Plus Environment

4

Page 5 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

The crystallization of polymers in the presence of carbon nanofillers has been reported in several papers: in addition to the above-mentioned work by Yudin et al.,13 we would like to refer to the experimental papers of Jia et al. and of the group of Dingemans,18 where the crystallization of a non-crystallizable polyimide was shown in the presence of CNTs. Furthermore, they reported that a crystallization resulted in a significant increase in the Young's modulus. In particular, the results reported by Hegde et al.18 emphasize the fundamental role of the CNTpolymer interface in providing higher values of the nanocomposites Young's modulus. An increase in the modulus due to crystallization caused by nanotubes and graphene (see Sharma et al., 201419and Sharma et al., 201520, respectively) was also observed for nanocomposites of poly vinyl alcohol. Changes in the properties of the polymer interacting with the carbon nanofillers were studied in Liu et al. using molecular dynamics simulations. A more significant increase in stiffness of the polymer was shown in the case of graphene, than in the case of the nanotube or fullerene, which is attributed to the formation of long fibrils of oriented polymer, in agreement with our previous study of an effect of carbon nanofiller curvature on crystallization of polyimides.16 The mechanical properties of graphene-epoxy nanocomposites were studied in the paper of Rahman by atomistic molecular dynamics simulations. A substantial increase in the Young's modulus upon addition of graphene to an epoxy matrix was shown. It was also demonstrated that a well-dispersed graphene possesses enhanced elastic modulus compared to the agglomerated graphene-epoxy system. The anisotropy of the polymer caused by crystallization and influencing its mechanical properties was studied in a number of experimental researches. The appearance of ordered structures of stacked parallel polymer chains oriented along CNTs (shish kebab) was shown for a number of polymers such as polyethylene21 and polylactic acid.22 Such ordering in polymer

ACS Paragon Plus Environment

5

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 32

systems can appear not only related to the presence of a carbon nanofiller but as a consequence of the tension applied as it was shown for polyethylene,23 PET24 and other polymers. A similar effect is strain (strain-induced crystallization) as it was shown for poly(keto-ether-imide) films, which is similar in chemical structure to polyimides.25 After a careful revision of the scientific literature, it is evident that the process of crystallization of polymers initiated by carbon nanofillers, as well as the study of the mechanical properties of nanocomposites with carbon fillers have been widely investigated both experimentally and theoretically using atomistic simulations. However, it should be noted that in all the above-mentioned papers on atomistic simulation in which the mechanical properties of polymer nanocomposites were investigated, the load was applied to the sample, which included both the polymer matrix and the nanoparticles. In our atomistic simulation we have separated the components after the interaction and we have calculated the mechanical properties of the ordered polymer without the nanofiller. In fact, we are convinced that the computation of the intrinsic mechanical properties of the polymer, ordered by the filler nanoparticles has a relevant interest. In particular, modeling such systems is important for a correct parameterization of calculations by finite elements. Obviously, such a study of the mechanical properties of nanocomposites separating the mechanical properties of the polymer matrix from those of a nanofiller is not possible experimentally but we demonstrate here that this problem can be solved via atomistic simulation.

ACS Paragon Plus Environment

6

Page 7 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 1. Chemical structure of the repeated unit of R-BAPB polyimide (upper part) and part of the infinite graphene sheet16 (lower part). The red and blue arrows indicate the vectors corresponding to the rigid fragment of the R-BAPB repeated unit chosen for structural analysis.

The choice of the R-BAPB-based graphene-filled nanocomposite as research target is determined by the following reasons: R-BAPB is a semi-crystalline thermoplastic polyimide and the incorporation of nanoparticles into crystallizable polymers further improves their mechanical properties26 increasing the material rigidity and preventing microscopic crack formation at the interface between crystalline and amorphous areas taking advantage of the formation of multiple randomly oriented nanoscale crystallites.18,27-29 The excellent thermal and mechanical properties of semi-crystalline thermoplastic polymers determine that these polymers and their nanocomposites are highly promising materials to be used in aerospace, automotive and other industries.4,18,28 For R-BAPB, it was shown experimentally that the addition of carbon nanofillers may significantly accelerate its crystallization rate compared to the rate of the non-filled polymer.13 To explain this phenomenon, in our previous studies16,17 we simulated the interaction of carbon nanofillers (graphene and nanotubes) and the R-BAPB matrix using atomisticallydetailed molecular dynamics. Notably, such simulations allow to study the initial steps of the crystallization of polymers such as polyolefines;9-11 obtaining results in line with experimental

ACS Paragon Plus Environment

7

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 32

data.30 In the case of less flexible and more complex heterocyclic polymers with a higher viscosity, the simulation of a crystallization process is hindered even when high-performance computational packages such as Gromacs,31,32 and modern supercomputers are used. However, during the microsecond-scale simulation we have observed a substantial number of linear fragments of repeated units of this polyimide oriented parallel to the nanofiller surface,16,17 which can be interpreted as a precursor stage of the crystallization.9,10 Graphene causes stronger ordering of R-BAPB than CNT; a greater number of chain segments are oriented parallel to the graphene plane than those which are oriented parallel to the CNT surface, which correlates well with experimental data.13 The presence of graphene affected the R-BAPB in all the volume of the sample considered. In fact, ordering involved not only the first sub-surface layer of R-BAPB but the whole volume of the polymer. This paper deals with the investigation of the influence of the structural ordering of the polymer at the interface with graphene on the mechanical properties of R-BAPB.

THEORETICAL METHODS A repeat unit of R-BAPB consists of dianhydride R (1,3-bis-(3',4-dicarboxyphenoxy)-benzene) and diamine BAPB (4,4'-bis-(4''-aminophenoxy)-diphenyl. Our study is based on our experience in atomistic simulation of polyimides,16,17,33-38 including R-BAPB. For the simulation, we used R-BAPB chains with a polymerization degree of Np = 8 in the samples of graphene-filled and non-filled R-BAPB samples containing 27 chains each. The atomistic simulation using the molecular-dynamics method was performed in the force field Gromos53a539 using the Gromacs package.31,32 The samples were obtained from the rarefied chain gas; the procedure is described in detail in our previous papers.16,17,33,37 The compaction was followed by a microsecond-long

ACS Paragon Plus Environment

8

Page 9 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

simulation. During the simulation, the average chain size reached constant values close to those calculated analytically.16,17,33 The simulations were carried out in the NpT ensemble. Constant temperature T = 600 K, that is 5K higher than the experimentally-determined value of the R-BAPB melting temperature,13,40 and pressure P = 1 bar were maintained using a Berendsen thermostat and barostat,41,42 with time constants of τT = 0.1 ps and τP = 0.5 ps, respectively. All bond lengths were kept constant with the LINCS algorithm.43 Atomistic partial charges were taken equal to zero since their negligence has a weak impact on the structural properties of the relatively weakly polarized R-BAPB, accelerating the polymer structural transformations by many times.16,17 During the simulation of the system with graphene, the mobility of its atoms was restricted by introducing additional potentials (position restraints) and a change of the periodic cell size in XY plane of graphene was prohibited in order to maintain the graphene sheet area.7 Nine configurations of the graphene-filled R-BAPB samples were chosen during productive run after each 50 ns and 11 configurations of the non-filled amorphous R-BAPB were chosen after each 100 ns. Starting from these configurations, the systems were cooled from 600 to 290 К which is much lower than the glass-transition point of R-BAPB.13 The research procedure was as follows: the polymer nanocomposite was equilibrated for 1 microsecond at T = 600 K and P = 1 bar and waited for the polymer matrix structural ordering near the graphene to occur.16,17 Then the system was cooled down to T = 290 K well below the glass transition temperature and hence all polymer degrees of freedom were arrested. Nine configurations of the system were selected each 100 ns from the final part of the MD trajectory to be cooled and then pulled. As it was shown in our previous papers,17,33,37 100 ns in enough to consider configurations of a nanoscale R-BAPB system as independent. Thus, a set of

ACS Paragon Plus Environment

9

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 32

independent quenched samples was obtained for further study of their mechanical properties. The cooling procedure consisted in gradual decrease of the temperature, with the temperature being reduced by 10 K at each step, followed by a 4 ns relaxation of the system.17,33,37 The size of the simulation cell was 6×6×6 nm. Then graphene atoms were removed simultaneously and the gap produced was eliminated by a short NpT run (1.5 ns at 50 bar along Z direction perpendicular to the plane where graphene laids, and 20 ns of equilibration at 1 bar), i.e. the two ordered surfaces were basically "glue" or pressed against each other to create a "bulk"-like environment, yet preserving the chain ordering created by the graphene. This was followed by the simulation of the uniaxial deforming of the structurally ordered polymer at the interface in various directions. For reference, mechanical properties of the non-filled amorphous R-BAPB sample were studied as well. Periodic boundary conditions were used at all simulation stages. The uniaxial deformation changes the periodic cell size at a constant rate, along the positive direction of one – X, Y or Z - of the reference axes, so that the isotropic Berendsen barostat with a time constant  =0.5 ps was replaced with the anisotropic Berendsen barostat with  =1 ps. The samples remained incompressible in the direction of the applied deformation, i.e. the compressibility of the system in this direction was set to zero. In the transverse direction the system compressibility was set to 4.5×10-10 Pa-1. Therefore, upon stretching, a simulation cell elongates in the direction of deformation and compresses in the directions transversal to deformation in response to the external pressure (1 bar). During deformation the values  ,  = , , of the pressure tensor and the simulation cell size  in the stretching direction were saved each 1 ps. The characteristics obtained were converted to the dependence of the stress σ on the relative strain  as

ACS Paragon Plus Environment

10

Page 11 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

 = −  = (  −  )/ 

(1)

where  is the simulation cell size prior to the deformation (t=0). The initial part of the () dependence shows a linear regime up to ~2% of the deformation , and the elastic modulus E is defined as  = 

(2)

The value of E was determined as a slope of the linearly approximated dependence σ(ε) in the linear viscoelasticity regime. The error bars in elasticity modulus calculations were computed as mean-square deviations from the average value of the elastic modulus, obtained by averaging over all samples and the three deformation directions. To study the mechanical properties of the amorphous and graphene-ordered R-BAPB samples, we carried out a uniaxial deformation of each of them with a rate of 1 m/s along the axes X, Y and Z, with the LINCS algorithm that maintains the covalent bond length constant switched off. The stress-strain curves obtained for each sample served as the basis for the calculation of the elasticity modulus and yield stress values. To determine the elasticity modulus, the initial linear section of the stress-strain curve with the relative strain of 0-2% was approximated by a linear function, and the elasticity modulus value was calculated as a tangent of the approximating straight line. The yield point was determined from the stress-strain curve as the offset yield peak, i.e. from the intersection height of the axis of ordinates and the straight line approximating the plateau of stress-strain curves.44The elasticity modulus and yield point values obtained for individual samples of the graphene-ordered and amorphous R-BAPB samples were approximated, with errors being evaluated as mean square deviations from the average value.

ACS Paragon Plus Environment

11

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 32

The presence of trustworthy distinctions between the average values was assessed using the Mann-Whitney criterion.45 In order to check out whether the removal of the graphene influences the properties of an ordered polymer matrix we have tested the MD-trajectory and found that the “glued” gap at the place where graphene was located does not become wider under pulling, see figure 2, where the gap is located at the left and right sides of density profiles. It is seen that polymer density at the gap increases under pulling. Thus we can conclude that the gap is not the “weak site” of the system and a comparably low resistance of the ordered polymer matrix to pulling in the Z direction is really due to its structural peculiarities, not due to the gap.

Figure 2. Density profiles along Z axis for R-BAPB with removed graphene before pulling (black curve) and at 10% strain (red curve).

RESULTS AND DISCUSSION Figure 3 shows the stress-strain curves for graphene-ordered and for amorphous R-BAPB polyimide samples averaged upon deforming along the axes X, Y and Z. According to Figure 3 results, the stress-strain curves for both the oriented and amorphous samples have a similar shape: initial growth until a the plateau is reached without any recognizable yield point. At the

ACS Paragon Plus Environment

12

Page 13 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

same time, the stress-strain curve for the oriented sample is located above the one for the amorphous sample along the entire strain length. Table 1 reports the results of the comparison of the elasticity modulus 〈〉 and the yield stress 〈 〉 values for amorphous and graphene-ordered R-BAPB samples, with errors bars being evaluated as mean square deviations from the average values obtained during the processing of individual stress-strain curves.

Figure 3. Stress-strain curves averaged over three directions (along the axes X, Y and Z) both for R-BAPB samples ordered by graphene and for amorphous R-BAPB samples.

Table 1. Comparison of mechanical properties (average values of elasticity modulus 〈〉 and yield stress values 〈 〉) of amorphous and graphene-ordered R-BAPB samples. 〈〉, GPa

〈 〉, MPa

amorphous

2.7±0.3

169±4

ordered by graphene

3.8±0.9

191±4

As shown in Table 1, the elasticity modulus and yield stress values for the graphene-ordered R-BAPB samples are higher than those for the amorphous one. Experiments also demonstrate3,4

ACS Paragon Plus Environment

13

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 32

that the addition of graphene can improve mechanical properties of the material. However, these experiments addressed the mechanical properties of nanocomposites filled with graphene whereas we simulate mechanical properties of the ordered polymer at the polymer-nanofiller interface showing that the polymer at the interface can significantly contribute to the mechanical properties of the nanocomposite. Therefore, such properties are clearly influenced by factors other than mechanical properties of the filler nanoparticles. Consequently, the specific behavior of the polymer at the interface should be taken into account when solving engineering tasks, for example, while using the finite elements simulations for nanocomposites.

Figure 4. Averaged stress-strain curves of (a) non-filled amorphous R-BAPB samples and (b) graphene-ordered R-BAPB samples, obtained after deforming along axes X (black curves), Y (red curves), and Z (blue curves)

In order to take a closer look at the influence of the R-BAPB structural ordering caused by graphene on its mechanical properties, we considered separately the stress-strain curves obtained as a result of deforming along the axes X, Y and Z (the graphene sheet was located in the XY plane). Figure 4 shows the curves averaged over all individual stress-strain curves, and

ACS Paragon Plus Environment

14

Page 15 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

demonstrates that all three stress-strain curves have the same shape as the averaged curve in Figure 3. Moreover, the stress-strain curve obtained upon deforming along the X axis is located higher than the curve obtained upon deforming along the Y axis which, in turn, is located above the curve obtained upon deforming along the Z axis. At the same time, in the case of the nonordered (amorphous) R-BAPB sample, the stress-strain curves obtained after deforming in various directions coincide. This fact confirms the significant anisotropy of mechanical properties of the polymer at the interface. Such shape of the stress-strain curves is an argument supporting the fact that the low resistance of the ordered polymer matrix to pulling in the Z direction is not caused by the gap produced by the graphene removal. In fact, in the case of a gap growth under pulling the stress should decrease, but this is not observed, see figure 4b. It is seen from Fig. 4 that the stress-strain curve for the amorphous polymer lays between the curves for Y and Z directions for the ordered polymer. This may be attributed to the amorphous polymer “avoiding” deformation paths where the already stretched segments should be additionally stretched as in the case of pulling along X direction. Table 2 features the elasticity modulus and yield stress values obtained during deforming along various axes averaged for 9 samples. The comparison of these results also confirms the anisotropy of the polymer mechanical properties at the interface: elasticity modulus values obtained by deforming along various axes differ (Mann-Whitney significance level45 p < 0.05), as well as the yield stress values (p < 0.01) show a similar behavior. On the contrary, the amorphous R-BAPB sample displays no anisotropy of mechanical properties.

ACS Paragon Plus Environment

15

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 32

Table 2. Comparison of mechanical properties of the graphene-ordered R-BAPB sample – average values of elasticity modulus 〈〉 and yield stress 〈 〉 obtained at deforming in various directions. amorphous R-BAPB

ordered R-BAPB

deforming  

!"

〈〉, GPa

〈 〉, MPa

〈〉, GPa

〈 〉, MPa

X

2.8±0.3

165±3

4.8±0.8

244±11

0.27±0.09

0.22±0.09

Y

2.7±0.2

169±4

3.8±0.9

190±6

0.17±0.08

0.20±0.10

Z

2.8±0.4

171±3

2.9±0.4

140±6

0.02±0.01

0.03±0.03

direction

To answer the question why the structural ordering causes the anisotropy of mechanical properties of R-BAPB, we investigated the relation between the computed mechanical properties and the orientation degree of rigid segments of R-BAPB chains (Fig. 1) along the coordinate axes. The segments are regarded as oriented along one of the axis if the angle between them and a given axis is less than a critical angle assumed as equal to 30 degrees (results for other critical angle values were qualitatively similar). We calculated the average proportion of segments  , oriented along each of the axes X, Y and Z, for two representative polymer segments and the results are given in Table 2. We can see that the fraction of segments oriented along the axis X is always higher than that of segments oriented along the axis Y, and that the fraction of segments oriented along the axis Y is much higher than that of segments oriented along the axis Z (p < # % & 0.01); that is  >  >  . Therefore, it is possible to establish a qualitative correlation

where most molecules are parallel to the XY plane and a stronger orientation of rigid fragments leads to a much stronger resistance of the sample in the same direction. Theses finding confirm

ACS Paragon Plus Environment

16

Page 17 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

that the structural ordering of R-BAPB fragments is responsible for the anisotropy of the mechanical properties of the polymer ordered at the matrix-graphene interface. Moreover, we can confirm that the origin of the polymer anisotropy at the matrix-nanoparticle interface is provided by the nanostructuration driven by the graphene structure. It is very well known that π-π interactions promote strong interactions among aromatic structures46 like those present in R-BAPB chains and in the graphene structure. Moreover, the π-π interactions of graphene with aromatic molecules have been well reported.47 A difference in tensile properties in X and Y is seen from Table 2 and Figure 4. It is unexpected because R-BAPB is ordered perpendicular to Z direction and thus a mechanical isotropy in X and Y is expected. One of the possible causes of this difference is an orientation effect of graphene which leads to an orientation of R-BAPB segments along X axis. Another possible cause is that an appearance of 3D ordering characteristic to crystals may occur even at the pre-crystallization stage leading to an anisotropy of mechanical properties along all three axis. This effect will be studied in our future simulations.

CONCLUSIONS The mechanical properties of a thermoplastic semi-crystalline polyimide R-BAPB, ordered by graphene were studied by atomistic molecular dynamics simulations. A novel method of modeling was developed, which consists in removing the graphene sheet before deformation, allowing the study of the intrinsic properties of the ordered polymer. The obtained results confirm that a pre-crystallization, also showed for polyolefins,9,10 affects the mechanical properties of the polyimide. This is crucial for the correct parameterization of the interface properties, which is necessary for simulation of the microscopic characteristics of

ACS Paragon Plus Environment

17

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 32

nanocomposites by finite elements methods48. A comparison of the results of the simulation of polymer nanocomposites by atomistic molecular dynamics with experimental data is always a non-trivial task. One of the few ways to compare is to study the dielectric relaxation by means of both simulation and experiments (see, eg, Wu et al.49 and Klonos et al.50); however, this study is devoted to a research of another kind. The results obtained in the study reported here can be compared with experimental studies of the mechanical properties of nanocomposites based on crystallizing polyimides filled with ordered graphene sheets. Such ordering can be achieved with the help of the application of a strong electric field.46 The method of studying the nanocomposite films using bugle test seems to be promising for this kind of comparison.51,52 The results can also be compared with those obtained in computer simulations of polymer nanocomposites filled with small graphene sheets.52 Therefore, the anisotropy of mechanical properties is relevant for practical applications of nanocomposites based on crystallizable polymers and graphene for micro- and nanoscale devices.53

ACKNOWLEDGMENT The study was carried out with the financial support from the Ministry of Education and Science of the Russian Federation under the Contract No. 14.Z50.31.0002 (megagrant of the Government of the Russian Federation according to the Resolution No. 220 of April 9, 2010). The simulations were carried out with the use of the computational resources of the Institute of Macromolecular Compounds, Russian Academy of Sciences, and the Lomonosov supercomputer at Moscow State University.

ACS Paragon Plus Environment

18

Page 19 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

REFERENCES

(1)

Stankovich, S.; Dikin, D. A.; Dommett, G. H.; Kohlhaas, K. M.; Zimney, E. J.; Stach, E. A.; Piner, R. D.; Nguyen, S. T.; Ruoff, R. S. Graphene-Based Composite Materials. Nature 2006, 442, 282-286.

(2)

Tang, H.; Ehlert, G. J.; Lin, Y.; Sodano, H. A. Highly Efficient Synthesis of Graphene Nanocomposites. Nano Lett. 2012, 12, 84-90.

(3)

Singh, V.; Joung, D.; Zhai, L.; Das, S.; Khondaker, S. I.; Seal, S. Graphene Based Materials: Past, Present and Future. Prog. Mater Sci. 2011, 56, 1178-1271.

(4)

Yoonessi, M.; Shi, Y.; Scheiman, D. A.; Lebron-Colon, M.; Tigelaar, D. M.; Weiss, R. A.; Meador, M. A. Graphene Polyimide Nanocomposites; Thermal, Mechanical, and High-Temperature Shape Memory Effects. ACS Nano 2012, 6, 7644-7655.

(5)

Zhang, T.; Xue, Q.; Zhang, S.; Dong, M. Theoretical Approaches to Graphene and Graphene-Based Materials. Nano Today 2012, 7, 180-200.

(6)

Zeng, Q. H.; Yua, A. B.; Lu, G. Q. Multiscale Modeling and Simulation of Polymer Nanocomposites. Prog. Polym. Sci. 2008, 33, 191-269.

ACS Paragon Plus Environment

19

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(7)

Page 20 of 32

Wang, L. Z; Duan, L. L. Isothermal crystallization of a single polyethylene chain induced by graphene: A molecular dynamics simulation. Comput. Theor. Chem. 2012, 1002, 5963.

(8)

Nie, Y. J; Zhang, R. J; Zheng, K. S; Zhou, Z. P. Nucleation details of nanohybrid shishkebabs in polymer solutions studied by molecular simulations. Polymer, 2015, 76, 1-7.

(9)

Anwar, M.; Turci, F.; Schilling, T. Crystallization Mechanism in Melts of Short n-alkane Chains. J. Chem. Phys. 2013, 139, 214904.

(10) Anwar, M.; Berryman, J. T.; Schilling, T. Crystal Nucleation Mechanism in Melts of Short Polymer Chains under Quiescent Conditions and under Shear Flow. J. Chem. Phys. 2014, 141, 124910. (11) Yang, J. S.; Yang, C. L.; Wang, M. S.; Chen, B. D.; Ma, X. G. Crystallization of Alkane Melts Induced by Carbon Nanotubes and Graphene Nanosheets: a Molecular Dynamics Simulation Study. Phys. Chem. Chem. Phys. 2011, 13, 15476-15482. (12) Xu, J. Z.; Chen, T.; Yang, C. L.; Li, Z. M.; Mao, Y. M.; Zeng, B. Q.; Hsiao, B. S. Isothermal Crystallization of Poly(L-lactide) Induced by Graphene Nanosheets and Carbon Nanotubes: A Comparative Study. Macromolecules 2010, 43(11), 5000-5008. (13) Yudin, V. E.; Feldman, A. Y.; Svetlichnyi, V. M.; Shumakov A. N.; Marom, G. Crystallization of R-BAPB Type Polyimide Modified by Carbon Nanoparticles. Compos. Sci. Technol. 2007, 67, 789-794.

ACS Paragon Plus Environment

20

Page 21 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(14) Rissanou, A. N.; Power, A. J.; Harmandaris, V. Structural and Dynamical Properties of Polyethylene/Graphene Nanocomposites through Molecular Dynamics Simulations Polymers 2015, 7(3), 390-417. (15) Bacova, P.; Rissanou, A. N.; Harmandaris, V. Edge-Functionalized Graphene as a Nanofiller: Molecular Dynamics Simulation Study. Macromolecules 2015, 48(24), 9024– 9038. (16) Falkovich, S. G.; Larin, S. V.; Lyulin, A. V.; Yudin, V. E.; Kenny, J. M.; Lyulin, S. V. Influence of the Carbon Nanofiller Surface Curvature on the Initiation of Crystallization in Thermoplastic Polymers. RSC Adv. 2014, 4, 48606-48612. (17) Larin, S. V.; Falkovich, S. G.; Nazarychev, V. M.; Gurtovenko, A. A.; Lyulin, A. V.; Lyulin, S. V. Molecular-Dynamics Simulation of Polyimide Matrix Pre-Crystallization Near the Surface of a Single-Walled Carbon Nanotube. RSC Adv. 2014, 4, 830-844. (18) Hegde, M.; Lafont, U.; Norder, B.; Picken, S. J.; Samulski, E. T.; Rubinstein, M.; Dingemans, T. Induced Crystallization in an Amorphous All-Aromatic Poly(ether imide). Macromolecules 2013, 46, 1492-1503. (19) Sharma, S. K.; Prakash, J.; Bahadur, J.; Sudarshan, K.; Maheshwari, P.; Mazumder, S.; Pujari, P. K. Investigation of nanolevel molecular packing and its role in thermomechanical properties of PVA-fMWCNT composites: positron annihilation and small angle X-ray scattering studies. Phys. Chem. Chem. Phys. 2014, 16(4), 1399-1408. (20) Sharma, S. K.; Prakash, J.; Pujari, P. K. Effects of the molecular level dispersion of graphene oxide on the free volume characteristics of poly(vinyl alcohol) and its impact

ACS Paragon Plus Environment

21

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 32

on the thermal and mechanical properties of their nanocomposites. Phys. Chem. Chem. Phys. 2015, 17(43), 29201-29209. (21) Liu, J.; Shen, J. X.; Zheng, Z. J.; Wu, Y. P.; Zhang, L. Q. Revealing the toughening mechanism of graphene-polymer nanocomposite through molecular dynamics simulation Nanotechnology 2015, 26(29), 291003. (22) Ning, N. Y.; Zhang, W.; Zhao, Y. S.; Luo, F.; Fu, Q. Nanohybrid shish kebab structure and

its

effect

on

mechanical

properties

in

poly(L-lactide)/carbon

nanotube

nanocomposite fibers. Polymer International 2012, 61(11), 1634-1639. (23) Somani, R. H.; Yang, L.; Zhu, L.; Hsiao, B. S. Flow-induced shish-kebab precursor structures in entangled polymer melts. Polymer 2005, 46(20), 8587-8623. (24) Gorlier, E.; Haudin, J. M.; Billon, N. Strain-induced crystallisation in bulk amorphous PET under uni-axial loading. Polymer 2001, 42(23), 9541-9549. (25) Hinkley, J. A.; Dezern, J. F.; Feuz, L.; Klinedinst, D. Uniaxial stretching of poly(ketoether-imide) films. Polym. Eng. Sci. 2004, 44(7), 1360-1367. (26) Haopeng, W.; Keum, J. K.; Hiltner, A.; Baer, E.; Freeman, B.; Rozanski, A.; Galeski, A. Confined Crystallization of Polyethylene Oxide in Nanolayer Assemblies. Science 2009, 323, 757-760. (27) Friedrich, K., Schlarb, A. K., Tribology of Polymeric Nanocomposites (Second Edition) Friedrich, K., Schlarb, A. K., Eds.; Elsevier: Oxford, U.K., 2013. (28) Huang, T.; Xin, Y.; Li, T.; Nutt, S.; Su, C.; Chen, H.; Liu, P.; Lai, Z. Modified Graphene/Polyimide Nanocomposites: Reinforcing and Tribological Effects ACS Appl. Mater. Interfaces 2013, 5, 4878-4891.

ACS Paragon Plus Environment

22

Page 23 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(29) Rafiee, M. A.; Rafiee, J.; Srivastava, I.; Wang, Z.; Song, H.; Yu, Z.-Z.; Koratkar, N. Fracture and Fatigue in Graphene Nanocomposites Small 2010, 6, 179-183. (30) Li, L.; Yang, Y.; Yang, G.; Chen, X.; Hsiao, B. S.; Chu, B.; Spanier, J. E.; Li, C. Y. Patterning Polyethylene Oligomers on Carbon Nanotubes Using Physical Vapor Deposition. Nano Lett. 2006, 6, 1007-1012. (31) Hess, C.; Kutzner, Spoel, D. van der; Lindahl, E. J. GROMACS 4: Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. J. Chem. Theory Comput. 2008, 4, 435-447. (32) Spoel, D. van der; Lindahl, E.; Hess, B.; Groenhoff, G.; Mark A. E.; Berendsen, H. J. C. GROMACS: Fast, Flexible, and Free. J. Comput. Chem. 2005, 26, 1701-1718. (33) Lyulin, S. V.; Larin, S. V.; Gurtovenko, A. A.; Lukasheva, N. V.; Yudin, V. E.; Svetlichnyj V. M.; Lyulin, A. V. Effect of the SO2 Group in the Diamine Fragment of Polyimides on Their Structural, Thermophysical, and Mechanical Properties. Polym. Sci. Ser. A. 2012, 54, 631-643. (34) Nazarychev, V. M.; Larin, S. V.; Lukasheva, N. V.; Glova, A. D.; Lyulin, S. V. Evaluation of the Characteristic Equilibration Times of Bulk Polyimides via Full-Atomic Computer Simulation. Polym. Sci. Ser. A. 2013, 55, 570-576. (35) Lyulin, S. V.; Larin, S. V.; Gurtovenko, A. A.; Nazarychev, V. M.; Falkovich, S. G.; Yudin, V. E.; Svetlichnyi, V. M.; Gofman I. V.; Lyulin, A. V. Thermal Properties of Bulk Polyimides: Insights from Computer Modeling versus Experiment. Soft Matter 2014, 10, 1224-1232.

ACS Paragon Plus Environment

23

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 32

(36) Falkovich, S. G.; Lyulin, S. V.; Nazarychev, V. M.; Larin, S. V.; Gurtovenko, A. A.; Lukasheva, N. V.; Lyulin, A. V. Influence of the Electrostatic Interactions on Thermophysical Properties of Polyimides: Molecular-Dynamics Simulations. J. Polym. Sci. Part B: Polym. Phys. 2014, 52, 640-646. (37) Lyulin, S. V. Gurtovenko, A. A.; Larin, S. V.; Nazarychev, V. M.; Lyulin, A. V. Microsecond

Atomic-Scale

Molecular

Dynamics

Simulations

of

Polyimides.

Macromolecules 2013, 46, 6357-6363. (38) Fal'kovich, S. G.; Larin, S. V..; Nazarychev, V. M..; Volgin, I. V..; Gurtovenko, A. A..; Lyulin, A. V.; Lyulin, S. V. Computer Simulation of the Heat-Resistant Polyimides ULTEM™ and EXTEM™ with the Use of GROMOS53a6 and AMBER99 Force Fields. Polym. Sci. Ser. A. 2014, 56, 558-567. (39) Oostenbrink, C.; Villa, A.; Mark, A.E.; van Gunsteren, W. F. A Biomolecular Force Field Based on the Free Enthalpy of Hydration and Solvation: The GROMOS Force-Field Parameter Sets 53A5 and 53A6. J. Comput. Chem. 2004, 25, 1656-1676. (40) Yudin, V. E.; Divoux, G. M.; Otaigbe, J. U.; Svetlichnyi, V. M. Synthesis and Rheological Properties of Oligoimide/Montmorillonite Nanocomposites. Polymer 2005, 46, 10866-10872. (41) H. J. C. Berendsen, in Computer Simulations in Material Science, ed. M. Meyer and V. Pontikis, Kluwer, Dordrecht, 1991.

ACS Paragon Plus Environment

24

Page 25 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(42) Berendsen, H. J. C.; Postma, J. P. M.; van Gunsteren, W. F.; DiNola, A.; Haak, J. R. Molecular-Dynamics with Coupling to an External Bath. J. Chem. Phys. 1984, 81, 36843690. (43) Hess, B.; Bekker, H.; Berendsen, H. J. C.; Fraaije., J. G. E. M. LINCS: a Linear Constraint Solver for Molecular Simulations. J. Comput. Chem. 1997, 18, 1463-1472. (44) Vorselaars, B.; Lyulin, A. V.; Michels, M. A. J. Microscopic Mechanisms of Strain Hardening in Glassy Polymers. Macromolecules. 2009, 42, 5829-5842. (45) Mann, H. B.; Whitney, D. R. On a Test of whether One of Two Random Variables is Stochastically Larger than the Other Ann. Math. Stat. 1947, 18, 50-60. (46) Yu, K.; Wang, P.; Lu, G.; Chen,K.; Bo, Z.; Chen, Z. Patterning Vertically Oriented Graphene Sheets for Nanodevice Applications. J. Phys. Chem. Lett. 2011, 2, 537-542. (47) Hunter, C.; Sanders J. K. M. The Nature of π-π Interactions J. Am. Chem. Soc. 1990, 112, 5525-5534. (48) Kontsos, A.; Cuadra, J. A. Multiscale Stochastic Finite Elements Modeling of Polymer Nanocomposites in Modeling and Prediction of Polymer Nanocomposite Properties. Mittal, V. (ed) Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim, 2013. (49) Wu, S. W.; Tang, Z. H.; Guo, B. C.; Zhang, L. Q.; Jia, D. M. Effects of interfacial interaction on chain dynamics of rubber/graphene oxide hybrids: a dielectric relaxation spectroscopy study. RSC Adv. 2013, 3(34), 14549-14559.

ACS Paragon Plus Environment

25

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 32

(50) Klonos, P.; Kripotou, S.; Kyritsis, A.; Papageorgiou, G. Z.; Bikiaris, D.; Gournis, D.; Pissis, P. Glass transition and segmental dynamics in poly(L-lactic acid)/graphene oxide nanocomposites. Thermochem. Acta 2015, 617, 44-53. (51) Zhang, J. J.; Sun, Y. Y.; Li, D. S.; Cao, Y.; Wang, Z.; Ma, J.; Zhao, G. Z. Modeling the mechanics of graphene-based polymer composite film measured by the bulge test. J. Phys. D Appl. Phys. 2015, 48(42), 425302. (52) Lin, C. S.; Horng, T. L.; Chen, J. H.; Chen, K. H.; Wu, J. J.; Chen, C. Y.; Ma, S. H. Mechanical Properties Measurement of Polymer Films by Bulge Test and Fringe Projection. Adv. Mater. Sci. Eng. 2014, 170279. (53) Roussou, R.; Karatasos, K. Graphene/poly(ethylene glycol) nanocomposites as studied by molecular dynamics simulations. Materials and Design 2016, 97, 163–174. (54) Yan, L. T.; Xie, X. M. Computational modeling and simulation of nanoparticle selfassembly in polymeric systems: Structures, properties and external field effects. Prog. Polym. Sci. 2013, 38(2), 369-405.

TOC GRAPHICS

(left) a snapshot of the first layer of R-BAPB polyimide near graphene (right) stress strain curves of R-BAPB, structurally ordered by the graphene

ACS Paragon Plus Environment

26

Page 27 of 32

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

27

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 1. Chemical structure of the repeated unit of R-BAPB polyimide (upper part) and part of the infinite graphene sheet7 (lower part). The red and blue arrows indicate the vectors corresponding to the rigid fragment of the R-BAPB repeated unit chosen for structural analysis. 81x40mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 28 of 32

Page 29 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 2. Density profiles along Z axis for R-BAPB with removed graphene before pulling (black curve) and at 10% strain (red curve) 82x57mm (300 x 300 DPI)

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 3. Stress-strain curves averaged over three directions (along the axes X, Y and Z) both for R-BAPB samples ordered by graphene and for amorphous R-BAPB samples 59x42mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 30 of 32

Page 31 of 32

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 4. Averaged stress-strain curves of (a) non-filled amorphous R-BAPB samples and (b) grapheneordered R-BAPB samples, obtained after deforming along axes X (black curves), Y (red curves), and Z (blue curves) 61x22mm (300 x 300 DPI)

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

TOC graphics 39x19mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 32 of 32