Dynamics of internally functionalized dendrimers

The internally functionalized dendrimers are novel polymers that differ from conventional dendrimers by having additional functional units which do not branch out further. We investigate the dynamics of these structures with the inclusion of local semiflexibility and analyze their eigenmodes. The functionalized units clearly manifest themselves leading to a group of eigenvalues which are not present for homogeneous dendrimers. This part of the spectrum reveals itself in the local relaxation, leading to a corresponding process in the imaginary part of the complex dielectric susceptibility.


Introduction
Dendrimers are hyperbranched polymers with a very high symmetry.Due to their perfectly branched structure, dendrimers are promising materials in the fields of chemistry, biology, and medicine, see e.g., ref. 1-3.While the conventional dendrimers were synthesized several decades ago, [4][5][6][7] recently new types of dendrimers have been created, see ref. 8 for review.][11][12][13][14] Their structure differs from homogeneous dendrimers: the internally functionalized dendrimers possess additional functional groups at each branching point (apart from the core).These functional units have distinct chemical properties so that they do not branch out further.As it has been highlighted in ref. 10, there are two major pathways to achieve functionalized dendrimers: premodification (i.e., functionalization before the construction of the macromolecule) or postmodification (i.e., functionalization after the construction of the macromolecule).With respect to these strategies one uses orthogonal protecting groups or groups that can be selectively activated. 10In particular, as has been recently shown, 13,14 the use of orthogonal click reactions based on epoxy-amine and thiol-ene chemistry is very efficient for the introduction of functional groups at each dendritic layer.The advantage of these compounds is that the attachment of drugs or imaging agents to the internal groups can lead to good solubility and biocompatibility properties in comparison to conventional dendrimers. 15Moreover, internally functionalized dendrimers can be used as multichromophoric light-harvesting systems, in which the internal functionalization helps to establish complex energy gradients. 11Furthermore, the layers with different polarities allow the internally functionalized dendrimers to operate as unimolecular catalytic nanosytems. 9n this paper we investigate the dynamics of internally functionalized dendrimers.For their description we use the model of semiflexible treelike polymers (STPs) which includes semiflexibility by restrictions on the bonds' orientation. 167][28] However, the STP framework leads to a higher computational effort stemming from the next-nearest neighboring interactions in the dynamical matrix.Nevertheless, due to the high symmetry of dendrimers, one can block-diagonalize the dynamical matrices very effectively as it has been shown for regular dendrimers. 29But, in contrast to regular semiflexible dendrimers, 29 for which the set of eigenmodes of ref. 30 can be readily transferred, for internally functionalized dendrimers the set of eigenmodes of ref. 30 is not straightforwardly applicable.The reason for this feature lies in the symmetry breaking due to the functional units, for which the corresponding degrees of freedom have to be consecutively included.These fundamental degrees of freedom for internally functionalized dendrimers lead to an additional group of eigenvalues, which is not present for regular dendrimers.The ensuing relaxation spectrum leaves its fingerprints on the dynamic properties of the macromolecules.Here we focus on the mechanical and dielectric relaxation, and investigate, in particular, the role of the functional beads.
The paper is structured as follows: in Section 2 we recall briefly the STP-model 16 and its formulas for dynamic quantities.In Section 3 we apply the STP-model to the internally functionalized dendrimers.Then, we introduce a complete set of eigenmodes which is used for the diagonalization of dynamic matrices.(The recursive construction of the resulting reduced matrices is relegated to Appendix A).We close Section 3 by analysis of the dynamics of internally functionalized dendrimers.Finally, Section 4 presents an overview of the most important results.

