Is the manifestation of the local dynamics in the spin – lattice NMR relaxation in dendrimers sensitive to excluded volume interactions ?

The effect of excluded volume (EV) interactions on the manifestation of the local dynamics in the spin– lattice NMR relaxation in dendrimers has been studied by using Brownian dynamics simulations. The study was motivated by the theory developed by Markelov et al., [J. Chem. Phys., 2014, 140, 244904] for a Gaussian dendrimer model without EV interactions. The theory connects the experimentally observed dependence of the spin–lattice relaxation rate 1/T1H on the location of NMR active groups with the restricted flexibility (semiflexibility) of dendrimers. Semiflexibility was introduced through the correlations between the orientations of different segments. However, these correlations exist even in flexible dendrimer models with EV interactions. We have simulated coarse-grained flexible and semiflexible dendrimer models with and without EV interactions. Every dendrimer segment consisted of two rigid bonds. Semiflexibility was introduced through a potential which restricts the fluctuations of angles between neighboring bonds but does not change orientational correlations in the EV model as compared to the flexible case. The frequency dependence of the reduced 1/T1H(oH) for segments and bonds belonging to different dendrimer shells was calculated. It was shown that the main effect of EV interactions consists of a much stronger contribution of the overall dendrimer rotation to the dynamics of dendrimer segments as compared to phantom models. After the exclusion of this contribution the manifestation of internal dynamics in spin–lattice NMR relaxation appears to be practically insensitive to EV interactions. For the flexible models, the position omax of the peak of the modified 1/T1H(oH) does not depend on the shell number. For semiflexible models, the maximum of 1/T1H(oH) for internal segments or bonds shifts to lower frequencies as compared to outer ones. The dependence of omax on the number of dendrimer shells appears to be universal for segments and bonds in dendrimer models with and without EV interactions.


