Next Article in Journal
Preparation of the Carbonized Zif−8@PAN Nanofiber Membrane for Cadmium Ion Adsorption
Next Article in Special Issue
Influence of Raster Pattern on Residual Stress and Part Distortion in FDM of Semi-Crystalline Polymers: A Simulation Study
Previous Article in Journal
Facile Electrodeposition of NiCo2O4 Nanosheets on Porous Carbonized Wood for Wood-Derived Asymmetric Supercapacitors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Star Polymers vs. Dendrimers: Studies of the Synthesis Based on Computer Simulations

by
Piotr Polanowski
1,
Krzysztof Hałagan
1 and
Andrzej Sikorski
2,*
1
Department of Molecular Physics, Lodz University of Technology, 90-924 Lodz, Poland
2
Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
*
Author to whom correspondence should be addressed.
Polymers 2022, 14(13), 2522; https://doi.org/10.3390/polym14132522
Submission received: 20 May 2022 / Revised: 18 June 2022 / Accepted: 19 June 2022 / Published: 21 June 2022
(This article belongs to the Special Issue Multiscale Simulation and Modeling in Polymers)

Abstract

:
A generic model was developed for studies of the polymerization process of regular branched macromolecules. Monte Carlo simulations were performed employing the Dynamic Lattice Liquid algorithm to study this process. A core-first methodology was used in a living polymerization of stars with up to 32 arms, and dendrimers consisted of 4-functional segments. The kinetics of the synthesis process for stars with different numbers of branches and dendrimers was compared. The size and structure of star-branched polymers and dendrimers during the synthesis were studied. The influence of the functionality of well-defined cores on the structure and on the dispersity of the system was also examined. The differences in the kinetics in the formation of both architectures, as well as changes to their structures, were described and discussed.

1. Introduction

Star-branched macromolecules are of interest mainly because their size is smaller than that of linear polymers of the same molecular masses, and their viscosity behaves differently from linear chains, as it depends on the arm mass exponentially [1,2,3,4,5,6]. There are two main groups of synthetic methods for star synthesis: core-first and arms-first. In the first method, a multifunctional core is prepared in the first stage. Synthesis can be carried out using a large amount of cross-linker, and a large amount of the high conversion monomer is added to start growth arms [7,8,9,10,11,12]. If the reactivity of the cross-linker is higher than that of the monomer, the cores are synthesized at the beginning, and the rest of the monomer forms the arms of the star [13]. In the arms-first method, the linear polymers are synthesized first, and then a cross-linker is added and star cores are formed. The effectiveness of the cross-linking process depends mainly on the concentration of the cross-linker and solvent molecules [7,11,14,15,16,17]. Experiments [7,14,18,19,20] and simulations [21,22] have shown that a considerable fraction of the cross-linker also forms intrachain bonds which do not contribute to star core formation.
The theoretical interest in star-branched polymers is also significant [23]. Computer simulations of the properties of star-branched macromolecules were carried out using various methods [2,24,25,26,27,28,29]. Molecular dynamics simulations of isolated stars with various numbers of arms (between 6 and 50) that were rather short (the length of an arm was between 10 and 50) [30,31] were carried out, while the Cooperative Motion Algorithm [32] was used for simulation of star melts with up to 32 arms. Radke simulated isolated stars with 3, 8 and 16 arms, looking for their static properties, including hydrodynamic radius, which was crucial for observations of these macromolecules using size exclusion chromatography [33,34].
A star-branched polymer can be treated as a model of the simplest branched chain, but dendrimers can also be treated as very regular macromolecules. Dendrimers are polymers that are completely and regularly branched [35]. They consist of repeating units emanating from a central branching point [36,37]. They are typically characterized by functionality, which is usually F = 3 or 5, the spacer m (the number of segments between branching points) and the number of generations g (the number of shells containing junctions).
The synthesis of dendrimers is performed according to two main schemes: divergent and convergent, with increasing and constant reactions per growth step, respectively [38,39,40]. Dendrimers are desired because of their potential applications in nanotechnology, drug delivery, supramolecular science, catalysis, as coatings and biosensors, etc. [35,37,38,41,42,43]. The properties of dendrimers differ significantly from those of their linear and star-branched counterparts [44]. Dendrimers with a small number of generations form an entangled melt, while those with larger g form densely packed structures [45]. The number of chain entanglements is smaller than in linear chain systems; thus, their influence on viscoelastic properties is rather small, and the viscosity of dendrimers was found to be a nonmonotonic function of molecular weight [46,47,48].
Theoretical treatment of dendrimers is rather easy due to their regular structure and was mainly carried out using self-consistent mean field and renormalization group calculations [49,50]. The other theoretical approaches are computer simulations [51,52] employing Monte Carlo [53,54,55,56,57,58,59,60] and Molecular Dynamics [61,62,63,64,65,66] methods. The structure and viscoelastic properties of dendrimers studied by computer simulations showed the impact of the spacer length on the polymer properties [67,68]. The structure of dendrimers requires further studies, as it has a great impact on their properties in practical applications such as containers for drugs and nanomaterials [35,36,69,70,71,72].
In this paper, a comparison of the polymerization process of star-branched polymers and dendrimers was simulated by means of Monte Carlo simulations. In contrast to previous studies where initiator, cross-linker and monomer molecules were present [54], the prefabricated, well-defined cores were used in the synthesis. The focus was on the course of the synthesis process and on the scaling laws concerning the dimensions of the chain. The polymerization process was studied on a generic lattice model using the Dynamic Lattice Liquid model. The polymerization process was performed in bulk.

2. Simulation Method