The model
The dynamics of internally functionalized dendrimers is modeled in the STP framework, 16 which we recall here briefly.
The STP model describes polymers as treelike structures of N beads that are connected by springs.Moreover, the bond variables (i.e.springs) {d a } are constrained to have: (i) fixed mean-square lengths, hd a Ád a i = l 2 ; (ii) two adjacent bonds d a and d b , connected by bead i, fulfilling hd a Ád b i = AEl 2 q i , where q i denotes the stiffness parameter of bead i, the plus sign applies for a head-to-tail orientation of As has been shown in ref. 16, the above listed restrictions (i)-(iii) on bonds {d a } lead to the potential energy V({d a }), where K is the spring constant that is related to temperature T and the mean-squared bond length l 2 by K = 3k B T/l 2 .In eqn (1)  the matrix W = (W ab ) is very sparse; its nonzero elements are either diagonal or related to adjacent bonds.It was shown in ref. 16 that the restriction (iii) follows from maximizing entropy under constraints (i) and (ii) and therefore no Lagrange multipliers associated with restriction (iii) are needed, i.e., the corresponding entries of W are zeros.We note that the vectors {r i } labeling the beads' positions are related to the bonds' vectors {d a } through where G is the so-called incidence matrix 31 and T denotes the transposition.Each line of G T is related to a bond and contains only two nonzero entries (+1 and À1) on the places related to the beads connected by the bond. 31Transforming eqn (1) to bead variables (using eqn (2)) yields the matrix A STP GWG T .Given that W couples only adjacent bonds and that the nonvanishing elements of G are related only to beads directly connected by a bond, the nonzero elements of A STP are either diagonal or related to the nearest neighboring (NN) and to the next-nearest neighboring (NNN) beads, vide infra.
The dynamics of a polymer is described by a set of 3N Langevin equations.The equation for the x-component of vector r i labeling the position of the ith bead reads Here t 0 = z/K, where z is the friction coefficient and K is the spring constant.Moreover, f ˜i(t) is the x-component of the stochastic force obeying white noise relations.Furthermore, the structure of a polymer (connectivity and local stiffness) is coded through the matrix A STP = (A STP ij ), whose elements are known in a closed form. 16To introduce them, Fig. 1 shows the NN and NNN of a bead i of a treelike structure.The NN beads of i are denoted by i k and the NNN ones by i k s .With this notation, the elements of the matrix A STP read as follows: (4) and Here, f i and f i k or q i and q i k denote the functionality or the stiffness of the beads i and i k , respectively.All other elements of the matrix A STP vanish.

Dynamic quantities
The properties of the matrix A STP are fundamental for the dynamics of STP. 32,33Using the eigenvalue spectrum of the matrix A STP one can compute the complex shear modulus G*(o) = G 0 (o) + iG 00 (o). 34G 0 (o) denotes the storage modulus and G 00 (o) the loss modulus.Using dimensionless variables one can obtain the reduced moduli 33 and The sum includes all non-vanishing eigenvalues {l j }.The eigenvalue l 1 = 0 and its corresponding eigenvector describe the translation of the whole macromolecule.
Using additionally the eigenvectors of the matrix A STP one can compute further quantities like the autocorrelation function M a 1 (t) hd a (t)Ád a (0)i/l 2 (9)   which turns out to read 27,35 M a 1 ðtÞ ¼ Here the matrix Q is constructed from the eigenvectors and the sum runs over all non-vanishing eigenvalues.As we proceed to show in Section 3.2 the eigenvalues {l j } are usually very degenerate for symmetric hyperbranched structures.Denoting the amount of different nonvanishing eigenvalues by w(G), eqn (10)  can be rewritten as with where the sum runs over all eigenvectors belonging to the eigenvalue l k .Using this coefficient one can determine through the Fourier-Laplace transform the complex dielectric susceptibility e*(o) = e 0 (o) À ie 00 (o) for any segment d a . 36The real and complex parts read 37 e 0 ðoÞ ¼ and respectively, where the sums run over the distinct nonvanishing eigenvalues.
3 Results and discussion