Introduction
Dendrimers are regular hierarchically branched polymers with a wide range of potential applications such as catalysts or nanocarriers for drugs, dye and metal nanoparticle delivery, [1][2][3][4] etc.The dynamical properties of dendrimers are of great importance for many practical applications.Experimentally, dendrimer dynamics is studied by dielectric relaxation, nuclear magnetic resonance (NMR), birefringence measurements and polarized luminescence.In particular, the NMR spin-lattice relaxation method allows us to obtain information about the orientational mobility of NMR active groups by measuring the spin-lattice relaxation rate, 1/T1 H .In the case when 1 H is used as a marker, 1/T1 H is connected with the spectral densities J(o H ) and J(2o H ) by the relation where o H is the angular frequency of the NMR spectrometer used.A 0 is a constant, which is insensitive to temperature or frequency and is determined by the type of NMR active group.Spectral densities can be calculated from the second order orientational autocorrelation function (ACF) P 2 (t) for an internuclear vector r JðoÞ ¼ 2 Due to the sensitivity of NMR spectra to the chemical environment it is possible to measure 1/T1 H for groups located in different dendrimer shells, and therefore to monitor their mobilities separately.Experimental data for many dendrimers [5][6][7][8][9] show a remarkable difference between the mobilities of inner and terminal segments.Similar results were obtained earlier by MD simulation (see ref .10).In spite of a large amount of experimental studies of NMR relaxation in dendrimers, the number of theoretical investigations in this field is insufficient.In contrast to NMR experiments, where the temperature dependence of 1/T1 H at a fixed frequency, o H , is studied, the theory usually considers the frequency dependence of 1/T1 H (o H ) at constant temperature.But its predictions can be compared with the results of FC relaxometry which covers a rather large frequency range (see, for example ref. 11).
In the majority of theoretical papers on NMR relaxation in dendrimers it is assumed that the dynamics of the NMR active groups reflects the orientation mobility of the segments of the dendrimer.This means that the internuclear vectors in these groups have a component directed along the segment.The theory of NMR relaxation in flexible dendrimers was developed on the basis of the Gaussian model of dendrimers (ref.12 and 13).In this model, the segments are represented by Gaussian subchains and excluded volume interactions are neglected (phantom model).The reduced dimensionless spin-lattice relaxation rate was considered: If the relaxation of ACF P 2 (t) is described by the simple exponential function with the relaxation time t, then the frequency dependence of [1/T1 H ](o H ) has a maximum at However, for dendrimers the time dependence of ACF P 2 (t) is not described by a single exponent.Nevertheless, the frequency dependence of [1/T1 H (o H )] for the Gaussian model has one maximum.The most interesting result for this model is that the position of this maximum on the frequency scale does not depend on the segment topological location inside a dendrimer.It means that the characteristic times appearing in the NMR relaxation should be the same for all NMR active groups.However, these conclusions are in contradiction with experimental data for many dendrimers.
To explain this disagreement, it was assumed 13 that the experimentally observed dependence of the NMR signal on the segment position inside a dendrimer is connected with a restricted flexibility (semiflexibility) of the dendrimer.Indeed, calculations have shown that the incorporation of semiflexibility into the Gaussian model of dendrimers leads to the dependence of o max on the number of dendrimer layers in the labeled segment. 13he Gaussian dendrimer model is a phantom model, i.e. non-bonded interactions are neglected.But in real dendrimers, the contribution of these interactions is significant.As a result, their size and internal structure differ qualitatively from those predicted by the Gaussian model.
Semiflexibility was incorporated into the Gaussian model through the correlations between the orientations of segments located in different dendrimer shells.
where d i and d j are unit vectors directed along the i-th and j-th segments which are separated from each other by k segments.The parameter of semiflexibility in this model, q, is an average cosine of the angle between neighboring segments.For the model without correlations (q = 0), the maxima of [1/T1 H (o H )] for segments belonging to different shells coincide.For the model with correlated segments, the position of the maximum shifts to lower frequencies for segments belonging to more inner layers.However, as it was shown by Kumar and Biswas, 14,15 the neglect of excluded volume (EV) interactions can be compensated by incorporating semiflexibility into the model, where the restrictions imposed on the segment direction and orientation exactly reproduce the features of EV interactions.
Therefore, the correlations between the orientations of segments are inherent even for dendrimers consisting of flexible elements, and thus there are no real dendrimers without such correlations.Due to excluded volume interactions, there are always correlations between the orientations of segments.Nevertheless, our previous simulation of simple coarse-grained dendrimer models with excluded volume interactions 12 shows that the positions of the maximum of 1/T1 H coincide for segments located in different dendrimer shells.It means that in this respect the flexible Gaussian and EV models behave similarly.The question arises as to whether this similarity remains when we speak about the effect of semiflexibility on the position of the [1/T1 H (o H )] maximum.
It should be noted that the dependence of the mobility of NMR active groups in poly-lysine dendrimers established experimentally and confirmed by full-atomic molecular dynamics simulations in our recent work 9 was interpreted as a manifestation of the restricted flexibility of these dendrimers.
However, a direct comparison of the effect of semiflexibility on the segmental mobility in phantom and excluded volume dendrimer models has not been performed up to now.In the present work, we will perform such a comparison by using computer simulation of simple coarse-grained models.We will show that the shift of the maximum of [1/T1 H (o H )] occurs due to the introduction of an additional angular potential into the EV dendrimer model even if the value of the correlation parameter q remains practically the same.
The remainder of this paper is organized as follows.In Section 2, we briefly describe the model and simulation details.Section 3 is devoted to our results and their discussion.Eventually, Section 4 contains the conclusions.