In the presented model of the macromolecular systems, coarse-grained representation was used. The second simplification was the limitation of degrees of freedom by the introduction of a lattice approximation, and a face-centered cubic lattice (fcc) was employed for this purpose. An element of the system like a monomer, a mer in a growing chain (a polymer bead) was allowed to occupy one lattice point. Thus, the size of all of the above objects treated as beads was assumed to be identical. The status of objects could change as a result of the polymerization: each reaction step between beads that were neighbors (located in a distance of lattice constant) was irreversible and resulted in the formation of an unbreakable bond. The attempt of the reaction and of the modification of chains’ conformations takes place simultaneously with the probability of the attachment of a monomer 0.02. The Dynamic Lattice Liquid (DLL) model was used in simulations in order to provide the dynamics to the modeled system [73] with time expressed in algorithm steps (Monte Carlo steps). This model was described in detail in previous publications [20,34,35,71,72]; therefore, only a short description highlighting its main features is given here. In one Monte Carlo step, a field of randomly chosen unit vectors representing motion attempts was generated and assigned to all objects in the system. Only those attempts were successful which coincided in such a way that the sum of displacements along a closed path (loop) that included more than two molecules was equal to zero (condition of continuity). Each bead was then moved to a neighboring lattice site along this loop. All beads that did not contribute to these correlated loops stayed in their previous positions. This model of molecular systems is provided with a dynamics consisting of local vibrations and occasional diffusion steps. Within a longer time interval, this kind of dynamics leads to displacements of individual beads along random walk trajectories with steps distributed randomly in time. The most important advantage of this method is that it naturally takes into account changes in the local mobility of the system components and steric hindrances resulting from polymerization.
The Monte Carlo box was 100 × 100 × 100, and thus the system consisted of 106 elements. The system was athermal (the excluded volume was the only interaction potential), and periodic boundary conditions were applied in all directions. At the beginning of a single simulation run, functional cores were randomly distributed in the system, while all the remaining lattice sites were filled with monomers while initiator molecules were omitted for simplicity. Simulated star-branched polymers were constructed using f = 4, 5, 6, 7, 8, 16 and 32 arms, and each system under consideration consisted of cores with the same functionality only. For stars with f = 4, a core consisted of one 4-functional monomer, but for a larger number of arms, next generations of 3-functional monomers must be added: for f from 5 to 8, the second generation of 3-functional monomers, for f = 16, the third generation, and for f = 32, the fourth generation. The concentrations of cores, i.e., the number of growing macromolecules, were selected in such a way as to obtain stars with approximately the same length of arms: the number of stars with f = 32 arms was twice that of stars with f = 16, four times that of stars with f = 8, etc. This made the concentration of cores 0.032% for stars with 4 arms, 0.016% for stars with 8 arms, etc. The concentration of multifunctional cores in the case of dendrimers’ synthesis was varied in the same way as for stars described above, i.e., 0.032, 0.01, etc. In all simulations of dendrimers, the spacer was set to 1, and the functionality of each monomer was 4. The length of arms in stars and the number of generations in dendrimers were not imposed on the model, and each simulation run lasted until exhaustion of the monomer molecules in the system (up to 106 Monte Carlo steps).

3. Results and Discussion

3.1. Polymerization Kinetics

The kinetics of the synthesis processes for systems containing stars and dendrimers was studied in terms of the conversion of monomers, i.e., the factor related to the growing of macromolecules along with elapsed time. Figure 1a shows the conversion of monomers as a function of time for star-branched polymers with various numbers of arms in the log–log scale. In this case, one can distinguish two clear regimes for all cases, i.e., for the numbers of arms from 4 to 32. Figure 1b presents the derivative of monomer conversion as a function of time, i.e., the reaction rate. In the first regime, the increase in monomer consumption is linear, which means that the total reaction rate is constant, as can be seen in Figure 1b. In the second regime, the conversion of the monomer rapidly slows down, which corresponds to a rapid linear decline in the reaction rate, as presented in Figure 1b. These dramatic changes take place for the monomer conversion around 0.75–0.80, which is marked by the vertical dashed line in Figure 1a. This point indicates the time moment for which free-growing star macromolecules (no limitation of surrounding monomers) begin to compete with each other for available monomer molecules. This implies that the system, which is relatively diluted at first, becomes dense in terms of polymer concentration, and the location of this point very weakly depends on the number of branches. Differences in monomer conversion and the reaction rate between all macromolecular systems considered are small, which can be explained by the fact that the number of growing chains does not change in all star systems under consideration.
Figure 2 shows the monomer conversion and the reaction rate as a function of time for dendrimers. The star-branched polymers and dendrimers cannot be compared directly by means of the number of branches, but they can be compared when systems contain the same number of macromolecules (the same number of multifunctional cores). Figure 2a shows that the conversion of monomers is much faster than in the case of star-branched polymers, and after the time near 103 Monte Carlo simulation steps, all monomer molecules reacted in all presented cases. The reaction rate presented in Figure 2b is slightly higher than for stars and increases in the initial period of time. This can be explained by the fact that for dendrimers, the number of active ends increases sharply as the dendrimer grows, while it remains constant in the case of stars. The overall monomer reaction rate for dendrimers slowly increases in the initial regime and then accelerates until the monomer is completely depleted.
Obviously, these differences in the monomer conversion for stars and dendrimers lead to differences in the growth dynamics and influence their scaling properties. Dynamic scaling of cluster size is usually discussed using the Smoluchowski coagulation equation [74]. It describes the cluster size distribution nl for dilute systems, where only binary collisions are relevant. For long time scales and large clusters, the number of clusters formed by l elements exhibits dynamic scaling, and the cluster size distribution approaches the limit form [75]:
nlN−2 Φ(N/l),
where N = N(t) is related to the number-average cluster size, and Φ(N/l) is a universally time-independent cluster size distribution, which characterizes the aggregation mechanism. Combining the solution of Equation (1) with the fractal scaling relationship leads to:
N(R) = k0Rdf,
where k0 is a prefactor, R is a characteristic size, df stands for fractal dimension and N(R) is the number of particles aggregated in the cluster. The mean cluster size defined in Equation (2) can be also obtained as a function of time [76]:
N(t) ≈ tzt1/(1−λ),
where λ is the van Dongen and Ernst homogeneity exponent [77] that describes different growth kinetics induced by the object–object interaction. It should take the value 0 for the diffusion-limited cluster aggregation and 1 for the reaction-limited cluster aggregation. In the presented case, the role of N(t) plays the number-averaged mass Mw of star or dendrimer (number of mer units).
Figure 3a presents the changes of the average mass Mw of the star-branched chain as a function of time t. Three regimes can be distinguished for all cases. The first is the initial one, where the growth of the star is slow. In the second regime, the mass increases significantly. In the third, the growth becomes slow again due to the high polymer density and the competition between stars in acquiring the remaining monomer. The influence of the arms number becomes visible if one considers the average mass of one arm Mw/f, which is presented in Figure 3b. Three regimes of behavior can also be distinguished in this case. In the first one, the increase of an average arm mass is slow, and a higher number of arms increase their mass. In the second regime, the arm mass increase does not depend on the number of arms and increases with time approximately as t0.84. This exponent is close to the reaction-limited cluster aggregation scenario characterized by exponent 1. In the third regime, the mass of arms stabilizes and does not depend on the number of arms—it scales with time as t0.045. In this regime, the mean size of an arm does not scale according to Equation (3) because the growth of the macromolecule takes place at very high concentrations, while Equation (3) describes the system at low and moderate concentrations.
The changes of dendrimer polymer mass over time are presented in Figure 4. The initial period is common for all cases with a slow increase of mass. Then, a rapid mass increase is observed, which lasts until the monomer is finally depleted. This increase is considerably stronger than in the case of star-branched polymers with the exponent of 2.30.
The kinetics of polymer growth during the polymerization process can be characterized by basic polymerization indicators such as Mw, weight average mass Pw and dispersity Pw/Mw. In this representation, the monomer conversion degree usually plays the role of non-linear time.
Figure 5 shows the number-averaged mass Mw, the weight-averaged mass Pw and the dispersity Pw/Mw for the star-branched systems under consideration in a log–log plot. After a short initial period, where the consumption of the monomer is low and stable, the increase is almost linearly proportional to the monomer conversion in a log–log scale. It was shown that in real experiments (ATRP polymerization), this dependency is the same [78]. The dispersity increases in the initial period, and the more arms the stars have, the lower the dispersity is. In general, the star-branched macromolecules with well-defined cores exhibit very small dispersion (below 1.07).
Figure 6 shows the number-averaged mass, the weight-averaged mass and the dispersity for the dendrimer systems under consideration in a log–log plot. In this case, the changes of Mw and Pw are similar to those for star-branched polymers, as discussed above. The dispersity almost does not depend on the number of macromolecules and increases with monomer conversion. This parameter is slightly higher than for star-branched polymers, and it remains rather low.
In order to obtain insight into details of dispersity values in both discussed cases, total mass distributions must be examined. Figure 7 presents the polymer mass distributions of macromolecules at the end of the polymerization process for all systems under consideration fitted with a Gaussian function. One can see that, in spite of small differences in the dispersity of stars and dendrimers (Figure 5 and Figure 6), the distributions of mass for dendrimers are considerably broader. For stars, the width of the distribution is almost constant, while for dendrimers it decreases when the number of macromolecules increases.