Analysis of A STP for SFDs
In this subsection we analyze the matrix A STP for semiflexible, internally functionalized dendrimers (SFDs).First of all, we recall the topology of these structures.The internally functionalized dendrimers differ from homogenous dendrimers by having additional functional beads which do not branch out further.The construction of a functionalized dendrimer is analogous to the homogeneous ones wherein the core bead is connected to f C beads and all the other beads to f beads.Therefore, the difference r = f À f C gives the number of functional beads attached to each bead except the core.Thus, setting f C = f leads to a homogeneous dendrimer of functionality f C .The number of beads of an SFD is hence given by the number of beads of a homogeneous dendrimer 38 plus the number of functional beads, Fig. 2 shows the topology of an internally functionalized dendrimer of functionalities f C = 3 and f = 4 and generation G = 3.
Based on the SFD structure, the matrix elements of A STP , which for an SFD we will call in the following by A SFD , can be readily constructed.We assume a homogeneous situation which means that all internal beads except the core have the stiffness value q.The core being different from all other beads has a stiffness value q C .First, we present the diagonal elements {A SFD ii }, then the NN elements A SFD For the diagonal elements one obtains eight different situations (for visualization, see Fig. 3, where the open circles mark the beads corresponding to the discussed matrix element): (a) If i is a peripheral or functional bead, one has f i = 1.Moreover, i has exactly one neighbor with functionality f (for G Z 2).Using eqn (4), the element A SFD ii , denoted by m 1 , equals (b) For G = 2, a bead in the first shell has functionality f and possesses ( f À 1) peripheral neighbors.The core is its inner neighbor with functionality f C .With this, the matrix element reads (c) For G 4 2, a bead adjacent to a peripheral bead has functionality f.Its ( f À 1) neighbors are peripheral and the remaining one inner neighbor has functionality f.It follows (d) A bead adjacent to the core has functionality f for G 4 2. Hence, one of its neighbors (namely, the core) has functionality f C , r neighbors have functionality 1 and ( f C À 1) neighbors have functionality f.Therefore, one has (e) A bead well inside the dendrimer of G 4 2 has functionality f and possesses f C neighbors with functionality f and r neighbors with functionality 1.Therefore, the matrix element reads (f) The core has functionality f C and (for G Z 2) f C neighbors with functionality f.This leads to (g) For G = 1, the core has f C neighbors with functionality 1, so the matrix element reads (h) For G = 1, a peripheral bead has one neighbor with functionality f C , therefore, one has Next, we turn to the NN elements, the corresponding interactions are depicted in Fig. 4 by wavy lines.There are four different elements:(i) If one of the beads is peripheral or functional, the other bead has functionality f (for G Z 2).
Eqn (5) then leads to the matrix element, denoted by n 1 , ( j) If one of the beads is the core, the other bead has functionality f (if G Z 2).Therefore, one gets (k) If both beads are inside the dendrimer, but none of them is the core, they have functionality f, which leads to (l) For G = 1, the element related to the core and to a peripheral bead reads The NNN elements A SFD ii ks depend only on the bead between the two considered beads, see eqn (6).Therefore, there are two possible distinct elements (the corresponding situations are exemplarily shown in Fig. 4): (m) If the middle bead is the core, eqn (6) leads to (n) If the middle bead is an inner bead with functionality f, one has All other elements of the matrix A SFD vanish.