Model and simulation method
We consider the coarse-grained model of dendrimers depicted in Fig. 1 with different numbers of generations and numbers of beads (Table 1).In the dendrimer models considered earlier in ref. 16 and 17 every segment connecting neighboring branching points consisted of only one bond.In our case, every segment contains two bonds.A similar dendrimer model was simulated in the papers of Karatasos et al. 18,19 This model allows monitoring the dendrimer dynamics not only on the segment scale but also on the scale of bonds.
The bonds are considered as springs with the potential where the value of the elasticity coefficient k bond = 1000.0kT.l 0 is the equilibrium bond length.Here and below, we will give all lengths in units l 0 and energies in units of the thermal energy k B T.
We consider two types of models with the topology shown in Fig. 1.The first type is a phantom (Ph) model, i.e. a model without the excluded volume interactions while the second model includes EV interactions between all beads: where e and s are the potential parameters.The presented form of the potential is different from the classical Lennard-Jones potential by term e.In such a form, the potential monotonically decreases to zero up to s Á ffiffi ffi 2 6 p .Thus, the potential U R is a purely repulsive potential that represents the short-ranged interactions between non-bonded monomers.It describes dendrimers under good solvent conditions and reproduces the main features of real dendrimers in solutions such as the molecular mass dependence of the gyration radius, density profile, terminal group distribution, etc.The potential parameters s = 0.8 and e = 1.0 are similar to those used in ref. 16.
Semiflexibility is introduced into the dendrimer model through the internal potential which restricts the angle g ij between neighboring bonds where the values of equilibrium angle g 0 = 109.5 and elasticity coefficient k ang = 10 were used.We intentionally have chosen such a potential so that the average cosine of the angle between neighboring segments will be equal to that in the flexible EV model.We will show that the introduction of such a potential will not change the mutual orientations of segments in the EV model remarkably in comparison with the flexible case.In the Ph model, this potential will induce correlations between the orientation of segments which are absent in the flexible Ph model.
In the previous work, 12 we used a free-draining dendrimer model where hydrodynamic interactions (HI) were neglected.As it was shown in ref. 12, the effect of HI interactions appears mainly in global motions, such as the rotation of the dendrimer as a whole.Their effect on internal modes is much weaker.As it was shown in our paper, 20 for a flexible dendrimer model with segments consisting of one bond the inclusion of HI does not change the conclusion about the independence of the position of the maximum of the reduced spin-lattice relaxation rate 1/T1 H (o H ) from the shell number.
The time evolution of particles in the BD method is described by coupled Langevin differential equations without an inertial term where x s the friction coefficient of particles, F(r i ) is the force affecting the i-th particle due to interactions with other particles and G is a random force where All results were obtained by averaging over the simulation time.In order to ensure that this time is sufficiently large we have studied the dynamics of terminal beads in our models.The density distribution of terminal beads relative to the center of mass (Fig. 4) shows that terminal groups are not only located at the dendrimer periphery but are also distributed along its volume.This well known effect is denoted in the literature as ''backfolding'' (see, for instance, ref. 22).This distribution is a dynamical one: terminal groups do not remain in definite positions but move to the center and back to the periphery.For an estimation of the time of this exchange we have calculated the mean squared displacement hDr 2 (t)i of terminal beads as a function of time (Fig. 5).hDr 2 (t)i increases initially with time and after some time t max reaches a plateau.Therefore, t max can be considered as the time of diffusion of terminal monomers from the inside to the outside of the dendrimer or backwards.In all cases, the time t max is much smaller than the simulation time t sim .Therefore, the dendrimer conformation changes many times during the simulation run and the t sim is sufficiently large.
The authors of ref.
12 considered the case when the internuclear vector of NMR active groups has a component directed along the segment and the orientational mobility of segments is displayed in the reorientation of these vectors.In our model, segments consist of two bonds (springs) (Fig. 1) and we will monitor not only the mobility of segments as a whole but also the mobility of constituent bonds.