3.2. Molecular Topology during the Polymerization Process

The basic quantities used to characterize the geometric properties of a macromolecule are the radius of gyration and the hydrodynamic radius. In the presented cases, for low dispersion, one can define the mean-squared radius of gyration < Rg2 > and the mean hydrodynamic radius < Rh > as follows:
< Rg2 > = 1/Mj R2g,j,
R2g,j = 1/Nji (ri,jrcm,j)2,
where M is the total number of macromolecular objects in the system (stars or dendrimers), Nj is the number of beads in a given macromolecule and rcm,j is the location of its center of mass.
< Rh > = 1/Mj Rh,j,
1/Rh,j = 1/2 ∑mn 1/rmn,j,
where rm,n is the distance between two different beads m and n in the molecule j.
The radius of gyration and hydrodynamic radius were calculated for both architectures studied according to Equations (4) and (5). The < Rg >/< Rh > ratios for both stars and dendrimers are presented in Figure 8. The results clearly reflect this ratio for stars, i.e., the ratio decreases with the number of arms, although it differs from the experimental data [79]. For small dendrimers, the ratio does not depend on the number of macromolecules in the system, while for larger ones it diverges.
The shape of macromolecules is usually described by means of the asphericity factor, which is defined as
A = λ 2 λ 1 2 + λ 3 λ 1 2 + λ 2 λ 3 2 λ 1 + λ 2 + λ 3 2
where λi are the components of the squared radius of gyration along the principal axes of inertia [80]. Figure 9a presents the mean asphericity factor A determined according to Equation (6). The shape of stars becomes more spherical when the number of arms increases, which is consistent with other findings [2]. In the case of dendrimers, there is a minimum on the curve, which was already found by experiments [40]. Polymer density profiles for both architectures are presented in Figure 9b. Polymer density decreases very slowly, which suggests the presence of a dense core and a rapid decrease for longer distances. In their theoretical considerations, Daoud and Cotton [81] predicted two regimes inside stars: an unswollen one with slope −1 and a swollen one with slope −4/3. One can conclude that it is impossible to localize the regions with such slopes on our curves.
Another important aspect of studying the topology of macromolecular objects is that their geometry changes with the synthesis progress with respect to both time and total mass (total number of beads N). Typically, for stars, the mass scaling is expressed in the following form:
< R2 > ≈ N2νN2/df,
where R is a characteristic dimension of the entire object (e.g., radius of gyration), ν the scaling exponent and df is the object fractal dimension.
Figure 10 shows the changes in the mean-squared radius of gyration as a function of star mass and time. A scaling behavior of < Rg2 > was found not to be related to the number of arms. The exponent 2ν decreases from 1.01 to 0.83, as the number of arms increases, which means that the fractal dimension becomes higher (from 1.98 to 2.41) with the number of arms. The exponent 2ν from Equation (7) for dense polymer melts is expected to be close to 1 [82], while lower values are predicted for dense collapsed structures. The size growth exponent is similar for all investigated numbers of arms.
In the case of dendrimers, the scaling of chain size is more complex. Using Flory theory, it was shown that under good solvent conditions, the radius of gyration should scale according to the following formula [50,53,54]:
RgN1/5 [(g + 1)m]2/3,
where g is the number of generations and m is the length of the spacer. Using simple arguments based on the properties of a geometric series (in investigated cases, the beads were 4-functional, with m = 1), it can be stated that the number of beads N corresponding to the generation g is given by the following formula [49,53]:
N = 1 + 2(3g+1 − 1),
By solving Equation (9) with respect to the number of generations g, one can easily find g corresponding to number (mass) N of the dendrimer′s elements.
g = log3((N + 1)/2) − 1,
Figure 11 shows the changes in the number of generations g calculated according to Equation (10) as a function of time of synthesis and the total mass of a macromolecules. In the initial period, the estimated number of generations does not depend on the number of growing macromolecules. Moreover, it is clear that the increased scales with time and the scaling exponent are found to be close to 0.47. Then, well below 103 Monte Carlo steps, g clearly depends on the number of cores (macromolecules) in the system. At the lowest initiator concentration (the smallest number of cores), almost 8 generation dendrimers were achieved, while for the highest number of cores, the number of generations does not exceed 6 (one must remember that the polymers studied are not the same length and the mean number of generations is studied here). The dependence of the number of generations on the total dendrimer mass is presented in Figure 11, showing no differences in the investigated cases.
To verify the validity of these theoretical considerations, a reduced radius of gyration Rg/((g + 1) − 3)2/5 was calculated based on Equation (8). Figure 12 presents it as a function of the total number of beads N. The scaling exponent was found to be 0.21 for all cases under consideration, and very close to the value 1/5 predicted by Flory theory given by Equation (8). This is consistent with the simulation results obtained for dendrimer lattice models at good solvent conditions [53,83]. The deviation from this value is observed when macromolecules become larger and interchain interactions become an important factor (as the polymer concentration approaches 1), which makes the dependence stronger. Figure 12 indicates that the size of dendrimers Rg2 does not scale ideally with the number of beads N but also with the topological parameter g. This means that dendrimers obtained in the simulated synthesis process do not possess a well-defined fractal dimension, and thus, are not true fractal objects (for a model dendrimer, the fractal dimension should be close to 3). This behavior is consistent with Flory’s theory and was observed in other simulation studies [53].
The visualization of a single macromolecule clearly shows the differences between the studied polymer architectures. Figure 13 presents examples of single star conformations with a small number of arms (f = 4) and a large number of arms (f = 32) and the corresponding dendrimer structures. One can observe the formation of a dense core for the star with a high number of arms. The difference between stars and dendrimers is also striking as the entire dendrimer must be considered as a densely packed globule.

