Effects of polymer concentration and chain length on aggregation in physically associating polymer solutions

The effects of polymer concentration and chain length on aggregation in associative polymer solutions, are studied using self-consistent e^ ld lattice model. Only two inhomogenous morphologies, i.e. microu^ ctuation homogenous (MFH) and micelle morphologies, are observed in the systems with differen...

Повний опис

Збережено в:
Бібліографічні деталі
Дата:2011
Автори: Han, X.-G., Zhang, X.-F., Ma, Y.-H., Zhang, C.-X., Guan, Y.-B.
Формат: Стаття
Мова:English
Опубліковано: Інститут фізики конденсованих систем НАН України 2011
Назва видання:Condensed Matter Physics
Онлайн доступ:https://nasplib.isofts.kiev.ua/handle/123456789/120038
Теги: Додати тег
Немає тегів, Будьте першим, хто поставить тег для цього запису!
Назва журналу:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Цитувати:Effects of polymer concentration and chain length on aggregation in physically associating polymer solutions / X.-G. Han, X.-F. Zhang, Y.-H. Ma, C.-X. Zhang, Y.-B. Guan // Condensed Matter Physics. — 2011. — Т. 14, № 4. — С. 43601:1-11. — Бібліогр.: 41 назв. — англ.

Репозитарії

