Theoretical and Experimental Insights into the Phase Transition of

Subscriber access provided by READING UNIV

Article

Theoretical and Experimental Insights into the Phase Transition of Rubber/Plastic Blends during Dynamic Vulcanization Zhaoyang Wei, Shangqing Li, Nanying Ning, Ming Tian, Liqun Zhang, and Jianguo Mi Ind. Eng. Chem. Res., Just Accepted Manuscript • DOI: 10.1021/acs.iecr.7b03452 • Publication Date (Web): 20 Oct 2017 Downloaded from http://pubs.acs.org on October 24, 2017

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.

Industrial & Engineering Chemistry Research 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 27

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

Industrial & Engineering Chemistry Research

Theoretical and Experimental Insights into the Phase Transition of Rubber/Plastic Blends during Dynamic Vulcanization Zhaoyang Wei,a,b Shangqing Li,a,b Nanying Ning,a,b Ming Tian,a,b,* Liqun Zhang,a,b and Jianguo Mia,*

a

State Key Laboratory of Organic−Inorganic Composites, Beijing University of Chemical Technology, Beijing 100029, China

b

Beijing Advanced Innovation Center for Soft Matter Science and Engineering, Beijing University of Chemical Technology, Beijing 100029, China

ABSTRACT: While the phase inversion phenomena in thermoplastic vulcanizates have been extensively investigated, the mechanism underlying the evolution from the sea−island to island−sea phase transition is far from being understood. In this work, such transition of the polyolefin elastomer (POE)/polypropylene (PP) blend was studied simultaneously by the density functional theory and experiments. Before vulcanization, the independent POE and PP phases were mixed into a homogenous blend. After vulcanization, the microstructure of POE was altered, leading to a metastable state of the blend. Thereafter, the blend decomposed into two new phases. Accordingly, the interfacial density changed with the sea−island morphology transition. According to the nucleation thermodynamics, the predicted nanoparticle size of POE nuclei is in general agreement with our experimental measurement. Thus, this work provides a novel thermodynamic explanation regarding the phase inversion in thermoplastic vulcanizates.

Corresponding author: [email protected]; [email protected]. 1

ACS Paragon Plus Environment


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 27

1. Introduction

Thermoplastic vulcanizates (TPVs) are well known as a group of high-performance thermoplastic elastomers.1-3 Involving the in-situ crosslinking of high rubber content and the simultaneous mixing with low thermoplastic content, the dynamic vulcanization (DV) is a reactive blending process above the melting temperature under high shear to realize the phase inversion. This class of polymer blends is of interest in applications ranging from structural components for aerospace to electronics packaging.4-6 The microstructures of TPVs, including particle sizes and morphologies of rubber phase, plays a remarkable role on their mechanical properties and processabilities.7-12 Therefore, understanding the molecular origins of these properties is becoming increasingly important for experiments aimed at identifying tailored materials for specific applications. Despite the apparent simplicity of the DV process, TPVs are highly complex and involve many chemical and physical processes. In the past few decades, extensive researches focused on the morphology evolution of TPVs. Some researchers have shown that the crosslinking of rubber phase leading to rapid enhancements of viscosity and shear stress; thus the continuous rubber phase elongates to balance the interfacial forces and breaks up into a multitude of particles dispersed in the plastic matrix, resulting in the phase inversion.13-21 Recently, we further confirmed that the particles are composed of agglomerates of nanoparticles in an ethylene propylene diene monomer (EPDM)/polypropylene (PP) TPV.22,23 However, the evolution mechanism is far from being understood. From a thermodynamic point of view, molecular recombination requires only modest thermodynamic costs, which is easy to occur. However, the recombination interface involves a large amount of

2


Page 3 of 27

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


thermodynamic cost. If the interface remains invariant, the cost of interfacial free energy should be compensated by other energy that facilitates the phase separation. The energy cost and compensation grow linearly with the surface area and volume, respectively. 24 This viewpoint is very common in the classical nucleation theory, which points out that the formation of a critical nucleus is necessary for the phase separation of supersaturated solutions.25 In order to understand the evaluation from the sea−island to island−sea structure of TPVs, the phase transition should be particularly illustrated. During a DV process, the structure of rubber undergoes variations, which plays a vital role in the phase inversion. Such inversion is a dynamic process accompanied by thermodynamic nucleation and phase transition. Developing a theoretical approach that can capture the essential mechanism of phase transition from a metastable to stable state is an attractive goal. Previous investigations employed the well-known Flory−Huggins theory to interpret the phase behavior in polymer blends.26-29 Good fittings to experimental data were achieved by the proper choice of the empirical Flory−Huggins parameter. However, the theory was restricted as a semi-empirical model since the effect of polymer chain conformation was overlooked. In the present work, the phase transition and interfacial evolution in polyolefin elastomer (POE)/PP TPV were investigated using the classical three-dimensional density functional theory (3D-DFT).30-33 Meanwhile, the corresponding experiments were conducted to assess the predicted results. In general, POE has relatively simple structural feature, which is particularly suitable for theoretical characterization and could serve as an excellent benchmark for the phase transition of TPVs. We aimed to (i) understand the molecular origin of the observed POE nanoparticles through their microstructure evolution and (ii) predict the variations of interfacial tension as well as the critical nucleation size. Seeking out a reasonable expression for the grand potential based on 3



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 27