4. Conclusions

The simulations of a generic model of star and dendrimer polymers were performed by means of the Monte Carlo method. The Dynamic Lattice Liquid model was employed for studies of the bulk polymerization process and characterization of structural properties. The model system contained well-defined cores where the growing of chains (the attaching of multifunctional monomers) takes place. The polymerization was carried out in bulk to the exhaustion of monomers with no solvent present.
The simulation results clearly indicated that the process of synthesis of highly branched macromolecules like stars and dendrimers using well-defined cores led to macromolecular structures with well-defined properties. It was shown that, for both macromolecular architectures and for all discussed cases (number of arms for stars and number of macromolecules for dendrimers), the dispersity was found to remain very low during the entire process of polymerization. Thus, this parameter was found to be considerably lower than for systems where cores were formed from initiator and cross-linker molecules simultaneously with the process of polymerization, i.e., attaching monomers to growing macromolecules. The mass of stars and dendrimers was rather comparable. The kinetics of the synthesis process for stars and dendrimers was compared, showing considerably faster polymerization of dendrimers, evidently because of a higher (and growing) number of points where a macromolecule could attach monomers. The reaction-limited cluster aggregation scenario was found to more suitable for stars. The size and structure of star-branched polymers and dendrimers during synthesis were found to be different, and the scaling of polymer size with its total number of beads was different too.

Author Contributions