Eigenmodes
In this subsection, we present a complete set of eigenvectors of the matrix A SFD .[40] Fig. 3 Schematic drawing represents the diagonal elements of the matrix A SFD .The core is depicted in orange, the other internal beads in black, the functional beads in red, and the peripheral beads in blue.Open circles mark the beads which are related to the elements {m i }.
Fig. 4 Schematic drawing of the nonvanishing nondiagonal elements of the matrix A SFD , where the corresponding interactions are indicated by wavy lines.The color scheme is as in Fig. 3.
This journal is © the Owner Societies 2016 In ref. 29 it was shown that the same construction of eigenvectors is valid also in the semiflexible case.Here, we extend these methods to examine SFDs of functionalities f and f C and generation G.
Using the symmetry and the SFD structure, the set of eigenvectors can be divided into (G + 1) groups.The first group consists of all eigenmodes in which only peripheral beads are moving.The second group consists of the eigenmodes in which the peripheral and their NN are moving.This scheme continues until the group (G + 1) is reached in which all the beads including the core are moving, see Fig. 6 and 7 for an illustration.
To be able to describe the eigenmotions of the SFD more accurately, we introduce the notation of the subwedge.A subwedge is a part of the dendrimer which starts with another bead rather than the core.Fig. 5 depicts subwedges Z (GZ) of different generations G Z .
As discussed in Section 2.1, the dynamics of SFDs is given by a set of Langevin equations, eqn (3).The corresponding set of homogeneous differential equations reads A SFD lk x k ðtÞ ¼ 0 for all l: By choosing convenient eigenmodes, the system of equations can be simplified to a few nontrivial equations.The first group of eigenmodes consists of eigenvectors which describe the movement of two NNN peripheral beads, see Fig. 6.Denoting the two moving beads as i and j, the amplitudes x l in eqn (30) read x l = 0 for i a l a j.
Here, x 1 labels the single nonvanishing amplitude.Using the respective matrix elements, the single nontrivial equation of motion reads All other equations of motion are trivial.For the NN and NNN of the beads i and j, the terms of the two moving beads cancel each other because they move with the same amplitude but in the opposite direction.This is sufficient to leave the rest of the dendrimers in rest as well.Therefore, there is one eigenvalue belonging to the first group of eigenmodes, denoted by l (1) .The described motion involves any two NNN peripheral beads.Each subwedge Z (2) contains ( f À 1) peripheral beads which leads to ( f À 2) different, linear independent eigenvectors.Given that the number of Z (2) subwedges is equal to f C ( f C À 1) GÀ2 , the eigenvalue l (1) for r r 1 is f C ( f C À 1) GÀ2 ( f À 2)-fold degenerate.
In the case of generation G = 1 one has to use different matrix elements which leads to the equation of motion Furthermore, having r Z 2, the movement of two functional NN beads (see Fig. 6(b)) exactly leads to the same equation of motion as eqn (33) and therefore to the same eigenvalue.Hence, the eigenvalue l (1) is additionally Altogether, the degeneracy of the eigenvalue connected to the first group is The second group consists of eigenvectors which describe the motion of two subwedges Z (2) against each other, see Fig. 7(a).Because of symmetry, all topologically equivalent beads of one subwedge move with the same amplitude, the respective beads of the other subwedge with the same amplitude in the opposite direction.There are two eigenvalues belonging to this group, denoted by l (2) 1 and l (2) 2 .Following an analogous discussion as above, the degeneracy of each of the eigenvalue is  The third group of eigenmodes describes the motion of two subwedges Z (3) , see Fig. 7(b).Here appears the first structural difference to the homogeneous dendrimer 29 because the subwedge Z (3) includes r nonperipheral functional beads.These beads are described by another matrix element than the internal beads in the same shell.Therefore, an extra equation of motion is required.This means that the change in the topology leads to an extra degree of freedom.Hence, there are four eigenvalues connected to the third group of eigenmodes instead of three as it is in the case of homogeneous dendrimers.The respective degeneracy is given by Now, one additional degree of freedom appears in each iteration from one group to the next one.That means, continuing the procedure leads to two new degrees of freedom for each transition, one for the internal beads of the new involved shell and one for its functionalized beads.Therefore, to the n-th group of eigenmodes (n = 2,. ..,G) belongs 2(n À 1) eigenvalues.The degeneracy is given by for n = 2,. .., (G À 1) and by D G = f C À 1 for the G-th group.
In the group (G + 1) all beads including the core are moving.All topologically equivalent beads move with the same amplitude which leads to (2G À 1) degrees of freedom and therefore to (2G À 1) nondegenerate eigenvalues.
To verify that this set of eigenmodes is complete one has to show that their number is equal to the number of beads N.
The number of eigenmodes is given by Here we used eqn (15) giving the number of beads of an SFD.

Reduced matrices
The knowledge of eigenmodes presented in Section 3.2 allows formulating the respective reduced Langevin equations for each group.The eigenvalues can be then determined by diagonalizing the corresponding reduced matrices.
In the first group of eigenmodes, two peripheral NNN beads move against each other.The Langevin equation is given by Hence, the respective eigenvalue reads for G Z 2. For G = 1, q has to be replaced by q C .The second group leads to the system of equations (G Z 3) For G = 2, the second equation of set ( 42) changes to The coefficient matrix for these two cases (G Z 3 and G = 2) reads Here and in the following, the first index in the last matrix entry has to be used in the case n o G and the second index in the case n = G, where n denotes the discussed group of eigenmodes.
For the third group of eigenmodes, one gets the reduced matrix Here, different from homogeneous dendrimers, one observes the additional row and column, as discussed in Section 3.2.
By following this procedure one can construct iteratively the matrices up to group G.These matrices are relegated to the Appendix.
Next, we discuss the group (G + 1).In this group, all beads including the core are mobile.All beads which are topologically equivalent move with the same amplitude.For a dendrimer of generation G = 1 the system of equations of motion for the second group reads Therefore, the corresponding reduced matrix is given by For generations two and three one gets and The B n -matrices for n 4 3 can be constructed iteratively from the matrices of an SFD of the previous generations.As for A n , these matrices can be found in the Appendix.