different interfacial microstructures depending on different interactions between POE and PP is the key issue of the DFT approach. Through minimization of the grand potential, the interfacial structures and properties were evaluated to decipher the phase inversion mechanism. Meanwhile, the biphasic morphology and actual sizes of rubber nanoparticles were observed in experiments. Finally, the theoretical predictions were tested with the experimental values to assess the reliability of the DFT model.

2. Experiment section 2.1 Materials Two types of commercially available polymers were used and are listed in Table 1. Dicumyl peroxide (DCP) was used as the crosslinking agent, which was purchased from Akzo Nobel Polymer Chemicals. A commercially available pentaerythritol tetrakis 3-(3,5-ditert-butyl-4-hydroxyphenyl) propionate (1010) was applied as the antioxidant. Cyclohexane and xylene were purchased from Beijing Chemical Works, China. Table 1. Characteristics of the materials used Polymer/Properties

POE

PP

Density (g/cm3)

0.86

0.91

MFI (g/10 min)

0.5a

0.5b

Tm (°C)

—

159

Octene content (wt %)

27

—

Mooney viscosity

33

—

4


Page 5 of 27

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


ML (1+4) 125 °C Grade Supplier

Engage8150

HP500D

DuPont Dow Elastomers Co., Ltd. (USA) DOW Co., Ltd (China)

MFI − Melt flow index; Tm − Melting point. a

at 190 °C and 2.16 kg, b at 230 °C and 2.16 kg.

2.2 Preparation of the samples The POE/PP TPV was prepared by melt mixing PP with POE in a Haake internal mixer (Haake Rheomix 600 OS, Thermo Fisher Scientific, America) equipped with two counter-rotating rotors. Firstly, the two composites of POE and PP with the mass ratio of 60/40 were fed into a two-roll mill at 180 °C. At initial stage, they formed an inhomogeneous premix. After long time mixing, the premix was transformed to a homogeneous blend. Meanwhile, antioxidant 1010 (0.35 wt % based on the POE) was added to prevent PP from aging. The sample with addition of DCP (1.5 wt % based on the POE) was then transferred to another two-roll mill at ambient temperature. Finally, the blend was fed into the Haake Rheomix under the condition of 170 °C and 100 rpm. Meanwhile, the blend decomposed into two new phases.

2.3 Interfacial tension measurement To determine the surface tension ( γ ), the contact angles for the samples in the form of films was measured using a Kruss instrument (model DCAT21; Germany). The interfacial tension ( γ 12 ) between PP and POE was therefore calculated from the harmonic mean equation15 5



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

γ 12 = γ 1 + γ 2 −

Page 6 of 27

4γ 1dγ 2d 4γ 1pγ 2p − γ 1d + γ 2d γ 1p + γ 2p

(1)

where γ 1 and γ 2 represent the surface tensions of POE and PP, respectively. γ d and γ p are the associated nonpolar and polar components. At high temperatures, the values of γ d and γ p were calculated from the χ value, which is the ratio of γ p to γ and does not change with temperature ( χ = γ p / γ , dχ / dT = 0 ).15, 16

2.4 Morphology investigation The morphologies of the samples were observed using a peak force atomic force microscopy (PF-AFM) (Nanoscope IIIa, Bruker Corporation, Karlsruhe, Germany). First of all, the samples were polished using a cryo-ultramicrotome (Leica EM UC7; Wetzlar, Germany) at −130 °C to obtain a flat surface for AFM characterization. The size of the individual POE rubber particles was also observed by the PF-AFM. The dynamically vulcanized POE/PP TPV was dissolved in hot xylene for 96 h. The solvent was refreshed every 24 h to wash off the dissolved PP. The obtained product was coated onto a clean coverslip using a spin coater (KW-4A; China). The obtained sample was then observed using PF-AFM. The POE particle size was counted with the Image-Pro Plus 4.5 software.

3. Theoretical section Within the 3D-DFT framework, the interactions between different sites in POE and PP molecules were described by the universal MARTINI force field,34 in which the groups of four or three units were represented by a single coarse-grained (CG) bead. For the ethylene and octene contents, we applied standard 4-to-1 mapping.

6


Page 7 of 27

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


For the propene content, however, we applied 3-to-1 mapping because the CH2 groups can be shared by the adjacent CG beads. Therefore, a POE molecule was described using a series of bonded C1 beads. Meanwhile, a PP molecule was simplified into a series of bonded SC1 beads. 34, 35 Figure 1 shows the Coarse-grained POE and PP. The POE and PP chains were composed of 200 C1 and 200 SC1 beads, respectively. Our calculations revealed that once the chain length reached 200, the effect of the chain length on the polymer structure was trivial. The pairwise nonbonded interactions were given by the Lennard−Jones (LJ) potential. As a consequence, the united intermolecular potentials were calculated by summing the pair potentials   σ αα ′ 12  σ αα ′ 6  ∑ 4ε αα ′   −   U (r ) =   r    r   ∞ 