Correlations between segment orientations
For the characterization of correlations between segments separated by s shells the average cosine between them was calculated: where d m is the end-to end vector for a segment belonging to the m-th dendrimer shell (Fig. 2).We calculated the number of shells from the outer (m = 0) to the inner layers.The dependence depicted in Fig. 3 can be presented in the scaling form The values of the correlation parameter q are shown in Table 2.
As expected, for the flexible Ph model q = 0, i.e. the orientations of segments are uncorrelated.The introduction of the potential (10) restricting the angle fluctuations between neighboring bonds leads to the appearance of correlations between the orientations of segments in this model.In the semiflexible phantom model, the value of the correlation parameter does not depend practically on the dendrimer generation number G.
But in the flexible EV model even without the angle potential (10) the correlation parameter q is rather large and its value increases with an increase of the generation number G from q = 0.2 for G = 2 to q = 0.5 for G = 5.The addition of the angle potential (10) to the EV model does not change the value of the correlation parameter q.
As a result, angular restrictions imposed by the angle potential (10) have only a minor effect on the equilibrium structural characteristics of a dendrimer such as the gyration radius and the core-terminal group distance as well as the density profile (see Appendix A).
Dendrimer rotation as a whole.In previous work, 12 only the contribution of internal modes to the orientational dynamics of segments was considered.However, in real dendrimers the dynamics of segments or bonds is determined not only by internal modes but also by the overall dendrimer rotation.The characteristic time of this rotation t rot (G) depends on the generation number G and can be estimated from the first order orientational ACF (see, for example ref.C rot (t) = he g (0)e g (t)i ( 16) for the unit vector Here, r t represents the position of an arbitrary terminal monomer and the vector r 0 represents the position of the dendrimer core.Averaging in ( 16) is performed over the simulation time and over all terminal monomers.Autocorrelation functions C rot for different dendrimer models are shown in Fig. 6.
Following the method described in ref. 12 we approximate C rot (t) using two exponents where the largest relaxation time is considered as t rot (G).The values of t rot (G) for dendrimer models are presented in the Table 3.
As expected, t rot (G) increases with an increase of G.This time is the largest characteristic time of a dendrimer of the generation G but it is much less than the simulation time t sim which is varied from 10 4 for G = 2 to 16 Â 10 4 for G = 5.This confirms again that the total simulation times used in these calculations are large enough to obtain the correct data.
The dependence of the orientational correlation time t rot on the overall number of beads N in a dendrimer is shown in Fig. 7.It can be presented in the scaling form The values of the exponent a are shown in Table 3.
Note that for the free-draining flexible EV model of a dendrimer with segments consisting of one bond the value of a = 1.6 was obtained earlier.This journal is © the Owner Societies 2016

Manifestation of the orientational dynamics of segments in NMR relaxation
The averaging in ( 19) is performed over the simulation time and over all segments belonging to the m-th shell.To compare the results obtained for different models the ACF P 2,s (t) for a single segment (not included in the dendrimer) consisting of three beads connected by two bonds was calculated P 2;s ðtÞ ¼ 3 2 Fig. 8 show the frequency dependence of [1/T1 H (o H )] m for flexible Ph and EV models for a single segment and for segments in shells with m = 0, 1 and 2 as well as the frequency dependence of [1/T1 H (o H )] m for a segment in the dendrimer with the ''frozen'' internal mobility.In this case, the ACF P 2,s (t) will be described by a simple exponent P 2 (t,G) with the relaxation time t rot ðGÞ 3 .
For the flexible Ph model (Fig. 8a) all curves have one maximum located practically at the same frequency for all m.Such a behavior is in agreement with predictions of the theory for the Gaussian dendrimer model. 13The position of this maximum is slightly shifted to lower frequencies as compared to that for a single segment.The curves of [1/T1 H (o H )] m become broader from the side of low frequencies by moving from m = 0 to m = 2.
In the flexible EV model, only the position of the maximum of [1/T1 H (o H )] 0 for terminal segments coincides with that for a single segment.The frequency dependence of [1/T1 H (o H )] m for internal segments not only becomes broader but its maxima shift to lower frequencies with an increase of m.
At first glance, such a behavior disagrees with that observed for the Ph model and with theoretical results for the Gaussian model. 13As we will show below, this difference is connected with a different contribution of global dendrimer motion to the segmental dynamics for these models.
The effect of the introduction of semiflexibility is the same for Ph and EV models (Fig. 8).The position of the maximum shifts to lower frequencies by moving from terminal to inner shells in agreement with the predictions of the Gaussian model.
Let us discuss now the origin of the difference between the results obtained for the flexible Ph and EV models.As mentioned above, the rate of reorientation of dendrimer segments is determined both by internal modes and by the overall dendrimer rotation.
Recall that in the theoretical work for the Gaussian model 13 only internal modes were taken into account.
For the separation of the contribution of the overall rotation we use the procedure proposed in previous work. 12g.6 The time dependence of ACF C rot for the Ph.(a) and EV (b) dendrimer models.
Table 3 The values of orientational correlation times for rotation of the dendrimer as a whole t rot (G) for different dendrimer models.Exponent a as fitted from We start consideration from the first order ACF: for segments belonging to the m-th shell in a dendrimer with the generation number G. The relaxation of P 1,s (t,G,m) is determined by internal modes and the overall dendrimer rotation.Following the method described in ref.
12 we assume that P 1,s (t,G,m) can be presented in the form: where the first term describes the contribution of internal modes and is the same for segments belonging to the shell with a number m in dendrimers with different G values.Functions P ˜1,s (t,m) and the factor a m for m = 0, 1 and 2 were obtained by fitting simulated P 1,s (t,G,m) for dendrimers with G = 2-5, 3-5 and 4-5, respectively (Table 4).
The table shows the values of a m for the models considered.As expected a m increases with m.For Ph models the contribution of the overall rotation to the orientational dynamics of segments is small in contrast to EV models, where this contribution reaches almost 50% for internal shells.We can introduce the function which characterizes the rate of reorientation of dendrimer segments only due to internal modes.Namely this function was considered in the work. 12For the calculation of the second order ACF P int 2,s (t,m) we use the relationship 23,24 Our calculations show that it is valid for segments both in flexible and semiflexible EV models.By using P int 2,s (t,m) the frequency dependence of [1/T1 H (o H )] m was calculated for both models (Fig. 9).The results depicted in Fig. 9 demonstrate not only a qualitative but also a quantitative agreement between the results obtained for the Ph and EV model.It means that the internal orientational dynamics of the segments which manifests in spin-lattice NMR relaxation is practically insensitive to excluded volume interactions.Returning to Fig. 9, we can conclude that the observed shift  In both models, the introduction of semiflexibility through the angular potential (10) leads to the shift of the maximum to lower frequencies by an increase of the shell number m in agreement with the theoretical predictions for the Gaussian dendrimer model. 13The absolute value of this shift is practically insensitive to excluded volume interactions if the contribution of only internal modes to the segment reorientation is considered.This value is determined by the angular potential (10) which restricts the fluctuations of angles between bonds, wherein the values of the correlation parameter q differ for Ph and EV models.

Orientational dynamics of bonds
Now we consider the situation where the internuclear vectors in NMR active groups are directed along bonds inside the dendrimer segments.Every segment consists of two bonds, and therefore every shell can be divided into two subshells numbering from the terminal s = 0 (Fig. 2).As in the case of segments we have calculated the frequency dependence of [1/T1 H (o H )] s for bonds belonging to a subshell with the number s in the dendrimer of generation G on the basis of the second order ACF P 2,b (t,G,s).
For comparison, we have also calculated the dependence of [1/T1 H (o H )] for a single bond (not included in the dendrimer) which consists of two beads connected by a rigid spring.In this case, the corresponding ACF P 2,b (t) is given by the expression: The analysis of the simulation results obtained for both models shows that the effect of the dendrimer generation number G on the frequency dependence of [1/T1 H (o H )] s for bonds is much weaker than for segments.For subshells s = 0-5 this effect disappears practically for dendrimers with generation number G Z 4. It means that the reorientation of bonds is determined mainly by internal modes.Therefore, it is not necessary to separate the contribution of the overall rotation as in the case of segments.Fig. 10 shows the calculated dependence of bonds in the models considered.In the same plot, we also show the dependence of o max,s (m) (reduced to o max,s (0)) for segments found from Fig. 11 where the effect of the overall dendrimer rotation was excluded.
We see that all points lie on a common line indicating a similar effect of semiflexibility on the internal orientational dynamics of segments and bonds which is insensitive to excluded volume interactions.Note that at the same potential (10) restricting the fluctuations of bond angles, the values of the correlation parameter q for Ph models are smaller than those for EV models (see Table 2).

Conclusions
The effect of excluded volume interactions on the manifestation of the local orientational dynamics of segments (or bonds) in spin-lattice NMR relaxation in dendrimers has been studied by using Brownian dynamics simulations.The study was stimulated by the theory developed in ref. 13 for the Gaussian dendrimer model without EV interactions.The theory predicts that for flexible dendrimers the position of the peak of the frequency dependence of the reduced spin-lattice relaxation rate [1/T1 H (o H )] does not depend on the proximity of dendrimer shells to the core by the assumption that the internuclear vectors have a component directed along segments.For semiflexible dendrimers, the peak of [1/T1 H (o H )] shifts to lower frequencies by moving from outer to inner segments.Semiflexibility was introduced into the theory 13 through the correlations between the orientations of different segments.However, these correlations exist even in flexible dendrimer models with EV interactions.
We have simulated coarse-grained flexible and semiflexible dendrimer models with generation number G = 2-5 with three functional branching points with and without EV interactions.Every dendrimer segment consisted of two rigid bonds.Semiflexibility was introduced through a potential which restricts the fluctuations of angles between neighboring bonds.The potential leads to the appearance of orientational correlations between segments in the Ph models but does not change these This journal is © the Owner Societies 2016 correlations in the EV models as compared to the flexible case.The frequency dependence of [1/T1 H (o H )] for segments and bonds belonging to different dendrimer shells was calculated.It was shown that the main effect of EV interactions consists of a much stronger contribution of the overall dendrimer rotation to the dynamics of dendrimer segments as compared to the phantom models.After the exclusion of this contribution the manifestation of internal segmental modes in spin-lattice NMR relaxation appears to be practically insensitive to EV interactions.For the flexible models the position o max of the peak of the modified [1/T1 H (o H )] does not depend on the shell number in spite of the orientational correlations between segments in the EV model.For semiflexible models the maximum of [1/T1 H (o H )] for internal segments shifts to lower frequencies as compared to the outer ones.The value of this shift is practically the same for Ph and EV models.For the orientational dynamics of bonds comprising the segments the contribution of the overall dendrimer rotation is much smaller than for segments even with the presence of EV interactions.The contribution of EV interactions is remarkable for dendrimers with segments consisting of one bond and becomes practically negligible for the reorientation of bonds inside two-bond segments.The average concentration of monomer units in a dendrimer decreases with the segment length l s as l s 4/5 at the same generation number under good solvent conditions.As a result, the effect of EV interactions will decrease with l s .It was shown in ref. 13 that the main contribution to the segmental orientational autocorrelation function (ACF) which is fundamental for the frequency dependent [1/T1 H (o H )] is made by local modes which are invariant to the shell number.The introduction of semiflexibility into this model leads to the shift of the maximum to lower frequencies for inner shells.It is explained by the increased contribution of modes with a larger scale (pulsation modes) to the ACF.We show that for flexible and semiflexible models the internal orientational dynamics of segments manifesting in spin-lattice NMR relaxation is not sensitive to the presence or absence of EV interactions.
Thus, our results give additional arguments for the interpretation 9 of the experimentally observed difference between the positions of [1/T1 H (o H )] maxima for terminal and inner segments in dendrimers as a manifestation of the restricted flexibility of these macromolecules.
The reported study was supported by the Supercomputing Center of Lomonosov Moscow State University. 25This work

¼ xl 0 2 k
a and G b are components of random force G.We use the simple algorithm proposed by Ermak and McCammon 21 to solve eqn (11) with a time unit equal to the characteristic time t 0 B T and time step Dt = 0.00005t 0 .The total simulation time was t sim from 2 Â 10 8 Dt = 10 4 t 0 (for the dendrimer of second generation, G = 2) to 32 Â 10 8 Dt = 16 Â 10 4 t 0 (G = 5).The relaxation time for a single bond (not included in the dendrimer) t b = t 0 /4.We measured the dynamic characteristics in t b .

Fig. 1
Fig. 1 Schematic model of a dendrimer of G = 3 generations and the core functionality 3. Dotted circles denote the boundaries of each generation layer numbered by m.Black, gray, and white circles (beads) represent the core monomer, branched monomers and linear monomers inside segments, respectively.The bonds (springs) connecting the beads are shown as solid lines. 12)

