环状聚合物及其对应的线性链熔体在启动剪切场下流变特性的分子动力学模拟研究
|
摘要
采用非平衡分子动力学方法研究了环状聚合物及其对应的线性链熔体在启动剪切场下的结构与流变特性.模拟结果显示:低剪切速率下(■<1×10~(-4)τ~(-1))环状链分子体系相比于同分子量的线性链体系并没有出现明显的过冲现象.该结果表明,在启动剪切场下环状分子与其对应的线性链相比较表现出了更弱的分子形变,同时模拟结果也与最近实验观察的结果一致.为了进一步探究这种现象背后的分子机理,在分子层面统计了不同流场强度下,链段的长度和取向角分布随着应变的演化,统计结果证明了环状聚合物分子链段弱的形变是导致其弱的剪切变稀和峰值应变的主要因素.本文还给出了过冲点和稳态在不同剪切速率下环状分子与其对应的线性链的流变特性(过冲应变、最大应力、最大黏性和稳态黏性)、结构和取向参数与维森伯格数(Wi_R)所满足的标度关系.
We analyze the structure and rheological properties of ring and linear polymers under shear byusing the non-equilibrium molecular dynamics simulation. The simulation results show that compared with the ring chains, the linear polymers do not present prominent stress over shoot phenomenon. Since the overshoot reflects the maximum flow-induced deformation of the polymer, this qualitative observation already implies that the ring experience less deformation than its linear precursor in simple shear flow. This is consistent with the recent experimental result. In order to further study the molecular mechanism of this phenomenon, the segmental structure and orientation angle distribution as a function of strain under the different Weissenberg numbers are given in this study. The weak overshoot of the stretching of the ring polymers proves that the weak shear thinning and peak strain are due to the weak deformation of the segment chain of the ring in the shear flow.The rheological properties of linear and ring system are extracted from the stress-strain curves, can be used further to analyze the data. The peak strain gmax as afunction of Wi_R follows a power-law with an exponent of 0.3 for linear polymer at Wi_R>1, however, for the ring system thepeak strain follows a power-law with an exponent of 0.1. The parameter h_(max)/h_(steady) is also the measure of the effective chain deformation at a steady state. The data show its progressive increase with Wi_R increasing, and follows a power-law with a scaling slop of 0.13 and 0.08 for linear and ring polymers, respectively. The peak stress smax as a function of Wi_R is also extracted from stress-strain curve. The two investigated systems both obey the scaling law with an exponent of 0.5. The normalized steady-state shear viscosity obeys a shear thinning slop of –0.86 for the linear polymer, the ring polymer obeysa shear thinning slop of –0.4. According to the gyration tensor and orientation angle, the power-law relationship between stretching and orientation is also given in this work.
We analyze the structure and rheological properties of ring and linear polymers under shear byusing the non-equilibrium molecular dynamics simulation. The simulation results show that compared with the ring chains, the linear polymers do not present prominent stress over shoot phenomenon. Since the overshoot reflects the maximum flow-induced deformation of the polymer, this qualitative observation already implies that the ring experience less deformation than its linear precursor in simple shear flow. This is consistent with the recent experimental result. In order to further study the molecular mechanism of this phenomenon, the segmental structure and orientation angle distribution as a function of strain under the different Weissenberg numbers are given in this study. The weak overshoot of the stretching of the ring polymers proves that the weak shear thinning and peak strain are due to the weak deformation of the segment chain of the ring in the shear flow.The rheological properties of linear and ring system are extracted from the stress-strain curves, can be used further to analyze the data. The peak strain gmax as afunction of Wi_R follows a power-law with an exponent of 0.3 for linear polymer at Wi_R>1, however, for the ring system thepeak strain follows a power-law with an exponent of 0.1. The parameter h_(max)/h_(steady) is also the measure of the effective chain deformation at a steady state. The data show its progressive increase with Wi_R increasing, and follows a power-law with a scaling slop of 0.13 and 0.08 for linear and ring polymers, respectively. The peak stress smax as a function of Wi_R is also extracted from stress-strain curve. The two investigated systems both obey the scaling law with an exponent of 0.5. The normalized steady-state shear viscosity obeys a shear thinning slop of –0.86 for the linear polymer, the ring polymer obeysa shear thinning slop of –0.4. According to the gyration tensor and orientation angle, the power-law relationship between stretching and orientation is also given in this work.
引文
[1]Li Y, Hsiao K W, Brockman C A, Yates D Y, RobertsonAnderson R M, Kornfield J A, San Francisco M J, Schroeder C M, Mc Kenna G B 2015 Macromolecules 48 5997
[2]Sanchez T, Kulic I M, Dogic Z 2010 Phys. Rev. Lett. 104098103
[3]Micheletti C, Orlandini E 2012 Macromolecules 45 2113
[4]Kapnistos M, Lang M, Vlassopoulos D, Pyckhout-Hintzen W,Richter D, Cho D, Chang T, Rubinstein M 2008 Nature Mater. 7 997
[5]Klein J 1986 Macromolecules 19 105
[6]Halverson J D, Grest G S, Grosberg A Y, Kremer K 2012Phys. Rev. Lett. 108 038301
[7]Hernández Cifre J G, Pamies R, López Martínez M C, García de la Torre J 2005 Polymer 46 267
[8]Yan Z C, Costanzo S, Jeong Y, Chang T, Vlassopoulos D2016 Macromolecules 49 1444
[9]Chen W, Zhang K, Liu L, Chen J, Li Y, An L 2017Macromolecules 50 1236
[10]Chen W, Zhao H, Liu L, Chen J, Li Y, An L 2015 Soft Matter11 5265
[11]Rubinstein M 1986 Phy. Rev. Lett. 57 3023
[12]Nechaev S K, Semenov A N, Koleva M K 1987 Physica A:Statist. Mech. Appl. 140 506
[13]Chen W, Li Y, Zhao H, Liu L, Chen J, An L 2015 Polymer 6493
[14]Cao J, Likhtman A E 2015 ACS Macro. Lett. 4 1376
[15]Luo C, Kr?ger M, Sommer J U 2017 Polymer 109 71
[16]Kremer K, Grest G S 1990 J. Chem. Phys. 92 5057
[17]Wang Z, Likhtman A E, Larson R G 2012 Macromolecules 453557
[18]Mc Leish T C B 2002 Adv. Phys. 51 1379
[19]Cao J, Likhtman A E 2012 Phy. Rev. Lett. 108 028302
[20]Doi Y, Matsubara K, Ohta Y, Nakano T, Kawaguchi D,Takahashi Y, Takano A, Matsushita Y 2015 Macromolecules48 3140
[21]Boukany P E, Wang S Q, Wang X 2009 J. Rheol. 53 617
[22]Jeong S, Kim J M, Baig C 2017 Macromolecules 50 3424
[2]Sanchez T, Kulic I M, Dogic Z 2010 Phys. Rev. Lett. 104098103
[3]Micheletti C, Orlandini E 2012 Macromolecules 45 2113
[4]Kapnistos M, Lang M, Vlassopoulos D, Pyckhout-Hintzen W,Richter D, Cho D, Chang T, Rubinstein M 2008 Nature Mater. 7 997
[5]Klein J 1986 Macromolecules 19 105
[6]Halverson J D, Grest G S, Grosberg A Y, Kremer K 2012Phys. Rev. Lett. 108 038301
[7]Hernández Cifre J G, Pamies R, López Martínez M C, García de la Torre J 2005 Polymer 46 267
[8]Yan Z C, Costanzo S, Jeong Y, Chang T, Vlassopoulos D2016 Macromolecules 49 1444
[9]Chen W, Zhang K, Liu L, Chen J, Li Y, An L 2017Macromolecules 50 1236
[10]Chen W, Zhao H, Liu L, Chen J, Li Y, An L 2015 Soft Matter11 5265
[11]Rubinstein M 1986 Phy. Rev. Lett. 57 3023
[12]Nechaev S K, Semenov A N, Koleva M K 1987 Physica A:Statist. Mech. Appl. 140 506
[13]Chen W, Li Y, Zhao H, Liu L, Chen J, An L 2015 Polymer 6493
[14]Cao J, Likhtman A E 2015 ACS Macro. Lett. 4 1376
[15]Luo C, Kr?ger M, Sommer J U 2017 Polymer 109 71
[16]Kremer K, Grest G S 1990 J. Chem. Phys. 92 5057
[17]Wang Z, Likhtman A E, Larson R G 2012 Macromolecules 453557
[18]Mc Leish T C B 2002 Adv. Phys. 51 1379
[19]Cao J, Likhtman A E 2012 Phy. Rev. Lett. 108 028302
[20]Doi Y, Matsubara K, Ohta Y, Nakano T, Kawaguchi D,Takahashi Y, Takano A, Matsushita Y 2015 Macromolecules48 3140
[21]Boukany P E, Wang S Q, Wang X 2009 J. Rheol. 53 617
[22]Jeong S, Kim J M, Baig C 2017 Macromolecules 50 3424