Conceptualization, P.P.; methodology, K.H. and A.S.; validation, K.H., A.S. and P.P.; formal analysis, K.H.; investigation, K.H., A.S. and P.P.; data curation, K.H. and P.P.; writing—original draft preparation, A.S., K.H. and P.P.; writing—review and editing, K.H. and A.S.; visualization, K.H.; supervision, A.S. and P.P.; project administration, A.S.; funding acquisition, A.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the Polish National Science Center grant UMO-2017/25/B/ST5/01970.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fetters, L.J.; Kiss, A.D.; Pearson, D.S.; Quack, G.F.; Vitus, F.J. Rheological behavior of star-shaped polymers. Macromolecules 1993, 26, 647–654. [Google Scholar] [CrossRef]
  2. Grest, G.S.; Fetters, L.J.; Huang, J.S.; Richter, D. Star polymers: Experiment, theory, and simulation. Adv. Chem. Phys. 1996, 94, 67–163. [Google Scholar]
  3. Vlassopoulos, D.; Fytas, G.; Pakula, T.; Roovers, J. Multiarm star polymers dynamics. J. Phys. Condens. Matter 2001, 13, R855–R876. [Google Scholar] [CrossRef]
  4. Seo, S.E.; Hawker, C.J. The beauty of branching in polymer science. Macromolecules 2020, 53, 3257–3261. [Google Scholar] [CrossRef]
  5. Wu, W.; Wang, W.; Li, J. Star polymers: Advances in biomedical applications. Prog. Polym. Sci. 2015, 46, 55–85. [Google Scholar] [CrossRef]
  6. Likos, C.N. Soft matter with soft particles. Soft Matter 2006, 2, 478–498. [Google Scholar] [CrossRef]
  7. Gao, H.F.; Matyjaszewski, K. Structural control in ATRP synthesis of star polymers using the arm-first method. Macromolecules 2006, 39, 3154–3160. [Google Scholar] [CrossRef]
  8. Gao, H.; Matyjaszewski, K. Low-polydispersity star polymers with core functionality by cross-linking macromonomers using functional ATRP initiators. Macromolecules 2007, 40, 399–401. [Google Scholar] [CrossRef]
  9. Gao, H.; Ohno, S.; Matyjaszewski, K. Low polydispersity star polymers via cross-linking macromonomers by ATRP. J. Am. Chem. Soc. 2006, 128, 15111–15113. [Google Scholar] [CrossRef]
  10. Gao, H.; Matyjaszewski, K. Synthesis of star polymers by a new “Core-First” method: Sequential polymerization of cross-linker and monomer. Macromolecules 2008, 41, 1118–1125. [Google Scholar] [CrossRef]
  11. Gao, H.; Matyjaszewski, K. Synthesis of functional polymers with controlled architecture by CRP of monomers in the presence of cross-linkers: From stars to gels. Prog. Polym. Sci. 2009, 34, 317–350. [Google Scholar] [CrossRef]
  12. Blencowe, A.; Tan, J.F.; Goh, T.K.; Qiao, G.G. Core cross-linked star polymers via controlled radical polymerization. Polymer 2009, 50, 5–32. [Google Scholar] [CrossRef] [Green Version]
  13. Gao, H.; Miasnikova, A.; Matyjaszewski, K. Effect of cross-linker reactivity on experimental gel points during ATRcP of monomer and cross-linker. Macromolecules 2008, 41, 7843–7849. [Google Scholar] [CrossRef]
  14. Rosselgong, J.; Armes, S.P.; Barton, W.R.S.; Price, D. Synthesis of branched methacrylic copolymers: Comparison between RAFT and ATRP and effect of varying the monomer concentration. Macromolecules 2010, 43, 2145–2156. [Google Scholar] [CrossRef]
  15. Liu, J.; Duong, H.; Whittaker, M.R.; Davis, T.P.; Boyer, C. Synthesis of functional core, star polymers via RAFT polymerization for drug delivery applications. Macromol. Theory Simul. 2012, 33, 760–766. [Google Scholar] [CrossRef]
  16. Ren, J.M.; McKenzie, T.G.; Fu, Q.; Wong, E.H.H.; Xu, J.; An, Z.; Shanmugam, S.; Davis, T.P.; Boyer, C.; Qiao, G.G. Star polymers. Chem. Rev. 2016, 116, 6743–6836. [Google Scholar] [CrossRef]
  17. Gao, H.; Matyjaszewski, K. Synthesis of miktoarm star polymers via ATRP using the “in-out” method: Determination of initiation efficiency of star macroinitiators. Macromolecules 2006, 39, 7216–7223. [Google Scholar] [CrossRef]
  18. Li, W.; Matyjaszewski, K. Star polymers via cross-linking amphiphilic macroinitiators by AGET ATRP in aqueous media. J. Am. Chem. Soc. 2009, 131, 10378–10379. [Google Scholar] [CrossRef]
  19. Li, W.; Yoon, J.A.; Zhong, M.; Matyjaszewski, K. Atom Transfer Radical Copolymerization of monomer and cross-linker under highly dilute conditions. Macromolecules 2011, 44, 3270–3275. [Google Scholar] [CrossRef]
  20. Rosselgong, J.; Armes, S.P. Quantification of intramolecular cyclization in branched copolymers by 1H NMR spectroscopy. Macromolecules 2012, 45, 2731–2737. [Google Scholar] [CrossRef]
  21. Polanowski, P.; Jeszka, J.K.; Li, W.; Matyjaszewski, K. Effect of dilution on branching and gelation in living copolymerization of monomer and divinyl cross-linker: Modeling using dynamic lattice liquid model (DLL) and Flory-Stockmayer (FS) model. Polymer 2011, 52, 5092–5101. [Google Scholar] [CrossRef]
  22. Polanowski, P.; Jeszka, J.K.; Matyjaszewski, K. Modeling of branching and gelation in living copolymerization of monomer and divinyl cross-linker using dynamic lattice liquid model (DLL) and Flory-Stockmayer model. Polymer 2010, 51, 6084–6092. [Google Scholar] [CrossRef]
  23. Freire, J.J. Conformational properties of branched polymers: Theory and simulations. Adv. Polym. Sci. 1999, 143, 35–112. [Google Scholar]
  24. Mavrantzas, V.G. Using Monte Carlo to simulate complex polymer systems: Recent progress and outlook. Front. Phys. 2021, 9, 661367. [Google Scholar] [CrossRef]
  25. Hsu, H.-P.; Nadler, W.; Grassberger, P. Scaling of star polymers with 1–80 arms. Macromolecules 2004, 37, 4658–4663. [Google Scholar] [CrossRef] [Green Version]
  26. Dasmahapatra, A.K.; Reddy, G.K. Conformational transition of telechlic star polymers. Polymer 2013, 54, 2392–2400. [Google Scholar] [CrossRef]
  27. Huissmann, S.; Blaak, R.; Likos, C.N. Star polymers in solvents of varying quality. Macromolecules 2009, 42, 2806–2816. [Google Scholar] [CrossRef]
  28. Holler, S.; Moreno, A.J.; Zamponi, M.; Bacova, P.; Willner, L.; Iatrou, H.; Falus, P.; Richter, D. The role of the functionality in the branch point motion in symmetric star polymers: A combined study by simulations and neutron spin echo spectroscopy. Macromolecules 2018, 51, 242–253. [Google Scholar] [CrossRef] [Green Version]
  29. Chremos, A.; Glynos, E.; Green, P.F. Structure and dynamics intra-molecular heterogenity od star polymer melts above glass transition temperature. J. Chem. Phys. 2015, 142, 044901. [Google Scholar] [CrossRef]
  30. Grest, G.S.; Kremer, K.; Witten, T.A. Structure of many-arm star polymers—A Molecular-Dynamics simulation. Macromolecules 1987, 20, 1376–1383. [Google Scholar] [CrossRef]
  31. Grest, G.S. Structure of many-arm star polymers in solvents of varying quality—A Molecular-Dynamics study. Macromolecules 1994, 27, 3493–3500. [Google Scholar] [CrossRef]
  32. Pakula, T. Static and dynamic properties of computer simulated melts of multiarm polymer stars. Comput. Theor. Polym. Sci. 1998, 8, 21–30. [Google Scholar] [CrossRef]
  33. Radke, W. Simulation of GPC-distribution coefficients of linear and star-shaped molecules in spherical pores. Macromol. Theory Simul. 2001, 10, 668–675. [Google Scholar] [CrossRef]
  34. Radke, W.; Gerber, J.; Wittmann, G. Simulation of GPC-distribution coefficients of linear and star-shaped molecules in spherical pores. 2. Comparison of simulation and experiment. Polymer 2003, 44, 519–525. [Google Scholar] [CrossRef]
  35. Astruc, D.; Boisselier, E.; Ornelas, C. Dendrimers designed for functions: From physical, photophysical, and supramolecular properties to applications in sensing, catalysis, molecular electronics, photonics, and nanomedicine. Chem. Rev. 2010, 110, 1857–1959. [Google Scholar] [CrossRef] [PubMed]
  36. Tomalia, D.A. Birth of a new macromolecular architecture: Dendrimers as quantized building blocks for nanoscale synthetic polymer chemistry. Prog. Polym. Sci. 2005, 30, 294–324. [Google Scholar] [CrossRef]
  37. Tomalia, D.A. Dendrons/Dendrimers: Quantized, nano-element like building blocks for soft-soft and soft-hard nano-compound synthesis. Soft Matter 2009, 6, 456–474. [Google Scholar] [CrossRef]
  38. Tomalia, D.A.; Frechet, J.M.J. Discovery of dendrimers and dendritic polymers: A brief historical perspective. J. Polym. Sci. Part A Polym. Chem. 2002, 40, 2719–2728. [Google Scholar] [CrossRef]
  39. Hawker, C.J.; Frechet, J.M.J. Preparation of polymers with controlled molecular architecture: A new convergent approach to dendritic macromolecules. J. Am. Chem. Soc. 1990, 112, 7638–7647. [Google Scholar] [CrossRef]
  40. Tomalia, D.A.; Christensen, J.B.; Boas, U. Dendrimers, Dendrons and dendritic Polymers: Discovery, Applications, and the Future; Cambridge University Press: Cambridge, UK, 2012. [Google Scholar]
  41. Patri, A.K.; Majoros, I.J.; Baker, J.R. Dendritic polymer macromolecular carriers for drug delivery. Curr. Opin. Chem. Biol. 2002, 6, 466–471. [Google Scholar] [CrossRef]
  42. Hourani, R.; Kakkar, A. Advances in the elegance of chemistry in designing dendrimers. Macromol. Rapid Commun. 2010, 31, 947–974. [Google Scholar] [CrossRef] [PubMed]
  43. Inoue, K. Functional dendrimers, hyperbranched and stars polymers. Prog. Polym. Sci. 2000, 25, 453–471. [Google Scholar] [CrossRef]
  44. Maiti, P.M.; Çağın, T.; Wang, G.; Goddard, W.A., III. Structure of PAMAM dendrimers: Generations 1 through 11. Macromolecules 2004, 37, 6236–6254. [Google Scholar] [CrossRef]
  45. Pilkington, G.A.; Pedersen, J.S.; Briscoe, W.H. Dendrimer nanofluids in the concentrated regime: From polymer melts to soft spheres. Langmuir 2015, 31, 3333–3342. [Google Scholar] [CrossRef]
  46. Frechet, J.M.J. Functional polymers and sendrimers: Reactivity, molecular architecture, and interfacial energy. Science 1994, 263, 1710–1715. [Google Scholar] [CrossRef]
  47. Mourey, T.H.; Turner, S.R.; Rubinstein, M.; Frechet, J.M.J.; Hawker, C.J.; Wooley, K.L. Unique behavior of dendritic macromolecules: Intrinsic viscosity of polyether Dendrimers. Macromolecules 1992, 25, 2401–2406. [Google Scholar] [CrossRef]
  48. Bosman, A.W.; Janssen, H.M.; Meijer, E.W. About dendrimers: Structure, physical properties, and applications. Chem. Rev. 1999, 99, 1665–1688. [Google Scholar] [CrossRef]
  49. La Ferla, R. Conformations and dynamics of dendrimers and cascade macromolecules. J. Chem. Phys. 1997, 106, 688–700. [Google Scholar] [CrossRef]
  50. Sheng, Y.J.; Jiang, S.Y.; Tsao, H.K. Radial size of a starburst dendrimer in solvents of varying quality. Macromolecules 2002, 35, 7865–7868. [Google Scholar] [CrossRef]
  51. Likos, C.N.; Ballauff, M. Equilibrium structure of dendrimer: Results and open questions. Top. Curr. Chem. 2005, 245, 239–252. [Google Scholar]
  52. Ballauff, M.; Likos, C.N. Dendrimers in solution: Insight from theory and simulation. Angew. Chem. Int. Ed. 2004, 43, 2998–3020. [Google Scholar] [CrossRef] [PubMed]
  53. Giupponi, G.; Buzza, D.M.A. Monte Carlo simulation of dendrimers in variable solvent quality. J. Chem. Phys. 2004, 120, 10290–10298. [Google Scholar] [CrossRef] [PubMed]
  54. Kłos, J.S.; Sommer, J.U. Properties of dendrimers with flexible spacer-chains: A Monte Carlo study. Macromolecules 2009, 42, 4878–4886. [Google Scholar] [CrossRef]
  55. Wawrzyńska, E.; Eisenhaber, S.; Parzuchowski, P.; Sikorski, A.; Zifferer, G. Simulation of hyperbranched polymers. 1. Properties of regular three generation dendrimers. Macromol. Theory Simul. 2014, 23, 288–299. [Google Scholar] [CrossRef]
  56. Wawrzyńska, E.; Sikorski, A.; Zifferer, G. Monte Carlo simulation studies of regular and irregular dendritic polymers. Macromol. Theory Simul. 2015, 24, 477–489. [Google Scholar] [CrossRef]
  57. Romiszowski, P.; Sikorski, A. Shape of star-branched polymers at various solvent conditions. A computer simulation study. J. Chem. Phys. 1998, 109, 6169–6174. [Google Scholar] [CrossRef]
  58. Sikorski, A. Properties of star-branched polymer chains. The application of the Replica Exchange Monte Carlo method. Macromolecules 2002, 35, 7132–7137. [Google Scholar] [CrossRef]
  59. Sikorski, A. Monte Carlo study of the collapse transition of flexible and semiflexible star-branched polymers. Polymer 1993, 34, 1271–1281. [Google Scholar] [CrossRef]
  60. Polanowski, P.; Jeszka, J.K.; Matyjaszewski, K. Synthesis of star polymers by “core-first” one-pot method via ATRP: Monte Carlo simulations. Polymer 2014, 55, 2552–2561. [Google Scholar] [CrossRef]
  61. Tian, W.D.; Ma, Y.Q. Coarse-grained molecular simulation of interacting dendrimers. Soft Matter 2011, 7, 500–505. [Google Scholar] [CrossRef]
  62. Kanchi, S.; Suresh, S.; Priyakumar, U.D.; Ayappa, K.G.; Maiti, P.K. Molecular Dynamics study of the structure, flexibility, and hydrophilicity of PETIM dendrimers: A comparison with PAMAM dendrimers. J. Phys. Chem. B 2015, 119, 12990–13001. [Google Scholar] [CrossRef] [PubMed]
  63. Karatasos, K.; Adolf, D.B.; Davies, G.R. Statics and dynamics of model dendrimers as studied by Molecular Dynamics simulations. J. Chem. Phys. 2001, 115, 5310–5318. [Google Scholar] [CrossRef] [Green Version]
  64. Kurbatov, A.O.; Balabaev, N.K.; Mazo, M.A.; Kramarenko, E.Y. Molecular Dynamics simulations of single siloxane dendrimers: Molecular structure and intramolecular mobility of terminal groups. J. Chen. Phys. 2018, 148, 014902. [Google Scholar] [CrossRef] [PubMed]
  65. Markelov, D.A.; Shishkin, A.N.; Matveev, V.V.; Penkova, A.V.; Lähderanta, E.; Chizhik, V.I. Orientational mobility in dendrimer melts: Molecular Dynamics simulations. Macromolecules 2018, 49, 9247–9257. [Google Scholar] [CrossRef]
  66. Khabaz, F.; Khare, R. Effect of chain architecture on the size, shape, and intrinsic viscosity of chains in polymer solutions: A molecular simulation study. J. Chem. Phys. 2014, 141, 214904. [Google Scholar] [CrossRef]
  67. Mansfield, M.L.; Klushin, L.I. Monte Carlo studies of dendrimer macromolecules. Macromolecules 1993, 26, 4262–4268. [Google Scholar] [CrossRef]
  68. Timoshenko, E.G.; Kuznetsov, Y.A.; Connolly, R. Conformations of dendrimers in dilute solution. J. Chem. Phys. 2002, 117, 9050–9062. [Google Scholar] [CrossRef] [Green Version]
  69. Kesharwani, P.; Jain, K.; Jain, N.K. Dendrimer as nanocarrier for drug delivery. Prog. Polym. Sci. 2014, 39, 268–307. [Google Scholar] [CrossRef]
  70. Mignani, S.; El Kazzouli, S.; Bousmina, M.; Majoral, J.P. Dendrimer space concept for innovative nanomedicine: A futuristic vision for medicinal chemistry. Prog. Polym. Sci. 2013, 38, 993–1008. [Google Scholar] [CrossRef]
  71. Polanowski, P.; Pakula, T. Studies of polymer conformation and dynamics in two dimensions using simulations based on the Dynamic Lattice Liquid (DLL) model. J. Chem. Phys. 2002, 117, 4022–4029. [Google Scholar] [CrossRef]
  72. Hałagan, K.; Banaszak, M.; Jung, J.; Polanowski, P.; Sikorski, A. The structure of opposing polymer brushes. Polymers 2021, 13, 4294. [Google Scholar] [CrossRef] [PubMed]
  73. Pakula, T. Simulation on the completely occupied lattices. In Simulation Methods for Polymers; Kotelyanskii, M., Theodorou, D.N., Eds.; Marcel Dekker: New York, NY, USA; Basel, Switzerland, 2004. [Google Scholar]
  74. von Smoluchowski, M. Versuch einer Mathematischen Theorie der Koagulations Kinetic Kolloider Lousungen. Phys. Z. Chem. 1917, 92, 129–168. [Google Scholar]
  75. Jullien, R.; Botet, R. Aggregation and Fractal Aggregate; World Scientific: Singapore, 1987. [Google Scholar]
  76. Tirado-Miranda, M.; Schmitt, A.; Callejas-Fernández, J.; Frenández-Barbero, A. Dynamic scaling and fractal structure of small colloidal clusters. Colloids Surf. A 2000, 162, 67–73. [Google Scholar] [CrossRef]
  77. Van Dongen, P.; Ernst, M.H. Dynamic scaling in the kinetics of clustering. Phys. Rev. Lett. 1985, 54, 1396–1399. [Google Scholar] [CrossRef] [PubMed]
  78. Rubinstein, M.; Colby, R.H. Polymer Physics; Oxford University Press: Oxford, UK, 2003. [Google Scholar]
  79. Krys, P.; Matyjaszewski, K. Kinetics of Atom Transfer Radical Polymerization. Eur. Polym. J. 2017, 89, 482–523. [Google Scholar] [CrossRef] [Green Version]
  80. Rudnick, J.; Gaspari, G. The asphericity of random walks. J. Phys. A Math. Gen. 1986, 19, L191–L193. [Google Scholar] [CrossRef]
  81. Daoud, M.; Cotton, J.P. Star shaped polymers: A model for the conformation and its concentration dependence. J. Phys. 1982, 43, 531–538. [Google Scholar] [CrossRef]
  82. Teraoka, I. Polymer Solutions: An Introduction to Physical Properties; Wiley: New York, NY, USA, 2002. [Google Scholar]
  83. Peng, X.; Cheng, L.; Jiang, B.; Ji, S.; Shen, X.-C.; Williams, D.R.M. The branching angle effect on the properties of rigid dendrimers studied by Monte Carlo simulation. J. Mol. Model. 2021, 27, 144. [Google Scholar] [CrossRef]