Fig. 2
Fig. 2 Schematic picture of a dendrimer branch.Vectors d m for segments belonging to the terminal (m = 0) and to the first two internal (m = 1, m = 2) shells (a) and vectors b s for bonds (springs) belonging to these shells (s = 0,. ..,5)(b) are shown.

Eqn ( 1 )
-(4) connect the frequency dependence of [1/T1 H (o H )] m or segment belonging to m-th shell in the dendrimer of the generation number G with the second order autocorrelation function (ACF) P 2,s (t,G,m) for this segment.

Fig. 4
Fig. 4 Density distribution of terminal beads relative to the center of mass.

Fig. 5
Fig.5The mean displacement hDr 2 (t)i of terminal beads as a function of time.

Fig. 3
Fig. 3 The dependence of mean hcos c(s)i for flexible and semiflexible Ph (a) and EV (b) models.

Fig. 9
Fig. 9 Frequency dependence of [1/T1 H (o H )] m calculated on the basis of (25) for segments belonging to different sub-shells for flexible (a) and semiflexible (b) Ph (dotted line) and EV (solid line) models.

[ 1 /
T1 H (o H )] s for bonds in flexible Ph and EV models for G = 5.It is seen that this dependence is very similar for flexible Ph and EV models.All curves of [1/T1 H (o H )] s corresponding to different subshells (from s = 0 to s = 5) have one maximum at o max,b (s).The values of o max,b (s) for different s values are close to each other and to the value of o max,b for a single bond.Similar to segments the introduction of the angle potential (10) for both models leads to the shift of the position of the maximum of [1/T1 H (o H )] s to lower frequencies by the increase of s i.e. by moving from terminal bonds to internal ones (Fig.10).Fig. 11 presents the dependence of the frequency o max,b (s) (reduced to o max,b (0)) corresponding to the maximum of [1/T1 H (o H )] s for

