1. Introduction
Currently, the theory of de Gennes, Doi and Edwards is considered to be the most advanced set of ideas in the physics of entangled polymers [
1,
2]. It replaces the complex system of linear polymer chains with a mean-field approach, employing the concept of a tube in order to represent the molecular environment of an entangled chain. Due to the simplifications it offers, the theory of de Gennes, Doi and Edwards, provides a powerful tool for applying rheology as an analytical method in industry, and the majority of research within the field of rheology follows this widely-accepted theoretical picture in order to describe the dynamics of macromolecules [
3,
4,
5,
6,
7,
8].
De Gennes [
1] was the first to simplify the complex many-body problem of a large number of interacting polymer chains diffusing in a melt. His strategy was to introduce a mean-field approach by reducing the many-body problem (all chains) to a one-body problem (“test” chain) in an average effective field. The topological constraints imposed by the surrounding chains on the “test” chain are assumed to be equivalent to the tube, in which the “test” chain is confined since it cannot move in the direction perpendicular to the tube axis. At equilibrium, the tube is assumed to be a random Gaussian walk consisting of
Z steps of length
a and average length
[
2]. Nevertheless, when a rapid strain is imposed on the melt, the tube is deformed and the assumption of the random walk is not valid anymore. The following relaxation towards equilibrium occurs via a change in the configuration of the tube which confines the “test” chain. By wriggling back and forth, the chain vacates parts of the original deformed tube and the stress linked to those portions is “forgotten”. In other words, the chain is forced to relax its stress by reptation, that is curvilinear Brownian motion along the tube axis. Therefore, the unreleased stress in the system,
, is proportional to
, the fraction of the initial tube that remains occupied by the chain at time
t.
Originally, the tube model was designed to describe the motion of a long polymer chain in a fixed cross-linked rubber network [
1]. Further development of the concept of reptation by Doi and Edwards [
2] led to a full theoretical framework which is a cornerstone of the field of rheology of monodisperse linear entangled melts. To explain and fix some drawbacks in the predicted results several relaxation mechanisms were subsequently added to the reptation theory. For instance, the time it takes for the chain to diffuse out of the initial stressed tube
was predicted to be proportional to the molecular mass
M as
, while experimental data showed a
dependency [
9]. The physical origin of this difference was explained later as some subsidiary relaxation process, which is faster than reptation [
10]. This process is called primitive-path or contour-length fluctuations (CLF), according to which, due to Brownian motion, the ends of the chain can wrinkle quite some distance inside the tube, thereby vacating part of the tube end. When one end is pulled out back, it occupies a new random tube segment and part of the stress associated with the end of the original tube is lost. Thus, the contour length of the chain
L fluctuates around the value
. Doi [
11,
12] performed the first calculation of CLF and analysed their effect upon reptation. The essential point of Doi’s arguments is that CLF reduces the reptation time
for medium-sized chains (number of entanglements
). His explanation is based on the idea that if the chain ends fluctuate rapidly over a distance
, the chain needs to reptate only over a distance
. As a result, the correction due to the contour length fluctuation is of the order of
and the impact of the CLF reduces as the chain length increases. Different estimates of the CLF effect were later obtained by des Cloizeaux [
13], Milner and McLeish [
3], and Likhtman and McLeish [
4]. An alternative view on the origin of 3.4 exponent for
was also proposed by Liu et al. [
14]. Their results indicate that the motions of surrounding chains have a significant influence on the latter dependency and the contribution of CLF to the
scaling is quite weak.
Another mechanism, which needs to be taken into account, is the constraint release (CR) process, which emerges from the mobility of the surrounding chains. It is also included in the original tube theory [
15]. Chains that form the tube of the “test” chain are moving as well in their own tubes. Thus, the tube itself evolves in time by releasing the topological constraints on the “test” chain and allowing the latter one to move on a large scale perpendicular to the tube axis. Marrucci [
16] proposed that the process of constraint release speeds up the relaxation of the melt, considering that the relaxed portion of the chain behaves as solvent for the unrelaxed part of the molecule. CR is then modeled through “dynamic tube dilation” (DTD) [
5,
15,
17] as well as “CR-Rouse motion” [
15,
18]. Within this framework, the relaxation modulus is given as
, where
is the fraction of the material which is not yet relaxed, and
is the volume fraction which corresponds to the tube diameter such that
. Under the assumption that an entanglement is formed by two chains, the dynamic dilution exponent
becomes 1. On the other hand, Colby and Rubinstein [
19] suggest that
should equal
by considering that an entanglement is formed by some fixed number of binary contacts between the chains. It is important to note here that constraint release effects are also present in monodisperse samples. Indeed, as it was recently shown by Matsumiya et al. [
20], viscoelastic relaxation process in the monodisperse system is accelerated by CR, and suppression of this mechanism is clearly observed if the chains are blended into a very long linear matrix. Thus, even for monodisperse samples, a good description of CR is needed in order to reach quantitative level of predictions.
All tube-based computational models that predict the linear viscoelasticity (LVE) agree that, in addition to the high frequency Rouse modes, the relaxation of linear polymer melts consists of the three main processes described above: reptation, contour length fluctuations and constraint release. For example, these mechanisms form the basis of the quantitative theory developed by Milner and McLeish (MM theory) for monodisperse linear polymers [
3]. In this model, a linear chain is represented as a two-arm star polymer able to reptate. In this way, relaxation starts with contour length fluctuations and retraction (activated CLF) of the chain ends, and at time
is followed by reptation along the unrelaxed part, as retraction is becoming slower than reptation. The Hierarchical model of Larson [
6] incorporates the physics of the MM theory and together with some additional features forms a model which can be applied to linear polymers as well as to mixtures of branched polymer molecules. This model was modified by Park [
21,
22,
23] to include the “early-time” fluctuations and some other refinements. The BoB model, which also uses the MM theory, was developed by Das et al. [
7] and incorporates modeling of branch on branch structures.
While all models use a common theoretical approach for the LVE of linear entangled polymers, they differ in some important assumptions [
24]. Notably, the definition of CLF in models which are based on the MM theory, contradicts the original idea developed by Doi [
2], according to which only the fast fluctuations at the end of the chain shorten the contour length of the tube. Therefore, although theoretically it is based on Doi’s assumption [
12], numerically the amount of the material that is being subtracted from the reptation time seems to be overestimated. The effect of this assumption is discussed in detail in
Section 4. On the other hand, as is discussed in
Section 5, the MM theory neglects the constraint release modes arising from reptation. In the present work, we investigate the importance of this modeling choice and its impact on the predictions of LVE by comparing the results obtained with the BoB model to a large set of experimental data. Tube models also require the definition of an expression for the relaxation modulus. The BoB and the Hierarchical models use the same expression as the one used in the MM theory. In
Section 6, we discuss the accuracy of this expression, and propose an alternative derivation of this relaxation modulus.
Another LVE theory was proposed by Likhtman and McLeish [
4] (LM theory). It involves all the relevant mechanisms considered above and does not include any new physical processes. This theory solves the existing mathematical models for the relaxation processes in a consistent and more accurate way. In particular, the expression for the amount of unrelaxed material at early times is derived by solving the one-chain problem combining analytical and stochastic approaches. The resulting function slightly enforces the influence of early fluctuations, as compared to the one used in the other models, and the transition between CLF and reptation regimes always happens at time
, where
is the Rouse time. As discussed in
Section 5, this estimate is much closer to the original idea of CLF, than the estimates employed in the BoB and Hierarchical models.
As a third approach, the time-marching algorithm (TMA) [
8] for monodisperse linear polymers is also a tube-based model, but in fact, it is quite different from the other models. For example, CLF and reptation processes act simultaneously, while the MM, Hierarchical and BoB models consider reptation and CLF as mutually exclusive processes with a transition time at
. From cursory analysis, there seems to be a substantial distinction between TMA and other models. For instance, reptation is often a slower process in TMA, however arm fluctuation can speed up the terminal relaxation indirectly through
, which is the fraction of material not yet relaxed by reptation or retraction. It must be noted that, in addition to these relaxation processes, all above models consider the high frequency Rouse relaxation in a similar way [
4].
Therefore, a first aim of this paper is to discuss the main differences between TMA approach and the models based on the MM theory, as well as the influence that each assumption has on the relaxation process. In particular, we would like to understand why all these models, despite their large differences, lead to similar predictions in case of linear polymers. To this end, we collected a large set of experimental data from literature, that we use, first, to discuss the material parameters of the different polymers, then to validate the different approaches. This comparison allows us to address two important questions: (1) what is the importance of CLF in the relaxation process of a linear chain? (2) how does it affect the terminal relaxation time of the chains? Based on our conclusions, we then propose to modify the TMA model. Another purpose of this work is to draw attention to the differences observed between the theory and the numerical implementation for the BoB and the Hierarchical models.
From a more general perspective, the objective of the present work is to revisit the existing tube models and propose improvements in order to be able to extend them towards a quantitative prediction of the viscoelastic properties of polydisperse polymers. Indeed, while this work only focuses on monodisperse linear polymers, it allows us to discuss in depth three issues which are intrinsic to the models, and which can have very large consequences on the prediction of the relaxation of more complex polymers.
The first question is related to the balance between reptation and fluctuation processes, which is of primary importance if we want to be able to predict the viscoelastic behavior of both linear and star-like polymers while keeping the same parameters [
25]. And this first step is of primary importance if we want to be able to predict the relaxation behavior of binary blends composed of linear and star-like polymers. As mentioned in ref. [
26], today, most of the tube models fail at doing so, which is probably the strongest limitation for applying such models to polydisperse samples.
The second question that we would like to address here is the too fast relaxation of very short polymer chains [
27], which leads to a wrong description of the plateau modulus [
14,
28,
29]. This requires a careful description of the chain dynamics at short times, in order to prevent the possible disentanglement of the chains before the Rouse relaxation of their entanglement segments.
The third question is related to the influence of constraint release (CR) process, and how this one is taken into account in the models. As already mentioned, its effect is quite large, even in the case of monodisperse linear polymers. By combining the present results obtained for monodisperse linear polymers with the ones coming from recent studies on the viscoelastic properties of binary blends of entangled linear polymers [
30,
31,
32,
33,
34], in which mechanisms such as the Constraint Release Rouse process and the dynamic tube dilution have been studied in detail, we believe all the ingredients will be there in order to extend these studies to the linear viscoelastic properties of polydisperse samples.
This paper is organized as follows. First, in
Section 2, we present the time-marching algorithm (TMA) for the case of linear macromolecules. Second, in
Section 3, we present the experimental data together with the material parameters used within the present study. In
Section 4 and
Section 5, we compare the theoretical predictions of TMA with LVE data of polybutadiene, polyisoprene and polystyrene polymer samples and propose a novel modification to the model. Then, in
Section 6, we discuss the main differences between the TMA and the BoB model, focusing on the way CLF is treated, and on the way the relaxation modulus is derived. Additionally, we propose changes to the BoB model, so that the implementation is consistent with the theory.
Section 7 compares the results of the TMA and BoB predictions for the various sets of experimental data, and finally conclusions are presented at the end of the article.
4. Renormalization of the Coordinate System
In our previous study [
25], we propose a renormalization of the segment coordinates in order to improve LVE predictions of the TMA model for short star-shaped polymers (arm length
). Indeed, we noticed that by time
, a large portion (
) of the shorter chains has already been relaxed due to CLF. A similar conclusion was made by van Ruymbeke et al. [
35] while investigating discrepancies in tube-model predictions regarding the plateau modulus of short linear chains. This fast relaxation has also been reported in several works, such as the work by Liu et al. [
14,
28] and by Hou et al. [
29] when comparing the LM theory to molecular dynamics (MD) simulations of bead-spring polymer melts. The early relaxation, observed at time smaller than
, predicted by the tube-based models contradicts the tube picture, since, at time
, chains should only relax by Rouse motion, and are unaware of the surrounding chains. However, classical tube-based models consider that relaxation by CLF starts immediately at
. This effect is especially important for short chains, as observed by Bacova et al. [
54] in their MD simulation data analyses of the entangled linear and star polymers. Furthermore, the latter study reveals that the primitive path near the chain ends is almost fully relaxed by simple Rouse dynamics at
, since the outer segments are weakly affected by the topological constraints. Therefore, the decay of CLF at time earlier than
should be omitted, since it is already taken into account by the high-frequency Rouse relaxation. Otherwise, the chain relaxation is underestimated and the level of the plateau modulus is predicted too low. This assumption strongly affects the relaxation of short molecules, and especially star-shaped molecules since segmental fluctuation times are exponentially dependent on the square of their depth and the unrelaxed fraction [
25].
In fact, many works on poorly entangled chains have shown that the dynamics of the outer segments should be treated differently than the dynamics of the middle part of a chain due to the variant nature of entanglements. It seems that the chain tips show a Rouse-like relaxation. Here, we briefly review these works, before discussing how tube models can be corrected in order to avoid double counting of the relaxation of the outer segments (i.e., by Rouse and CLF). This issue was well described by Kavassalis and Noolandi [
55,
56,
57], who presented a packing entanglement theory, where the loose chain ends are not expected to contribute to the formation of effective entanglements. High mobility of the chain ends prevents them from forming long-living connections with other chains. According to their model, an effective entanglement is defined by the presence of a certain number of non-tail neighboring polymer segments; the concentration of the latter segments scales as
. The average tail length is
. Thus, as chains get shorter, more and more dangling ends are excluded from contributing to entanglements, and the entanglement spacing
varies with molecular weight
M. To support this assumption, Rathgeber et al. [
58] used neutron spin-echo spectroscopy (NSE) to investigate on a microscopic length scale the influence of topological constraints on the single-chain dynamics. To this end, they varied the mass fraction of long tracer chains in a short-chain matrix. As a result, the authors found that the number of end segments relative to the number of inner segments increases with decreasing tracer fraction, thus weakening the topological constraints. Consequently, the value of
and the concentration of longer chains are directly proportional, since shorter chains create fewer entanglements. Taking chain-end effects into account, Rathgeber [
58] obtained an excellent description of the NSE data.
In addition to rheology and NSE measurements, nuclear magnetic resonance (NMR) T
relaxation time studies [
59,
60] can provide information on the chain mobility. Most NMR experiments measure segmental dynamics ultimately in terms of an orientation autocorrelation function [
61]. Highly mobile polymer chain ends exhibit an isotropic motion (fast reorientation), which shows a slow transverse relaxation with a long corresponding relaxation time T
. At the same time, the entanglement network induces an anisotropy of chain motion on the central segments. Segments located near chain ends are expected to be subject to molecular motions different than those in the central block, provided that the molecular mass of the chain is above the critical value
, in order to ensure entangled dynamics. As far as NMR is concerned, Kimmich et al. [
62] proposed that the outer entanglement lengths of a reptating chain behave in a similar fashion to a Rouse chain not subjected to topological constraints. Indeed, through analysis of bi-exponential fits to experimental data, they showed that the effective molecular weight of the end performing Rouse motion is independent of the molecular weight of the chain, and approximately equal to
. This is consistent with the fact that unentangled linear chains with a molecular weight approaching
, display Rouse dynamics.
Moreover, also based on NMR, Guillermo and Addad [
63] showed that it is possible to represent a chain as two end-submolecules, independent of
M, and one inner part with clearly distinct dynamics. The idea that the end sections are undergoing fast Rouse-like dynamics, whereas the central section is undergoing slower dynamics is also confirmed by Ries et al. [
64] and by Chavez and Saalwachter [
61]. The latter research showed that during the constraint Rouse regime the chain is less restricted than predicted by the tube model. Wang et al. [
65] argue that this conclusion reflects the importance of ruling out the chain end contributions. Nevertheless, today, the main problem is that the transition between isotropic and anisotropic parts cannot be clearly defined [
60,
66]. Furthermore, it is not completely clear to which extent the constrained Rouse motion is observed by NMR as isotropic or not. However, it is clear from the above that it is important to separate contribution of end monomers to the relaxation mechanism and include it to the current TMA model, even if quantitative description of chain-end effect still needs to be refined by NMR experiments or computer simulations.
A first solution was proposed by van Ruymbeke et al. [
35], in order to slow down the CLF process: (1) adding to the chain an abstract sub-chain of length
, thus forming an “equivalent” chain, and, (2) rescaling the coordinate system. As a result, the fluctuation time of the original chain end is equal to
in the new coordinate system of the “equivalent” molecule. The obtained result is in good agreement with the experimental data for short linear chains. Nevertheless, this procedure is rather a numerical technique for being consistent with the tube theory. Instead, to attain consistency with our previous works, we propose here to account for the role of chain ends in the relaxation process through three considerations: (1) effective entanglements are only formed by the chain segments with relaxation time higher than
, in other words, chain ends are not part of the entanglement network due to their high mobility, (2) they relax via free Rouse-like motion, and therefore, (3) they are neglected in the calculations of the relaxation of the chain by fluctuations or reptation. To this end, we introduced [
25] a new coordinate system with the origin
localized at the molecular segment
, which has a relaxation time
. This correction is in qualitative agreement with the approach of van Ruymbeke et al. [
35], while it includes the specific influence of chain ends on the entanglement network as the chain length is shortened. In this way, the latter effect makes our approach consistent with the tube theory.
The above mentioned discrepancies in tube-model predictions regarding the plateau modulus of short linear chains are illustrated in
Figure 2. In this figure, the original TMA predictions (dash-dotted line) of the storage modulus are compared to experimental results for short PBD polymers. The results are obtained with material parameters from
Table 4. Whereas the agreement between theory and experiments is relatively good, the plateau modulus is poorly described for really short chains, like PBD8 and PBD21, due to the fact that the chain starts to relax their tube before
. Hence, using the new coordinate system in order to prevent this too fast relaxation seems important.
In our work on star-shaped polymers [
25], we showed that part of the tube end of length
is relaxed by time
via the free Rouse motion. This means that the relaxation of the ends is already included in high-frequency Rouse relaxation defined by Equation (6) and must be excluded from the determination of the relaxation function
in Equation (3). Here, we use the same result, which is in good agreement with the estimation provided by the theory of Kavassalis and Noolandi [
55,
56,
57]. This requires working in a new coordinate system, as already mentioned.
Indeed, once relaxation reaches
, the original tube with equilibrium length
will already lose its ends up to segment
, and the molecular weight between entanglements needs to be renormalized as
. The segmental Rouse time is also rescaled as
and the Rouse relaxation process is postponed till
, such that segments up to
are relaxed. Thus,
represents the deepest segment which relaxes by high-frequency Rouse relaxation. Segments after
are relaxing in a tube. Therefore, after a time
, only the central part of the chain with molecular weight of
is confined in the tube. This means that the relaxation function
is calculated over the interval
. Based on the above arguments, following on previous work [
25], retraction and reptation should be calculated in a new coordinate system with the origin
, localized at the tube segment
, which has a relaxation time
in the original coordinate system. However, to keep the calculations simple and due to negligible difference in predictions made using
and
, we start the new coordinate system at
as indicated in
Figure 3. The same holds for
, where
is used instead.
Following the classical tube theory, early fluctuations are not influenced by dilution. The relaxation time
is found for
, but only values from
are valid for the calculations. Thus, for the new coordinate system Equation (11) should be rescaled as
The above renormalization ensures that at time
the tube in the new coordinate system is yet to start relaxing. Additionally, it only affects the properties of the short chains (
). Results obtained with the new coordinate system are shown in
Figure 2 (dashed line). The agreement between the experimental plateau modulus of short linear chains and the theoretical one is certainly improved. It must be noted that in our initial work on linear binary blends [
34] the prefactor 0.5 in the fluctuations process of TMA could have been avoided by considering the coordinate transformation under discussion, along with considering the effects of deep retraction process.
5. Comprehension of CLF in the Tube Models
Besides the importance of avoiding CLF before
, another important question to address is the transition between fluctuations and reptation processes, and how it is taken into account in the different models. To this end, we first discuss the comparison between the experimental and theoretical PBD viscoelastic data presented in
Figure 2. For the theoretical predictions of
and
calculated by the TMA model with the coordinate renormalization as described in
Section 4, we see that the modeling results (dashed line) correctly predict the plateau modulus, but show a too slow relaxation in the terminal zone for the shorter polymers.
This deviation can be attributed to the fact that the influence of CLF process on reptation is neglected. Indeed, although the TMA model incorporates the retraction of the chain ends, this effect is not accounted for as tube length fluctuations that speed up reptation. This contradicts, the original idea behind CLF [
2,
11,
12], which states that only in the absence of fluctuations the disengagement time is estimated as the time
necessary for the chain to move by distance
L. In case of fluctuations, the chain can disengage from the tube when it moves the distance
. However, only under the condition that the chain ends are fluctuating rapidly over the distance
. Consequently, the disengagement time should be
[
2]. While, in TMA, the reptation time in Equation (14) involves the entire chain length
L, most models consider a linear chain being relaxed by reptation, when its center of mass diffuses along the tube over a shorter distance
. In the case of MM-based theories,
is designed as the whole fraction of the tube that did not survive from the full arm retraction by that time. On the other hand, the Likhtman–McLeish (LM) approach has a special scaling form for the effective tube length:
Introducing CLF into reptation allowed explaining the experimentally observed dependency
[
9]. Doi suggested [
12] that the effective length,
, decreases with the chain length as
and
X being a constant larger than
. Stephanou predicted [
67] the value of
X from the primitive path analysis to be close to
.
In
Figure 4 we plot the estimations of the effective fraction of the tube
participating in reptation versus
Z for polymers with a molecular weight up to 400 entanglements. As expected, the Likhtman–McLeish results are close to Doi’s, where the chain ends move quickly on the time scale of the Rouse time
. On the other hand, the primitive path (PP) simulations reveal a much weaker effect of CLF. Similar observations were done by Likhtman in his MD simulations of polymer melts [
68]: chains with the length of
Z = 7–15 had tube shortening by only
, instead of 32–42%, that we determined here, based on the BoB model. Conversely, according to the predictions of BoB, CLF has a much stronger impact on reptation than it was assumed initially. Prior to giving any details, we note that the BoB model correctly estimates that part of the chain which relaxes by retraction rather than by reptation (based on survival probabilities). However, subtracting this length from the length to be relaxed by reptation is, to our opinion, a too strong assumption.
At the other extreme, TMA, by excluding completely the influence of the fluctuation process on the reptation time, misses a quite important fundamental concept in its modeling. In order to include this effect, we need first to highlight an important point: only the segments relaxing fast in comparison to the reptation should contribute to the shortening of the tube during diffusion along the tube. Therefore, one should answer the question: which portion of the linear chain is fluctuating rapidly enough to be considered negligible for reptation? Similarly to the way we defined a possible solvent effect in reptation process [
8], we argue that it is estimated as the part of the tube which has relaxed due to fluctuations by time
. In the other words, a chain segment can take advantage of the widened effective tube, if its constraint release characteristic time resulting from CLF is on the same time scale as the reptation time. After calculating values of
for a wide range of
Z we follow the approach of Likhtman [
4] and conclude that the effective length of the primitive path can be well approximated by the empirical equation
In
Figure 4, the corresponding results appear as TMA2 in the new coordinate system and is very similar to the results obtained by Stephanou [
67]. Based on this result, a new reptation time in the new coordinate system can be defined as
, where
is the initial reptation time calculated in the old coordinate system according to Equation (14). In order to be consistent in a comparison of the part relaxed by CLF with the ones found with the other models, we can account for the amount of the material relaxed before
, and use the old coordinate system. In such a case, the tube is predicted to be shortened by CLF (TMA1 in
Figure 4), as
with the reptation time
. The lower limit
for CLF influence is obtained empirically and plays a role only for very short chains. The results obtained by using this formula in TMA can be seen in
Figure 2, where PBD experimental data are compared to the corresponding predictions. Now TMA properly captures the behavior of linear polymers. Note that the effect of CLF on reptation time is important only for polymers with molecular weight up to 50 entanglements. For higher molecular weights it becomes negligible.
6. Comparison between the TMA and BoB Models
In this section, we discuss important differences between the TMA and BoB models. This includes designing the expression used for describing the relaxation modulus, the influence of retraction process, as well as the way reptation and CLF are interconnected.
6.1. Relaxation Modulus
At first, we intend to show how the stress relaxation formula used in the BoB and the Hierarchical models is derived from the continuous Maxwell model. This allows us, in a second step, to show how the CR process is taken into account in this expression. The relaxation function for the viscoelastic fluid
is described by a continuous distribution of Maxwell modes, called the relaxation spectrum
.
Rheological relaxation processes tend to occur over so long periods of time that relaxation processes are usually modeled using logarithmic time. Nevertheless, in the derivations further on we use time on a linear scale. The exponential function is called the kernel. The relaxation spectrum
is independent of time and is the fundamental description of the system. It can be calculated from the Alfrey approximation [
69], where a set of exponents is replaced by a set of step functions. Simply,
at time
and
at time
. Equation (19) is then rewritten as
so that
It should be highlighted that the integral in Equation (19), with the relaxation spectrum from Equation (21), represents the stress decay through time and is always defined in the interval []. This means that a single Maxwell element is defined as an amount of the stress decrease at relaxation time . Then defines the fraction of Maxwell’s elements with relaxation times between and , namely, an amount of stress relaxed within that time interval.
Using Equation (21) with
, together with Equation (19), the BoB and the Hierarchical models, define the relaxation modulus as:
where
is the amount of unrelaxed material at time
and
is the volume fraction which corresponds to the tube diameter (“supertube” fraction). Both models start from
and at each computational step, they increase the time by
, up to
. Here,
is the reptation time, and
is the total number of steps. This means that although the simulation time is longer than
, the calculation of
and
stops after
steps. Thus, this model step is based on the assumption that reptation happens instantaneously at
and
becomes zero at the last computational step. Later, we show that this approximation leads to an overestimation of the reptation time.
Since the modeling of BoB is applied in the discrete time domain, all variables are calculated at each time step. The value of the integral in Equation (22) is approximated numerically by computing the product of the value of the integrand at the end-point of the interval by the interval length. Given the values at each time step we can approximate the integrand using a linear spline and
[
7]. Applying the numerical approximation of the derivative product rule, the integral over the interval
becomes
where
Thus, the numerical approximation of Equation (22) for the BoB model becomes
In order to separate contributions of retraction and reptation in the relaxation modulus, we rewrite Equation (22) for the time interval
such that only the reptation integrand from
till
is numerically approximated
The reptation integrand has been simplified due to the above mentioned assumption, which leads to for .
The resulting Equation (26) is nothing more than the expression for the relaxation modulus in the Milner–McLeish model. The first term is a portion of the stress relaxing by arm retraction and has a form identical to the one used in other models, including the TMA. However, the second term which is responsible for the reptation differs from the one in TMA. In the Milner–McLeish theory, this has the following form
where
is the fraction of the stress which remains after the relaxation by retraction and is identical to
in Equation (26). The sum over weighted relaxation times in Equation (27) is heavily dominated by the longest relaxation time
, and thus, can be simplified. Looking at Equation (27), we see that it excludes the contribution of constraint release to reptation, as it is explicitly stated in Milner–McLeish paper [
3]. That is why the exponent has the form
, instead of
, when constraint release is considered. While, the BoB and the Hierarchical models assume constraint release always to be active, as suggested in Equation (22), they remove its effect in the derivation of Equation (26) and use
for the reptation part .
In fact, the idea of reptation being an instantaneous process, i.e., for , appears to be a strong approximation. Indeed, at time , the stress of the segments relaxing by reptation only drops to of its initial value, meaning that should not be 0. Since no CCR regime is considered for monodisperse chains, also means that . Thus, Equation (19) with the relaxation spectrum from Equation (21) can not be used, as the maximum relaxation time of the chosen Maxwell elements should not be limited by and should be defined over the time interval []. If instead of the integral form for the relaxation modulus, one just uses the equivalent form , then the idea of , simply leads to , at . The latter conclusion is incorrect, as stress relaxes exponentially over a long period of time.
In order to resolve this issue, the approximation
can be used in a slightly different context. With the introduction of constraint release to the tube model, it is difficult to determine Maxwell elements for Equation (20), such that the maximum relaxation time is
. Instead, we define
being the amount of material with relaxation time larger than
t. Then,
and
. The resulting expression for the relaxation modulus becomes
where the second term represents relaxation via reptation. Recall that
, if CCR is inactive and consequently,
This form of relaxation modulus includes the double reptation expression, where the unrelaxed portion of the entanglements contributes to the modulus in proportion to the square of the remaining fraction of entanglements.
Keeping Equation (29) aside, we can now show that similar expression can be drawn from Equations (19) and (21) assuming that a part of the tube
is relaxed mostly by retraction up to time
and the remaining part
by reptation with the survival function given by
. Knowing that
and neglecting CCR (
), Equation (22) becomes
Note that
is constant for
and
in that interval. Thus, although the second integral is responsible for relaxation via reptation for
, integration can start from 0 rather than from
. The second term of Equation (30), can then be simplified to
, if we replace the exponential with a step function. Similarly to the
derived above, this expression reproduces the concept of double reptation. In reality, reptation should start earlier than
, when
. Empirical estimations show that
and the segment with this relaxation time is close to
, as calculated by the BoB model.
Based on this analysis, our present concern is that the BoB model fails to properly capture reptation behaviour. Due to the strong assumption, i.e., , the term in Equation (26) misses a factor of 2 in the exponent compared to Equation (29) (simplified version), and consequently, the relaxation time in the BoB and the Hierarchical models is smaller than the original by approximately a factor of 2.
Interestingly, the too slow reptation term seems to be well compensated by the too strong effect of CLF (see
Figure 4). There seems however to be no justification regarding the hypothesis that deep fluctuations of chain ends speed up the reptation process.
6.2. Arm Retraction
If we deal with the linear chains as two-arm stars, then the contour-length fluctuations of the ends are the same as the retracting motion of a star arm. The retraction dynamics of the TMA and the BoB models is already compared in our previous work [
25], where it was shown that both models successfully predict the linear viscoelastic data of star-like molecules. However, the difference in the definition of the prefactor in the first-passage time of the activated fluctuation
,
together with some small algorithmic discrepancies lead to slightly different material parameters for the two models (
). Here, the best fit to the experimental data was obtained with the same entanglement spacing parameter
for the TMA and BoB predictions (see
Table 4). It means that the deep fluctuations of linear chains in the TMA model are slightly faster than in the BoB model. Due to that fact, we expect to see a slower relaxation of shorter linear polymers (
).
6.3. Cumulative Effect of Reptation and CLF
All models, except TMA, assume that at some time, retraction of the linear polymer will be interrupted, and stress relaxation is completed by reptation using the theoretical framework developed by Doi and Edwards [
2]. Similarly, the Hierarchical and the BoB models consider retraction and reptation as mutually exclusive processes with a transition around reptation time
. On the other hand, the TMA model assumes both processes to take place simultaneously and only the long chains show this obvious crossover behaviour. At first sight, it seems that these choices are the cause of the alternative definitions of reptation time in each model. However, we showed in the previous section that the reason is different. In the present work, TMA shortens the contour length by excluding the part of the chain that has relaxed by
, whereas other models exclude much more—everything that has relaxed by
(see
Figure 4).
Naturally, it is interesting to check the viscoelastic behavior when the long chains (
) have this implicit transition between reptation and retraction in the TMA model. At each time step we calculate the probabilities
and
that the segment
x survives from relaxation by fluctuation and by reptation. Furthermore, at any time
t we can determine the so called relaxation front, a segment
, which separates relaxed and unrelaxed segments. This segment has a relaxation time
and the fraction of material not yet relaxed by reptation or retraction can be found as
. Once the reptation time
is reached, the relaxation front is at segment
. In
Figure 5 and
Figure 6, we plot the relaxation probabilities
and
for segments
at the time
they relaxed, and
and
for segments
at the reptation time
estimated by Equation (14). Thus, the observed decrease of
with the depth of the segment,
x, is due to the fact that for
, the time at which this probability is determined increases with
x. By comparing
to
one can notice that both processes cross over around segment
. This indicates that it is more probable that segments
are relaxed by fluctuations rather than by reptation, and segments
are relaxed primarily via reptation, as
. This estimation of transition between relaxation process is identical to the one from the BoB model (see
Figure 4). Fluctuations stop at segment
with relaxation time
and reptation takes over for segments
.
If we analyse the short chains (see
Figure 7), there is no clear separation between fluctuations and reptation because the probability of being relaxed by retraction
for all segments. However, for
the chance to be relaxed by reptation is very low, same as for the long chains. The rest of the segments are relaxed via reptation and retraction, as
.
In fact, it is easy to show the numerical equivalence of the TMA and the BoB models, especially for long chains (). To illustrate this conclusion, we add several features from the BoB model to the TMA model and compare the results. First, we modify the TMA algorithm (TMA) in order to allow reptation to start only after the arm fluctuation is finished. This transition from fluctuation process to reptation takes place at segment . In other words, the joint survival probability becomes , if and , if , where . It must be noted that in the TMA is calculated at each time step, so we need to determine a priori. The TMA model which accounts for this first modification is called TMA. In addition, we examine how the choice of the relaxation modulus equation used in the BoB model influences the relaxation process. Therefore, we introduce one extra modification in our model, called TMA, where the reptation time and of TMA are now determined according to BoB’s formulations, in addition to preventing fluctuations and reptation process to take place simultaneously. This version of the TMA model is almost identical to the BoB model, except for the transition between the early and the deep fluctuations and calculation of the relaxation time for the latter one, as was explained before.
Figure 8 shows a comparison between TMA (solid line), TMA
(dotted line) and TMA
(dashed line) predictions for Auhl’s PI samples from
Table 2. Both TMA and TMA
models predict nearly the same results for samples with high number of entanglements. The explanation goes as follows: by using the joint survival probability
, TMA introduces a smooth transition between retraction and reptation processes. The domination of the arm retraction (
for all unrelaxed segments
x) is followed by a period when reptation starts to take over; namely
tends to
for segments of the relaxation front. Finally, a chain will complete its relaxation by reptation, rather than by retraction, since
for unrelaxed segments. In general, this rule holds for linear chains with a molecular weight of 50 entanglements or higher.
A different behavior is observed for shorter chains, where TMA shows a slower decay of and in the terminal zone. Although TMA and TMA have the same reptation time, there is no actual transition between processes in the TMA model, as in case of the long chains. Relaxation of segments after is also governed by those two processes, rather than by pure reptation, because probabilities of being relaxed by reptation or retraction are almost equivalent for those segments. This explains the difference in the obtained results of the two models, as in TMA relaxation ends with pure reptation.
Interestingly, comparing these data with the TMA model allows us to illustrate how the two main approximations discussed above, i.e., a reptation process two times slower and an overestimated influence of CLF on reptation, compensate each other in the BoB model. Indeed, BoB can be viewed as the TMA model with the reptation time . Nevertheless, samples having a high molecular weight show slightly slower terminal relaxation with BoB. Indeed, when for short chains, TMA still has retraction, which implicitly reduces reptation time to the level of . For very long chains the reptation relaxation time ratio is slightly less than 2 and leads to longer terminal relaxation in the BoB model.
In the BoB model, the missing constraint release in reptation is compensated by choosing a smaller value for by excluding the part of the chain relaxed by deep fluctuations, such that a relationship between TMA and BoB reptation times of is found. Here, we argue against this hypothesis. Indeed, during retraction, one end of the chain will be first dragged inside the tube over a maximum distance by forming some big loops between entanglements. However, at the same time, the other end will be pulled out of the tube by reptation over the same distance and will first stretch those loops rather than translating them over the distance , so that the fluctuating tip of the chain will stay fixed. Therefore, on average, the equilibrium length of the chain stays constant. Only the fast fluctuations (compared to the reptation time) should speed up reptation. Furthermore, deep retraction cannot be considered as a fast process compared to reptation.
7. Comparison of Theoretical and Experimental LVE Data
First, we summarize in
Table 5 the important features of the TMA, the modified TMA and the Bob which were discussed in this paper.
The modified TMA and the BoB models with the input parameters presented in
Table 4 are used to compute
. The comparison between experimental and theoretical data for PI samples is shown in
Figure 9 and
Figure 10.
The overall agreement between the Auhl’s PI data and the more advanced version of the TMA model is remarkable. The BoB model exhibits a deviation from the experimental plateau modulus for the short chains. This is explained by accounting twice the effect of chain ends relaxation as discussed in
Section 4. Also, the predictions of BoB for short chains (
) are relaxing a bit slower than the ones from TMA. We identified the nature of this discrepancy previously: the largest part of the chain is relaxed by retraction, and BoB would need larger
to coincide with TMA predictions [
25]. For middle range molecular weight samples, the TMA and the BoB models are found to yield similar results when using the same set of input parameters.
All the PBD data are predicted with the same values for the input parameters
and
. Recall that Likhtman and McLeish [
4] had to set the corresponding values differently for Baumgaertel’s [
37] and Lee’s [
38] data. The results in
Figure 11,
Figure 12,
Figure 13,
Figure 14 and
Figure 15 show the same quality of predictions as for PI samples. All deviations between the experimental data and the BoB model are already explained above.
Finally, in
Figure 16 and
Figure 17 we compare our predictions for linear viscoelasticity with the experimental data of linear PS samples coming from several research groups.
Again, the predictions made by the TMA model are in a good agreement with all data, while the BoB model overpredicts the relaxation of really short chains together with a clear deviation from the plateau modulus. The data from Graessley [
47] and Schausberger [
46] is also considered in the work of Likhtman and McLeish [
4]. To fit the data accurately, they had to adjust the values of
,
and
, separately for each set. On the contrary, we use the same value of
and
for the calculations employing these samples.
There are several comments regarding the comparison of the data. First, the shape of the relaxation spectra
and
predicted by TMA is in good agreement with the experimental data. Second, the agreement between TMA and BoB is quite reasonable, but with some deviations. The most pronounced disagreement is observed for the short-chain samples, where the BoB model overestimates the terminal relaxation time, and predictions are slightly shifted downwards. Presumably, this effect can be attributed to the difference in the coordinate systems used in both models, since TMA
shows a better fit. As discussed in
Section 4, the value of
Z for short chains is somewhat lower due to the renormalization. Note that the BoB model cannot predict linear and star polymers with the same
value (
). Thus, BoB predictions could lead to a better agreement for these chains, if the values of
for the star-shaped polymers were used. This fix is consistent with the fact that the relaxation of the short linear polymers is mostly governed by fluctuations. Third, the distinction between BoB and TMA is expected to be more pronounced for highly entangled samples (
). Nonetheless, the lack of experimental data prevent this comparison.
8. Conclusions
Although most of tube-based models produce similar quantitative results for the monodisperse linear viscoelastic behavior of linear macromolecules, there are some crucial differences between them; especially between TMA and the others. In this paper, we compared the theoretical frameworks for several well-known models: Milner–McLeish, Likhtman–McLeish, the Hierarchical model, the BoB model and the TMA model. Our comparison revealed that the most pronounced distinction between these models is the way the dynamic modulus is determined, and as a consequence, the relaxation time . Additionally, the description of the transition between relaxation processes also plays an important role in the predictions, especially for shorter polymers.
Our analyses indicate that the BoB and the Hierarchical models use the same expression for the relaxation modulus as MM theory, i.e., excluding the effect of constraint release modes arising from reptation of other chains. This assumption is compensated in the MM theory, the BoB and the Hierarchical models by including the possible deep retraction of the chain extremities into CLF calculations and by removing this CLF from the reptation time. In this work, we showed that only fast fluctuations of chain ends should speed up reptation, as initially proposed in the tube theory.
On the other hand, we discussed the consequences of an incorrect approximation applied in the derivation of the relaxation modulus and proposed several alternative ways to obtain the correct expression for calculation. We also showed that, the crude assumption made in the BoB model does not dramatically affect the quality of the results for linear chains. However, one can wonder what would be their impact in the case of binary blends or for very polydisperse samples. Furthermore, modeling of branched polymers might also be affected.
Based on the comparison to experimental data, we added several modifications to the time-marching algorithm (TMA) for the case of monodisperse linear molecules. Thanks to renormalization of the coordinate system, we avoid considering the segments which relax before in further calculations. This procedure helped us to improve the performance of the TMA model and to eliminate the discrepancy of tube model predictions for describing the plateau modulus for short chains. Furthermore, we included the influence of the fluctuation process on reptation time in the modeling. This quite important fundamental concept was missing in the TMA model. The amount of material excluded from the reptation process is consistent with the original CLF concept and agrees with values extracted from MD simulations.
We validated the modified TMA model on a wide range of nearly monodisperse polyisoprene, polystyrene and polybutadiene linear melt samples. All predictions are in good agreement with the experimental data. Notably, same entanglement molecular weight
could be used for each chemical group, contrarily to the work of Likhtman and McLeish [
4], who used different values of
within a chemical group. Furthermore, we were able to use the same parameters as the ones used to describe the viscoelastic properties of start polymers. Therefore, this gives us a solid base to extend tube models to polydisperse blends of linear and star polymers. Validation of the BoB model was performed on the same set of experimental data. The numerical predictions of the BoB model also fit the data quite well, apart from the deviations in the case of short chains. The latter observation can be attributed to the fact that in the BoB model, stress relaxation by contour length fluctuations is calculated starting at time
.