Figure 1. Conversion of a monomer (a) and its derivative (b) as a function of time during polymerization of star-branched polymers. The number of arms is given in the (b) legend.
Figure 1. Conversion of a monomer (a) and its derivative (b) as a function of time during polymerization of star-branched polymers. The number of arms is given in the (b) legend.
Polymers 14 02522 g001
Figure 2. Conversion of a monomer (a) and its derivative (b) as a function of time during polymerization of dendrimers. The color indicates the number of macromolecules in the system (see text for details).
Figure 2. Conversion of a monomer (a) and its derivative (b) as a function of time during polymerization of dendrimers. The color indicates the number of macromolecules in the system (see text for details).
Polymers 14 02522 g002
Figure 3. Average total mass (a) and average mass of arm (b) of star-branched polymer as a function of time. The number of arms is given in the legend.
Figure 3. Average total mass (a) and average mass of arm (b) of star-branched polymer as a function of time. The number of arms is given in the legend.
Polymers 14 02522 g003
Figure 4. The average total mass of the dendrimer as a function of time. The color indicates the number of macromolecules in the system (see text for details).
Figure 4. The average total mass of the dendrimer as a function of time. The color indicates the number of macromolecules in the system (see text for details).
Polymers 14 02522 g004
Figure 5. Number-averaged mass (upper panel), weight-averaged mass (middle panel) and dispersity (lower panel) of star-branched polymers as a function of the monomer conversion. The number of arms is given in the legend.
Figure 5. Number-averaged mass (upper panel), weight-averaged mass (middle panel) and dispersity (lower panel) of star-branched polymers as a function of the monomer conversion. The number of arms is given in the legend.
Polymers 14 02522 g005
Figure 6. Number-averaged mass (upper panel), weight-averaged mass (middle panel) and dispersity (lower panel) for dendrimers as a function of the monomer conversion. The color indicates the number of macromolecules in the system (see text for details).
Figure 6. Number-averaged mass (upper panel), weight-averaged mass (middle panel) and dispersity (lower panel) for dendrimers as a function of the monomer conversion. The color indicates the number of macromolecules in the system (see text for details).
Polymers 14 02522 g006
Figure 7. The mass distribution for star-branched polymers (open symbols) and dendrimers (solid symbols). The lines indicate Gaussian fits.
Figure 7. The mass distribution for star-branched polymers (open symbols) and dendrimers (solid symbols). The lines indicate Gaussian fits.
Polymers 14 02522 g007
Figure 8. The ratio < Rg >/< Rh > as a function of total number of beads N for stars (a) and dendrimers (b). The color indicates the number of macromolecules in the system (see text for details).
Figure 8. The ratio < Rg >/< Rh > as a function of total number of beads N for stars (a) and dendrimers (b). The color indicates the number of macromolecules in the system (see text for details).
Polymers 14 02522 g008
Figure 9. Asphericity as a function of polymer mass (a) and density profiles (b) for star-branched chains and dendrimers (see text for details).
Figure 9. Asphericity as a function of polymer mass (a) and density profiles (b) for star-branched chains and dendrimers (see text for details).
Polymers 14 02522 g009
Figure 10. The mean-squared radius of gyration < Rg2 > of star-branched chains as a function of total number of beads N (a) and as a function of time (b). The colors indicate the number of arms.
Figure 10. The mean-squared radius of gyration < Rg2 > of star-branched chains as a function of total number of beads N (a) and as a function of time (b). The colors indicate the number of arms.
Polymers 14 02522 g010
Figure 11. The number of dendrimer generations g as a function of time (a) and total number of beads N (b). The color indicates the number of macromolecules in the system (see text for details).
Figure 11. The number of dendrimer generations g as a function of time (a) and total number of beads N (b). The color indicates the number of macromolecules in the system (see text for details).
Polymers 14 02522 g011
Figure 12. The reduced radius of gyration for dendrimers Rg/((g + 1) − 3)2/5 as a function of total number of beads N (a) and time (b). The color indicates the number of macromolecules in the system (see text for details).
Figure 12. The reduced radius of gyration for dendrimers Rg/((g + 1) − 3)2/5 as a function of total number of beads N (a) and time (b). The color indicates the number of macromolecules in the system (see text for details).
Polymers 14 02522 g012
Figure 13. Snapshots of randomly selected star-branched polymers and dendrimers. The stars consisted of f = 4 (a) and f = 32 (b). The dendrimers (c,d) represent systems containing the same number of cores (macromolecules) as those with stars presented in the upper row.
Figure 13. Snapshots of randomly selected star-branched polymers and dendrimers. The stars consisted of f = 4 (a) and f = 32 (b). The dendrimers (c,d) represent systems containing the same number of cores (macromolecules) as those with stars presented in the upper row.
Polymers 14 02522 g013
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Polanowski, P.; Hałagan, K.; Sikorski, A. Star Polymers vs. Dendrimers: Studies of the Synthesis Based on Computer Simulations. Polymers 2022, 14, 2522. https://doi.org/10.3390/polym14132522

AMA Style

Polanowski P, Hałagan K, Sikorski A. Star Polymers vs. Dendrimers: Studies of the Synthesis Based on Computer Simulations. Polymers. 2022; 14(13):2522. https://doi.org/10.3390/polym14132522

Chicago/Turabian Style

Polanowski, Piotr, Krzysztof Hałagan, and Andrzej Sikorski. 2022. "Star Polymers vs. Dendrimers: Studies of the Synthesis Based on Computer Simulations" Polymers 14, no. 13: 2522. https://doi.org/10.3390/polym14132522

APA Style

Polanowski, P., Hałagan, K., & Sikorski, A. (2022). Star Polymers vs. Dendrimers: Studies of the Synthesis Based on Computer Simulations. Polymers, 14(13), 2522. https://doi.org/10.3390/polym14132522

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop