3.1.1. Analytical Approaches
The first investigation of the intramolecular dynamics of dendrimers was carried out by La Ferla [
19] within the framework of the Rouse-Zimm approach. The modelled dendrimers had a ternary core (
f = 3) and binary dendra (
m = 2), as shown in
Figure 1, resulting in a variable number of bonds
p between consecutive branch points, while the generations
g were numbered from 0 onward. The intramolecular dynamics were described on the basis of the Langevin equation. As for the chosen conformational model, which determines the intramolecular elastic forces acting on each bead, La Ferla adopted a freely rotating model, whereby the connected springs associated with each chain segment and which have one bead in common exhibit fixed bond angles that may be different depending on whether the bead is the central core, another branching point, or is part of a linear portion; the simplest freely jointed model was recovered by setting the average cosines of all the bond angles to 0. Only topologically short-ranged interactions were included, while long-range excluded volume interactions were neglected. The scalar products between the bond vectors determine the intramolecular elastic forces through an incidence (or connectivity) matrix, depending only on the molecular topology. From a conformational viewpoint, La Ferla was able to obtain analytical results for the radius of gyration
in the freely jointed dendrimer that were equivalent to results obtainable from a random walk model, as a function of
f,
m,
g and the number of bonds
p between adjacent branch points [
19]. Here, we are only interested in the dynamical properties, starting from the spectrum of relaxation times and the intrinsic viscosity [η]. The relaxation spectrum shows a large degeneracy in the relaxation times τ
p, which become a universal function at high
g if they are plotted as a function of
p/
N (see
Section 2.1). Considering dendrimers having
p = 1, La Ferla was also able to determine the degeneracy of the collective modes as a function of
f,
m,
g for dendrimers, that is, based only on the dendrimer topology. The intrinsic viscosity [η] can thus be obtained through Equation (1), as it is mainly controlled by the longest relaxation times of the collective modes. Unfortunately, La Ferla did not report the
g-dependence of [η] for a comparison with experimental data, but since [η] is proportional to
/
M based on the known coefficients, where
Rη is the viscosimetric radius and
M is the molar mass, he showed that the calculated
Rη did not follow a power-law dependence on
M, unlike what is found for linear chains and star polymers. Moreover, the
Rη values, calculated for dendrimers with
f = 3 and
m = 2, exceeded the experimental values, with higher values corresponding to higher generation
g. The reason for this discrepancy was tentatively attributed to a small incomplete branching at high
g of the experimental samples, which was also implied by the small but non-negligible polydispersity, which would affect the viscosimetric radius much more that the radius of gyration.
The spectrum of relaxation times and the intrinsic viscosity were later investigated while also considering the effect of good solvent expansion [
20] in view of the large covolume effects that are present in the sterically crowded dendrimers. Intramolecular expansion was calculated through self-consistent minimization of the intramolecular free energy, accounting for the configurational entropy loss experienced by the swollen molecule and for the covolume two-body interactions. In turn, the latter interactions are calculated on the basis of the pairwise contact probability among the beads, assuming a Gaussian distribution of the (perturbed) inter-bead distances. Both these free-energy contributions can be written in terms of the scalar products among all the bond vectors connecting the beads for the topology of interest [
20,
21,
22]. As a result, the equilibrium state is determined, which corresponds to the optimal compromise between the repulsive covolume interactions, which tend to non-affinely swell the molecule, and the elastic penalty opposing it as a result of the entropic configurational loss. This approach led to an asymptotically finite expansion dictated by the finite, and quite small, span of each dendron. Moreover, the optimized scalar products among the bond vectors also yield the elastic forces acting on the beads and the inter-bead distances required to account for the hydrodynamic interaction. However, while the agreement with the experimental values for the PPI dendrimers were quite satisfactory for
, the viscosimetric radius obtained from the calculated [η] was found to exceed the experimental data at high values of
g, which is in keeping with the La Ferla’s results. Moreover, the intrinsic viscosity [η] did not show any maximum if plotted as a function of
g, and only the phantom molecule (which corresponds to a random walk) hinted at the possible presence of such a maximum, and only for unrealistically high
g values. In the same paper [
20], it was also shown, however, that slightly attractive pairwise interactions among the beads, resisted only by configurational entropy, could indeed lead to a maximum in [η]. Interestingly, such weak attractive interactions, corresponding to a slightly negative binary cluster that is integral among the beads [
13], are actually required in order to achieve an unperturbed Θ state, whereby the temperature produces a vanishing second virial coefficient between the molecules that is measured, for instance, in osmometry and in light scattering experiments [
23]. Because of the molecular topology, in branched systems, the interactions between two molecules also entail repulsive three-body interactions among the beads; in fact, two beads of different molecules cannot freely approach one another because of the covolume of a topologically neighboring third bead. Such intermolecular repulsion must be compensated by a slightly negative binary cluster that is integral to achieving the Θ state [
23]. Therefore, it can be concluded that the maximum [η] obtained in Reference [
20] corresponds to what was calculated for the Θ state. One general issue related to this point is that, in overcrowded systems such as dendrimers at high values of
g, the excluded volume effects are effectively screened out by the locally dense environment, similar to what is found in polymer melts. As for the spectrum of relaxation times, this has been analyzed in more detail in a subsequent paper, discussed below in connection with the calculations of other dynamical observables [
24], where the degeneracy of the relaxation times, related to the topological symmetry of the molecules alone was fully determined as a function of
g,
f,
m and
p.
The analysis of the normal modes of motion of dendrimers and the multiplicity of relaxation times related to molecular symmetry was investigated independently by Cai and Chen [
25] with respect to the Rouse free-draining limit, that is, neglecting hydrodynamic interactions. Moreover, the intramolecular excluded-volume interactions were ignored, so that the molecule could be conformationally described as an appropriate random walk. In this way, many quantities of interest can be analytically calculated, such as
, for instance [
19]. More interestingly, in this way, the dynamical problem can be solved quasi-analytically, and the normal modes of motion can be obtained explicitly together with the relaxation times (or at least the collective ones) for the low-generation dendrimers, and then, by extension, for higher
g values. The degeneracy of the relaxation times was largely determined in this way (see also Reference [
24] for a more systematic list), while the intrinsic viscosity could also be analytically determined as a function of the generation,
g, in consideration of dendrimers with a tri-functional core (
f = 3) and binary dendra (
m = 2); in keeping with the other investigations mentioned above, [η] was found to increase monotonically with
g, with no indication of any maximum [
25]. In a subsequent paper [
26], Cai and Chen investigated the intrinsic viscosity [η] of dendrimers in more detail, adopting the variational approach of Fixman [
27] to tackle this problem, the exact solution of which is a formidable task, usually requiring some approximation such as the pre-averaged approximation. Since the effect of this approximation on the calculation of [η] is unknown, but can be large in dense systems, Fixman proposed two different approximations that could bracket the “true” values, rather than the pre-averaged one, which provides an overestimation of [η]. Cai and Chen adopted the Rotne-Prager tensor to model hydrodynamic interaction, since it is more accurate that the Oseen tensor, but this introduced a further parameter—hydrodynamic bead diameter. If this value is equal to the excluded-volume bead diameter used to account for the bead covolume, no maximum can be obtained for [η] for either the upper or the lower limit. If, however, the hydrodynamic bead diameter is close to the segment length, while the covolume bead diameter is very small (thus approaching the Θ state), a shallow maximum can indeed be obtained. It can be noted that this procedure is quite ad hoc, and no clear molecular basis for this difference can be given on physical grounds. This problem was investigated again by the same group with Monte Carlo simulations [
28], as described in the next section.
The theoretical study of dendrimer dynamics was subsequently carried out by Biswas’ group employing the standard Rouse-Zimm approach with hydrodynamic interaction and using the pre-averaged Oseen tensor [
29]. Biswas et al. adopted a semiflexible conformational model by imposing restrictions on the direction and orientations of consecutive bond vectors, somehow accounting locally for the excluded-volume effects. The conformational model yields equilibrium quantities such as
, for instance, and determines the elastic force terms in the Langevin equation. The restrictions on the direction and orientation of the bond vectors are applied by fixing the bond angle values formed by the vectors stemming from a branch point in polar coordinates using appropriate spherical harmonics, thus describing the local correlation between consecutive bonds of different generations (since
p = 1, see
Figure 1) and allowing for more compressed or more expanded conformations. The elastic forces are then evaluated based on the average scalar products between the bond vectors, so that finally the spectrum of the relaxation times and of the normal modes of motion can be obtained. While the multiplicity of the relaxation times is not affected, the first result is that the rigidity constraints of the semiflexible dendrimers produces longer relaxation times for the local modes, but does not greatly affect the relaxation times of the collective modes. Interestingly, the resulting intrinsic viscosity shows a well-defined maximum when plotted as a function of
g, this maximum being present at
g = 7 quite independently of stiffness (generations are numbered from
g = 1 onward), and being only slightly shifted to
g = 6–7 with more compressed conformations. It should be pointed out that such
g values (calculations were carried out up to
g = 10) are quite high, producing a very high intramolecular density that is only achieved for PAMAM dendrimers that have a relatively long and flexible spacer between adjacent branch points and an EDA central unit with two functional groups, such that a value of
g = 10 was realized. It is noteworthy that a pronounced maximum was also obtained in the same paper [
29] for the model adopted by La Ferla [
19], carrying out calculation for up to
g = 10 (La Ferla considered the
g = 6 case, at most); furthermore, in this case, a maximum of [η] was obtained for
g = 7. It may therefore be concluded that in other cases, too (for instance in References [
20,
25,
28]), such a maximum could indeed be found at
g values higher than those considered. Therefore, the existence of a maximum of [η] as a function of
g could be a general feature of current approaches that are close to the ideal Θ state, even though the quantitative agreement would often be quite poor, at least at this location.
Biswas’ group also investigated the dynamical properties of randomly hyperbranched polymers along the same lines [
30,
31]. The hyperbranched polymers were closely similar to dendrimers, but were built using a growth algorithm starting from a three-functional core (
f = 3) and binary branches (
m = 2) after each bond. However, at each generation, one end unit was randomly selected as a dead end, wherefrom no further branching (or linear bonding) took place [
30], while a limited flexibility was introduced in the same way as previously described in References [
29,
32]. In the second paper on hyperbranched polymers, the elastic forces were accounted for through the elastic springs associated with each bond and through a local excluded volume term between beads belonging to the same or to the adjacent shell only with two different parameters [
31]. It may be pointed out again that this description would not account for the long-range excluded volume effects of linear polymers, since it would simply renormalize the segment length of an equivalent random walk chain. However, it could be used for dendrimers, assuming that the covolume effects are effectively screened out in dense systems. The dynamical equations, which account for hydrodynamic interactions through the pre-averaged Oseen tensor, yield the spectrum of relaxation times and the normal modes of motion for different values of the excluded volume parameters. The intrinsic viscosity [η] is then obtained as a function of
g; a clear maximum is seen for
g = 5 in the case of regular dendrimers with a significant excluded volume interaction, while in the analogous randomly hyperbranched polymer the maximum takes place for
g = 6 and is less pronounced with a weaker decrease of [η] at higher generations. Moreover, in the hyperbranched polymers, the magnitude of [η] increases with decreasing strength of the excluded volume interaction, ultimately approaching the values of star polymers with the same number of beads, with a monotonous increase of [η] at an increasing molar mass.
Dendrimer dynamics were also discussed in a further paper by Biswas’ group [
33], again using a pre-averaged hydrodynamic interaction and accounting for the excluded volume between the nearest non-bonded beads only. The strength of this interaction was characterized by a parameter derived from Flory’s mean-field approximation through minimization of the molecular free energy accounting for the configurational entropy and an excluded volume term estimated in the mean-field approximation through the volume pervaded by the polymer (a sphere having a radius equal to
). Because of that, the excluded-volume parameter depends on the number of beads, on the number of the nearest non-bonded interactions and on the sum of the distances among the bead pairs, hence on the dendrimer generation. An alternative approach adopted in the same paper consists of adopting the geometrical procedure proposed by the same group [
31] in terms of two parameters for beads belonging to the same or to the adjacent generation only (see before). With either method, in order to account for the excluded-volume parameter, a maximum of [η] is obtained for relatively small
g values. In fact, using the geometrical procedure with a small excluded volume parameter, the maximum can be seen at
g = 6, with a small shift at higher
g for a weaker excluded value strength; when using the parameter obtained with Flory’s mean-field approximation, it takes place at
g = 4 (the generations are numbered from
g = 1 onward). By tuning the excluded-volume parameters of the geometrical approach, Biswas et al. [
33] were also able to satisfactorily reproduce a large number of experimental and simulation results on the intrinsic viscosity of dendrimers.
Figure 2 summarizes the experimental results mentioned in the introduction and the fitting results of Biswas et al. [
33].
As a final theoretical approach, we should also mention an entirely different method for calculating the intrinsic viscosity of dendrimers based on a sort of mesoscopic molecular description, rather than a fully microscopic one as used before. The theory was first proposed in terms of a simple two-zone model of dendrimers [
34]. In this model, the dendrimer is described as being formed of a dense core impermeable to the flow field, where the solvent molecules are effectively trapped (thus adopting the Einstein result for the viscosity of a suspension of spheres [
13]), and a thin outer region where the solvent drains freely (i.e., in the Rouse limit). Assuming spherical symmetry and a Gaussian radial density profile, the dependence of [η] on
g could be calculated in quantitative agreement with the experimental results for PBzE dendrimers, in particular with respect to the presence of the maximum and its position. The theory and the general mesoscopic molecular picture of polymers was subsequently developed in greater detail [
35], introducing a drag function in order to describe the local volume fraction of the solvent flowing with the polymer, and a drainage function related to the solvent that moves through the polymer. Both functions are calculated from a spherical symmetric density profile, which is now determined on the basis of appropriate Monte Carlo simulations for the topologies of interest using the Lennard-Jones potential for the non-bonded interactions. The results of these calculations favorably compare semi-quantitatively with experimental data from polymers of different topologies (linear, ring, star polymers) and for the dendrimers (PAMAM, PBzE, PPI) with only a few adjustable parameters. We note that the connection of this model with a more microscopic one is quite unclear, and would require a more detailed investigation, as also indicated at the end of this review.
3.1.2. Simulation Methods
The first simulation study investigating the intrinsic viscosity of dendrimers and of hyperbranched polymers was probably carried out in 1998 using Metropolis Monte Carlo simulations [
36]. The polymers were built starting from a trifunctional B
3 core by sequentially adding AB
2 monomers assuming definite reaction probabilities at the ends of each unit. By tuning these a priori probabilities, both dendrimers and hyperbranched polymers with some dead ends could be obtained. The chosen monomers corresponded to a branched alkane, described on the basis of the Rotational Isomeric Scheme (RIS) [
37] for the rotations around single bonds, but they effectively adopted a phantom-chain model (bond crossing was not forbidden). The excluded-volume interactions were neglected, as well as the hydrodynamic interactions. The intrinsic viscosity was then calculated through the Flory-Fox equation [η] ~
/
N, with
being the radius of gyration obtained by the simulations and
N being the total number of monomers. Quite surprisingly, in view of the above-mentioned approximations and the inadequacy of using the Flory-Fox equation noted in other studies [
34,
38], this approach yielded a maximum of [η] for
g = 3 for regular dendrimers (
g was numbered from 0 onward), in a range consistent with the experimental data. A similar maximum was also predicted for hyperbranched molecules, with the maximum progressively shifting to a somewhat larger molecular weight corresponding to g ≈ 4–5 of the respective dendrimers with a decreasing branching pattern, and with larger [η] values than the corresponding dendrimers [
36].
The calculation of the spectrum of relaxation times and of intrinsic viscosity was later investigated in more detail by Cai and Chen [
28], with an analysis of the approximations involved in the analytical studies. The dendrimer conformation was modelled by standard Monte Carlo simulations in continuous space for a freely jointed model with
g = 3 and
m = 2 (see
Figure 1) and with a single bond between adjacent branch points in order to calculate averages of interest, including those yielding the elastic potential and the elements of the hydrodynamic tensor, both in the simple Oseen form and in the more sophisticated Rotne-Prager form. The dynamical equations were then solved numerically, producing the spectrum of relaxation times and the normal modes of motion. The relaxation times showed the expected degeneracies independently of the strength of the excluded volume interaction. As for the intrinsic viscosity, the pre-averaged Oseen tensor led to values that somewhat exceeded both the upper and the lower bounds calculated using Fixman’s variational method with the Rotne-Prager tensor [
27]. In any case, no maximum of [η] was obtained as a function of
g using the same bead diameter for the hydrodynamic interaction and for the excluded volume [
28].
A more complete simulation study of dendrimer dynamics was later performed by Mansfield in a seminal paper [
39], adopting lattice Monte Carlo simulations for the dendrimers. The issue of the simulation of the transport properties, including the intrinsic viscosity in particular, was dealt with using an analogy with electrostatic and random-walk statistics, so that the Navier-Stokes equation was transformed into the Laplace equation, while the intrinsic viscosity was proportional to the trace of the polarizability tensor, even though the proportionality constant depends on the shape of the body. It is not clear, however, how much this analogy depends on the step length of the random walk, and on the assumption of a continuous solvent medium, as adopted in the Rouse-Zimm approach. The simulations were carried out for dendrimers on a diamond lattice, assuming that seven bonds were present between adjacent branch points under excluded volume conditions, and showed a well-defined maximum of the intrinsic viscosity for
g = 6 [
39]. Unfortunately, the applied methodology violated the detailed balance [
40], potentially casting some doubt on the calculated properties. Further simulations by Mansfield and Jeong [
41] took this issue into account and overcame the original problem by using a different criterion for accepting the trial moves. The result was that the former Mansfield’s results were significantly affected by the detailed-balance violation for phantom-chain dendrimers, but were qualitatively valid provided the exclude volume interactions were accounted for, so that in this case a sharp maximum of [η] was calculated for
g = 6–7. No direct comparison with the experimental data was attempted, however, as the main emphasis was on the conformational equilibrium properties.
Another simulation study of the dynamical properties of dendrimers was carried out by Lyulin et al. using Brownian Dynamics simulations of dendrimers under a shear flow [
42]. Again adopting a dendrimer model with a trifunctional core (
f = 3) and binary dendra (
m = 2) comprising either one or two bonds between adjacent branch points (
p = 1 or 2), the equations of motion accounted for the hydrodynamic interaction with the Rotne-Prager tensor. The interaction potential among the beads was modelled on the basis of Lennard-Jones potential, with parameters that reproduced the unperturbed Θ state in linear chains, while a rigid constraint was applied to the bond lengths, and the solvent velocity in one direction was given in terms of an applied shear rate
. The intrinsic viscosity [η], calculated from its definition as the ratio between the appropriate non-diagonal component of the stress tensor and the applied shear rate in the limit
→0, displayed a shallow maximum for
g = 4–5. At a significantly higher shear rate, a shear thinning behavior was obtained, such that [η] was found to decrease with an increase of
according to the power-law
−1/3, with a concomitant significant increase in the molecular size in the shear direction, as expected. Hyperbranched polymers were also investigated using the same Brownian Dynamics simulations [
43], with the polymers being built by using the same algorithm as in Reference [
36]. A maximum of the intrinsic viscosity was again observed for dendrimers and for hyperbranched polymers with a large branching degree, while the shear thinning behavior was again detected at a large shear rate, although the simulations did not allow to detect a clear power-law dependence of [η] from
.
The intrinsic viscosity of dendrimers has also been modelled using Molecular Dynamics methods [
44], adopting a coarse-grained model of dendrimers in explicit solvent, modelled as Lennard-Jones particles. Dendrimers were described up to
g = 7 with a trifunctional core (
f = 3), and binary dendra (
m = 2), comprising a single harmonic spring between consecutive branch points (
p = 1), while the non-bonded beads interacted through Lennard-Jones potential. The shear viscosity was calculated with the Green-Kubo formula as the integral over time of the stress autocorrelation function, using the off-diagonal terms of the stress tensor, which in turn were obtained from the momenta of the particles, and on the interparticle forces and distances. The intrinsic viscosity thus calculated, which is in principle quite computationally demanding, showed a clear maximum at
g = 5, and in general compared very well with the Brownian Dynamics simulations [
42].
A different approach was later subsequently by Freire et al. [
45], who used coarse-grained Monte Carlo simulations in continuous space for dendrimers with a single core or a central unit with two functional groups and tri-functional branch points, comprising a further bead between adjacent branch points of the dendra so as to avoid bond crossing. Furthermore, the beads interacted through a hard-sphere potential with an appropriate bead diameter, which was indirectly related to the solvent quality. The starting conformations for the simulations were obtained through a lengthy procedure involving a preliminary atomistic Molecular Dynamics simulation using a standard force-field (the CVFF force field) for a few specific dendrimers [
45]. A careful ad hoc strategy was used in these runs carried out in vacuo, but eventually the modelled
reproduced the experimental value. In this way, a reasonable set of conformations was achieved, wherefrom a systematic coarse-grained sampling of relevant atoms (in particular the branch points) provided the starting coordinates of the beads for the Metropolis Monte Carlo runs. These simulations, which were much more efficient than the atomistic Molecular Dynamics runs for sampling the phase space of the system, were performed by small random displacements of a randomly selected bead, provided the resulting bond lengths were consistent with the distribution inferred from the atomistic Molecular Dynamics runs. These simulations yielded the conformational averages (for instance, those involved in the pre-averaged hydrodynamic interaction) required to calculate the intrinsic viscosity according to the lower bound of the variational procedure introduced by Fixman [
27], assuming that a frictional radius of each bead had been appropriately chosen. The intrinsic viscosity plotted as a function of
g did indeed show a maximum, which was quantitatively consistent with the experimental data [
45]. It should be noted, however, that the ratio between the bead friction radius and the bead hard-sphere radius required for dendrimers appears to be much larger than that required in linear chains, which suggests a different non-draining behavior. This point would require a more detailed investigation, in our opinion; in fact, with these dense systems, the assumption of the solvent treated as a continuum that pervades the inner part of the molecule may be incorrect, while a different treatment of a few discrete molecules would be required, with a physically different approach.
Subsequently, Freire et al. improved their model somewhat by applying a few corrections of varying importance [
46]. The first correction, and probably the most important one, was the addition of the contribution of an individual friction bead, which could be important for small molecules. This term amounted to assuming that the friction forces are distributed on the bead surface, rather than acting on the bead centers. The second correction amounted to incorporating a realistic distribution of the internal angle formed by the bonds connecting adjacent branch points. The third correction involved adopting a more realistic distribution of the distances between the adjacent branch points in order to better account for the presence of the solvent molecules. The simulation results again provided the observed maxima of [η] as a function of
g, and were generally in very good agreement with the experimental data. Moreover, the optimized bead friction radius was generally more reasonable than what was required by Reference [
45], even though its value, which was larger than that of the hard-sphere radius, was still not fully clear. Subsequently, these parameters were further revised [
47], but the relative size of these two radii was not qualitatively affected. It should also be noted that in a subsequent paper, Freire et al. [
48] also investigated the effect of pre-averaging on the calculated intrinsic viscosity obtained adopting the same model as in Reference [
46]. The pre-averaged approximation turned out to be inaccurate, especially at high
g, producing results that were too large by a factor greater than 2 at the highest values of
g (
g = 7).
In a later simulation study [
49], Brownian Dynamics simulations were used to investigate the dynamics of dendrimers adopting a finite-extensible nonlinear elastic (i.e., non-Hookean) FENE potential for the springs connecting the adjacent branch points. Excluded-volume interactions were ignored, while the hydrodynamic interaction was modelled through the Rotne-Prager tensor. Each spring was assumed to stand for a given number of statistical segments (or better, Kuhn steps), which were entered as a parameter in the expression of the FENE potential, so that different spacers between the branch points could be effectively modelled. The viscosity was then calculated as the ratio between the applied stress tensor and the shear rate
, so that [η] was obtained within the limit
→0. The intrinsic viscosity plotted vs. the generation was found to show a maximum at
g = 3, which is a bit less than, although still close to, what was observed experimentally. It is also interesting to point out that a significant effect of shear rate on viscosity was detected upon increasing
. In particular, the intrinsic viscosity first exhibited a minor decrease when
≈ 1 (here,
is a characteristic relaxation time derived from the zero-shear intrinsic viscosity), then a minor increase with a small shear thickening effect within a decade in the applied frequency, and finally a pronounced shear thinning effect at lager
with a power-law scaling behavior so that [η] ~
−α with α ≈ 0.55–0.6. Interestingly, the same pattern of shear thinning, followed by shear thickening and then by a much larger shear thinning was previously predicted for heavily branched star polymers with long arm lengths using a different approach [
50]. In the latter study, an analytical approach made it possible to describe the behavior of star polymers under shear, considering both the excluded-volume and the pre-averaged hydrodynamic interaction through the stochastic Langevin equation under the constraint of a constant contour length, somehow consistent with, but different in detail from, the approach of Reference [
49]. These changes of the intrinsic viscosity with increasing applied
, starting roughly at the same reduced shear rate
≈ 1, where
is the longest relaxation time of the chain under a very small
(→0), were attributed to the incipient arm unwinding, followed by a change in the draining regime due to the arm deformation, and eventually to a decrease of the hydrodynamic interaction due to the full arm deformation of the arms in the flow direction with an asymptotic power-law
−α dependence with α = 2/3 [
50]. Such physical effects are also at play in the case of dendrimers, as also briefly noted in Reference [
49].