r ≥ σ αα ′

(2)

r < σ αα ′

where ε αα ′ and σ αα ′ are LJ potential parameters. The levels of interactions between different CG sites are given in Table 2. The cross-interaction parameters were estimated using the mixing rules with ε αα ′ = ε α ε α ′ and σ αα ′ = (σ α + σ α ′ ) / 2 .

7



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 27

Figure 1. United atoms to CG maps of POE and PP (red and blue circles). C1 and SC1 are standard CG types. Table 2. CG potential interaction parameters34, 35 Site

σ (Å)

ε (KJ/mol)

C1

4.7

3.5

SC1

4.3

2.6

The grand potential Ω[{ρα (r)}] of polymer blend system, can be expressed as

Ω[{ρα (r)}] = kBT ∑ ∫ drρα (r) ln ( ρα (r) ) − 1 + F ex [{ρα (r)}] − kBT ∑ ∫ ρα (r) µα dr α

(3)

α

where F ex [{ρα (r )}] is the excess free energy contributions. α denotes a site in polymer molecules, kB is

the Boltzmann constant and T is the temperature. µα represents the bulk chemical potential of component

α . The excess free energy is presented as follows 8


Page 9 of 27

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


F ex [{ρα (r )}] = F hs [{ρα (r )}] + F att [{ρα (r )}] + F chain [{ρα (r )}]

(4)

where three terms on the right-hand side represent hard-sphere repulsion, dispersive attraction, and chain connectivity, respectively. The contribution of hard-sphere repulsion can be described by the fundamental measure theory36

F hs [{ρα (r )}] = kBT ∫ dr Φ hs [nm (r )]

(5)

hs where Φ [ nm (r )] is the free energy density, which stems from the modified fundamental measure theory37

and includes both the scalar and vector contributions



n1n2 − nV1 ⋅ nV2



1 − n3

Φ hs [ nm (r ) ] =  −n0 ln (1 − n3 ) +

+

1 

 n3 ln (1 − n3 ) +

36π 

 n23 − 3n2nV2 ⋅ nV2   2  n33 (1 − n3 )   n32

(6)

where nm (r) represents the weighted densities with m = 0,1, 2, 3, V1 ,V2 . The details of nm (r) can be found in the literature.37 For the attractive contribution, F att [{ρα (r)}] can be described using the weighted density approximation method38

F att [{ρα (r)}] = kBT ∑ ∫ ( ρα (r ){F1[ ρα (r )] + F2 [ ρα (r )]})dr

(7)

α

where the expressions of F1[ ρα (r)] and F2 [ ρα (r)] can be given by the first-order mean spherical approximation expansion38

  1 + z1,αα ' Rαα ' z 'R ' F1[ ρα (r)] = −2πρα (r)β ∑∑ xα xα ' εαα ' k1,αα '  G0,αα ' ( z1,αα ' )e 1,αα αα −  z1,2αα '  α α'   1 + z2,αα ' Rαα ' z 'R ' − k2,αα '  G0,αα ' ( z2,αα ' )e 2,αα αα −  z2,2 αα ' 

 3  + 8πρα (r)β ∑∑ xα xα ' εαα ' Rαα 'I αα ' ,∞  ' α α 

3 −8πρα (r)β ∑∑ xα xα ' εαα ' g0,αα ' ( Rαα ' ) Rαα 'I αα ' ,1

α

   

α'

9


(8)


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

(

F2[ρα (r)] =−πρα (r)β∑∑xα xα' εαα' k1,αα' G1,αα' (z1,αα' )e 1,αα  α α' z

' Rαα'

) −k ( G 2,αα

'

Page 10 of 27 z

' Rαα'

(z )e 2,αα 2,αα 2,αα' '

3 −4πρα (r)β∑∑xα xα' εαα' g1,αα' (Rαα' )Rαα 'I ' αα ,1

)

(9)

α α'

with 12

1σ ' =  αα 9  Dαα '

 1σ '  −  αα 3  Dαα ' 

1σ ' Iαα ' ,1 =  αα 9  Dαα '

 1σ '  −  αα 3  Dαα ' 

Iαα ' ,∞

12

6

  ,  6

 2  σ αα '  +   9  Dαα '

  

(10)

3

where xα is mole fraction of component α , z1,αα ' = 2.9637 / σ αα ' , z2,αα ' = 14.0167 / σ αα ' . G0 and G1 are the radial distribution functions at hard-sphere contact. The details of calculation of F1[ ρα (r)] and F2 [ ρα (r)] can be seen elsewhere. 38, 39 The weighted density ρα (r) is defined as att ρα (r ) = ∑ ∫ ρα (r')ωαα ( r − r' )dr' , '

(11)

'

α'

att att att att with the weight function ωαα ' (r ) = uαα ' ( r ) / ∫ uαα ' ( r )dr . Here uαα ′ (r) is the LJ attractive interaction potential