Fig. 11
Fig. 11 The dependence of the positions.o max,b (s)/o max,b(s=0) of the maximum of [1/T1 H (o H )] s for bonds and o max,s (m)/o max,s (m = 0) of the maximum of [1/T1 H (o H )] m for segments for flexible and semiflexible Ph and EV models.The figure consists of two abscissa axes: at the bottom for the bonds, at the top for the segments.The segment with m and the bond with s = 2m + 1 originate from the same branching point.

Fig. 12 2 1 2
Fig.12The square root of the mean-squared gyration radius R g 21 2 (a)and of the mean-squared distance R 2 1 2 between the core and terminal beads (b) for the dendrimer models considered as a function of the number N of beads.

Fig. 13
Fig. 13 Density profiles for flexible and semiflexible EV dendrimer models.

Table 1
Characteristics of simulated dendrimers.G -number of generations, N -number of beads.N t -number of terminal beads

Table 2
The values of the correlation parameter q 5 a Fig. 7 The dependence of the orientational correlation time t rot of the dendrimer rotation as a whole on the overall number of beads N for the models considered.This journal is © the Owner Societies 2016 Phys.Chem.Chem.Phys., 2016, 18, 24307--24317 | 24313

Table 4
Values of a m characterizing the contribution of the overall rotation for the dendrimer models considered T1 H (o H )] m to lower frequencies for inner segments in the flexible EV model is connected with a much stronger contribution of the overall dendrimer rotation to the reorientation of inner segments as compacted in the Ph model.