Dynamics of SFDs
We start the discussion of dynamic properties of an SFD by the analysis of the eigenvalue spectrum of A SFD .Using the reduced matrices of Section 3.3, one can easily compute the eigenvalue spectrum of a very large SFD.
Fig. 8(a) shows the eigenvalue spectra of an SFD of generation G = 10.The core's functionality is f C = 3, its stiffness parameter is q C = 0.49.The functionality of the inner beads is f = f C + r, their stiffness . Such choice of q ensures that both stiffness parameters reach the two limiting cases at the same time. 33,41,42One can observe in Fig. 8(a) that the eigenvalue 1 1 þ q takes a larger part of the spectrum with increasing amount of functional beads.This is obvious since this eigenvalue describes the motion of two peripheral beads against each other.The number of such beads increases for growing r.On the left-hand side of this plateau there is a difference between the cases r = 0 and r a 0. While the plateau of eigenvalue 1 1 þ q for r = 0 is preceded by a jump from another (smaller) plateau, for r a 0, the transition to this plateau is rather steady.The corresponding in-between eigenvalues (which do not exist for the homogeneous dendrimer, r = 0) originate from the additional degrees of freedom due to the functional beads.As discussed in Sections 3.2 and 3.3 these degrees of freedom appear starting from the third group, i.e. they exist for dendrimers of G Z 3.
Based on the spectra of Fig. 8(a) we compute the reduced loss moduli [G 00 (o)] which we show in Fig. 8(b).(We focus on the loss modulus because it typically displays a richer structure than the corresponding storage modulus. 33) As can be inferred from Fig. 8(b), the increasing influence of the eigenvalue 1 1 þ q manifests itself in increasing significance of the left maximum.The growth of the jump after this eigenvalue with a larger amount of functional beads squeezes the region of higher eigenvalues; consequently leading to a shift of the right maximum to higher frequencies.Now, holding the functionality of the core fixed and varying the amount of functional beads leads to a different functionality of all other beads.Thus, it is also interesting to compare a homogeneous dendrimer of functionality f = 4 (r = 0) and stiffness parameter q = 0.32 with an SFD of functionalities f C = 3 and f = 4 and stiffness parameters q C = 0.48 and q = 0.32.Therefore, we have a comparable setup concerning the functionalities and stiffness, while the number of beads of the two dendrimers differs.Fig. 9(a) shows the eigenvalue spectra of the considered dendrimers (and the corresponding SFD spectrum, rescaled with the number of beads).One can see that the spectra look similar after rescaling.The significant difference is the absence of a jump for SFDs in front of the longest plateau as discussed above.There are also minor differences in some of the eigenvalues and in their degeneracy, but they are rather insignificant.Therefore, these two spectra lead to very similar loss moduli as shown in Fig. 9(b).This feature can be traced back to the fact that the periphery of the dendrimers contains a large part of the molecular mass, so that it influences very much the mechanical relaxation, what has been also observed for heterogeneous dendrimers. 43ince one can find a homogeneous counterpart to SFDs which lead to similar mechanical relaxation moduli, we are going to look at local properties of such counterparts.The relaxation of single segments is described by the correlation function M a 1 (t) (see eqn (10)).This quantity requires the knowledge of coefficients involving eigenvectors.Although one can do rigorous calculations of these coefficients by means of projection operator 44 or reduced description 42 techniques, here for investigation of the local properties it is enough to perform brute-force computations for an SFD of a smaller generation.Hence in the following we discuss dendrimers of generation G = 5 and the setup of Fig. 9.
Fig. 10(a) shows the autocorrelation function M a 1 for different segments of a homogeneous dendrimer and of an SFD as a function of time.Due to the symmetry, for homogeneous dendrimers the segments belonging to the same shell g have the same dynamics.For SFDs there are two types of segments (functionalized and non-functionalized) in each shell g 4 1, see Fig. 2. In order to have a direct correspondence to the homogeneous dendrimer, we consider for the SFD non-functionalized segments only.One can see that the relaxation of the segments of the homogeneous dendrimer is slower than that of the corresponding ones of the SFD.The difference is smaller for the more peripheral segments.To understand this behaviour we look at the amplitudes C a k of the relaxation decays given by eqn (12) for different relaxation times t k = t 0 /l k , see Fig. 10(b).As is known, 27 homogeneous semiflexible dendrimers have one dominating coefficient C a k for each shell.The same behaviour we find for the SFD, as can be observed in Fig. 10(b).However, in the case of SFDs the dominating coefficient has a smaller influence than for the homogeneous dendrimer.Moreover, for SFDs there are some additional amplitudes in the region between t 0 and 10t 0 .These amplitudes arise from the eigenvalues in front of the plateau 1 1 þ q in the eigenvalue spectrum, as it was observed in Fig. 8. Furthermore, for each shell g 4 1 the most significant relaxation time of SFDs is smaller than that of the homogeneous dendrimer.These two facts lead to faster decaying correlation functions for higher g.Now, using the amplitude C a k one can compute the real part e 0 (o) and the imaginary part e 00 (o) of the complex dielectric susceptibility, see eqn (13) and ( 14).Here we focus on the e 00 (o) since it displays more emphasized differences.Fig. 11 shows the e 00 (o) computed based on the amplitudes of Fig. 10(b).As can be inferred from the figure, e 00 (o) is dominated by two processes.The right maximum originates from the large eigenvalues, i.e. from the relaxation times t k o t 0 .This maximum is less pronounced since the corresponding amplitudes C a k are small in this region.The left maximum originates mostly from the dominating, very large amplitude.Moreover, for lower g the left maximum appears for the SFD at a bit higher frequencies than for the corresponding shells of the homogeneous dendrimer.This indicates that the viscosity and the gyration radius of the SFD is a bit lower than those of the corresponding homogeneous dendrimer, because the behavior of these characteristics is dominated by lower eigenvalues, see the ESI † for details.Now, for functionalized dendrimers there is an additional shoulder in the intermediate frequency region, which is significantly developed especially for internal shells (i.e. for lower g).This additional process in the e 00 (o) arises from the relaxation times between t 0 and 10t 0 , see Fig. 10(b).Thus the presence of functional beads is clearly manifested in the local dynamics of the non-functionalized segments on intermediate frequencies.We note that for fully-flexible structures the presence of the additional process in e 00 (o) is hardly pronounced, i.e., the semiflexibility highlights the presence of functionalized units, see the ESI † for details.We close this section by discussing the role of the dendrimers' generation G.In Fig. 12 we display the autocorrelation function M a 1 (t) and the corresponding imaginary part of the complex dielectric susceptibility e 00 (o) for different shells g of an SFD of generations G = 5 and G = 7.As it was observed for homogeneous dendrimers, 27 for higher g the curves overlap each other: the corresponding segments feel the size of the branches which originated from them, but not the total size of the dendrimer, i.e., G.The behavior of the core segments is somewhat different.Here the difference between the core beads and all other beads is also revealed.In particular, the similarity of neighborhood for core segments of both dendrimers (G = 5 and G = 7) is reflected in the initial decay of M a 1 (t) or at higher frequencies of e 00 (o).