between any two sites of species

α and α ′ .

Chain connectivity is considered as a series of bonding segments by providing a label to each segment and allowing the segments to be exclusively integrated with their particular matching segments.40 Therefore, the contribution of chain connectivity to the free energy functional is presented as follows N

F chain [{ ρ α (r )}] = k BT ∫ d r ′∑ ρ α (r ′) × α =1

 X Aα (r ′) 1  α +   ln X A (r ′) − 2 2 A ∈Γ ( α ) 

∑

(12)

where the first summation is for the overall segments α , the second one is for the overall association sites on segment α as Γ(α ) , and the variable X Aα stand for the fraction of segments α that are not bonded at association site A. The details of the above functional derivatives can be seen elsewhere. 41, 42 In fact, more

10


Page 11 of 27

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


association schemes, in which polymer molecules may involve multiple segments or multiple sites, may lead to more complex polymer structures.43 Here the multiple associating segments of POE were determined by the degree of crosslinking as follows: (1) Measure the volume swell ratio Q in the experiment. The details of measurement of Q can be seen our previous study23; (2) Estimate the molecular weight M c between the crosslinks using the Flory-Rehner equation44; (3) Calculate the multiple associating segments with n = M n / M c , where M n is the number average molecular weight of PEO. The equilibrium density profiles can be obtained by ∂Ω / ∂ρα ( r ) = 0 . Then, the Euler–Lagrange Equation can be generalized as

 δ ( F hs [{ρα (r )}] + F att [{ρα (r )}] + F chain [{ρα (r )}])   ρα (r ) = exp  βµα − β   δρα (r )  

(13)

The interfacial tension γ of POE−PP was calculated from the grand potential, which can be written as

γ=

∆Ω Ω [ ρ (r)] − Ω ( ρb ) = A A

(14)

where Ω ( ρ b ) is the grand potential of the bulk phase, A is the interfacial area. According to the liquid–liquid phase equilibrium thermodynamics, the chemical potential of each component has to be the same in different phases. For a binary mixture, the equations for calculating the chemical potentials can be written as  ∂µ1  d µ1      ∂P  =  d µ   ∂µ 2  2  ∂P

∂µ1   dp  ∂y1    =  ∂µ 2     dy1  ∂y1 

(15)

where d µi ( i = 1, 2) means the differential of µi . 11



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 27

The iteration method for obtaining the phase equilibrium at given temperature and pressure is as follows：(1) Set the premix as 1.0 mole, and the mole quantity of POE-rich phase is LPOE-rich = 0.995 ; (2) Assume different mole

fractions

of

POE

in

POE-rich

and

PP-rich

phases

( x POE

and

y POE )

to

satisfy

z POE = x POE LPOE-rich + y POE (1 − LPOE-rich ) , where z POE is the constant mole fraction of POE in the premix; (3) Substitute the two values into eq.15, solve the equation to obtain new x′POE and y ′POE ; (4) Repeat step 3 and estimate the chemical potentials of POE and PP in the two phases. (5) Evaluate the chemical potentials. If POE-rich PP-rich POE-rich PP-rich µPOE − µPOE > 10 −3 , or µPP − µPP > 10 −3 , set a new mole quantity of POE-rich phase with

Lnew POE-rich = LPOE-rich − 0.005 ; (6) Perform the iteration procedure including the steps of 2, 3, and 4 until POE-rich PP-rich POE-rich PP-rich µPOE − µPOE < 10−3 and µPP − µPP < 10−3 . Accordingly, the phase equilibrium at given temperature

and pressure are determined. Then change the temperature to obtain a new equilibrium. As such, the temperature-composition phase diagram is finally determined. According to nucleation thermodynamics, nucleation occurred under a supersaturated condition. A growing embryo should overcome its free energy barrier before reaching a stable state. When the free energy reaches its maximum, the droplet size is defined as the critical size ( d c ), and the maximum is the energy barrier. The supersaturation ratio can be defined as S = ρsuper / ρpoor , where ρsuper and ρ poor denote the averaged and bulk densities of POE in the POE-poor phase, respectively. The averaged density was obtained by integrating the density distribution in the nucleation region. The initial droplet density profile can be written as

 ρrich

ρ (r ) = 

 ρsuper

0 ≤ r ≤ R0

(16)

r > R0

12


Page 13 of 27

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


for the iteration calculations. Here R0 is the initial radius of the nucleus, ρ rich is the bulk density of POE in POE-rich phase. In order to avoid the influence of the edge effect, the interfacial densities were calculated on a 3D grid of

256 × 256 × 256 points. 0.1σ was set as the grid spacing in each direction. The outmost site of the surface was arranged in the z coordinate of the polymer surface sites and was located at the point where z = 0 . The traditional Picard-type iterative method was employed to calculate the density profiles. For description of interfacial density distributions in theoretical model, the external force fields given by the pure PP and POE/PP blend were constructed using molecular dynamics simulations to mimic the real surfaces. In these simulations, all molecule parameters were derived from the MARTINI force field.

34, 35

More