Digital Library of Periodicals of National Academy of Sciences of Ukraine
id nasplib_isofts_kiev_ua-123456789-120038
record_format dspace
spelling nasplib_isofts_kiev_ua-123456789-1200382025-02-09T21:51:17Z Effects of polymer concentration and chain length on aggregation in physically associating polymer solutions Вплив концентрацiї полiмера i довжини ланцюга на агрегацiю у фiзично асоцiйованих розчинах полiмерiв Han, X.-G. Zhang, X.-F. Ma, Y.-H. Zhang, C.-X. Guan, Y.-B. The effects of polymer concentration and chain length on aggregation in associative polymer solutions, are studied using self-consistent e^ ld lattice model. Only two inhomogenous morphologies, i.e. microu^ ctuation homogenous (MFH) and micelle morphologies, are observed in the systems with different chain lengths. The temperatures at which the above two inhomogenous morphologies r^ st appear, which are denoted by TMFH and Tm, respectively, are dependent on polymer concentration and chain length. The variation of the logarithm of critical MFH concentration with the logarithm of chain length full^ s a linear-t^ ting relationship with a slope equaling −1. Furthermore, the variation of the average volume fraction of stickers at the micellar core (AVFSM) with polymer concentration and chain length is focused in the system at Tm. It is founded by calculations that the above behavior of AVFSM, is explained in terms of intrachain and interchain associations. Використовуючи ґраткову модель самоузгодженого поля, вивчається вплив концентрацiї полiмера i довжини ланцюга на агрегацiю в асоцiативних полiмерних розчинах. У системах з рiзною довжиною ланцюга спостерiгається тiльки двi неоднорiднi морфологiї, а саме мiкрофлуктуацiйна однорiдна (MFH) i мiцелярна морфологiї. Температури, при яких вище згаданi двi морфологiї виникають впер-ше i позначаються як TMFH i Tm, вiдповiдно, є незалежними вiд концентрацiї полiмера i довжини ланцюга. Змiна логарифма критичної концентрацiї MFH зi змiною логарифма довжини ланцюга задовiльняє спiввiдношенню лiнiйного допасування з нахилом рiвним −1. Крiм того, змiна середньої об’ємної фракцiї стiкерiв при мiцелярному корi (AVFSM) з концентрацiєю полiмера i довжиною ланцюга сфокусована в системi при Tm. Знайдено шляхом розрахункiв, що вище згадана поведiнка AVFSM пояснюється в термiнах iнтраланцюгових та iнтерланцюгових асоцiацiй. This research is supported by the Innovation Fund of Inner Mongolia University of Science and Technology (Grant No. 2010NC065) the High Performance Computers of Inner Mongolia University of Science and Technology. 2011 Article Effects of polymer concentration and chain length on aggregation in physically associating polymer solutions / X.-G. Han, X.-F. Zhang, Y.-H. Ma, C.-X. Zhang, Y.-B. Guan // Condensed Matter Physics. — 2011. — Т. 14, № 4. — С. 43601:1-11. — Бібліогр.: 41 назв. — англ. 1607-324X PACS: 61.25.Hp, 87.15.Nr, 82.60.Fa DOI:10.5488/CMP.14.43601 arXiv:1202.4595 https://nasplib.isofts.kiev.ua/handle/123456789/120038 en Condensed Matter Physics application/pdf Інститут фізики конденсованих систем НАН України
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection DSpace DC
language English
description The effects of polymer concentration and chain length on aggregation in associative polymer solutions, are studied using self-consistent e^ ld lattice model. Only two inhomogenous morphologies, i.e. microu^ ctuation homogenous (MFH) and micelle morphologies, are observed in the systems with different chain lengths. The temperatures at which the above two inhomogenous morphologies r^ st appear, which are denoted by TMFH and Tm, respectively, are dependent on polymer concentration and chain length. The variation of the logarithm of critical MFH concentration with the logarithm of chain length full^ s a linear-t^ ting relationship with a slope equaling −1. Furthermore, the variation of the average volume fraction of stickers at the micellar core (AVFSM) with polymer concentration and chain length is focused in the system at Tm. It is founded by calculations that the above behavior of AVFSM, is explained in terms of intrachain and interchain associations.
format Article
author Han, X.-G.
Zhang, X.-F.
Ma, Y.-H.
Zhang, C.-X.
Guan, Y.-B.
spellingShingle Han, X.-G.
Zhang, X.-F.
Ma, Y.-H.
Zhang, C.-X.
Guan, Y.-B.
Effects of polymer concentration and chain length on aggregation in physically associating polymer solutions
Condensed Matter Physics
author_facet Han, X.-G.
Zhang, X.-F.
Ma, Y.-H.
Zhang, C.-X.
Guan, Y.-B.
author_sort Han, X.-G.
title Effects of polymer concentration and chain length on aggregation in physically associating polymer solutions
title_short Effects of polymer concentration and chain length on aggregation in physically associating polymer solutions
title_full Effects of polymer concentration and chain length on aggregation in physically associating polymer solutions
title_fullStr Effects of polymer concentration and chain length on aggregation in physically associating polymer solutions
title_full_unstemmed Effects of polymer concentration and chain length on aggregation in physically associating polymer solutions
title_sort effects of polymer concentration and chain length on aggregation in physically associating polymer solutions
publisher Інститут фізики конденсованих систем НАН України
publishDate 2011
url https://nasplib.isofts.kiev.ua/handle/123456789/120038
citation_txt Effects of polymer concentration and chain length on aggregation in physically associating polymer solutions / X.-G. Han, X.-F. Zhang, Y.-H. Ma, C.-X. Zhang, Y.-B. Guan // Condensed Matter Physics. — 2011. — Т. 14, № 4. — С. 43601:1-11. — Бібліогр.: 41 назв. — англ.
series Condensed Matter Physics
work_keys_str_mv AT hanxg effectsofpolymerconcentrationandchainlengthonaggregationinphysicallyassociatingpolymersolutions
AT zhangxf effectsofpolymerconcentrationandchainlengthonaggregationinphysicallyassociatingpolymersolutions
AT mayh effectsofpolymerconcentrationandchainlengthonaggregationinphysicallyassociatingpolymersolutions
AT zhangcx effectsofpolymerconcentrationandchainlengthonaggregationinphysicallyassociatingpolymersolutions
AT guanyb effectsofpolymerconcentrationandchainlengthonaggregationinphysicallyassociatingpolymersolutions
AT hanxg vplivkoncentraciípolimeraidovžinilancûganaagregaciûufizičnoasociiovanihrozčinahpolimeriv
AT zhangxf vplivkoncentraciípolimeraidovžinilancûganaagregaciûufizičnoasociiovanihrozčinahpolimeriv
AT mayh vplivkoncentraciípolimeraidovžinilancûganaagregaciûufizičnoasociiovanihrozčinahpolimeriv
AT zhangcx vplivkoncentraciípolimeraidovžinilancûganaagregaciûufizičnoasociiovanihrozčinahpolimeriv
AT guanyb vplivkoncentraciípolimeraidovžinilancûganaagregaciûufizičnoasociiovanihrozčinahpolimeriv
first_indexed 2025-12-01T04:04:57Z
last_indexed 2025-12-01T04:04:57Z
_version_ 1850277267956039680
fulltext Condensed Matter Physics, 2011, Vol. 14, No 4, 43601: 1–11 DOI: 10.5488/CMP.14.43601 http://www.icmp.lviv.ua/journal Effects of polymer concentration and chain length on aggregation in physically associating polymer solutions X.-G. Han1∗, X.-F. Zhang1, Y.-H. Ma1, C.-X. Zhang2, Y.-B. Guan2 1 Inner Mongolia Key Laboratory for Utilization of Bayan Obo Multi-Metallic Resources: Elected State Key Laboratory, School of Mathematics, Physics and Biological engineering, Inner Mongolia University of Science and Technology, 014010 Baotou, China 2 Department of Physics, Jilin University, 130021 Changchun, China Received July 2, 2011, in final form September 19, 2011 The effects of polymer concentration and chain length on aggregation in associative polymer solutions, are studied using self-consistent field lattice model. Only two inhomogenous morphologies, i.e. microfluctuation homogenous (MFH) and micelle morphologies, are observed in the systems with different chain lengths. The temperatures at which the above two inhomogenous morphologies first appear, which are denoted by TMFH and Tm, respectively, are dependent on polymer concentration and chain length. The variation of the logarithm of critical MFH concentration with the logarithm of chain length fulfils a linear-fitting relationship with a slope equaling −1. Furthermore, the variation of the average volume fraction of stickers at the micellar core (AVFSM) with polymer concentration and chain length is focused in the system at Tm. It is founded by calculations that the above behavior of AVFSM, is explained in terms of intrachain and interchain associations. Key words: concentration, chain length, aggregation, associative polymer PACS: 61.25.Hp, 87.15.Nr, 82.60.Fa 1. Introduction Physically, associating polymers are polymer chains containing a small fraction of attractive groups along the backbones. The attractive groups, for example, solvophobic groups, tend to form physical links which can play an important role in reversible junctions between different polymer chains. The junctions can be broken and recombined frequently on experimental time scales. This property of junctions makes associative polymer solutions behave reversibly when ambient condi- tions, such as temperature and concentrations, change. This tunable characteristic of the system produces extensive applications [1–4] that possess great potential as smart materials [5–7]. The attractive groups (also called sticker monomers) drive the self-assembly of the polymers, leading to the formation of polymeric micelles in physically associating polymer solutions (PAPSs). In telechelic [8, 9] and multiblock [10] associative polymers, flower micelles were observed. This aggregation, as well as their ability to form bridges between micelles, profoundly affect their macro- scopic properties, particularly their rheological behavior, which is required in many applications. In telechelic associative polymers, the effects of architectural parameters of polymers have been assessed, such as chain length, end-group length, and chemical composition [11–15]. In multi- block associative polymers, the effect of chain architecture of polymer on the property of aggre- gates [16, 17] was studied, which suggested that chain architecture can be an important factor in controlling macroscopic properties of the systems. However, the above studies of multiblock PAPSs were carried out in two dimensions and at low concentrations. It is well known that self-consistent field theory (SCFT), as a mean-field theory, has been applied to the study of a great deal of problems in polymeric systems [18–21]. Recently, SCFT has ∗E-mail: xghan0@163.com c© X.-G. Han, X.-F. Zhang, Y.-H. Ma, C.-X. Zhang, Y.-B. Guan, 2011 43601-1 http://dx.doi.org/10.5488/CMP.14.43601 http://www.icmp.lviv.ua/journal X.-G. Han et al. been applied to the study of the properties of micelles in polymer solutions [22–24]. In the previous paper [10], given a fixed chain length, we focused on the thermodynamic properties and structure transitions in PAPSs. The microfluctuation homogenous (MFH) morphology, which corresponds to the onset of gelation [10, 25], and micelle morphology were observed. The degrees of aggregations of the above two morphologies are very different. The volume fraction of stickers in micelle morphology is much bigger than that in MFH morphology. In this work, the property of the aggregation in PAPAs is studied using self-consistent field lattice model. The temperatures at which MFH and micelle morphologies first appear, which are denoted by TMFH and Tm, respectively, and the average volume fraction of stickers at the micellar core, 〈φs(rco)〉, are calculated. Such calculations are carried out for different chain lengths and polymer concentrations. It is found that TMFH, Tm and 〈φ(rco)〉 are dependent on chain length and polymer concentration, and the relationship of 〈φs(rco)〉 with chain length and polymer concentration is explained in terms of intrachain and interchain associations. 2. Theory We consider a system of incompressible PAPSs, where nP polymers, each of which is composed of Nst segments of sticker monomer type (attractive group) and Nns segments of nonsticky monomer type, are distributed over a lattice. Each sticker monomer is a regularly placed apart l monomer along the chain backbone. The degree of polymerization of chain is N = Nst +Nns. In addition to polymer monomers, nh solvent molecules are placed on the vacant lattice sites. Sticker, nonsticky monomers and solvent molecules are of the same size and each occupies one lattice site. The total number of lattice sites is NL = nh + nPN . Nearest neighbor pairs of stickers have attractive interaction −ǫ with ǫ > 0, which is the only non-bonded interaction in the present system. The interaction energy is expressed as: U = − ǫ 2 ∑ r ∑ r′ φ̂st(r)φ̂st(r ′), (2.1) where ∑ r means the summation over all the lattice sites r and ∑ r′ means the summation over the nearest neighbor sites of r. φ̂st(r) = ∑ j ∑ s∈st δr,rj,s is the volume fraction of stickers on site r, where j and s are the indexes of chain and monomer of a polymer, respectively. s ∈ st means that the sth monomer belongs to sticker monomer type. In this simulation, however, instead of directly using the exact expression of the nearest neighbor interaction for stickers, we introduce a local concentration approximation for the non-bonded interaction [10, 26]. ∑ r′ φ̂st(r ′) in equation (2.1) is replaced with z φ̂st(r), where z is the coordination number of the lattice used. Within this approximation, the interaction energy is expressed as: U kBT = −χ ∑ r φ̂st(r)φ̂st(r), (2.2) where χ is the Flory-Huggins interaction parameter in the solutions, which is equal to z 2kBT ǫ. We perform the SCFT calculations in the canonical ensemble, and the field-theoretic free energy F is defined as follows: F [ω+, ω−] kBT = ∑ r { 1 4χ ω2 −(r)− ω+(r) } − nP lnQP[ωst, ωns]− nh lnQh[ωh], (2.3) where Qh is the partition function of a solvent molecule subject to the field ωh(r) = ω + (r), which is defined as Qh = 1 nh ∑ r exp [− ωh(r)]. QP is the partition function of a noninteraction polymer chain subject to the fields ωst(r) = ω + (r) − ω − (r) and ωns(r) = ω + (r), which act on sticker and nonsticky segments, respectively. Following the scheme by Schentiens and Leermakers [27], QP is expressed as QP = 1 NL 1 z ∑ rN ∑ αN GαN(r,N |1), where rN and αN denote the position and orientation of the Nth segment of the chain, respectively. ∑ rN ∑ αN means the summation over 43601-2 Effects of polymer concentration and chain length on aggregation all the possible positions and orientations of the Nth segment of the chain. Gαs(r, s|1) is the end segment distribution function of the sth segment of the chain, which is evaluated from the following recursive relation: Gαs(r, s|1) = G(r, s) ∑ r′ s−1 ∑ αs−1 λ αs−αs−1 rs−r′ s−1 Gαs−1(r′, s− 1|1), (2.4) where G(r, s) is the free segment weighting factor and is expressed as G(r, s) = { exp[−ωst(rs )], s ∈ st ; exp[−ωns(rs )], s ∈ ns . The initial condition is Gα1(r, 1|1) = G(r, 1) for all the values of α1. In the above expression, the values of λ αs−αs−1 rs−r′s−1 depend on the chain model used. We assume that λ αs−αs−1 rs−r′s−1 = { 0, αs = αs−1 ; 1/(z − 1), otherwise . This means that the chain is described as a random walk without the possibility of direct back- folding. Although self-intersections of a chain are not permitted, the excluded volume effect is sufficiently taken into account [28]. Another end segment distribution function Gαs(r, s|N) is eval- uated from the following recursive relation: Gαs(r, s|N) = G(r, s) ∑ r′ s+1 ∑ αs+1 λ αs+1−αs r′ s+1 −rs Gαs+1(r′, s+ 1|N), (2.5) with the initial condition GαN(r,N |N) = G(r,N) for all the values of αN. Using the expressions of the end segment distribution functions, the single-segment probability distribution function P (1)(r, s) and the two-segment probability distribution function P (2)(r1, s1; r2, s2) of the chain can be defined as follows: P (1)(r, s) = 1 zNLQP ∑ r′s ∑ αs Gαs(r′, s|1)Gαs(r′, s|N) G(r′, s) δr′s,r , (2.6) which is the normalized probability that the monomer s of the chain is on the lattice site r; P (2)(r1, s1; r2, s2) = 1 zNLQP ∑ r ′ s1 ∑ αs1 ∑ r ′ s2 ∑ αs2 Gαs1 (r ′ , s1|1)δr′s1 ,r1 × G(r ′ , s1; r ′ , s2)G αs2 (r ′ , s2|N)δr′s2 ,r2 (2.7) and G(r, s1; r, s2) = ∑ rs1+1 ∑ αs1+1 ..... ∑ rs2−1 ∑ αs2−1 { s2−1∏ s=s1+1 λ αs−αs−1 rs−rs−1 G(r, s) } λ αs2 −αs2−1 rs2−rs2−1 (for s2 > s1) yield the probability that the monomers s1 and s2 of the chain are on the lattice sites r1 and r2, respectively. It can be verified that ∑ r P (1)(r, s) = 1, and ∑ r2 P (2)(r1, s1; r2, s2) = P (1)(r1, s1). Equation (2.3) can be considered as the alternative form of the self-consistent field free energy functional for incompressible polymer solutions [29]. When a local concentration approximation for the non-bonded interaction is introduced, the SCFT description of lattice model for PAPSs presented in this work is basically equivalent to that of the “Gaussian thread model” chain for the similar polymer solutions [29]. The related details are presented in [10]. Minimization of the free energy function F with ω − (r) and ω + (r) leads to the following saddle point equations: ω − (r) = 2χφst(r), (2.8) 43601-3 X.-G. Han et al. φst(r) + φns(r) + φh(r) = 1, (2.9) where φst(r) = 1 NL 1 z nP QP ∑ s∈st ∑ αs Gαs(r, s|1)Gαs(r, s|N) G(r, s) (2.10) and φns(r) = 1 NL 1 z nP QP ∑ s∈ns ∑ αs Gαs(r, s|1)Gαs(r, s|N) G(r, s) (2.11) are the average numbers of sticker and nonsticky segments at r, respectively, and φh(r) = 1 NL nh Q h exp [−ωh(r)] is the average numbers of solvent molecules at r. The saddle point is calculated using the pseudo-dynamical evolution process presented by Fredrickson et al. [30]: ωnew + (r) = ωold + (r) + λ+(φst(r) + φns(r) + φh(r)− 1), (2.12) ωnew − (r) = ωold − (r) + λ−(φst(r) − ω − (r) 2χ ). (2.13) The calculation is initiated from appropriately randomly-chosen fields ω + (r) and ω−(r), and stopped when the change of free energy F between two successive iterations is reduced to the needed precision. The resulting configuration is taken as a saddle point configuration. By comparing the free energies of the saddle point configurations obtained from different initial fields, the relative stability of the observed morphologies can be assessed. 3. Result and discussion In our studies, the properties of associating polymers depend on three tunable molecular pa- rameters: χ (The Flory-Huggins interaction parameter), N (Chain length) and l (The number of nonsticky monomers between two neighboring stickers along the backbone, l equals 9 in this paper). The calculations are performed in three-dimensional simple cubic lattice with periodic boundary condition. The aggregation behavior is first sketched using the lattice with the size NL = 263. Then, the obtained results are verified using larger size lattices. The results presented below are obtained from the lattice with NL = 403. Three different morphologies, i.e., the homogenous, micro-fluctuation homogenous and micelle morphologies, are observed in PAPSs. By comparing the relative stability of the observed states, the phase diagram is constructed. Figure 1 shows the phase diagram of the systems with different chain lengths N . In this study, when N is changed, only MFH morphology and micelles are observed as inhomogeneous morpholo- gies. The structural morphology of MFH morphology does not change, and micellar shape remains spherelike (see figure 2). For N = 21, the χ value on micellar boundary (∼ 1/Tm) increases with decreasing φ̄P. Approaching to critical micelle concentration (φ̄CMC = 0.04)1, micellar boundary abruptly becomes steep. The χ value on MFH boundary (∼ 1/TMFH) also goes up with the decrease in φ̄P, which resembles the behavior on the micellar boundary. The critical MFH concentration (φ̄CFC = 0.37) is much higher than that of micellar morphology. Therefore, the MFH boundary intersects with the micellar boundary at φ̄CFC. When N is increased, at fixed φ̄P, the χ value on micellar boundary shifts slightly to larger value, and the χ value on MFH boundary decreases markedly, which is different from that on micellar boundary. φ̄CFC also decreases with the increase in N . Figure 3 shows the variations of the logarithm of φ̄CFC with the logarithm of N . It is seen that the straight line with a slope equaling −1 fits the results for all chain lengths considered in this study. The appearance of MFH morphology is considered as the onset of gelation [10, 25]. Therefore, 1Critical micelle concentration (CMC) is generally considered as the minimum concentration of the micellar appearance. In this study, the concentration which corresponds to the χ abrupt increase on micellar boundary is regarded as CMC, which is denoted by φ̄CMC. When N = 21, φ̄CMC = 0.04. When N is increased φ̄CMC decreases. For N = 101, the φ̄CMC is not observed till φ̄P = 0.015. 43601-4 Effects of polymer concentration and chain length on aggregation 0.00 0.15 0.30 0.45 0.60 0.75 0.90 0 2 4 6 8 P Figure 1. (Color online) The phase diagram for the systems with different chain lengths N . The boundaries of MFH and micelle morphologies are obtained. The red open and solid squares, green open and solid triangles, blue open and solid diamonds, magenta open and solid pen- tagons correspond to the boundaries of MFH and micelle morphologies for N = 21, 41, 81, 101, respectively. 0 1 2 0.00 0.02 0.04 0.06 0.08 0.8 0.9 1.0 1.1 0.0 0.2 0.4 0.6 0.8 1.0 ( ) st r r ns Figure 2. (Color online) The variations of the volume fractions of sticker and nonsticker compo- nents of polymers in micelle morphology with r when φ̄P = 0.1 and χ = 5.5 in the systems with N = 21 and 101. r is the distance from the sticker-rich core. The red solid and open squares, blue solid and open triangles denote the volume fractions of stickers and nonstickers for N = 101 and 21, respectively. critical MFH concentration φ̄CFC should correspond to the critical gelation concentration C∗, which is analyzed in terms of the concept of chain overlap. The critical gelation concentration C∗, at which quasi-ideal coils begin to overlap, the pervaded volume of one another is related to the chain length as N = KC−2 ∗ , where K is a constant [31–33]. The corresponding slope of the critical concentration C∗ at quasi-ideal coil is −2. Similarly, at the excluded volume chain, the corresponding quantity is −5/4. The slope of fitting straight line at associative polymer chain is the smallest of the above three cases. It is demonstrated that the physically associating polymer chain in solution should be elongated compared with the excluded volume chain and quasi-ideal coil. In order to measure the degree of aggregation of stickers, an order-parameter-type variable [10] is used in this work. It is expressed as: Φ({φst(r)}) = 1 NL ∑ r (φst(r)− φ̄st) 2 = 1 NL ∑ r φ2 st(r) − φ̄2 st , (3.1) 43601-5 X.-G. Han et al. 1.20 1.35 1.50 1.65 1.80 1.95 2.10 -1.20 -1.05 -0.90 -0.75 -0.60 -0.45 -0.30 log CFC log N Figure 3. The logarithm of critical MFH concentration φ̄CFC a function of the logarithm of chain length N , which fulfils a linear-fitting relationship with the slope equaling -1. 0.00 0.75 1.50 2.25 3.00 3.75 4.50 0.000 0.015 0.030 0.045 0.060 0.075 Figure 4. (Color online) The variation of the order-parameter-type variable Φ with χ near the HS-MFH and MFH-micelle transition points at φ̄P = 0.8 in the systems with different chain length N . The red squares and blue triangles correspond to N = 21 and N = 81, respectively. With the increase in χ, there are the appearances of two nonzero steps on the above two curves, which denote the HS-MFH and MFH-micelle transitions, respectively. where Φ is determined by the distribution of the volume fraction of stickers {φst(r)}, which is a function of χ and φ̄P. For homogenous solutions (HS), Φ is equal to zero. Figure 4 shows the variations of Φ in MFH and micelle morphologies with an increasing χ when φ̄P = 0.8 near HS- MFH and MFH-micelle transition points, in the systems with N = 21 and N = 81. It is seen that the degree of the aggregation of stickers increases with an increasing χ at fixed N . Compared with N = 81, although the value of Φ at fixed χ above MFH-micelle transition point at N = 21 is larger, the variations of Φ on the above two cases with χ are similar. However, near HS-MFH transition point, the variation of Φ with χ at N = 21 is obviously different from that of N = 81. The value of Φ at N = 21 increases more rapidly than that of N = 81 when χ is increased. The heat capacity CV is proportional to the first derivative of Φ with respect to temperature [10, 34, 35]. When N is increased, the maximum of the peak of CV decreases, and the corresponding half-width increases near the HS-MFH transition. The shape of the peak of CV near MFH-micelle transition point does not practically change with the increase in N (not shown). It is demonstrated that the change of N has a greater effect on HS-MFH transition than that on MFH-micelle transition. The aggregation number of micelles, which is used to account for some properties, for example, 43601-6 Effects of polymer concentration and chain length on aggregation 0.00 0.15 0.30 0.45 0.60 0.75 0.90 0.75 0.80 0.85 0.90 0.95 1.00 S co r P Figure 5. (Color online) The average volume fractions of stickers at micellar core as a function of φ̄P in micellar boundary systems with different chain lengths. The red squares, green triangles, blue diamonds and magenta pentagons correspond to the boundaries of micelle morphology for N = 21, 41, 81, 101, respectively. the rheological behavior of associative polymers, is a major focus in many experimental [36–38] and theoretic [39] studies. In this paper, the average volume fraction of stickers at micellar cores 〈φs(rco)〉 is similar to the average aggregation number of the micelle. Figure 5 shows the variations of 〈φs(rco)〉 on micellar boundary with φ̄P and N in the systems. Being given a fixed chain length, when φ̄P is increased, 〈φs(rco)〉 decreases. When N is increased, at fixed φ̄P, the variation of 〈φs(rco)〉 with φ̄P, does not practically change for φ̄P 6 0.4, on the other hand, when 0.4 < φ̄P 6 0.9, 〈φs(rco)〉 increases. The decreasing tendency of 〈φs(rco)〉 with φ̄P becomes gentle with the increase in N at high polymer concentrations. It is shown that the effect of the increase in N on 〈φs(rco)〉 is related to φ̄P. In order to interpret the effects of N and φ̄P on the aggregation of stickers in micelle morphology, we evaluate the probabilities that a polymer chain forms intrachain and interchain associations in the system using an approach similar to the one presented in Refs. [10, 40, 41]. We suppose that there are no other sticker aggregates in the micellar system except the micellar cores because the volume fraction of stickers at micellar core approaches 1, and that of the rest of the system is less than 10−1. A sticker in a particular chain can form an intrachain association, as well as interchain association. Ignoring the probabilities that more than two stickers of a definite chain are attached to a micellar core, the conditional probability that the sticker s1 concerns the intrachain association, provided that the sticker s1 is at the micelle core rco, can be expressed as: ploop(rco, s1) = 1 P (1)(rco, s1) ∑ s2∈st,s2 6=s1 P (2)(rco, s1; rco, s2), (3.2) where ∑ s2∈st,s2 6=s1 means the summation over all the stickers of a polymer chain except the s1th one, P (1)(rco, s1) and P (2)(rco, s1; rco, s2) are the single-segment and two-segment probability distribution functions of a chain, respectively. Then, 1−ploop(rco, s1) is the conditional probability that the sticker s1 is linked with those belonging to other chains when the sticker s1 is at rco, and Plk(s1) = ∑ rco P (1)(rco, s1) · ⌈1− ploop(rco, s1)⌉ is the probability that a sticker s1 of a chain is related to interchain association, where ∑ rco means the summation over all the micellar cores of the system. The summation of Plk(s1) over all the stickers in a chain, 〈n lk 〉 = ∑ s1,s1∈st Plk(s1), can be viewed as the average sticker num- ber from a particular polymer chain linked with other chains by sticker aggregates. The total number of stickers participating in interchain association in the system is expressed as ninte = 43601-7 X.-G. Han et al. nP〈nlk 〉. The total number of stickers belonging to intrachain association is defined as nintr = (1/Nst)[ ∑ s1,s1∈st{(1/Nm){ ∑ rco ploop(rco, s1)}}]·〈φs(rco)〉Nm, where Nm denotes the micellar num- ber in the system. 0.00 0.15 0.30 0.45 0.60 0.75 0.90 0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 0.6 0.8 1.0 int r fint e f P Figure 6. (Color online) The variations of the average fractions of intrachain and interchain associations at a micellar core with φ̄P in the same systems presented in figure 5, denoted by fintr and finte, respectively. The red open and solid squares, green open and solid triangles, blue open and solid diamonds, magenta open and solid pentagons correspond to fintr and finte for N = 21, 41, 81, 101, respectively. Figure 6 shows the variations of the average fractions of intrachain and interchain associations within a micelle, denoted by fintr = nintr/Nm and finte = ninte/Nm, respectively, with φ̄P and N in the systems identical with those presented in figure 5. Being given a fixed chain length, when φ̄P is increased, fintr does decrease markedly first, and then the decreasing tendency of fintr becomes gentle when φ̄P exceeds some value. Except for the case of N = 21, finte at fixed N increases markedly first, and then increases gently with the increase in φ̄P, which is contrary to fintr. At N = 21, finte increases to φ̄P = 0.6, then decreases slightly with the increase in φ̄P. When N is increased, fintr at fixed φ̄P increases, and corresponding finte decreases. When φ̄P is increased the differences of fintr and finte in various N become small. In other words, the effect of chain length on fintr and finte at low concentrations is more pronounced than that at high concentrations. The relative magnitude of fintr and finte at fixed N is dependent on φ̄P. Except for N = 21, when φ̄P is low, finte at fixed N is smaller than corresponding fintr. For N = 101, fintr equals 1, and corresponding finte equals =0 at φ̄P = 0.015 and χ = 6.0. The micelle is absolutely aggregated by intrachain association. When φ̄P exceeds some value, which is related to N , finte at fixed N is larger than corresponding fintr. When N = 21, finte is always larger than fintr in the considered range of φ̄P. In light of the variations of 〈φs(rco)〉 with φ̄P and N in the systems, the effect of chain length on aggregation of stickers is more pronounced at high concentrations. However, seen from respec- tive value of fintr and finte at fixed φ̄P, the increase of chain length is more favorable to the change of intrachain and interchain associations at low concentrations, which is different from the above conclusions drawn in light of 〈φs(rco)〉. It is shown that [see figure 6], being given a fixed chain length, although fintr is smaller than the corresponding finte at high concentrations, the decrease of 〈φs(rco)〉 with the increase in φ̄P is similar to that of fintr, instead of finte, in the considered range of φ̄P. When N is increased, fintr at fixed φ̄P goes up. It is notable that the effect of chain length on fintr at low concentrations is much more pronounced than that at high concentrations. However, the increase of fintr at fixed φ̄P does not result in the dramati- cal change of 〈φs(rco)〉 at low concentrations. At high concentrations (φ̄P > 0.4), on the con- trary, 〈φs(rco)〉 with a big chain length has a larger value. It is reasonable that 〈φs(rco)〉 is also affected by finte. When φ̄P is low, finte at fixed φ̄P with a long chain length is much smaller than that with a short chain. Therefore, the increase of fintr at fixed φ̄P is canceled by the de- 43601-8 Effects of polymer concentration and chain length on aggregation crease of corresponding finte, thus increasing in N . With the increase in φ̄P, finte at fixed φ̄P goes up markedly. The difference of finte among different N becomes small. The contribution of fintr to 〈φs(rco)〉 resulting from the increase in N begins to appear. It is demonstrated that the effects of chain length and polymer concentration on 〈φs(rco)〉 do interact each other. The in- crease of chain length is favorable to intrachain association, and retains interchain association, which is contrary to the effect of the increase in φ̄P. When φ̄P is increased, the effect of chain length on intrachain and interchain associations is weakened. In other words, the polymer con- centration and chain length simultaneously control the formations of intrachain and interchain associations. 4. Conclusion and summary Using the self-consistent field lattice model, the effects of polymer concentration φ̄P and chain length N on aggregation in physically associating polymer solutions are studied. When N is changed, only two inhomogenous aggregates, i.e., the MFH and micelle morphologies, are ob- served in PAPSs. When φ̄P is decreased, being given a fixed N , the χ values on MFH and micellar boundaries increase. At fixed φ̄P, the increase in N remarkably decreases the χ value on MFH boundary, but slightly increases the χ value on micellar boundary. The logarithm of critical MFH concentration as a function of the logarithm of N fulfils a fitting straight line with a slope equal- ing -1, which demonstrates that the associating polymer chain in solution should be elongated compared with an excluded volume chain. Furthermore, on micellar boundary, the average volume fraction of stickers at a micellar core, 〈φs(rco)〉, which is similar to the average aggregation number, decreases at fixed N when φ̄P is increased. With the increase in N , on the other hand, 〈φs(rco)〉, at fixed φ̄P, does not practically change when φ̄P 6 0.4. When φ̄P > 0.4, on the other hand, the increase in N causes the increase of 〈φs(rco)〉 at fixed φ̄P. There is found a decreasing tendency of variation of 〈φs(rco)〉 at fixed N , with φ̄P is determined by intrachain association, instead of interchain association, and the magnitude of 〈φs(rco)〉 with different N is affected by intrachain and interchain associations. At fixed φ̄P, when the difference of the contribution of interchain as- sociation to 〈φs(rco)〉 between different N goes down to some extent, the corresponding effect of intrachain association on 〈φs(rco)〉 begins to appear. The architectural parameters of polymer are important to the properties of PAPSs. In this paper, only the effect of chain length of the polymer is investigated. When the other parameters, for example, the nonsticky monomer number between two neighboring stickers l, are changed, the aggregation behavior is different. When l is increased, the χ values on MFH and micellar boundaries rise. The increase of l has a different effect on the specific heat peaks for HS-MFH and MFH-micelle transitions. Systemical studies are presented in our subsequent work. Our calculations are within the mean field framework, which may also be taken as a starting point to the further study, i.e., the direct calculation of the partition function of the system using complex Langevin simulations [29]. Acknowledgements This research is supported by the Innovation Fund of Inner Mongolia University of Science and Technology (Grant No. 2010NC065) the High Performance Computers of Inner Mongolia University of Science and Technology. References 1. Clark A.H., Ross-Murphy S.B., Structural and mechanical properties of biopolymer gels, In: Advances in Polymer Science, 1987, 83, 57; doi:10.1007/BFb0023332. 2. Slater G.W., Kist T.B.L., Ren H., Drouin G., Electrophoresis, 1998, 19, No. 10, 1525; doi:10.1002/elps.1150191003. 3. Tonge S.R., Tighe B.J., Adv. Drug Delivery Rev., 2001, 53, No. 1, 109; doi:10.1016/S0169-409X(01)00223-X. 43601-9 http://dx.doi.org/10.1007/BFb0023332 http://dx.doi.org/10.1002/elps.1150191003 http://dx.doi.org/10.1016/S0169-409X(01)00223-X X.-G. Han et al. 4. Taylor K.C., Nasr-El-Din H. A., J. Petrol. Sci. Eng., 1998, 19, No. 3–4, 265; doi:10.1016/S0920-4105(97)00048-X. 5. Brunsveld L., Folmer B.J.B., Meijer E.W., Sijbesma R.P., Chem. Rev., 2001, 101, No. 12, 4071; doi:10.1021/cr990125q. 6. De Greef T. F. A., Meijer E. W., Nature, 2008, 453, No. 7192, 171; doi:10.1038/453171a. 7. Cordier P., Tournilhac F., Soulié-Ziakovic C., Leibler L., Nature, 2008 451, 977; doi:10.1038/nature06669. 8. Alami E., Almgren M., Brown W., Macromolecules, 1996, 29, No. 6, 2229; doi:10.1021/ma951174h. 9. Chassenieux C., Nicolai T., Durand D., Macromolecules, 1997, 30, No. 17, 4952; doi:10.1021/ma970354j. 10. Han X.-G., Zhang C.-X., J. Chem. Phys., 2010, 132, No. 16, 164905; doi:10.1063/1.3400648. 11. Chassenieux C., Nicolai T., Durand D., Macromolecules, 1998, 31, No. 12, 4035; doi:10.1021/ma9802273. 12. Beaudoin E., Borisov O., Lapp A., Billon L., Hiorns R.C., Franc J., Macromolecules, 2002, 35, No. 19, 7436; doi:10.1021/ma011027l. 13. Beaudoin E., Gourierb C., Hiornsc R.C., Franc J., J. Colloid Interface Sci., 2002, 251, No. 2, 398; doi:10.1006/jcis.2002.8402. 14. Lafleche F, Durand D., Nicolai T., Macromolecules, 2003, 36, No. 4, 1331; doi:10.1021/ma021076d. 15. Lafleche F., Nicolai T., Durand D., Macromolecules, 2003, 36, No. 4, 1341; doi:10.1021/ma0210776. 16. Balazs A.C., Gempe M.C., Lentvorski A.P., J. Chem. Phys., 1991, 95, No. 11, 8467; doi:10.1063/1.461276. 17. Brown G., Chakrabarti A., J. Chem. Phys., 1992, 96, No. 10, 3251; doi:10.1063/1.461969. 18. Orland H., Schick M., Macromolecules, 1996, 29, No. 2, 713; doi:10.1021/ma9508461. 19. Matsen M.W., Schick M., Phys. Rev. Lett., 1994, 72, No. 16, 2660; doi:10.1103/PhysRevLett.72.2660. 20. Tang P., Qiu F., Zhang H., Yang Yu., Phys. Rev. E, 2004, 69, No. 3, 031803; doi:10.1103/PhysRevE.69.031803. 21. He X., Liang H., Huang L., Pan C., J. Phys. Chem. B, 2004, 108, No. 5, 1731; doi:10.1021/jp0359337. 22. Cavallo A., Muller M., Binder K., Macromolecules, 2006, 39, No. 26, 9539; doi:10.1021/ma061493g. 23. Jelinek K., Limpouchova Z., Uhlik F., Prochazka K., Macromolecules, 2007, 40, No.21, 7656; doi:10.1021/ma070928c. 24. Charlaganov M., Borisov O.V., Leermakers F.A.M., Macromolecules, 2008, 41, No. 10, 3668; doi:10.1021/ma800130q. 25. Kumar S. K., Doughlas J. F., Phys. Rev. Lett., 2001, 87, No. 5, 188301; doi:10.1103/PhysRevLett.87.188301. 26. Chen J.-Z., Zhang Ch.-X., Sun Zh.-Y., Zheng Y.-S., An L.-J., J. Chem. Phys., 2006, 124, No. 10, 104907; doi:10.1063/1.2176619. 27. Leermakers F.A.M., Scheutjens J.M.H.M., J. Chem. Phys., 1988, 89, No. 5, 3264; doi:10.1063/1.454931. 28. Medvedevskikh Y.G. Condens. Matter Phys., 2001, 4, No. 2, 209. 29. Fredrickson G.H., The Equilibrium Theory of Inhomogenous Polymers, Clarendon Press, Oxford, 2005; doi:10.1093/acprof:oso/9780198567295.001.0001. 30. Alexander-Katz A., Fredrickson G.H., Macromolecules, 2002, 35, 16; doi:10.1021/ma011515t. 31. Daoud M., Jannink G., J. Phys. (Paris), 1976, 37, 973; doi:10.1051/jphys:01976003707-8097300. 32. Daoud M., Cotton J.P., Farnoux B., Jannink G., Sarma G., Benoit H., Duplessix C., Picot C., de Gennes P.G., Macromolecules, 1975, 8, No. 6, 804; doi:10.1021/ma60048a024. 33. Klein J., Macromolecules, 1978, 11, No. 5, 852; doi:10.1021/ma60065a002. 34. Dudowicz J., Freed K. F., Douglas J. F., J. Chem. Phys., 1999, 111, No. 15, 7116; doi:10.1063/1.480004. 35. Kennedy J., Wheeler J. C., J. Chem. Phys., 1983, 78, No. 2, 953; doi:10.1063/1.444800. 36. Elliot P.T., Xing L., Wetzel W.H., Glass J.E., Macromolecules, 2003, 36, No. 22, 8449; doi:10.1021/ma020166f. 37. Vorobyova O., Lau W., Winnik M. A., Langmuir, 2001, 17, No. 5, 1357; doi:10.1021/la0011026. 38. Serero Y., Aznar R., Porte G., Berret J.-F., Calvet D., Collet A., Viguier M., Phys. Rev. Lett., 1998, 81, No. 25, 5584; doi:10.1103/PhysRevLett.81.5584. 39. Meng. X.-X., Russel W. B., Macromolecules, 2005, 38, No. 2, 593; doi:10.1021/ma048968t. 40. Matsen M. W., Thompson R. B., J. Chem. Phys., 1999, 111, No. 15, 7139; doi:10.1063/1.480006. 41. Bras R. E., Shull K. R., Macromolecules, 2009, 42, No. 21, 8513; doi:10.1021/ma9013687. 43601-10 http://dx.doi.org/10.1016/S0920-4105(97)00048-X http://dx.doi.org/10.1021/cr990125q http://dx.doi.org/10.1038/453171a http://dx.doi.org/10.1038/nature06669 http://dx.doi.org/10.1021/ma951174h http://dx.doi.org/10.1021/ma970354j http://dx.doi.org/10.1063/1.3400648 http://dx.doi.org/10.1021/ma9802273 http://dx.doi.org/10.1021/ma011027l http://dx.doi.org/10.1006/jcis.2002.8402 http://dx.doi.org/10.1021/ma021076d http://dx.doi.org/10.1021/ma0210776 http://dx.doi.org/10.1063/1.461276 http://dx.doi.org/10.1063/1.461969 http://dx.doi.org/10.1021/ma9508461 http://dx.doi.org/10.1103/PhysRevLett.72.2660 http://dx.doi.org/10.1103/PhysRevE.69.031803 http://dx.doi.org/10.1021/jp0359337 http://dx.doi.org/10.1021/ma061493g http://dx.doi.org/10.1021/ma070928c http://dx.doi.org/10.1021/ma800130q http://dx.doi.org/10.1103/PhysRevLett.87.188301 http://dx.doi.org/10.1063/1.2176619 http://dx.doi.org/10.1063/1.454931 http://dx.doi.org/10.1093/acprof:oso/9780198567295.001.0001 http://dx.doi.org/10.1021/ma011515t http://dx.doi.org/10.1051/jphys:01976003707-8097300 http://dx.doi.org/10.1021/ma60048a024 http://dx.doi.org/10.1021/ma60065a002 http://dx.doi.org/10.1063/1.480004 http://dx.doi.org/10.1063/1.444800 http://dx.doi.org/10.1021/ma020166f http://dx.doi.org/10.1021/la0011026 http://dx.doi.org/10.1103/PhysRevLett.81.5584 http://dx.doi.org/10.1021/ma048968t http://dx.doi.org/10.1063/1.480006 http://dx.doi.org/10.1021/ma9013687 Effects of polymer concentration and chain length on aggregation Вплив концентрацiї полiмера i довжини ланцюга на агрегацiю у фiзично асоцiйованих розчинах полiмерiв К.Ґ. Ган1, К.-Ф. Жанг1, Й.-Г. Ма1, С.-К. Жан2, Й.-Б. Ґуан2 1 Школа математики, фiзики i бiотехнологiй, Унiверситет науки i технологiй Внутрiшньої Монголiї, 014010 Баоту, Китай 2 Фiзичний факультет, Унiверситет Джiлiн, 130021 Чанґчун, Китай Використовуючи ґраткову модель самоузгодженого поля, вивчається вплив концентрацiї полiмера i довжини ланцюга на агрегацiю в асоцiативних полiмерних розчинах. У системах з рiзною довжиною ланцюга спостерiгається тiльки двi неоднорiднi морфологiї, а саме мiкрофлуктуацiйна однорiдна (MFH) i мiцелярна морфологiї. Температури, при яких вище згаданi двi морфологiї виникають впер- ше i позначаються як TMFH i Tm, вiдповiдно, є незалежними вiд концентрацiї полiмера i довжини ланцюга. Змiна логарифма критичної концентрацiї MFH зi змiною логарифма довжини ланцюга за- довiльняє спiввiдношенню лiнiйного допасування з нахилом рiвним −1. Крiм того, змiна середньої об’ємної фракцiї стiкерiв при мiцелярному корi (AVFSM) з концентрацiєю полiмера i довжиною лан- цюга сфокусована в системi при Tm. Знайдено шляхом розрахункiв, що вище згадана поведiнка AVFSM пояснюється в термiнах iнтраланцюгових та iнтерланцюгових асоцiацiй. Ключовi слова: концентрацiя, довжина ланцюга, агрегацiя, асоцiативний полiмер 43601-11 Introduction Theory Result and discussion Conclusion and summary