Conclusions
In this work we studied the dynamics of semiflexible, internally functionalized dendrimers (SFDs).Based on their symmetry we have block-diagonalized the dynamical matrix stemming from the STP framework.The ensuing reduced matrices grow linearly with the generation G of the dendrimer, whose molecular mass has an exponential dependence on G.The reduced scheme provides not only the computational effort, but also gives a better understanding   of the SFD's eigenmotions.Starting from generation G = 3 the SFD eigenmodes become distinct from those of a homogeneous dendrimer.These differences come from the functional beads inside the SFD, resulting in the group of new relaxation times close to the segmental relaxation time.While the mechanical relaxation of SFDs is mostly determined by the peripheral shells, and hence one can find for each SFD a homogeneous counterpart with a similar behavior, the local dynamic properties allow uncovering the presence of functional units.These lead to an additional process (related to the group of new relaxation times close to the segmental relaxation time) in the imaginary part of the complex dielectric susceptibility in the intermediate frequency region.

A Appendix
Here we present reduced matrices A n and B n for n 4 3.
The fourth and fifth groups of eigenmodes lead to the matrices and  where the matrix A ˜nÀ1 is obtained by crossing out the last row and the last column of the matrix A nÀ1 .
This journal is © the Owner Societies 2016 The B-matrices of generations four and five read and where the matrix B ˜nÀ1 is obtained by crossing out the last row and the last column of the matrix B nÀ1 and replacing its last entry (i.e.(B nÀ1 ) 2nÀ4,2nÀ4 ) by m 5 + ( f C À 2)r 2 .
d a and d b and the minus sign for all other configurations of d a and d b ; (iii) for two nonadjacent bonds d a and d c , connected by the path (b 1 , b 2 ,. .., b kÀ1 , b k ), one has hd a Ád a i = hd a Ád b 1 ihd b 1 Ád b 2 iÁ Á Áhd b k Ád c il À2k .