importantly, the force field parameters and contents of POE and PP in the simulations and theoretical calculations were exactly the same to ensure their consistency. First of all, the elongated polymer chains with extremely low density were placed in a large periodic box to obtain the initial configuration. Then a uniform dispersion could be obtained by running the NVT ensemble for a long time. The system was then compressed to the actual density by the NPT ensembles. At atmospheric pressure, the NVT simulations run for at least 6 ns to ensure system balance. Finally, the structural information of equilibrated system was collected by running of (10-20)×106 NVT simulation steps at 443 K. During all simulation runs, 1.0 fs was chosen as the time step for the multiple time step integrator.45 A Nose-Hoover thermostat at the coupling frequency of 0.02 fs-1 was used to perform the isothermal calculations. The pressure was controlled by the Nose-Hoover barostat.46 A SHAKE algorithm was used to constrain the bond lengths.47 As a consequence, the external force fields could be represented as the slices of the equilibrium systems. 13



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 27

3. Results and discussion The initial morphology of the premix was firstly determined, at which PP (denoted as island) was dispersed in POE (denoted as sea). Figure 2 shows the spatial density profile of POE on the PP surface to display the interfacial microstructure. The POE molecules display pronounced and disordered packing on the PP surface. The strongest interactions, which were ascribed to the close rows of the PP chains, resulted in the highest peaks at the slots of the surface. On the other hand, one can see that both ordered and disordered voids were dispersed on the rough PP surface. These voids were induced by the synergistic effect of chemical characteristics and topological structures. When the POE molecules were in contact with the PP matrix, some of the voids were occupied by POE molecules because their sizes were approximately equal.

Figure 2. Two-dimensional cut of the spatial density distribution of the sea (POE) on the island (PP) surface. The distribution is cut at y = 0 .

14


Page 15 of 27

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


Due to the similarity in the chemical structure of POE and PP, they can form a compatible blend. 48, 49 From the thermodynamics point of view, the molecular weight of POE increases with the vulcanization of POE, which reduces the mixing entropy and leads to the phase decomposition.2 Figure 3 presents the calculated temperature−composition phase diagram for the vulcanized POE/PP at constant pressure. In descending order of temperature, the POE/PP blend revealed a complex phase behavior involving a lower critical solution temperature (LCST). The vertical dashed line indicated the mole fraction of POE in the homogenous blend in the current experiment. At the experimental temperature of phase separation, the POE-rich phase (denoted as island) and the PP-rich phase (denoted as sea) were formed. Therefore, the composition of the two phases was determined based on the requirement of equal chemical potential and pressure in the coexistent phases.

Figure 3. Temperature−composition phase diagram of vulcanized POE (1) and PP (2) at constant pressure. The vertical dashed line denotes the mole fraction of POE in the blend. x POE is the mole fraction of POE in the POE-rich phase, and y POE is the mole fraction of POE in the PP-rich phase.

15



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 27

After phase separation at 443 K, the initial sea−island morphology in the premix had been transferred to the new island−sea one. Accordingly, the newly interfacial microstructure was determined. Structure analysis is an effective tool to interpret a wide range of interfacial phenomenon. Figure 4 shows the corresponding density distributions of the sea (PP-rich) on the island (POE-rich) surface. One can immediately see the obvious multilayer adsorptions of POE and PP in the POE-rich region. The density distributions close to the interface are high peaks, which can be attributed to the enhancement of interfacial adsorption. Outside the interface, the density distributions vary from peak to trough, implying the transformation of local adsorption/desorption. The comparison of the subfigures 4a and 4b implies that more adsorption peaks emerge in Figure 4a, since the attractions of POE-POE are stronger than those of POE-PP,.

Figure 4. Two-dimensional cuts of the spatial density distributions of the sea (PP-rich phase) containing POE (a) and PP (b) on the island (POE-rich phase) surface. The distributions are cut at y = 0 .

The above density fluctuations correlate to interfacial energy or interfacial tension, which was calculated by eq 14. It is unambiguous that the variations of density profiles, including the height and width of the fluctuations, affect the interfacial tension. Table 3 presents the experimental and calculated results. For the 16


Page 17 of 27

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


premixed system, the calculated value generally agrees with the experimental one, which indicates that the theory provides quantitatively reasonable description of the interfacial microstructures and properties of POE/PP blend. Table 3. Interfacial tensions of POE/PP blend

γ (mN/m) Species Theory

Experiment

POE/PP (premixed)

0.47

0.50 ± 0.05

POE/PP (vulcanized)

0.34

—

Figure 5 shows the experimental phase images of POE/PP blend. The lighter and darker regions represent the POE and PP phases, respectively. The initial sea−island morphology in the premix can be observed in Figure 5a. As expected, the inhomogeneous PP domains, including large or small PP regions, were dispersed in the POE matrix. In Figure 5b, one can see a new island−sea structure after the DV disposal, indicating the occurrence of liquid-liquid phase decomposition. Since the POE-rich domains were dispersed in the PP-rich phase, and some PP crystals were located in the crosslinked POE phase, the morphology of TPV product was quite complex.50