Fig. 1
Fig.1Nearest (e.g., i k ) and next-nearest (e.g., i k s ) neighbors of a bead i of a treelike network.

Fig. 2
Fig. 2 Schematic drawing of an internally functionalized dendrimer of functionalities f C = 3 and f = 4 and generation G = 3.The different shells g = 1, g = 2, and g = 3 are indicated by dashed circles.The functional beads are depicted in red.

Fig. 6
Fig. 6 (a) A sketch of a possible motion related to the first group.(b) A sketch of a possible motion related to the first group involving two functional beads of an SFD of functionalities f C = 3 and f = 5.The moving beads i and j are highlighted by color.See text for details.

Fig. 7
Fig. 7 (a) A sketch of a possible motion related to the second group.(b) A sketch of a possible motion related to the third group.The involved subwedges are colored in green and orange.See text for details.

Fig. 9
Fig. 9 (a) Eigenvalue spectra of a homogeneous dendrimer (f = 4, q = 0.32) and an SFD (f C = 3, q C = 0.48, f = 4, q = 0.32) of generation G = 10.The green line represents the rescaled spectrum of the SFD.(b) Reduced loss moduli computed based on the spectra of figure (a).See text for details.

Fig. 8
Fig. 8 (a) Eigenvalue spectra and (b) the corresponding reduced loss moduli of an SFD of generation G = 10 and functionality f C = 3 with different amounts of functional beads, r = f À f C , see text for details.In (a) the highlighted area indicates the eigenvalues of an SFD (r Z 1) that originate from the additional degrees of freedom due to the functional beads.

Fig. 11
Fig. 11 Imaginary part of the complex dielectric susceptibility corresponding to the structures of Fig. 10.The lines and the color code are the same as in Fig. 10.

Fig. 12 (
Fig. 12 (a) Autocorrelation function M a 1 (t) and (b) the corresponding imaginary part of the complex dielectric susceptibility e 00 (o) for functionalized (f C = 3, f = 4) dendrimers of generations G = 5 (dashed lines) and G = 7 (solid lines).The dendrimers' shells are numbered by g in an ascending order, starting from the core shell, as indicated in Fig. 2.

Fig. 10 (
Fig. 10 (a) Autocorrelation function M a 1 (t) and (b) the corresponding amplitude C a k of the relaxation decays of different segments of a homogeneous (f = 4) and a functionalized (f C = 3, f = 4) dendrimer of generation G = 5.In (a) the solid lines are related to segments of the homogeneous dendrimer, while the dashed lines represent the results of the corresponding non-functional segments of the SFD.In (b) the dots indicate the amplitude C a k for the homogeneous dendrimer and the triangles for the SFD.The dendrimers' shells are numbered by g in an ascending order, starting from the core shell, as indicated in Fig. 2.
higher groups the matrices can be constructed iteratively.The matrix A n is a 2(n À 1) Â 2(n À 1)-matrix and it reads higher generations, the matrix can be given in an iterative form.The matrix B n (n = 6,. .., G) is of the dimension (2n À 1) Â (2n À 1) and reads