17



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 27

Figure 5. AFM phase images of the POE/PP blend taken during an (a) premixed and (b) dynamic vulcanized. The lighter and darker regions represent the POE and PP phases, respectively.

Figure 6 displays the formed POE nuclei from the theoretical calculations and experimental measurements. The size was determined using the interfacial density distributions. During the nucleation process, a droplet that is growing in size should overcome its free energy barrier before reaching a stable state in the system. The critical size ( d c ) was defined as the droplet size at the point where the free energy needed for the formation of POE-rich droplet is maximum. When an embryo has a size less than the critical value, any increase of the size is accompanied by additional free-energy consumption, which means the embryo would shrink. On the other hand, when an embryo has a size greater than the critical value, any increase in size is accompanied by the depressed free-energy, and the embryo grows. Figure 6a shows an AFM image of the POE/PP sample, which was coated on a coverslip. To characterize the individual POE-rich nanoparticle, the PP content was washed off with hot xylene. The lighter and darker regions represent the glass substrate and the crosslinked POE-rich phase, 18


Page 19 of 27

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


respectively. The minimum size of the POE-rich nanoparticles is 20 nm. Figure 6b displays the calculated excess free energy curve during the formation of POE-rich droplets. it is shown that the size of POE-rich droplets is 22 nm. In general, the theoretical prediction is in consistent with the experimental results. The reason lies in the fully homogeneous mixing, such that the thermodynamic bases of the theoretical and experimental phase transitions are the same. The differences between the theoretical prediction and experiment measurement could be caused by (a) the lack of details of topologic structures due to the coarse-grained disposal of various groups and (b) the coalescence of nanoparticles during the DV process.

Figure 6. (a) AFM phase image of the dissolved the POE/PP blend. (b) Dimensionless free energy barrier of the growing POE-rich droplets.

4. Conclusion The 3D-DFT calculations and AFM measurements were combined to investigate the phase separation of POE/PP blend during DV. The phase equilibrium and interfacial properties were quantitatively predicted to 19



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 20 of 27

illustrate the morphology transition from sea−island to island−sea in the blend. The results indicate that the phase transition is induced by thermodynamic driving force which is developed during the dynamic crosslinking process. In other words, the dynamic crosslinking causes the microstructure variations of POE, which results in the reduced interfacial tension, and finally the liquid-liquid phase transition. The observed POE nanoparticles nucleated into new phases, and the nucleus size was demonstrated by the 3D-DFT approach. The agreement between the theory and experiment enables a fundamental understanding of the phase inversion of the TPV.

Acknowledgments The financial support from the National Basic Research Program of China (Grant No. 2015CB654700 (2015CB674704)) and the National Natural Science Foundation of China (Nos. 51525301 and 51521062) is gratefully acknowledged. We also express our sincere thanks to the CHEMCLOUDCOMPUTING of Beijing University of Chemical Technology.

20


Page 21 of 27

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


References (1) Babu, R. R.; Naskar, K. Recent Developments on Thermoplastic Elastomers by Dynamic Vulcanization. In

Advanced Rubber Composites; Heinrich, G., Ed.; Springer Berlin Heidelberg: Berlin, Heidelberg, 2011; pp 219−247. (2) Abolhasani, M. M.; Zarejousheghani, F.; Naebe, M.; Guo, Q. Does dynamic vulcanization induce phase separation? Soft Matter 2014, 10, 5550−5558. (3) Duin, M. Recent developments for EPDM-based thermoplastic vulcanizates. Macromol. Symp. 2006, 233, 11−16. (4) Li, Y.; Wu, X.; Song, J.; Li, J.; Shao, Q.; Cao, N.; Lu, N.; Guo, Z., Reparation of recycled acrylonitrilebutadiene-styrene by pyromellitic dianhydride: Reparation performance evaluation and property analysis.

Polymer 2017, 124, 41−47. (5)Zhang, K.; Li, G.-H.; Feng, L.-M.; Wang, N.; Guo, J.; Sun, K.; Yu, K.-X.; Zeng, J.-B.; Li, T.; Guo, Z.; Wang, M.,

Ultralow

percolation

threshold

and

enhanced

electromagnetic

interference

shielding

in

poly(l-lactide)/multi-walled carbon nanotube nanocomposites with electrically conductive segregated networks.

J. Mater. Chem. C 2017, 5, 9359−9369. (6)Zhang, X.; Gao, H.; Guo, M.; Li, G.; Liu, Y.; Li, D., A study on key technologies of unmanned driving.

CAAI Transactions on Intelligence Technology 2016, 1, 4−13. (7) Mousa, A. The effect of dynamic vulcanization on the mechanical properties of EPDM/PP thermoplastic elastomers. Int. J. Polymer. Mater. 2005, 54, 619−631.

21



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 27

(8) Banerjee, S. S.; Bhowmick, A. K. Tailored nanostructured thermoplastic elastomers from polypropylene and fluoroelastomer: morphology and functional properties. Ind. Eng. Chem. Res. 2015, 54, 8137−8146. (9) Naderi, G.; Lafleur, P. G.; Dubois, C. Microstructure‐properties correlations in dynamically vulcanized nanocomposite thermoplastic elastomers based on PP/EPDM. Polym. Eng. Sci. 2007, 47, 207−217. (10) l'Abee, R. M. A.; van Duin, M.; Spoelstra, A. B.; Goossens, J. G. P. The rubber particle size to control the properties-processing balance of thermoplastic/cross-linked elastomer blends. Soft Matter 2010, 6, 1758−1768. (11) Liang, B.; Ji, J. M.; Hong, X. D. Study on mechanical properties of EPDM/HPVC prepared by dynamic vulcanization. Adv. Mater. Res. 2011, 299, 751−754. (12) Tang, Y.; Lu, K.; Cao, X.; Li, Y. Nanostructured thermoplastic vulcanizates by selectively cross-linking a thermoplastic blend with similar chemical structures. Ind. Eng. Chem. Res. 2013, 52, 12613−12621. (13)Ellul, M. D.; Tsou, A. H.; Hu, W. Crosslink densities and phase morphologies in thermoplastic vulcanizates.

Polymer 2004, 45, 3351−3358. (14) Coran, A. Y.; Patel, R. Rubber-thermoplastic compositions. Part I. EPDM-polypropylene thermoplastic vulcanizates. Rubber Chem. Technol. 1980, 53, 141−150. (15) Wu, S. Formation of dispersed phase in incompatible polymer blends: Interfacial and rheological effects.

Polym. Eng. Sci. 1987, 27, 335−343. (16) Goharpey, F.; Katbab, A.; Nazockdast, H. Mechanism of morphology development in dynamically cured EPDM/PP TPEs. I. Effects of state of cure. J. Appl. Polym. Sci. 2001, 81, 2531−2544.

22


Page 23 of 27

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


(17) Goharpey, F.; Katbab, A.; Nazockdast, H. Formation of rubber particle agglomerates during morphology development in dynamically crosslinked EPDM/PP thermoplastic elastomers. Part 1: effects of processing and polymer structural parameters. Rubber Chem. Technol. 2003, 76, 239−252. (18) Goharpey, F.; Nazockdast, H.; Katbab, A. Relationship between the rheology and morphology of dynamically vulcanized thermoplastic elastomers based on EPDM/PP. Polym. Eng. Sci. 2005, 45, 84−94. (19) Martin, G.; Barrès, C.; Sonntag, P.; Garois, N.; Cassagnau, P. Morphology development in thermoplastic vulcanizates (tpv): Dispersion mechanisms of a pre-crosslinked epdm phase. Eur. Polym. J. 2009, 45, 3257−3268. (20) Antunes, C. F.; Machado, A.; Van Duin, M. Morphology development and phase inversion during dynamic vulcanisation of EPDM/PP blends. Eur. Polym. J. 2011, 47, 1447−1459. (21) Antunes, C.; Van Duin, M.; Machado, A. Effect of crosslinking on morphology and phase inversion of EPDM/PP blends. Mater. Chem. Phys. 2012, 133, 410−418. (22) Wu, H.; Tian, M.; Zhang, L.; Tian, H.; Wu, Y.; Ning, N. New understanding of microstructure formation of the rubber phase in thermoplastic vulcanizates (TPV). Soft Matter 2014, 10, 1816−1822. (23) Wu, H.; Tian, M.; Zhang, L.; Tian, H.; Wu, Y.; Ning, N.; Chan, T. W. New understanding of morphology evolution of thermoplastic vulcanizate (TPV) during dynamic vulcanization. ACS Sustainable Chem. Eng. 2014,

3, 26−32. (24) Chandler, D. Interfaces and the driving force of hydrophobic assembly. Nature 2005, 437, 640−647. (25) Barrat, J. L.; Hansen, J. P. Basic Concepts for Simple and Complex Liquids; Cambridge University Press, 2003. 23



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 27

(26) Flory, P. J. Thermodynamics of High Polymer Solutions. J. Chem. Phys. 1942, 10, 51−61. (27) Huggins, M. L. Some Properties of Solutions of Long-chain Compounds. J. Phys. Chem. 1942, 46, 151−158. (28)Gai, J.-G.; Li, H.-L.; Schrauwen, C.; Hu, G.-H., Dissipative particle dynamics study on the phase morphologies of the ultrahigh molecular weight polyethylene/polypropylene/poly(ethylene glycol) blends.

Polymer 2009, 50, 336−346. (29)Gai, J.-G.; Hu, G.-H.; Li, H.-L.; Zhu, S.-P.; Hoppe, S., Dissipative Particle Dynamics and Flory−Huggins Theories for Predicting the Rheological Behavior of Ultrahigh Molecular Weight Polyethylene Blends. Ind. Eng.

Chem. Res. 2010, 49, 11369−11379. (30) Douglas Frink, L. J.; Salinger, A. G. Two- and Three-Dimensional Nonlocal Density Functional Theory for Inhomogeneous Fluids. J. Comput. Phys. 2000, 159, 407−424. (31) Ravikovitch, P. I.; Neimark, A. V. Density Functional Theory of Adsorption in Spherical Cavities and Pore Size Characterization of Templated Nanoporous Silicas with Cubic and Three-Dimensional Hexagonal Structures. Langmuir 2002, 18, 1550−1560. (32)Chang, Y.; Mi, J.; Zhong, C. Density functional theory for carbon dioxide crystal. J. Chem. Phys. 2014, 140, 204706. (33) Knepley, M. G.; Karpeev, D. A.; Davidovits, S.; Eisenberg, R. S.; Gillespie, D. An efficient algorithm for classical density functional theory in three dimensions: Ionic solutions. J. Chem. Phys. 2010, 132, 124101. (34) Marrink, S. J.; Risselada, H. J.; Yefimov, S.; Tieleman, D. P.; de Vries, A. H. The MARTINI Force Field: Coarse Grained Model for Biomolecular Simulations. J. Phys. Chem. B 2007, 111, 7812−7824. 24


Page 25 of 27

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


(35) Panizon, E.; Bochicchio, D.; Monticelli, L.; Rossi, G. MARTINI Coarse-Grained Models of Polyethylene and Polypropylene. J. Phys. Chem. B 2015, 119, 8209−8216. (36) Rosenfeld, Y. Free-energy model for the inhomogeneous hard-sphere fluid mixture and density-functional theory of freezing. Phys. Rev. Lett. 1989, 63, 980−983. (37) Yu, Y.-X.; Wu, J. Structures of hard-sphere fluids from a modified fundamental-measure theory. J. Chem.

Phys. 2002, 117, 10156. (38) Soon-Chul, K.; Song Hi, L. A density functional perturbative approach for simple fluids: the structure of a nonuniform Lennard-Jones fluid at interfaces. J. Phys.: Condens. Matter 2004, 16, 6365. (39) Tang, Y.; Tong, Z.; Lu, B. C. Y. Analytical equation of state based on the Ornstein-Zernike equation. Fluid

Phase Equilib. 1997, 134, 21−42. (40) Tripathi, S.; Chapman, W. G. Microstructure and thermodynamics of inhomogeneous polymer blends and solutions. Phys. Rev. Lett. 2005, 94, 087801. (41) Jain, S.; Dominik, A.; Chapman, W. G. Modified interfacial statistical associating fluid theory: A perturbation density functional theory for inhomogeneous complex fluids. J. Chem. Phys. 2007, 127, 244904. (42) Jain, S.; Jog, P.; Weinhold, J.; Srivastava, R.; Chapman, W. G. Modified interfacial statistical associating fluid theory: Application to tethered polymer chains. J. Chem. Phys. 2008, 128, 154910. (43) Bymaster, A.; Chapman, W. G. An i SAFT density functional theory for associating polyatomic molecules.

J. Phys. Chem. B 2010, 114, 12298−12307. (44)Dijkstra, D. J.; Hoogsteen, W.; Pennings, A. J., Cross-linking of ultra-high molecular weight polyethylene in the melt by means of electron beam irradiation. Polymer 1989, 30, 866−873. 25



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 27

(45) Curro, J. G.; Frischknecht, A. L. The structure of poly(ethylene oxide) liquids: comparison of integral equation theory with molecular dynamics simulations and neutron scattering. Polymer 2005, 46, 6500−6506. (46) Martyna, G. J.; Tuckerman, M. E.; Tobias, D. J.; Klein, M. L. Explicit reversible integrators for extended systems dynamics. Mol. Phys. 1996, 87, 1117−1157. (47) Ryckaert, J.-P.; Ciccotti, G.; Berendsen, H. J. C. Numerical integration of the cartesian equations of motion of a system with constraints: molecular dynamics of n-alkanes. J. Comput. Phys. 1977, 23, 327−341. (48)McNally, T.; McShane, P.; Nally, G. M.; Murphy, W. R.; Cook, M.; Miller, A., Rheology, phase morphology, mechanical, impact and thermal properties of polypropylene/metallocene catalysed ethylene 1-octene copolymer blends. Polymer 2002, 43, 3785−3793. (49)Da Silva, A. L. N.; Rocha, M. C. G.; Coutinho, F. M. B.; Bretas, R.; Scuracchio, C., Rheological, mechanical, thermal, and morphological properties of polypropylene/ethylene-octene copolymer blends. J. Appl.

Polym. Sci. 2000, 75, 692−704. (50)Tian, M.; Li, T.; Zhang, L.; Tian, H.; Wu, Y.; Ning, N., Interfacial crystallization and its mechanism in in-situ dynamically vulcanized iPP/POE blends. Polymer 2014, 55, 3068−3074.

26


Page 27 of 27


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

For Table of Contents use only

Theoretical and Experimental Insights into the Phase Transition of Rubber/Plastic Blends during Dynamic Vulcanization Zhaoyang Wei, Shangqing Li, Nanying Ning, Ming Tian, Liqun Zhang, and Jianguo Mi

27


Theoretical and Experimental Insights into the Phase Transition of

Recommend Documents